The Phase Behavior and Organization of Sphingomyelin/Cholesterol Membranes: A Deuterium NMR Study

Cholesterol and Lipid Rafts in the Biogenesis of Amyloid-β Protein and Alzheimer's Disease

Cholesterol has been conjectured to be a modulator of the amyloid cascade, the mechanism that produces the amyloid-β (Aβ) peptides implicated in the onset of Alzheimer's disease. We propose that cholesterol impacts the genesis of Aβ not through direct interaction with proteins in the bilayer, but indirectly by inducing the liquid-ordered phase and accompanying liquid-liquid phase separations, which partition proteins in the amyloid cascade to different lipid domains and ultimately to different endocytotic pathways. We explore the full process of Aβ genesis in the context of liquid-ordered phases induced by cholesterol, including protein partitioning into lipid domains, mechanisms of endocytosis experienced by lipid domains and secretases, and pH-controlled activation of amyloid precursor protein secretases in specific endocytotic environments. Outstanding questions on the essential role of cholesterol in the amyloid cascade are identified for future studies.

Biophysical investigation of liposome systems decorated with bioconjugated copolymers in the presence of amantadine

Liposome-based technologies derived from lipids and polymers (e.g., PEGylated liposomes) have been recognized because of their applications in nanomedicine. However, since such systems represent myriad challenges and may promote immune responses, investigation of new biomaterials is mandatory. Here, we report on a biophysical investigation of liposomes decorated with bioconjugated copolymers in the presence (or absence) of amantadine (an antiviral medication). First, copolymers of poly(N,N-dimethylacrylamide-co-fluoresceinacrylate-co-acrylic acid-N-succinimide ester)-block-poly(N-isopropylacrylamide) (PDMA-b-PNIPAM) containing a fluorescence label were biofunctionalized with short peptides that resemble the sequence of the loops 220 and 130 of the binding receptor of the hemagglutinin (HA) protein of the influenza A virus. Then, the bioconjugated copolymers were self-assembled along with liposomes composed of 1,2 dimyristoyl-sn-glycero-3-phosphocholine, sphingomyelin, and cholesterol (MSC). These biohybrid systems, with and without amantadine, were systematically characterized using differential scanning calorimetry (DSC), dynamic light scattering (DLS), and cryogenic transmission electron microscopy (cryoTEM). Finally, the systems were tested in an in vitro study to evaluate cytotoxicity and direct immunofluorescence in Madin Darbin Canine Kidney (MDCK) cells. The biohybrid systems displayed long-term stability, thermo-responsiveness, hydrophilic-hydrophobic features, and fluorescence properties and were presumable endowed with cell targeting properties intrinsically integrated into the amino acid sequences of the utilized peptides, which indeed turn them into promising nanodevices for biomedical applications.

CGCompiler: Automated Coarse-Grained Molecule Parametrization via Noise-Resistant Mixed-Variable Optimization

Coarse-grained force fields (CG FFs) such as the Martini model entail a predefined, fixed set of Lennard-Jones parameters (building blocks) to model virtually all possible nonbonded interactions between chemically relevant molecules. Owing to its universality and transferability, the building-block coarse-grained approach has gained tremendous popularity over the past decade. The parametrization of molecules can be highly complex and often involves the selection and fine-tuning of a large number of parameters (e.g., bead types and bond lengths) to optimally match multiple relevant targets simultaneously. The parametrization of a molecule within the building-block CG approach is a mixed-variable optimization problem: the nonbonded interactions are discrete variables, whereas the bonded interactions are continuous variables. Here, we pioneer the utility of mixed-variable particle swarm optimization in automatically parametrizing molecules within the Martini 3 coarse-grained force field by matching both structural (e.g., RDFs) as well as thermodynamic data (phase-transition temperatures). For the sake of demonstration, we parametrize the linker of the lipid sphingomyelin. The important advantage of our approach is that both bonded and nonbonded interactions are simultaneously optimized while conserving the search efficiency of vector guided particle swarm optimization (PSO) methods over other metaheuristic search methods such as genetic algorithms. In addition, we explore noise-mitigation strategies in matching the phase-transition temperatures of lipid membranes, where nucleation and concomitant hysteresis introduce a dominant noise term within the objective function. We propose that noise-resistant mixed-variable PSO methods can both improve and automate parametrization of molecules within building-block CG FFs, such as Martini.

