Journal
JOURNAL OF CHEMICAL INFORMATION AND MODELING
Volume 53, Issue 10, Pages 2659-2670Publisher
AMER CHEMICAL SOC
DOI: 10.1021/ci400312v
Keywords
-
Categories
Funding
- NIH [CA135380, AA0197461]
- INGEN grant from the Lilly Endowment, Inc.
- National Science Foundation
- Shared University Research grants from IBM, Inc.
Ask authors/readers for more resources
End-point free energy calculations using MM-GBSA and MM-PBSA provide a detailed understanding of molecular recognition in protein-ligand interactions. The binding free energy can be used to rank-order protein-ligand structures in virtual screening for compound or target identification. Here, we carry out free energy calculations for a diverse set of 11 proteins bound to 14 small molecules using extensive explicit-solvent MD simulations. The structure of these complexes was previously solved by crystallography and their binding studied with isothermal titration calorimetry (ITC) data enabling direct comparison to the MM-GBSA and MM-PBSA calculations. Four MM-GBSA and three MM-PBSA calculations reproduced the ITC free energy within 1 kcal.mol(-1) highlighting the challenges in reproducing the absolute free energy from end-point free energy calculations. MM-GBSA exhibited better rank-ordering with a Spearman rho of 0.68 compared to 0.40 for MM-PBSA with dielectric constant (epsilon = 1). An increase in epsilon resulted in significantly better rank-ordering for MM-PBSA (rho = 0.91 for epsilon = 10), but larger epsilon significantly reduced the contributions of electrostatics, suggesting that the improvement is due to the nonpolar and entropy components, rather than a better representation of the electrostatics. The SVRKB scoring function applied to MD snapshots resulted in excellent rank-ordering (rho = 0.81). Calculations of the configurational entropy using normal-mode analysis led to free energies that correlated significantly better to the ITC free energy than the MD-based quasi-harmonic approach, but the computed entropies showed no correlation with the ITC entropy. When the adaptation energy is taken into consideration by running separate simulations for complex, apo, and ligand (MM-PBSA(ADAPT)), there is less agreement with the ITC data for the individual free energies, but remarkably good rank-ordering is observed (rho = 0.89). Interestingly, filtering MD snapshots by prescoring protein-ligand complexes with a machine learning-based approach (SVMSP) resulted in a significant improvement in the MM-PBSA results (epsilon = 1) from rho = 0.40 to rho = 0.81. Finally, the nonpolar components of MM-GBSA and MM-PBSA, but not the electrostatic components, showed strong correlation to the ITC free energy; the computed entropies did not correlate with the ITC entropy.
Authors
I am an author on this paper
Click your name to claim this paper and add it to your profile.
Reviews
Recommended
No Data Available