High-Performance Drug Discovery: Computational Screening by Combining Docking and Molecular Dynamics Simulations

Virtual compound screening using molecular docking is widely used in the discovery of new lead compounds for drug design. However, this method is not completely reliable and therefore unsatisfactory. In this study, we used massive molecular dynamics simulations of protein-ligand conformations obtained by molecular docking in order to improve the enrichment performance of molecular docking. Our screening approach employed the molecular mechanics/Poisson-Boltzmann and surface area method to estimate the binding free energies. For the top-ranking 1,000 compounds obtained by docking to a target protein, approximately 6,000 molecular dynamics simulations were performed using multiple docking poses in about a week. As a result, the enrichment performance of the top 100 compounds by our approach was improved by 1.6–4.0 times that of the enrichment performance of molecular dockings. This result indicates that the application of molecular dynamics simulations to virtual screening for lead discovery is both effective and practical. However, further optimization of the computational protocols is required for screening various target proteins.     

基于处方挖掘与分子动力学模拟筛选严重急性呼吸综合征冠状病毒2潜在抑制剂分子的研究

目的 从中药筛选具有潜在抑制严重急性呼吸综合征冠状病毒2 （SARS-CoV-2） 活性的成分，进一步从原子水平揭 示其抑制SARS-CoV-2 表面刺突蛋白（S 蛋白） 受体结合域（RBD） 与血管紧张素转化酶2 （ACE2） 结合的内在机制。 方法 检索新型冠状病毒（简称“新冠肺炎”） 治疗中药处方，构建“新冠肺炎中药候选活性成分数据库”。用具有ACE2 抑制活性的小分子化合物构建HipHop药效团模型，并对“新冠肺炎中药候选活性成分数据库”中活性成分筛选。采用分子 对接和分子动力学模拟方法研究候选活性成分与ACE2 的结合方式及其对SARS-CoV-2 S 蛋白与ACE2 识别的影响。 结果 本文通过中药处方挖掘和分子动力学模拟，从143 个新冠肺炎治疗中药处方中筛选出10 种可与SARS-CoV-2 S 蛋白/ 人源ACE2 识别位点结合的中药成分。其中，枇杷叶主要活性成分23-trans-p-coumaryhormentic acid 与ACE2 具有最高的亲和 力，且23-trans-p-coumaryhormentic acid 的结合可有效阻断SARS-CoV-2 S蛋白与宿主细胞ACE2 的结合。结论 本文通过虚 拟筛选发现了SARS-CoV-2 潜在抑制剂分子23-trans-p-coumaryhormentic acid，同时从原子水平预测了其抑制SARS-CoV-2 S 蛋白与ACE2 结合的内在机制，这将为SARS-CoV-2 特异性抗病毒药物的研发提供理论依据。     

Key Features for Designing Phosphodiesterase-5 Inhibitors

Abstract

Phosphodiesterase superfamily is the key regulator of 3′,5′-cyclic guanosine monophosphate (cGMP) decomposition in human body. Phosphodiesterase-5 (PDE-5) inhibitors, sildenafil, vardenafil and tadalafil, are well known oral treatment for males with erectile dysfunction. To investigate the inhibitory effects of traditional Chinese medicine (TCM) compounds to PDE-5, we performed both ligand-based and structure-based studies on this topic. Comparative molecular field analysis (CoMFA) and comparative molecular similarity indices analysis (CoMSIA) studies were conducted to construct three dimensional quantitative structure-activity relationship (3D-QSAR) models of series of known PDE-5 inhibitors. The predictive models had cross-validated, q², and non cross-validated coefficient, r², values of 0.791 and 0.948 for CoMFA and 0.724 and 0.908 for CoMSIA. These two 3D-QSAR models were used to predict activity of TCM compounds. Docking simulations were performed to further analyze the binding mode of training set and TCM compounds. A putative binding model was proposed based on CoMFA and CoMSIA contour maps and docking simulations; formation of pi-stacking, water bridge and specific hydrogen bonding were deemed important interactions between ligands and PDE-5. Of our TCM compounds, engeletin, satisfied our binding model, and hence, emerged as PDE-5 inhibitor candidate.

