Line Tension Controls Liquid-Disordered + Liquid-Ordered Domain Size Transition in Lipid Bilayers

Membrane Bending Moduli of Coexisting Liquid Phases Containing Transmembrane Peptide

A number of highly curved membranes in vivo, such as epithelial cell microvilli, have the relatively high sphingolipid content associated with "raft-like" composition. Given the much lower bending energy measured for bilayers with "nonraft" low sphingomyelin and low cholesterol content, observing high curvature for presumably more rigid compositions seems counterintuitive. To understand this behavior, we measured membrane rigidity by fluctuation analysis of giant unilamellar vesicles. We found that including a transmembrane helical GWALP peptide increases the membrane bending modulus of the liquid-disordered (Ld) phase. We observed this increase at both low-cholesterol fraction and higher, more physiological cholesterol fraction. We find that simplified, commonly used Ld and liquid-ordered (Lo) phases are not representative of those that coexist. When Ld and Lo phases coexist, GWALP peptide favors the Ld phase with a partition coefficient of 3-10 depending on mixture composition. In model membranes at high cholesterol fractions, Ld phases with GWALP have greater bending moduli than the Lo phase that would coexist.

Exploring the impact of proteins on the line tension of a phase-separating ternary lipid mixture

The separation of lipid mixtures into thermodynamically stable phase-separated domains is dependent on lipid composition, temperature, and system size. Using molecular dynamics simulations, the line tension between thermodynamically stable lipid domains formed from ternary mixtures of di-C16:0 PC:di-C18:2 PC:cholesterol at 40:40:20 mol. % ratio was investigated via two theoretical approaches. The line tension was found to be 3.1 ± 0.2 pN by capillary wave theory and 4.7 ± 3.7 pN by pressure tensor anisotropy approaches for coarse-grained models based on the Martini force field. Using an all-atom model of the lipid membrane based on the CHARMM36 force field, the line tension was found to be 3.6 ± 0.9 pN using capillary wave theory and 1.8 ± 2.2 pN using pressure anisotropy approaches. The discrepancy between estimates of the line tension based on capillary wave theory and pressure tensor anisotropy methods is discussed. Inclusion of protein in Martini membrane lipid mixtures was found to reduce the line tension by 25%-35% as calculated by the capillary wave theory approach. To further understand and predict the behavior of proteins in phase-separated membranes, we have formulated an analytical Flory-Huggins model and parameterized it against the simulation results. Taken together these results suggest a general role for proteins in reducing the thermodynamic cost associated with domain formation in lipid mixtures and quantifies the thermodynamic driving force promoting the association of proteins to domain interfaces.

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.

Scaling of lipid membrane rigidity with domain area fraction

Biological membranes are highly heterogeneous in composition which in turn leads to local variations in the physical properties. Here we quantify how heterogeneity in stiffness determines the effective bending modulus, κ_eff, of model phase-separated membranes with coexisting soft fluid and rigid gel domains. We find that the temperature- and composition- dependent trends in membrane rigidity collapse onto a single curve, such that κ_eff directly scales with the area fraction of the rigid gel domains. Using no adjustable parameters, the measurements are found to agree with theoretical predictions for inhomogenous membranes and indicate that κ_eff is sensitive to the lateral distribution of the rigid phase within the membrane. This key finding confirms that the properties of heterogeneous membranes can be quantitatively predicted if the area fraction and properties of the individual phases are known.

Lowering line tension with high cholesterol content induces a transition from macroscopic to nanoscopic phase domains in model biomembranes

Chemically simplified lipid mixtures are used here as models of the cell plasma membrane exoplasmic leaflet. In such models, phase separation and morphology transitions controlled by line tension in the liquid-disordered (Ld) + liquid-ordered (Lo) coexistence regime have been described [1]. Here, we study two four-component lipid mixtures at different cholesterol fractions: brain sphingomyelin (BSM) or 1,2-distearoyl-sn-glycero-3-phosphocholine (DSPC)/1,2-dioleoyl-sn-glycero-3-phosphocholine (DOPC)/1-palmitoyl-2-oleoyl-sn-glycero-3-phosphocholine (POPC)/cholesterol (Chol). On giant unilamellar vesicles (GUVs) display a nanoscopic-to-macroscopic transition of Ld + Lo phase domains as POPC is replaced by DOPC, and this transition also depends on the cholesterol fraction. Line tension decreases with increasing cholesterol mole fractions in both lipid mixtures. For the ternary BSM/DOPC/Chol mixture, the published phase diagram [19] requires a modification to show that when cholesterol mole fraction is >~0.33, coexisting phase domains become nanoscopic.