Eukaryotic Cell Membranes: Structure, Composition, Research Methods and Computational Modelling

This paper deals with the problems encountered in the study of eukaryotic cell membranes. A discussion on the structure and composition of membranes, lateral heterogeneity of membranes, lipid raft formation, and involvement of actin and cytoskeleton networks in the maintenance of membrane structure is included. Modern methods for the study of membranes and their constituent domains are discussed. Various simplified models of biomembranes and lipid rafts are presented. Computer modelling is considered as one of the most important methods. This is stated that from the study of the plasma membrane structure, it is desirable to proceed to the diverse membranes of all organelles of the cell. The qualitative composition and molar content of individual classes of polar lipids, free sterols and proteins in each of these membranes must be considered. A program to create an open access electronic database including results obtained from the membrane modelling of individual cell organelles and the key sites of the membranes, as well as models of individual molecules composing the membranes, has been proposed.

Critical Role of Molecular Packing in Lo Phase Membrane Solubilization

Membrane solubilization induced by Triton X-100 (TX-100) was investigated. Different membrane compositions and phase states were studied along the detergent titration. Expected solubilization profiles were obtained but new information is provided. The fluorescence of nitrobenzoxadiazole (NBD)-labeled lipids indicates that the liquid-ordered (Lo)/liquid-disordered (Ld) phase coexistence is barely unaffected at sub-solubilizing detergent concentrations and highlights the vesicle-to-micelle transition. Moreover, the location of the NBD group in the bilayer emphasizes a detergent-membrane interaction in the case of the insoluble Lo phase membrane. It has also been shown that the molecular packing of the membrane loosens in the presence of TX-100, regardless of the solubilization profile. Motivated by studies on GPMVs, the solubilization of less ordered Lo phase membranes was considered in order to improve the effect of molecular packing on the extent of solubilization. Membranes composed of SM and Chol in an equimolar ratio doped with different amounts of PC were studied. The more ordered the Lo phase membrane is in the absence of detergent, the less likely it is to be solubilized. Furthermore, and in contrast to what is observed for membranes exhibiting an Lo/Ld phase coexistence, a very small decrease in the molecular packing of the Lo phase membrane radically modifies the extent of solubilization. These results have implications for the reliability of TX-100 insolubility as a method to detect ordered domains.

OxPAPC stabilizes liquid-ordered domains in biomimetic membranes

Long-chain polyunsaturated fatty acids (PUFAs) are prone to nonenzymatic oxidation in response to differing environmental stressors and endogenous cellular sources. There is increasing evidence that phospholipids containing oxidized PUFA acyl chains control the inflammatory response. However, the underlying mechanism(s) of action by which oxidized PUFAs exert their functional effects remain unclear. Herein, we tested the hypothesis that replacement of 1-palmitoyl-2-arachidonyl-phosphatidylcholine (PAPC) with oxidized 1-palmitoyl-2-arachidonyl-phosphatidylcholine (oxPAPC) regulates membrane architecture. Specifically, with solid-state 2H NMR of biomimetic membranes, we investigated how substituting oxPAPC for PAPC modulates the molecular organization of liquid-ordered (Lo) domains. 2H NMR spectra for bilayer mixtures of 1,2-dipalmitoylphosphatidylcholine-d62 (an analog of DPPC deuterated throughout sn-1 and -2 chains) and cholesterol to which PAPC or oxPAPC was added revealed that replacing PAPC with oxPAPC disrupted molecular organization, indicating that oxPAPC does not mix favorably in a tightly packed Lo phase. Furthermore, unlike PAPC, adding oxPAPC stabilized 1,2-dipalmitoylphosphatidylcholine-d6-rich/cholesterol-rich Lo domains formed in mixtures with 1,2-dioleoylphosphatidylcholine while decreasing the molecular order within 1,2-dioleoylphosphatidylcholine-rich liquid-disordered regions of the membrane. Collectively, these results suggest a mechanism in which oxPAPC stabilizes Lo domains-by disordering the surrounding liquid-disordered region. Changes in the structure, and thereby functionality, of Lo domains may underly regulation of plasma membrane-based inflammatory signaling by oxPAPC.

