Structure of sphingomyelin bilayers: a simulation study

Unravelling the influence of surface lipids on the structure, dynamics and interactome of high-density lipoproteins

Surface lipids influence the biological activities of high-density lipoproteins (HDLs) but their species-specific effects on HDL structure, dynamics, and surface interactome has remained unclear. Building upon the five-lipid species HDL models developed and characterised in previous work, representative models of the major HDL subpopulations found in human plasma containing apolipoprotein A-I (apoA-I) have been studied using molecular dynamics simulation to describe their varying degrees of surface lipidome complexity. Specifically, two additional sets of representative HDL subpopulation particles were developed, one with sphingomyelin (SM) and the other with SM, phosphatidylethanolamine, phosphatidylinositol, and ceramide in quantities reflecting average levels characterised for HDL subpopulations derived from normolipidemic patients. These lipid species were assessed in terms of HDL size, morphology, dynamics, and overall interactome. The findings reveal that the presence of a representative SM fraction marginally enhanced HDL interfacial curvature and surface monolayer rigidity, manifesting in tighter phospholipid packing and slower surface lipid dynamics relative to SM-deficient HDL models. Furthermore, the presence of SM resulted in a reduction in the solvent exposure of core lipids and cholesterol molecules, whilst also enhancing apolipoprotein conformational flexibility and its overall twisting across the HDL surface. The hydrophobicity of apoA-I-bound lipid patches and the proportion of apoA-I hydrophobic surface area is enhanced by the overall lipidation of apoA-I irrespective of lipid composition. These findings offer new insights into how the surface lipid composition of different HDL subpopulations can significantly impact the overall interactome of HDL particles, potentially influencing subpopulation-specific biological functions like lipid scavenging and receptor interactions.

Lipid nanoparticle mRNA systems containing high levels of sphingomyelin engender higher protein expression in hepatic and extra-hepatic tissues

Lipid nanoparticles (LNPs) for delivery of mRNA usually contain ionizable lipid/helper lipid/cholesterol/PEG-lipid in molar ratios of 50:10:38.5:1.5, respectively. These LNPs are rapidly cleared from the circulation following intravenous (i.v.) administration, limiting uptake into other tissues. Here, we investigate the properties of LNP mRNA systems prepared with high levels of "helper" lipids such as 1,2-distearoyl-sn-glycero-3-phosphorylcholine (DSPC) or N-(hexadecanoyl)-sphing-4-enine-1-phosphocholine (egg sphingomyelin [ESM]). We show that LNP mRNAs containing 40 mol % DSPC or ESM have a unique morphology with a small interior "solid" core situated in an aqueous compartment that is bounded by a lipid bilayer. The encapsulated mRNA exhibits enhanced stability in the presence of serum. LNP mRNA systems containing 40 mol % DSPC or ESM exhibit significantly improved transfection properties in vitro compared with systems containing 10 mol % DSPC or ESM. When injected i.v., LNP mRNAs containing 40 mol % ESM exhibit extended circulation lifetimes compared with LNP mRNA systems containing 10 mol % DSPC, resulting in improved accumulation in extrahepatic tissues. Systems containing 40 mol % ESM result in significantly improved gene expression in spleen and bone marrow as well as liver post i.v. injection compared with 10 mol % DSPC LNP mRNAs. We conclude that LNP mRNAs containing high levels of helper lipid provide a new approach for transfecting hepatic and extrahepatic tissues.

Study of lipid heterogeneity on bilayer membranes using molecular dynamics simulations

Human cell membranes consist of various lipids that are essential for their structure and function. It typically comprises phosphatidylcholine (POPC), phosphatidylethanolamine (POPE), phosphatidylserine (POPS), sphingomyelin (PSM), and cholesterol (CHL). Several experimental and computational techniques have been employed to characterize the composition of human cell membranes, however, CHL enriched membrane is still not clearly understood through these techniques. Molecular dynamics simulation results illustrated the biophysical properties of heterogeneous membranes based on the lipid composition as well as the concentration of lipids, exclusively for CHL and PSM. Herein, we have investigated the structure-function relationships of lipids comparatively to delineate the effect of heterogeneity on the biophysical properties of different membranes. It has been observed that the significant fraction of CHL (i.e., ~33% in ternary, ~25% in quaternary, and ~16% in senary type bilayers) in combination with other lipids introduced compactness, and increased the thickness of the membrane. The analysis of lipid mass density stated that the density of lipid head group, phosphate_, and glycerol-ester in presence of CHL with or without PSM is an underlying reason for membrane ordering. Results also revealed that the presence of POPI and POPS are the reasons for an adequate drop in the ordering of lipid chain, particularly on POPE chain. The self-interaction of CHL, PSM, POPE and the interaction of CHL and POPC with POPE seem to determine the structure and function of the heterogeneous membrane. Our findings provide a qualitative understanding of the effect of membrane heterogeneity on the physiological properties of membranes. The structures inspected in this study would help to select the heterogeneous bilayer model to mimic the human cell membranes to analyse or characterize the membrane-associated phenomena.

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.

