Accuracy Of Virtual Screening Research Articles

How to Improve Docking Accuracy of AutoDock4.2: A Case Study Using Different Electrostatic Potentials

Molecular docking, which is the indispensable emphasis in predicting binding conformations and energies of ligands to receptors, constructs the high-throughput virtual screening available. So far, increasingly numerous molecular docking programs have been released, and among them, AutoDock 4.2 is a widely used docking program with exceptional accuracy. It has heretofore been substantiated that the calculation of partial charge is very fundamental for the accurate conformation search and binding energy estimation. However, no systematic comparison of the significances of electrostatic potentials on docking accuracy of AutoDock 4.2 has been determined. In this paper, nine different charge-assigning methods, including AM1-BCC, Del-Re, formal, Gasteiger-Hückel, Gasteiger-Marsili, Hückel, Merck molecular force field (MMFF), and Pullman, as well as the ab initio Hartree-Fock charge, were sufficiently explored for their molecular docking performance by using AutoDock4.2. The results clearly demonstrated that the empirical Gasteiger-Hückel charge is the most applicable in virtual screening for large database; meanwhile, the semiempirical AM1-BCC charge is practicable in lead compound optimization as well as accurate virtual screening for small databases.

A Comparative Assessment of Ranking Accuracies of Conventional and Machine-Learning-Based Scoring Functions for Protein-Ligand Binding Affinity Prediction

Accurately predicting the binding affinities of large sets of protein-ligand complexes efficiently is a key challenge in computational biomolecular science, with applications in drug discovery, chemical biology, and structural biology. Since a scoring function (SF) is used to score, rank, and identify drug leads, the fidelity with which it predicts the affinity of a ligand candidate for a protein's binding site has a significant bearing on the accuracy of virtual screening. Despite intense efforts in developing conventional SFs, which are either force-field based, knowledge-based, or empirical, their limited ranking accuracy has been a major roadblock toward cost-effective drug discovery. Therefore, in this work, we explore a range of novel SFs employing different machine-learning (ML) approaches in conjunction with a variety of physicochemical and geometrical features characterizing protein-ligand complexes. We assess the ranking accuracies of these new ML-based SFs as well as those of conventional SFs in the context of the 2007 and 2010 PDBbind benchmark data sets on both diverse and protein-family-specific test sets. We also investigate the influence of the size of the training data set and the type and number of features used on ranking accuracy. Within clusters of protein-ligand complexes with different ligands bound to the same target protein, we find that the best ML-based SF is able to rank the ligands correctly based on their experimentally determined binding affinities 62.5 percent of the time and identify the top binding ligand 78.1 percent of the time. For this SF, the Spearman correlation coefficient between ranks of ligands ordered by predicted and experimentally determined binding affinities is 0.771. Given the challenging nature of the ranking problem and that SFs are used to screen millions of ligands, this represents a significant improvement over the best conventional SF we studied, for which the corresponding ranking performance values are 57.8 percent, 73.4 percent, and 0.677.

HLA‐DP2 binding prediction by molecular dynamics simulations

Major histocompatibility complex (MHC) II proteins bind peptide fragments derived from pathogen antigens and present them at the cell surface for recognition by T cells. MHC proteins are divided into Class I and Class II. Human MHC Class II alleles are grouped into three loci: HLA-DP, HLA-DQ, and HLA-DR. They are involved in many autoimmune diseases. In contrast to HLA-DR and HLA-DQ proteins, the X-ray structure of the HLA-DP2 protein has been solved quite recently. In this study, we have used structure-based molecular dynamics simulation to derive a tool for rapid and accurate virtual screening for the prediction of HLA-DP2-peptide binding. A combinatorial library of 247 peptides was built using the "single amino acid substitution" approach and docked into the HLA-DP2 binding site. The complexes were simulated for 1 ns and the short range interaction energies (Lennard-Jones and Coulumb) were used as binding scores after normalization. The normalized values were collected into quantitative matrices (QMs) and their predictive abilities were validated on a large external test set. The validation shows that the best performing QM consisted of Lennard-Jones energies normalized over all positions for anchor residues only plus cross terms between anchor-residues.

Local kernel canonical correlation analysis with application to virtual drug screening

