Biological lipid membranes are generally asymmetric, not only with respect to the composition of the two membrane leaflets but also with respect to the state of mechanical stress on the two sides. Computer simulations of such asymmetric membranes pose unique challenges with respect to the choice of boundary conditions and ensemble in which such simulations are to be carried out. Here, we demonstrate an alternative to the usual choice of fully periodic boundary conditions: The membrane is only periodic in one direction, with free edges running parallel to the single direction of periodicity. In order to maintain bilayer asymmetry under these conditions, nanoscale "sticky tapes" are adhered to the membrane edges in order to prevent lipid flip-flop across the otherwise open edge. In such semi-periodic simulations, the bilayer is free to choose both its area and mean curvature, allowing for minimization of the bilayer elastic free energy. We implement these principles in a highly coarse-grained model and show how even the simplest examples of such simulations can reveal useful membrane elastic properties, such as the location of the monolayer neutral surface.
Read full abstract