Membrane elasticity modulated by cholesterol in model of porcine eye lens-lipid membrane

Experimental evidence shows that the eye lens loses its elasticity dramatically with age. It has also been reported that the cholesterol (Chol) content in the eye lens fiber cell plasma membrane increases significantly with age. High Chol content leads to the formation of cholesterol bilayer domains (CBDs) in the lens membrane. The role of high Chol associated with lens elasticity is unclear. The purpose of this research is to investigate the membrane elasticity of the model of porcine lens-lipid (MPLL) membrane with increasing Chol content to elucidate the role of high Chol in lens membrane elasticity. In this study, we used atomic force microscopy (AFM) to study the mechanical properties (breakthrough force and area compressibility modulus (K_A)) of the MPLL membrane with increasing Chol content where K_A is the measure of membrane elasticity. We varied Chol concentration in Chol/MPLL membrane from 0 to ∼71 mol%. Supported Chol/MPLL membranes were prepared by fusion of small unilamellar vesicles (SUVs) on top of a flat mica surface. SUVs of the Chol/MPLL lipid mixture were prepared with the rapid solvent exchange method followed by probe-tip sonication. For the Chol/MPLL mixing ratio of 0, AFM image showed the formation of two distinct phases of the membrane, i.e., liquid-disordered phase (l_d) and solid-ordered phase (s_o) membrane. However, with Chol/MPLL mixing ratio of 0.5 and above, only liquid-ordered phase (l_o) membrane was formed. Also, two distinct breakthrough forces corresponding to l_d and s_o were observed for Chol/MPLL mixing ratio of 0, whereas only one breakthrough force was observed for membranes with Chol/MPLL mixing ratio of 0.5 and above. No significant difference in the membrane surface roughness was measured with increasing Chol content for these membranes; however, breakthrough force and K_A for l_o membrane increased when Chol/MPLL mixing ratio was increased from 0.5 to 1. Interestingly above the Chol/MPLL mixing ratio of 1, both breakthrough force and K_A decreased, indicating the formation of CBDs. Furthermore, these results showed that membrane elasticity increases at high Chol content, suggesting that high Chol content in lens membrane might be responsible for maintaining lens membrane elasticity.

Raman scattering in protonated and deuterated 1,2-dipalmitoyl-sn-glycero-3-phosphatidylcholine (DPPC): Indicators of conformational and lateral orders

The use of deuterocarbons is an effective method in the Raman spectroscopy of multicomponent lipid materials and biological samples. Here, Raman spectra of hydrated multilamellar vesicles of 1,2-dipalmitoyl-sn-glycero-3-phosphocholine (DPPC), its deuterated analog 1,2-dipalmitoyl-d62-sn-glycero-3-phosphocholine (DPPC_d62), and DPPC-DPPC_d62 mixtures were studied in a wide temperature range to specify the Raman indicators of conformational and lateral orders. The temperature dependence of the 985 cm^-1 line in the deuterated phospholipid unequivocally indicates that this line corresponds to the CC stretching vibrations of deuterated hydrocarbon chains in the all-trans conformation. It was also concluded that the ratio of Raman intensities at the maximum of the peak of the symmetric CD₂ stretching and at a maximum near 2168 cm^-1 reflects the conformational order of the hydrocarbon chain and can be used for an evaluation of the fraction of the all-trans sequences. The frequency of the symmetric CD₂ stretching peak is sensitive to the phase state (gel or fluid) but has a low sensitivity to the partial conformational disordering within the gel phase. The Raman study of DPPC-DPPC_d62 mixtures reveals that the lateral order contributes to the ratio of intensities of the antisymmetric and symmetric CH₂ stretching peaks as a prefactor enhancing the effect of conformational ordering.

Molecular Organization of a Raft-like Domain in a Polyunsaturated Phospholipid Bilayer: A Supervised Machine Learning Analysis of Molecular Dynamics Simulations

