RESUMO
Most biological membranes are curved, and both lipids and proteins play a role in generating curvature. For any given membrane shape and composition, it is not trivial to determine whether lipids are laterally distributed in a homogeneous or inhomogeneous way, and whether the inter-leaflet distribution is symmetric or not. Here we present a simple computational tool that allows to predict the preference of any lipid type for membranes with positive vs. negative curvature, for any given value of curvature. The tool is based on molecular dynamics simulations of tubular membranes with hydrophilic pores. The pores allow spontaneous, barrierless flip-flop of most lipids, while also preventing differences in pressure between the inner and outer water compartments and minimizing membrane asymmetric stresses. Specifically, we provide scripts to build and analyze the simulations. We test the tool by performing simulations on simple binary lipid mixtures, and we show that, as expected, lipids with negative intrinsic curvature distribute to the tubule inner leaflet, the more so when the radius of the tubular membrane is small. Compared to other existing computational methods, relying on membrane buckles and tethers, our method is based on spontaneous inter-leaflet transport of lipids, and therefore allows to explore lipid distribution in asymmetric membranes. The method can easily be adapted to work with any molecular dynamics code and any force field.