Autodock Vina与Discovery Studio在虚拟筛选耐药蛋白抑制剂中的比较

随着大量与细菌耐药相关的基因的发现和其表达蛋白结构的成功测定,从已有的化合物中通过计算机模拟方法筛选对耐药蛋白靶点有作用的候选化合物,成为了药物发现的一个标准途径。虚拟筛选在耐药基因抑制剂的发现中可以提高效率、降低实验成本。本文介绍了Autodock Vina和Discovery Studio在基于分子对接法的虚拟筛选中的使用,并对比分析其对β-内酰胺酶活性位点的筛选结果。希望通过这种比较促进虚拟筛选在药物设计领域中的应用,提高耐药基因抑制剂的发现速度。     

Emerging frontiers in virtual drug discovery: From quantum mechanical methods to deep learning approaches

Virtual screening-based approaches to discover initial hit and lead compounds have the potential to reduce both the cost and time of early drug discovery stages, as well as to find inhibitors for even challenging target sites such as protein–protein interfaces. Here in this review, we provide an overview of the progress that has been made in virtual screening methodology and technology on multiple fronts in recent years. The advent of ultra-large virtual screens, in which hundreds of millions to billions of compounds are screened, has proven to be a powerful approach to discover highly potent hit compounds. However, these developments are just the tip of the iceberg, with new technologies and methods emerging to propel the field forward. Examples include novel machine-learning approaches, which can reduce the computational costs of virtual screening dramatically, while progress in quantum-mechanical approaches can increase the accuracy of predictions of various small molecule properties.     

Molecular modeling of drug-DNA interactions: virtual screening to structure-based design

Virtual ligand screening (VLS) and structure-based design are strategies that have been routinely used for the development of pharmaceuticals, particularly those targeting enzymes and other protein targets. In recent years, an increased understanding of the role played by nucleic acids in biological systems made DNA an alternative candidate for the development of new drugs. This review highlights some successful applications of molecular modeling in virtual ligand screening and structure-based design of organic and inorganic molecules that target non-canonical nucleic acid structures such as G-quadruplex and triplex DNA.     

Evaluation of molecular model-based discovery of ecto-5′-nucleotidase inhibitors on the basis of X-ray structures

The enzyme ecto-5′-nucleotidase (e5NT, CD73), a metallophosphoesterase, is a critical component of adenosine metabolism and signaling and implicated in different disease states. Therefore, attempts have been made to discover inhibitors of e5NT. For example, a virtual screening study using a molecular model of the enzyme has led to the identification of a new series of sulfonamide-containing e5NT inhibitors. The recent availability of several X-ray structures of human e5NT in complex with inhibitors has made it possible to re-evaluate this model building and virtual screening effort. We have assessed the quality of the model in detail and analyzed the question why it was possible to identify a new series of inhibitors on the basis of model-based docking calculations. The model utilized for virtual screening was found to be topologically correct and approach experimental accuracy in the active site region. Two key features within the active site were identified as major determinants for the successfully identification of inhibitors. Taken together, the results rationalize the computer-aided discovery of sulfonamide inhibitors of e5NT and provide further support for the use of carefully built protein models for virtual screening.     

Critical assessment of the automated AutoDock as a new docking tool for virtual screening

A major problem in virtual screening concerns the accuracy of the binding free energy between a target protein and a putative ligand. Here we report an example supporting the outperformance of the AutoDock scoring function in virtual screening in comparison to the other popular docking programs. The original AutoDock program is in itself inefficient to be used in virtual screening because the grids of interaction energy have to be calculated for each putative ligand in chemical database. However, the automation of the AutoDock program with the potential grids defined in common for all putative ligands leads to more than twofold increase in the speed of virtual database screening. The utility of the automated AutoDock in virtual screening is further demonstrated by identifying the actual inhibitors of various target enzymes in chemical databases with accuracy higher than the other docking tools including DOCK and FlexX. These results exemplify the usefulness of the automated AutoDock as a new promising tool in structure-based virtual screening.     

Scoring functions and enrichment: a case study on Hsp90

Background  

The need for fast and accurate scoring functions has been driven by the increased use of in silico virtual screening twinned with high-throughput screening as a method to rapidly identify potential candidates in the early stages of drug development. We examine the ability of some the most common scoring functions (GOLD, ChemScore, DOCK, PMF, BLEEP and Consensus) to discriminate correctly and efficiently between active and non-active compounds among a library of ~3,600 diverse decoy compounds in a virtual screening experiment against heat shock protein 90 (Hsp90).     

