Interleaflet Interaction and Asymmetry in Phase Separated Lipid Bilayers: Molecular Dynamics Simulations

Efficient calculation of orientation-dependent lipid dynamics from membrane simulations

Molecular dynamics simulations of lipid membranes have become increasingly impactful in biophysics because they offer atomistic resolution of structural fluctuations in relation to their functional outputs. Yet quantitative characterization of multiscale processes is a formidable challenge due to the distribution of motions that evade analysis of discrete simulation data. Here we investigate the efficient calculation of CH bond relaxation rates from membrane simulations. Widely used computational approaches offer numerical simplicity but fall short of capturing crucial aspects of the orientation dependence of the dynamics. To circumvent this problem, we introduced a robust framework based on liquid crystal theory which considers explicitly the CH bond motions with respect to the director axis (bilayer normal). Analysis of the orientation dependence of the dynamics shows excellent agreement with experiment, illustrating how the ordering potential affects the calculated relaxation rates. Furthermore, a fit-based resampling of the autocorrelation function of the bond fluctuations validates the new approach for low-temporal resolution data. The recovered relaxation rates indicate that at short timescales, both with and without cholesterol, the local motions of CH bonds describe the bilayer microviscosity and resemble liquid hydrocarbons. Our results establish the critical role of the orientational anisotropy in analysis of membrane simulations, explain fundamental aspects of lipid dynamics, and provide guidelines for extracting information that can be compared to experimental data.

The rise of FTIR spectroscopy in the characterization of asymmetric lipid membranes

In contrast to symmetric unilamellar liposomes (sLUVs) prepared from a mixture of different lipids, asymmetric ones (aLUVs) with different lipid composition in the inner and outer membrane leaflets are more suitable model systems of eukaryotic plasma membranes. However, apart from the challenging preparation of asymmetric liposomes and small amounts of obtained asymmetric unilamellar liposomes (aLUVs), a major drawback is the qualitative characterization of asymmetry, as each of the techniques used so far has certain limitations. In this regard, we prepared aLUVs composed dominantly of DPPC(out)/DPPS(in) lipids and, along with 1H NMR and DSC characterization, we showed for the first time how FTIR spectroscopy can be used in the presence of (a)symmetry between DPPC/DPPS lipid bilayers. Using second derivative FTIR spectra we demonstrated not only that the hydration of lipids glycerol backbone and choline moiety of DPPC differs in s/aLUVs, but in addition that the lateral interactions between hydrocarbon chains during the phase change display different trend in s/aLUVs. Molecular dynamics simulations confirmed different chain ordering and packing between s/a bilayers, with a significant influence of temperature, i.e. membrane phase.

The Membrane-Mediated Interaction of Liquid-Ordered Lipid Domains in the Presence of Amphipathic Peptides

The lipid membranes of living cells are composed of a large number of lipid types and can undergo phase separation with the formation of nanometer-scale liquid-ordered lipid domains, also called rafts. Raft coalescence, i.e., the fusion of lipid domains, is involved in important cell processes, such as signaling and trafficking. In this work, within the framework of the theory of elasticity of lipid membranes, we explore how amphipathic peptides adsorbed on lipid membranes may affect the domain-domain fusion processes. We show that the elastic deformations of lipid membranes drive amphipathic peptides to the boundary of lipid domains, which leads to an increase in the average energy barrier of the domain-domain fusion, even if the surface concentration of amphipathic peptides is low and the domain boundaries are only partially occupied by the peptides. This inhibition of the fusion of lipid domains may lead to negative side effects of using amphipathic peptides as antimicrobial agents.

Membrane Domain Anti-Registration Induces an Intrinsic Transmembrane Potential

Plasma membrane segregation into various nanoscale membrane domains is driven by distinct interactions between diverse lipids and proteins. Among them, liquid-ordered (Lo) membrane domains are defined as "lipid rafts" and liquid-disordered (Ld) ones as "lipid non-rafts". Using model membrane systems, both intra-leaflet and inter-leaflet dynamics of these membrane domains are widely studied. Nevertheless, the biological impact of the latter, which is accompanied by membrane domain registration/anti-registration, is far from clear. Hence, in this work, we studied the biological relevance of the membrane domain anti-registration using both all-atom molecular dynamics (MD) simulations and confocal fluorescence microscopy. All-atom MD simulations suggested an intrinsic transmembrane potential for the case of the membrane anti-registration (Lo/Ld). Meanwhile, confocal fluorescence microscopy experiments of HeLa and 293T cell lines indicated that membrane cholesterol depletion could significantly alter the transmembrane potential of cells. Considering differences in the cholesterol content between Lo and Ld membrane domains, our confocal fluorescence microscopy experiments are consistent with our all-atom MD simulations. In short, membrane domain anti-registration induces local membrane asymmetry and, thus, an intrinsic transmembrane potential.

Building Asymmetric Lipid Bilayers for Molecular Dynamics Simulations: What Methods Exist and How to Choose One?

The compositional asymmetry of biological membranes has attracted significant attention over the last decade. Harboring more differences from symmetric membranes than previously appreciated, asymmetric bilayers have proven quite challenging to study with familiar concepts and techniques, leaving many unanswered questions about the reach of the asymmetry effects. One particular area of active research is the computational investigation of composition- and number-asymmetric lipid bilayers with molecular dynamics (MD) simulations. Offering a high level of detail into the organization and properties of the simulated systems, MD has emerged as an indispensable tool in the study of membrane asymmetry. However, the realization that results depend heavily on the protocol used for constructing the asymmetric bilayer models has sparked an ongoing debate about how to choose the most appropriate approach. Here we discuss the underlying source of the discrepant results and review the existing methods for creating asymmetric bilayers for MD simulations. Considering the available data, we argue that each method is well suited for specific applications and hence there is no single best approach. Instead, the choice of a construction protocol-and consequently, its perceived accuracy-must be based primarily on the scientific question that the simulations are designed to address.

Phosphatidylinositol 4,5-Bisphosphate Sensing Lipid Raft via Inter-Leaflet Coupling Regulated by Acyl Chain Length of Sphingomyelin

