Glycolipid Membranes through Atomistic Simulations: Effect of Glucose and Galactose Head Groups on Lipid Bilayer Properties

Lipid/water interface of galactolipid bilayers in different lyotropic liquid-crystalline phases

In this study, carried out using computational methods, the organisation of the lipid/water interface of bilayers composed of galactolipids with both α-linolenoyl acyl chains is analysed and compared in three different lyotropic liquid-crystalline phases. These systems include the monogalactosyldiglyceride (MGDG) and digalactosyldiglyceride (DGDG) bilayers in the lamellar phase, the MGDG double bilayer during stalk phase formation and the inverse hexagonal MGDG phase. For each system, lipid-water and direct and water-mediated lipid-lipid interactions between the lipids of one bilayer leaflet and those of two apposing leaflets at the onset of new phase (stalk) formation, are identified. A network of interactions between DGDG molecules and its topological properties are derived and compared to those for the MGDG bilayer.

The Exploration of the Thermococcus barophilus Lipidome Reveals the Widest Variety of Phosphoglycolipids in Thermococcales

One of the most distinctive characteristics of archaea is their unique lipids. While the general nature of archaeal lipids has been linked to their tolerance to extreme conditions, little is known about the diversity of lipidic structures archaea are able to synthesize, which hinders the elucidation of the physicochemical properties of their cell membrane. In an effort to widen the known lipid repertoire of the piezophilic and hyperthermophilic model archaeon Thermococcus barophilus, we comprehensively characterized its intact polar lipid (IPL), core lipid (CL), and polar head group compositions using a combination of cutting-edge liquid chromatography and mass spectrometric ionization systems. We tentatively identified 82 different IPLs based on five distinct CLs and 10 polar head group derivatives of phosphatidylhexoses, including compounds reported here for the first time, e.g., di-N-acetylhexosamine phosphatidylhexose-bearing lipids. Despite having extended the knowledge on the lipidome, our results also indicate that the majority of T. barophilus lipids remain inaccessible to current analytical procedures and that improvements in lipid extraction and analysis are still required. This expanded yet incomplete lipidome nonetheless opens new avenues for understanding the physiology, physicochemical properties, and organization of the membrane in this archaeon as well as other archaea.

The Anti-Inflammatory Effect of Nanoarchaeosomes on Human Endothelial Cells

Archaebacterias are considered a unique source of novel biomaterials of interest for nanomedicine. In this perspective, the effects of nanoarchaeosomes (ARC), which are nanovesicles prepared from polar lipids extracted from the extreme halophilic Halorubrum tebenquinchense, on human umbilical vein endothelial cells (HUVEC) were investigated in physiological and under inflammatory static conditions. Upon incubation, ARC (170 nm mean size, −41 mV ζ) did not affect viability, cell proliferation, and expression of intercellular adhesion molecule-1 (ICAM-1) and E-selectin under basal conditions, but reduced expression of both molecules and secretion of IL-6 induced by lypopolysaccharide (LPS), Pam3CSK4 or Escherichia coli. Such effects were not observed with TNF-α or IL-1β stimulation. Interestingly, ARC significantly decreased basal levels of von Willebrand factor (vWF) and levels induced by all stimuli. None of these parameters was altered by liposomes of hydrogenated phosphatidylcholine and cholesterol of comparable size and concentration. Only ARC were endocytosed by HUVEC and reduced mRNA expression of ICAM-1 and vWF via NF-ĸB and ERK1/2 in LPS-stimulated cells. This is the first report of the anti-inflammatory effect of ARC on endothelial cells and our data suggest that its future use in vascular disease may hopefully be of particular interest.

Macrophage apoptosis using alendronate in targeted nanoarchaeosomes