The Alzheimer's disease amyloid-β peptide affects the size-dynamics of raft-mimicking Lo domains in GM1-containing lipid bilayers

Alzheimer's disease (AD) is characterized by the overproduction of the amyloid-β peptide (Aβ) which forms fibrils under the influence of raft microdomains containing the ganglioside GM1. Raft-mimicking artificial liquid ordered (Lo) domains containing GM1 enhance amyloid-β polymerization. Other experiments suggest that Aβ binds preferably to the non-raft liquid disordered (Ld) phase rather than to the Lo phase in the presence of GM1. Here, the interaction of Aβ(1-42) with GM1-containing biphasic Lo-Ld giant vesicles was investigated. Fluorescence colocalisation experiments confirm that Aβ(1-42) binds preferentially to the Ld phase. The effect of Aβ(1-42) on the Lo-Ld size dynamics was studied using photoinduced spinodal decomposition which mimics the nanodomain-microdomain raft coalescence. Aβ affects the kinetics of the coarsening phase and the size of the resulting microdomains. The effect depends on which phase is in a majority: when the Lo microdomains are formed inside an Ld phase, their growth rate becomes slower and their final size smaller in the presence of Aβ(1-42), whereas when the Ld microdomains are formed inside an Lo phase, the growth rate becomes faster and the final size larger. Fluorimetric measurements on large vesicles using the probe Laurdan indicate that Aβ(1-42) binding respectively increases or decreases the packing of the Ld phase in the presence or absence of GM1. The differential effects of Aβ on spinodal decomposition are accordingly interpreted as resulting from distinct effects of the peptide on the Lo-Ld line tension modulated by GM1. Such modulating effect of Aβ on domain dynamics could be important for lipid rafts in signaling disorders in AD as well as in Aβ fibrillation.

Regimes of Complex Lipid Bilayer Phases Induced by Cholesterol Concentration in MD Simulation

Cholesterol is essential to the formation of phase-separated lipid domains in membranes. Lipid domains can exist in different thermodynamic phases depending on the molecular composition and play significant roles in determining structure and function of membrane proteins. We investigate the role of cholesterol in the structure and dynamics of ternary lipid mixtures displaying phase separation using molecular dynamics simulations, employing a physiologically relevant span of cholesterol concentration. We find that cholesterol can induce formation of three regimes of phase behavior: 1) miscible liquid-disordered bulk, 2) phase-separated, domain-registered coexistence of liquid-disordered and liquid-ordered domains, and 3) phase-separated, domain-antiregistered coexistence of liquid-disordered and newly identified nanoscopic gel domains composed of cholesterol threads we name "cholesterolic gel" domains. These findings are validated and discussed in the context of current experimental knowledge, models of cholesterol spatial distributions, and models of ternary lipid-mixture phase separation.

A Rationale for Mesoscopic Domain Formation in Biomembranes

Cell plasma membranes display a dramatically rich structural complexity characterized by functional sub-wavelength domains with specific lipid and protein composition. Under favorable experimental conditions, patterned morphologies can also be observed in vitro on model systems such as supported membranes or lipid vesicles. Lipid mixtures separating in liquid-ordered and liquid-disordered phases below a demixing temperature play a pivotal role in this context. Protein-protein and protein-lipid interactions also contribute to membrane shaping by promoting small domains or clusters. Such phase separations displaying characteristic length-scales falling in-between the nanoscopic, molecular scale on the one hand and the macroscopic scale on the other hand, are named mesophases in soft condensed matter physics. In this review, we propose a classification of the diverse mechanisms leading to mesophase separation in biomembranes. We distinguish between mechanisms relying upon equilibrium thermodynamics and those involving out-of-equilibrium mechanisms, notably active membrane recycling. In equilibrium, we especially focus on the many mechanisms that dwell on an up-down symmetry breaking between the upper and lower bilayer leaflets. Symmetry breaking is an ubiquitous mechanism in condensed matter physics at the heart of several important phenomena. In the present case, it can be either spontaneous (domain buckling) or explicit, i.e., due to an external cause (global or local vesicle bending properties). Whenever possible, theoretical predictions and simulation results are confronted to experiments on model systems or living cells, which enables us to identify the most realistic mechanisms from a biological perspective.

Tuning Length Scales of Small Domains in Cell-Derived Membranes and Synthetic Model Membranes

