Peptide folding kinetics from replica exchange molecular dynamics.

Area of Science:

• Computational chemistry
• Biophysics
• Molecular dynamics

Background:

• Determining protein folding and unfolding rates is crucial for understanding protein function.
• Traditional molecular dynamics simulations often struggle to capture rare folding events due to timescale limitations.

Purpose of the Study:

• To develop and validate a method for extracting accurate kinetic information, specifically protein folding and unfolding rates, from replica-exchange molecular dynamics (REMD) simulations.
• To demonstrate the applicability of this method to a helical peptide in explicit water.

Main Methods:

• Utilizing short, continuous trajectory segments between replica exchanges in REMD simulations.
• Estimating short-time propagators in conformation space.
• Constructing a master equation from these propagators.
• Applying the method to a helical peptide in explicit water.

Main Results:

• Accurate rates for local transitions between conformational states and global folding/unfolding events were determined.
• REMD simulations with 5 ps exchange times yielded accurate rates in the 1-100 ns range.
• The obtained rates showed excellent agreement with results from long equilibrium molecular dynamics simulations.

Conclusions:

• Replica-exchange molecular dynamics (REMD) is a viable and accurate method for determining protein folding and unfolding kinetics.
• The developed approach enables efficient extraction of kinetic information from REMD trajectories.
• This method can significantly advance the study of protein dynamics and folding pathways.