Lipid interactions of an actinoporin pore-forming oligomer

The actinoporins are cytolytic toxins produced by sea anemones. Upon encountering a membrane, preferably containing sphingomyelin, they oligomerize and insert their N-terminal helix into the membrane, forming a pore. Whether sphingomyelin is specifically recognized by the protein or simply induces phase coexistence in the membrane has been debated. Here, we perform multi-microsecond molecular dynamics simulations of an octamer of fragaceatoxin C, a member of the actinoporin family, in lipid bilayers containing either pure 1,2-dioleoyl-sn-glycero-3-phosphocholine (DOPC) or a 1:1 mixture of DOPC and palmitoyl sphingomyelin (PSM). The complex is highly stable in both environments, with only slight fraying of the inserted helices near their N-termini. Analyzing the structural parameters of the mixed membrane in the course of the simulation, we see signs of a phase transition for PSM in the inner leaflet of the bilayer. In both leaflets, cross-interactions between lipids of different type decrease over time. Surprisingly, the aromatic loop thought to be responsible for sphingomyelin recognition interacts more with DOPC than PSM by the end of the simulation. These results support the notion that the key membrane property that actinoporins recognize is lipid phase coexistence.

Protective Role of Sphingomyelin in Eye Lens Cell Membrane Model against Oxidative Stress

In the eye lens cell membrane, the lipid composition changes during the aging process: the proportion of sphingomyelins (SM) increases, that of phosphatidylcholines decreases. To investigate the protective role of the SMs in the lens cell membrane against oxidative damage, analytical techniques such as electrochemistry, high-resolution mass spectrometry (HR-MS), and atomic force microscopy (AFM) were applied. Supported lipid bilayers (SLB) were prepared to mimic the lens cell membrane with different fractions of PLPC/SM (PLPC: 1-palmitoyl-2-linoleoyl-sn-glycero-3-phosphocholine). The SLBs were treated with cold physical plasma. A protective effect of 30% and 44% in the presence of 25%, and 75% SM in the bilayer was observed, respectively. PLPC and SM oxidation products were determined via HR-MS for SLBs after plasma treatment. The yield of fragments gradually decreased as the SM ratio increased. Topographic images obtained by AFM of PLPC-bilayers showed SLB degradation and pore formation after plasma treatment, no degradation was observed in PLPC/SM bilayers. The results of all techniques confirm the protective role of SM in the membrane against oxidative damage and support the idea that the SM content in lens cell membrane is increased during aging in the absence of effective antioxidant systems to protect the eye from oxidative damage and to prolong lens transparency.

Additional contributions to elastic energy of lipid membranes: Tilt-curvature coupling and curvature gradient

Lipid bilayer membranes under biologically relevant conditions are flexible thin laterally fluid films consisting of two unimolecular layers (monolayers) each about 2 nm thick. On spatial scales much larger than the bilayer thickness, the membrane elasticity is well determined by its shape. The classical Helfrich theory considers the membrane as an elastic two-dimensional (2D) film, which has no particular internal structure. However, various local membrane heterogeneities can result in a lipids tilt relative to the membrane surface normal. On the basis of the classical elasticity theory of 3D bodies, Hamm and Kozlov [Eur. Phys. J. E 3, 323 (2000)10.1007/s101890070003] derived the most general energy functional, taking into account the tilt and lipid monolayer curvature. Recently, Terzi and Deserno [J. Chem. Phys. 147, 084702 (2017)10.1063/1.4990404] showed that Hamm and Kozlov's derivation was incomplete because the tilt-curvature coupling term had been missed. However, the energy functional derived by Terzi and Deserno appeared to be unstable, thereby being invalid for applications that require minimizations of the overall energy of deformations. Here, we derive a stable elastic energy functional, showing that the squared gradient of the curvature was missed in both of these works. This change in the energy functional arises from a more accurate consideration of the transverse shear deformation terms and their influence on the membrane stability. We also consider the influence of the prestress terms on the stability of the energy functional, and we show that it should be considered small and the effective Gaussian curvature should be neglected because of the stability requirements. We further generalize the theory, including the stretching-compressing deformation modes, and we provide the geometrical interpretation of the terms that were previously missed by Hamm and Kozlov. The physical consequences of the new terms are analyzed in the case of a membrane-mediated interaction of two amphipathic peptides located in the same monolayer. We also provide the expression for director fluctuations, comparing it with that obtained by Terzi and Deserno.

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 dynamics study of lipid bilayers modeling outer and inner leaflets of plasma membranes of mouse hepatocytes. I. Differences in physicochemical properties between the two leaflets

