EnCurv: Simple Technique of Maintaining Global Membrane Curvature in Molecular Dynamics Simulations
Semen Yesylevskyy, Himanshu Khandelia
Abstract
The EnCurv method for maintaining membrane curvature in molecular dynamics simulations is introduced. The method allows maintaining any desired curvature in a sector of lipid membrane bent in a single plane without adding any unphysical interactions into the system and without restrictions on lateral and transversal lipid diffusion and distribution. The current implementation is limited to the membranes curved in a single plane but generalization to arbitrary curvature and membrane topology is possible. The method is simple, easy to implement, and scales linearly with the system size. EnCurv is agnostic to the force field, simulation parameters, and membrane composition. The proof of principle implementation (https://github.com/yesint/EnCurv) is compatible with the majority of modern simulation packages and shows consistent results on the model systems.