On the Use of Enveloping Distribution Sampling (EDS) to Compute Free Enthalpy Differences between Different Conformational States of Molecules: Application to 310-, α-, and π-Helices.

Enveloping distribution sampling (EDS) is a powerful method to compute relative free energies from simulation. So far, the EDS method has only been applied to alchemical free energy differences, i.e., between different Hamiltonians defining different systems, and not yet to obtain free energy differences between different conformations or conformational states of a system. In this article, we extend the EDS formalism such that it can be applied to compute free energy differences of different conformations and apply it to compute the relative free enthalpy ΔG of 310-, α-, and π-helices of an alanine deca-peptide in explicit water solvent. The resulting ΔG values are compared to those obtained by standard thermodynamic integration (TI) and from so-called end-state simulations. A TI simulation requires the definition of a λ-dependent pathway which in the present case is based on hydrogen bonds of the different helical conformations. The values of ⟨(∂VTI)/(∂λ)⟩λ show a sharp change for a particular range of λ values, which is indicative of an energy barrier along the pathway, which lowers the accuracy of the resulting ΔG value. In contrast, in a two-state EDS simulation, an unphysical reference-state Hamiltonian which connects the parts of conformational space that are relevant to the different end states is constructed automatically; that is, no pathway needs to be defined. In the simulation using this reference state, both helices were sampled, and many transitions between them occurred, thus ensuring the accuracy of the resulting free enthalpy difference. According to the EDS simulations, the free enthalpy differences of the π-helix and the 310-helix versus the α-helix are 5 kJ mol(-1) and 47 kJ mol(-1), respectively, for an alanine deca-peptide in explicit SPC water solvent using the GROMOS 53A6 force field. The EDS method, which is a particular form of umbrella sampling, is thus applicable to compute free energy differences between conformational states as well as between systems and has definite advantages over the traditional TI and umbrella sampling methods to compute relative free energies.

[1]  Wilfred F van Gunsteren,et al.  Simple, Efficient, and Reliable Computation of Multiple Free Energy Differences from a Single Simulation: A Reference Hamiltonian Parameter Update Scheme for Enveloping Distribution Sampling (EDS). , 2009, Journal of chemical theory and computation.

[2]  H. Meirovitch.  Recent developments in methodologies for calculating the entropy and free energy of biological systems by computer simulation. , 2007, Current opinion in structural biology.

[3]  H. Berendsen,et al.  Molecular dynamics with coupling to an external bath , 1984 .

[4]  Charles L Brooks,et al.  Protein and peptide folding explored with molecular simulations. , 2002, Accounts of chemical research.

[5]  Wilfred F. van Gunsteren,et al.  A generalized reaction field method for molecular dynamics simulations , 1995 .

[6]  Jun Wang,et al.  Determination of conformational free energies of peptides by multidimensional adaptive umbrella sampling. , 2006, The Journal of chemical physics.

[7]  Charles L. Brooks,et al.  λ‐Dynamics free energy simulation methods , 2009, J. Comput. Chem..

[8]  M. Gilson,et al.  Calculation of protein-ligand binding affinities. , 2007, Annual review of biophysics and biomolecular structure.

[9]  Markus Christen,et al.  The GROMOS software for biomolecular simulation: GROMOS05 , 2005, J. Comput. Chem..

[10]  D. Kofke.  Free energy methods in molecular simulation , 2005 .

[11]  Wilfred F van Gunsteren,et al.  Free energies of ligand binding for structurally diverse compounds. , 2005, Proceedings of the National Academy of Sciences of the United States of America.

[12]  J. G. Powles,et al.  Non-destructive molecular-dynamics simulation of the chemical potential of a fluid , 1982 .

[13]  J. Kirkwood.  Statistical Mechanics of Fluid Mixtures , 1935 .

[14]  X. Daura,et al.  Reversible peptide folding in solution by molecular dynamics simulation. , 1998, Journal of molecular biology.

