The analysis of the structural organization of lipid bilayers is generally performed across the direction normal to the bilayer/water interface, whereas the surface properties of the bilayer at the interface with water are often neglected. Here, we present PackMem, a bioinformatic tool that performs a topographic analysis of the bilayer surface from various molecular dynamics simulations. PackMem unifies and rationalizes previous analyses based on a Cartesian grid. The grid allows identification of surface regions defined as lipid-packing defects where lipids are loosely packed, leading to cavities in which aliphatic carbons are exposed to the solvent, either deep inside or close to the membrane surface. Examples are provided to show that the abundance of lipid-packing defects varies according to the temperature and to the bilayer composition. Because lipid-packing defects control the adsorption of peripheral proteins with hydrophobic insertions, PackMem is instrumental for us to understand and quantify the adhesive properties of biological membranes as well as their response to mechanical perturbations such as membrane deformation.

Phospholipid membranes form cellular barriers but need to be flexible enough to divide by fission. Phospholipids generally contain a saturated fatty acid (FA) at position sn1 whereas the sn2-FA is saturated, monounsaturated or polyunsaturated. Our understanding of the impact of phospholipid unsaturation on membrane flexibility and fission is fragmentary. Here, we provide a comprehensive view of the effects of the FA profile of phospholipids on membrane vesiculation by dynamin and endophilin. Coupled to simulations, this analysis indicates that: (i) phospholipids with two polyunsaturated FAs make membranes prone to vesiculation but highly permeable; (ii) asymmetric sn1-saturated-sn2-polyunsaturated phospholipids provide a tradeoff between efficient membrane vesiculation and low membrane permeability; (iii) When incorporated into phospholipids, docosahexaenoic acid (DHA; omega-3) makes membranes more deformable than arachidonic acid (omega-6). These results suggest an explanation for the abundance of sn1-saturated-sn2-DHA phospholipids in synaptic membranes and for the importance of the omega-6/omega-3 ratio on neuronal functions.