Implementation of Telescoping Boxes in Adaptive Steered Molecular Dynamics.

Area of Science:

• Computational chemistry
• Molecular dynamics simulations
• Biophysics

Background:

• Steered molecular dynamics (SMD) accelerates simulations of long-time processes like protein unfolding.
• Extracting generalizable information from nonequilibrium SMD simulations remains challenging.
• Jarzynski's equality enables free energy calculation but requires slow-converging nonequilibrium trajectory sampling.

Purpose of the Study:

• To reduce the computational cost associated with large solvent boxes in adaptive steered molecular dynamics (ASMD).
• To optimize solvent molecule usage across different simulation stages.
• To maintain accuracy while improving efficiency in molecular dynamics simulations.

Main Methods:

• Introduction of telescoping box schemes within adaptive steered molecular dynamics (ASMD).
• Adjustment of solvent box sizes between simulation stages to vary solvent molecule numbers.
• Benchmarking on an α-helical peptide (Ala30) to assess potential of mean force and hydrogen bonds.

Main Results:

• Telescoping box schemes in ASMD significantly reduce computational cost.
• The method introduces minimal numerical error compared to standard simulations.
• Optimized solvent box sizes lead to substantial savings in computational resources.

Conclusions:

• Telescoping box schemes offer an efficient approach for ASMD simulations.
• This technique effectively addresses the high computational cost of large solvent boxes.
• The method provides a viable strategy for accurate and cost-effective molecular dynamics studies.