Alchemical free energy calculations and multiple conformational substates.

Thermodynamic integration (TI) was combined with (adaptive) umbrella sampling to improve the convergence of alchemical free energy simulations in which multiple conformational substates are present. The approach, which we refer to as non-Boltzmann TI (NBTI), was tested by computing the free energy differences between three five-atomic model systems, as well as the free energy difference of solvation between leucine and asparagine. In both cases regular TI failed to give converged results, whereas the NBTI results were free from hysteresis and had standard deviations well below +/-0.7 kcal/mole. We also present theoretical considerations that make it possible to compute free energy differences between simple molecules, such as the five-atomic model systems, by numerical integration of the partition functions at the respective end points.

[1]  T. Straatsma,et al.  Separation‐shifted scaling, a new scaling method for Lennard‐Jones interactions in thermodynamic integration , 1994 .

[2]  James Andrew McCammon,et al.  Ligand-receptor interactions , 1984, Comput. Chem..

[3]  Michael R. Shirts,et al.  Extremely precise free energy calculations of amino acid side chain analogs: Comparison of common molecular mechanics force fields for proteins , 2003 .

[4]  Kentaro Shimizu,et al.  Molecular Dynamics Simulations in Aqueous Solution: Application to Free Energy Calculation of Oligopeptides , 1998 .

[5]  Wei Yang,et al.  Free energy simulations: use of reverse cumulative averaging to determine the equilibrated region and the time required for convergence. , 2004, The Journal of chemical physics.

[6]  Overcoming barriers in macromolecular simulations: non-Boltzmann thermodynamic integration , 1998 .

[7]  W. Im,et al.  Generalized solvent boundary potential for computer simulations , 2001 .

[8]  Y. Sugita,et al.  Replica-exchange molecular dynamics method for protein folding , 1999 .

[9]  Y. Okamoto,et al.  Molecular dynamics, Langevin, and hybrid Monte Carlo simulations in multicanonical ensemble , 1996, physics/9710018.

[10]  William L. Jorgensen,et al.  Free Energies of Hydration and Pure Liquid Properties of Hydrocarbons from the OPLS All-Atom Model , 1994 .

[11]  P. Kollman.  Advances and Continuing Challenges in Achieving Realistic and Predictive Simulations of the Properties of Organic and Biological Molecules , 1996 .

[12]  M. Karplus,et al.  Molecular recognition in proteins. Simulation analysis of substrate binding by a tyrosyl-tRNA synthetase mutant. , 1994, Journal of molecular biology.

[13]  M. Gilson,et al.  The statistical-thermodynamic basis for computation of binding affinities: a critical review. , 1997, Biophysical journal.

[14]  Peter A. Kollman,et al.  The overlooked bond‐stretching contribution in free energy perturbation calculations , 1991 .

[15]  B. Berg,et al.  Multicanonical algorithms for first order phase transitions , 1991 .

[16]  Charles L. Brooks,et al.  λ‐dynamics: A new approach to free energy calculations , 1996 .

[17]  Christopher J. Woods,et al.  Fluoride-selective binding in a new deep cavity calix[4]pyrrole: experiment and theory. , 2002, Journal of the American Chemical Society.

[18]  G. Torrie,et al.  Nonphysical sampling distributions in Monte Carlo free-energy estimation: Umbrella sampling , 1977 .

[19]  Benoît Roux,et al.  Extension to the weighted histogram analysis method: combining umbrella sampling with free energy calculations , 2001 .

[20]  S J Wodak,et al.  Contribution of the hydrophobic effect to protein stability: analysis based on simulations of the Ile-96----Ala mutation in barnase. , 1991, Proceedings of the National Academy of Sciences of the United States of America.

[21]  T. P. Straatsma,et al.  Treatment of rotational isomers in free energy evaluations. Analysis of the evaluation of free energy differences by molecular dynamics simulations of systems with rotational isomeric states , 1989 .

[22]  D. Beveridge,et al.  Free energy via molecular simulation: applications to chemical and biomolecular systems. , 1989, Annual review of biophysics and biophysical chemistry.

[23]  T. P. Straatsma,et al.  Treatment of rotational isomers in free energy calculations. II. Molecular dynamics simulation study of 18‐crown‐6 in aqueous solution as an example of systems with large numbers of rotational isomeric states , 1989 .

[24]  P. Weiner,et al.  Computer Simulation of Biomolecular Systems , 1997 .

[25]  Vijay S. Pande,et al.  Screen Savers of the World Unite! , 2000, Science.

[26]  Peter A. Kollman,et al.  FREE ENERGY CALCULATIONS : APPLICATIONS TO CHEMICAL AND BIOCHEMICAL PHENOMENA , 1993 .

[27]  Stefan Boresch,et al.  Absolute Binding Free Energies: A Quantitative Approach for Their Calculation , 2003 .

[28]  Christopher J. Woods,et al.  Enhanced configurational sampling in binding free-energy calculations , 2003 .

[29]  P M Cullis,et al.  Affinities of amino acid side chains for solvent water. , 1981, Biochemistry.

[30]  W. L. Jorgensen,et al.  Comparison of simple potential functions for simulating liquid water , 1983 .

