Efficient molecular dynamics using geodesic integration and solvent-solute splitting.

Area of Science:

• Computational chemistry
• Molecular dynamics simulations
• Biophysics

Background:

• Langevin dynamics is crucial for simulating molecular motion.
• Holonomic constraints are often necessary for biomolecular simulations.
• Existing methods face accuracy and efficiency limitations.

Purpose of the Study:

• To develop a novel, accurate, and efficient integrator for Langevin dynamics with holonomic constraints.
• To enhance the sampling efficiency of molecular dynamics simulations for solvated biomolecules.

Main Methods:

• Decomposition of Langevin dynamics into geodesic flow, constrained impulse, and constrained diffusion.
• Strategic ordering of system components for integrator development.
• Integration of geodesic method with solvent-solute force splitting.

Main Results:

• A new integrator demonstrates an order of magnitude improvement in configurational average accuracy.
• Large simulation steps (≥ 8 fs) are feasible for solvated biomolecules.
• Molecular dynamics sampling efficiency is approximately doubled without altering diffusion rates.

Conclusions:

• The presented method offers a significant advancement in constrained molecular dynamics.
• The approach is readily implementable in standard simulation software.
• This facilitates more efficient and accurate biomolecular simulations.