This Article Full Text Full Text (PDF) Alert me when this article is cited Alert me if a correction is posted Services Similar articles in this journal Similar articles in PubMed Alert me to new issues of the journal Download to citation manager Citing Articles Citing Articles via HighWire Citing Articles via Google Scholar Google Scholar Articles by Rhee, Y. M. Articles by Pande, V. S. Search for Related Content PubMed PubMed Citation Articles by Rhee, Y. M. Articles by Pande, V. S.

Multiplexed-Replica Exchange Molecular Dynamics Method for Protein Folding Simulation

Department of Chemistry, Stanford University, Stanford California 94305-5080

Correspondence: Address reprint requests to Vijay S. Pande, E-mail: pande{at}stanford.edu.

Simulating protein folding thermodynamics starting purely from a protein sequence is a grand challenge of computational biology. Here, we present an algorithm to calculate a canonical distribution from molecular dynamics simulation of protein folding. This algorithm is based on the replica exchange method where the kinetic trapping problem is overcome by exchanging noninteracting replicas simulated at different temperatures. Our algorithm uses multiplexed-replicas with a number of independent molecular dynamics runs at each temperature. Exchanges of configurations between these multiplexed-replicas are also tried, rendering the algorithm applicable to large-scale distributed computing (i.e., highly heterogeneous parallel computers with processors having different computational power). We demonstrate the enhanced sampling of this algorithm by simulating the folding thermodynamics of a 23 amino acid miniprotein. We show that better convergence is achieved compared to constant temperature molecular dynamics simulation, with an efficient scaling to large number of computer processors. Indeed, this enhanced sampling results in (to our knowledge) the first example of a replica exchange algorithm that samples a folded structure starting from a completely unfolded state.

This article has been cited by other articles:

				M. B. Kubitzki and B. L. de Groot Molecular Dynamics Simulations Using Temperature-Enhanced Essential Dynamics Replica Exchange Biophys. J., June 15, 2007; 92(12): 4262 - 4270. [Abstract] [Full Text] [PDF]

				J. Juraszek and P. G. Bolhuis Sampling the multiple folding mechanisms of Trp-cage in explicit solvent PNAS, October 24, 2006; 103(43): 15859 - 15864. [Abstract] [Full Text] [PDF]

				Y. Mu, L. Nordenskiold, and J. P. Tam Folding, Misfolding, and Amyloid Protofibril Formation of WW Domain FBP28 Biophys. J., June 1, 2006; 90(11): 3983 - 3992. [Abstract] [Full Text] [PDF]

				V. Cojocaru, R. Klement, and T. M. Jovin Loss of G-A base pairs is insufficient for achieving a large opening of U4 snRNA K-turn motif Nucleic Acids Res., June 13, 2005; 33(10): 3435 - 3446. [Abstract] [Full Text] [PDF]

				Y. Sugita and Y. Okamoto Molecular Mechanism for Stabilizing a Short Helical Peptide Studied by Generalized-Ensemble Simulations with Explicit Solvent Biophys. J., May 1, 2005; 88(5): 3180 - 3190. [Abstract] [Full Text] [PDF]

				E. J. Sorin and V. S. Pande Exploring the Helix-Coil Transition via All-Atom Equilibrium Ensemble Simulations Biophys. J., April 1, 2005; 88(4): 2472 - 2493. [Abstract] [Full Text] [PDF]

				G. S. Jas and K. Kuczera Equilibrium Structure and Folding of a Helix-Forming Peptide: Circular Dichroism Measurements and Replica-Exchange Molecular Dynamics Simulations Biophys. J., December 1, 2004; 87(6): 3786 - 3798. [Abstract] [Full Text] [PDF]

				E. Metwally, H. A. Ismail, J. S. Davison, and R. Mathison A Tree-Based Algorithm for Determining the Effects of Solvation on the Structure of Salivary Gland Tripeptide NH3+-D-PHE-D-GLU-GLY-COO- Biophys. J., September 1, 2003; 85(3): 1503 - 1511. [Abstract] [Full Text] [PDF]

				E. J. Sorin, Y. M. Rhee, B. J. Nakatani, and V. S. Pande Insights into Nucleic Acid Conformational Dynamics from Massively Parallel Stochastic Simulations Biophys. J., August 1, 2003; 85(2): 790 - 803. [Abstract] [Full Text] [PDF]