Mixing properties of sphingomyelin ceramide bilayers: a simulation study

Interactions between the Nicotinic and Endocannabinoid Receptors at the Plasma Membrane

Compartmentalization, together with transbilayer and lateral asymmetries, provide the structural foundation for functional specializations at the cell surface, including the active role of the lipid microenvironment in the modulation of membrane-bound proteins. The chemical synapse, the site where neurotransmitter-coded signals are decoded by neurotransmitter receptors, adds another layer of complexity to the plasma membrane architectural intricacy, mainly due to the need to accommodate a sizeable number of molecules in a minute subcellular compartment with dimensions barely reaching the micrometer. In this review, we discuss how nature has developed suitable adjustments to accommodate different types of membrane-bound receptors and scaffolding proteins via membrane microdomains, and how this "effort-sharing" mechanism has evolved to optimize crosstalk, separation, or coupling, where/when appropriate. We focus on a fast ligand-gated neurotransmitter receptor, the nicotinic acetylcholine receptor, and a second-messenger G-protein coupled receptor, the cannabinoid receptor, as a paradigmatic example.

The role of the 'sphingoid motif' in shaping the molecular interactions of sphingolipids in biomembranes

Sphingolipids can be differentiated from other membrane lipids by the distinctive chemistry of the sphingoid long chain base (LCB), which is generated by the condensation of an amino acid (normally but not always serine) and a fatty acyl CoA (normally palmitoyl CoA) by the pyridoxal phosphate-dependent enzyme, serine palmitoyl transferase (SPT). The first five carbon atoms of the sphingoid LCB, herein defined as the 'sphingoid motif', are largely responsible for the unique chemical and biophysical properties of sphingolipids since they can undergo a relatively large number (compared to other lipid species) of molecular interactions with other membrane lipids, via hydrogen-bonding, charge-pairing, hydrophobic and van der Waals interactions. These interactions are responsible, for instance, for the association of sphingolipids with cholesterol in the membrane lipid bilayer. Here, we discuss some of the unique properties of this sphingoid motif, and in addition to outlining how this structural motif drives intra-bilayer interactions, discuss the atomic details of the interactions with two critical players in the biosynthetic pathway, namely SPT, and the ceramide transport protein, CERT. In the former, the selectivity of sphingolipid synthesis relies on a hydrogen bond interaction between Lys379 of SPTLC2 and the l-serine sidechain hydroxyl moiety. In the latter, the entire sphingoid motif is stereoselectively recognized by a hydrogen-bonding network involving all three sphingoid motif heteroatoms. The remarkable selectivity of these interactions, and the subtle means by which these interactions are modified and regulated in eukaryotic cells raises a number of challenging questions about the generation of these proteins, and of their interactions with the sphingoid motif in evolutionary history.

Update of the CHARMM36 United Atom Chain Model for Hydrocarbons and Phospholipids

Accurate lipid force field (FF) parameters used in molecular dynamics (MD) simulations are crucial for understanding the properties of lipid-containing systems and biological processes related to lipids. The last update of the CHARMM36 united atom chain model (C36UA) was in 2013 [Lee, S. J. Phys. Chem. B 2014, 118, 547 556]; it utilized CHARMM36 (C36) lipid FF parameters for headgroups and OPLS-UA Lennard-Jones (LJ) parameters for tails. Simulations with the FF were able to reproduce many experimental observables of lipid bilayers accurately, but to be more applicable for a wide range of lipids, additional FF parameter optimization was needed. In this work, we present an update of the model, named C36UAr. The parameterization included the LJ parameters for hydrocarbons and related dihedrals. Bulk liquid properties (density, heat of vaporization, isothermal compressibility, and diffusion constant) of model compounds were used as targets for the LJ parameter fitting, and dihedrals were fit to either quantum mechanical (QM) or potential of mean force (PMF) calculations using C36. Thermodynamic reweighting was used to further improve the parameters. Bilayer simulations of various lipid headgroups (phosphatidylcholine, phosphatidylethanolamine, and phosphatidylglycerol) and tails (saturated, monounsaturated, and polyunsaturated) were performed to validate the model, and significant improvements were seen in bilayer properties, including surface area, membrane thicknesses, NMR deuterium order parameters, and density profiles. C36UAr was also compared to the hydrogen mass repartitioning (HMR) method. The high accuracy and competitive efficiency shown in this study make C36UAr one of the best choices for studies of membrane structure and membrane-associated proteins.

