Solvent Interaction Energy Calculations on Molecular Dynamics Trajectories: Increasing the Efficiency Using Systematic Frame Selection
- 28 August 2011
- journal article
- research article
- Published by American Chemical Society (ACS) in Journal of Chemical Information and Modeling
- Vol. 51 (10), 2680-2689
- https://doi.org/10.1021/ci200191m
Abstract
End-point methods such as linear interaction energy (LIE) analysis, molecular mechanics generalized Born solvent-accessible surface (MM/GBSA), and solvent interaction energy (SIE) analysis have become popular techniques to calculate the free energy associated with protein–ligand binding. Such methods typically use molecular dynamics (MD) simulations to generate an ensemble of protein structures that encompasses the bound and unbound states. The energy evaluation method (LIE, MM/GBSA, or SIE) is subsequently used to calculate the energy of each member of the ensemble, thus providing an estimate of the average free energy difference between the bound and unbound states. The workflow requiring both MD simulation and energy calculation for each frame and each trajectory proves to be computationally expensive. In an attempt to reduce the high computational cost associated with end-point methods, we study several methods by which frames may be intelligently selected from the MD simulation including clustering and address the question of how the number of selected frames influences the accuracy of the SIE calculations.Keywords
This publication has 45 references indexed in Scilit:
- Computational Determination of Binding Structures and Free Energies of Phosphodiesterase-2 with Benzo[1,4]diazepin-2-one DerivativesThe Journal of Physical Chemistry B, 2010
- Modeling Protein−Ligand Binding by Mining MinimaJournal of Chemical Theory and Computation, 2010
- Absolute binding free energy calculations: On the accuracy of computational scoring of protein–ligand interactionsProteins-Structure Function and Bioinformatics, 2010
- Predictions of Binding for Dopamine D2 Receptor Antagonists by the SIE MethodJournal of Chemical Information and Modeling, 2009
- Development and testing of a general amber force fieldJournal of Computational Chemistry, 2004
- A point‐charge force field for molecular mechanics simulations of proteins based on condensed‐phase quantum mechanical calculationsJournal of Computational Chemistry, 2003
- Fast, efficient generation of high‐quality atomic charges. AM1‐BCC model: II. Parameterization and validationJournal of Computational Chemistry, 2002
- Enhanced docking with the mining minima optimizer: Acceleration and side-chain flexibilityJournal of Computational Chemistry, 2002
- Asparagine and glutamine: using hydrogen atom contacts in the choice of side-chain amide orientation 1 1Edited by J. ThorntonJournal of Molecular Biology, 1999
- Design of Novel, Potent, Noncovalent Inhibitors of Thrombin with Nonbasic P-1 Substructures: Rapid Structure−Activity Studies by Solid-Phase SynthesisJournal of Medicinal Chemistry, 1998