Using this study as an example, we demonstrated that docking should be conducted for qualitative purposes, such as identifying protein characteristics, rather than for quantitative analyses that rank compound efficacy based on results of scoring functions. Prediction of compound activity should be reserved for QSAR analyses, and scoring functions and docking scores should be used for preliminary screening of TCM database (http://tcm.cmu.edu.tw/index.php).     

Shape-based virtual screening,docking, and molecular dynamics simulations to identify Mtb-ASADH inhibitors

Aspartate β-semialdehyde dehydrogenase (ASADH) is a key enzyme for the biosynthesis of essential amino acids and several important metabolites in microbes. Inhibition of ASADH enzyme is a promising drug target strategy against Mycobacterium tuberculosis (Mtb). In this work, in silico approach was used to identify potent inhibitors of Mtb-ASADH. Aspartyl β-difluorophosphonate (β-AFP), a known lead compound, was used to understand the molecular recognition interactions (using molecular docking and molecular dynamics analysis). This analysis helped in validating the computational protocol and established the participation of Arg99, Glu224, Cys130, Arg249, and His256 amino acids as the key amino acids in stabilizing ligand–enzyme interactions for effective binding, an essential feature is H-bonding interactions with the two arginyl residues at the two ends of the ligand. Best binding conformation of β-AFP was selected as a template for shape-based virtual screening (ZINC and NCI databases) to identify compounds that competitively inhibit the Mtb-ASADH. The top rank hits were further subjected to ADME and toxicity filters. Final filter was based on molecular docking analysis. Each screened molecule carries the characteristics of the highly electronegative groups on both sides separated by an average distance of 6?Å. Finally, the best predicted 20 compounds exhibited minimum three H-bonding interactions with Arg99 and Arg249. These identified hits can be further used for designing the more potent inhibitors against ASADH family. MD simulations were also performed on two selected compounds (NSC4862 and ZINC02534243) for further validation. During the MD simulations, both compounds showed same H-bonding interactions and remained bound to key active residues of Mtb-ASADH.     

Key features for designing phosphodiesterase-5 inhibitors

Phosphodiesterase superfamily is the key regulator of 3',5'-cyclic guanosine monophosphate (cGMP) decomposition in human body. Phosphodiesterase-5 (PDE-5) inhibitors, sildenafil, vardenafil and tadalafil, are well known oral treatment for males with erectile dysfunction. To investigate the inhibitory effects of traditional Chinese medicine (TCM) compounds to PDE-5, we performed both ligand-based and structure-based studies on this topic. Comparative molecular field analysis (CoMFA) and comparative molecular similarity indices analysis (CoMSIA) studies were conducted to construct three dimensional quantitative structure-activity relationship (3D-QSAR) models of series of known PDE-5 inhibitors. The predictive models had cross-validated, q(2), and non cross-validated coefficient, r(2), values of 0.791 and 0.948 for CoMFA and 0.724 and 0.908 for CoMSIA. These two 3D-QSAR models were used to predict activity of TCM compounds. Docking simulations were performed to further analyze the binding mode of training set and TCM compounds. A putative binding model was proposed based on CoMFA and CoMSIA contour maps and docking simulations; formation of pi-stacking, water bridge and specific hydrogen bonding were deemed important interactions between ligands and PDE-5. Of our TCM compounds, engeletin, satisfied our binding model, and hence, emerged as PDE-5 inhibitor candidate. Using this study as an example, we demonstrated that docking should be conducted for qualitative purposes, such as identifying protein characteristics, rather than for quantitative analyses that rank compound efficacy based on results of scoring functions. Prediction of compound activity should be reserved for QSAR analyses, and scoring functions and docking scores should be used for preliminary screening of TCM database (http://tcm.cmu.edu.tw/index.php).     

Discovery of promising FtsZ inhibitors by E-pharmacophore, 3D-QSAR,molecular docking study,and molecular dynamics simulation

Filamentous temperature-sensitive protein Z (FtsZ), playing a key role in bacterial cell division, is regarded as a promising target for the design of antimicrobial agent. This study is looking for potential high-efficiency FtsZ inhibitors. Ligand-based pharmacophore and E-pharmacophore, virtual screening and molecular docking were used to detect promising FtsZ inhibitors, and molecular dynamics simulation was used to study the stability of protein-ligand complexes in this paper. Sixty-three inhibitors from published literatures with pIC₅₀ ranging from 2.483 to 5.678 were collected to develop ligand-based pharmacophore model. 4DXD bound with 9PC was selected to develop the E-pharmacophore model. The pharmacophore models validated by test set method and decoy set were employed for virtual screening to exclude inactive compounds against ZINC database. After molecular docking, ADME analysis, IFD docking and MM-GBSA, 8 hits were identified as potent FtsZ inhibitors. A 50?ns molecular dynamics simulation was implemented on the compounds to assess the stability between potent inhibitors and FtsZ. The results indicated that the candidate compounds had a high docking score and were strongly combined with FtsZ by forming hydrogen bonding interactions with key amino acid residues, and van der Waals forces and hydrophobic interactions had significant contribution to the stability of the binding. Molecular dynamics simulation results showed that the protein-ligand compounds performed well in both the stability and flexibility of the simulation process.     

Locking the 150-Cavity Open: In Silico Design and Verification of Influenza Neuraminidase Inhibitors

Neuraminidase (NA) of influenza is a key target for virus infection control and the recently discovered open 150-cavity in group-1 NA provides new opportunity for novel inhibitors design. In this study, we used a combination of theoretical methods including fragment docking, molecular linking and molecular dynamics simulations to design ligands that specifically target at the 150-cavity. Through in silico screening of a fragment compound library on the open 150-cavity of NA, a few best scored fragment compounds were selected to link with Zanamivir, one NA-targeting drug. The resultant new ligands may bind both the active site and the 150-cavity of NA simultaneously. Extensive molecular dynamics simulations in explicit solvent were applied to validate the binding between NA and the designed ligands. Moreover, two control systems, a positive control using Zanamivir and a negative control using a low-affinity ligand 3-(p-tolyl) allyl-Neu5Ac2en (ETT, abbreviation reported in the PDB) found in a recent experimental work, were employed to calibrate the simulation method. During the simulations, ETT was observed to detach from NA, on the contrary, both Zanamivir and our designed ligand bind NA firmly. Our study provides a prospective way to design novel inhibitors for controlling the spread of influenza virus.     

Discovery of benzimidazole derivatives as novel multi-target EGFR, VEGFR-2 and PDGFR kinase inhibitors

Multi-target EGFR, VEGFR-2 and PDGFR inhibitors are highly useful anticancer agents with improved therapeutic efficacies. In this work, we used two virtual screening methods, support vector machines (SVM) and molecular docking, to identify a novel series of benzimidazole derivatives, 2-aryl benzimidazole compounds, as multi-target EGFR, VEGFR-2 and PDGFR inhibitors. 2-Aryl benzimidazole compounds were synthesized and their biological activities against a tumor cell line HepG-2 and specific kinases were evaluated. Among these compounds, compounds 5a and 5e exhibited high cytotoxicity against HepG-2 cells with IC?? values at ～2 μM. Further kinase assay study showed that compound 5a have good EGFR inhibitory activity and moderate VEGFR-2 and PDGFR inhibitory activities, while 5e have moderate EGFR inhibitory activity and slightly weaker VEGFR-2 and PDGFR inhibitory activities. Molecular docking analysis suggested that compound 5a more tightly interacts with EGFR and PDGFR than compound 5e. Our study discovered a novel series of benzimidazole derivatives as multi-target EGFR, VEGFR-2 and PDGFR kinases inhibitors.     

Structure-based virtual screening of novel tubulin inhibitors and their characterization as anti-mitotic agents

Microtubule cytoskeletons are involved in many essential functions throughout the life cycle of cells, including transport of materials into cells, cell movement, and proper progression of cell division. Small compounds that can bind at the colchicine site of tubulin have drawn great attention because these agents can suppress or inhibit microtubule dynamics and tubulin polymerization. To find novel tubulin polymerization inhibitors as anti-mitotic agents, we performed a virtual screening study of the colchicine binding site on tubulin. Novel tubulin inhibitors were identified and characterized by their inhibitory activities on tubulin polymerization in vitro. The structural basis for the interaction of novel inhibitors with tubulin was investigated by molecular modeling, and we have proposed binding models for these hit compounds with tubulin. The proposed docking models were very similar to the binding pattern of colchicine or podophyllotoxin with tubulin. These new hit compound derivatives exerted growth inhibitory effects on the HL60 cell lines tested and exhibited strong cell cycle arrest at G2/M phase. Furthermore, these compounds induced apoptosis after cell cycle arrest. In this study, we show that the validated derivatives of compound 11 could serve as potent lead compounds for designing novel anti-cancer agents that target microtubules.     

Discovery of chromenes as inhibitors of macrophage migration inhibitory factor

Macrophage migration inhibitory factor (MIF) is an essential signaling cytokine with a key role in the immune system. Binding of MIF to its molecular targets such as, among others, the cluster of differentiation 74 (CD74) receptor plays a key role in inflammatory diseases and cancer. Therefore, the identification of MIF binding compounds gained importance in drug discovery. In this study, we aimed to discover novel MIF binding compounds by screening of a focused compound collection for inhibition of its tautomerase enzyme activity. Inspired by the known chromen-4-one inhibitor Orita-13, a focused collection of compounds with a chromene scaffold was screened for MIF binding. The library was synthesized using versatile cyanoacetamide chemistry to provide diversely substituted chromenes. The screening provided inhibitors with IC₅₀’s in the low micromolar range. Kinetic evaluation suggested that the inhibitors were reversible and did not bind in the binding pocket of the substrate. Thus, we discovered novel inhibitors of the MIF tautomerase activity, which may ultimately support the development of novel therapeutic agents against diseases in which MIF is involved.     

De novo design of VEGFR-2 tyrosine kinase inhibitors based on a linked-fragment approach

Vascular endothelial growth factor receptor-2 (VEGFR-2) tyrosine kinase inhibitors have been demonstrated to possess substantial antitumor activity. VEGFR-2 tyrosine kinase inhibitors are crucial for development of antitumor drugs. Based on the crystal structure of VEGFR-2 tyrosine kinase, a linked-fragment strategy was employed to design novel VEGFR-2 tyrosine kinase inhibitors, and 1000 compounds were generated in this process. Absorption, distribution, metabolism, excretion and toxicity (ADMET) were used to screen the 1000 compounds, and 59 compounds were acceptable. Scaffold hopping was then used for further screening, and only four compounds were obtained in this way. Then, the binding energy of the four molecules to VEGFR-2 tyrosine kinase was calculated using molecular docking, and their values were found to be lower than that of Sorafenib. Finally, molecular dynamics simulations were performed on the complex of the compound with the lowest binding energy with VEGFR-2 tyrosine kinase, and the binding model was analyzed. At the end, four chemical entities with novel structures were obtained, and were suggested for experimental testing in future studies.     

Identification of novel scaffolds to inhibit human mitotic kinesin Eg5 targeting the second allosteric binding site using in silico methods

Human mitotic kinesins are potential anticancer drug targets because of their essential role in mitotic cell division. The kinesin Eg5 (Kinesin-5, kif11) has gained much attention in this regard and has many inhibitors in different phases of clinical trials. All drug candidates considered for Eg5 so far binds to the binding site (Site 1) formed by the loop L5, helices α2 and α3 and are uncompetitive to ATP/ADP. Recently, it has been reported that Eg5 also has a second binding site (Site 2) formed by helices α4 and α6. In the current work, we have screened the compounds in the diversity set-III from National Cancer Institute (NCI) and Zinc database to identify potential inhibitors for Eg5 that specifically binds to the site 2. The compounds were ranked based on the glide extra precision docking scores and the top ranked compounds were found to have pyridazine scaffold. The top five compounds were further evaluated for other drug like properties. Stability of protein-ligand complexes were analyzed using molecular dynamic simulations. Our studies suggest that pyridazine analogs have good MDCK, permeability properties and high binding affinity to the human Eg5.     

Modeling of human CCR5 as Target for HIV-I and Virtual Screening with Marine Therapeutic compounds

Marine derivatives are of great pharmaceutical interest as inhibitory compound and search of bioactive compounds from Marine organism which is relatively new to medicinal chemistry. Our main aim in the study is to screen possible inhibitors against CCR5 which acts as co-receptor M-tropic HIV-1, through virtual screening of 122 Marine derived compounds from various organisms known to have biological activity. Homology Model of CCR5 was constructed using MODELLER and the Model was energy minimized and validated using PROCHECK to obtain a stable structure, which was further used for virtual screening of Marine derived compounds through molecular Docking studies using GOLD. The Docked complexes were validated and Enumerated based on the GOLD Scoring function to pick out the best Marine inhibitor based on GOLD score. Thus from the entire 122 Marine compounds which were Docked, we got best 4 of them with optimal GOLD Score. (LAMIVUDINE: 45.0218, BATZELLINE-D: 44.3852.ACYCLOVIR: 43.1362 and THIIOACETAMIDE: 42.7412) Further the Complexes were analyzed through LIGPLOT for their interaction for the 4 best docked Marine compounds. Thus from the Complex scoring and binding ability its deciphered that these Marine compounds could be promising inhibitors for M-tropic HIV-1 using CCR5 as Drug target yet pharmacological studies have to confirm it.     

Integrating docking and molecular dynamics approaches for a series of proline-based 2,5-diketopiperazines as novel αβ-tubulin inhibitors

In this work, docking tools were utilized in order to study the binding properties of more than five hundred of proline-based 2,5-diketopiperazine in the binding site of αβ-tubulin. Results revealed that 20 compounds among them showed lower binding energies in comparison with Tryprostatin-A, a well known tubulin inhibitor and therefore could be potential inhibitors of tubulin. However, the precise evaluation of binding poses represents the similar binding modes for all of these compounds and Tryprostatin-A. Finally, the best docked complex was subjected to a 25 ns molecular dynamics simulation to further validate the proposed binding mode of this compound.     

In silico multi-filter screening approaches for developing novel beta-secretase inhibitors

A large database of chemical structures was screened for potential inhibitors of β-secretase was carried out using in silico multi-filter techniques. Substructure screening, computer-aided ligand docking, binding free energy calculations, and partial interaction energy analyses were performed successively to identify chemical compounds which could serve as different scaffolds from known β-secretase inhibitors for future drug design. We showed that our in silico multi-filter screening retrieved all known inhibitors from the compound database investigated, which suggests that the other compounds identified as inhibitors by this computerized screening process are potential β-secretase inhibitors.     

Identification of novel CDK2 inhibitors by a multistage virtual screening method based on SVM,pharmacophore and docking model

Abstract

Cyclin-dependent kinase 2 (CDK2) is the family of Ser/Thr protein kinases that has emerged as a highly selective with low toxic cancer therapy target. A multistage virtual screening method combined by SVM, protein-ligand interaction fingerprints (PLIF) pharmacophore and docking was utilised for screening the CDK2 inhibitors. The evaluation of the validation set indicated that this method can be used to screen large chemical databases because it has a high hit-rate and enrichment factor (80.1% and 332.83 respectively). Six compounds were screened out from NCI, Enamine and Pubchem database. After molecular dynamics and binding free energy calculation, two compounds had great potential as novel CDK2 inhibitors and they also showed selective inhibition against CDK2 in the kinase activity assay.     

Novel virtual lead identification in the discovery of hematopoietic cell kinase (HCK) inhibitors: application of 3D QSAR and molecular dynamics simulation

High level of hematopoietic cell kinase (Hck) is associated with drug resistance in chronic myeloid leukemia. Additionally, Hck activity has also been connected with the pathogenesis of HIV-1 and chronic obstructive pulmonary disease. In this study, three-dimensional (3D) QSAR pharmacophore models were generated for Hck based on experimentally known inhibitors. A best pharmacophore model, Hypo1, was developed with high correlation coefficient (0.975), Low RMS deviation (0.60) and large cost difference (49.31), containing three ring aromatic and one hydrophobic aliphatic feature. It was further validated by the test set (r?=?0.96) and Fisher’s randomization method (95%). Hypo 1 was used as a 3D query for screening the chemical databases, and the hits were further screened by applying Lipinski’s rule of five and ADMET properties. Selected hit compounds were subjected to molecular docking to identify binding conformations in the active site. Finally, the appropriate binding modes of final hit compounds were revealed by molecular dynamics (MD) simulations and free energy calculation studies. Hence, we propose the final three hit compounds as virtual candidates for Hck inhibitors.     

Exploration of Virtual Candidates for Human HMG-CoA Reductase Inhibitors Using Pharmacophore Modeling and Molecular Dynamics Simulations

3-hydroxy-3-methylglutaryl coenzyme A reductase (HMGR) is a rate-controlling enzyme in the mevalonate pathway which involved in biosynthesis of cholesterol and other isoprenoids. This enzyme catalyzes the conversion of HMG-CoA to mevalonate and is regarded as a drug target to treat hypercholesterolemia. In this study, ten qualitative pharmacophore models were generated based on chemical features in active inhibitors of HMGR. The generated models were validated using a test set. In a validation process, the best hypothesis was selected based on the statistical parameters and used for virtual screening of chemical databases to find novel lead candidates. The screened compounds were sorted by applying drug-like properties. The compounds that satisfied all drug-like properties were used for molecular docking study to identify their binding conformations at active site of HMGR. The final hit compounds were selected based on docking score and binding orientation. The HMGR structures in complex with the hit compounds were subjected to 10 ns molecular dynamics simulations to refine the binding orientation as well as to check the stability of the hits. After simulation, binding modes including hydrogen bonding patterns and molecular interactions with the active site residues were analyzed. In conclusion, four hit compounds with new structural scaffold were suggested as novel and potent HMGR inhibitors.     

Identification of potential PKC inhibitors through pharmacophore designing, 3D-QSAR and molecular dynamics simulations targeting Alzheimer’s disease

Protein kinases are ubiquitously expressed as Serine/Threonine kinases, and play a crucial role in cellular activities. Protein kinases have evolved through stringent regulation mechanisms. Protein kinases are also involved in tauopathy, thus are important targets for developing Anti-Alzheimer’s disease compounds. Structures with an indole scaffold turned out to be potent new leads. With the aim of developing new inhibitors for human protein kinase C, here we report the generation of four point 3D geometric featured pharmacophore model. In order to identify novel and potent PKCθ inhibitors, the pharmacophore model was screened against 80,000,00 compounds from various chemical databases such as., ZINC, SPEC, ASINEX, which resulted in 127 compound hits, and were taken for molecular docking filters (HTVS, XP docking). After in-depth analysis of binding patterns, induced fit docking (flexible) was employed for six compounds along with the cocrystallized inhibitor. Molecular docking study reveals that compound 6F found to be tight binder at the active site of PKCθ as compared to the cocrystal and has occupancy of 90 percentile. MM-GBSA also confirmed the potency of the compound 6F as better than cocrystal. Molecular dynamics results suggest that compound 6F showed good binding stability of active sites residues similar to cocrystal 7G compound. Present study corroborates the pharmacophore-based virtual screening, and finds the compound 6F as a potent Inhibitor of PKC, having therapeutic potential for Alzheimer’s disease. Worldwide, 46.8 million people are believed to be living with Alzheimer’s disease. When elderly population increases rapidly and neurodegenerative burden also increases in parallel, we project the findings from this study will be useful for drug developing efforts targeting Alzheimer’s disease.     

The identification of novel,high affinity AQP9 inhibitors in an intracellular binding site

Abstract

Background: The involvement of aquaporin (AQP) water and small solute channels in the etiology of several diseases, including cancer, neuromyelitis optica and body fluid imbalance disorders, has been suggested previously. Furthermore, results obtained in a mouse model suggested that AQP9 function contributes to hyperglycemia in type-2 diabetes. In addition, the physiological role of several AQP family members remains poorly understood. Small molecule inhibitors of AQPs are therefore desirable to further study AQP physiological and pathophysiological functions. Methods: The binding of recently established AQP9 inhibitors to a homology model of AQP9 was investigated by molecular dynamics simulations and molecular docking. Putative inhibitor binding sites identified with this procedure were modified by site-directed mutagenesis. Active compounds were measured in a mammalian cell water permeability assay of mutated AQP9 isoforms and tested for changes in inhibitory effects. Controls: Three independent cell lines were established for each mutated AQP9 isoform and functionality of mutant isoforms was established. Principal findings: We have identified putative binding sites of recently established AQP9 inhibitors. This information facilitated successful identification of novel AQP9 inhibitors with low micromolar IC50 values in a cell based assay by in silico screening of a compound library targeting specifically this binding site. Significance: We have established a successful strategy for AQP small molecule inhibitor identification. AQP inhibitors may be relevant as experimental tools, to enhance our understanding of AQP function, and in the treatment of various diseases.