Simulations of proteins with inhomogeneous degrees of freedom: The effect of thermostats

Proteins with a long flexible polymeric tail attached at their N‐ or C‐terminus are studied using molecular dynamics (MD) simulations of a coarse‐grained model for protein folding where the temperature is regulated by either the Berendsen or the Langevin thermostat. These thermostats show different abilities to regulate the temperature of these systems that include flexible and more rigid regions. In the simulations with the Berendsen thermostat, the flexible tail is significantly hotter than the protein, both in its folded and unfolded states. Upon weakening the strength of the Berendsen thermostat, the temperature gradient between the fast and the slow degrees of freedom is significantly decreased, yet linkage between the temperatures of the flexible tail and the protein remains. The Langevin thermostat is proven to regulate the temperature of these inhomogenous systems reliably, without discriminating between the slow and fast degrees of freedom. The Langevin thermostat is less sensitive than is the Berendsen thermostat to the strength of the coupling between the protein system and the thermal bath. Our study calls for special care in choosing the thermostat for MD simulations of systems with inhomogenous degrees of freedom. Using the Berendsen thermostat with strong coupling would result in mistaken thermodynamic descriptions of such systems. © 2008 Wiley Periodicals, Inc. J Comput Chem 2008

[1]  J. Onuchic,et al.  Equilibrium and kinetic folding pathways of a TIM barrel with a funneled energy landscape. , 2005, Biophysical journal.

[2]  M. Karplus,et al.  Advances in chemical physics, volume 71: Proteins: A theoretical perspective of dynamics, structure, and thermodynamics , 2006 .

[3]  H. Dyson,et al.  Coupling of folding and binding for unstructured proteins. , 2002, Current opinion in structural biology.

[4]  Kenneth M. Merz,et al.  Application of the Nosé−Hoover Chain Algorithm to the Study of Protein Dynamics , 1996 .

[5]  Effect of the Berendsen thermostat on the dynamical properties of water , 2004, cond-mat/0403458.

[6]  Berend Smit,et al.  Understanding molecular simulation: from algorithms to applications , 1996 .

[7]  J. Onuchic,et al.  Topological and energetic factors: what determines the structural details of the transition state ensemble and "en-route" intermediates for protein folding? An investigation for small globular proteins. , 2000, Journal of molecular biology.

[8]  V. Uversky.  Intrinsically Disordered Proteins , 2000 .

[9]  Tetsuya Morishita,et al.  Fluctuation formulas in molecular-dynamics simulations with the weak coupling heat bath , 2000 .

[10]  S. Plotkin.  Speeding protein folding beyond the Gō model: How a little frustration sometimes helps , 2000, Proteins.

[11]  H. Dyson,et al.  Intrinsically unstructured proteins and their functions , 2005, Nature Reviews Molecular Cell Biology.

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

[13]  Shoji Takada,et al.  How protein thermodynamics and folding mechanisms are altered by the chaperonin cage: Molecular simulations , 2003, Proceedings of the National Academy of Sciences of the United States of America.

[14]  John Karanicolas,et al.  Improved Gō-like models demonstrate the robustness of protein folding mechanisms towards non-native interactions. , 2003, Journal of molecular biology.

[15]  Peter G Wolynes,et al.  The energy landscape of modular repeat proteins: topology determines folding mechanism in the ankyrin family. , 2005, Journal of molecular biology.

[16]  Cecilia Clementi,et al.  The experimental folding landscape of monomeric lactose repressor, a large two-domain protein, involves two kinetic intermediates. , 2005, Proceedings of the National Academy of Sciences of the United States of America.

[17]  Roland Winter,et al.  The amino-terminal PrP domain is crucial to modulate prion misfolding and aggregation. , 2005, Biophysical journal.

[18]  Andrea Amadei,et al.  Dynamical and Statistical Mechanical Characterization of Temperature Coupling Algorithms , 2002 .

[19]  D. Thirumalai,et al.  Molecular crowding enhances native state stability and refolding rates of globular proteins. , 2005, Proceedings of the National Academy of Sciences of the United States of America.

[20]  R. Swendsen,et al.  THE weighted histogram analysis method for free‐energy calculations on biomolecules. I. The method , 1992 .

[21]  S. Harvey,et al.  The flying ice cube: Velocity rescaling in molecular dynamics leads to violation of energy equipartition , 1998, J. Comput. Chem..

[22]  Peter G Wolynes,et al.  A survey of flexible protein binding mechanisms and their transition states using native topology based energy landscapes. , 2005, Journal of molecular biology.

[23]  B. Brooks,et al.  Constant pressure molecular dynamics simulation: The Langevin piston method , 1995 .

[24]  Tamar Schlick,et al.  Molecular Modeling and Simulation: An Interdisciplinary Guide , 2010 .

[25]  P. Hünenberger.  Thermostat Algorithms for Molecular Dynamics Simulations , 2005 .

[26]  P. Tompa.  Intrinsically unstructured proteins. , 2002, Trends in biochemical sciences.

[27]  Yaakov Levy,et al.  The folding energy landscape of the dimerization domain of Escherichia coli Trp repressor: a joint experimental and theoretical investigation. , 2006, Journal of molecular biology.

[28]  P. Wolynes,et al.  The experimental survey of protein-folding energy landscapes , 2005, Quarterly Reviews of Biophysics.

[29]  Tongye Shen,et al.  φ-Value Analysis of Apo-Azurin Folding: Comparison between Experiment and Theory† , 2006 .

[30]  M. Karplus,et al.  Proteins: A Theoretical Perspective of Dynamics, Structure, and Thermodynamics , 1988 .

[31]  José N Onuchic,et al.  The native energy landscape for interleukin-1beta. Modulation of the population ensemble through native-state topology. , 2005, Journal of molecular biology.