Molecular Structure of Sphingomyelin in Fluid Phase Bilayers Determined by the Joint Analysis of Small-Angle Neutron and X-ray Scattering Data

We have determined the fluid bilayer structure of palmitoyl sphingomyelin (PSM) and stearoyl sphingomyelin (SSM) by simultaneously analyzing small-angle neutron and X-ray scattering data. Using a newly developed scattering density profile (SDP) model for sphingomyelin lipids, we report structural parameters including the area per lipid, total bilayer thickness, and hydrocarbon thickness, in addition to lipid volumes determined by densitometry. Unconstrained all-atom simulations of PSM bilayers at 55 °C using the C36 CHARMM force field produced a lipid area of 56 Ų, a value that is 10% lower than the one determined experimentally by SDP analysis (61.9 Ų). Furthermore, scattering form factors calculated from the unconstrained simulations were in poor agreement with experimental form factors, even though segmental order parameter (S_CD) profiles calculated from the simulations were in relatively good agreement with S_CD profiles obtained from NMR experiments. Conversely, constrained area simulations at 61.9 Ų resulted in good agreement between the simulation and experimental scattering form factors, but not with S_CD profiles from NMR. We discuss possible reasons for the discrepancies between these two types of data that are frequently used as validation metrics for molecular dynamics force fields.

Theoretical investigations into the intermolecular hydrogen-bonding interactions of N-(hydroxymethyl)acetamide dimers

The structures of the N-(hydroxymethyl)acetamide (model molecule of ceramide) dimers have been fully optimized at B3LYP/6-311++G** level. The intermolecular hydrogen bonding interaction energies have been calculated using the B3LYP/6-311++G**, B3LYP/6-311++G(2df,2p), MP2(full)/6-311++G** and MP2(full)/6-311++G(2df,2p) methods, respectively. The results show that the O-H···O, N-H···O, O-H···N, and C-H···O hydrogen bonding interactions could exist in N-(hydroxymethyl)acetamide dimers, and the O-H···O, N-H···O, and O-H···N hydrogen bonding interactions could be stronger than C-H···O. The three-dimensional network structure formed by ceramide molecules through intermolecular hydrogen bonding interactions may be the main reason why the stratum corneum of skin could prevent foreign substances from entering our body, as is in accordance with the experimental results. The stability of hydrogen-bonding interactions follow the order of (a) > (b) ≈ (c) > (d) > (e) ≈ (f) > (g) > (h). The analyses of the energy decomposition, frequency, atoms in molecules (AIM), natural bond orbital (NBO), and electron density shift are used to further reveal the nature of the complex formation. In the range of 263.0-328.0 K, the complex is formed via an exothermic reaction, and the solvent with lower temperature and dielectric constant is favorable to this process. Graphical abstract The structures and the O-H···O=C, N-H···O=C and C-H···O=C H-bonding interactions in the N-(hydroxymethyl)acetamide (model molecule of ceramide) dimers were investigated using the B3LYP and MP2(full) methods.

Anandamide Revisited: How Cholesterol and Ceramides Control Receptor-Dependent and Receptor-Independent Signal Transmission Pathways of a Lipid Neurotransmitter

Anandamide is a lipid neurotransmitter derived from arachidonic acid, a polyunsaturated fatty acid. The chemical differences between anandamide and arachidonic acid result in a slightly enhanced solubility in water and absence of an ionisable group for the neurotransmitter compared with the fatty acid. In this review, we first analyze the conformational flexibility of anandamide in aqueous and membrane phases. We next study the interaction of the neurotransmitter with membrane lipids and discuss the molecular basis of the unexpected selectivity of anandamide for cholesterol and ceramide from among other membrane lipids. We show that cholesterol behaves as a binding partner for anandamide, and that following an initial interaction mediated by the establishment of a hydrogen bond, anandamide is attracted towards the membrane interior, where it forms a molecular complex with cholesterol after a functional conformation adaptation to the apolar membrane milieu. The complex is then directed to the anandamide cannabinoid receptor (CB1) which displays a high affinity binding pocket for anandamide. We propose that cholesterol may regulate the entry and exit of anandamide in and out of CB1 by interacting with low affinity cholesterol recognition sites (CARC and CRAC) located in transmembrane helices. The mirror topology of cholesterol binding sites in the seventh transmembrane domain is consistent with the delivery, extraction and flip-flop of anandamide through a coordinated cholesterol-dependent mechanism. The binding of anandamide to ceramide illustrates another key function of membrane lipids which may occur independently of protein receptors. Interestingly, ceramide forms a tight complex with anandamide which blocks the degradation pathway of both lipids and could be exploited for anti-cancer therapies.

