Energetic analysis of fragment docking and application to structure-based pharmacophore hypothesis generation

We have developed a method that uses energetic analysis of structure-based fragment docking to elucidate key features for molecular recognition. This hybrid ligand- and structure-based methodology uses an atomic breakdown of the energy terms from the Glide XP scoring function to locate key pharmacophoric features from the docked fragments. First, we show that Glide accurately docks fragments, producing a root mean squared deviation (RMSD) of <1.0 Å for the top scoring pose to the native crystal structure. We then describe fragment-specific docking settings developed to generate poses that explore every pocket of a binding site while maintaining the docking accuracy of the top scoring pose. Next, we describe how the energy terms from the Glide XP scoring function are mapped onto pharmacophore sites from the docked fragments in order to rank their importance for binding. Using this energetic analysis we show that the most energetically favorable pharmacophore sites are consistent with features from known tight binding compounds. Finally, we describe a method to use the energetically selected sites from fragment docking to develop a pharmacophore hypothesis that can be used in virtual database screening to retrieve diverse compounds. We find that this method produces viable hypotheses that are consistent with known active compounds. In addition to retrieving diverse compounds that are not biased by the co-crystallized ligand, the method is able to recover known active compounds from a database screen, with an average enrichment of 8.1 in the top 1% of the database.     

Lessons for fragment library design: analysis of output from multiple screening campaigns

Over the past 8 years, we have developed, refined and applied a fragment based discovery approach to a range of protein targets. Here we report computational analyses of various aspects of our fragment library and the results obtained for fragment screening. We reinforce the finding of others that the experimentally observed hit rate for screening fragments can be related to a computationally defined druggability index for the target. In general, the physicochemical properties of the fragment hits display the same profile as the library, as is expected for a truly diverse library which probes the relevant chemical space. An analysis of the fragment hits against various protein classes has shown that the physicochemical properties of the fragments are complementary to the properties of the target binding site. The effectiveness of some fragments appears to be achieved by an appropriate mix of pharmacophore features and enhanced aromaticity, with hydrophobic interactions playing an important role. The analysis emphasizes that it is possible to identify small fragments that are specific for different binding sites. To conclude, we discuss how the results could inform further development and improvement of our fragment library. Electronic supplementary material  The online version of this article (doi:) contains supplementary material, which is available to authorized users.     

Virtual screening and rational drug design method using structure generation system based on 3D-QSAR and docking

An efficient virtual and rational drug design method is presented. It combines virtual bioactive compound generation with 3D-QSAR model and docking. Using this method, it is possible to generate a lot of highly diverse molecules and find virtual active lead compounds. The method was validated by the study of a set of anti-tumor drugs. With the constraints of pharmacophore obtained by DISCO implemented in SYBYL 6.8, 97 virtual bioactive compounds were generated, and their anti-tumor activities were predicted by CoMFA. Eight structures with high activity were selected and screened by the 3D-QSAR model. The most active generated structure was further investigated by modifying its structure in order to increase the activity. A comparative docking study with telomeric receptor was carried out, and the results showed that the generated structures could form more stable complexes with receptor than the reference compound selected from experimental data. This investigation showed that the proposed method was a feasible way for rational drug design with high screening efficiency.     

Design of compound libraries for fragment screening

Approaches to the design of libraries for fragment screening are illustrated with reference to a 20 k generic fragment screening library and a 1.2 k generic NMR screening library. Tools and methods for library design that have been developed within AstraZeneca are described, including Foyfi fingerprints and the Flush program for neighborhood characterization. It will be shown how Flush and the BigPicker, which selects maximally diverse sets of compounds, are used to apply the Core and Layer method for library design. Approaches to partitioning libraries into cocktails are also described.     

Distilling the essential features of a protein surface for improving protein-ligand docking,scoring, and virtual screening

For the successful identification and docking of new ligands to a protein target by virtual screening, the essential features of the protein and ligand surfaces must be captured and distilled in an efficient representation. Since the running time for docking increases exponentially with the number of points representing the protein and each ligand candidate, it is important to place these points where the best interactions can be made between the protein and the ligand. This definition of favorable points of interaction can also guide protein structure-based ligand design, which typically focuses on which chemical groups provide the most energetically favorable contacts. In this paper, we present an alternative method of protein template and ligand interaction point design that identifies the most favorable points for making hydrophobic and hydrogen–bond interactions by using a knowledge base. The knowledge-based protein and ligand representations have been incorporated in version 2.0 of SLIDE and resulted in dockings closer to the crystal structure orientations when screening a set of 57 known thrombin and glutathione S–transferase (GST) ligands against the apo structures of these proteins. There was also improved scoring enrichment of the dockings, meaning better differentiation between the chemically diverse known ligands and a 15,000-molecule dataset of randomly-chosen small organic molecules. This approach for identifying the most important points of interaction between proteins and their ligands can equally well be used in other docking and design techniques. While much recent effort has focused on improving scoring functions for protein-ligand docking, our results indicate that improving the representation of the chemistry of proteins and their ligands is another avenue that can lead to significant improvements in the identification, docking, and scoring of ligands.(These authors contributed equally to this work)     

Identification of novel nt-MGAM inhibitors for potential treatment of type 2 diabetes: Virtual screening,atom based 3D-QSAR model,docking analysis and ADME study

