How cholesterol is distributed between monolayers in asymmetric lipid membranes

Thermal Stability of bivalent cation/phosphoinositide domains in model membranes

As key mediators in a wide array of signaling events, phosphoinositides (PIPs) orchestrate the recruitment of proteins to specific cellular locations at precise moments. This intricate spatiotemporal regulation of protein activity often necessitates the localized enrichment of the corresponding PIP. We investigate the extent and thermal stabilities of phosphatidylinositol-4-phosphate (PI(4)P), phosphatidylinositol-4,5-bisphosphate (PI(4,5)P2 and phosphatidylinositol-3,4,5-trisphosphate (PI(3,4,5)P3) clusters with calcium and magnesium ions. We observe negligible or minimal clustering of all examined PIPs in the presence of Mg2+ ions. While PI(4)P shows in the presence of Ca2+ no clustering, PI(4,5)P2 forms with Ca2+ strong clusters that exhibit stablity up to at least 80°C. The extent of cluster formation for the interaction of PI(3,4,5)P3 with Ca2+ is less than what was observed for PI(4,5)P2, yet we still observe some clustering up to 80°C. Given that cholesterol has been demonstrated to enhance PIP clustering, we examined whether bivalent cations and cholesterol synergistically promote PIP clustering. We found that the interaction of Mg2+ or Ca2+ with PI(4)P remains extraordinarily weak, even in the presence of cholesterol. In contrast, we observe synergistic interaction of cholesterol and Ca2+ with PI(4,5)P2. Also, in the presence of cholesterol, the interaction of Mg2+ with PI(4,5)P2 remains weak. PI(3,4,5)P3 does not show strong clustering with cholesterol for the experimental conditions of our study and the interaction with Ca2+ and Mg2+ was not influenced by the presence of cholesterol.

Biomembranes balance many types of leaflet asymmetries

Many biological membranes host different lipid species in their two leaflets. Since their spontaneous curvatures are typically not the same, this compositional asymmetry generally entails bending torques, which can be counteracted by differential stress-the difference between the two leaflet tensions. This stress, in turn, can affect elastic parameters or phase behavior of the membrane or each individual leaflet, or push easily flippable species, especially cholesterol, from the compressed leaflet into the tense leaflet. In short, breaking the symmetry of a single observable (to wit: composition), essentially breaks all other symmetries as well, with many potentially interesting consequences. This brief report examines the elastic aspects of this interplay, focusing on some elementary conditions of mechanical and thermodynamic equilibrium, but also shows how this poses novel questions that we are only beginning to appreciate.

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.

Interdependence of cholesterol distribution and conformational order in lipid bilayers

We show, via molecular simulations, that not only does cholesterol induce a lipid order, but the lipid order also enhances cholesterol localization within the lipid leaflets. Therefore, there is a strong interdependence between these two phenomena. In the ordered phase, cholesterol molecules are predominantly present in the bilayer leaflets and orient themselves parallel to the bilayer normal. In the disordered phase, cholesterol molecules are mainly present near the center of the bilayer at the midplane region and are oriented orthogonal to the bilayer normal. At the melting temperature of the lipid bilayers, cholesterol concentration in the leaflets and the bilayer midplane is equal. This result suggests that the localization of cholesterol in the lipid bilayers is mainly dictated by the degree of ordering of the lipid bilayer. We validate our findings on 18 different lipid bilayer systems, obtained from three different phospholipid bilayers with varying concentrations of cholesterol. To cover a large temperature range in simulations, we employ the Dry Martini force field. We demonstrate that the Dry and the Wet Martini (with polarizable water) force fields produce comparable results.

Quantitative relationship between cholesterol distribution and ordering of lipids in asymmetric lipid bilayers

The plasma membrane of eukaryotic cells is known to be compositionally asymmetric. Certain phospholipids, such as sphingomyelin and phosphatidylcholine species, are predominantly localized in the outer leaflet, while phosphatidylethanolamine and phosphatidylserine species primarily reside in the inner leaflet. While phospholipid asymmetry between the membrane leaflets is well established, there is no consensus about cholesterol distribution between the two leaflets. We have performed a systematic study, via molecular simulations, of how the spatial distribution of cholesterol molecules in different "asymmetric" lipid bilayers are affected by the lipids' backbone, head-type, unsaturation, and chain-length by considering an asymmetric bilayer mimicking the plasma membrane lipids of red blood cells, as well as seventeen other asymmetric bilayers comprising of different lipid types. Our results reveal that the distribution of cholesterol in the leaflets is solely a function of the extent of ordering of the lipids within the leaflets. The ratio of the amount of cholesterol matches the ratio of lipid order in the two leaflets, thus providing a quantitative relationship between the two. These results are understood by the observation that asymmetric bilayers with equimolar amount of lipids in the two leaflets develop tensile and compressive stresses due to differences in the extent of lipid order. These stresses are alleviated by the transfer of cholesterol from the leaflet in compressive stress to the one in tensile stress. These findings are important in understanding the biology of the cell membrane, especially with regard to the composition of the membrane leaflets.

Long chain sphingomyelin depletes cholesterol from the cytoplasmic leaflet in asymmetric lipid membranes

The transbilayer distribution of cholesterol (CHL) in complex asymmetric lipid membranes remains controversial, with contrasting investigations suggesting that there is more CHL either in the exoplasmic, outer leaflet (OL) or the cytoplasmic, inner leaflet (IL) depending on cell type or model, membrane composition, and method of investigation. Here, we launch systematic coarse-grained molecular dynamics simulations to investigate the impact of the sphingomyelin (SM) acyl chain length upon CHL distribution in asymmetric lipid membrane mixtures which account for the variation of the most abundant headgroups and acyl chain unsaturation in the two membrane leaflets. We find that there is always more CHL in the OL, but longer chain SM depletes more CHL from the IL than short chain SM in simple membrane mixtures containing SM and 16 : 0, 18 : 1 phospholipids. The difference between longer and shorter chain SM is neutralised in a more complex asymmetric membrane, where there are more saturated tails in the outer leaflet. We propose that interdigitation of long-chain SM into the opposing IL pushes cytoplasmic CHL towards the OL, but higher chain saturation of the outer leaflet compensates for the effect of SM chain length.

Long acyl chain sphingomyelin and saturated phospholipid tails in the outer membrane leaflet deplete cholesterol from the inner leaflet in mammalian membranes.

Thermal conductivity of the cell membrane in the presence of cholesterol and amyloid precursor protein

The cell membrane is responsible for the transportation of heat between inside and outside the cell. Whether the thermal properties of the cell membrane are affected by the cholesterol concentration or the membrane proteins has not been investigated so far. Although the experimental measurement of the membrane thermal conductivity was not available until very recently, computational methods have been widely used for this purpose. In this study, we carry out molecular dynamics simulations to investigate the relation between the concentration of cholesterol and the thermal conductivity of a model membrane. Our results suggest an increase in the membrane thermal conductivity upon increasing the concentration of cholesterol in the membrane. Moreover, we find that the asymmetric distribution of cholesterol in the two membrane leaflets decreases thermal conductivity. We also find a rectification effect when heat flows in opposite directions through a model membrane decorated with the amyloid precursor protein. The results of this study apply to the advancement of selective treatment methods, as well as the development of new materials such as biological rectifiers.

The ABCs of Sterol Transport

Cholesterol homeostasis and trafficking are critical to the maintenance of the asymmetric plasma membrane of eukaryotic cells. Disruption or dysfunction of cholesterol trafficking leads to numerous human diseases. ATP-binding cassette (ABC) transporters play several critical roles in this process, and mutations in these sterol transporters lead to disorders such as Tangier disease and sitosterolemia. Biochemical and structural information on ABC sterol transporters is beginning to emerge, with published structures of ABCA1 and ABCG5/G8; these two proteins function in the reverse cholesterol transport pathway and mediate the efflux of cholesterol and xenosterols to high-density lipoprotein and bile salt micelles, respectively. Although both of these transporters belong to the ABC family and mediate the efflux of a sterol substrate, they have many distinct differences. Here, we summarize the current understanding of sterol transport driven by ABC transporters, with an emphasis on these two extensively characterized transporters.

Probing cholesterol binding and translocation in P-glycoprotein

P-glycoprotein (Pgp) is a biomedically important member of the ABC transporter superfamily that mediates multidrug resistance in various cancer types. Substrate binding and transport in Pgp are modulated by the presence of cholesterol in the membrane. Structural information on cholesterol binding sites and mechanistic details of its redistribution are, however, largely unknown. In this study, a set of 40 independent molecular dynamics (MD) simulations of Pgp embedded in cholesterol-rich lipid bilayers are reported, totaling 8 μs, enabling extensive sampling of cholesterol-protein interactions in Pgp. Clustering analyses of the ensemble of cholesterol molecules (∼5740) sampled around Pgp in these simulations reveal specific and asymmetric cholesterol-binding regions formed by the transmembrane (TM) helices TM1-6 and TM8. Notably, not all the putative cholesterol binding sites identified by MD can be predicted by the primary sequence based cholesterol-recognition amino acid consensus (CRAC) or inverted CRAC (CARC) motifs, an observation that we attribute to inadequacy of these motifs to account for binding sites formed by remote amino acids in the sequence that can still be spatially adjacent to each other. Binding of cholesterol to Pgp occurs more frequently through its rough β-face formed by the two protruding methyl groups, whereas the opposite smooth α-face prefers packing alongside the membrane lipids. One full and two partial cholesterol flipping events between the two leaflets of the bilayer mediated by the surface of Pgp are also captured in these simulations. All flipping events are observed in a region formed by helices TM1, TM2, and TM11, featuring two full and two partial CRAC/CARC motifs, with Tyr49 and Tyr126 identified as key residues interacting with cholesterol during this event. Our study is the first to report direct observation of unconventional cholesterol translocation on the surface of Pgp, providing a secondary transport model for the known flippase activity of ABC exporters of cholesterol. This article is part of a Special Issue entitled: Molecular biophysics of membranes and membrane proteins.

Curvature increases permeability of the plasma membrane for ions, water and the anti-cancer drugs cisplatin and gemcitabine

In this work the permeability of a model asymmetric plasma membrane, for ions, water and the anti-cancer drugs cisplatin and gemcitabine is studied by means of all-atom molecular dynamics simulations. It is shown for the first time that permeability of the highly curved membrane increases from one to three orders of magnitude upon membrane bending depending on the compound and the sign of curvature. Our results suggest that the membrane curvature could be an important factor of drug translocation through the membrane.

Raman Spectroscopy Study of Curvature-Mediated Lipid Packing and Sorting in Single Lipid Vesicles

Cellular plasma membrane deformability and stability is important in a range of biological processes. Changes in local curvature of the membrane affect the lateral movement of lipids, affecting the biophysical properties of the membrane. An integrated holographic optical tweezers and Raman microscope was used to investigate the effect of curvature gradients induced by optically stretching individual giant unilamellar vesicles (GUVs) on lipid packing and lateral segregation of cholesterol in the bilayer. The spatially resolved Raman analysis enabled detection of induced phase separation and changes in lipid ordering in individual GUVs. Using deuterated cholesterol, the changes in lipid ordering and phase separation were linked to lateral sorting of cholesterol in the stretched GUVs. Stretching the GUVs in the range of elongation factors 1-1.3 led to an overall decrease in cholesterol concentration at the edges compared to the center of stretched GUVs. The Raman spectroscopy results were consistent with a model of the bilayer accounting for cholesterol sorting in both bilayer leaflets, with a compositional asymmetry of 0.63 ± 0.04 in favor of the outer leaflet. The results demonstrate the potential of the integrated holographic optical tweezers-Raman technique to induce deformations to individual lipid vesicles and to simultaneously provide quantitative and spatially resolved molecular information. Future studies can extend to include more realistic models of cell membranes and potentially live cells.

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.

The asymmetry of plasma membranes and their cholesterol content influence the uptake of cisplatin

The composition of the plasma membrane of malignant cells is thought to influence the cellular uptake of cisplatin and to take part in developing resistance to this widespread anti-cancer drug. In this work we study the permeation of cisplatin through the model membranes of normal and cancer cells using molecular dynamics simulations. A special attention is paid to lipid asymmetry and cholesterol content of the membranes. The loss of lipid asymmetry, which is common for cancer cells, leads to a decrease in their permeability to cisplatin by one order of magnitude in comparison to the membranes of normal cells. The change in the cholesterol molar ratio from 0% to 33% also decreases the permeability of the membrane by approximately one order of magnitude. The permeability of pure DOPC membrane is 5-6 orders of magnitude higher than one of the membranes with realistic lipid composition, which makes it as an inadequate model for the studies of drug permeability.

Transverse distribution of plasma membrane bilayer cholesterol: Picking sides

The transverse asymmetry (sidedness) of phospholipids in plasma membrane bilayers is well characterized, distinctive, actively maintained and functionally important. In contrast, numerous studies using a variety of techniques have concluded that plasma membrane bilayer cholesterol is either mostly in the outer leaflet or the inner leaflet or is fairly evenly distributed. Sterols might simply partition according to their differing affinities for the asymmetrically disposed phospholipids, but some studies have proposed that it is actively transported to the outer leaflet. Other work suggests that the sterol is enriched in the inner leaflet, driven by either positive interactions with the phosphatidylethanolamine on that side or by its exclusion from the outer leaflet by the long chain sphingomyelin molecules therein. This uncertainty raises three questions: is plasma membrane cholesterol sidedness fixed in a given cell or cell type; is it generally the same among mammalian species; and does it serve specific physiological functions? This review grapples with these issues.

Molecular mechanism of vesicle division induced by coupling between lipid geometry and membrane curvatures

We investigated the effects of lipid geometry on vesicle division using coarse grained molecular dynamics simulations. When the vesicle is composed of zero and negative spontaneous curvature lipids (ZSLs and NSLs), the difference in their molecular spontaneous curvatures destabilizes the neck of the limiting shape vesicle. In the vesicle division pathway, the neck developed into the stalk intermediates. The stalk was broken when the NSLs were expelled from the stalk. Free energy analysis shows that the coupling between the lipid geometry and the Gaussian rigidity is responsible for the observed vesicle division.

The influence of curvature on the properties of the plasma membrane. Insights from atomistic molecular dynamics simulations

In this work we conduct a systematic analysis of the influence of curvature on various properties of a realistic model of mammalian plasma membrane with asymmetric lipid content of monolayers and a realistic concentration of cholesterol. In order to do this we developed new technique for simulating membranes with the global membrane curvature restricted to any desirable value while keeping free lateral diffusion of the lipids and without introducing artifacts or perturbing the membrane structure. We show that the thickness of the hydrophobic core of the concave monolayer decreases by approximately 1.3 Å in comparison to that of the flat membrane, while the thickness of the convex monolayer does not change. The order parameter of the lipid tails decreases significantly in the certain layers of the curved membrane. The area per lipid increases in the convex monolayer and decreases in the concave monolayer. The cholesterol inclination angle decreases when the curvature of a particular monolayer changes from concave to convex. The amount of cholesterol in the minor fraction located between the membrane leaflets is zero in the membrane with positive curvature and increases to 1.7% in the flat membrane and to 2.5% in the membrane with negative curvature.

Molecular dynamics exploration of poration and leaking caused by Kalata B1 in HIV-infected cell membrane compared to host and HIV membranes

The membrane disruption activities of kalata B1 (kB1) were investigated using molecular dynamics simulations with membrane models. The models were constructed to mimic the lipid microdomain formation in membranes of HIV particle, HIV-infected cell, and host cell. The differences in the lipid ratios of these membranes caused the formation of liquid ordered (lo) domains of different sizes, which affected the binding and activity of kB1. Stronger kB1 disruptive activity was observed for the membrane with small sized lo domain. Our results show that kB1 causes membrane leaking without bilayer penetration. The membrane poration mechanism involved in the disorganization of the lo domain and in cholesterol inter-leaflet translocation is described. This study enhances our understanding of the membrane activity of kB1, which may be useful for designing novel and potentially therapeutic peptides based on the kB1 framework.

Interaction of C60 fullerenes with asymmetric and curved lipid membranes: a molecular dynamics study

Interaction of fullerenes with asymmetric and curved DOPC/DOPS bicelles is studied by means of coarse-grained molecular dynamics simulations. The effects caused by asymmetric lipid composition of the membrane leaflets and the curvature of the membrane are analyzed. It is shown that the aggregates of fullerenes prefer to penetrate into the membrane in the regions of the moderately positive mean curvature. Upon penetration into the hydrophobic core of the membrane fullerenes avoid the regions of the extreme positive or the negative curvature. Fullerenes increase the ordering of lipid tails, which are in direct contact with them, but do not influence other lipids significantly. Our data suggest that the effects of the membrane curvature should be taken into account in the studies concerning permeability of the membranes to fullerenes and fullerene-based drug delivery systems.

Asymmetric lipid membranes: towards more realistic model systems

Despite the ubiquity of transbilayer asymmetry in natural cell membranes, the vast majority of existing research has utilized chemically well-defined symmetric liposomes, where the inner and outer bilayer leaflets have the same composition. Here, we review various aspects of asymmetry in nature and in model systems in anticipation for the next phase of model membrane studies.

Cholesterol behavior in asymmetric lipid bilayers: insights from molecular dynamics simulations

Asymmetric lipid composition of the cell membranes plays an important role in the multitude of important biological functions. Much less is known, however, about the distribution and dynamics of cholesterol in asymmetric biological membranes. In this work we show how this issue could be addressed computationally by molecular dynamics simulations. The influence of the lipid head group charge, acyl chain saturation, spontaneous membrane curvature and the surface tension of the membrane on cholesterol distribution in asymmetric lipid bilayers is investigated. Four asymmetric bilayers containing DOPC, DOPS, DSPC, or DSPS lipids, were simulated on the time scale extended to tens of microseconds. We show that cholesterol strongly prefers anionic lipids to neutral and saturated lipid tails to unsaturated with distribution ratio ~0.4-0.6. Multiple flip-flop transitions of cholesterol were observed directly and their mean times range from 350 to 2,000 ns. It was shown that the distribution of cholesterol in the asymmetric bilayer depends not only on the type of lipids but also on the local membrane curvature and the surface tension. The geometric shape of spontaneously curved asymmetric bilayer changes dramatically in the presence of cholesterol. The membrane curvature becomes less homogeneous with large patches of flattened regions interleaved by rather sharp bends.

Cholesterol, sphingolipids, and glycolipids: what do we know about their role in raft-like membranes?

Lipids rafts are considered to be functional nanoscale membrane domains enriched in cholesterol and sphingolipids, characteristic in particular of the external leaflet of cell membranes. Lipids, together with membrane-associated proteins, are therefore considered to form nanoscale units with potential specific functions. Although the understanding of the structure of rafts in living cells is quite limited, the possible functions of rafts are widely discussed in the literature, highlighting their importance in cellular functions. In this review, we discuss the understanding of rafts that has emerged based on recent atomistic and coarse-grained molecular dynamics simulation studies on the key lipid raft components, which include cholesterol, sphingolipids, glycolipids, and the proteins interacting with these classes of lipids. The simulation results are compared to experiments when possible.

Lipid clustering correlates with membrane curvature as revealed by molecular simulations of complex lipid bilayers

Cell membranes are complex multicomponent systems, which are highly heterogeneous in the lipid distribution and composition. To date, most molecular simulations have focussed on relatively simple lipid compositions, helping to inform our understanding of in vitro experimental studies. Here we describe on simulations of complex asymmetric plasma membrane model, which contains seven different lipids species including the glycolipid GM3 in the outer leaflet and the anionic lipid, phosphatidylinositol 4,5-bisphophate (PIP₂), in the inner leaflet. Plasma membrane models consisting of 1500 lipids and resembling the in vivo composition were constructed and simulations were run for 5 µs. In these simulations the most striking feature was the formation of nano-clusters of GM3 within the outer leaflet. In simulations of protein interactions within a plasma membrane model, GM3, PIP₂, and cholesterol all formed favorable interactions with the model α-helical protein. A larger scale simulation of a model plasma membrane containing 6000 lipid molecules revealed correlations between curvature of the bilayer surface and clustering of lipid molecules. In particular, the concave (when viewed from the extracellular side) regions of the bilayer surface were locally enriched in GM3. In summary, these simulations explore the nanoscale dynamics of model bilayers which mimic the in vivo lipid composition of mammalian plasma membranes, revealing emergent nanoscale membrane organization which may be coupled both to fluctuations in local membrane geometry and to interactions with proteins.

Cell membranes play important roles in vivo both in shielding the cell interior from the surrounding environment and in cell function through lipid components of the membrane having roles in controlling protein function, cell signaling etc. We employ molecular dynamics simulations to explore the behavior of biologically realistic membrane models. Our simulations reveal nano-domain clustering of the glycolipid GM3 and to a lesser extent of the anionic lipid phosphatidylinositol 4,5-bisphophate (PIP₂). When including transmembrane proteins we are able to observe preferential interactions of known regulatory lipids (e.g. GM3, PIP₂ and cholesterol) with the proteins. Membrane curvature is shown to be coupled to the local lipid composition, suggestive of a link between lipid nano-domains and membrane geometry.

Influence of different membrane environments on the behavior of cholesterol

Our results show the distribution of cholesterol between stress-free and stressed membranes or between the inner leaflet and the outer leaflet of curved membrane.

Cholesterol induces uneven curvature of asymmetric lipid bilayers

A remarkable flexibility is observed in biological membranes, which allows them to form the structures of different curvatures. We addressed the question of intrinsic ability of phospholipid membranes to form highly curved structures and the role of cholesterol in this process. The distribution of cholesterol in the highly curved asymmetric DOPC/DOPS lipid bilayer was investigated by the coarse-grained molecular dynamics simulations in the membrane patches with large aspect ratio. It is shown that cholesterol induces uneven membrane curvature promoting the formation of extended flattened regions of the membrane interleaved by sharp bends. It is shown that the affinity of cholesterol to anionic DOPS or neutral DOPC lipids is curvature dependent. The cholesterol prefers DOPS to DOPC in either planar or highly curved parts of the membrane. In contrast, in the narrow interval of moderate membrane curvatures this preference is inverted. Our data suggest that there is a complex self-consistent interplay between the membrane curvature and cholesterol distribution in the asymmetric lipid bilayers. The suggested new function of cholesterol may have a biological relevance.