Self-consistent calculation of protein folding pathways.

Area of Science:

• Computational biology
• Biophysics
• Molecular dynamics

Background:

• Protein folding is crucial for biological function.
• Simulating conformational transitions is computationally intensive.
• Accurate simulation methods are needed to understand protein dynamics.

Purpose of the Study:

• To introduce an efficient iterative algorithm for simulating protein folding.
• To directly yield reaction pathways and estimate reaction coordinates.
• To validate the algorithm's accuracy and computational efficiency.

Main Methods:

• Developed an iterative algorithm based on the Langevin equation.
• Derived a self-consistent stochastic equation of motion.
• Validated against all-atom molecular dynamics (MD) simulations on a supercomputer.
• Applied the algorithm to a 130-amino acid protein.

Main Results:

• The algorithm efficiently simulates protein folding and conformational transitions.
• It directly provides reaction pathways and a stochastic estimate of the reaction coordinate.
• Validation against MD simulations confirmed the approach's accuracy.
• Demonstrated computational efficiency on a larger protein system.

Conclusions:

• The iterative algorithm offers an efficient and accurate method for simulating protein dynamics.
• This approach can accelerate the study of protein folding pathways.
• The method holds promise for exploring complex conformational changes in proteins.