Nanoarchaeosomes are non-hydrolysable nanovesicles made of archaeolipids, naturally functionalised with ligand for scavenger receptor class 1. We hypothesized that nitrogenate bisphosphonate alendronate (ALN) loaded nanoarchaeosomes (nanoarchaeosomes(ALN)) may constitute more efficient macrophage targeted apoptotic inducers than ALN loaded nanoliposomes (nanoliposomes (ALN)). To that aim, ALN was loaded in cholesterol containing (nanoARC-chol(ALN)) or not (nanoARC(ALN)) nanoarchaeosomes. Nanoarchaeosomes(ALN) (220-320 nm sized, ~ -40 mV ξ potential, 38-50 μg ALN/mg lipid ratio) displayed higher structural stability than nanoliposomes(ALN) of matching size and ξ potential, retaining most of ALN against a 1/200 folds dilution. The cytotoxicity of nanoARC(ALN) on J774A.1 cells, resulted > 30 folds higher than free ALN and nanoliposomes(ALN) and was reduced by cholesterol in nanoARC-chol(ALN). Devoid of ALN, nanoARC-chol was non-cytotoxic, exhibited pronounced anti-inflammatory activity on J774.1 cells, strongly reducing reactive oxygen species (ROS) and IL-6 induced by LPS. Nanoarchaeosomes bilayer extensively interacted with serum proteins but resulted refractory to phospholipases. Upon J774A.1 cells uptake, nanoarchaeosomes induced cytoplasmic acid vesicles, reduced the mitochondrial membrane potential by 20-40 % without consuming ATP neither damaging lysosomes and increasing pERK. Refractory to chemoenzymatic attacks, either void or drug loaded, nanoarchaeosomes induced either anti-inflammation or macrophages apoptosis, constituting promising targeted nanovesicles for multiple therapeutic purposes.

Effects of lipid heterogeneity on model human brain lipid membranes

Cell membranes naturally contain a heterogeneous lipid distribution. However, homogeneous bilayers are commonly preferred and utilised in computer simulations due to their relative simplicity, and the availability of lipid force field parameters. Recently, experimental lipidomics data for the human brain cell membranes under healthy and Alzheimer's disease (AD) conditions were investigated, since disruption to the lipid composition has been implicated in neurodegenerative disorders, including AD [R. B. Chan et al., J. Biol. Chem., 2012, 287, 2678-2688]. In order to observe the effects of lipid complexity on the various bilayer properties, molecular dynamics simulations were used to study four membranes with increasing heterogeneity: a pure POPC membrane, a POPC and cholesterol membrane in a 1 : 1 ratio (POPC-CHOL), and to our knowledge, the first realistic models of a healthy brain membrane and an Alzheimer's diseased brain membrane. Numerous structural, interfacial, and dynamical properties, including the area per lipid, interdigitation, dipole potential, and lateral diffusion of the two simple models, POPC and POPC-CHOL, were analysed and compared to those of the complex brain models consisting of 27 lipid components. As the membranes gain heterogeneity, a number of alterations were found in the structural and dynamical properties, and more significant differences were observed in the lateral diffusion. Additionally, we observed snorkeling behaviour of the lipid tails that may play a role in the permeation of small molecules across biological membranes. In this work, atomistic description of realistic brain membrane models is provided, which can add insight towards the permeability and transport pathways of small molecules across these membrane barriers.

Supramolecular organization of α-galactosylceramide in pure dispersions and in cationic DODAB bilayers

α-galactosylceramide (α-GalCer; KRN7000) strongly stimulates NKT cells. The structures of α-GalCer assemblies and of cationic DODAB bilayers containing α-GalCer were investigated by differential scanning calorimetry (DSC) and electron spin resonance (ESR) spectroscopy. Assemblies of α-GalCer have a very tightly packed gel phase, causing spin labels to cluster and display spin exchange interactions. An endothermic phase transition is observed by DSC, leading to a fluid phase. This phase transition peak disappears upon mixing with DODAB, showing that up to 9 mol% α-GalCer is miscible with the cationic lipid. ESR spectra show that α-GalCer decreases DODAB gel phase packing, resulting in a decrease of gel-fluid transition temperature and cooperativity in DSC thermograms of mixed bilayers. In contrast, α-GalCer increases the rigidity of the fluid phase. These effects are probably due to the conformation of the rigid amide bond that connects the phytosphingosine base of α-GalCer to its long and saturated acyl chain. Possibly, α-GalCer adopts a V-shaped conformation because of the perpendicular orientation of the amide bond towards the axes of the hydrocarbon chains. Apparently, the effect of the amide bond configuration is a key structural feature for the interaction between ceramide-based glycolipids and DODAB molecules, since we have previously reported a similar decrease of gel phase packing and increase in fluid phase rigidity for DODAB bilayers containing C24:1β-glucosylceramide. Since the structure of delivery systems is critical to the biological activity of α-GalCer, this work certainly contributes to the planning and development of novel immunotherapeutic tools.

