Computer-aided Drug Design Approaches Research Articles

A computational biology approach for the identification of potential SARS-CoV-2 main protease inhibitors from natural essential oil compounds.

Background: Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) has fomented a climate of fear worldwide due to its rapidly spreading nature, and high mortality rate. The World Health Organization (WHO) declared it as a global pandemic on 11th March, 2020. Many endeavors have been made to find appropriate medications to restrain the SARS CoV-2 infection from spreading but there is no specific antiviral therapy to date. However, a computer-aided drug design approach can be an alternative to identify probable drug candidates within a short time. SARS-CoV-2 main protease is a proven drug target, and it plays a pivotal role in viral replication and transcription. Methods: In this study, we identified a total of 114 essential oil compounds as a feasible anti-SARS-CoV-2 agent from several online reservoirs. These compounds were screened by incorporating ADMET profiling, molecular docking, and 50 ns of molecular dynamics simulation to identify potential drug candidates against the SARS-CoV-2 main protease. The crystallized SARS-CoV-2 main protease structure was collected from the RCSB PDB database (PDB ID 6LU7). Results: According to the results of the ADMET study, none of the compounds have any side effects that could reduce their druglikeness or pharmacokinetic properties. Out of 114 compounds, we selected bisabololoxide B, eremanthin, and leptospermone as our top drug candidates based on their higher binding affinity scores, and strong interaction with the Cys 145-His 41 catalytic dyad. Finally, the molecular dynamics simulation was implemented to evaluate the structural stability of the ligand-receptor complex. MD simulations disclosed that all the hits showed conformational stability compared to the positive control α-ketoamide. Conclusions: Our study showed that the top three hits might work as potential anti-SARS-CoV-2 agents, which can pave the way for discovering new drugs, but for experimental validation, they will require more in vivo trials.

Cheminformatics-based analysis identified (Z)-2-(2,5-dimethoxy benzylidene)-6-(2-(4-methoxyphenyl)-2-oxoethoxy) benzofuran-3(2H)-one as an inhibitor of Marburg replication by interacting with NP

The highly pathogenic Marburg virus (MARV) is a member of the Filoviridae family, a non-segmented negative-strand RNA virus. This article represents the computer-aided drug design (CADD) approach for identifying drug-like compounds that prevent the MARV virus disease by inhibiting nucleoprotein, which is responsible for their replication. This study used a wide range of in silico drug design techniques to identify potential drugs. Out of 368 natural compounds, 202 compounds passed ADMET, and molecular docking identified the top two molecules (CID: 1804018 and 5280520) with a high binding affinity of −6.77 and −6.672 kcal/mol, respectively. Both compounds showed interactions with the common amino acid residues SER_216, ARG_215, TYR_135, CYS_195, and ILE_108, which indicates that lead compounds and control ligands interact in the common active site/catalytic site of the protein. The negative binding free energies of CID: 1804018 and 5280520 were −66.01 and −31.29 kcal/mol, respectively. Two lead compounds were re-evaluated using MD modeling techniques, which confirmed CID: 1804018 as the most stable when complexed with the target protein. PC3 of the (Z)-2-(2,5-dimethoxybenzylidene)-6-(2-(4-methoxyphenyl)-2-oxoethoxy) benzofuran-3(2H)-one (CID: 1804018) was 8.74 %, whereas PC3 of the 2′-Hydroxydaidzein (CID: 5280520) was 11.25 %. In this study, (Z)-2-(2,5-dimethoxybenzylidene)-6-(2-(4-methoxyphenyl)-2-oxoethoxy) benzofuran-3(2H)-one (CID: 1804018) unveiled the significant stability of the proteins' binding site in ADMET, Molecular docking, MM-GBSA and MD simulation analysis studies, which also showed a high negative binding free energy value, confirming as the best drug candidate which is found in Angelica archangelica which may potentially inhibit the replication of MARV nucleoprotein.

Revolutionizing Drug Discovery: A Comprehensive Review of Computer-Aided Drug Design Approaches

