Accelerating Simulations

We have implemented the virtual sites technique for lipid membranes, 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. The hydrogen virtual sites still have non-bonded interactions, and thus the method is more “accurate” than the united force field method, where information on hydrogen atoms is lost completely.

The method performs well for DPPC, POPC and DOPC membranes. We are currently working on application of the method to other force fields and lipids.

The positions of the H atoms on the CH2 group of a lipid tail can be calculated from the positions of the three nearest heavy atoms using simple geometric rules.

Relevant Publications:

Loubet, B., Kopec, W., and Khandelia, H. (2014) Accelerating All-Atom MD Simulations of Lipids Using a Modified Virtual-Sites Technique, J Chem Theory Comput 10, 5690-5695.