Naturally Derived Membrane Lipids Impact Nanoparticle-Based Messenger RNA Delivery

Introduction

Lipid based nanoparticles (LNPs) are clinically successful vectors for hepatic delivery of nucleic acids. These systems are being developed for non-hepatic delivery of mRNA for the treatment of diseases like cystic fibrosis or retinal degeneration as well as infectious diseases. Localized delivery to the lungs requires aerosolization. We hypothesized that structural lipids within LNPs would provide features of integrity which can be tuned for attributes required for efficient hepatic and non-hepatic gene delivery. Herein, we explored whether naturally occurring lipids that originate from the cell membrane of plants and microorganisms enhance mRNA-based gene transfection in vitro and in vivo and whether they assist in maintaining mRNA activity after nebulization.

Methods

We substituted DSPC, a structural lipid used in a conventional LNP formulation, to a series of naturally occurring membrane lipids. We measured the effect of these membrane lipids on size, encapsulation efficiency and their impact on transfection efficiency. We further characterized LNPs after nebulization and measured whether they retained their transfection efficiency.

Results

One plant-derived structural lipid, DGTS, led to a significant improvement in liver transfection of mRNA. DGTS LNPs had similar transfection ability when administered in the nasal cavity to conventional LNPs. In contrast, we found that DGTS LNPs had reduced transfection efficiency in cells pre-and post-nebulization while maintaining size and encapsulation similar to DSPC LNPs.

Conclusions

We found that structural lipids provide differential mRNA-based activities in vitro and in vivo which also depend on the mode of administration. Understanding influence of structural lipids on nanoparticle morphology and structure can lead to engineering potent materials for mRNA-based gene therapy applications.

Lipid remodeling under acidic conditions and its interplay with low Pi stress in Arabidopsis

Here we show that accumulation of galactose-containing lipids in plastid membranes in shoots and the other membranes in roots maintains Arabidopsis growth under acidic stress and acidic phosphate deficiency. Soil acidification and phosphate deficiency are closely related to each other in natural environments. In addition to the toxicity of high proton concentrations, acid soil can lead to imbalances of ion availability and nutritional deficiencies, including inorganic phosphate (Pi). Among plants, activation of non-phosphorus-containing galactolipid, digalactosyldiacylglycerol (DGDG), synthesis concomitant with phospholipid degradation, namely membrane lipid remodeling, is crucial for coping with Pi starvation. However, regulation mechanisms of membrane lipid composition during acidic stress have not been clarified. Here, we investigated lipid metabolism in Arabidopsis thaliana grown under acidic stress with or without Pi. Under Pi-sufficient acidic conditions, DGDG was increased in shoot membranes, and some Pi starvation-responsive genes that are involved in lipid remodeling were upregulated without reducing Pi content in leaves. In contrast, under acidic Pi deficiency, membrane lipid remodeling in roots was partially repressed at a lower external pH. Nevertheless, phenotypic comparison between wild type and the double mutant of MGD2/3, which are responsible for DGDG accumulation during Pi starvation, indicated that the complete absence of lipid remodeling in roots resulted in a loss of tolerance to Pi deficiency rather specifically under acidic conditions. This result suggested important physiological roles of galactolipid-enriched membranes under acidic Pi deficiency.

Structural characterization of cationic DODAB bilayers containing C24:1 β-glucosylceramide