Phase behavior of palmitoyl and egg sphingomyelin

Despite the biological significance of sphingomyelins (SMs), there is far less structural information available for SMs compared to glycerophospholipids. Considerable confusion exists in the literature regarding even the phase behavior of SM bilayers. This work studies both palmitoyl (PSM) and egg sphingomyelin (ESM) in the temperature regime from 3 °C to 55 °C using X-ray diffraction and X-ray diffuse scattering on hydrated, oriented thick bilayer stacks. We observe clear evidence for a ripple phase for ESM in a large temperature range from 3 °C to the main phase transition temperature (T_M) of ∼38 °C. This unusual stability of the ripple phase was not observed for PSM, which was in a gel phase at 3 °C, with a gel-to-ripple transition at ∼24 °C and a ripple-to-fluid transition at ∼41 °C. We also report structural results for all phases. In the gel phase at 3 °C, PSM has chains tilted by ∼30° with an area/lipid ∼45 Ų as determined by wide angle X-ray scattering. The ripple phases for both PSM and ESM have temperature dependent ripple wavelengths that are ∼145 Å near 30 °C. In the fluid phase, our electron density profiles combined with volume measurements allow calculation of area/lipid to be ∼64 Ų for both PSM and ESM, which is larger than that from most of the previous molecular dynamics simulations and experimental studies. Our study demonstrates that oriented lipid films are particularly well-suited to characterize ripple phases since the scattering pattern is much better resolved than in unoriented samples.

Analytical Techniques in Lipidomics: State of the Art

Current studies related to lipid identification and determination, or lipidomics in biological samples, are one of the most important issues in modern bioanalytical chemistry. There are many articles dedicated to specific analytical strategies used in lipidomics in various kinds of biological samples. However, in such literature, there is a lack of articles dedicated to a comprehensive review of the actual analytical methodologies used in lipidomics. The aim of this article is to characterize the lipidomics methods used in modern bioanalysis according to the methodological point of view: (1) chromatography/separation methods, (2) spectroscopic methods and (3) mass spectrometry and also hyphenated methods. In the first part, we discussed thin layer chromatography (TLC), high-pressure liquid chromatography (HPLC), gas chromatography (GC) and capillary electrophoresis (CE). The second part includes spectroscopic techniques such as Raman spectroscopy (RS), Fourier transform infrared spectroscopy (FT-IR) and nuclear magnetic resonance (NMR). The third part is a synthetic review of mass spectrometry, matrix-assisted laser desorption/ionization (MALDI), hyphenated methods, which include liquid chromatography-mass spectrometry (LC-MS), gas chromatography-mass spectrometry (GC-MS) and also multidimensional techniques. Other aspects are the possibilities of the application of the described methods in lipidomics studies. Due to the fact that the exploration of new methods of lipidomics analysis and their applications in clinical and medical studies are still challenging for researchers working in life science, we hope that this review article will be very useful for readers.

A Coarse-Grained Model of Stratum Corneum Lipids: Free Fatty Acids and Ceramide NS

Ceramide (CER)-based biological membranes are used both experimentally and in simulations as simplified model systems of the skin barrier. Molecular dynamics studies have generally focused on simulating preassembled structures using atomistically detailed models of CERs, which limit the system sizes and time scales that can practically be probed, rendering them ineffective for studying particular phenomena, including self-assembly into bilayer and lamellar superstructures. Here, we report on the development of a coarse-grained (CG) model for CER NS, the most abundant CER in human stratum corneum. Multistate iterative Boltzmann inversion is used to derive the intermolecular pair potentials, resulting in a force field that is applicable over a range of state points and suitable for studying ceramide self-assembly. The chosen CG mapping, which includes explicit interaction sites for hydroxyl groups, captures the directional nature of hydrogen bonding and allows for accurate predictions of several key structural properties of CER NS bilayers. Simulated wetting experiments allow the hydrophobicity of CG beads to be accurately tuned to match atomistic wetting behavior, which affects the whole system, since inaccurate hydrophobic character is found to unphysically alter the lipid packing in hydrated lamellar states. We find that CER NS can self-assemble into multilamellar structures, enabling the study of lipid systems more representative of the multilamellar lipid structures present in the skin barrier. The coarse-grained force field derived herein represents an important step in using molecular dynamics to study the human skin barrier, which gives a resolution not available through experiment alone.

