Verifying trivial parallelization of multicanonical molecular dynamics for conformational sampling of a polypeptide in explicit water

Multicanonical molecular dynamics (McMD) simulation is an enhanced conformational sampling method. For more accelerated sampling, we examined a trajectory-parallelization method, which trivially integrates multiple McMD trajectories into one. This method decreases communications between computing nodes largely. We applied this method to conformational sampling of a seven-residue peptide in explicit water. Physical quantities converged well with those from a long McMD simulation in a wide temperature range. Thus, the integration can generate a canonical ensemble without performing the long simulation. The current method is useful for folding study of long polypeptides (proteins) in explicit water with massive parallelism.

[1]  Gary P. Morriss,et al.  The isothermal/isobaric molecular dynamics ensemble , 1983 .

[2]  Haruki Nakamura,et al.  Energy landscape of a β-hairpin peptide in explicit water studied by multicanonical molecular dynamics , 2001 .

[3]  D. Teplow,et al.  Kinetic Studies of Amyloid β-Protein Fibril Assembly , 2002, The Journal of Biological Chemistry.

[4]  Junichi Higo,et al.  AMBER-based hybrid force field for conformational sampling of polypeptides , 2005 .

[5]  Haruki Nakamura,et al.  Conformational sampling of a 40-residue protein consisting of α and β secondary-structure elements in explicit solvent , 2007 .

[6]  Haruki Nakamura,et al.  Energy landscape of a peptide consisting of α‐helix, 310‐helix, β‐turn, β‐hairpin, and other disordered conformations , 2001 .

[7]  B. Berg,et al.  Multicanonical algorithms for first order phase transitions , 1991 .

[8]  Joan-Emma Shea,et al.  Folding of the 25 residue Abeta(12-36) peptide in TFE/water: temperature-dependent transition from a funneled free-energy landscape to a rugged one. , 2007, The journal of physical chemistry. B.

[9]  Haruki Nakamura,et al.  Protein‐inhibitor flexible docking by a multicanonical sampling: Native complex structure with the lowest free energy and a free‐energy barrier distinguishing the native complex from the others , 2007, Proteins.

[10]  A. Kidera,et al.  Multicanonical Ensemble Generated by Molecular Dynamics Simulation for Enhanced Conformational Sampling of Peptides , 1997 .

[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]  T. Hirokawa,et al.  Protein-segment universe exhibiting transitions at intermediate segment length in conformational subspaces , 2008, BMC Structural Biology.

[13]  K. Hukushima,et al.  Exchange Monte Carlo Method and Application to Spin Glass Simulations , 1995, cond-mat/9512035.

[14]  Haruki Nakamura,et al.  Conformational transition states of a β‐hairpin peptide between the ordered and disordered conformations in explicit water , 2002, Protein science : a publication of the Protein Society.

[15]  W. L. Jorgensen,et al.  Comparison of simple potential functions for simulating liquid water , 1983 .

[16]  A. Pohorille,et al.  The development/application of a ‘minimalist’ organic/biochemical molecular mechanic force field using a combination of ab initio calculations and experimental data , 1997 .

[17]  W. Goddard,et al.  Atomic level simulations on a million particles: The cell multipole method for Coulomb and London nonbond interactions , 1992 .

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

[19]  Haruki Nakamura,et al.  Presto(protein Engineering Simulator): A Vectorized Molecular Mechanics Program for Biopolymers , 1992, Comput. Chem..

[20]  Haruki Nakamura,et al.  Flexible docking of a ligand peptide to a receptor protein by multicanonical molecular dynamics simulation , 1997 .