The effect of 5 mol%, 9 mol%, and 16 mol% of C24:1 β-glucosylceramide (βGlcCer) on the structure of cationic DODAB bilayers was investigated by means of differential scanning calorimetry (DSC), electron spin resonance (ESR) spectroscopy and fluorescence microscopy. βGlcCer is completely miscible with DODAB at all fractions tested, since no domains were observed in fluorescence microscopy or ESR spectra. The latter showed that βGlcCer destabilized the gel phase of DODAB bilayers by decreasing the gel phase packing. As a consequence, βGlcCer induced a decrease in the phase transition temperature and cooperativity of DODAB bilayers, as seen in DSC thermograms. ESR spectra also showed that βGlcCer induced an increase in DODAB fluid phase order and/or rigidity. Despite their different structures, a similar effect of loosening the gel phase packing and turning the fluid phase more rigid/organized has also been observed when low molar fractions of cholesterol were incorporated in DODAB bilayers. The structural characterization of mixed membranes made of cationic lipids and glucosylceramides may be important for developing novel immunotherapeutic tools such as vaccine adjuvants.

Exploring Dynamics and Structure of Biomolecules, Cryoprotectants, and Water Using Molecular Dynamics Simulations: Implications for Biostabilization and Biopreservation

Successful stabilization and preservation of biological materials often utilize low temperatures and dehydration to arrest molecular motion. Cryoprotectants are routinely employed to help the biological entities survive the physicochemical and mechanical stresses induced by cold or dryness. Molecular interactions between biomolecules, cryoprotectants, and water fundamentally determine the outcomes of preservation. The optimization of assays using the empirical approach is often limited in structural and temporal resolution, whereas classical molecular dynamics simulations can provide a cost-effective glimpse into the atomic-level structure and interaction of individual molecules that dictate macroscopic behavior. Computational research on biomolecules, cryoprotectants, and water has provided invaluable insights into the development of new cryoprotectants and the optimization of preservation methods. We describe the rapidly evolving state of the art of molecular simulations of these complex systems, summarize the molecular-scale protective and stabilizing mechanisms, and discuss the challenges that motivate continued innovation in this field.

Make It Simple: (SR-A1+TLR7) Macrophage Targeted NANOarchaeosomes

Hyperhalophilic archaebacteria exclusively produce sn2,3 diphytanylglycerol diether archaeolipids, unique structures absent in bacteria and eukaryotes. Nanovesicles made of archaeolipids known as nanoarchaeosomes (nanoARC), possess highly stable bilayers, some of them displaying specific targeting ability. Here we hypothesize that nanoARC made from Halorubrum tebenquichense archaebacteria, may constitute efficient carriers for the TLR7 agonist imiquimod (IMQ). NanoARC-IMQ takes advantage of the intense interaction between IMQ and the highly disordered, poorly fluid branched archaeolipid bilayers, rich in archaeol analog of methyl ester of phosphatidylglycerophosphate (PGP-Me), a natural ligand of scavenger receptor A1 (SR-A1). This approach lacks complex manufacture steps required for bilayers labeling, enabling future analytical characterization, batch reproducibility, and adaptation to higher scale production. SR-A1 mediated internalization of particulate material is mostly targeted to macrophages and is extensive because it is not submitted to a negative feedback. A massive and selective intracellular delivery of IMQ may concentrate its effect specifically into the endosomes, where the TLR7 is expressed, magnifying its immunogenicity, at the same time reducing its systemic bioavailability, and therefore it's in vivo adverse effects. NanoARC-IMQ (600-900 nm diameter oligolamellar vesicles of ~-43 mV Z potential) were heavily loaded with IMQ at ~44 μg IMQ/mg phospholipids [~20 folds higher than the non-SR-A1 ligand soyPC liposomes loaded with IMQ (LIPO-IMQ)]. In vitro, nanoARC-IMQ induced higher TNF-α and IL-6 secretion by J774A1 macrophages compared to same dose of IMQ and same lipid dose of LIPO-IMQ. In vivo, 3 subcutaneous doses of nanoARC-IMQ+ 10 μg total leishmania antigens (TLA) at 50 μg IMQ per Balb/C mice, induced more pronounced DTH response, accompanied by a nearly 2 orders higher antigen-specific systemic IgG titers than IMQ+TLA and LIPO-IMQ. The isotype ratio of nanoARC-IMQ+TLA remained ~0.5 indicating, the same as IMQ+TLA, a Th2 biased response distinguished by a pronounced increase in antibody titers, without negative effects on splenocytes lymphoproliferation, with a potential CD8+LT induction 10 days after the last dose. Overall, this first approach showed that highly SR-A1 mediated internalization of heavily loaded nanoARC-IMQ, magnified the effect of IMQ on TLR7 expressing macrophages, leading to a more intense in vivo immune response.