[31]  Thomas Simonson,et al.  Continuum Treatment of Long-Range Interactions in Free Energy Calculations. Application to Protein−Ligand Binding. , 1997 .

[32]  Peter A. Kollman,et al.  Solvation Free Energies of the Nucleic Acid Bases , 1996 .

[33]  G. Ciccotti,et al.  Numerical Integration of the Cartesian Equations of Motion of a System with Constraints: Molecular Dynamics of n-Alkanes , 1977 .

[34]  T. P. Straatsma,et al.  Treatment of rotational isomeric states. III. The use of biasing potentials , 1994 .

[35]  Stefan Boresch,et al.  The Role of Bonded Energy Terms in Free Energy Simulations - Insights from Analytical Results , 2002 .

[36]  D. Herschbach,et al.  Molecular Partition Functions in Terms of Local Properties , 1959 .

[37]  Benoît Roux,et al.  Free Energy Simulations: Thermodynamic Reversibility and Variability , 2000 .

[38]  A. Mark,et al.  Avoiding singularities and numerical instabilities in free energy calculations based on molecular simulations , 1994 .

[39]  Thomas Simonson,et al.  Free energy simulations come of age: protein-ligand recognition. , 2002, Accounts of chemical research.

[40]  Thomas Simonson,et al.  Free energy of particle insertion , 1993 .

[41]  Alessandra Villa,et al.  Sampling and convergence in free energy calculations of protein-ligand interactions: The binding of triphenoxypyridine derivatives to factor Xa and trypsin , 2003, J. Comput. Aided Mol. Des..

[42]  Lennart Nilsson,et al.  Effects of base substitutions in an RNA hairpin from molecular dynamics and free energy simulations. , 2003, Biophysical journal.

[43]  Wei Yang,et al.  Generalized ensembles serve to improve the convergence of free energy simulations , 2003 .

[44]  J. Mongan,et al.  Accelerated molecular dynamics: a promising and efficient simulation method for biomolecules. , 2004, The Journal of chemical physics.

[45]  M Karplus,et al.  Binding free energies and free energy components from molecular dynamics and Poisson-Boltzmann calculations. Application to amino acid recognition by aspartyl-tRNA synthetase. , 2001, Journal of molecular biology.

[46]  Jonathan W. Essex,et al.  The development of replica-exchange-based free-energy methods , 2003 .

[47]  Grubmüller,et al.  Predicting slow structural transitions in macromolecular systems: Conformational flooding. , 1995, Physical review. E, Statistical physics, plasmas, fluids, and related interdisciplinary topics.

[48]  M. Karplus,et al.  Simulation of activation free energies in molecular systems , 1996 .

[49]  A T Brünger,et al.  Proline cis‐trans isomerization in staphylococcal nuclease: Multi‐substate free energy perturbation calculations , 1995, Protein science : a publication of the Protein Society.

[50]  B. Brooks,et al.  Self-guided Langevin dynamics simulation method , 2003 .

[51]  M Mezei,et al.  Free Energy Simulations a , 1986, Annals of the New York Academy of Sciences.

[52]  David S. Saxon,et al.  Diffuse Surface Optical Model for Nucleon-Nuclei Scattering , 1954 .

[53]  U. Hansmann.  Parallel tempering algorithm for conformational studies of biological molecules , 1997, physics/9710041.

[54]  Stefan Boresch,et al.  THE ROLE OF BONDED TERMS IN FREE ENERGY SIMULATIONS : 1. THEORETICAL ANALYSIS , 1999 .

[55]  Hoover,et al.  Canonical dynamics: Equilibrium phase-space distributions. , 1985, Physical review. A, General physics.

[56]  M L Lamb,et al.  Investigations of neurotrophic inhibitors of FK506 binding protein via Monte Carlo simulations. , 1998, Journal of medicinal chemistry.

[57]  Martin Karplus,et al.  Probability Distributions for Complex Systems: Adaptive Umbrella Sampling of the Potential Energy , 1998 .

[58]  T. Simonson.  Electrostatic Free Energy Calculations for Macromolecules: A Hybrid Molecular Dynamics/Continuum Electrostatics Approach , 2000 .

[59]  Y. Sugita,et al.  Multidimensional replica-exchange method for free-energy calculations , 2000, cond-mat/0009120.

[60]  Peter A. Kollman,et al.  Use of Locally Enhanced Sampling in Free Energy Calculations: Testing and Application to the α → β Anomerization of Glucose , 1998 .

[61]  C. Brooks,et al.  Conformational flexibility in free energy simulations , 1989 .

[62]  Alexander D. MacKerell,et al.  All-atom empirical potential for molecular modeling and dynamics studies of proteins. , 1998, The journal of physical chemistry. B.

[63]  Benoît Roux,et al.  Hydration of Amino Acid Side Chains: Nonpolar and Electrostatic Contributions Calculated from Staged Molecular Dynamics Free Energy Simulations with Explicit Water Molecules , 2004 .

[64]  Stefan Boresch,et al.  The Role of Bonded Terms in Free Energy Simulations. 2. Calculation of Their Influence on Free Energy Differences of Solvation , 1999 .