Phosphatidylinositol 4,5-bisphosphate (PIP₂) is an important molecule located at the inner leaflet of cell membrane, where it serves as anchoring sites for a cohort of membrane-associated molecules and as a broad-reaching signaling intermediate. The lipid raft is thought as the major platform recruiting proteins for signal transduction and also known to mediate PIP₂ accumulation across the membrane. While the significance of this cross-membrane coupling is increasingly appreciated, it remains unclear whether and how PIP₂ senses the dynamic change of the ordered lipid domains over the packed hydrophobic core of the bilayer. Herein, by means of molecular dynamic simulation, we reveal that inner PIP₂ molecules can sense the outer lipid domain via inter-leaflet coupling, and the coupling manner is dictated by the acyl chain length of sphingomyelin (SM) partitioned to the lipid domain. Shorter SM promotes membrane domain registration, whereby PIP₂ accumulates beneath the domain across the membrane. In contrast, the anti-registration is thermodynamically preferred if the lipid domain has longer SM due to the hydrophobic mismatch between the corresponding acyl chains in SM and PIP₂. In this case, PIP₂ is expelled by the domain with a higher diffusivity. These results provide molecular insights into the regulatory mechanism of correlation between the outer lipid domain and inner PIP₂, both of which are critical components for cell signal transduction.

Characterizing Oscillatory and Excitability Regimes in a Protein-Free Lipid Membrane

Observations of electric potential oscillations in artificial lipid bilayers near the order-disorder transition indicate the existence of a stable limit cycle and, therefore, the possibility of producing excitable signals close to the bifurcation. We present a theoretical investigation of membrane oscillatory and excitability regimes induced by an increase in ion permeability at the order-disorder transition. The model accounts for the coupled effects of state-dependent permeability, membrane charge density, and hydrogen ion adsorption. A bifurcation diagram shows a transition between fixed-point and limit cycle solutions, enabling both oscillatory and excitability responses at different values of the acid association parameter. Oscillations are identified in terms of the membrane state, electric potential difference, and ion concentration near the membrane. The emerging voltage and time scales agree with measurements. Excitability is demonstrated by applying an external electric current stimulus, and the emerging signals display a threshold response and the appearance of repetitive signals upon using a long-lasting stimulus. The approach highlights the important role of the order-disorder transition, enabling membrane excitability in the absence of specialized proteins.

Lysolipids regulate raft size distribution

The lipid matrix of cellular membranes, directly and indirectly, regulates many vital functions of the cell. The diversity of lipids in membranes leads to the formation of ordered domains called rafts, which play a crucial role in signal transduction, protein sorting and other cellular processes. Rafts are believed to impact the development of different neurodegenerative diseases, such as Alzheimer’s, Parkinson’s, Huntington’s ones, amyotrophic lateral sclerosis, some types of cancer, etc. These diseases correlate with the change in the membrane lipid composition resulting from an oxidative stress, age-related processes, dysfunction of proteins, and many others. In particular, a lot of studies report a significant rise in the level of lysolipids. Physicochemical properties of rafts are determined by membrane composition, in particular, by the content of lysolipids. Lysolipids may thus regulate raft-involving processes. However, the exact mechanism of such regulation is unknown. Although studying rafts in vivo still seems to be rather complicated, liquid-ordered domains are well observed in model systems. In the present study, we used atomic force microscopy (AFM) to examine how lysophospholipids influence the liquid-ordered domains in model ternary membranes. We demonstrated that even a small amount of lysolipids in a membrane significantly impacts domain size depending on the saturation of the lysolipid hydrocarbon tails and the amount of cholesterol. The mixture with the bigger relative fraction of cholesterol was more susceptible to the action of lysolipids. This data helped us to generalize our previous theoretical model of the domain size regulation by lipids with particular molecular shape expanding it to the case of lysolipids and dioleoylglycerol.

Small ionic radii limit time step in Martini 3 molecular dynamics simulations

Among other improvements, the Martini~3 coarse-grained force field provides a more accurate description of the solvation of protein pockets and channels through the consistent use of various bead types and sizes. Here, we show that the representation of Na+ and Cl- ions as ``tiny' (TQ5) beads limits the accessible time step to 25~fs. By contrast, with Martini~2, time steps of 30-40~fs were possible for lipid bilayer systems without proteins. This limitation is relevant for systems that require long equilibration times. We derive a quantitative kinetic model of time-integration instabilities in molecular dynamics as a function of time step, ion concentration and mass, system size, and simulation time. We demonstrate that ion-water interactions are the main source of instability at physiological conditions, followed closely by ion-ion interactions. We show that increasing the ionic masses makes it possible to use time steps up to 40~fs with minimal impact on static equilibrium properties and on dynamical quantities such as lipid and ion diffusion coefficients. Increasing the size of the bead representing the ions (and thus changing their hydration) also permits longer time steps. For a soluble protein, we find that increasing the mass of tiny beads also on the protein permits simulations with 30-fs time step. The use of larger time steps in Martini~3 results in a more efficient exploration of configuration space. The kinetic model of MD simulation crashes can be used to determine the maximum allowed time step upfront for an efficient use of resources and whenever sampling efficiency is critical.

Order vs. Disorder: Cholesterol and Omega-3 Phospholipids Determine Biomembrane Organization

Lipid structural diversity strongly affects biomembrane chemico-physical and structural properties in addition to membrane-associated events. At high concentrations, cholesterol increases membrane order and rigidity, while polyunsaturated lipids are reported to increase disorder and flexibility. How these different tendencies balance in composite bilayers is still controversial. In this study, electron paramagnetic resonance spectroscopy, small angle neutron scattering, and neutron reflectivity were used to investigate the structural properties of cholesterol-containing lipid bilayers in the fluid state with increasing amounts of polyunsaturated omega-3 lipids. Either the hybrid 1-stearoyl-2-docosahexaenoyl-sn-glycero-3-phosphocholine or the symmetric 1,2-docosahexaenoyl-sn-glycero-3-phosphocholine were added to the mixture of the naturally abundant 1-palmitoyl-2-oleyl-sn-glycero-3-phosphocholine and cholesterol. Our results indicate that the hybrid and the symmetric omega-3 phospholipids affect the microscopic organization of lipid bilayers differently. Cholesterol does not segregate from polyunsaturated phospholipids and, through interactions with them, is able to suppress the formation of non-lamellar structures induced by the symmetric polyunsaturated lipid. However, this order/disorder balance leads to a bilayer whose structural organization cannot be ascribed to either a liquid ordered or to a canonical liquid disordered phase, in that it displays a very loose packing of the intermediate segments of lipid chains.

