Analyzing the robustness of the MM/PBSA free energy calculation method: Application to DNA conformational transitions
Authors:
- 1 February 2011
- journal article
- research article
- Published by Wiley in Journal of Computational Chemistry
- Vol. 32 (7) , 1431-1440
Abstract
The ability to predict and characterize free energy differences associated with conformational equilibria or the binding of biomolecules is vital to understanding the molecular basis of many important biological functions. As biological studies focus on larger molecular complexes and properties of the genome, proteome, and interactome, the development and characterization of efficient methods for calculating free energy becomes increasingly essential. The aim of this study is to examine the robustness of the end‐point free energy method termed the molecular mechanics Poisson‐Boltzmann solvent accessible surface area (MM/PBSA) method. Specifically, applications of MM/PBSA to the conformational equilibria of nucleic acid (NA) systems are explored. This is achieved by comparing A to B form DNA conformational free energy differences calculated using MM/PBSA with corresponding free energy differences determined with a more rigorous and time‐consuming umbrella sampling algorithm. In addition, the robustness of NA MM/PBSA calculations is also evaluated in terms of the sensitivity towards the choice of force field and the choice of solvent model used during conformational sampling. MM/PBSA calculations of the free energy difference between A‐form and B‐form DNA are shown to be in very close agreement with the PMF result determined using an umbrella sampling approach. Further, it is found that the MM/PBSA conformational free energy differences were also in agreement using either the CHARMM or AMBER force field. The influence of ionic strength on conformational stability was particularly insensitive to the choice of force field. Finally, it is also shown that the use of a generalized Born implicit solvent during conformational sampling results in free energy estimates that deviate slightly from those obtained using explicitly solvated MD simulations in these NA systems. © 2011 Wiley Periodicals, Inc. J Comput Chem, 2011Keywords
This publication has 43 references indexed in Scilit:
- Poisson-Boltzmann Calculations of Nonspecific Salt Effects on Protein-Protein Binding Free EnergiesBiophysical Journal, 2007
- Ligand Affinities Predicted with the MM/PBSA Method: Dependence on the Simulation Method and the Force FieldJournal of Medicinal Chemistry, 2006
- Implicit Solvent Simulations of DNA and DNA−Protein Complexes: Agreement with Explicit Solvent vs ExperimentThe Journal of Physical Chemistry B, 2006
- Balancing an accurate representation of the molecular surface in generalized born formalisms with integrator stability in molecular dynamics simulationsJournal of Computational Chemistry, 2006
- The Dynamics of the B−A Transition of Natural DNA Double HelicesJournal of the American Chemical Society, 2005
- DNA Structure: What's in Charge?Journal of Molecular Biology, 2000
- A-form Conformational Motifs in Ligand-bound DNA StructuresJournal of Molecular Biology, 2000
- Calculation and analysis of low frequency normal modes for DNAJournal of Computational Chemistry, 1997
- CHARMM: A program for macromolecular energy, minimization, and dynamics calculationsJournal of Computational Chemistry, 1983
- Polymorphism of DNA double helicesJournal of Molecular Biology, 1980