Abstract: Computer-Aided Drug Design (CADD) has significantly advanced the drug discovery process, offering tools to enhance efficiency and reduce costs. This review explores essential CADD methodologies, including molecular docking, virtual screening, ADMET profiling, homology modeling, and Quantitative Structure-Activity Relationship (QSAR) models. Molecular docking predicts interactions between drugs and targets, while virtual screening evaluates large compound libraries to identify promising candidates. ADMET profiling assesses pharmacokinetic and toxicological properties early in development. Homology modeling constructs three-dimensional protein models to aid target identification, and QSAR models predict biological activities based on chemical structures. These integrated approaches streamline drug development, providing a robust framework for modern pharmaceutical research.

Targeting small druggable compounds against 3RZE histamine H1 receptor as potential of anti-allergic drug applying molecular modeling approach

BackgroundAllergic disorders, prevalent global health concerns, afflict a substantial portion of the world’s population. These maladies result from an exaggerated immune system response to ordinarily innocuous substances, such as pollen, dust mites, and specific dietary components. Clinical manifestations of this heightened immune response include itching, swelling, and respiratory impairment, often accompanied by releasing mediators like histamine. The pathophysiological mechanisms of allergy disorders are intricate, arising from a complex interplay between genetic and environmental factors. While clinical presentations may vary, all allergy conditions share a common foundation in the dysregulated immune response to allergens.ResultThe current aim of this study was to identify innovative anti-allergic agents capable of inhibiting histamine and effectively mitigating allergic reactions by utilizing the computer-aided drug design approach by discovery studio (DS) 2022 v 23.1.1 package. The overarching aim was identifying potential drug candidates targeting the active site within the histamine H1 receptor complex; therefore, a collection of 4000 small druggable compounds was curated from ZINC, PubChem, and DRUG BANK databases sources. Four compounds appeared as promising candidates after assessing docking scores and binding energies. Notably, Compound ID 34154, recognized as tymazoline, showed the highest affinity for the H1 receptor of 3RZE, suggesting it may be the most promising choice for more research. Further chemoinformatic and ADMET (absorption, distribution, metabolism, excretion, and toxicity) analyses were conducted to assess the drug-like qualities of this chosen molecule. In addition, bioisosteric substitution techniques were employed to enhance tymazoline’s ADMET characteristics.ConclusionTymazoline shows strong binding affinity with 3RZE and verified all the drug-likeness criteria to inhibit the allergic disorders. Furthermore, molecular dynamics (MD) studies corroborated tymazoline’s potential as an anti-allergic agent, demonstrating contact between the ligand and the receptor that is well defined and stable.

Identifying p56lck SH2 Domain Inhibitors Using Molecular Docking and In Silico Scaffold Hopping

Bacterial infections are the second-leading cause of death, globally. The prevalence of antibacterial resistance has kept the demand strong for the development of new and potent drug candidates. It has been demonstrated that Src protein tyrosine kinases (TKs) play an important role in the regulation of inflammatory responses to tissue injury, which can trigger the onset of several severe diseases. We carried out a search for novel Src protein TK inhibitors, commencing from reported highly potent anti-bacterial compounds obtained using the Mannich reaction, using a combination of e-pharmacophore modeling, virtual screening, ensemble docking, and core hopping. The top-scoring compounds from ligand-based virtual screening were modified using protein structure-based design approaches, and their binding to the Src homology-2 domain of p56lck TK was predicted using ensemble molecular docking. We have prepared a database of 202 small molecules and have identified six novel top hits that can be subjected to further investigation. We have also performed in silico ADMET property prediction for the hit compounds. This combined computer-aided drug design approach can serve as a starting point for identifying novel TK inhibitors that could be further subjected to in vitro studies and validation of antimicrobial activity.

In Silico Design of Potential Small-Molecule Antibiotic Adjuvants against Salmonella typhimurium Ortho Acetyl Sulphydrylase Synthase to Address Antimicrobial Resistance.