Numerous health benefits are associated with omega-3 polyunsaturated fatty acids (n-3 PUFA) consumed in fish oils. An understanding of the mechanism remains elusive. The plasma membrane as a site of action is the focus in this study. With large-scale all-atom MD simulations run on a model membrane (1050 lipid molecules), we observed the evolution over time (6 μs) of a circular (raft-like) domain composed of N-palmitoylsphingomyelin (PSM) and cholesterol embedded into a surrounding (non-raft) patch composed of polyunsaturated 1-palmitoyl-2-docosahexaenoylphosphatylcholine (PDPC) (1:1:1 mol). A supervised machine learning algorithm was developed to characterize the migration of each lipid based on molecular conformation and the local environment. PDPC molecules were seen to infiltrate the ordered raft-like domain in a small amount, while a small concentration of PSM and cholesterol molecules was seen to migrate into the disordered non-raft region. Enclosing the raft-like domain, a narrow (∼2 nm in width) interfacial zone composed of PDPC, PSM, and cholesterol that buffers the substantial difference in order (ΔS_CD ≈ 0.12) between raft-like and non-raft environments was seen to form. Our results suggest that n-3 PUFA regulate the architecture of lipid rafts enriched in sphingolipids and cholesterol with a minimal effect on order within their interior in membranes.

CHO/LY-B cell growth under limiting sphingolipid supply: Correlation between lipid composition and biophysical properties of sphingolipid-restricted cell membranes

Sphingolipids (SL) are ubiquitous in mammalian cell membranes, yet there is little data on the behavior of cells under SL-restriction conditions. LY-B cells derive from a CHO linein whichserine palmitoyl transferase (SPT), thus de novo SL synthesis, is suppressed, while maintaining the capacity of taking up and metabolizing exogenous sphingoid bases from the culture medium. In this study, LY-B cells were adapted to grow in a fetal bovine serum (FBS)-deficient medium to avoid external uptake of lipids. The lowest FBS concentration that allowed LY-B cell growth, though at a slow rate, under our conditions was 0.04%, that is, 250-fold less than the standard (10%) concentration. Cells grown under limiting SL concentrations remained viable for at least 72 hours. Enriching with sphingomyelin the SL-deficient medium allowed the recovery of growth rates analogous to those of control LY-B cells. Studies including whole cells, plasma membrane preparations, and derived lipid vesicles were carried out. Laurdan fluorescence was recorded to measure membrane molecular order, showing a significant decrease in the rigidity of LY-B cells, not only in plasma membrane but also in whole cell lipid extract, as a result of SL limitation in the growth medium. Plasma membrane preparations and whole cell lipid extracts were also studied using atomic force microscopy in the force spectroscopy mode. Force measurements demonstrated that lower breakthrough forces were required to penetrate samples obtained from SL-poor LY-B cells than those obtained from control cells. Mass-spectroscopic analysis was also a helpful tool to understand the rearrangement undergone by the LY-B cell lipid metabolism. The most abundant SL in LY-B cells, sphingomyelin, decreased by about 85% as a result of SL limitation in the medium, the bioactive lipid ceramide and the ganglioside precursor hexosylceramide decreased similarly, together with cholesterol. Quantitative SL analysis showed that a 250-fold reduction in sphingolipid supply to LY-B cells led only to a sixfold decrease in membrane sphingolipids, underlining the resistance to changes in composition of these cells. Plasma membrane compositions exhibited similar changes, at least qualitatively, as the whole cells with SL restriction. A linear correlation was observed between the sphingomyelin concentration in the membranes, the degree of lipid order as measured by laurdan fluorescence, and membrane breakthrough forces assessed by atomic force microscopy. Smaller, though significant, changes were also detected in glycerophospholipids under SL-restriction conditions.

Melatonin Alters Fluid Phase Coexistence in POPC/DPPC/Cholesterol Membranes