虚拟筛选与新药发现

虚拟筛选是创新药物研究的新方法和新技术,近年来引起了研究机构和制药公司的高度重视,并且已经成为一种与高通量筛选互补的实用化工具,加入到了创新药物研究的工作流程(pipeline)中。本文介绍国际上虚拟筛选及其在创新药物发现中应用的研究进展,特别介绍了我国这方面研究的状况。     

Protein flexibility in ligand docking and virtual screening to protein kinases

The main complicating factor in structure-based drug design is receptor rearrangement upon ligand binding (induced fit). It is the induced fit that complicates cross-docking of ligands from different ligand-receptor complexes. Previous studies have shown the necessity to include protein flexibility in ligand docking and virtual screening. Very few docking methods have been developed to predict the induced fit reliably and, at the same time, to improve on discriminating between binders and non-binders in the virtual screening process.We present an algorithm called the ICM-flexible receptor docking algorithm (IFREDA) to account for protein flexibility in virtual screening. By docking flexible ligands to a flexible receptor, IFREDA generates a discrete set of receptor conformations, which are then used to perform flexible ligand-rigid receptor docking and scoring. This is followed by a merging and shrinking step, where the results of the multiple virtual screenings are condensed to improve the enrichment factor. In the IFREDA approach, both side-chain rearrangements and essential backbone movements are taken into consideration, thus sampling adequately the conformational space of the receptor, even in cases of large loop movements.As a preliminary step, to show the importance of incorporating protein flexibility in ligand docking and virtual screening, and to validate the merging and shrinking procedure, we compiled an extensive small-scale virtual screening benchmark of 33 crystal structures of four different protein kinases sub-families (cAPK, CDK-2, P38 and LCK), where we obtained an enrichment factor fold-increase of 1.85±0.65 using two or three multiple experimental conformations. IFREDA was used in eight protein kinase complexes and was able to find the correct ligand conformation and discriminate the correct conformations from the “misdocked” conformations solely on the basis of energy calculation. Five of the generated structures were used in the small-scale virtual screening stage and, by merging and shrinking the results with those of the original structure, we show an enrichment factor fold increase of 1.89±0.60, comparable to that obtained using multiple experimental conformations.Our cross-docking tests on the protein kinase benchmark underscore the necessity of incorporating protein flexibility in both ligand docking and virtual screening. The methodology presented here will be extremely useful in cases where few or no experimental structures of complexes are available, while some binders are known.     

Novel and potent NPY5 receptor antagonists derived from virtual screening and iterative parallel chemistry design

In the quest for NPY5 receptor antagonists a virtual screening approach yielded a novel and potent hit class from a limited compound selection. The tight and seamless integration between virtual screening and rapid parallel chemistry within the framework of the Roche Lead Generation unit led in only two rounds of iterative chemistry optimisation to a much broader understanding of the factors which influence the potency of the thiazole hit class.     

Novel GSK-3beta inhibitors from sequential virtual screening

Glycogen synthase kinase-3 (GSK-3beta) has been emerging as a key therapeutic target for type-2 diabetics, Alzheimer's disease, cancer, and chronic inflammation. For the purpose of finding biologically active and novel compounds and providing new idea for drug-design, we performed virtual screening using commercially available database. Three-dimensional common feature pharmacophore model was developed by using HipHop program provided in Catalyst software and it was used as a query for screening database. Recursive partitioning (RP) model was developed as a filtering system, which was able to classify active and inactive compounds. Eventually, a sequential virtual screening procedure (SQSP) was conducted by applying the common feature pharmacophore and RP model in succession to discover novel potent GSK-3beta inhibitors. The final 56 hit compounds were carefully selected considering predicted docking mode in crystal structures. Subsequent enzyme assay for human GSK-3beta protein confirmed that three compounds of these hit compounds exhibit micromolar inhibitory activity. Here, we report novel hit compounds and their binding mode in the active site of GSK-3beta crystal structure.     

Identification of potential inhibitors for HCV NS3 genotype 4a by combining protein–ligand interaction fingerprint, 3D pharmacophore,docking, and dynamic simulation