Outer and inner leaflets of plasma cell membranes have different lipid compositions, and the membrane properties of each leaflet can differ from each other significantly due to these composition differences. However, because of the experimental difficulty in measuring the membrane properties for each leaflet separately, the differences are not well understood at a molecular level. In this study, we constructed two lipid bilayer systems, modeling outer and inner leaflets of plasma membranes of mouse hepatocytes based on experimental composition data. The ion concentration in the interlamellar water phase was also set to match the concentration in extra- and intracellular fluids. The differences in physical properties between the outer and inner leaflets of mouse hepatocyte cell membrane models were investigated by performing 1.2 μs-long all-atomistic molecular dynamics calculations under physiological temperature and pressure conditions (310.15 K and 1 atm). The calculated electron density profiles along the bilayer normal for each model bilayer system captured well the asymmetric feature of the experimental electron density profile across actual cell plasma membranes, indicating that our procedure of modeling the outer and inner leaflets of the cell plasma membranes was satisfactory. We found that compared to the outer leaflet model, the inner leaflet model had a very bulky and soft structure in the lateral direction. To confirm the differences, membrane fluidity was measured from the lateral diffusivity and relaxation times. The fluidity was significantly higher in the inner leaflet model than in the outer leaflet model. We also discuss two topics that are of wide interest in biology, i.e., the interdigitation of acyl tails of lipid molecules between two monolayers and the lateral concentration fluctuation of lipid molecules in the bilayers.

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.

Normal and shear forces between boundary sphingomyelin layers under aqueous conditions

Sphingomyelin is one of the predominant phospholipid groups in synovial joints, where lipids have been strongly implicated in the boundary lubrication of articular cartilage; however, little attention has been paid to its lubrication behavior. In this study, we demonstrate that sphingomyelin is an excellent boundary lubricant by measuring the normal and shear forces between sphingomyelin-layer-coated surfaces with a surface force balance under aqueous conditions. Slightly negatively-charged egg sphingomyelin vesicles were adsorbed on mica either by calcium bridging or by charge screening with high concentration monovalent salt. The normal force profiles between opposing egg sphingomyelin layers (vesicles or bilayers) show long-ranged weak repulsion and short-ranged strong repulsion on approaching. Friction coefficients, calculated from the highest load, were (7.2 ± 1.7) × 10-4 at contact stresses of 9.1 ± 0.7 MPa across 0.3 mM liposome dispersion in 0.03 mM Ca2+, and (0.8-3.5) × 10-3 at contact stresses of 7.6 ± 0.8 MPa across 0.3 mM liposome dispersion in 150 mM NaNO3. Similar or slightly lower friction coefficients of (5.3 ± 0.8) × 10-4 at 9.8 ± 0.2 MPa were obtained by replacing the liposome dispersion in 0.03 mM Ca2+ by water. Such low friction coefficients, attributed to the hydration lubrication mechanism, are comparable to those of phosphatidylcholine lipids, which have been widely recognized as excellent aqueous biolubricants. Therefore, we believe that sphingomyelin, in parallel with phosphatidylcholine, contributes to the remarkably good boundary lubrication in synovial joints.

Sphingomyelins Prevent Propagation of Lipid Peroxidation-LC-MS/MS Evaluation of Inhibition Mechanisms

Free radical driven lipid peroxidation is a chain reaction which can lead to oxidative degradation of biological membranes. Propagation vs. termination rates of peroxidation in biological membranes are determined by a variety of factors including fatty acyl chain composition, presence of antioxidants, as well as biophysical properties of mono- or bilayers. Sphingomyelins (SMs), a class of sphingophospholipids, were previously described to inhibit lipid oxidation most probably via the formation of H-bond network within membranes. To address the “antioxidant” potential of SMs, we performed LC-MS/MS analysis of model SM/glycerophosphatidylcholine (PC) liposomes with different SM fraction after induction of radical driven lipid peroxidation. Increasing SM fraction led to a strong suppression of lipid peroxidation. Electrochemical oxidation of non-liposomal SMs eliminated the observed effect, indicating the importance of membrane structure for inhibition of peroxidation propagation. High resolution MS analysis of lipid peroxidation products (LPPs) observed in in vitro oxidized SM/PC liposomes allowed to identify and relatively quantify SM- and PC-derived LPPs. Moreover, mapping quantified LPPs to the known pathways of lipid peroxidation allowed to demonstrate significant decrease in mono-hydroxy(epoxy) LPPs relative to mono-keto derivatives in SM-rich liposomes. The results presented here illustrate an important property of SMs in biological membranes, acting as “biophysical antioxidant”. Furthermore, a ratio between mono-keto/mono-hydroxy(epoxy) oxidized species can be used as a marker of lipid peroxidation propagation in the presence of different antioxidants.

Fat SIRAH: Coarse-Grained Phospholipids To Explore Membrane-Protein Dynamics