Micron-scale, coexisting liquid-ordered (L_o) and liquid-disordered (L_d) phases are straightforward to observe in giant unilamellar vesicles (GUVs) composed of ternary lipid mixtures. Experimentally, uniform membranes undergo demixing when temperature is decreased: domains subsequently nucleate, diffuse, collide, and coalesce until only one domain of each phase remains. The sizes of these two domains are limited only by the size of the system. Under different conditions, vesicles exhibit smaller-scale domains of fixed sizes, leading to the question of what sets the length scale. In membranes with excess area, small domains are expected when coarsening is hindered or when a microemulsion or modulated phase arises. Here, we test predictions of how the size, morphology, and fluorescence levels of small domains vary with the membrane's temperature, tension, and composition. Using GUVs and cell-derived giant plasma membrane vesicles, we find that 1) the characteristic size of domains decreases when temperature is increased or membrane tension is decreased, 2) stripes are favored over circular domains for lipid compositions with low energy per unit interface, 3) fluorescence levels are consistent with domain registration across both monolayer leaflets of the bilayer, and 4) small domains form in GUVs composed of lipids both with and without ester-linked lipids. Our experimental results are consistent with several elements of current theories for microemulsions and modulated phases and inconsistent with others, suggesting a motivation to modify or enhance current theories.

Docosahexaenoic acid regulates the formation of lipid rafts: A unified view from experiment and simulation

Docosahexaenoic acid (DHA, 22:6) is an n-3 polyunsaturated fatty acid (n-3 PUFA) that influences immunological, metabolic, and neurological responses through complex mechanisms. One structural mechanism by which DHA exerts its biological effects is through its ability to modify the physical organization of plasma membrane signaling assemblies known as sphingomyelin/cholesterol (SM/chol)-enriched lipid rafts. Here we studied how DHA acyl chains esterified in the sn-2 position of phosphatidylcholine (PC) regulate the formation of raft and non-raft domains in mixtures with SM and chol on differing size scales. Coarse grained molecular dynamics simulations showed that 1-palmitoyl-2-docosahexaenoylphosphatylcholine (PDPC) enhances segregation into domains more than the monounsaturated control, 1-palmitoyl-2-oleoyl-phosphatidylcholine (POPC). Solid state ²H NMR and neutron scattering experiments provided direct experimental evidence that substituting PDPC for POPC increases the size of raft-like domains on the nanoscale. Confocal imaging of giant unilamellar vesicles with a non-raft fluorescent probe revealed that POPC had no influence on phase separation in the presence of SM/chol whereas PDPC drove strong domain segregation. Finally, monolayer compression studies suggest that PDPC increases lipid-lipid immiscibility in the presence of SM/chol compared to POPC. Collectively, the data across model systems provide compelling support for the emerging model that DHA acyl chains of PC lipids tune the size of lipid rafts, which has potential implications for signaling networks that rely on the compartmentalization of proteins within and outside of rafts.

Frontiers of Complex Disease Mechanisms: Membrane Surface Tension May Link Genotype to Phenotype in Glaucoma

Although many monogenic diseases are understood based upon structural changes of gene products, less progress has been made concerning polygenic disease mechanisms. This article presents a new interdisciplinary approach to understand complex diseases, especially their genetic polymorphisms. I focus upon primary open angle glaucoma (POAG). Although elevated intraocular pressure (IOP) and oxidative stress are glaucoma hallmarks, the linkages between these factors and cell death are obscure. Reactive oxygen species (ROS) promote the formation of oxidatively truncated phosphoglycerides (OTP), free fatty acids, lysophosphoglycerides, oxysterols, and other chemical species that promote membrane disruption and decrease membrane surface tension. Several POAG-linked gene polymorphisms identify proteins that manage damaged lipids and/or influence membrane surface tension. POAG-related genes expected to participate in these processes include: ELOVL5, ABCA1, APOE4, GST, CYP46A1, MYOC, and CAV. POAG-related gene products are expected to influence membrane surface tension, strength, and repair. I propose that heightened IOP overcomes retinal ganglion cell (RGC) membrane compressive strength, weakened by damaged lipid accumulation, to form pores. The ensuing structural failure promotes apoptosis and blindness. The linkage between glaucoma genotype and phenotype is mediated by physical events. Force balancing between the IOP and compressive strength regulates pore nucleation; force balancing between pore line tension and membrane surface tension regulates pore growth. Similar events may contribute to traumatic brain injury, Alzheimer's disease, and macular degeneration.

Membrane perturbing properties of toxin mycolactone from Mycobacterium ulcerans