HCV NS3 protease domain has been one of the most attractive targets for developing new drugs for HCV infection and many drugs were successfully developed, but all of them were designed for targeting HCV genotype 1 infection. HCV genotype 4a dominant in Egypt has paid less attention. Here, we describe our protocol of virtual screening in identification of novel potential potent inhibitors for HCV NS3 of genotype 4a using homology modeling, PLIF (protein–ligand interaction fingerprint), docking, pharmacophore, and dynamic simulation. A high-quality 3D model of HCV NS3 protease of genotype 4a was constructed using crystal structure of HCV NS3 protease of genotype 1b (PDB ID: 4u01) as a template. PLIF was generated using five crystal structures of HCV NS3 (PDB ID: 4u01, 3kee, 4ktc, 4i33, and 5epn) which revealed the most important residues and their interactions with the co-crystalized ligands. A 3D pharmacophore model consisting of six features was developed from the generated PLIF data and then used as a screening filter for 11,244 compounds. Only 423 compounds passed the pharmacophore filter and entered the docking-based virtual screening stage. The highest ranked five hits from docking result (compound (C1–C5)) were selected for further analysis. They exhibited stronger interaction and higher binding affinity than HCV NS3 protease ligands. Dynamic simulation of the protein–best lead complex was performed to validate and augment the virtual screening results and it showed that these compounds have a strong binding affinity and could be very effective in treating HCV genotype 4a infections.     

Prediction of molecular targets of cancer preventing flavonoid compounds using computational methods

Plant-based polyphenols (i.e., phytochemicals) have been used as treatments for human ailments for centuries. The mechanisms of action of these plant-derived compounds are now a major area of investigation. Thousands of phytochemicals have been isolated, and a large number of them have shown protective activities or effects in different disease models. Using conventional approaches to select the best single or group of best chemicals for studying the effectiveness in treating or preventing disease is extremely challenging. We have developed and used computational-based methodologies that provide efficient and inexpensive tools to gain further understanding of the anticancer and therapeutic effects exerted by phytochemicals. Computational methods involving virtual screening, shape and pharmacophore analysis and molecular docking have been used to select chemicals that target a particular protein or enzyme and to determine potential protein targets for well-characterized as well as for novel phytochemicals.     

How to choose relevant multiple receptor conformations for virtual screening: a test case of Cdk2 and normal mode analysis

Better treatment of protein flexibility is essential in structure-based drug design projects such as virtual screening and protein-ligand docking. Diversity in ligand-binding mechanisms and receptor conformational changes makes it difficult to treat dynamic features of the receptor during the docking simulation. Thus, the use of pregenerated multiple receptor conformations is applied today in virtual screening studies. However, generation of a small relevant set of receptor conformations remains challenging. To address this problem, we propose a new protocol for the generation of multiple receptor conformations via normal mode analysis and for the selection of several receptor conformations suitable for docking/virtual screening. We validated this protocol on cyclin-dependent kinase 2, which possesses a binding site located at the interface between two subdomains and is known to undergo significant conformational changes in the active site region upon ligand binding. We believe that the suggested rules for the choice of suitable receptor conformations can be applied to other targets when dealing with in silico screening on flexible receptors.     

Towards structure-based protein drug design

Structure-based drug design for chemical molecules has been widely used in drug discovery in the last 30 years. Many successful applications have been reported, especially in the field of virtual screening based on molecular docking. Recently, there has been much progress in fragment-based as well as de novo drug discovery. As many protein-protein interactions can be used as key targets for drug design, one of the solutions is to design protein drugs based directly on the protein complexes or the target structure. Compared with protein-ligand interactions, protein-protein interactions are more complicated and present more challenges for design. Over the last decade, both sampling efficiency and scoring accuracy of protein-protein docking have increased significantly. We have developed several strategies for structure-based protein drug design. A grafting strategy for key interaction residues has been developed and successfully applied in designing erythropoietin receptor-binding proteins. Similarly to small-molecule design, we also tested de novo protein-binder design and a virtual screen of protein binders using protein-protein docking calculations. In comparison with the development of structure-based small-molecule drug design, we believe that structure-based protein drug design has come of age.     

The use of chemical design tools to transform proteomics data into drug candidates

The Human Genome Project has fueled the massive information-driven growth of genomics and proteomics and promises to deliver new insights into biology and medicine. Since proteins represent the majority of drug targets, these molecules are the focus of activity in pharmaceutical and biotechnology organizations. In this article, we describe the processes by which computational drug design may be used to exploit protein structural information to create virtual small molecules that may become novel medicines. Experimental protein structure determination, site exploration, and virtual screening provide a foundation for small molecule generation in silico, thus creating the bridge between proteomics and drug discovery.     

