Gaussian Accelerated Molecular Dynamics in OpenMM.

Area of Science:

• Computational chemistry and biophysics
• Molecular dynamics simulations
• Biomolecular modeling

Background:

• Gaussian accelerated molecular dynamics (GaMD) is a powerful computational method for enhanced sampling and free energy calculations.
• Efficient simulation techniques are crucial for understanding complex biomolecular processes like protein and RNA folding.
• Implementing advanced simulation methods in widely used packages like OpenMM broadens their accessibility.

Purpose of the Study:

• To implement Gaussian accelerated molecular dynamics (GaMD) within the OpenMM simulation package.
• To validate the accuracy and efficiency of the OpenMM GaMD implementation using model systems.
• To demonstrate the utility of GaMD in OpenMM for studying RNA folding and ligand-RNA interactions.

Main Methods:

• Implementation of GaMD algorithms in the OpenMM molecular dynamics simulation framework.
• Validation using alanine dipeptide simulations, comparing GaMD results with conventional molecular dynamics (cMD).
• Application to complex systems: RNA tetraloop folding and ligand binding to HIV-1 Tar RNA.

Main Results:

• GaMD simulations achieved comparable free energy profiles to significantly longer cMD simulations for alanine dipeptide.
• The OpenMM GaMD implementation successfully captured folding pathways for hyperstable RNA tetraloops (UUCG, GCAA, CUUG).
• GaMD simulations accurately modeled the binding of the rbt203 ligand to HIV-1 Tar RNA, highlighting key electrostatic interactions.

Conclusions:

• The OpenMM implementation of GaMD provides an efficient and accurate tool for enhanced sampling and free energy calculations.
• GaMD in OpenMM is effective for studying critical electrostatic interactions in biomolecular systems, including RNA folding and ligand binding.
• This integration expands the applicability of GaMD for simulating proteins, RNA, and other biomolecules.