The capability to handle highly heterogeneous molecular assemblies in a consistent manner is among the greatest challenges faced when deriving simulation parameters. This is particularly the case for coarse-grained (CG) simulations in which chemical functional groups are lumped into effective interaction centers for which transferability between different chemical environments is not guaranteed. Here, we introduce the parametrization of a set of CG phospholipids compatible with the latest version of the SIRAH force field for proteins. The newly introduced lipid species include different acylic chain lengths and partial unsaturation, as well as polar and acidic head groups that show a very good reproduction of structural membrane determinants, such as areas per lipid, thickness, order parameter, etc., and their dependence with temperature. Simulation of membrane proteins showed unprecedented accuracy in the unbiased description of the thickness-dependent membrane-protein orientation in systems where this information is experimentally available (namely, the SarcoEndoplasmic Reticulum Calcium-SERCA-pump and its regulator Phospholamban). The interactions that lead to this faithful reproduction can be traced down to the single amino acid-lipid interaction level and show full agreement with biochemical data present in the literature. Finally, the present parametrization is implemented in the GROMACS and AMBER simulation packages facilitating its use by a wide portion of the biocomputing community.

Lipid-Surrounding Water Molecules Probed by Time-Resolved Emission Spectra of Laurdan

The hydration states of the interfacial region of lipid bilayers were investigated on the basis of the time-resolved emission spectra (TRES) analysis of 6-lauroyl-2-dimethylamino naphthalene (Laurdan), a common fluorescence probe used to analyze membrane hydration. TRES derived from long and short lifetime components were extracted from samples of different lipid species: 1,2-dipalmitoyl- sn-glycero-3-phosphocholine (DPPC), 1,2-dioleoyl- sn-glycero-3-phosphocholine (DOPC), d- erythro- N-palmitoyl-sphingosylphosphorylcholine (PSM), and a DOPC/PSM binary bilayer system. Neither lifetime component (short or long) corresponded with the hydration properties; the short lifetime component of DOPC (1.97 ns) exhibited a peak at 440 nm, and the long lifetime components of DPPC and PSM (7.76 and 7.77 ns, respectively) exhibited peaks at the same wavelength. This similarity arose from the competition between the collisional quenching and the hydration effects of water molecules. Herein, this phenomenon was investigated using a plot of the lifetime τ and the peak position λ (τ vs λ plot), simultaneously visualizing both effects by deconvoluting the TRES. On the basis of collisional quenching theory, the distribution of the water population per lipid (water map) was generated. According to this theory, the τ vs λ plot was applied to the water map and the calculation of the number of water molecules per lipid, which is consistent with previous reports. This approach provides novel insights for the analysis of molecular hydration states using the fluorescence of Laurdan.

Mechanical properties of bilayers containing sperm sphingomyelins and ceramides with very long-chain polyunsaturated fatty acids

Sphingomyelins (SM) and ceramides (Cer) with very long chain polyunsaturated fatty acids (V) are important components of spermatozoa membranes. In this study, the mechanical properties of bilayers of SM and Cer with nonhydroxy (n-V) and 2-hydroxy (h-V) fatty acid (30:5) were studied by molecular dynamics simulation at different temperatures and in the presence and the absence of salt. From our results, it was evidenced how n-V SM and h-V SM bilayers showed similar behavior. When n-V Cer was added to a h-V SM bilayer, the Gaussian curvature modulus and E_curve of binary bilayers decreased. This variation in the mechanical properties of the bilayer can be associated with an incipient step during the fecundation process.

A simple model of planar membrane: An integral equation investigation

A simple model of a lipid membrane, a binary mixture of saturated lipids and unsaturated lipids, was studied using an integral equation theory. The planar membrane is modeled as mixture of linear and bent molecules in two-dimensional space, and site-site radial distribution function, Kirkwood-Buff (KB) integral and related quantities were computed over the whole range of the molar fraction to understand their mixing behavior. We found that a close packing of linear molecules is enhanced as the fraction of bent molecules increases, but a long range correlation between the linear molecules is weakened. A high concentration of linear molecules promotes the demixing of linear molecules and bent molecules, and enhances the long range correlation between molecules. This implies that, the higher the concentration of linear molecules, the larger clusters tend to be formed. © 2018 Wiley Periodicals, Inc.

Atomistic Picture of Fluorescent Probes with Hydrocarbon Tails in Lipid Bilayer Membranes: An Investigation of Selective Affinities and Fluorescent Anisotropies in Different Environmental Phases