Enthalpic and entropic contributions to interleaflet coupling drive domain registration and antiregistration in biological membrane

Biological membrane is a complex self-assembly of lipids, sterols, and proteins organized as a fluid bilayer of two closely stacked lipid leaflets. Differential molecular interactions among its diverse constituents give rise to heterogeneities in the membrane lateral organization. Under certain conditions, heterogeneities in the two leaflets can be spatially synchronized and exist as registered domains across the bilayer. Several contrasting theories behind mechanisms that induce registration of nanoscale domains have been suggested. Following a recent study showing the effect of position of lipid tail unsaturation on domain registration behavior, we decided to develop an analytical theory to elucidate the driving forces that create and maintain domain registry across leaflets. Towards this, we formulated a Hamiltonian for a stacked lattice system where site variables capture the lipid molecular properties such as the position of unsaturation and various other interactions that could drive phase separation and interleaflet coupling. We solve the Hamiltonian using Monte Carlo simulations and create a complete phase diagram that reports the presence or absence of registered domains as a function of various Hamiltonian parameters. We find that the interleaflet coupling should be described as a competing enthalpic contribution due to interaction of lipid tail termini, primarily due to saturated-saturated interactions, and an interleaflet entropic contribution from overlap of unsaturated tail termini. A higher position of unsaturation is seen to provide weaker interleaflet coupling. Thermodynamically stable nanodomains could also be observed for certain points in the parameter space in our bilayer model, which were further verified by carrying out extended Monte Carlo simulations. These persistent noncoalescing registered nanodomains close to the lower end of the accepted nanodomain size range also point towards a possible "nanoscale" emulsion description of lateral heterogeneities in biological membrane leaflets.

Hydrophobic Mismatch Controls the Mode of Membrane-Mediated Interactions of Transmembrane Peptides

Various cellular processes require the concerted cooperative action of proteins. The possibility for such synchronization implies the occurrence of specific long-range interactions between the involved protein participants. Bilayer lipid membranes can mediate protein-protein interactions via relatively long-range elastic deformations induced by the incorporated proteins. We considered the interactions between transmembrane peptides mediated by elastic deformations using the framework of the theory of elasticity of lipid membranes. An effective peptide shape was assumed to be cylindrical, hourglass-like, or barrel-like. The interaction potentials were obtained for membranes of different thicknesses and elastic rigidities. Cylindrically shaped peptides manifest almost neutral average interactions-they attract each other at short distances and repel at large ones, independently of membrane thickness or rigidity. The hourglass-like peptides repel each other in thin bilayers and strongly attract each other in thicker bilayers. On the contrary, the barrel-like peptides repel each other in thick bilayers and attract each other in thinner membranes. These results potentially provide possible mechanisms of control for the mode of protein-protein interactions in membrane domains with different bilayer thicknesses.

Characteristic lengths of transmembrane peptides controlling their tilt and lateral distribution between membrane domains

Lipids and proteins of plasma membranes of eukaryotic cells are supposed to form protein-lipid domains, characterized by a different molecular order, bilayer thickness, and elastic parameters. Several mechanisms of preferable distribution of transmembrane proteins to the ordered or disordered membrane domains have been revealed. The mismatch between the length of the protein transmembrane domain and hydrophobic thickness of the lipid bilayer is considered to be an important driving force of protein lateral sorting. Utilizing the continuum theory of elasticity, we analyzed optimal configurations and preferable membrane domains for single-pass transmembrane peptides of various hydrophobic lengths and effective molecular shapes. We obtained that short transmembrane peptides stand perpendicularly to the membrane plane. The exceedance of a certain characteristic length leads to the tilt of the peptide. This length depends on the bilayer thickness. Thus, in the membrane with coexisting ordered (thicker) and disordered (thinner) phases tilting of the peptide in each phase is governed by its individual characteristic length. The lateral distribution of the peptides between ordered and disordered membrane domains is shown to be described by two additional characteristic lengths. The exceedance of the smaller one drives the peptide towards a more ordered and thicker membrane, while the exceedance of the larger characteristic length switches the preferable membrane domain from ordered and thicker to disordered and thinner. Thus, membrane proteins with long enough transmembrane domains are predicted to accumulate in the thinner disordered membrane as compared to the thicker ordered bilayer. For hourglass-like and barrel-like shaped transmembrane peptides the specific regime of sorting was obtained: the peptides distributed almost equally between the phases in a wide range of peptide lengths. This finding allowed explaining the experimental data on lateral distribution of transmembrane peptide tLAT.

Dimensional Changes in Lipid Rafts from Human Brain Cortex Associated to Development of Alzheimer's Disease. Predictions from an Agent-Based Mathematical Model

Alzheimer’s disease (AD) is a neurodegenerative disease caused by abnormal functioning of critical physiological processes in nerve cells and aberrant accumulation of protein aggregates in the brain. The initial cause remains elusive—the only unquestionable risk factor for the most frequent variant of the disease is age. Lipid rafts are microdomains present in nerve cell membranes and they are known to play a significant role in the generation of hallmark proteinopathies associated to AD, namely senile plaques, formed by aggregates of amyloid β peptides. Recent studies have demonstrated that human brain cortex lipid rafts are altered during early neuropathological phases of AD as defined by Braak and Braak staging. The lipid composition and physical properties of these domains appear altered even before clinical symptoms are detected. Here, we use a coarse grain molecular dynamics mathematical model to predict the dimensional evolution of these domains using the experimental data reported by our group in human frontal cortex. The model predicts significant size and frequency changes which are detectable at the earliest neuropathological stage (ADI/II) of Alzheimer’s disease. Simulations reveal a lower number and a larger size in lipid rafts from ADV/VI, the most advanced stage of AD. Paralleling these changes, the predictions also indicate that non-rafts domains undergo simultaneous alterations in membrane peroxidability, which support a link between oxidative stress and AD progression. These synergistic changes in lipid rafts dimensions and non-rafts peroxidability are likely to become part of a positive feedback loop linked to an irreversible amyloid burden and neuronal death during the evolution of AD neuropathology.

Sphingomyelinase-Mediated Multitimescale Clustering of Ganglioside GM1 in Heterogeneous Lipid Membranes