The structure and biophysical properties of lipid membranes are important for cellular functions in health and disease. In Alzheimer's disease, the neuronal membrane is a target for toxic amyloid-β (Aβ). Melatonin is an important pineal gland hormone that has been shown to protect against Aβ toxicity in cellular and animal studies, but the molecular mechanism of this protection is not fully understood. Melatonin is a small membrane-active molecule that has been shown to interact with model lipid membranes and alter the membrane biophysical properties, such as membrane molecular order and dynamics. This effect of melatonin has been previously studied in simple model bilayers with one or two lipid components. To make it more relevant to neuronal membranes, we used a more complex ternary lipid mixture as our membrane model. In this study, we used ²H-NMR to investigate the effect of melatonin on the phase behavior of 1,2-dipalmitoyl-sn-glycero-3-phosphocholine (DPPC), 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphocholine (POPC), and cholesterol lipid membranes. We used deuterium-labeled POPC-d₃₁ and DPPC-d₆₂,separately to probe the changes in hydrocarbon chain order as a function of temperature and melatonin concentration. We find that POPC/DPPC/cholesterol at molar proportions of 3:3:2 is close to liquid-disordered/liquid-ordered phase separation and that melatonin can induce phase separation in these ternary mixtures by preferentially incorporating into the disordered phase and increasing its level of disorder. At 5 mol% melatonin, we observed phase separation in samples with POPC-d₃₁, but not with DPPC-d₆₂, whereas at 10 mol% melatonin, phase separation was observed in both samples with either POPC-d₃₁ or DPPC-d₆₂. These results indicate that melatonin can have a strong effect on membrane structure and physical properties, which may provide some clues to understanding how melatonin protects against Aβ, and that choice of chain perdeuteration is an important consideration from a technical point of view.

Structure of supported DPPC/cholesterol bilayers studied via X-ray reflectivity

The electron density profile of bilayers of DPPC/cholesterol mixtures supported on semiconductor grade silicon substrates were studied with the objective of determining how the proximity of a solid interface modifies the phase diagram of mixed bilayers. The bilayers were studied in situ immersed in water via synchrotron X-ray reflectivity (XRR). Measurements were performed as a function of temperature through the main phase transition and cholesterol mole fractions up to 40%. Analysis of XRR yields the bilayer thickness, roughness and leaflet asymmetry. We find that the structure of the pure DPPC bilayers in the gel phase is in agreement with previous X-ray measurements of supported bilayers deposited via vesicle fusion and multilamellar vesicles but show more clearly defined features than measurements made on films formed using Langmuir-Blodget Langmuir-Shaffer (LB) deposition. Examination of bilayer thickness vs. temperature shows that the melting temperature for supported bilayers is shifted upwards by approximately 4 °C relative to multilamellar vesicles and that the melting temperature decreases with increasing cholesterol content up to 20%. For pure DPPC bilayers the leaflets melt in two stages with the distal leaflet melting first. For cholesterol concentrations of 10% and 20% there is no clear indication of separate melting. For 33% and 40% cholesterol content no clear transition is seen in the bilayer thickness, but an abrupt change in roughness indicates possible microdomain formation in the 40% cholesterol sample.

Improved Characterization of Raft-Mimicking Phase-Separation Phenomena in Lipid Bilayers Using Laurdan Fluorescence with Log-Normal Multipeak Analysis

The study of biomimetic model membrane systems undergoing liquid-ordered (Lo)-liquid-disordered (Ld) phase separation using spectroscopic methods has played an important role in understanding the properties of lipid rafts in plasma membranes. In particular, the membrane-associated fluorescence probe Laurdan has proved to be a very efficient reporter of Lo-Ld phase separation in lipid bilayers using the general polarization (GP) parameter. A limitation of the GP approach is that it monitors only global average packing so that the contribution of each phase remains undetermined. The decomposition of Laurdan emission spectra has been proposed as an additional approach to overcoming this limitation. Here, further developments of this method for the study of Lo-Ld phase separation are described here for Laurdan in sphingomyelin-phosphatidylcholine-cholesterol large unilamellar vesicles. Lipid compositions corresponding to homogeneous Lo or Ld phases as well as undergoing thermally induced Lo-Ld phase separation were investigated. In addition, the occurrence of phase separation was checked by the fluorescence imaging of giant unilamellar vesicles. Decomposition into three log-normal components is used to show that an intermediate energy component is specifically associated with the occurrence of the Lo phase, with a small contribution from this component occurring above the phase-separation temperature being attributable to phase fluctuations. The ratio R_X of the relative area of this intermediate-energy peak to that of the low-energy peak is shown to provide a straightforward index of Lo-Ld phase separation as a function of temperature, which is occasionally more sensitive than GP. It is also shown that R_X can be used in conjunction with GP to gain further insight into Lo-Ld, the phase-separation processes. This latter feature is illustrated by the influence of the alcohol butanol on the Lo-Ld phase separation in sphingomyelin-phosphatidylcholine-cholesterol bilayers by showing that the effect of the alcohol occurs specifically at the onset of the phase separation, indicating a line tension mechanism. It is proposed that the three components of log-normal decomposition approaching Laurdan emission spectra provide a useful improvement for characterizing Lo-Ld phase-separation phenomena.