The inhibition of O-acetyl sulphydrylase synthase isoforms has been reported to represent a promising approach for the development of antibiotic adjuvants. This occurs via the organism developing an unpaired oxidative stress response, causing a reduction in antibiotic resistance in vegetative and swarm cell populations. This consequently increases the effectiveness of conventional antibiotics at lower doses. This study aimed to predict potential inhibitors of Salmonella typhimurium ortho acetyl sulphydrylase synthase (StOASS), which has lower binding energy than the cocrystalized ligand pyridoxal 5 phosphate (PLP), using a computer-aided drug design approach including pharmacophore modeling, virtual screening, and in silico ADMET (Absorption, Distribution, Metabolism, Excretion, and Toxicity) evaluation. The screening and molecular docking of 4254 compounds obtained from the PubChem database were carried out using AutoDock vina, while a post-screening analysis was carried out using Discovery Studio. The best three hits were compounds with the PubChem IDs 118614633, 135715279, and 155773276, possessing binding affinities of -9.1, -8.9, and -8.8 kcal/mol, respectively. The in silico ADMET prediction showed that the pharmacokinetic properties of the best hits were relatively good. The optimization of the best three hits via scaffold hopping gave rise to 187 compounds, and they were docked against StOASS; this revealed that lead compound 1 had the lowest binding energy (-9.3 kcal/mol) and performed better than its parent compound 155773276. Lead compound 1, with the best binding affinity, has a hydroxyl group in its structure and a change in the core heterocycle of its parent compound to benzimidazole, and pyrimidine introduces a synergistic effect and consequently increases the binding energy. The stability of the best hit and optimized compound at the StOASS active site was determined using RMSD, RMSF, radius of gyration, and SASA plots generated from a molecular dynamics simulation. The MD simulation results were also used to monitor how the introduction of new functional groups of optimized compounds contributes to the stability of ligands at the target active site. The improved binding affinity of these compounds compared to PLP and their toxicity profile, which is predicted to be mild, highlights them as good inhibitors of StOASS, and hence, possible antimicrobial adjuvants.

Exploring Lead-Like Molecules of Traditional Chinese Medicine for Treatment Quest against Aliarcobacter butzleri: In Silico Toxicity Assessment, Dynamics Simulation, and Pharmacokinetic Profiling.