Several signaling processes in the plasma membrane are intensified by ceramides that are formed by sphingomyelinase-mediated hydrolysis of sphingomyelin. These ceramides trigger clustering of signaling-related biomolecules, but how they concentrate such biomolecules remains unclear. Here, the spatiotemporal localization of ganglioside GM1, a glycolipid receptor involved in signaling, during sphingomyelinase-mediated hydrolysis is described. Real-time visualization of the dynamic remodeling of the heterogeneous lipid membrane that occurs due to sphingomyelinase action is used to examine GM1 clustering, and unexpectedly, it is found that it is more complex than previously thought. Specifically, lipid membranes generate two distinct types of condensed GM1: 1) rapidly formed but short-lived GM1 clusters that are formed in ceramide-rich domains nucleated from the liquid-disordered phase; and 2) late-onset yet long-lasting, high-density GM1 clusters that are formed in the liquid-ordered phase. These findings suggest that multiple pathways exist in a plasma membrane to synergistically facilitate the rapid amplification and persistence of signals.

LiPyphilic: A Python Toolkit for the Analysis of Lipid Membrane Simulations

Molecular dynamics simulations are now widely used to study emergent phenomena in lipid membranes with complex compositions. Here, we present LiPyphilic-a fast, fully tested, and easy-to-install Python package for analyzing such simulations. Analysis tools in LiPyphilic include the identification of cholesterol flip-flop events, the classification of local lipid environments, and the degree of interleaflet registration. LiPyphilic is both force field- and resolution-agnostic, and by using the powerful atom selection language of MDAnalysis, it can handle membranes with highly complex compositions. LiPyphilic also offers two on-the-fly trajectory transformations to (i) fix membranes split across periodic boundaries and (ii) perform nojump coordinate unwrapping. Our implementation of nojump unwrapping accounts for fluctuations in the box volume under the NPT ensemble-an issue that most current implementations have overlooked. The full documentation of LiPyphilic, including installation instructions and links to interactive online tutorials, is available at https://lipyphilic.readthedocs.io/en/latest.

Model Membrane Systems Used to Study Plasma Membrane Lipid Asymmetry

It is well known that the lipid distribution in the bilayer leaflets of mammalian plasma membranes (PMs) is not symmetric. Despite this, model membrane studies have largely relied on chemically symmetric model membranes for the study of lipid–lipid and lipid–protein interactions. This is primarily due to the difficulty in preparing stable, asymmetric model membranes that are amenable to biophysical studies. However, in the last 20 years, efforts have been made in producing more biologically faithful model membranes. Here, we review several recently developed experimental and computational techniques for the robust generation of asymmetric model membranes and highlight a new and particularly promising technique to study membrane asymmetry.

Investigating lipid headgroup composition within epithelial membranes: a systematic review

Membrane lipid composition is often quoted within the literature, but with very little insight into how or why these compositions vary when compared to other biological membranes. One prominent area that lacks understanding in terms of rationale for lipid variability is the human gastro-intestinal tract (GIT). We have carried out a comprehensive systematic literature search to ascertain the key lipid components of epithelial membranes, with a particular focus on addressing the human GIT and to use compositional data to understand structural aspects of biological membranes. Both bacterial outer membranes and the human erythrocyte membrane were used as a comparison for the mammalian [epithelial] membranes and to understand variations in lipid presence. We show that phosphatidylcholine (PC) lipid types tend to dominate (33%) with phosphatidylethanolamines (PE) and cholesterol having very similar abundances (25 and 23% respectively). This systematic review presents a detailed insight into lipid headgroup composition and roles in various membrane types, with a summary of the distinction between the major lipid bilayer forming lipids and how peripheral lipids regulate charge and fluidity. The variety of lipids present in biological membranes is discussed and rationalised in terms function as well as cellular position.

Transmembrane potential of physiologically relevant model membranes: Effects of membrane asymmetry

Transmembrane potential difference (V_m) plays important roles in regulating various biological processes. At the macro level, V_m can be experimentally measured or calculated using the Nernst or Goldman-Hodgkin-Katz equation. However, the atomic details responsible for its generation and impact on protein and lipid dynamics still need to be further elucidated. In this work, we performed a series of all-atom molecular dynamics (MD) simulations of symmetric model membranes of various lipid compositions and cation contents to evaluate the relationship between membrane asymmetry and V_m. Specifically, we studied the impact of the asymmetric distribution of POPS (1-palmitoyl-2-oleoyl-sn-glycero-3-phospho-l-serine), PIP2 (phosphatidylinositol 4,5-bisphosphate), as well as Na⁺ and K⁺ on V_m using atomically detailed MD simulations of symmetric model membranes. The results suggest that, for an asymmetric POPC-POPC/POPS bilayer in the presence of NaCl, the presence of the monovalent anionic lipid POPS in the inner leaflet polarizes the membrane (ΔV_m < 0). Intriguingly, replacing a third of the POPS lipids by the polyvalent anionic signaling lipid PIP2 counteracts this effect, resulting in a smaller negative membrane potential. We also found that replacing Na⁺ ions in the inner region by K⁺ depolarizes the membrane (ΔV_m > 0). These divergent effects arise from variations in the strength of cation-lipid interactions and are correlated with changes in lipid chain order and head-group orientation.

Role of Lipids in Morphogenesis of T-Cell Microvilli

T cells communicate with the environment via surface receptors. Cooperation of surface receptors regulates T-cell responses to diverse stimuli. Recently, finger-like membrane protrusions, microvilli, have been demonstrated to play a role in the organization of receptors and, hence, T-cell activation. However, little is known about the morphogenesis of dynamic microvilli, especially in the cells of immune system. In this review, I focus on the potential role of lipids and lipid domains in morphogenesis of microvilli. Discussed is the option that clustering of sphingolipids with phosphoinositides at the plasma membrane results in dimpling (curved) domains. Such domains can attract phosphoinositide-binding proteins and stimulate actin cytoskeleton reorganization. This process triggers cortical actin opening and bundling of actin fibres to support the growing of microvilli. Critical regulators of microvilli morphogenesis in T cells are unknown. At the end, I suggest several candidates with a potential to organize proteins and lipids in these structures.

Investigation of the domain line tension in asymmetric vesicles prepared via hemifusion

