Using Free Energy of Binding Calculations to Improve the Accuracy of Virtual Screening Predictions

Overview

Journal J Chem Inf Model

Publisher American Chemical Society

Specialties Chemistry
Medical Informatics

Date 2011 Jun 24

PMID 21696204

Citations 9

Authors

Robert D Malmstrom

Stanley J Watowich

Affiliations

Soon will be listed here.

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.

Citing Articles

Virtual screening: hope, hype, and the fine line in between.

Nada H, Meanwell N, Gabr M Expert Opin Drug Discov. 2025; 20(2):145-162.

PMID: 39862145 PMC: 11844436. DOI: 10.1080/17460441.2025.2458666.

Synthesis, molecular docking, and binding Gibbs free energy calculation of β-nitrostyrene derivatives: Potential inhibitors of SARS-CoV-2 3CL protease.

Jia Z, Lan X, Lu K, Meng X, Jing W, Jia S J Mol Struct. 2023; 1284:135409.

PMID: 36993878 PMC: 10033154. DOI: 10.1016/j.molstruc.2023.135409.

Free Energy Calculations Using the Movable Type Method with Molecular Dynamics Driven Protein-Ligand Sampling.

Liu W, Liu Z, Liu H, Westerhoff L, Zheng Z J Chem Inf Model. 2022; 62(22):5645-5665.

PMID: 36282990 PMC: 9709919. DOI: 10.1021/acs.jcim.2c00278.

The anti-inflammatory potential of cefazolin as common gamma chain cytokine inhibitor.

Zyzynska-Granica B, Trzaskowski B, Dutkiewicz M, Zegrocka-Stendel O, Machcinska M, Bocian K Sci Rep. 2020; 10(1):2886.

PMID: 32076052 PMC: 7031511. DOI: 10.1038/s41598-020-59798-3.

PIP3-binding proteins promote age-dependent protein aggregation and limit survival in C. elegans.

Ayyadevara S, Balasubramaniam M, Johnson J, Alla R, Mackintosh S, Shmookler Reis R Oncotarget. 2016; 7(31):48870-48886.

PMID: 27429199 PMC: 5226477. DOI: 10.18632/oncotarget.10549.

References

Pang Y, Perola E, Xu K, Prendergast F . EUDOC: a computer program for identification of drug interaction sites in macromolecules and drug leads from chemical databases. J Comput Chem. 2002; 22(15):1750-1771. DOI: 10.1002/jcc.1129. View

Morton A, Baase W, Matthews B . Energetic origins of specificity of ligand binding in an interior nonpolar cavity of T4 lysozyme. Biochemistry. 1995; 34(27):8564-75. DOI: 10.1021/bi00027a006. View

Collins M, Hummer G, Quillin M, Matthews B, Gruner S . Cooperative water filling of a nonpolar protein cavity observed by high-pressure crystallography and simulation. Proc Natl Acad Sci U S A. 2005; 102(46):16668-71. PMC: 1283839. DOI: 10.1073/pnas.0508224102. View

Ravindranathan K, Mandiyan V, Ekkati A, Bae J, Schlessinger J, Jorgensen W . Discovery of novel fibroblast growth factor receptor 1 kinase inhibitors by structure-based virtual screening. J Med Chem. 2010; 53(4):1662-72. PMC: 2842983. DOI: 10.1021/jm901386e. View

MacKerell A, Bashford D, Bellott M, Dunbrack R, Evanseck J, Field M . All-atom empirical potential for molecular modeling and dynamics studies of proteins. J Phys Chem B. 2014; 102(18):3586-616. DOI: 10.1021/jp973084f. View

Garcia-Sosa A, Hetenyi C, Maran U . Drug efficiency indices for improvement of molecular docking scoring functions. J Comput Chem. 2009; 31(1):174-84. DOI: 10.1002/jcc.21306. View

Deng Y, Roux B . Computations of standard binding free energies with molecular dynamics simulations. J Phys Chem B. 2009; 113(8):2234-46. PMC: 3837708. DOI: 10.1021/jp807701h. View

Kramer B, Rarey M, Lengauer T . Evaluation of the FLEXX incremental construction algorithm for protein-ligand docking. Proteins. 1999; 37(2):228-41. DOI: 10.1002/(sici)1097-0134(19991101)37:2<228::aid-prot8>3.0.co;2-8. View

Tomlinson S, Malmstrom R, Watowich S . New approaches to structure-based discovery of dengue protease inhibitors. Infect Disord Drug Targets. 2009; 9(3):327-43. DOI: 10.2174/1871526510909030327. View

10.

Aqvist J, Marelius J . The linear interaction energy method for predicting ligand binding free energies. Comb Chem High Throughput Screen. 2002; 4(8):613-26. DOI: 10.2174/1386207013330661. View

11.

Mobley D, Chodera J, Dill K . On the use of orientational restraints and symmetry corrections in alchemical free energy calculations. J Chem Phys. 2006; 125(8):084902. PMC: 3583553. DOI: 10.1063/1.2221683. View

12.

Kang N, Lee G, Kim C, Bae M, Kim I, Cho Y . Identification of small molecules that inhibit GSK-3beta through virtual screening. Bioorg Med Chem Lett. 2008; 19(2):533-7. DOI: 10.1016/j.bmcl.2008.10.120. View

13.

Lee K, Kim J, Jeong K, Lee K, Lee Y, Song J . Structure-based virtual screening of Src kinase inhibitors. Bioorg Med Chem. 2009; 17(8):3152-61. DOI: 10.1016/j.bmc.2009.02.054. View

14.

Moustakas D, Lang P, Pegg S, Pettersen E, Kuntz I, Brooijmans N . Development and validation of a modular, extensible docking program: DOCK 5. J Comput Aided Mol Des. 2006; 20(10-11):601-19. DOI: 10.1007/s10822-006-9060-4. View

15.

Wang J, Wolf R, Caldwell J, Kollman P, Case D . Development and testing of a general amber force field. J Comput Chem. 2004; 25(9):1157-74. DOI: 10.1002/jcc.20035. View

16.

Deng Y, Roux B . Calculation of Standard Binding Free Energies: Aromatic Molecules in the T4 Lysozyme L99A Mutant. J Chem Theory Comput. 2015; 2(5):1255-73. DOI: 10.1021/ct060037v. View

17.

MacKerell Jr A, Feig M, Brooks 3rd C . Extending the treatment of backbone energetics in protein force fields: limitations of gas-phase quantum mechanics in reproducing protein conformational distributions in molecular dynamics simulations. J Comput Chem. 2004; 25(11):1400-15. DOI: 10.1002/jcc.20065. View

18.

Charifson P, Corkery J, Murcko M, Walters W . Consensus scoring: A method for obtaining improved hit rates from docking databases of three-dimensional structures into proteins. J Med Chem. 1999; 42(25):5100-9. DOI: 10.1021/jm990352k. View

19.

Mobley D, Graves A, Chodera J, McReynolds A, Shoichet B, Dill K . Predicting absolute ligand binding free energies to a simple model site. J Mol Biol. 2007; 371(4):1118-34. PMC: 2104542. DOI: 10.1016/j.jmb.2007.06.002. View

20.

Deng Y, Roux B . Computation of binding free energy with molecular dynamics and grand canonical Monte Carlo simulations. J Chem Phys. 2008; 128(11):115103. DOI: 10.1063/1.2842080. View