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

Citations Over Time

Abstract

We present two new implementations of the virtual sites technique which completely suppresses the degrees of freedom of the hydrogen atoms in a lipid bilayer allowing for an increased time step of 5 fs in all-atom simulations of the CHARMM36 force field. One of our approaches uses the derivation of the virtual sites used in GROMACS while the other uses a new definition of the virtual sites of the CH2 groups. Our methods is tested on a DPPC (no unsaturated chain), a POPC (one unsaturated chain), and a DOPC (two unsaturated chains) lipid bilayers. We calculate various physical properties of the membrane of our simulations with and without virtual sites and explain the differences and similarity observed. The best agreements are obtained for the GROMACS original virtual sites on the DOPC bilayer where we get an area per lipid of 67.3 ± 0.3 Å(2) without virtual sites and 67.6 ± 0.3 Å(2) with virtual sites. In conclusion the virtual-sites technique on lipid membranes is a powerful simulation tool, but it should be used with care. The procedure can be applied to other force fields and lipids in a straightforward manner.

Related Papers

• → Effects of a Carane Derivative Local Anesthetic on a Phospholipid Bilayer Studied by Molecular Dynamics Simulation(2003)31 cited
• → Effects of imidazolium-based ionic surfactants on the size and dynamics of phosphatidylcholine bilayers with saturated and unsaturated chains(2015)7 cited
• → Concentration effect of cimetidine with POPC bilayer: a molecular dynamics simulation study(2016)7 cited
• → Aggregation of Lysozyme in the Presence of a Mixed Bilayer of POPC and POPG(2021)3 cited
• → Molecular Dynamics Simulations of Structural Changes of Lipid Bilayers Under High Pressure(2021)