The plasma membrane (PM) is asymmetric in lipid composition. The distinct and characteristic lipid compositions of the exoplasmic and cytoplasmic leaflets lead to different lipid-lipid interactions and physical-chemical properties in each leaflet. The exoplasmic leaflet possesses an intrinsic ability to form coexisting ordered and disordered fluid domains, whereas the cytoplasmic leaflet seems to form a single fluid phase. To better understand the interleaflet interactions that influence domains, we compared asymmetric model membranes that capture salient properties of the PM with simpler symmetric membranes. Using asymmetric giant unilamellar vesicles (aGUVs) prepared by hemifusion with a supported lipid bilayer, we investigate the domain line tension that characterizes the behavior of coexisting ordered + disordered domains. The line tension can be related to the contact perimeter of the different phases. Compared to macroscopic phase separation, the appearance of modulated phases was found to be a robust indicator of a decrease in domain line tension. Symmetric GUVs of 1,2-distearoyl-sn-glycero-3-phosphocholine (DSPC)/1,2-dioleoyl-sn-glycero-3-phosphocholine (DOPC)/1-palmitoyl-2-oleoyl-glycero-3-phosphocholine (POPC)/cholesterol (chol) were formed into aGUVs by replacing the GUV outer leaflet with DOPC/chol = 0.8/0.2 in order to create a cytoplasmic leaflet model. These aGUVs revealed lower line tension for the ordered + disordered domains of the exoplasmic model leaflet.

Domain formation in bicomponent vesicles induced by composition-curvature coupling

Lipid vesicles composed of a mixture of two types of lipids are studied by intensive Monte Carlo numerical simulations. The coupling between the local composition and the membrane shape is induced by two different spontaneous curvatures of the components. We explore the various morphologies of these biphasic vesicles coupled to the observed patterns such as nano-domains or labyrinthine mesophases. The effect of the difference in curvatures, the surface tension, and the interaction parameter between components is thoroughly explored. Our numerical results quantitatively agree with the previous analytical results obtained by Gueguen et al. [Eur. Phys. J. E 37, 76 (2014)] in the disordered (high temperature) phase. Numerical simulations allow us to explore the full parameter space, especially close to and below the critical temperature, where analytical results are not accessible. Phase diagrams are constructed and domain morphologies are quantitatively studied by computing the structure factor and the domain size distribution. This mechanism likely explains the existence of nano-domains in cell membranes as observed by super-resolution fluorescence microscopy.

Polycation-Anionic Lipid Membrane Interactions

Natural or synthetic polycations are used as biocides or as drug/gene carriers. Understanding the interactions between these macromolecules and cell membranes at the molecular level is therefore of great importance for the design of effective polymer biocides or biocompatible polycation-based delivery systems. Until now, details of the processes at the interface between polycations and biological systems have not been fully recognized. In this study, we consider the effect of strong polycations with quaternary ammonium groups on the properties of anionic lipid membranes that we use as a model system for protein-free cell membranes. For this purpose, we employed experimental measurements and atomic-scale molecular dynamics (MD) simulations. MD simulations reveal that the polycations are strongly hydrated in the aqueous phase and do not lose the water shell after adsorption at the bilayer surface. As a result of strong hydration, the polymer chains reside at the phospholipid headgroup and do not penetrate to the acyl chain region. The polycation adsorption involves the formation of anionic lipid-rich domains, and the density of anionic lipids in these domains depends on the length of the polycation chain. We observed the accumulation of anionic lipids only in the leaflet interacting with the polymer, which leads to the formation of compositionally asymmetric domains. Asymmetric adsorption of the polycation on only one leaflet of the anionic membrane strongly affects the membrane properties in the polycation-membrane contact areas: (i) anionic lipid accumulates in the region near the adsorbed polymer, (ii) acyl chain ordering and lipid packing are reduced, which results in a decrease in the thickness of the bilayer, and (iii) polycation-anionic membrane interactions are strongly influenced by the presence and concentration of salt. Our results provide an atomic-scale description of the interactions of polycations with anionic lipid bilayers and are fully supported by the experimental data. The outcomes are important for understanding the correlation of the structure of polycations with their activity on biomembranes.

Modulated and spiral surface patterns on deformable lipid vesicles

We investigate the behavior of two-dimensional systems that exhibit a transition between homogeneous and spatially inhomogeneous phases, which have spherical topology, and whose mechanical properties depend on the local value of the order parameter. One example of such a system is multicomponent lipid bilayer vesicles, which serve as a model to study cellular membranes. Under certain conditions, such bilayers separate into coexisting liquid-ordered and liquid-disordered regions. When arranged into the shape of small vesicles, this phase coexistence can result in spatial patterns that are more complex than the basic two-domain configuration encountered in typical bulk systems. The difference in bending rigidity between the liquid-ordered and liquid-disordered regions couples the shape of the vesicle to the local composition. We show that this interplay gives rise to a rich phase diagram that includes homogeneous, separated, and axisymmetric modulated phases that are divided by regions of spiral patterns in the surface morphology.

A combined molecular/continuum-modeling approach to predict the small-angle neutron scattering of curved membranes

This paper develops a framework to compute the small-angle neutron scattering (SANS) from highly curved, dynamically fluctuating, and potentially inhomogeneous membranes. This method is needed to compute the scattering from nanometer-scale membrane domains that couple to curvature, as predicted by molecular modeling. The detailed neutron scattering length density of a small planar bilayer patch is readily available via molecular dynamics simulation. A mathematical, mechanical transformation of the planar scattering length density is developed to predict the scattering from curved bilayers. By simulating a fluctuating, curved, surface-continuum model, long time- and length-scales can be reached while, with the aid of the planar-to-curved transformation, the molecular features of the scattering length density can be retained. A test case for the method is developed by constructing a coarse-grained lipid vesicle following a protocol designed to relieve both the osmotic stress inside the vesicle and the lipid-number stress between the leaflets. A question was whether the hybrid model would be able to replicate the scattering from the highly deformed inner and outer leaflets of the small vesicle. Matching the scattering of the full (molecular vesicle) and hybrid (continuum vesicle) models indicated that the inner and outer leaflets of the full vesicle were expanded laterally, consistent with previous simulations of the Martini forcefield that showed thinning in small vesicles. The vesicle structure is inconsistent with a zero-tension leaflet deformed by a single set of elastic parameters, and the results show that this is evident in the scattering. The method can be applied to translate observations of any molecular model's neutron scattering length densities from small patches to large length and timescales.