By reverting to spectroscopy, changes in the biological environment of a fluorescent probe can be monitored and the presence of various phases of the surrounding lipid bilayer membranes can be detected. However, it is currently not always clear in which phase the probe resides. The well-known orange 1,1'-dioctadecyl-3,3,3',3'-tetramethylindodicarbo-cyanine perchlorate (DiI-C18(5)) fluorophore, for instance, and the new, blue BODIPY (4,4-difluoro-4-bora-3 a,4 a-diaza- s-indacene) derivative were experimentally seen to target and highlight identical parts of giant unilamellar vesicles of various compositions, comprising mixtures of dipalmitoylphosphatidylcholine (DPPC), dioleoylphosphatidylcholine (DOPC), sphingomyelin (SM), and cholesterol (Chol). However, it was not clear which of the coexisting membrane phases were visualized (Bacalum et al., Langmuir. 2016, 32, 3495). The present study addresses this issue by utilizing large-scale molecular dynamics simulations and the z-constraint method, which allows evaluating Gibbs free-energy profiles. The current calculations give an indication why, at room temperature, both BODIPY and DiI-C18(5) probes prefer the gel (S_o) phase in DOPC/DPPC (2:3 molar ratio) and the liquid-ordered (L_o) phase in DOPC/SM/Chol (1:2:1 molar ratio) mixtures. This study highlights the important differences in orientation and location and therefore in efficiency between the probes when they are used in fluorescence microscopy to screen various lipid bilayer membrane phases. Dependent on the lipid composition, the angle between the transition-state dipole moments of both probes and the normal to the membrane is found to deviate clearly from 90°. It is seen that the DiI-C18(5) probe is located in the headgroup region of the SM/Chol mixture, in close contact with water molecules. A fluorescence anisotropy study also indicates that DiI-C18(5) gives rise to a distinctive behavior in the SM/Chol membrane compared to the other considered membranes. The latter behavior has not been seen for the studied BODIPY probe, which is located deeper in the membrane.

Sphingomyelin ability to act as chiral selector using nanodisc electrokinetic chromatography

Understanding how stereochemistry affects interactions with cell membranes is important for effective drug development. Chirality has been shown to greatly effect pharmaceutical distribution and metabolism within the cell. However it has been thought that interactions with, and passive diffusion through, the membrane are not stereochemically selective. Various studies have produced conflicting results regarding whether interactions with lipid bilayers are or can be stereoselective. In the current work, stereoselective interactions between a pair of atropisomers, R-(+)/(S)-(-) 1,1'-Bi-2-naphthol, and sphingomyelin nanodisc bilayers, are demonstrated. This is accomplished using nanodisc electrokinetic chromatography, demonstrating that this approach is sensitive to subtle differences in affinity between small molecule probes and lipid bilayers. Using the same approach, no evidence of stereoselectivity was observed using enantiomer or diastereomer probes of varied chemistry and structure.

Atomistic Simulations of Electroporation of Model Cell Membranes

Electroporation is a phenomenon that modifies the fundamental function of the cell since it perturbs transiently or permanently the integrity of its membrane. Today, this technique is applied in fields ranging from biology and biotechnology to medicine, e.g., for drug and gene delivery into cells, tumor therapy, etc., in which it made it to preclinical and clinical treatments. Experimentally, due to the complexity and heterogeneity of cell membranes, it is difficult to provide a description of the electroporation phenomenon in terms of atomically resolved structural and dynamical processes, a prerequisite to optimize its use. Atomistic modeling in general and molecular dynamics (MD) simulations in particular have proven to be an effective approach for providing such a level of detail. This chapter provides the reader with a comprehensive account of recent advances in using such a technique to complement conventional experimental approaches in characterizing several aspects of cell membranes electroporation.

Determination of lipid bilayer affinities and solvation characteristics by electrokinetic chromatography using polymer-bound lipid bilayer nanodiscs

Styrene-maleic acid polymer-bound lipid bilayer nanodiscs have been investigated and characterized by electrokinetic chromatography. Linear solvation energy relationship analysis was employed to characterize the changes in solvation environment of nanodiscs of varied belt to lipid ratio, belt polymer chemistry and molecular weight, and lipid composition. Increases in the lipid to belt polymer ratio resulted in smaller, more cohesive nanodiscs with greater electrophoretic mobility. Nanodisc structures with belt polymers of different chemistry and molecular weight were compared and showed only minor changes in solvent characteristics and selectivity consistent with changes in structure of the lipid bilayer. Seven phospholipid and sphingomyelin nanodiscs of different lipid composition were characterized. Changes in lipid head group structure had a significant effect on bilayer-solute interactions. In most cases, changes in alkyl tail structure had no discernible effect on solvation environment aside from those explained by changes in the gel-liquid transition temperature. Comparison to vesicles of similar lipid composition show only minor differences in solvation environment, likely due to differences in lipid composition and bilayer curvature. Together these results provide evidence that the dominant solute-nanodisc interactions are with the lipid bilayer and that head group chemistry has a greater impact on bilayer-solute interactions than alkyl tail or belt polymer structure. Nanodisc electrokinetic chromatography is demonstrated to allow characterization of solute interactions with lipid bilayers of varied composition.

Effect of Polycation Structure on Interaction with Lipid Membranes

