Efficient Computation of Small-Molecule Configurational Binding Entropy and Free Energy Changes by Ensemble Enumeration
Name
Silver-2013-Efficient Computation.pdf
Size
3.18 MB
Format
Adobe PDF
Checksum (MD5)
8925a2ad847085c65e46486b0bd6d353
Author(s) • • • • • • • •
Silver, Nathaniel W.
King, Bracken Matheny
Nalam, Madhavi N. L.
Cao, Hong
Ali, Akbar
Kiran Kumar Reddy, G. S.
Rana, Tariq M.
Schiffer, Celia A.
Tidor, Bruce
Date Issued
November 2013
Journal
Journal of Chemical Theory and Computation
Publisher
American Chemical Society (ACS)
Citation
Silver, Nathaniel W., Bracken M. King, Madhavi N. L. Nalam, Hong Cao, Akbar Ali, G. S. Kiran Kumar Reddy, Tariq M. Rana, Celia A. Schiffer, and Bruce Tidor. “Efficient Computation of Small-Molecule Configurational Binding Entropy and Free Energy Changes by Ensemble Enumeration.” Journal of Chemical Theory and Computation 9, no. 11 (November 12, 2013): 5098–5115. © 2013 American Chemical Society.
Version
Final published version
Abstract
Here we present a novel, end-point method using the dead-end-elimination and A* algorithms to efficiently and accurately calculate the change in free energy, enthalpy, and configurational entropy of binding for ligand–receptor association reactions. We apply the new approach to the binding of a series of human immunodeficiency virus (HIV-1) protease inhibitors to examine the effect ensemble reranking has on relative accuracy as well as to evaluate the role of the absolute and relative ligand configurational entropy losses upon binding in affinity differences for structurally related inhibitors. Our results suggest that most thermodynamic parameters can be estimated using only a small fraction of the full configurational space, and we see significant improvement in relative accuracy when using an ensemble versus single-conformer approach to ligand ranking. We also find that using approximate metrics based on the single-conformation enthalpy differences between the global minimum energy configuration in the bound as well as unbound states also correlates well with experiment. Using a novel, additive entropy expansion based on conditional mutual information, we also analyze the source of ligand configurational entropy loss upon binding in terms of both uncoupled per degree of freedom losses as well as changes in coupling between inhibitor degrees of freedom. We estimate entropic free energy losses of approximately +24 kcal/mol, 12 kcal/mol of which stems from loss of translational and rotational entropy. Coupling effects contribute only a small fraction to the overall entropy change (1–2 kcal/mol) but suggest differences in how inhibitor dihedral angles couple to each other in the bound versus unbound states. The importance of accounting for flexibility in drug optimization and design is also discussed.
MIT Department
Massachusetts Institute of Technology. Computer Science and Artificial Intelligence Laboratory
Massachusetts Institute of Technology. Department of Biological Engineering
Massachusetts Institute of Technology. Department of Chemistry
Massachusetts Institute of Technology. Department of Electrical Engineering and Computer Science
Terms of Use
Article is made available in accordance with the publisher's policy and may be subject to US copyright law. Please refer to the publisher's site for terms of use.
Persistent DSpace Link
DOI of Published Version
https://doi.org/10.1021/ct400383v