Mycolactone is the exotoxin produced by Mycobacterium ulcerans and is the virulence factor behind the neglected tropical disease Buruli ulcer. The toxin has a broad spectrum of biological effects within the host organism, stemming from its interaction with at least two molecular targets and the inhibition of protein uptake into the endoplasmic reticulum. Although it has been shown that the toxin can passively permeate into host cells, it is clearly lipophilic. Association with lipid carriers would have substantial implications for the toxin’s distribution within a host organism, delivery to cellular targets, diagnostic susceptibility, and mechanisms of pathogenicity. Yet the toxin’s interactions with, and distribution in, lipids are unknown. Herein we have used coarse-grained molecular dynamics simulations, guided by all-atom simulations, to study the interaction of mycolactone with pure and mixed lipid membranes. Using established techniques, we calculated the toxin’s preferential localization, membrane translocation, and impact on membrane physical and dynamical properties. The computed water-octanol partition coefficient indicates that mycolactone prefers to be in an organic phase rather than in an aqueous environment. Our results show that in a solvated membrane environment the exotoxin mainly localizes in the water-membrane interface, with a preference for the glycerol moiety of lipids, consistent with the reported studies that found it in lipid extracts of the cell. The calculated association constant to the model membrane is similar to the reported association constant for Wiskott-Aldrich syndrome protein. Mycolactone is shown to modify the physical properties of membranes, lowering the transition temperature, compressibility modulus, and critical line tension at which pores can be stabilized. It also shows a tendency to behave as a linactant, a molecule that localizes at the boundary between different fluid lipid domains in membranes and promotes inter-mixing of domains. This property has implications for the toxin’s cellular access, T-cell immunosuppression, and therapeutic potential.

Mycolactone is a macrolide exotoxin secreted by Mycobacterium ulcerans, which causes a skin disease called Buruli ulcer, a neglected emerging disease. It is the third most common mycobacterial disease after tuberculosis and leprosy. Studies have shown how mycolactone plays a pivotal role in Buruli ulcer pathogenesis, and identified it as an attractive therapeutic target. This multifunctional cytotoxin exerts multiple local and global responses, including ulcerative, analgesic, and anti-inflammatory effects. Prompted by its lipid-like structure, we used extensive multi-resolution simulations to probe mycolactone’s interactions with model membranes. Our results suggest that mycolactone is sequestered in membranes where it alters several dynamical, physical, and mechanical properties. It also behaves as a linactant, localizing at the interface between lipid domains and decreasing the inter-domain line tension. Our results shed light on how mycolactone permeates host cell membranes and is distributed between lipid and aqueous environments. These findings have significant implications for the toxin’s distribution in the host environment and mechanisms of pathogenicity. Understanding the toxin’s distribution and mechanism of trafficking will have ramifications for targeted diagnostics, therapeutic approaches, and our understanding of Buruli ulcer pathogenesis.

Composition Fluctuations in Lipid Bilayers

Cell membranes contain multiple lipid and protein components having heterogeneous in-plane (lateral) distribution. Nanoscale rafts are believed to play an important functional role, but their phase state—domains of coexisting phases or composition fluctuations—is unknown. As a step toward understanding lateral organization of cell membranes, we investigate the difference between nanoscale domains of coexisting phases and composition fluctuations in lipid bilayers. We simulate model lipid bilayers with the MARTINI coarse-grained force field on length scales of tens of nanometers and timescales of tens of microseconds. We use a binary and a ternary mixture: a saturated and an unsaturated lipid, or a saturated lipid, an unsaturated lipid, and cholesterol, respectively. In these mixtures, the phase behavior can be tuned from a mixed state to a coexistence of a liquid-crystalline and a gel, or a liquid-ordered and a liquid-disordered phase. Transition from a two-phase to a one-phase state is achieved by raising the temperature and adding a hybrid lipid (with a saturated and an unsaturated chain). We analyze the evolution of bilayer properties along this transition: domains of two phases transform to fluctuations with local ordering and compositional demixing. Nanoscale domains and fluctuations differ in several properties, including interleaflet overlap and boundary length. Hybrid lipids show no enrichment at the boundary, but decrease the difference between the coexisting phases by ordering the disordered phase, which could explain their role in cell membranes.

Complex biomembrane mimetics on the sub-nanometer scale

Biomimetic lipid vesicles are indispensable tools for gaining insight into the biophysics of cell physiology on the molecular level. The level of complexity of these model systems has steadily increased, and now spans from domain-forming lipid mixtures to asymmetric lipid bilayers. Here, we review recent progress in the development and application of elastic neutron and X-ray scattering techniques for studying these systems in situ and under physiologically relevant conditions on the nanometer to sub-nanometer length scales. In particular, we focus on: (1) structural details of coexisting liquid-ordered and liquid-disordered domains, including their thickness and lipid packing mismatch as a function of a size transition from nanoscopic to macroscopic domains; (2) membrane-mediated protein partitioning into lipid domains; (3) the role of the aqueous medium in tuning interactions between membranes and domains; and (4) leaflet-specific structure in asymmetric bilayers and passive lipid flip-flop.