Glucosylceramide modifies the LPS-induced inflammatory response in macrophages and the orientation of the LPS/TLR4 complex in silico

Toll-like receptor 4 (TLR4) is activated by bacterial lipopolysaccharide (LPS), which drives the production of proinflammatory cytokines. Earlier studies have indicated that cholesterol- and glycosphingolipid-rich subregions of the plasma membrane (lipid domains) are important for TLR4-mediated signaling. We report that inhibition of glucosylceramide (GluCer) synthase, which resulted in decreased concentrations of the glycosphingolipid GluCer in lipid domains, reduced the LPS-induced inflammatory response in both mouse and human macrophages. Atomistic molecular dynamics simulations of the TLR4 dimer complex (with and without LPS in its MD-2 binding pockets) in membranes (in the presence and absence of GluCer) showed that: (1) LPS induced a tilted orientation of TLR4 and increased dimer integrity; (2) GluCer did not affect the integrity of the LPS/TLR4 dimer but reduced the LPS-induced tilt; and (3) GluCer increased electrostatic interactions between the membrane and the TLR4 extracellular domain, which could potentially modulate the tilt. We also showed that GCS inhibition reduced the interaction between TLR4 and the intracellular adaptor protein Mal. We conclude that the GluCer-induced effects on LPS/TLR4 orientation may influence the signaling capabilities of the LPS/TLR4 complex by affecting its interaction with downstream signaling proteins.

Membrane property and biofunction of phospholiposome incorporated with anomeric galactolipids

There has been increasing interest in the construction of liposomes containing a targeting reagent for target-specific drug delivery. Glycoconjugates that can be recognized by transmembrane glycoprotein receptors have been extensively used to form glyco-liposomal drug carriers. However, the impact of anomerism, which is a common identity of natural glycoconjugates, on the glyco-liposomal properties has been hardly probed in previous studies. Here we investigate the liposomal properties of phospholipid incorporated with a pair of anomeric galactolipids. The anomeric galacto-liposomes are characterized and their membrane fluidity, thermo-stability, DNA condensation efficiency and fluorescence leakage are comparatively tested. The in vitro cellular internalization effect of the galacto-liposomes is also demonstrated. This study suggests that anomerism might give distinct impact on the membrane properties and even biofunctions of glyco-liposomes.

Computer modelling studies of the bilayer/water interface

This review summarises high resolution studies on the interface of lamellar lipid bilayers composed of the most typical lipid molecules which constitute the lipid matrix of biomembranes. The presented results were obtained predominantly by computer modelling methods. Whenever possible, the results were compared with experimental results obtained for similar systems. The first and main section of the review is concerned with the bilayer-water interface and is divided into four subsections. The first describes the simplest case, where the interface consists only of lipid head groups and water molecules and focuses on interactions between the lipid heads and water molecules; the second describes the interface containing also mono- and divalent ions and concentrates on lipid-ion interactions; the third describes direct inter-lipid interactions. These three subsections are followed by a discussion on the network of direct and indirect inter-lipid interactions at the bilayer interface. The second section summarises recent computer simulation studies on the interactions of antibacterial membrane active compounds with various models of the bacterial outer membrane. This article is part of a Special Issue entitled: Biosimulations edited by Ilpo Vattulainen and Tomasz Róg.

Molecular simulations of glycolipids: Towards mammalian cell membrane models

Glycolipids are key components of mammalian cell membranes, influencing a diverse range of cellular functions. For example, a number of receptor tyrosine kinases, including the epidermal growth factor receptor (EGFR), are allosterically regulated by the glycolipid monosialodihexosylganglioside (GM3). Recent advances in molecular dynamics methods, especially the development of coarse-grained models, have enabled simulations of increasingly complex models of cell membranes. We demonstrate these methodological developments via a case study of a coarse-grained model for the ganglioside GM3. This glycolipid is included in simulations of a mixed lipid bilayer model reflecting the compositional complexity of a mammalian cell membrane. The resultant membrane model is used to simulate the interactions of GM3 with the transmembrane domain of the EGFR.