In this study, a virtual screening procedure was applied to identify new potential nt-MGAM inhibitors as a possible medication for type 2 diabetes. To this aim, a series of salacinol analogues were first investigated by docking analysis for their binding to the X-ray structure of the biological target nt-MGAM. Key interactions for ligand binding into the receptor active site were identified which shared common features to those found for other known inhibitors, which strengthen the results of this study. 3D QSAR model was then built and showed to be statistically significant and with a good predictive power for the training (R² = 0.99, SD = 0.17, F = 555.3 and N = 27) and test set (Q² = 0.81, Pearson(r) = 0.92, RMSE = 0.52, N = 08). The model was then used to virtually screen the ZINC database with the aim of identifying novel chemical scaffolds as potential nt-MGAM inhibitors. Further, in silico predicted ADME properties were investigated for the most promising molecules. The outcome of this investigation sheds light on the molecular characteristics of the binding of salacinol analogues to nt-MGAM enzyme and identifies new possible inhibitors which have the potential to be developed into drugs, thus significantly contributing to the design and optimization of therapeutic strategies against type 2 diabetes.     

Virtual screening of B-Raf kinase inhibitors: A combination of pharmacophore modelling,molecular docking, 3D-QSAR model and binding free energy calculation studies

B-Raf kinase has been identified as an important target in recent cancer treatment. In order to discover structurally diverse and novel B-Raf inhibitors (BRIs), a virtual screening of BRIs against ZINC database was performed by using a combination of pharmacophore modelling, molecular docking, 3D-QSAR model and binding free energy (ΔG_bind) calculation studies in this work. After the virtual screening, six promising hit compounds were obtained, which were then tested for inhibitory activities of A375 cell lines. In the result, five hit compounds show good biological activities (IC₅₀ < 50 μM). The present method of virtual screening can be applied to find structurally diverse inhibitors, and the obtained five structurally diverse compounds are expected to develop novel BRIs.     

Critical evaluation of search algorithms for automated molecular docking and database screening

The DOCK program explores possible orientations of a molecule within a macromolecular active site by superimposing atoms onto precomputed site points. Here we compare a number of different search methods, including an exhaustive matching algorithm based on a single docking graph. We evaluate the performance of each method by screening a small database of molecules to a variety of macromolecular targets. By varying the amount of sampling, we can monitor the time convergence of scores and rankings. We not only show that the site point–directed search is tenfold faster than a random search, but that the single graph matching algorithm boosts the speed of database screening up to 60-fold. The new algorithm, in fact, outperforms the bipartite graph matching algorithm currently used in DOCK. The results indicate that a critical issue for rapid database screening is the extent to which a search method biases run time toward the highest-ranking molecules. The single docking graph matching algorithm will be incorporated into DOCK version 4.0. © 1997 John Wiley & Sons, Inc. J Comput Chem 18: 1175–1189     

AutoDock Vina: Improving the speed and accuracy of docking with a new scoring function,efficient optimization,and multithreading

AutoDock Vina, a new program for molecular docking and virtual screening, is presented. AutoDock Vina achieves an approximately two orders of magnitude speed‐up compared with the molecular docking software previously developed in our lab (AutoDock 4), while also significantly improving the accuracy of the binding mode predictions, judging by our tests on the training set used in AutoDock 4 development. Further speed‐up is achieved from parallelism, by using multithreading on multicore machines. AutoDock Vina automatically calculates the grid maps and clusters the results in a way transparent to the user. © 2009 Wiley Periodicals, Inc. J Comput Chem 2010     

Deep Learning and Structure-Based Virtual Screening for Drug Discovery against NEK7: A Novel Target for the Treatment of Cancer

NIMA-related kinase7 (NEK7) plays a multifunctional role in cell division and NLRP3 inflammasone activation. A typical expression or any mutation in the genetic makeup of NEK7 leads to the development of cancer malignancies and fatal inflammatory disease, i.e., breast cancer, non-small cell lung cancer, gout, rheumatoid arthritis, and liver cirrhosis. Therefore, NEK7 is a promising target for drug development against various cancer malignancies. The combination of drug repurposing and structure-based virtual screening of large libraries of compounds has dramatically improved the development of anticancer drugs. The current study focused on the virtual screening of 1200 benzene sulphonamide derivatives retrieved from the PubChem database by selecting and docking validation of the crystal structure of NEK7 protein (PDB ID: 2WQN). The compounds library was subjected to virtual screening using Auto Dock Vina. The binding energies of screened compounds were compared to standard Dabrafenib. In particular, compound 762 exhibited excellent binding energy of −42.67 kJ/mol, better than Dabrafenib (−33.89 kJ/mol). Selected drug candidates showed a reactive profile that was comparable to standard Dabrafenib. To characterize the stability of protein–ligand complexes, molecular dynamic simulations were performed, providing insight into the molecular interactions. The NEK7–Dabrafenib complex showed stability throughout the simulated trajectory. In addition, binding affinities, pIC50, and ADMET profiles of drug candidates were predicted using deep learning models. Deep learning models predicted the binding affinity of compound 762 best among all derivatives, which supports the findings of virtual screening. These findings suggest that top hits can serve as potential inhibitors of NEK7. Moreover, it is recommended to explore the inhibitory potential of identified hits compounds through in-vitro and in-vivo approaches.