Accelerating All-Atom MD Simulations of Lipids Using a Modified Virtual-Sites Technique

Area of Science:

• Computational chemistry
• Molecular dynamics simulations
• Biophysics

Background:

• All-atom simulations of lipid bilayers are crucial for understanding membrane properties.
• Current methods face limitations in computational efficiency due to small time steps.
• The virtual sites technique offers a way to reduce degrees of freedom and increase simulation speed.

Purpose of the Study:

• To develop and evaluate new implementations of the virtual sites technique for lipid bilayers.
• To assess the impact of virtual sites on the accuracy of simulated membrane properties.
• To enable larger time steps in all-atom simulations using the CHARMM36 force field.

Main Methods:

• Two novel virtual sites implementations were developed: one based on GROMACS derivation and another with a new CH2 group definition.
• Simulations were performed on dipalmitoylphosphatidylcholine (DPPC), 1-palmitoyl-2-oleoylphosphatidylcholine (POPC), and 1,2-dioleoyl-3-phosphatidylcholine (DOPC) lipid bilayers.
• Key physical properties of the membranes were calculated and compared between simulations with and without virtual sites.

Main Results:

• Both virtual sites methods allowed for an increased time step of 5 fs.
• The GROMACS-derived virtual sites on the DOPC bilayer showed excellent agreement with results without virtual sites (Area per lipid: 67.3 ± 0.3 Ų vs. 67.6 ± 0.3 Ų).
• Minor differences in calculated properties were observed, highlighting the accuracy of the virtual sites approach.

Conclusions:

• The virtual sites technique is a powerful tool for accelerating lipid membrane simulations.
• Careful implementation and validation are necessary for reliable results.
• The presented methods are adaptable to other force fields and lipid types.