•

We review recent progress in molecular dynamics simulations of glycolipids.

•

A coarse-grained model of the ganglioside GM3 is described.

•

The GM3 model is used in coarse-grained simulations of a mammalian cell membrane.

•

We describe the interactions of GM3 with the transmembrane domain of the EGFR.

A computer model of a polyunsaturated monogalactolipid bilayer

1,2-di-O-acyl-3-O-β-D-galactopyranosyl-sn-glycerol (MGDG) is the main lipid component of thylakoid membranes of higher plants and algae. This monogalactolipid is thought of as a non-bilayer lipid but actually it can form both lamellar and nonlamellar phases. In this study, molecular dynamics (MD) simulations of the fully hydrated di-18:3 MGDG bilayer in the lamellar phase were carried out at 310 and 295 K for 200 and 450 ns, respectively, using the GROMACS 4 software package and OPLS-AA force field. At both temperatures, the lamellar phase of the systems was stable. The pure di-18:3 MGDG bilayer is the first step towards creating a computer model of the lipid matrix of the thylakoid membrane and the main aim of this study was to validate the computer model of di-18:3 MGDG in the bilayer and also to assess the properties of the bilayer. However, only a few of the predicted properties could be compared with those derived experimentally and in other MD simulations because of insufficient amount of such data. Thus, direct validation of the MGDG bilayer proved difficult. Therefore, in the validation process also an indirect approach was used, in which a computer model of the 1,2-dioleoyl-sn-glycero-3-phosphatidylcholine (DOPC) bilayer simulated at the same temperatures using the same force field as the MGDG bilayer was assessed. Successful validation of the DOPC bilayer parameterized in the OPLS-AA force field and similar properties of the MGDG molecules in the pure 18:3 MGDG and in binary 18:3 MGDG-PC bilayers indicate that the computer model of the MDGD molecule is faithful and the MGDG bilayer is representative on the time scales covered in these MD simulations.

Phase Transition of Glycolipid Membranes Studied by Coarse-Grained Simulations

Glycolipids are important components of biological membranes. High concentrations of glycolipids are particularly found in lipid rafts, which take part in many physiological phenomena. This different partitioning and interaction pattern of glycolipids in the membrane as compared to those of phospholipids are likely due to their different chemical structures: the polar regions of glycosphingolipids can be even larger than for their hydrophobic moieties, giving rise to a rich conformational landscape. Here we study the influence of glycosphingolipids galactosylceramide (GCER) and monosialotetrahexosylganglioside (GM1) on the structural and thermodynamic properties of a phospholipid (DPPC) bilayer. Using the method of coarse-grained molecular dynamics simulation we show that both glycolipids increase the phase-transition temperature of phospholipid membranes and that the extent of this increase depends on the headgroup size and structure. GM1 shows a strong tendency to form mixed clusters with phospholipids, thereby stabilizing the membrane. In contrast, GCER is dispersed in the membrane. By occupying the interstitial space between phospholipids it causes a tighter packing of the lipids in the membrane.

Synthetic biology outside the cell: linking computational tools to cell-free systems

As mathematical models become more commonly integrated into the study of biology, a common language for describing biological processes is manifesting. Many tools have emerged for the simulation of in vivo synthetic biological systems, with only a few examples of prominent work done on predicting the dynamics of cell-free synthetic systems. At the same time, experimental biologists have begun to study dynamics of in vitro systems encapsulated by amphiphilic molecules, opening the door for the development of a new generation of biomimetic systems. In this review, we explore both in vivo and in vitro models of biochemical networks with a special focus on tools that could be applied to the construction of cell-free expression systems. We believe that quantitative studies of complex cellular mechanisms and pathways in synthetic systems can yield important insights into what makes cells different from conventional chemical systems.

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.

Conformational dynamics of dry lamellar crystals of sugar based lipids: an atomistic simulation study

