Multiple branched adaptive steered molecular dynamics.

Area of Science:

• Computational chemistry
• Biophysics
• Molecular dynamics simulations

Background:

• Steered molecular dynamics (SMD) with Jarzynski's equality is crucial for calculating free energy profiles in biological systems like protein folding and ligand binding.
• Standard SMD methods often suffer from rare events dominating averages, necessitating extensive computations.
• Adaptive steered molecular dynamics (ASMD) was developed to improve efficiency by selecting significant trajectories.

Purpose of the Study:

• To introduce and evaluate a novel, highly efficient multiple branching adaptive steered molecular dynamics (MB-ASMD) method.
• To enhance the sampling of important trajectories and eliminate non-contributing segments in molecular dynamics simulations.
• To improve the accuracy and reduce the computational cost of free energy profile calculations.

Main Methods:

• Implementation of a multiple branching strategy within the adaptive steered molecular dynamics framework.
• Selection of multiple configurations along the reaction coordinate based on work trajectory distributions.
• Application of MB-ASMD to analyze hydrogen bond breaking in decaalanine peptide during helix-to-coil transition.

Main Results:

• MB-ASMD achieves more complete enhanced sampling of critical trajectories while discarding irrelevant ones.
• The method demonstrates faster convergence of potential of mean force (PMF) compared to ASMD, even at high pulling speeds.
• MB-ASMD significantly reduces the likelihood of simulations getting trapped in non-significant pathways and yields less noisy results.

Conclusions:

• MB-ASMD offers superior efficiency and accuracy for free energy calculations in complex biological systems.
• The multiple branching approach enhances sampling completeness and computational savings over existing ASMD methods.
• This novel technique provides a powerful tool for investigating molecular mechanisms in biophysics and computational chemistry.