CHARMM all-atom additive force field for sphingomyelin: elucidation of hydrogen bonding and of positive curvature

The C36 CHARMM lipid force field has been extended to include sphingolipids, via a combination of high-level quantum mechanical calculations on small molecule fragments, and validation by extensive molecular dynamics simulations on N-palmitoyl and N-stearoyl sphingomyelin. NMR data on these two molecules from several studies in bilayers and micelles played a strong role in the development and testing of the force field parameters. Most previous force fields for sphingomyelins were developed before the availability of the detailed NMR data and relied on x-ray diffraction of bilayers alone for the validation; these are shown to be too dense in the bilayer plane based on published chain order parameter data from simulations and experiments. The present simulations reveal O-H:::O-P intralipid hydrogen bonding occurs 99% of the time, and interlipid N-H:::O=C (26-29%, depending on the lipid) and N-H:::O-H (17-19%). The interlipid hydrogen bonds are long lived, showing decay times of 50 ns, and forming strings of lipids, and leading to reorientational correlation time of nearly 100 ns. The spontaneous radius of curvature for pure N-palmitoyl sphingomyelin bilayers is estimated to be 43-100 Å, depending on the assumptions made in assigning a bending constant; this unusual positive curvature for a two-tailed neutral lipid is likely associated with hydrogen bond networks involving the NH of the sphingosine group.

Pressure-temperature phase behavior of mixtures of natural sphingomyelin and ceramide extracts