The rational design of a glycolipid application (e.g. drug delivery) with a tailored property depends on the detailed understanding of its structure and dynamics. Because of the complexity of sugar stereochemistry, we have undertaken a simulation study on the conformational dynamics of a set of synthetic glycosides with different sugar groups and chain design, namely dodecyl β-maltoside, dodecyl β-cellobioside, dodecyl β-isomaltoside and a C12C10 branched β-maltoside under anhydrous conditions. We examined the chain structure in detail, including the chain packing, gauche/trans conformations and chain tilting. In addition, we also investigated the rotational dynamics of the headgroup and alkyl chains. Monoalkylated glycosides possess a small amount of gauche conformers (∼20%) in the hydrophobic region of the lamellar crystal (LC) phase. In contrast, the branched chain glycolipid in the fluid Lα phase has a high gauche population of up to ∼40%. Rotational diffusion analysis reveals that the carbons closest to the headgroup have the highest correlation times. Furthermore, its value depends on sugar type, where the rotational dynamics of an isomaltose was found to be 11-15% and more restrained near the sugar, possibly due to the chain disorder and partial inter-digitation compared to the other monoalkylated lipids. Intriguingly, the present simulation demonstrates the chain from the branched glycolipid bilayer has the ability to enter into the hydrophilic region. This interesting feature of the anhydrous glycolipid bilayer simulation appears to arise from a combination of lipid crowding and the amphoteric nature of the sugar headgroups.

Refined OPLS all-atom force field for saturated phosphatidylcholine bilayers at full hydration

We report parametrization of dipalmitoyl-phosphatidylcholine (DPPC) in the framework of the Optimized Parameters for Liquid Simulations all-atom (OPLS-AA) force field. We chose DPPC as it is one of the most studied phospholipid species and thus has plenty of experimental data necessary for model validation, and it is also one of the highly important and abundant lipid types, e.g., in lung surfactant. Overall, PCs have not been previously parametrized in the OPLS-AA force field; thus, there is a need to derive its bonding and nonbonding parameters for both the polar and nonpolar parts of the molecule. In the present study, we determined the parameters for torsion angles in the phosphatidylcholine and glycerol moieties and in the acyl chains, as well the partial atomic charges. In these calculations, we used three methods: (1) Hartree-Fock (HF), (2) second order Møller-Plesset perturbation theory (MP2), and (3) density functional theory (DFT). We also tested the effect of the polar environment by using the polarizable continuum model (PCM), and for acyl chains the van der Waals parameters were also adjusted. In effect, six parameter sets were generated and tested on a DPPC bilayer. Out of these six sets, only one was found to be able to satisfactorily reproduce experimental data for the lipid bilayer. The successful DPPC model was obtained from MP2 calculations in an implicit polar environment (PCM).

Deuterium NMR Investigation of the lyotropic phases of alkyl β-glycoside/D2O systems

We have investigated the phase behavior of four glycosides (βC8OGlc, βC8SGlc, βC10OGlc, βC8OGal) in water and D2O by optical polarizing microscopy and deuterium NMR. Previously published phase diagrams were evaluated by deuterium NMR, via monitoring D2O spectra, and confirmed the presence of the hexagonal, bicontinuous cubic, and lamellar phases in these glycosides. We have also shown the presence of the gel phase in (βC10OGlc) and observed the extensive supercooling of the lamellar phase to temperatures well below the Kraft line. While the main features of the phase diagrams were confirmed, some phase boundaries were found to be slightly different. Magnetically aligned spectra were also observed for relatively dilute samples for the hexagonal phase (βC8OGlc and βC8OGal) and the lamellar phase (βC8SGlc and βC10OGlc). The average number of bound water molecules per headgroup in the lamellar phase for the glycosides was determined by the systematic measurement of the quadrupolar splitting of D2O over a wide range of values of the (glycoside/water) molar ratio. The number of water molecules bound to the headgroup was found on average to be about 1.6-1.7 water molecules with no significant differences in this value for the different glycosides (and over the temperature range investigated), indicating that the bound water content is predominately influenced by the number of hydroxyl groups of the headgroup only. However, this bound water content of only 1.6-1.7 water molecules per sugar headgroup is surprisingly low, suggesting strong intermolecular interactions of the OH groups of headgroup sugars. The results are in line with computational results reported earlier for the octyl-β-glucoside and β-galactoside, which show the presence of strong intralayer hydrogen bonding.

Molecular dynamics simulations of membranes composed of glycolipids and phospholipids

