Spontaneous curvature generation by peptides in asymmetric bilayers

Peptides and proteins play crucial roles in membrane remodeling by inducing spontaneous curvature. However, extracting spontaneous curvatures from simulations of asymmetric bilayers is challenging because differential stress (i.e., the difference of the leaflet surface tensions) arising from leaflet area strains can vary substantially among initial conditions. This study investigates peptide-induced spontaneous curvatureδc 0 p in asymmetric bilayers consisting of a single lipid type and a peptide confined to one leaflet;δc 0 p is calculated from the Helfrich equation using the first moment of the lateral pressure tensor and an alternative expression using the differential stress. It is shown that differential stress introduced during initial system generation is effectively relaxed by equilibrating using P21 periodic boundary conditions, which allows lipids to switch leaflets across cell boundaries and equalize their chemical potentials across leaflets. This procedure leads to robust estimates ofδc 0 p for the systems simulated, and is recommended when equality of chemical potentials between the leaflets is a primary consideration.

The Extended Eighth-Shell method for periodic boundary conditions with rotational symmetry

The Eighth-Shell method for parallelization of molecular dynamics simulations has previously been shown to be the most optimal for efficiency at large process counts. However, in its current formulation only the P1 space group is supported for periodic boundary conditions (PBC) and thus reflection and/or rotational crystal symmetries are not supported. In this work, we outline the development and implementation of the Extended Eighth-Shell (EES) method that allows rotational symmetry by using an extended import region compared to the ES method. It simulates only the asymmetric unit and communicates coordinates and forces with images that correspond to P2₁ PBC. The P2₁ PBC has application in lipid bilayer simulations as it can be used to allow lipids to switch leaftlets, thus rapidly balancing the chemical potential difference between the two layers. Our results show that the EES method scales efficiently over large number of processes and can be used for simulations with P2₁ symmetry in an orthorhombic crystal.