Characterization of conformational equilibria through Hamiltonian and temperature replica‐exchange simulations: Assessing entropic and environmental effects

Molecular dynamics simulations based on the replica‐exchange framework (REMD) are emerging as a useful tool to characterize the conformational variability that is intrinsic to most chemical and biological systems. In this work, it is shown that a simple extension of the replica‐exchange method, known as Hamiltonian REMD, greatly facilitates the characterization of conformational equilibria across large energetic barriers, or in the presence of substantial entropic effects, overcoming some of the difficulties of REMD based on temperature alone. In particular, a comparative assessment of the HREMD and TREMD approaches was made, through computation of the gas‐phase free‐energy difference between the so‐called D2d and S4 states of tetrabutylammonium (TBA), an ionic compound of frequently used in biophysical studies of ion channels. Taking advantage of the greater efficiency of the HREMD scheme, the conformational equilibrium of TBA was characterized in a variety of conditions. Simulation of the gas‐phase equilibrium in the 100–300 K range allowed us to compute the entropy difference between these states as well as to describe its temperature dependence. Through HREMD simulations of TBA in a water droplet, the effect of solvation on the conformational equilibrium was determined. Finally, the equilibrium of TBA in the context of a simplified model of the binding cavity of the KcsA potassium channel was simulated, and density maps for D2d and S4 states analogous to those derived from X‐ray crystallography were constructed. Overall, this work illustrates the potential of the HREMD approach in the context of computational drug design, ligand‐receptor structural prediction and more generally, molecular recognition, where one of the most challenging issues remains to account for conformational flexibility as well for the solvation and entropic effects thereon. © 2007 Wiley Periodicals, Inc. J Comput Chem, 2007

[1]  C. Brooks Computer simulation of liquids , 1989 .

[2]  P. Focia,et al.  Structural basis of TEA blockade in a model potassium channel , 2005, Nature Structural &Molecular Biology.

[3]  Michael W Deem,et al.  Parallel tempering: theory, applications, and new perspectives. , 2005, Physical chemistry chemical physics : PCCP.

[4]  Jyoti Chattopadhyaya,et al.  Computational and NMR study of quaternary ammonium ion conformations in solution , 2002 .

[5]  Yuji Sugita,et al.  Replica-exchange multicanonical algorithm and multicanonical replica-exchange method for simulating systems with rough energy landscape , 2000, cond-mat/0009119.

[6]  B. Berne,et al.  Replica exchange with solute tempering: a method for sampling biological systems in explicit water. , 2005, Proceedings of the National Academy of Sciences of the United States of America.

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

[8]  Michael R. Shirts,et al.  Atomistic protein folding simulations on the submillisecond time scale using worldwide distributed computing. , 2003, Biopolymers.

[9]  Edward Lyman,et al.  A Second Look at Canonical Sampling of Biomolecules using Replica Exchange Simulation. , 2006, Journal of chemical theory and computation.

[10]  S. Takada,et al.  On the Hamiltonian replica exchange method for efficient sampling of biomolecular systems: Application to protein structure prediction , 2002 .

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

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

[13]  Benoît Roux,et al.  Extracellular Blockade of K+ Channels by Tea , 2001, The Journal of general physiology.

[14]  David Chandler,et al.  Transition path sampling: throwing ropes over rough mountain passes, in the dark. , 2002, Annual review of physical chemistry.

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

[16]  Dusanka Janezic,et al.  Harmonic analysis of large systems. I. Methodology , 1995, J. Comput. Chem..

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

[18]  I. Tavernelli,et al.  A Novel Hamiltonian Replica Exchange MD Protocol to Enhance Protein Conformational Space Sampling. , 2006, Journal of chemical theory and computation.

[19]  Preston Moore,et al.  Metadynamics as a tool for exploring free energy landscapes of chemical reactions. , 2006, Accounts of chemical research.

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

[21]  M. Karplus,et al.  CHARMM: A program for macromolecular energy, minimization, and dynamics calculations , 1983 .