Ceramides are a group of sphingolipids that act as highly important signaling molecules in a variety of cellular processes including differentiation and apoptosis. The predominant in vivo synthetic pathway for ceramide formation is via sphingomyelinase catalyzed hydrolysis of sphingomyelin. The biochemistry of this essential pathway has been studied in detail; however, there is currently a lack of information on the structural behavior of sphingomyelin- and ceramide-rich model membrane systems, which is essential for developing a bottom-up understanding of ceramide signaling and platform formation. We have studied the lyotropic phase behavior of sphingomyelin-ceramide mixtures in excess water as a function of temperature (30-70 °C) and pressure (1-200 MPa) by small- and wide-angle X-ray scattering. At low ceramide concentrations the mixtures form the ripple gel phase (P(β)') below the gel transition temperature for sphingomyelin, and this observation has been confirmed by atomic force microscopy. Formation of the ripple gel phase can also be induced at higher temperatures via the application of hydrostatic pressure. At high ceramide concentration an inverse hexagonal phase (HII) is formed coexisting with a cubic phase.

Coarse-grain simulations of skin ceramide NS with newly derived parameters clarify structure of melted phase

Ceramides are lipids that are involved in numerous biologically important structures (e.g., the stratum corneum and ceramide-rich platforms) and processes (e.g., signal transduction and membrane fusion), but their behavior is not fully understood. We report coarse-grain force field parameters for N-lignocerylsphingosine (ceramide NS, also known as ceramide 2) that are consistent with the Martini force field. These parameters were optimized for simulations in the gel phase and validated against atomistic simulations. Coarse-grained simulations with our parameters provide areas per lipid, membrane thicknesses, and electron density profiles that are in good agreement with atomistic simulations. Properties of the simulated membranes are compared with available experimental data. The obtained parameters were used to model the phase behavior of ceramide NS as a function of temperature and hydration. At low water content and above the main phase transition temperature, the bilayer melts into an irregular phase, which may correspond to the unstructured melted-chain phase observed in X-ray diffraction experiments. The developed parameters also reproduce the extended conformation of ceramide, which may occur in the stratum corneum. The parameters presented herein will facilitate studies on important complex functional structures such as the uppermost layer of the skin and ceramide-rich platforms in phospholipid membranes.

Cholesterol stimulates and ceramide inhibits Sticholysin II-induced pore formation in complex bilayer membranes

The pore forming capacity of Sticholysin II (StnII; isolated from Stichodactyla helianthus) in bilayer membranes containing 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphocholine (POPC), palmitoylsphingomyelin (PSM) and either cholesterol or palmitoyl ceramide (PCer) has been examined. The aim of the study was to elucidate how the presence of differently ordered PSM domains affected StnII oligomerization and pore formation. Cholesterol is known to enhance pore formation by StnII, and our results confirmed this and provide kinetic information for the process. The effect of cholesterol on bilayer permeabilization kinetics was concentration-dependent. In the concentration regime used (2.5-10nmol cholesterol in POPC:PSM 80:20 by nmol), cholesterol also increased the acyl chain order in the fluid PSM domain and thus decreased bilayer fluidity, suggesting that fluidity per se was not responsible for cholesterol's effect. Addition of PCer (2.5-10nmol) to the POPC:PSM (80:20 by nmol) bilayers attenuated StnII-induced pore formation, again in a concentration-dependent fashion. This addition also led to the formation of a PCer-rich gel phase. Addition of cholesterol to PCer-containing membranes could partially reduce the inhibitory effect of PCer on StnII pore formation. We conclude that the physical state of PSM (as influenced by either cholesterol or PCer) affected StnII binding and pore formation under the conditions examined.

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.

Simulating G protein-coupled receptors in native-like membranes: from monomers to oligomers

G protein-coupled receptors (GPCRs) are one of the most relevant superfamilies of transmembrane proteins as they participate in an important variety of biological events. Recently, the scientific community is witnessing an advent of a GPCR crystallization age along with impressive improvements achieved in the field of computer simulations during the last two decades. Computer simulation techniques such as molecular dynamics (MD) simulations are now frequent tools to study the dynamic behavior of GPCRs and, more importantly, to model the complex membrane environment where these proteins spend their lifetime. Thanks to these tools, GPCRs can be simulated not only longer but also in a more "physiological" fashion. In this scenario, scientists are taking advantage of such advances to approach certain phenomena such as GPCR oligomerization occurring only at timescales not reachable until now. Thus, despite current MD simulations having important limitations today, they have become an essential tool to study key biophysical properties of GPCRs and GPCR oligomers.

Ceramide-lipid interactions studied by MD simulations and solid-state NMR

Ceramides play a key modulatory role in many cellular processes, which results from their effect on the structure and dynamics of biological membranes. In this study, we investigate the influence of C16-ceramide (C16) on the biophysical properties of DMPC lipid bilayers using solid-state NMR and atomistic molecular dynamics (MD) simulations. MD simulations and NMR measurements were carried out for a pure DMPC bilayer and for a 20% DMPC-C16 mixture. Calculated key structural properties, namely area per lipid, chain order parameters, and mass density profiles, indicate that C16 has an ordering effect on the DMPC bilayer. Furthermore, the simulations predict that specific hydrogen-bonds form between DMPC and C16 molecules. Multi-nuclear solid-state NMR was used to verify these theoretical predictions. Chain order parameters extracted from (13)C(1)H dipole couplings were measured for both lipid and ceramide and follow the trend suggested by the MD simulations. Furthermore, (1)H-MAS NMR experiments showed a direct contact between ceramide and lipids.

Simulation study of the structure and phase behavior of ceramide bilayers and the role of lipid head group chemistry

Ceramides are known to be a key component of the stratum corneum, the outermost protective layer of the skin that controls barrier function. In this work, molecular dynamics simulations are used to examine the behavior of ceramide bilayers, focusing on non-hydroxy sphingosine (NS) and non-hydroxy phytosphingosine (NP) ceramides. Here, we propose a modified version of the CHARMM force field for ceramide simulation, which is directly compared to the more commonly used GROMOS-based force field of Berger (Biophys. J. 1997, 72); while both force fields are shown to closely match experiment from a structural standpoint at the physiological temperature of skin, the modified CHARMM force field is better able to capture the thermotropic phase transitions observed in experiment. The role of ceramide chemistry and its impact on structural ordering is examined by comparing ceramide NS to NP, using the validated CHARMM-based force field. These simulations demonstrate that changing from ceramide NS to NP results in changes to the orientation of the OH groups in the lipid headgroups. The arrangement of OH groups perpendicular to the bilayer normal for ceramide NP, verse parallel for NS, results in the formation of a distinct hydrogen bonding network, that is ultimately responsible for shifting the gel-to-liquid phase transition to higher temperature, in direct agreement with experiment.

Effects of sphingomyelin headgroup size on interactions with ceramide

Sphingomyelins (SMs) and ceramides are known to interact favorably in bilayer membranes. Because ceramide lacks a headgroup that could shield its hydrophobic body from unfavorable interactions with water, accommodation of ceramide under the larger phosphocholine headgroup of SM could contribute to their favorable interactions. To elucidate the role of SM headgroup for SM/ceramide interactions, we explored the effects of reducing the size of the phosphocholine headgroup (removing one, two, or three methyls on the choline moiety, or the choline moiety itself). Using differential scanning calorimetry and fluorescence spectroscopy, we found that the size of the SM headgroup had no marked effect on the thermal stability of ordered domains formed by SM analog/palmitoyl ceramide (PCer) interactions. In more complex bilayers composed of a fluid glycerophospholipid, SM analog, and PCer, the thermal stability and molecular order of the laterally segregated gel domains were roughly identical despite variation in SM headgroup size. We suggest that that the association between PCer and SM analogs was stabilized by ceramide's aversion for disordered phospholipids, by interfacial hydrogen bonding between PCer and the SM analogs, and by attractive van der Waals' forces between saturated chains of PCer and SM analogs.

Cholesteryl phosphocholine--a study on its interactions with ceramides and other membrane lipids

We prepared cholesteryl phosphocholine (CholPC) by chemical synthesis and studied its interactions with small (ceramide and cholesterol) and large headgroup (sphingomyelin (SM) and phosphatidylcholine) colipids in bilayer membranes. We established that CholPC could form bilayers (giant uni- and multilamellar vesicles, as well as extruded large unilamellar vesicles) with both ceramides and cholesterol (initial molar ratio 1:1). The extruded bilayers appeared to be fluid, although highly ordered, even when the ceramide had an N-linked palmitoyl acyl chain. In binary systems containing CholPC and either palmitoyl SM or 1,2-dipalmitoyl-sn-glycero-3-phospholine, CholPC markedly destabilized the gel phase of the respective large headgroup lipid. In 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphocholine (POPC) bilayers, CholPC was much less efficient than cholesterol in ordering the acyl chains. In complex bilayers containing POPC and cholesterol or palmitoyl ceramide, CholPC appeared to prefer interacting with the small headgroup lipids over POPC. When the degree of order in CholPC/PCer bilayers was compared to Chol/PSM bilayers, CholPC/PCer bilayers were more disordered (based on DPH anisotropy). This finding may result from different headgroup orientation and dynamics in CholPC and PSM. Our results overall can be understood if one takes into account the molecular shape of CholPC (large polar headgroup and modest size hydrophobic part) when interpreting molecular interactions between small and large headgroup colipids. The results are also consistent with the proposed umbrella model" for explaining cholesterol/colipid interactions.

Insights into sphingolipid miscibility: separate observation of sphingomyelin and ceramide N-acyl chain melting

Ceramide produced from sphingomyelin in the plasma membrane is purported to affect signaling through changes in the membrane's physical properties. Thermal behavior of N-palmitoyl sphingomyelin (PSM) and N-palmitoyl ceramide (PCer) mixtures in excess water has been monitored by ²H NMR spectroscopy and compared to differential scanning calorimetry (DSC) data. The alternate use of either perdeuterated or proton-based N-acyl chain PSM and PCer in our ²H NMR studies has allowed the separate observation of gel-fluid transitions in each lipid in the presence of the other one, and this in turn has provided direct information on the lipids' miscibility over a wide temperature range. The results provide further evidence of the stabilization of the PSM gel state by PCer. Moreover, overlapping NMR and DSC data reveal that the DSC-signals parallel the melting of the major component (PSM) except at intermediate (20 and 30 mol %) fractions of PCer. In such cases, the DSC endotherm reports on the presumably highly cooperative melting of PCer. Up to at least 50 mol % PCer, PSM and PCer mix ideally in the liquid crystalline phase; in the gel phase, PCer becomes incorporated into PSM:PCer membranes with no evidence of pure solid PCer.

Another Piece of the Membrane Puzzle: Extending Slipids Further

To be able to model complex biological membranes in a more realistic manner, the force field Slipids (Stockholm lipids) has been extended to include parameters for sphingomyelin (SM), phosphatidylglycerol (PG), phosphatidylserine (PS) lipids, and cholesterol. Since the parametrization scheme was faithful to the scheme used in previous editions of Slipids, all parameters are consistent and fully compatible. The results of careful validation of a number of key structural properties for one and two component lipid bilayers are in excellent agreement with experiments. Potentials of mean force for transferring water across binary mixtures of lipids and cholesterol were also computed in order to compare water permeability rates to experiments. In agreement with experimental and simulation studies, it was found that the permeability and partitioning of water is affected by cholesterol in lipid bilayers made of saturated lipids to the largest extent. With the extensions of Slipids presented here, it is now possible to study complex systems containing many different lipids and proteins in a fully atomistic resolution in the isothermic-isobaric (NPT) ensemble, which is the proper ensemble for membrane simulations.