Dataset Information

Using free energy of binding calculations to improve the accuracy of virtual screening predictions.

ABSTRACT: Virtual screening of small molecule databases against macromolecular targets was used to identify binding ligands and predict their lowest energy bound conformation (i.e., pose). AutoDock4-generated poses were rescored using mean-field pathway decoupling free energy of binding calculations and evaluated if these calculations improved virtual screening discrimination between bound and nonbound ligands. Two small molecule databases were used to evaluate the effectiveness of the rescoring algorithm in correctly identifying binders of L99A T4 lysozyme. Self-dock calculations of a database containing compounds with known binding free energies and cocrystal structures largely reproduced experimental measurements, although the mean difference between calculated and experimental binding free energies increased as the predicted bound poses diverged from the experimental poses. In addition, free energy rescoring was more accurate than AutoDock4 scores in discriminating between known binders and nonbinders, suggesting free energy rescoring could be a useful approach to reduce false positive predictions in virtual screening experiments.

SUBMITTER: Malmstrom RD

PROVIDER: S-EPMC3145410 | biostudies-literature | 2011 Jul

REPOSITORIES: biostudies-literature

ACCESS DATA

Publications

Using free energy of binding calculations to improve the accuracy of virtual screening predictions.

Malmstrom Robert D RD Watowich Stanley J SJ

Journal of chemical information and modeling 20110711 7

Virtual screening of small molecule databases against macromolecular targets was used to identify binding ligands and predict their lowest energy bound conformation (i.e., pose). AutoDock4-generated poses were rescored using mean-field pathway decoupling free energy of binding calculations and evaluated if these calculations improved virtual screening discrimination between bound and nonbound ligands. Two small molecule databases were used to evaluate the effectiveness of the rescoring algorithm ...[more]

PMID: 21696204

Similar Datasets

Project description:We describe binding free energy calculations in the D3R Grand Challenge 2015 for blind prediction of the binding affinities of 180 ligands to Hsp90. The present D3R challenge was built around experimental datasets involving Heat shock protein (Hsp) 90, an ATP-dependent molecular chaperone which is an important anticancer drug target. The Hsp90 ATP binding site is known to be a challenging target for accurate calculations of ligand binding affinities because of the ligand-dependent conformational changes in the binding site, the presence of ordered waters and the broad chemical diversity of ligands that can bind at this site. Our primary focus here is to distinguish binders from nonbinders. Large scale absolute binding free energy calculations that cover over 3000 protein-ligand complexes were performed using the BEDAM method starting from docked structures generated by Glide docking. Although the ligand dataset in this study resembles an intermediate to late stage lead optimization project while the BEDAM method is mainly developed for early stage virtual screening of hit molecules, the BEDAM binding free energy scoring has resulted in a moderate enrichment of ligand screening against this challenging drug target. Results show that, using a statistical mechanics based free energy method like BEDAM starting from docked poses offers better enrichment than classical docking scoring functions and rescoring methods like Prime MM-GBSA for the Hsp90 data set in this blind challenge. Importantly, among the three methods tested here, only the mean value of the BEDAM binding free energy scores is able to separate the large group of binders from the small group of nonbinders with a gap of 2.4 kcal/mol. None of the three methods that we have tested provided accurate ranking of the affinities of the 147 active compounds. We discuss the possible sources of errors in the binding free energy calculations. The study suggests that BEDAM can be used strategically to discriminate binders from nonbinders in virtual screening and to more accurately predict the ligand binding modes prior to the more computationally expensive FEP calculations of binding affinity.

Project description:The absolute (standard) binding free energy of eight FK506-related ligands to FKBP12 is calculated using free energy perturbation molecular dynamics (FEP/MD) simulations with explicit solvent. A number of features are implemented to improve the accuracy and enhance the convergence of the calculations. First, the absolute binding free energy is decomposed into sequential steps during which the ligand-surrounding interactions as well as various biasing potentials restraining the translation, orientation, and conformation of the ligand are turned "on" and "off." Second, sampling of the ligand conformation is enforced by a restraining potential based on the root mean-square deviation relative to the bound state conformation. The effect of all the restraining potentials is rigorously unbiased, and it is shown explicitly that the final results are independent of all artificial restraints. Third, the repulsive and dispersive free energy contribution arising from the Lennard-Jones interactions of the ligand with its surrounding (protein and solvent) is calculated using the Weeks-Chandler-Andersen separation. This separation also improves convergence of the FEP/MD calculations. Fourth, to decrease the computational cost, only a small number of atoms in the vicinity of the binding site are simulated explicitly, while all the influence of the remaining atoms is incorporated implicitly using the generalized solvent boundary potential (GSBP) method. With GSBP, the size of the simulated FKBP12/ligand systems is significantly reduced, from approximately 25,000 to 2500. The computations are very efficient and the statistical error is small ( approximately 1 kcal/mol). The calculated binding free energies are generally in good agreement with available experimental data and previous calculations (within approximately 2 kcal/mol). The present results indicate that a strategy based on FEP/MD simulations of a reduced GSBP atomic model sampled with conformational, translational, and orientational restraining potentials can be computationally inexpensive and accurate.

Dataset Information

Using free energy of binding calculations to improve the accuracy of virtual screening predictions.

Publications

Using free energy of binding calculations to improve the accuracy of virtual screening predictions.

Similar Datasets

OmicsDI is part of the ELIXIR infrastructure

Tweets