Toward canonical ensemble distribution from self-guided Langevin dynamics simulation.

This work derives a quantitative description of the conformational distribution in self-guided Langevin dynamics (SGLD) simulations. SGLD simulations employ guiding forces calculated from local average momentums to enhance low-frequency motion. This enhancement in low-frequency motion dramatically accelerates conformational search efficiency, but also induces certain perturbations in conformational distribution. Through the local averaging, we separate properties of molecular systems into low-frequency and high-frequency portions. The guiding force effect on the conformational distribution is quantitatively described using these low-frequency and high-frequency properties. This quantitative relation provides a way to convert between a canonical ensemble and a self-guided ensemble. Using example systems, we demonstrated how to utilize the relation to obtain canonical ensemble properties and conformational distributions from SGLD simulations. This development makes SGLD not only an efficient approach for conformational searching, but also an accurate means for conformational sampling.

[1]  Xiongwu Wu,et al.  Enhancing systematic motion in molecular dynamics simulation , 1999 .

[2]  Bernard R. Brooks,et al.  Open Science Grid Study of the Coupling between Conformation and Water Content in the Interior of a Protein , 2008, J. Chem. Inf. Model..

[3]  M. Karplus,et al.  Self-guided enhanced sampling methods for thermodynamic averages , 2003 .

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

[5]  Mark A Olson,et al.  Protein Folding Simulations Combining Self-Guided Langevin Dynamics and Temperature-Based Replica Exchange. , 2010, Journal of chemical theory and computation.

[6]  Shen-Shu Sung,et al.  Constraint dynamics algorithm for simulation of semiflexible macromolecules , 1998 .

[7]  Bernard R Brooks,et al.  Backbone relaxation coupled to the ionization of internal groups in proteins: a self-guided Langevin dynamics study. , 2008, Biophysical journal.

[8]  Jianpeng Ma,et al.  CHARMM: The biomolecular simulation program , 2009, J. Comput. Chem..

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

[10]  Bernard R Brooks,et al.  Isotropic periodic sum: a method for the calculation of long-range interactions. , 2005, The Journal of chemical physics.

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

[12]  Xiongwu Wu,et al.  Self-Guided Molecular Dynamics Simulation for Efficient Conformational Search , 1998 .

[13]  Donald Hamelberg,et al.  A statistical analysis of the precision of reweighting-based simulations. , 2008, The Journal of chemical physics.

[14]  Cheng-I Lee,et al.  Characterizing the denatured state of human prion 121-230. , 2010, Biophysical chemistry.