Near-field studies of anisotropic variations and temperature-induced structural changes in a supported single lipid bilayer

Temperature-controlled polarization modulation near-field scanning optical microscopy measurements of a single supported L_{β^{'}} 1,2-dipalmitoyl-sn-glycero-3-phosphocholine (DPPC) lipid bilayer are presented. The effective retardance (S=2π(n_{e}-n_{o})t/λ, where t is the thickness of the bilayer and λ is the wavelength of light used) and the direction of the projection of the acyl chains (θ) were measured simultaneously. We demonstrate how one is able to align the system over the sample and measure a relative retardance ΔS, a crucial step in performing temperature-controlled experiments. Maps of ΔS and θ, with a lateral resolution on the order of ∼100 nm are presented, highlighting variations deriving from changes in the average molecular orientation across a lipid membrane at room temperature. A discussion of how this information can be used to map the average three-dimensional orientation of the molecules is presented. From ΔS and the known thickness of the membrane t the birefringence (n_{e}-n_{o}) of the bilayer was determined. Temperature-controlled measurements yielded a change of ΔS∼(3.8±0.3) mrad at the main transition temperature (T_{m}∼41^{∘}C) of a single planar bilayer of DPPC, where the membrane transitioned between the gel L_{β^{'}} to liquid disorder L_{α} state. The result agrees well with previous values of (n_{e}-n_{o}) in the L_{β^{'}} phase and translates to an assumed average acyl chain orientation relative to the membrane normal (〈ϕ〉∼32^{∘}) when T<T_{m} and 0^{∘} when T>T_{m}. Evidence of super heating and cooling are presented. A discussion on how the observed behavior as T_{m} is approached, could relate to the existence of varying microconfigurations within the lipid bilyer are presented. This conversation includes ideas from a Landau-Ginzburg picture of first-order phase transitions in nematic-to-isotropic systems.

Ordering of Saturated and Unsaturated Lipid Membranes near Their Phase Transitions Induced by an Amphiphilic Cyclodextrin and Cholesterol

When inserted in membranes of dimyristoyl phosphatidylcholine (DMPC), methylated β-cyclodextrins with one (TrimβMLC) or two (TrimβDLC) lauryl acyl chains grafted onto the hydrophilic cavity exert a "cholesterol-like ordering effect", by straightening the acyl chains in the fluid phase at temperatures near the chain melting transition. This effect may be related to pretransitional events such as the "anomalous swelling" known to occur with saturated phosphatidylcholine membranes. To investigate this model, order profiles and bilayer thicknesses of DMPC and unsaturated 1-palmitoyl-2-oleoyl-phosphatidylcholine (POPC) membranes containing amphiphilic cyclodextrins or cholesterol were determined by deuterium NMR. The pure lipid membranes display both a qualitatively similar chain ordering upon cooling in the fluid phase, more important at the chain extremity, which gets more pronounced near their fluid-to-gel transitions. Both membranes show a bilayer thickness increase by ∼0.5 Å just above their transition, as observed previously with saturated phosphatidylcholines of various chain lengths. Membrane-insertion of 5% TrimβMLC or cholesterol induces an important ordering of the DMPC acyl chains just above the transition, which is also more pronounced at the chain extremity. There is an additional increase of the bilayer thickness, most probably due to a deep insertion of these amphiphilic molecules, facilitated by increased bilayer softness in the anomalous swelling regime. These effects are more important with TrimβMLC than with cholesterol. By contrast, no enhanced acyl chain ordering was observed when approaching the transition of TrimβMLC-containing POPC membranes, as a possible consequence of an eventual lack of anomalous swelling in unsaturated lipid membranes. Insertion of higher concentrations of TrimβMLC was found to induce a magnetic orientation of the DMPC membranes in the fluid phase with 10% of this derivative, coupled with the appearance of a broad isotropic component when the concentration is raised to 20%. No membrane orientation or isotropic component was detected with TrimβMLC-containing POPC membranes.

Impact of Acyl Chain Mismatch on the Formation and Properties of Sphingomyelin-Cholesterol Domains

Lateral segregation and the formation of lateral domains are well-known phenomena in ternary lipid bilayers composed of an unsaturated (low gel-to-liquid phase transition temperature (T_m)) phospholipid, a saturated (high-T_m) phospholipid, and cholesterol. The formation of lateral domains has been shown to be influenced by differences in phospholipid acyl chain unsaturation and length. Recently, we also showed that differential interactions of cholesterol with low- and high-T_m phospholipids in the bilayer can facilitate phospholipid segregation. Now, we have investigated phospholipid-cholesterol interactions and their role in lateral segregation in ternary bilayers composed of different unsaturated phosphatidylcholines (PCs) with varying acyl chain lengths, N-palmitoyl-D-erythro-sphingomyelin (PSM), and cholesterol. Using deuterium NMR spectroscopy, we determined how PSM was influenced by the acyl chain composition in surrounding PC environments and correlated this with the affinity of cholestatrienol (a fluorescent cholesterol analog) for PSM in the different PC environments. Results from a combination of time-resolved fluorescence measurements of trans-parinaric acid and Förster resonance energy transfer experiments showed that the relative affinity of cholesterol for phospholipids determined the degree to which the sterol promoted domain formation. From Förster resonance energy transfer, deuterium NMR, and differential scanning calorimetry results, it was clear that cholesterol also influenced both the thermostability of the domains and the degree of order in and outside the PSM-rich domains. The results of this study have shown that the affinity of cholesterol for both low-T_m and high-T_m phospholipids and the effects of low- and high-T_m phospholipids on each other influence both lateral structure and domain properties in complex bilayers. We envision that similar effects also contribute to lateral heterogeneity in even more complex biological membranes.

Lamellar Phases Composed of Phospholipid, Cholesterol, and Ceramide, as Studied by 2H NMR

Sphingolipids constitute a significant fraction of cellular plasma membrane lipid content. Among sphingolipids, ceramide levels are usually very low. However, in some cell processes like apoptosis, cell membrane ceramide levels increase markedly because of the activation of enzymes like acid sphingomyelinase. This increase can change the physical state of the membrane by promoting molecular order and inducing solid-ordered (So) phase domains. This effect has been observed in a previous 2H NMR study on membranes consisting of palmitoyl sphingomyelin (PSM) and palmitoyl ceramide (PCer). Cholesterol (Chol), too, is present at high concentrations in mammalian plasma membranes and has a favorable interaction with sphingomyelin (SM), together forming domains in the liquid-ordered phase in model membranes. There are reports that Chol is able to displace ceramide (Cer) in SM bilayers and abolish the So phase domains formed by SM:Cer. This ability of Chol appears to be concentration dependent; in membranes with low Chol and high Cer contents, So phase domains rich in Cer coexist with the continuous fluid phase of the membrane. Here, we studied the effect of increasing PCer concentration in PSM:Chol bilayers, using 2H NMR. Chol:PCer mole ratios were 3:1, 3:2, and 3:3, at a fixed 7:3 phospholipid:cholesterol mol ratio. Both PSM and PCer were monitored in separate samples for changes in their physical state by introducing a perdeuterated palmitoyl chain in either molecule. Moreover, the effect of replacing PSM with DPPC was investigated to test the impact on membrane phase behavior of replacing the sphingosine with a palmitoylated glycerol backbone. We found that PCer can increase acyl chain order in both PSM:Chol and DPPC:Chol bilayers. Especially in bilayers with Chol:PCer 1:1 molar ratios, PCer induces highly stable So phase domains in both PSM and DPPC bilayers near 37°C. However, PCer has a more pronounced ordering effect on PSM compared to DPPC bilayers.

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.

How to Determine Lipid Interactions in Membranes from Experiment Through the Ising Model

The determination and the meaning of interactions in lipid bilayers are discussed and interpreted through the Ising model. Originally developed to understand phase transitions in ferromagnetic systems, the Ising model applies equally well to lipid bilayers. In the case of a membrane, the essence of the Ising model is that each lipid is represented by a site on a lattice and that the interaction of each site with its nearest neighbors is represented by an energy parameter ω. To calculate the thermodynamic properties of the system, such as the enthalpy, the Gibbs energy, and the heat capacity, the partition function is derived. The calculation of the configurational entropy factor in the partition function, however, requires approximations or the use of Monte Carlo (MC) simulations. Those approximations are described. Ultimately, MC simulations are used in combination with experiment to determine the interaction parameters ω in lipid bilayers. Several experimental approaches are described, which can be used to obtain interaction parameters. They include nearest-neighbor recognition, differential scanning calorimetry, and Förster resonance energy transfer. Those approaches are most powerful when used in combination of MC simulations of Ising models. Lipid membranes of different compositions are discussed, which have been studied with these approaches. They include mixtures of cholesterol, saturated (ordered) phospholipids, and unsaturated (disordered) phospholipids. The interactions between those lipid species are examined as a function of molecular properties such as the degree of unsaturation and the acyl chain length. The general rule that emerges is that interactions between different lipids are usually unfavorable. The exception is that cholesterol interacts favorably with saturated (ordered) phospholipids. However, the interaction of cholesterol with unsaturated phospholipids becomes extremely unfavorable as the degree of unsaturation increases.

Cholesterol Effects on the Physical Properties of Lipid Membranes Viewed by Solid-state NMR Spectroscopy

In this chapter, we review the physical properties of lipid/cholesterol mixtures involving studies of model membranes using solid-state NMR spectroscopy. The approach allows one to quantify the average membrane structure, fluctuations, and elastic deformation upon cholesterol interaction. Emphasis is placed on understanding the membrane structural deformation and emergent fluctuations at an atomistic level. Lineshape measurements using solid-state NMR spectroscopy give equilibrium structural properties, while relaxation time measurements study the molecular dynamics over a wide timescale range. The equilibrium properties of glycerophospholipids, sphingolipids, and their binary and tertiary mixtures with cholesterol are accessible. Nonideal mixing of cholesterol with other lipids explains the occurrence of liquid-ordered domains. The entropic loss upon addition of cholesterol to sphingolipids is less than for glycerophospholipids, and may drive formation of lipid rafts. The functional dependence of ²H NMR spin-lattice relaxation (R _1Z) rates on segmental order parameters (S _CD) for lipid membranes is indicative of emergent viscoelastic properties. Addition of cholesterol shows stiffening of the bilayer relative to the pure lipids and this effect is diminished for lanosterol. Opposite influences of cholesterol and detergents on collective dynamics and elasticity at an atomistic scale can potentially affect lipid raft formation in cellular membranes.

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.

Heat Capacity of DPPC/Cholesterol Mixtures: Comparison of Single Bilayers with Multibilayers and Simulations

The excess heat capacity (Δ C _p) of mixtures of dipalmitoylphosphatidylcholine (DPPC) and cholesterol (Chol) is examined in detail in large unilamellar vesicles (LUVs), both experimentally, using differential scanning calorimetry (DSC), and theoretically, using a three-state Ising model. The model postulates that DPPC can access three conformational states: gel, liquid-disordered (L_d), and liquid-ordered (L_o). The L_o state, however, is only available if coupled with interaction with an adjacent Chol. Δ C _p was calculated using Monte Carlo simulations on a lattice and compared to experiment. The DSC results in LUVs are compared with literature data on multilamellar vesicles (MLVs). The enthalpy change of the complete phase transition from gel to L_d is identical in LUVs and MLVs, and the melting temperatures ( T_m) are similar. However, the DSC curves in LUVs are significantly broader, and the maxima of Δ C _p are accordingly smaller. The parameters in the Ising model were chosen to match the DSC curves in LUVs and the nearest-neighbor recognition (NNR) data. The model reproduces the NNR data very well. It also reproduces the phase transition in DPPC, the freezing point depression induced by Chol, and the broad component of Δ C _p in DPPC/Chol LUVs. However, there is a sharp component, between 5 and 15 mol % Chol, that the model does not reproduce. The broad component of Δ C _p becomes dominant as Chol concentration increases, indicating that it involves melting of the L_o phase. Because the simulations reproduce this component, the conclusions regarding the nature of the phase transition at high Chol concentrations and the structure of the L_o phase are important: there is no true phase separation in DPPC/Chol LUVs. There are large domains of gel and L_o phase coexisting below T_m of DPPC, but above T_m the three states of DPPC are mixed with Chol, although clusters persist.