Microscopic Picture of Calcium-Assisted Lipid Demixing and Membrane Remodeling Using Multiscale Simulations

The specificity of anionic phospholipids-calcium ion interaction and lipid demixing has been established as a key regulatory mechanism in several cellular signaling processes. The mechanism and implications of this calcium-assisted demixing have not been elucidated from a microscopic point of view. Here, we present an overview of atomic interactions between calcium and phospholipids that can drive nonideal mixing of lipid molecules in a model lipid bilayer composed of zwitterionic (1-palmitoyl-2-oleoyl-sn-glycero-3-phosphocholine (POPC)) and anionic (1-palmitoyl-2-oleoyl-sn-glycero-3-phospho-l-serine (POPS)) lipids with computer simulations at multiple resolutions. Lipid nanodomain formation and growth were driven by calcium-enabled lipid bridging of the charged phosphatidylserine (PS) headgroups, which were favored against inter-POPS dipole interactions. Consistent with several experimental studies of calcium-associated membrane sculpting, our analyses also suggest modifications in local membrane curvature and cross-leaflet couplings as a response to such induced lateral heterogeneity. In addition, reverse mapping to a complementary atomistic description revealed structural insights in the presence of anionic nanodomains, at timescales not accessed by previous computational studies. This work bridges information across multiple scales to reveal a mechanistic picture of calcium ion's impact on membrane biophysics.

Simulations of Asymmetric Membranes Illustrate Cooperative Leaflet Coupling and Lipid Adaptability

Biological membranes are composed of lipid bilayers that are often asymmetric with regards to the lipid composition and/or aqueous solvent they separate. Studying lipid asymmetry both experimentally and computationally is challenging. Molecular dynamics simulations of lipid bilayers with asymmetry are difficult due to finite system sizes and time scales accessible to simulations. Due to the very slow flip-flop rate for phospholipids, one must first choose how many lipids are on each side of the bilayer, but the resulting bilayer may be unstable (or metastable) due to differing tensile and compressive forces between leaflets. Here we use molecular dynamics simulations to investigate a number of different asymmetric membrane systems, both with atomistic and coarse-grained models. Asymmetries studied include differences in number of lipids, lipid composition (unsaturated and saturated tails and different headgroups), and chemical gradients between the aqueous phases. Extensive analysis of the bilayers' properties such as area per lipid, density, and lateral pressure profiles are used to characterize bilayer asymmetry. We also address how cholesterol (which flip-flops relatively quickly) influences membrane asymmetries. Our results show how each leaflet is influenced by the other and can mitigate the structural changes to the bilayer overall structure. Cholesterol can respond to changes in bilayer asymmetry to alleviate some of the effect on the bilayer structure, but that will alter its leaflet distribution, which in turn affects its chemical potential. Ionic imbalances are shown to have a modest change in bilayer structure, despite large changes in the electrostatic potential. Bilayer asymmetry can also induce a modest electrostatic potential across the membrane. Our results highlight the importance of membrane asymmetry on bilayer properties, the influence of lipid headgroups, tails and cholesterol on asymmetry, and the ability of lipids to adapt to different environments.

Pivotal Role of Interdigitation in Interleaflet Interactions: Implications from Molecular Dynamics Simulations

The asymmetric lipid composition in plasma membranes within the inner leaflet is not typically suitable for domain formation. Thus elucidation of the likelihood of the formation or stability of a raft-like domain in the inner leaflet is necessary. Herein we investigated the phase behavior of asymmetric membranes using coarse-grained molecular dynamics simulations. The lipid leaflet comprising dioleoylphosphatidylcholine (DOPC) and cholesterol (Chol) does not typically show well-developed domains in symmetric bilayer membranes; however, it does separate into liquid ordered (L_o) and liquid disordered (L_d) phases when the opposing leaflet containing sphingomyelin (SM), DOPC, and Chol demonstrates domain formation. We determine that interdigitated acyl chains modulated the partitioning of Chol in the opposing leaflet, resulting in phase separation. Similarly, the acyl chain length of SM within the opposing leaflet affected the phase behavior of the leaflet. Our results reveal the crucial role of interdigitation in determining the phase status in asymmetric membranes.

Carotenoids promote lateral packing and condensation of lipid membranes

Carotenoids are pigment molecules that protect biomembranes against degradation and may be involved in the formation of functional bacterial membrane microdomains. Little is known on whether different types of carotenoids have different effects on the membrane or if there is any concentration dependence of these effects. In this work, we present results from molecular dynamics simulations of phospholipid bilayers containing different amounts of either β-carotene or zeaxanthin. Both β-carotene and zeaxanthin show the ability to laterally condense the membrane lipids and reduce their inter-leaflet interactions. With increasing concentrations, both carotenoids increase the bilayer thickness and rigidity. The results reveal that carotenoids have similar effects to cholesterol on regulating the behavior of fluid-phase membranes, suggesting that they could function as sterol substitutes and confirming their potential role in the formation of functional membrane domains.

Interleaflet Coupling of Lipid Nanodomains - Insights From in vitro Systems

The plasma membrane is a complex system, consisting of two layers of lipids and proteins compartmentalized into small structures called nanodomains. Despite the asymmetric composition of both leaflets, coupling between the layers is surprisingly strong. This can be evidenced, for example, by recent experimental studies performed on phospholipid giant unilamellar vesicles showing that nanodomains formed in the outer layer are perfectly registered with those in the inner leaflet. Similarly, microscopic phase separation in one leaflet can induce phase separation in the opposing leaflet that would otherwise be homogeneous. In this review, we summarize the current theoretical and experimental knowledge that led to the current view that domains are – irrespective of their size – commonly registered across the bilayer. Mechanisms inducing registration of nanodomains suggested by theory and calculations are discussed. Furthermore, domain coupling is evidenced by experimental studies based on the sparse number of methods that can resolve registered from independent nanodomains. Finally, implications that those findings using model membrane studies might have for cellular membranes are discussed.

Ordered Lipid Domains Assemble via Concerted Recruitment of Constituents from Both Membrane Leaflets

