Significance: The multifaceted involvement of GM1 as a ligand in many cellular functions has been well recognized. We find that GM1 readily desorbs from the membrane of cell-sized model biomimetic systems (giant unilamellar vesicles). The desorption is unbalanced, creating an asymmetry between the bilayer leaflets. This results in reshaping weakly curved membranes into nanotubular invaginations stabilized by the membrane spontaneous curvature, which we quantify experimentally. Computer simulations confirm the experimental results. Uncovering the role of GM1 as a fine regulator of membrane curvature broadens our perspective on its important function in reshaping neuronal membranes and emphasizes that GM1 desorption can strongly affect the cell membrane morphology.

Currently, there is no comprehensive model for the dynamics of cellular membranes. The understanding of even the basic dynamic processes, such as lateral diffusion of lipids, is still quite limited. Recent studies of one-component membrane systems have shown that instead of single-particle motions, the lateral diffusion is driven by a more complex, concerted mechanism for lipid diffusion (E. Falck et al., J. Am. Chem. Soc., 2008, 130, 44–45), where a lipid and its neighbors move in unison in terms of loosely defined clusters. In this work, we extend the previous study by considering the concerted lipid diffusion phenomena in many-component raft-like membranes. This nature of diffusion phenomena emerge in all the cases we have considered, including both atom-scale simulations of lateral diffusion within rafts and coarse-grained MARTINI simulations of diffusion in membranes characterized by coexistence of raft and non-raft domains. The data allows us to identify characteristic time scales for the concerted lipid motions, which turn out to range from hundreds of nanoseconds to several microseconds. Further, we characterize typical length scales associated with the correlated lipid diffusion patterns and find them to be about 10 nm, or even larger if weak correlations are taken into account. Finally, the concerted nature of lipid motions is also found in dissipative particle dynamics simulations of lipid membranes, clarifying the role of hydrodynamics (local momentum conservation) in membrane diffusion phenomena.

Biomembranes are built up from lipid bilayers with two leaflets that typically differ in their lipid composition. Each lipid molecule stays within one leaflet of the bilayer before it undergoes a transition or flip-flop to the other leaflet. The corresponding flip-flop times are very different for different lipid species and vary over several orders of magnitude. Here, we use molecular dynamics simulations to elucidate the consequences of this separation of time scales for compositionally asymmetric bilayers. We first study bilayers with two lipid components that do not undergo flip-flops on the accessible time scales. In such a situation, one must distinguish a bilayer state in which both leaflets have the same preferred area from another state in which each leaflet is tensionless. However, when we add a third lipid component that undergoes frequent flip-flops, the bilayer relaxes towards the state with tensionless leaflets, not to the state with equal preferred leaflet areas. Furthermore, we show that bilayers with compositional asymmetry acquire a significant spontaneous curvature even if both leaflets are tensionless. Our results can be extended to lipid bilayers with a large number of lipid components provided at least one of these components undergoes frequent flip-flops. For cellular membranes containing lipid pumps, the leaflet tensions also depend on the rates of protein-induced flip-flops.