Structure based virtual screening to identify inhibitors against MurE Enzyme of Mycobacterium tuberculosis using AutoDock Vina

The Mur E enzyme of Mur pathway of Mycobacterium tuberculosis is an attractive drug target as it is unique to bacteria and is absent in mammalian cells. The virtual screening of large libraries of drug like molecules against a protein target is a common strategy used to identify novel inhibitors. However, the method has a large number of pitfalls, with large variations in accuracy caused in part by inaccurate protocols, use of improper standards and libraries, and system dependencies such as the potential for nonspecific docking from large active-site cavities. The screening of drug-like small molecules from diversity sets can, however, be used to short-list potential fragments as building blocks to generate leads with improved specificity. We describe a protocol to implement this strategy, which involves an analysis of the active site and known inhibitors to identify orthospecific determinants, virtual screening of a drug-like diversity library to identify potential drug primitives, and inspection of the potential docked fragments for both binding potential and toxicity. The protocol is implemented on the M.tb Mur E protein which has a large active site with poor enrichment of known positives and a set of drug-like molecules that meets this criteria is presented for further analysis.

Abbreviations

MTB - Mycobacterium tuberculosis, NCI - National Cancer Institute, PDB - Protein Databank.     

GABPalpha regulates Oct-3/4 expression in mouse embryonic stem cells

There is a dire need for novel therapeutics to treat the virulent malarial parasite, Plasmodium falciparum. Recently, the X-ray crystal structure of enoyl-acyl carrier protein reductase (ENR) in complex with triclosan has been determined and provides an opportunity for the rational design of novel inhibitors targeting the active site of ENR. Here, we report the discovery of several compounds by virtual screening and their experimental validation as high potency PfENR inhibitors.     

Structure Based Screening for Inhibitory Therapeutics of CTLA-4 Unveiled New Insights About Biology of ACTH

Although the biology of adrenocorticotropic hormone (ACTH) protein has already been scrutinized, some functional aspects of its biology are yet to be elucidated in the context of immunological disorders. In this regard, virtual screening of a compound library was performed against the structure of Cytotoxic T-Lymphocyte Associated Protein-4 (CTLA-4) (assessed both spatially and energetically) to discover novel biological functions for ACTH. The results of virtual screening and the MD simulation demonstrated that DB01284 has high binding energy along with proper interaction orientation against CTLA-4 (FG loop) by a clamp like structure. The employed methodology was checked using confirmatory control analyses. Intriguingly, DB01284 belongs to Tetracosactide (already prescribed protein drug for clinical conditions) which is the N-terminal region of ACTH. This is the first study to reveal that ACTH protein binds to the same amino acids of CTLA-4 (FG-loop) as B7 and anti-CTLA-4 antibody binds. In light of this finding, the molecular mechanism of ACTH function in patients suffering from Cushing’s Syndrom and the immunological bases for ACTH therapy of multiple sclerosis (MS) patients could be further delineated. Moreover, this finding suggests that ACTH could also act to block CTLA-4 in the context of anticancer immune check point blockade.

     

Efficient molecular docking of NMR structures: application to HIV-1 protease

Docking ligands into an ensemble of NMR conformers is essential to structure-based drug discovery if only NMR structures are available for the target. However, sequentially docking ligands into each NMR conformer through standard single-receptor-structure docking, referred to as sequential docking, is computationally expensive for large-scale database screening because of the large number of NMR conformers involved. Recently, we developed an efficient ensemble docking algorithm to consider protein structural variations in ligand binding. The algorithm simultaneously docks ligands into an ensemble of protein structures and achieves comparable performance to sequential docking without significant increase in computational time over single-structure docking. Here, we applied this algorithm to docking with NMR structures. The HIV-1 protease was used for validation in terms of docking accuracy and virtual screening. Ensemble docking of the NMR structures identified 91% of the known inhibitors under the criterion of RMSD < 2.0 A for the best-scored conformation, higher than the average success rate of single docking of individual crystal structures (66%). In the virtual screening test, on average, ensemble docking of the NMR structures obtained higher enrichments than single-structure docking of the crystal structures. In contrast, docking of either the NMR minimized average structure or a single NMR conformer performed less satisfactorily on both binding mode prediction and virtual screening, indicating that a single NMR structure may not be suitable for docking calculations. The success of ensemble docking of the NMR structures suggests an efficient alternative method for standard single docking of crystal structures and for considering protein flexibility.