Interaction of polycations with lipid membranes is a very important issue in many biological and medical applications such as gene delivery or antibacterial usage. In this work, we address the influence of hydrophobic substitution of strong polycations containing quaternary ammonium groups on the polymer-zwitterionic membrane interactions. In particular, we focus on the polymer tendency to adsorb on or/and incorporate into the membrane. We used complementary experimental and computational methods to enhance our understanding of the mechanism of the polycation-membrane interactions. Polycation adsorption on liposomes was assessed using dynamic light scattering (DLS) and zeta potential measurements. The ability of the polymers to form hydrophilic pores in the membrane was evaluated using a calcein-release method. The polymer-membrane interaction at the molecular scale was explored by performing atomistic molecular dynamics (MD) simulations. Our results show that the length of the alkyl side groups plays an essential role in the polycation adhesion on the zwitterionic surface, while the degree of substitution affects the polycation ability to incorporate into the membrane. Both the experimental and computational results show that the membrane permeability can be dramatically affected by the amount of alkyl side groups attached to the polycation main chain.

Molecular dynamics study of lipid bilayers modeling the plasma membranes of mouse hepatocytes and hepatomas

Molecular dynamics (MD) calculations of lipid bilayers modeling the plasma membranes of normal mouse hepatocytes and hepatomas in water have been performed under physiological isothermal-isobaric conditions (310.15 K and 1 atm). The changes in the membrane properties induced by hepatic canceration were investigated and were compared with previous MD calculations included in our previous study of the changes in membrane properties induced by murine thymic canceration. The calculated model membranes for normal hepatocytes and hepatomas comprised 23 and 24 kinds of lipids, respectively. These included phosphatidylcholine, phosphatidylethanolamine, phosphatidylserine, phosphatidylinositol, sphingomyelin, lysophospholipids, and cholesterol. We referred to previously published experimental values for the mole fraction of the lipids adopted in the present calculations. The calculated structural and dynamic properties of the membranes such as lateral structure, order parameters, lateral self-diffusion constants, and rotational correlation times all showed that hepatic canceration causes plasma membranes to become more ordered laterally and less fluid. Interestingly, this finding contrasts with the less ordered structure and increased fluidity of plasma membranes induced by thymic canceration observed in our previous MD study.

Membrane stiffness is modified by integral membrane proteins

The ease with which a cell membrane can bend and deform is important for a wide range of biological functions. Peripheral proteins that induce curvature in membranes (e.g. BAR domains) have been studied for a number of years. Little is known, however, about the effect of integral membrane proteins on the stiffness of a membrane (characterised by the bending rigidity, Kc). We demonstrate by computer simulation that adding integral membrane proteins at physiological densities alters the stiffness of the membrane. First we establish that the coarse-grained MARTINI forcefield is able to accurately reproduce the bending rigidity of a small patch of 1500 phosphatidyl choline lipids by comparing the calculated value to both experiment and an atomistic simulation of the same system. This enables us to simulate the dynamics of large (ca. 50 000 lipids) patches of membrane using the MARTINI coarse-grained description. We find that altering the lipid composition changes the bending rigidity. Adding integral membrane proteins to lipid bilayers also changes the bending rigidity, whilst adding a simple peripheral membrane protein has no effect. Our results suggest that integral membrane proteins can have different effects, and in the case of the bacterial outer membrane protein, BtuB, the greater the density of protein, the larger the reduction in stiffness.

An overview of molecular dynamics simulations of oxidized lipid systems, with a comparison of ELBA and MARTINI force fields for coarse grained lipid simulations

Biological membranes and model lipid systems containing high amounts of unsaturated lipids and sterols are subject to chemical and/or photo-induced lipid oxidation, which leads to the creation of exotic oxidized lipid products (OxPLs). OxPLs are known to have significant physiological impact in cellular systems and also affect physical properties of both biological and model lipid bilayers. In this paper we (i) provide a perspective on the existing literature on simulations of lipid bilayer systems containing oxidized lipid species as well as the main related experimental results, (ii) describe our new data of all-atom and coarse-grained simulations of hydroperoxidized lipid monolayer and bilayer systems and (iii) provide a comparison of the MARTINI and ELBA coarse grained force fields for lipid bilayer systems. We show that the better electrostatic treatment of interactions in ELBA is able to resolve previous conflicts between experiments and simulations. This article is part of a Special Issue entitled: Biosimulations edited by Ilpo Vattulainen and Tomasz Róg.

A molecular insight into the electro-transfer of small molecules through electropores driven by electric fields

The transport of chemical compounds across the plasma membrane into the cell is relevant for several biological and medical applications. One of the most efficient techniques to enhance this uptake is reversible electroporation. Nevertheless, the detailed molecular mechanism of transport of chemical species (dyes, drugs, genetic materials, …) following the application of electric pulses is not yet fully elucidated. In the past decade, molecular dynamics (MD) simulations have been conducted to model the effect of pulsed electric fields on membranes, describing several aspects of this phenomenon. Here, we first present a comprehensive review of the results obtained so far modeling the electroporation of lipid membranes, then we extend these findings to study the electrotransfer across lipid bilayers subject to microsecond pulsed electric fields of Tat11, a small hydrophilic charged peptide, and of siRNA. We use in particular a MD simulation protocol that allows to characterize the transport of charged species through stable pores. Unexpectedly, our results show that for an electroporated bilayer subject to transmembrane voltages in the order of 500mV, i.e. consistent with experimental conditions, both Tat11 and siRNA can translocate through nanoelectropores within tens of ns. We discuss these results in comparison to experiments in order to rationalize the mechanism of drug uptake by cells. This article is part of a Special Issue entitled: Biosimulations edited by Ilpo Vattulainen and Tomasz Róg.