Lipid membranes composed of 1,2-di-(9Z,12Z,15Z)-octadecatrienoyl-3-O-β-D-galactosyl-sn-glycerol or monogalactosyldiacylglycerol (MGDG) and 1,2-dimyristoyl-sn-glycero-3-phosphocholine (DMPC) were studied by means of molecular dynamics (MD) computer simulations. Three lipid compositions were considered: 0%, 20%, and 45% MGDG (by mole) denoted as MG-0, MG-20, and MG-45, respectively. The article is focused on the calculation of NMR dipolar interactions, which were confronted with previously reported experimental couplings. Dynamical processes and orientational distributions relevant for the averaging of dipolar interactions were evaluated. Furthermore, several parameters important for characterization of the bilayer structure, molecular organization, and dynamics were investigated. In general, only a minor change in DMPC properties was observed upon the increased MGDG/DMPC ratio, whereas properties related to MGDG undergo a more pronounced change. This effect was ascribed to the fact that DMPC is a bilayer (L(α)) forming lipid, whereas MGDG prefers a reverse hexagonal (H(II)) arrangement.

Study of PEGylated lipid layers as a model for PEGylated liposome surfaces: molecular dynamics simulation and Langmuir monolayer studies

We have combined Langmuir monolayer film experiments and all-atom molecular dynamics (MD) simulation of a bilayer to study the surface structure of a PEGylated liposome and its interaction with the ionic environment present under physiological conditions. Lipids that form both gel and liquid-crystalline membranes have been used in our study. By varying the salt concentration in the Langmuir film experiment and including salt at the physiological level in the simulation, we have studied the effect of salt ions present in the blood plasma on the structure of the poly(ethylene glycol) (PEG) layer. We have also studied the interaction between the PEG layer and the lipid bilayer in both the liquid-crystalline and gel states. The MD simulation shows two clear results: (a) The Na(+) ions form close interactions with the PEG oxygens, with the PEG chains forming loops around them and (b) PEG penetrates the lipid core of the membrane for the case of a liquid-crystalline membrane but is excluded from the tighter structure of the gel membrane. The Langmuir monolayer results indicate that the salt concentration affects the PEGylated lipid system, and these results can be interpreted in a fashion that is in agreement with the results of our MD simulation. We conclude that the currently accepted picture of the PEG surface layer acting as a generic neutral hydrophilic polymer entirely outside the membrane, with its effect explained through steric interactions, is not sufficient. The phenomena we have observed may affect both the interaction between the liposome and bloodstream proteins and the liquid-crystalline-gel transition and is thus relevant to nanotechnological drug delivery device design.

Structure and Dynamics of Glycosphingolipids in Lipid Bilayers: Insights from Molecular Dynamics Simulations

Glycolipids are important constituents of biological membranes, and understanding their structure and dynamics in lipid bilayers provides insights into their physiological and pathological roles. Experimental techniques have provided details into their behavior at model and biological membranes; however, computer simulations are needed to gain atomic level insights. This paper summarizes the insights obtained from MD simulations into the conformational and orientational dynamics of glycosphingolipids and their exposure, hydration, and hydrogen-bonding interactions in membrane environment. The organization of glycosphingolipids in raft-like membranes and their modulation of lipid membrane structure are also reviewed.

Mitochondrial membranes with mono- and divalent salt: changes induced by salt ions on structure and dynamics

We employ atomistic simulations to consider how mono- (NaCl) and divalent (CaCl(2)) salt affects properties of inner and outer membranes of mitochondria. We find that the influence of salt on structural properties is rather minute, only weakly affecting lipid packing, conformational ordering, and membrane electrostatic potential. The changes induced by salt are more prominent in dynamical properties related to ion binding and formation of ion-lipid complexes and lipid aggregates, as rotational diffusion of lipids is slowed down by ions, especially in the case of CaCl(2). In the same spirit, lateral diffusion of lipids is slowed down rather considerably for increasing concentration of CaCl(2). Both findings for dynamic properties can be traced to the binding of ions with lipid head groups and the related changes in interaction patterns in the headgroup region, where the binding of Na(+) and Ca(2+) ions is clearly different. The role of cardiolipins in these phenomena turns out to be important.