Polymer dynamics in torsion space

The large scale motion of proteins, or covalently bonded polymers in general, is governed by the dynamics of the torsion angles, with bond lengths and bond angles kept approximately constant. In the present work, the Lagrangian equations of torsion motion are derived for a general macromolecule. The dynamics is implemented numerically for a test protein, using the velocity Verlet method as the integrator. The results indicate time steps of up to about 30 fs can be used for short time (up to at least 20 ps) simulations, before the dynamics and energy start to differ significantly from results obtained with smaller time steps. For longer time simulations, up to 1000 ps, a time step of 10 fs is relatively safe. © 2010 Wiley Periodicals, Inc. J Comput Chem, 2010

[1]  H. Scheraga,et al.  Variable step molecular dynamics: An exploratory technique for peptides with fixed geometry , 1990 .

[2]  M. Fixman.  Simulation of polymer dynamics. II. Relaxation rates and dynamic viscosity , 1978 .

[3]  Arnold Neumaier,et al.  Molecular Modeling of Proteins and Mathematical Prediction of Protein Structure , 1997, SIAM Rev..

[4]  J. H. Weiner,et al.  Brownian dynamics study of a polymer chain of linked rigid bodies , 1979 .

[5]  Abhinandan Jain,et al.  A fast recursive algorithm for molecular dynamics simulation , 1993 .

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

[7]  Harold A. Scheraga,et al.  Analysis of the Contribution of Internal Vibrations to the Statistical Weights of Equilibrium Conformations of Macromolecules , 1969 .

[8]  P. Chacón,et al.  Thorough validation of protein normal mode analysis: a comparative study with essential dynamics. , 2007, Structure.

[9]  Kai Nordlund,et al.  Analytical interatomic potential for modeling nonequilibrium processes in the W–C–H system , 2005 .

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

[11]  Ilpo Vattulainen,et al.  Multiscale modeling of emergent materials: biological and soft matter. , 2009, Physical chemistry chemical physics : PCCP.

[12]  A. Mazur.  Quasi‐Hamiltonian equations of motion for internal coordinate molecular dynamics of polymers , 1997 .

[13]  Tirion,et al.  Large Amplitude Elastic Motions in Proteins from a Single-Parameter, Atomic Analysis. , 1996, Physical review letters.

[14]  D. Brenner,et al.  Empirical potential for hydrocarbons for use in simulating the chemical vapor deposition of diamond films. , 1990, Physical review. B, Condensed matter.

[15]  W. F. van Gunsteren,et al.  Effect of constraints on the dynamics of macromolecules , 1982 .

[16]  Pablo Echenique,et al.  Quantum mechanical calculation of the effects of stiff and rigid constraints in the conformational equilibrium of the alanine dipeptide , 2006, J. Comput. Chem..

[17]  J. Pesonen.  Vibrational coordinates and their gradients: A geometric algebra approach , 2000 .

[18]  Brenner.  Erratum: Empirical potential for hydrocarbons for use in simulating the chemical vapor deposition of diamond films , 1992, Physical review. B, Condensed matter.

[19]  Wilfred F van Gunsteren,et al.  Biomolecular modeling: Goals, problems, perspectives. , 2006, Angewandte Chemie.

[20]  H. Scheraga,et al.  Macromolecular conformational dynamics in torsional angle space , 1998 .

[21]  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..

[22]  M. Fixman,et al.  Simulation of polymer dynamics. I. General theory , 1978 .

[23]  Kai Nordlund,et al.  Simulations of cementite: An analytical potential for the Fe-C system , 2009 .

[24]  Tosiyuki Noguti,et al.  Dynamics of Native Globular Proteins in Terms of Dihedral Angles , 1983 .