Premium
Calculation of ligand binding free energies from molecular dynamics simulations
Author(s) -
Marelius J.,
Hansson T.,
Åqvist J.
Publication year - 1998
Publication title -
international journal of quantum chemistry
Language(s) - English
Resource type - Journals
SCImago Journal Rank - 0.484
H-Index - 105
eISSN - 1097-461X
pISSN - 0020-7608
DOI - 10.1002/(sici)1097-461x(1998)69:1<77::aid-qua10>3.0.co;2-2
Subject(s) - free energy perturbation , chemistry , molecular dynamics , thermodynamic integration , computational chemistry , scaling , umbrella sampling , binding energy , ligand (biochemistry) , intermolecular force , chemical physics , molecule , atomic physics , physics , receptor , geometry , mathematics , organic chemistry , biochemistry
A recently developed method for predicting binding affinities in ligand–receptor complexes, based on interaction energy averaging and conformational sampling by molecular dynamics simulation, is presented. Polar and nonpolar contributions to the binding free energy are approximated by a linear scaling of the corresponding terms in the average intermolecular interaction energy for the bound and free states of the ligand. While the method originally assumed the validity of electrostatic linear response, we show that incorporation of systematic deviations from linear response derived from free energy perturbation calculations enhances the accuracy of the approach. The method is applied to complexes of wild‐type and mutant human dihydrofolate reductases with 2,4‐diaminopteridine and 2,4‐diaminoquinazoline inhibitors. It is shown that a binding energy accuracy of about 1 kcal/mol is attainable even for multiply ionized compounds, such as methotrexate, for which electrostatic interactions energies are very large. © 1998 John Wiley & Sons, Inc. Int J Quant Chem 69: 77–88, 1998