Aliarcobacter butzleri is a Gram-negative, curved or spiral-shaped, microaerophilic bacterium and causes human infections, specifically diarrhea, fever, and sepsis. The research objective of this study was to employ computer-aided drug design techniques to identify potential natural product inhibitors of a vital enzyme in this bacterium. The pyrimidine biosynthesis pathway in its core genome fraction is crucial for its survival and presents a potential target for novel therapeutics. Hence, novel small molecule inhibitors were identified (from traditional Chinese medicinal (TCM) compound library) against it, which may be used for possible curbing of infection by A. butzleri. Methods. A comprehensive subtractive genomics approach was utilized to identify a key enzyme (orotidine-5'-phosphate decarboxylase) cluster conserved in the core genome fraction of A. butzleri. It was selected for inhibitor screening due to its vital role in pyrimidine biosynthesis. TCM library (n > 36,000 compounds) was screened against it using pharmacophore model based on orotidylic acid (control), and the obtained lead-like molecules were subjected to structural docking using AutoDock Vina. The top-scoring compounds, ZINC70454134, ZINC85632684, and ZINC85632721, underwent further scrutiny via a combination of physiological-based pharmacokinetics, toxicity assessment, and atomic-scale dynamics simulations (100 ns). Among the screened compounds, ZINC70454134 displayed the most favorable characteristics in terms of binding, stability, absorption, and safety parameters. Overall, traditional Chinese medicine (TCM) compounds exhibited high bioavailability, but in diseased states (cirrhosis, renal impairment, and steatosis), there was a significant decrease in absorption, Cmax, and AUC of the compounds compared to the healthy state. Furthermore, MD simulation demonstrated that the ODCase-ZINC70454134 complex had a superior overall binding affinity, supported by PCA proportion of variance and eigenvalue rank analysis. These favorable characteristics underscore its potential as a promising drug candidate. The computer-aided drug design approach employed for this study helped expedite the discovery of antibacterial compounds against A. butzleri, offering a cost-effective and efficient approach to address infection by it. It is recommended that ZINC70454134 should be considered for further experimental analysis due to its indication as a potential therapeutic agent for combating A. butzleri infections. This study provides valuable insights into the molecular basis of biophysical inhibition of A. butzleri through TCM compounds.

Targeting with Structural Analogs of Natural Products the Purine Salvage Pathway in Leishmania (Leishmania) infantum by Computer-Aided Drug-Design Approaches.

Visceral Leishmaniasis (VL) has a high death rate, with 500,000 new cases and 50,000 deaths occurring annually. Despite the development of novel strategies and technologies, there is no adequate treatment for the disease. Therefore, the purpose of this study is to find structural analogs of natural products as potential novel drugs to treat VL. We selected structural analogs from natural products that have shown antileishmanial activities, and that may impede the purine salvage pathway using computer-aided drug-design (CADD) approaches. For these, we started with the vastly studied target in the pathway, the adenine phosphoribosyl transferase (APRT) protein, which alone is non-essential for the survival of the parasite. Keeping this in mind, we search for a substance that can bind to multiple targets throughout the pathway. Computational techniques were used to study the purine salvage pathway from Leishmania infantum, and molecular dynamic simulations were used to gather information on the interactions between ligands and proteins. Because of its low homology to human proteins and its essential role in the purine salvage pathway proteins network interaction, the findings further highlight the significance of adenylosuccinate lyase protein (ADL) as a therapeutic target. An analog of the alkaloid Skimmianine, N,N-diethyl-4-methoxy-1-benzofuran-6-carboxamide, demonstrated a good binding affinity to APRT and ADL targets, no expected toxicity, and potential for oral route administration. This study indicates that the compound may have antileishmanial activity, which was granted in vitro and in vivo experiments to settle this finding in the future.

Design of new small molecules derived from indolin-2-one as potent TRKs inhibitors using a computer-aided drug design approach

Tropomyosin receptor kinase (TRKs) enzymes are responsible for cancers associated with the neurotrophic tyrosine kinase receptor gene fusion and are identified as effective targets for anticancer drug discovery. A series of small-molecule indolin-2-one derivatives showed remarkable biological activity against TRKs enzymatic activity. These small molecules could have an excellent profile for pharmaceutical application in the treatment of cancers caused by TRKs activity. The aim of this study is to modify the structure of these molecules to obtain new molecules with improved TRK inhibitory activity and pharmacokinetic properties favorable to the design of new drugs. Based on these series, we carried out a 3D-QSAR study. As a result, robust and reliable CoMFA and CoMSIA models are developed and applied to the design of 11 new molecules. These new molecules have a biological activity superior to the most active molecule in the starting series. The eleven designed molecules are screened using drug-likeness, ADMET proprieties, molecular docking, and MM-GBSA filters. The results of this screening identified the T1, T3, and T4 molecules as the best candidates for strong inhibition of TRKs enzymatic activity. In addition, molecular dynamics simulations are performed for TRK free and complexed with ligands T1, T3, and T4 to evaluate the stability of ligand-protein complexes over the simulation time. On the other hand, we proposed experimental synthesis routes for these newly designed molecules. Finally, the designed molecules T1, T2, and T3 have great potential to become reliable candidates for the conception of new drug inhibitors of TRKs. Communicated by Ramaswamy H. Sarma

Targeting Leishmania infantum Mannosyl-oligosaccharide glucosidase with natural products: potential pH-dependent inhibition explored through computer-aided drug design.

Visceral Leishmaniasis (VL) is a serious public health issue, documented in more than ninety countries, where an estimated 500,000 new cases emerge each year. Regardless of novel methodologies, advancements, and experimental interventions, therapeutic limitations, and drug resistance are still challenging. For this reason, based on previous research, we screened natural products (NP) from Nuclei of Bioassays, Ecophysiology, and Biosynthesis of Natural Products Database (NuBBEDB), Mexican Compound Database of Natural Products (BIOFACQUIM), and Peruvian Natural Products Database (PeruNPDB) databases, in addition to structural analogs of Miglitol and Acarbose, which have been suggested as treatments for VL and have shown encouraging action against parasite's N-glycan biosynthesis. Using computer-aided drug design (CADD) approaches, the potential inhibitory effect of these NP candidates was evaluated by inhibiting the Mannosyl-oligosaccharide Glucosidase Protein (MOGS) from Leishmania infantum, an enzyme essential for the protein glycosylation process, at various pH to mimic the parasite's changing environment. Also, computational analysis was used to evaluate the Absorption, Distribution, Metabolism, Excretion, and Toxicity (ADMET) profile, while molecular dynamic simulations were used to gather information on the interactions between these ligands and the protein target. Our findings indicated that Ocotillone and Subsessiline have potential antileishmanial effects at pH 5 and 7, respectively, due to their high binding affinity to MOGS and interactions in the active center. Furthermore, these compounds were non-toxic and had the potential to be administered orally. This research indicates the promising anti-leishmanial activity of Ocotillone and Subsessiline, suggesting further validation through in vitro and in vivo experiments.

Insilico exploration C. koseri ATP synthase inhibitors by pharmacophore-based virtual screening, molecular docking and MD simulation.

Citrobacter koseri is a gram-negative rod that causes infections in people who have significant comorbidities and are immunocompromised. Antibiotic-resistant strains are becoming more common, which complicates infection treatment and highlights the need for innovative, effective drugs to fight these resistant strains. The enzyme complex ATP synthase participates in the adenosine triphosphate (ATP) synthesis, the fundamental energy currency of cells. This study used Computer-Aided Drug Design approaches to identify potential inhibitors of C. koseri ATP synthase. SWISS-MODEL was used to predict the 3D structure of C. koseri ATP synthase. A ligand-based pharmacophore model was developed using chemical features of ampicillin. Following ligand-based virtual screening across nine databases, the 2043 screened hits were docked to the ATP synthase active site using the standard precision mode of the glide tool. Based on their binding affinities, the top ten compounds were selected for additional investigation. The binding affinities of the chosen compounds ranged from -10.021 to -8.452 kcal/mol. The top four compounds (PubChem-25230613, PubChem-74936833, CHEMBL263035, PubChem-44208924) with the best ADMET characteristics and binding modes were chosen. Thus, the feasible binding mechanisms of the selected compounds were subjected to stability analysis using the MD Simulation study, which revealed the compounds' stability as potent inhibitors within the protein binding pocket. This computational approach provides important insights into the rational design of novel therapeutics and emphasizes the importance of targeting essential metabolic pathways when combating antibiotic-resistant pathogens. Future experimental validation and optimization of the identified inhibitors is required to determine their efficacy and safety profiles for clinical use.

Dual Targeting of DNA Damage Response Proteins Implicated in Cancer Radioresistance.

Ionizing radiation can induce different types of DNA lesions, leading to genomic instability and ultimately cell death. Radiation therapy or radiotherapy, a major modality in cancer treatment, harnesses the genotoxic potential of radiation to target and destroy cancer cells. Nevertheless, cancer cells have the capacity to develop resistance to radiation treatment (radioresistance), which poses a major obstacle in the effective management of cancer. It has been shown that administration of platinum-based drugs to cancer patients can increase tumor radiosensitivity, but despite this, it is associated with severe adverse effects. Several lines of evidence support that activation of the DNA damage response and repair machinery in the irradiated cancer cells enhances radioresistance and cellular survival through the efficient repair of DNA lesions. Therefore, targeting of key DNA damage repair factors would render cancer cells vulnerable to the irradiation effects, increase cancer cell killing, and reduce the risk of side effects on healthy tissue. Herein, we have employed a computer-aided drug design approach for generating ab initio a chemical compound with drug-like properties potentially targeting two proteins implicated in multiple DNA repair pathways. The findings of this study could be taken into consideration in clinical decision-making in terms of co-administering radiation with DNA damage repair factor-based drugs.

Identification of potential inhibitors targeting Ebola virus VP35 protein: a computational strategy

Ebola virus (EBOV) poses a severe threat as a highly infectious pathogen, causing devastating hemorrhagic fever in both humans and animals. The EBOV virus VP35 protein plays a crucial role in viral replication and exhibits the ability to suppress the host interferon cascade, leading to immune system depletion. As a potential drug target, VP35 protein inhibition holds promise for combating EBOV. To discover new drug candidates, we employed a computer-aided drug design approach, focusing on compounds capable of inhibiting VP35 protein replication. In this connection, a pharmacophore model was generated using molecular interactions between the VP35 protein and its inhibitor. ZINC and Cambridge database were screened using validated pharmacophore model. Further the compounds were filtered based on Lipinski’s rule of five and subjected to MD simulation and relative binding free energy calculation. Six compounds manifest a significant docking score and strong binding interaction towards VP35 protein. MD simulations further confirmed the remarkable stability of these six complexes. Relative binding free energy calculations also showed significant ΔG value in the range of −132.3 and −49.3 kcal/mol. This study paves the way for further optimization of these compounds as potential inhibitors of VP35, facilitating subsequent experimental in vitro studies. Communicated by Ramaswamy H. Sarma

Antibreast cancer activities of phytochemicals from Anonna muricata using computer-aided drug design (CADD) approach

Abstract Antibreast cancer activities of 131 phytochemicals from Annona muricata (Soursop) were investigated against human placental aromatase (PDB ID: 3S7S), a prominent target receptor in breast cancer therapy using computer aided-drug design approach. An antibreast cancer drug (tamoxifen) was used for comparison. The result of this work flourishes caffeoquinic acid (−8.4 kcal/mol), roseoside (−8.3 kcal/mol), chlorogenic acid (−8.2 kcal/mol), feruloylglycoside (−8.1 kcal/mol), citroside A (−8.0 kcal/mol), and coreximine (−7.8 kcal/mol), as probable inhibitors of human placental aromatase. This is due to their excellent binding affinities (ΔG), coupled with outstanding druglike, absorption, distribution, metabolism, excretion, and toxicity profiles, bioavailability and oral-bioavailability properties, and the interactions of important residues with the active pocket of human placental aromatase. All the results obtained were similar to that of the standards tamoxifen (−8.0 kcal/mol) but could be better when optimized. Thus, lead optimization, molecular dynamics, and in vivo investigations are thereby recommended on the identified potent compounds in the quest of developing new therapeutic agents against breast cancer.

An Exploration of the Inhibitory Mechanism of Rationally Screened Benzofuran-1,3,4-Oxadiazoles and-1,2,4-Triazoles as Inhibitors of NS5B RdRp Hepatitis C Virus through Pharmacoinformatic Approaches.

Benzofuran, 1,3,4-oxadiazole, and 1,2,4-triazole are privileged heterocyclic moieties that display the most promising and wide spectrum of biological activities against a wide variety of diseases. In the current study, benzofuran-1,3,4-oxadiazole BF1-BF7 and benzofuran-1,2,4-triazole compounds BF8-BF15 were tested against HCV NS5B RNA-dependent RNA polymerase (RdRp) utilizing structure-based screening via a computer-aided drug design (CADD) approach. A molecular docking approach was applied to evaluate the binding potential of benzofuran-appended 1,3,4-oxadiazole and 1,2,4-triazole BF1-BF15 molecules. Benzofuran-1,3,4-oxadiazole scaffolds BF1-BF7 showed lesser binding affinities (-12.63 to -14.04 Kcal/mol) than benzofuran-1,2,4-triazole scaffolds BF8-BF15 (-14.11 to -16.09 Kcal/mol) against the HCV NS5B enzyme. Molecular docking studies revealed the excellent binding affinity scores exhibited by benzofuran-1,2,4-triazole structural motifs BF-9 (-16.09 Kcal/mol), BF-12 (-15.75 Kcal/mol), and BF-13 (-15.82 Kcal/mol), respectively, which were comparatively better than benzofuran-based HCV NS5B inhibitors' standard reference drug Nesbuvir (-15.42 Kcal/mol). A molecular dynamics simulation assay was also conducted to obtain valuable insights about the enzyme-compounds interaction profile and structural stability, which indicated the strong intermolecular energies of the BF-9+NS5B complex and the BF-12+NS5B complex as per the MM-PBSA method, while the BF-12+NS5B complex was the most stable system as per the MM-GBSA calculation. The drug-likeness and ADMET studies of all the benzofuran-1,2,4-triazole derivatives BF8-BF15 revealed that these compounds possessed good medicinal chemistry profiles in agreement with all the evaluated parameters for being drugs. The molecular docking affinity scores, MM-PBSA/MM-GBSA and MD-simulation stability analysis, drug-likeness profiling, and ADMET study assessment indicated that N-4-fluorophenyl-S-linked benzofuran-1,2,4-triazole BF-12 could be a future promising anti-HCV NS5B RdRp inhibitor therapeutic drug candidate that has a structural agreement with the Nesbuvir standard reference drug.

Design, Synthesis, Evaluation and Molecular Dynamics Simulation of Dengue Virus NS5-RdRp Inhibitors.

Dengue virus (DENV) is a major mosquito-borne human pathogen in tropical countries; however, there are currently no targeted antiviral treatments for DENV infection. Compounds 27 and 29 have been reported to be allosteric inhibitors of DENV RdRp with potent inhibitory effects. In this study, the structures of compounds 27 and 29 were optimized using computer-aided drug design (CADD) approaches. Nine novel compounds were synthesized based on rational considerations, including molecular docking scores, free energy of binding to receptor proteins, predicted Absorption, Distribution, Metabolism, Excretion, and Toxicity (ADMET) parameters, structural diversity, and feasibility of synthesis. Subsequently, the anti-DENV activity was assessed. In the cytopathic effect (CPE) assay conducted on BHK-21 cells using the DENV2 NGC strain, both SW-b and SW-d demonstrated comparable or superior activity against DENV2, with IC50 values of 3.58 ± 0.29 μM and 23.94 ± 1.00 μM, respectively, compared to that of compound 27 (IC50 = 19.67 ± 1.12 μM). Importantly, both SW-b and SW-d exhibited low cytotoxicity, with CC50 values of 24.65 μmol and 133.70 μmol, respectively, resulting in selectivity indices of 6.89 and 5.58, respectively. Furthermore, when compared to the positive control compound 3'-dATP (IC50 = 30.09 ± 8.26 μM), SW-b and SW-d displayed superior inhibitory activity in an enzyme inhibitory assay, with IC50 values of 11.54 ± 1.30 μM and 13.54 ± 0.32 μM, respectively. Molecular dynamics (MD) simulations elucidated the mode of action of SW-b and SW-d, highlighting their ability to enhance π-π packing interactions between benzene rings and residue W795 in the S1 fragment, compared to compounds 27 and 29. Although the transacylsulphonamide fragment reduced the interaction between T794 and NH, it augmented the interaction between R729 and T794. In summary, our study underscores the potential of SW-b and SW-d as allosteric inhibitors targeting the DENV NS5 RdRp domain. However, further in vivo studies are warranted to assess their pharmacology and toxicity profiles.

An Updated Review on Developing Small Molecule Kinase Inhibitors Using Computer-Aided Drug Design Approaches.

Small molecule kinase inhibitors (SMKIs) are of heightened interest in the field of drug research and development. There are 79 (as of July 2023) small molecule kinase inhibitors that have been approved by the FDA and hundreds of kinase inhibitor candidates in clinical trials that have shed light on the treatment of some major diseases. As an important strategy in drug design, computer-aided drug design (CADD) plays an indispensable role in the discovery of SMKIs. CADD methods such as docking, molecular dynamic, quantum mechanics/molecular mechanics, pharmacophore, virtual screening, and quantitative structure-activity relationship have been applied to the design and optimization of small molecule kinase inhibitors. In this review, we provide an overview of recent advances in CADD and SMKIs and the application of CADD in the discovery of SMKIs.

The importance of good practices and false hits for QSAR-driven virtual screening real application: a SARS-CoV-2 main protease (Mpro) case study

Computer-Aided Drug Design (CADD) approaches, such as those employing quantitative structure-activity relationship (QSAR) methods, are known for their ability to uncover novel data from large databases. These approaches can help alleviate the lack of biological and chemical data, but some predictions do not generate sufficient positive information to be useful for biological screenings. QSAR models are often employed to explain biological data of chemicals and to design new chemicals based on their predictions. In this review, we discuss the importance of data set size with a focus on false hits for QSAR approaches. We assess the challenges and reliability of an initial in silico strategy for the virtual screening of bioactive molecules. Lastly, we present a case study reporting a combination approach of hologram-based quantitative structure-activity relationship (HQSAR) models and random forest-based QSAR (RF-QSAR), based on the 3D structures of 25 synthetic SARS-CoV-2 Mpro inhibitors, to virtually screen new compounds for potential inhibitors of enzyme activity. In this study, optimal models were selected and employed to predict Mpro inhibitors from the database Brazilian Compound Library (BraCoLi). Twenty-four compounds were then assessed against SARS-CoV-2 Mpro at 10 µM. At the time of this study (March 2021), the availability of varied and different Mpro inhibitors that were reported definitely affected the reliability of our work. Since no hits were obtained, the data set size, parameters employed, external validations, as well as the applicability domain (AD) could be considered regarding false hits data contribution, aiming to enhance the design and discovery of new bioactive molecules.

Computer-aided drug design approaches applied to screen natural product’s structural analogs targeting arginase in Leishmania spp

Introduction: Leishmaniasis is a disease with high mortality rates and approximately 1.5 million new cases each year. Despite the new approaches and advances to fight the disease, there are no effective therapies. Methods: Hence, this study aims to screen for natural products' structural analogs as new drug candidates against leishmaniasis. We applied Computer-aided drug design (CADD) approaches, such as virtual screening, molecular docking, molecular dynamics simulation, molecular mechanics–generalized Born surface area (MM–GBSA) binding free estimation, and free energy perturbation (FEP) aiming to select structural analogs from natural products that have shown anti-leishmanial and anti-arginase activities and that could bind selectively against the Leishmania arginase enzyme. Results: The compounds 2H-1-benzopyran, 3,4-dihydro-2-(2-methylphenyl)-(9CI), echioidinin, and malvidin showed good results against arginase targets from three parasite species and negative results for potential toxicities. The echioidinin and malvidin ligands generated interactions in the active center at pH 2.0 conditions by MM-GBSA and FEP methods. Conclusions: This work suggests the potential anti-leishmanial activity of the compounds and thus can be further in vitro and in vivo experimentally validated.

Medicinal and toxicological investigation of some common NSAIDs; A computer-aided drug design approach

Nonsteroidal anti-inflammatory drugs (NSAIDs) are commonly used for the treatment of pain, fever, and inflammation. Their consumption rises some crucial side effects like gastrointestinal, cardiovascular, and renal injury because of the non-selective reduction of prostaglandin synthesis, with differences in the extent of inhibition of the cyclooxygenase. Attempts have been taken for a comparative biochemical, medicinal, and toxicological investigation to raise awareness of the consumption of NSAIDs. In-silico methods were incorporated to investigate their physicochemical, spectral, medicinal, pharmacological, and toxicological properties. Molecular docking and non-bonding calculations were performed against the human prostaglandin synthase protein to reveal their binding affinity and interactions with the amino acid residues of the receptor protein. Further, molecular dynamics simulation was conducted to validate the binding mode and drug-protein stability at the inhibiting binding pocket. The physical and chemical analysis supports their geometries and spectral results confirming the presence of essential functional groups at the core structure. From ADMET and PASS prediction, all the drugs are non-carcinogenic (except IBP), and nonselective NSAIDs showed relatively higher adverse effects compared to selective agents. Based on the above studies, this report can be helpful to understand more deeply the medicinal, and toxicological effects of the reported NSAIDs.