United Atom Lipid Parameters for Combination with the Optimized Potentials for Liquid Simulations All-Atom Force Field

We have developed a new united-atom set of lipid force field parameters for dipalmitoylphosphatidylcholine (DPPC) lipid bilayers that can be combined with the all-atom optimized potentials for liquid simulations (OPLS-AA) protein force field. For this, all torsions have been refitted for a nonbonded 1-4 scale factor of 0.5, which is the standard in OPLS-AA. Improved van der Waals parameters have been obtained for the acyl lipid tails by matching simulation results of bulk pentadecane against recently improved experimental measurements. The charge set has been adjusted from previous lipid force fields to allow for an identical treatment of the alkoxy ester groups. This reduces the amount of parameters required for the model. Simulation of DPPC bilayers in the tension-free NPT ensemble at 50 °C gives the correct area per lipid of 62.9 ± 0.1 Å(2), which compares well with the recently refined experimental value of 63.0 Å(2). Electron density profiles and deuterium order parameters are similarly well reproduced. The new parameters will allow for improved simulation results in microsecond scale peptide partitioning simulations, which have proved problematic with prior parametrizations.

The trans-Golgi SNARE syntaxin 10 is required for optimal development of Chlamydia trachomatis

Chlamydia trachomatis, an obligate intracellular pathogen, grows inside of a vacuole, termed the inclusion. Within the inclusion, the organisms differentiate from the infectious elementary body (EB) into the reticulate body (RB). The RB communicates with the host cell through the inclusion membrane to obtain the nutrients necessary to divide, thus expanding the chlamydial population. At late time points within the developmental cycle, the RBs respond to unknown molecular signals to redifferentiate into infectious EBs to perpetuate the infection cycle. One strategy for Chlamydia to obtain necessary nutrients and metabolites from the host is to intercept host vesicular trafficking pathways. In this study we demonstrate that a trans-Golgi soluble N-ethylmaleimide–sensitive factor attachment protein (SNARE), syntaxin 10, and/or syntaxin 10-associated Golgi elements colocalize with the chlamydial inclusion. We hypothesized that Chlamydia utilizes the molecular machinery of syntaxin 10 at the inclusion membrane to intercept specific vesicular trafficking pathways in order to create and maintain an optimal intra-inclusion environment. To test this hypothesis, we used siRNA knockdown of syntaxin 10 to examine the impact of the loss of syntaxin 10 on chlamydial growth and development. Our results demonstrate that loss of syntaxin 10 leads to defects in normal chlamydial maturation including: variable inclusion size with fewer chlamydial organisms per inclusion, fewer infectious progeny, and delayed or halted RB-EB differentiation. These defects in chlamydial development correlate with an overabundance of NBD-lipid retained by inclusions cultured in syntaxin 10 knockdown cells. Overall, loss of syntaxin 10 at the inclusion membrane negatively affects Chlamydia. Understanding host machinery involved in maintaining an optimal inclusion environment to support chlamydial growth and development is critical toward understanding the molecular signals involved in successful progression through the chlamydial developmental cycle.

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.

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.

Atomistic Monte Carlo simulation of lipid membranes

Biological membranes are complex assemblies of many different molecules of which analysis demands a variety of experimental and computational approaches. In this article, we explain challenges and advantages of atomistic Monte Carlo (MC) simulation of lipid membranes. We provide an introduction into the various move sets that are implemented in current MC methods for efficient conformational sampling of lipids and other molecules. In the second part, we demonstrate for a concrete example, how an atomistic local-move set can be implemented for MC simulations of phospholipid monomers and bilayer patches. We use our recently devised chain breakage/closure (CBC) local move set in the bond-/torsion angle space with the constant-bond-length approximation (CBLA) for the phospholipid dipalmitoylphosphatidylcholine (DPPC). We demonstrate rapid conformational equilibration for a single DPPC molecule, as assessed by calculation of molecular energies and entropies. We also show transition from a crystalline-like to a fluid DPPC bilayer by the CBC local-move MC method, as indicated by the electron density profile, head group orientation, area per lipid, and whole-lipid displacements. We discuss the potential of local-move MC methods in combination with molecular dynamics simulations, for example, for studying multi-component lipid membranes containing cholesterol.

CHARMM36 united atom chain model for lipids and surfactants