Lipid rafts serve as anchoring platforms for membrane proteins. Thus far they escaped direct observation by light microscopy due to their small size. Here we used differently colored dyes as reporters for the registration of both ordered and disordered lipids from the two leaves of a freestanding bilayer. Photoswitchable lipids dissolved or reformed the domains. Measurements of domain mobility indicated the presence of 120 nm wide ordered and 40 nm wide disordered domains. These sizes are in line with the predicted roles of line tension and membrane undulation as driving forces for alignment.

Elastic deformations mediate interaction of the raft boundary with membrane inclusions leading to their effective lateral sorting

Liquid-ordered lipid domains represent a lateral inhomogeneity in cellular membranes. These domains have elastic and physicochemical properties different from those of the surrounding membrane. In particular, their thickness exceeds that of the disordered membrane. Thus, elastic deformations arise at the domain boundary in order to compensate for the thickness mismatch. In equilibrium, the deformations lead to an incomplete register of monolayer ordered domains: the elastic energy is minimal if domains in opposing monolayers lie on the top of each other, and their boundaries are laterally shifted by about 3 nm. This configuration introduces a region, composed of one ordered and one disordered monolayers, with an intermediate bilayer thickness. Besides, a jump in a local monolayer curvature takes place in this intermediate region, concentrating here most of the elastic stress. This region can participate in a lateral sorting of membrane inclusions by offering them an optimal bilayer thickness and local curvature conditions. In the present study, we consider the sorting of deformable lipid inclusions, undeformable peripheral and deeply incorporated peptide inclusions, and undeformable transmembrane inclusions of different molecular geometry. With rare exceptions, all types of inclusions have an affinity to the ordered domain boundary as compared to the bulk phases. The optimal lateral distribution of inclusions allows relaxing the elastic stress at the boundary of domains.

The Effect of Transmembrane Protein Shape on Surrounding Lipid Domain Formation by Wetting

Signal transduction through cellular membranes requires the highly specific and coordinated work of specialized proteins. Proper functioning of these proteins is provided by an interplay between them and the lipid environment. Liquid-ordered lipid domains are believed to be important players here, however, it is still unclear whether conditions for a phase separation required for lipid domain formation exist in cellular membranes. Moreover, membrane leaflets are compositionally asymmetric, that could be an obstacle for the formation of symmetric domains spanning the lipid bilayer. We theoretically show that the presence of protein in the membrane leads to the formation of a stable liquid-ordered lipid phase around it by the mechanism of protein wetting by lipids, even in the absence of conditions necessary for the global phase separation in the membrane. Moreover, we show that protein shape plays a crucial role in this process, and protein conformational rearrangement can lead to changes in the size and characteristics of surrounding lipid domains.

Role of Transmembrane Proteins for Phase Separation and Domain Registration in Asymmetric Lipid Bilayers

It is well known that the formation and spatial correlation of lipid domains in the two apposed leaflets of a bilayer are influenced by weak lipid–lipid interactions across the bilayer’s midplane. Transmembrane proteins span through both leaflets and thus offer an alternative domain coupling mechanism. Using a mean-field approximation of a simple bilayer-type lattice model, with two two-dimensional lattices stacked one on top of the other, we explore the role of this “structural” inter-leaflet coupling for the ability of a lipid membrane to phase separate and form spatially correlated domains. We present calculated phase diagrams for various effective lipid–lipid and lipid–protein interaction strengths in membranes that contain a binary lipid mixture in each leaflet plus a small amount of added transmembrane proteins. The influence of the transmembrane nature of the proteins is assessed by a comparison with “peripheral” proteins, which result from the separation of one single integral protein into two independent units that are no longer structurally connected across the bilayer. We demonstrate that the ability of membrane-spanning proteins to facilitate domain formation requires sufficiently strong lipid–protein interactions. Weak lipid–protein interactions generally tend to inhibit phase separation in a similar manner for transmembrane as for peripheral proteins.

Intrinsic Curvature-Mediated Transbilayer Coupling in Asymmetric Lipid Vesicles

We measured the effect of intrinsic lipid curvature, J0, on structural properties of asymmetric vesicles made of palmitoyl-oleoyl-phosphatidylethanolamine (POPE; J0<0) and palmitoyl-oleoyl-phosphatidylcholine (POPC; J0∼0). Electron microscopy and dynamic light scattering were used to determine vesicle size and morphology, and x-ray and neutron scattering, combined with calorimetric experiments and solution NMR, yielded insights into leaflet-specific lipid packing and melting processes. Below the lipid melting temperature we observed strong interleaflet coupling in asymmetric vesicles with POPE inner bilayer leaflets and outer leaflets enriched in POPC. This lipid arrangement manifested itself by lipids melting cooperatively in both leaflets, and a rearrangement of lipid packing in both monolayers. On the other hand, no coupling was observed in vesicles with POPC inner bilayer leaflets and outer leaflets enriched in POPE. In this case, the leaflets melted independently and did not affect each other’s acyl chain packing. Furthermore, we found no evidence for transbilayer structural coupling above the melting temperature of either sample preparation. Our results are consistent with the energetically preferred location of POPE residing in the inner leaflet, where it also resides in natural membranes, most likely causing the coupling of both leaflets. The loss of this coupling in the fluid bilayers is most likely the result of entropic contributions.

Computational Modeling of Realistic Cell Membranes

Cell membranes contain a large variety of lipid types and are crowded with proteins, endowing them with the plasticity needed to fulfill their key roles in cell functioning. The compositional complexity of cellular membranes gives rise to a heterogeneous lateral organization, which is still poorly understood. Computational models, in particular molecular dynamics simulations and related techniques, have provided important insight into the organizational principles of cell membranes over the past decades. Now, we are witnessing a transition from simulations of simpler membrane models to multicomponent systems, culminating in realistic models of an increasing variety of cell types and organelles. Here, we review the state of the art in the field of realistic membrane simulations and discuss the current limitations and challenges ahead.

Multiscale Simulations of Biological Membranes: The Challenge To Understand Biological Phenomena in a Living Substance