[15]  Philippe H Hünenberger,et al.  Ball-and-Stick Local Elevation Umbrella Sampling: Molecular Simulations Involving Enhanced Sampling within Conformational or Alchemical Subspaces of Low Internal Dimensionalities, Minimal Irrelevant Volumes, and Problem-Adapted Geometries. , 2010, Journal of chemical theory and computation.

[16]  W. L. Jorgensen.  The Many Roles of Computation in Drug Discovery , 2004, Science.

[17]  Stefan Bruckner,et al.  Unorthodox uses of Bennett's acceptance ratio method , 2009, J. Comput. Chem..

[18]  Thomas B. Woolf,et al.  Improving the efficiency and reliability of free energy perturbation calculations using overlap sampling methods , 2004, J. Comput. Chem..

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

[20]  Wilfred F. van Gunsteren,et al.  Basic ingredients of free energy calculations: A review , 2009, J. Comput. Chem..

[21]  Charles H. Bennett,et al.  Efficient estimation of free energy differences from Monte Carlo data , 1976 .

[22]  Wilfred F van Gunsteren,et al.  Comparison of enveloping distribution sampling and thermodynamic integration to calculate binding free energies of phenylethanolamine N-methyltransferase inhibitors. , 2011, The Journal of chemical physics.

[23]  Alan E. Mark,et al.  Estimating the Relative Free Energy of Different Molecular States with Respect to a Single Reference State , 1996 .

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

[25]  Pengyu Y. Ren,et al.  Ion solvation thermodynamics from simulation with a polarizable force field. , 2003, Journal of the American Chemical Society.

[26]  Wilfred F. van Gunsteren,et al.  Computation of free energy , 2002 .

[27]  H. Berendsen,et al.  THERMODYNAMICS OF CAVITY FORMATION IN WATER - A MOLECULAR-DYNAMICS STUDY , 1982 .

[28]  W. V. van Gunsteren,et al.  Prediction of folding equilibria of differently substituted peptides using one-step perturbation. , 2010, Journal of the American Chemical Society.

[29]  William L Jorgensen,et al.  Perspective on Free-Energy Perturbation Calculations for Chemical Equilibria. , 2008, Journal of chemical theory and computation.

[30]  Wilfred F van Gunsteren,et al.  Enveloping distribution sampling: a method to calculate free energy differences from a single simulation. , 2007, The Journal of chemical physics.

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

[32]  M Karplus,et al.  Ion transport in a model gramicidin channel. Structure and thermodynamics. , 1991, Biophysical journal.

[33]  Eric Vanden-Eijnden,et al.  Some recent techniques for free energy calculations , 2009, J. Comput. Chem..

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

[35]  R. Zwanzig.  High‐Temperature Equation of State by a Perturbation Method. I. Nonpolar Gases , 1954 .

[36]  Krzysztof Kuczera,et al.  Conformational Free Energy Surfaces of Ala10 and Aib10 Peptide Helices in Solution , 2001 .

[37]  Régis Pomès,et al.  Enhancing the accuracy, the efficiency and the scope of free energy simulations. , 2005, Current opinion in structural biology.

[38]  K. Shing,et al.  The chemical potential in dense fluids and fluid mixtures via computer simulation , 1982 .

[39]  Wilfred F. van Gunsteren,et al.  Comparison of three enveloping distribution sampling Hamiltonians for the estimation of multiple free energy differences from a single simulation , 2009, J. Comput. Chem..

[40]  Philippe H. Hünenberger,et al.  Using the local elevation method to construct optimized umbrella sampling potentials: Calculation of the relative free energies and interconversion barriers of glucopyranose ring conformers in water , 2010, J. Comput. Chem..

[41]  Wilfred F van Gunsteren,et al.  Multiple free energies from a single simulation: extending enveloping distribution sampling to nonoverlapping phase-space distributions. , 2008, The Journal of chemical physics.

[42]  Chris Oostenbrink,et al.  A biomolecular force field based on the free enthalpy of hydration and solvation: The GROMOS force‐field parameter sets 53A5 and 53A6 , 2004, J. Comput. Chem..

[43]  David A. Kofke,et al.  Appropriate methods to combine forward and reverse free-energy perturbation averages , 2003 .