Molecular simulations of lipids and surfactants require accurate parameters to reproduce and predict experimental properties. Previously, a united atom (UA) chain model was developed for the CHARMM27/27r lipids (Hénin, J., et al. J. Phys. Chem. B. 2008, 112, 7008-7015) but suffers from the flaw that bilayer simulations using the model require an imposed surface area ensemble, which limits its use to pure bilayer systems. A UA-chain model has been developed based on the CHARMM36 (C36) all-atom lipid parameters, termed C36-UA, and agreed well with bulk, lipid membrane, and micelle formation of a surfactant. Molecular dynamics (MD) simulations of alkanes (heptane and pentadecane) were used to test the validity of C36-UA on density, heat of vaporization, and liquid self-diffusion constants. Then, simulations using C36-UA resulted in accurate properties (surface area per lipid, X-ray and neutron form factors, and chain order parameters) of various saturated- and unsaturated-chain bilayers. When mixed with the all-atom cholesterol model and tested with a series of 1,2-dimyristoyl-sn-glycero-3-phosphocholine (DMPC)/cholesterol mixtures, the C36-UA model performed well. Simulations of self-assembly of a surfactant (dodecylphosphocholine, DPC) using C36-UA suggest an aggregation number of 53 ± 11 DPC molecules at 0.45 M of DPC, which agrees well with experimental estimates. Therefore, the C36-UA force field offers a useful alternative to the all-atom C36 lipid force field by requiring less computational cost while still maintaining the same level of accuracy, which may prove useful for large systems with proteins.

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.

Calculating the Bending Modulus for Multicomponent Lipid Membranes in Different Thermodynamic Phases

We establish a computational approach to extract the bending modulus, KC , for lipid membranes from relatively small-scale molecular simulations. Fluctuations in the splay of individual pairs of lipids faithfully inform on KC in multicomponent membranes over a large range of rigidities in different thermodynamic phases. Predictions are validated by experiments even where the standard spectral analysis-based methods fail. The local nature of this method potentially allows its extension to calculations of KC in protein-laden membranes.

Stiffened lipid platforms at molecular force foci. Proc Natl Acad Sci U S A. 2013;110:4886-4892. [PMID: 23476066 DOI: 10.1073/pnas.1302018110]

How mechanical forces are sensed remains largely mysterious. The forces that gate prokaryotic and several eukaryotic channels were found to come from the lipid membrane. Our survey of animal cells found that membrane force foci all have cholesterol-gathering proteins and are reinforced with cholesterol. This result is evident in overt force sensors at the tips of stereocilia for vertebrate hearing and the touch receptor of Caenorhabditis elegans and mammalian neurons. For less specialized cells, cadherins sustain the force between neighboring cells and integrins between cells and matrix. These tension bearers also pass through and bind to a cholesterol-enriched platform before anchoring to cytoskeleton through other proteins. Cholesterol, in alliance with sphingomyelin and specialized proteins, enforces a more ordered structure in the bilayer. Such a stiffened platform can suppress mechanical noise, redirect, rescale, and confine force. We speculate that such platforms may be dynamic. The applied force may allow disordered-phase lipids to enter the platform-staging channel opening in the thinner mobile neighborhood. The platform may also contain specialized protein/lipid subdomains enclosing mechanosensitive channels to open with localized tension. Such a dynamic stage can mechanically operate structurally disparate channels or enzymes without having to tie them directly to cadherin, integrin, or other protein tethers.

Molecular dynamics study of lipid bilayers modeling the plasma membranes of normal murine thymocytes and leukemic GRSL cells

Molecular dynamics (MD) calculations for the plasma membranes of normal murine thymocytes and thymus-derived leukemic GRSL cells in water have been performed under physiological isothermal-isobaric conditions (310.15K and 1 atm) to investigate changes in membrane properties induced by canceration. The model membranes used in our calculations for normal and leukemic thymocytes comprised 23 and 25 kinds of lipids, respectively, including phosphatidylcholine, phosphatidylethanolamine, phosphatidylserine, phosphatidylinositol, sphingomyelin, lysophospholipids, and cholesterol. The mole fractions of the lipids adopted here were based on previously published experimental values. Our calculations clearly showed that the membrane area was increased in leukemic cells, and that the isothermal area compressibility of the leukemic plasma membranes was double that of normal cells. The calculated membranes of leukemic cells were thus considerably bulkier and softer in the lateral direction compared with those of normal cells. The tilt angle of the cholesterol and the conformation of the phospholipid fatty acid tails both showed a lower level of order in leukemic cell membranes compared with normal cell membranes. The lateral radial distribution function of the lipids also showed a more disordered structure in leukemic cell membranes than in normal cell membranes. These observations all show that, for the present thymocytes, the lateral structure of the membrane is considerably disordered by canceration. Furthermore, the calculated lateral self-diffusion coefficient of the lipid molecules in leukemic cell membranes was almost double that in normal cell membranes. The calculated rotational and wobbling autocorrelation functions also indicated that the molecular motion of the lipids was enhanced in leukemic cell membranes. Thus, here we have demonstrated that the membranes of thymocyte leukemic cells are more disordered and more fluid than normal cell membranes.