Drug discovery is the process of identifying compounds which have potentially meaningful biological activity. A major challenge that arises is that the number of compounds to search over can be quite large, sometimes numbering in the millions, making experimental testing intractable. For this reason computational methods are employed to filter out those compounds which do not exhibit strong biological activity. This filtering step, also called virtual screening reduces the search space, allowing for the remaining compounds to be experimentally tested.In this paper we propose several novel approaches to the problem of virtual screening based on Canonical Correlation Analysis (CCA) and on a kernel-based extension. Spectral learning ideas motivate our proposed new method called Indefinite Kernel CCA (IKCCA). We show the strong performance of this approach both for a toy problem as well as using real world data with dramatic improvements in predictive accuracy of virtual screening over an existing methodology.

On the Integration of Pharmacophore Model and Molecular Docking Method

Structure-based molecular docking has been recognized as the most effective virtual screening technology. Two aspects decide the accuracy of virtual screening result: one is the simulation degree to the interaction of ligands and acceptors, the other is the quantitative evaluation technology to the result of interaction. The pharmacophore model is an important method of drug design, which can be used for the discovery of lead compounds. The development of scoring-technique based on the pharmacophore is an important supplement and improvement to the quantitative evaluation technology of the molecular docking result, which can enhance the accuracy and reliability of virtual screening based on the molecular docking technology, and impulse the development of new drug molecular design technique. Therefore, the Molecular Docking method of integration pharmacophore model is necessary for researching the creation of new drugs. This paper focuses on researching following content: 1) Pharmacophore of digital processing; explore the use of appropriate mathematical model to describe the specific target pharmacophore. 2) The integration of molecular docking algorithm; the development of the combination of molecular docking algorithm based on the established mathematical model of pharmacophore. 3) Establishing the verification of the reliability of the methods; through screening the results of the known activity of the compound-library, testing the reliability of the method.

FRED Pose Prediction and Virtual Screening Accuracy

Results of a previous docking study are reanalyzed and extended to include results from the docking program FRED and a detailed statistical analysis of both structure reproduction and virtual screening results. FRED is run both in a traditional docking mode and in a hybrid mode that makes use of the structure of a bound ligand in addition to the protein structure to screen molecules. This analysis shows that most docking programs are effective overall but highly inconsistent, tending to do well on one system and poorly on the next. Comparing methods, the difference in mean performance on DUD is found to be statistically significant (95% confidence) 61% of the time when using a global enrichment metric (AUC). Early enrichment metrics are found to have relatively poor statistical power, with 0.5% early enrichment only able to distinguish methods to 95% confidence 14% of the time.

Electron Density Fingerprints (EDprints): Virtual Screening Using Assembled Information of Electron Density

We have designed a method to encode properties related to the electron densities of molecules (calculated (1)H and (13)C NMR shifts and atomic partial charges) in molecular fingerprints (EDprints). EDprints was evaluated in terms of their retrospective virtual screening accuracy against the Directory of Useful Decoys (DUD) and compared to the established ligand-based similarity search methods MOLPRINT 2D and FCFP-4. Although there are no significant differences in the overall virtual screening accuracies of the three methods, specific examples highlight interesting differences between the new EDprints fingerprint method and the atom-centered circular fingerprint methods of MOLPRINT 2D and FCFP-4. On one hand, EDprints similarity searches can be biased by the molecular protonation state, especially when reference ligands contain multiple ionizable groups. On the other hand, EDprints models are more robust toward subtle rearrangements of chemical groups and more suitable for screening against reference molecules with fused ring systems than MOLPRINT 2D and FCFP-4. EDprints is furthermore the fastest method under investigation in comparing fingerprints (average 56-233-fold increase in speed), which makes it highly suitable for all-against-all similarity searches and for repetitive virtual screening against large chemical databases of millions of compounds.

Comparison of Several Molecular Docking Programs: Pose Prediction and Virtual Screening Accuracy