Biological membranes are tricky to investigate. They are complex in terms of molecular composition and structure, functional over a wide range of time scales, and characterized by nonequilibrium conditions. Because of all of these features, simulations are a great technique to study biomembrane behavior. A significant part of the functional processes in biological membranes takes place at the molecular level; thus computer simulations are the method of choice to explore how their properties emerge from specific molecular features and how the interplay among the numerous molecules gives rise to function over spatial and time scales larger than the molecular ones. In this review, we focus on this broad theme. We discuss the current state-of-the-art of biomembrane simulations that, until now, have largely focused on a rather narrow picture of the complexity of the membranes. Given this, we also discuss the challenges that we should unravel in the foreseeable future. Numerous features such as the actin-cytoskeleton network, the glycocalyx network, and nonequilibrium transport under ATP-driven conditions have so far received very little attention; however, the potential of simulations to solve them would be exceptionally high. A major milestone for this research would be that one day we could say that computer simulations genuinely research biological membranes, not just lipid bilayers.

Understanding the Role of Lipids in Signaling Through Atomistic and Multiscale Simulations of Cell Membranes

Cell signaling controls essentially all cellular processes. While it is often assumed that proteins are the key architects coordinating cell signaling, recent studies have shown more and more clearly that lipids are also involved in signaling processes in a number of ways. Lipids do, for instance, act as messengers, modulate membrane receptor conformation and dynamics, and control membrane receptor partitioning. Further, through structural modifications such as oxidation, the functions of lipids as part of signaling processes can be modified. In this context, in this article we discuss the understanding recently revealed by atomistic and coarse-grained computer simulations of nanoscale processes and underlying physicochemical principles related to lipids’ functions in cellular signaling.

Experimental Evidence of the Existence of Interleaflet Coupled Nanodomains: An MC-FRET Study

Plasma membranes of living cells are compartmentalized into small submicroscopic structures (nanodomains) having potentially relevant biological functions. Despite this, structural features of these nanodomains remain elusive, primarily due to the difficulties in characterizing such small dynamic entities. It is unclear whether nanodomains found in the upper bilayer leaflet are transversally registered with those found in the lower leaflet. Experiments performed on larger microscopic domains indicate that the coupling between the leaflets is strong, forcing the domains to be in perfect registration, but can the same thing be said about the biologically more relevant nanodomains? This work provides experimental evidence that even small nanodomains of variable sizes between 10 and 160 nm are interleaflet coupled. Importantly, the alternative scenarios of partially registered, independent, or antiregistered nanodomains could be excluded.

Molecular Dynamics Simulations Reveal Leaflet Coupling in Compositionally Asymmetric Phase-Separated Lipid Membranes

The eukaryotic plasma membrane has an asymmetric distribution of its component lipids. Rafts that result from liquid-liquid phase separation are a feature of its exoplasmic leaflet, but how these exoplasmic leaflet domains are coupled to the cytoplasmic leaflet is not understood. These rafts can be studied in model membranes of three-component mixtures that produce coexisting liquid ordered (Lo) and liquid disordered (Ld) domains. We conducted all-atom molecular dynamics simulations of compositionally asymmetric lipid bilayers that reflect a more realistic model of the plasma membrane. One leaflet contained phase-separated domains with phosphatidylcholine and cholesterol, representing the exoplasmic leaflet, whereas the other contained phosphatidylethanolamine, phosphatidylserine, and cholesterol, which are the predominant components of the cytoplasmic leaflet. Inspired by findings of domain alignment across the two leaflets in compositionally symmetric model membranes, we examined the coupling between the two leaflets to see how the single-phase cytoplasmic leaflet would respond to phase separation in the other leaflet and if information could be communicated across the membrane. We found the region of the single-phase leaflet apposing the Lo domain to be slightly more ordered and thicker than the region apposing the Ld domain. The region across from the Lo domain is somewhat enriched in cholesterol and significantly depleted of polyunsaturated lipids.

Gramicidin Increases Lipid Flip-Flop in Symmetric and Asymmetric Lipid Vesicles

Unlike most transmembrane proteins, phospholipids can migrate from one leaflet of the membrane to the other. Because this spontaneous lipid translocation (flip-flop) tends to be very slow, cells facilitate the process with enzymes that catalyze the transmembrane movement and thereby regulate the transbilayer lipid distribution. Nonenzymatic membrane-spanning proteins with unrelated primary functions have also been found to accelerate lipid flip-flop in a nonspecific manner and by various hypothesized mechanisms. Using deuterated phospholipids, we examined the acceleration of flip-flop by gramicidin channels, which have well-defined structures and known functions, features that make them ideal candidates for probing the protein-membrane interactions underlying lipid flip-flop. To study compositionally and isotopically asymmetric proteoliposomes containing gramicidin, we expanded a recently developed protocol for the preparation and characterization of lipid-only asymmetric vesicles. Channel incorporation, conformation, and function were examined with small angle x-ray scattering, circular dichroism, and a stopped-flow spectrofluorometric assay, respectively. As a measure of lipid scrambling, we used differential scanning calorimetry to monitor the effect of gramicidin on the melting transition temperatures of the two bilayer leaflets. The two calorimetric peaks of the individual leaflets merged into a single peak over time, suggestive of scrambling, and the effect of the channel on the transbilayer lipid distribution in both symmetric 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphocholine and asymmetric 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphocholine/1,2-dimyristoyl-sn-glycero-3-phosphocholine vesicles was quantified from proton NMR measurements. Our results show that gramicidin increases lipid flip-flop in a complex, concentration-dependent manner. To determine the molecular mechanism of the process, we used molecular dynamics simulations and further computational analysis of the trajectories to estimate the extent of membrane deformation. Together, the experimental and computational approaches were found to constitute an effective means for studying the effects of transmembrane proteins on lipid distribution in both symmetric and asymmetric model membranes.

SPICA Force Field for Lipid Membranes: Domain Formation Induced by Cholesterol

Heterogeneity is essential for multicomponent lipid membranes. Especially, sterol-induced domain formation in membranes has recently attracted attention because of its biological importance. To investigate such membrane domains at the molecular level, coarse-grained molecular dynamics (CG-MD) simulations are a promising approach since they allow one to consider the temporal and spatial scales involved in domain formation. In this work, we present a new CG force field, named SPICA, which can accurately predict domain formation within various lipids in membranes. The SPICA force field was developed as an extension of a previous CG model, known as SDK (Shinoda-DeVane-Klein), in which membrane properties such as tension, elasticity, and structure are well reproduced. By examining domain formation in a series of ternary lipid bilayers, we observed a separation into liquid-ordered and liquid-disordered phases fully consistent with experimental observations. Importantly, it is shown that the SPICA force field can detect the different phase behavior that results from subtle differences in the lipid composition of the bilayer.