Exploring the protein folding free energy landscape: coupling replica exchange method with P3ME/RESPA algorithm.

Area of Science:

• Computational chemistry
• Biophysics
• Molecular dynamics simulations

Background:

• Protein folding is crucial for biological function but challenging to simulate.
• Efficiently sampling the free energy landscape of protein folding requires advanced computational methods.

Purpose of the Study:

• To develop and apply a highly parallel replica exchange method (REM) coupled with a new molecular dynamics algorithm (P3ME/RESPA).
• To efficiently sample the protein folding free energy landscape for beta-hairpin and Trp-cage proteins.

Main Methods:

• Utilized a highly parallel replica exchange method (REM).
• Coupled REM with the particle-particle particle-mesh Ewald (P3ME)/RESPA molecular dynamics algorithm.
• Performed all-atom simulations with the OPLSAA force field and explicit solvent model for up to 64 replicas.
• Simulated beta-hairpin and Trp-cage protein folding over a wide temperature range.

Main Results:

• Successfully sampled the free energy landscape for protein folding.
• Overcame free energy barriers at low temperatures by combining trajectories in temperature and configurational space.
• Revealed detailed folding mechanisms, intermediate structures, and thermodynamic properties for both proteins.
• Determined temperature dependencies for protein folding.

Conclusions:

• The proposed REM-P3ME/RESPA method is efficient for sampling protein folding free energy landscapes.
• Large-scale parallel simulations provide detailed insights into protein folding mechanisms and thermodynamics.
• The method is applicable to various protein systems, including beta-hairpin and Trp-cage.