Molecular docking programs are widely used modeling tools for predicting ligand binding modes and structure based virtual screening. In this study, six molecular docking programs (DOCK, FlexX, GLIDE, ICM, PhDOCK, and Surflex) were evaluated using metrics intended to assess docking pose and virtual screening accuracy. Cognate ligand docking to 68 diverse, high-resolution X-ray complexes revealed that ICM, GLIDE, and Surflex generated ligand poses close to the X-ray conformation more often than the other docking programs. GLIDE and Surflex also outperformed the other docking programs when used for virtual screening, based on mean ROC AUC and ROC enrichment values obtained for the 40 protein targets in the Directory of Useful Decoys (DUD). Further analysis uncovered general trends in accuracy that are specific for particular protein families. Modifying basic parameters in the software was shown to have a significant effect on docking and virtual screening results, suggesting that expert knowledge is critical for optimizing the accuracy of these methods.

Molecular modeling of the second extracellular loop of G‐protein coupled receptors and its implication on structure‐based virtual screening

The current study describes the validation of high-throughput modeling procedures for the construction of the second extracellular loop (ecl2) of all nonolfactory human G Protein-coupled receptors. Our modeling flowchart is based on the alignment of essential residues determining the particular ecl2 fold observed in the bovine rhodopsin (bRho) crystal structure. For a set of GPCR targets, the dopamine D2 receptor (DRD2), adenosine A3 receptor (AA3R), and the thromboxane A2 receptor (TA2R), the implications of including ecl2 atomic coordinates is evaluated in terms of structure-based virtual screening accuracy: the suitability of the 3D models to distinguish between known antagonists and randomly chosen decoys using automated docking approaches. The virtual screening results of different models describing increasingly exhaustive receptor representations (seven helices only, seven helices and ecl2 loop, full model) have been compared. Explicit modeling of the ecl2 loop was found to be important in only one of three test cases whereas a loopless model was shown to be accurate enough in the two other receptors. An exhaustive comparison of ecl2 loops of 365 receptors to that of bRho suggests that explicit ecl2 loop modeling should be reserved to receptors where loop building can be guided by experimental restraints.

Catalytic Site Prediction and Virtual Screening of Cytochrome P450 2D6 Substrates by Consideration of Water and Rescoring in Automated Docking

Automated docking strategies successfully applied to binding mode predictions of ligands in Cyt P450 crystal structures in an earlier study (de Graaf et al. J. Med. Chem. 2005, 7, 2308-2318) were used for the catalytic site prediction (CSP) of 65 substrates in a CYP2D6 homology model. The consideration of water molecules at predicted positions in the active site and the rescoring of pooled docking poses from four different docking programs (AutoDock, FlexX, GOLD-Goldscore, and GOLD-Chemscore) with the SCORE scoring function enabled the successful prediction of experimentally reported sites of catalysis of more than 80% of the substrates. Three docking algorithms (FlexX, GOLD-Goldscore, and GOLD-Chemscore) were subsequently used in combination with six scoring functions (Chemscore, DOCK, FlexX, GOLD, PMF, and SCORE) to assess the ability of docking-based virtual screening methods to prioritize known CYP2D6 substrates seeded into a drug-like chemical database (in the absence and presence of active-site water molecules). Finally, the optimal docking strategy in terms of virtual screening accuracy, GOLD-Chemscore with the consideration of active-site water (60% of known substrates recovered in the top 5% of the ranked drug-like database), was verified experimentally; it was successfully used to identify high-affinity CYP2D6 ligands among a larger proprietary database.

Enhancing the accuracy of virtual screening: molecular dynamics with quantum-refined force fields

A methodology aimed at improving the accuracy of current docking-scoring procedures is proposed, and validated through detailed tests of its performance in predicting the activity of HIV-1 protease inhibitors. This methodology is based on molecular dynamics simulations using a force field whose effective charges are refined by means of a novel procedure that relies on quantum-mechanical calculations and preserves the internal consistency of the parameterization scheme.

Pharmacophore-based molecular docking to account for ligand flexibility.

Rapid computational mining of large 3D molecular databases is central to generating new drug leads. Accurate virtual screening of large 3D molecular databases requires consideration of the conformational flexibility of the ligand molecules. Ligand flexibility can be included without prohibitively increasing the search time by docking ensembles of precomputed conformers from a conformationally expanded database. A pharmacophore-based docking method whereby conformers of the same or different molecules are overlaid by their largest 3D pharmacophore and simultaneously docked by partial matches to that pharmacophore is presented. The method is implemented in DOCK 4.0.