Modification of Lipid Bilayer Structure by Diacylglycerol: A Comparative Study of Diacylglycerol and Cholesterol

A Membrane-Assisted Mechanism for the Release of Ceramide from the CERT START Domain

Ceramide transfer protein CERT is the mediator of nonvesicular transfer of ceramide from the ER to Golgi. In CERT, START is the domain responsible for the binding and transport of ceramide. A wealth of structural data has revealed a helix-grip fold surrounding a large hydrophobic cavity holding the ceramide. Yet, little is known about the mechanisms by which START releases the ceramide through the polar region and into the packed environment of cellular membranes. As such events do not lend themselves easily to experimental investigations, we used multiple unbiased microsecond-long molecular simulations. We propose a membrane-assisted mechanism in which the membrane acts as an allosteric effector initiating the release of ceramide and where the passage of the ceramide acyl chains is facilitated by the intercalation of a single phosphatidylcholine lipid in the cavity, practically greasing the ceramide way out. We verify using free energy calculation and experimental lipidomics data that CERT forms stable complexes with phosphatidylcholine lipids, in addition to ceramide, thus providing validation for the proposed mechanism.

Phase-Separated Lipid-Based Nanoparticles: Selective Behavior at the Nano-Bio Interface

The membrane-protein interface on lipid-based nanoparticles influences their in vivo behavior. Better understanding may evolve current drug delivery methods toward effective targeted nanomedicine. Previously, the cell-selective accumulation of a liposome formulation in vivo is demonstrated, through the recognition of lipid phase-separation by triglyceride lipases. This exemplified how liposome morphology and composition can determine nanoparticle-protein interactions. Here, the lipase-induced compositional and morphological changes of phase-separated liposomes-which bear a lipid droplet in their bilayer- are investigated, and the mechanism upon which lipases recognize and bind to the particles is unravelled. The selective lipolytic degradation of the phase-separated lipid droplet is observed, while nanoparticle integrity remains intact. Next, the Tryptophan-rich loop of the lipase is identified as the region with which the enzymes bind to the particles. This preferential binding is due to lipid packing defects induced on the liposome surface by phase separation. In parallel, the existing knowledge that phase separation leads to in vivo selectivity, is utilized to generate phase-separated mRNA-LNPs that target cell-subsets in zebrafish embryos, with subsequent mRNA delivery and protein expression. Together, these findings can expand the current knowledge on selective nanoparticle-protein communications and in vivo behavior, aspects that will assist to gain control of lipid-based nanoparticles.

Insertion Depth Modulates Protein Kinase C-δ-C1b Domain Interactions with Membrane Cholesterol as Revealed by MD Simulations

Protein kinase C delta (PKC-δ) is an important signaling molecule in human cells that has both proapoptotic as well as antiapoptotic functions. These conflicting activities can be modulated by two classes of ligands, phorbol esters and bryostatins. Phorbol esters are known tumor promoters, while bryostatins have anti-cancer properties. This is despite both ligands binding to the C1b domain of PKC-δ (δC1b) with a similar affinity. The molecular mechanism behind this discrepancy in cellular effects remains unknown. Here, we have used molecular dynamics simulations to investigate the structure and intermolecular interactions of these ligands bound to δC1b with heterogeneous membranes. We observed clear interactions between the δC1b-phorbol complex and membrane cholesterol, primarily through the backbone amide of L250 and through the K256 side-chain amine. In contrast, the δC1b-bryostatin complex did not exhibit interactions with cholesterol. Topological maps of the membrane insertion depth of the δC1b-ligand complexes suggest that insertion depth can modulate δC1b interactions with cholesterol. The lack of cholesterol interactions suggests that bryostatin-bound δC1b may not readily translocate to cholesterol-rich domains within the plasma membrane, which could significantly alter the substrate specificity of PKC-δ compared to δC1b-phorbol complexes.

Phase-Separated Liposomes Hijack Endogenous Lipoprotein Transport and Metabolism Pathways to Target Subsets of Endothelial Cells In Vivo

Plasma lipid transport and metabolism are essential to ensure correct cellular function throughout the body. Dynamically regulated in time and space, the well-characterized mechanisms underpinning plasma lipid transport and metabolism offers an enticing, but as yet underexplored, rationale to design synthetic lipid nanoparticles with inherent cell/tissue selectivity. Herein, a systemically administered liposome formulation, composed of just two lipids, that is capable of hijacking a triglyceride lipase-mediated lipid transport pathway resulting in liposome recognition and uptake within specific endothelial cell subsets is described. In the absence of targeting ligands, liposome-lipase interactions are mediated by a unique, phase-separated ("parachute") liposome morphology. Within the embryonic zebrafish, selective liposome accumulation is observed at the developing blood-brain barrier. In mice, extensive liposome accumulation within the liver and spleen - which is reduced, but not eliminated, following small molecule lipase inhibition - supports a role for endothelial lipase but highlights these liposomes are also subject to significant "off-target" by reticuloendothelial system organs. Overall, these compositionally simplistic liposomes offer new insights into the discovery and design of lipid-based nanoparticles that can exploit endogenous lipid transport and metabolism pathways to achieve cell selective targeting in vivo.

Insights into the behavior of unsaturated diacylglycerols in mixed lipid bilayers in relation to protein kinase C activation-A molecular dynamics simulation study

The lipid second messenger diacylglycerol (DAG) is known for its involvement in many types of cellular signaling, especially as an endogenous agonist for protein kinase C (PKC). Evidence has emerged that the degree of saturation of the DAG molecules can affect PKC activation. DAG molecules with different acyl chain saturation have not only been observed to induce varying extents of PKC activation, but also to express selectivity towards different PKC isozymes. Both qualities are important for precise therapeutic activation of PKC; understanding DAG behavior at the molecular level in different environments has much potential in the development of drugs to target PKC. We used molecular dynamics simulations to study the behavior of two different unsaturated DAG species in lipid environments with varying degrees of unsaturation. We focus on phosphatidylethanolamine (PE) instead of phosphatidylcholine (PC) to more accurately model the relevant biomembranes. The effect of cholesterol (CHOL) on these systems was also explored. We found that both high level of unsaturation in the acyl chains of the DAG species and presence of CHOL in the surrounding membrane increase DAG molecule availability at the lipid-water interface. This can partially explain the previously observed differences in PKC activation strength and specificity, the complete mechanism is, however, likely to be more complex. Our simulations coupled with the current understanding of lipids highlight the need for more simulations of biologically accurate lipid environments in order to determine the correct correlations between molecular mechanisms and biological behavior when studying PKC activation.

Interaction between glycolipid MPIase and proteinaceous factors during protein integration into the cytoplasmic membrane of E. coli

Protein integration into biomembranes is an essential biological phenomenon common to all organisms. While various factors involved in protein integration, such as SRP, SecYEG and YidC, are proteinaceous, we identified a glycolipid named MPIase (Membrane Protein Integrase), which is present in the cytoplasmic membrane of E. coli. In vitro experiments using inverted membrane vesicles prepared from MPIase-depleted strains, and liposomes containing MPIase showed that MPIase is required for insertion of a subset of membrane proteins, which has been thought to be SecYEG-independent and YidC-dependent. Also, SecYEG-dependent substrate membrane proteins require MPIase in addition. Furthermore, MPIase is also essential for insertion of proteins with multiple negative charges, which requires both YidC and the proton motive force (PMF). MPIase directly interacts with SecYEG and YidC on the membrane. MPIase not only cooperates with these factors but also has a molecular chaperone-like function specific to the substrate membrane proteins through direct interaction with the glycan chain. Thus, MPIase catalyzes membrane insertion by accepting nascent membrane proteins on the membrane through its chaperone-like function, i.e., direct interaction with the substrate proteins, and then MPIase functionally interacts with SecYEG and YidC for substrate delivery, and acts with PMF to facilitate and complete membrane insertion when necessary. In this review, we will outline the mechanisms underlying membrane insertion catalyzed by MPIase, which cooperates with proteinaceous factors and PMF.

Biogenesis and Breakdown of Lipid Droplets in Pathological Conditions

Lipid droplets (LD) have long been considered as mere fat drops; however, LD have lately been revealed to be ubiquitous, dynamic and to be present in diverse organelles in which they have a wide range of key functions. Although incompletely understood, the biogenesis of eukaryotic LD initiates with the synthesis of neutral lipids (NL) by enzymes located in the endoplasmic reticulum (ER). The accumulation of NL leads to their segregation into nanometric nuclei which then grow into lenses between the ER leaflets as they are further filled with NL. The lipid composition and interfacial tensions of both ER and the lenses modulate their shape which, together with specific ER proteins, determine the proneness of LD to bud from the ER toward the cytoplasm. The most important function of LD is the buffering of energy. But far beyond this, LD are actively integrated into physiological processes, such as lipid metabolism, control of protein homeostasis, sequestration of toxic lipid metabolic intermediates, protection from stress, and proliferation of tumours. Besides, LD may serve as platforms for pathogen replication and defense. To accomplish these functions, from biogenesis to breakdown, eukaryotic LD have developed mechanisms to travel within the cytoplasm and to establish contact with other organelles. When nutrient deprivation occurs, LD undergo breakdown (lipolysis), which begins with the LD-associated members of the perilipins family PLIN2 and PLIN3 chaperone-mediated autophagy degradation (CMA), a specific type of autophagy that selectively degrades a subset of cytosolic proteins in lysosomes. Indeed, PLINs CMA degradation is a prerequisite for further true lipolysis, which occurs via cytosolic lipases or by lysosome luminal lipases when autophagosomes engulf portions of LD and target them to lysosomes. LD play a crucial role in several pathophysiological processes. Increased accumulation of LD in non-adipose cells is commonly observed in numerous infectious diseases caused by intracellular pathogens including viral, bacterial, and parasite infections, and is gradually recognized as a prominent characteristic in a variety of cancers. This review discusses current evidence related to the modulation of LD biogenesis and breakdown caused by intracellular pathogens and cancer.

Rotational decoupling between the hydrophilic and hydrophobic regions in lipid membranes

Cells use homeostatic mechanisms to ensure an optimal composition of distinct types of lipids in cellular membranes. The hydrophilic region of biological lipid membranes is mainly composed of several types of phospholipid headgroups that interact with incoming molecules, nanoparticles, and viruses, whereas the hydrophobic region consists of a distribution of acyl chains and sterols affecting membrane fluidity/rigidity related properties and forming an environment for membrane-bound molecules such as transmembrane proteins. A fundamental open question is to what extent the motions of these regions are coupled and, consequently, how strongly the interactions of phospholipid headgroups with other molecules depend on the properties and composition of the membrane hydrophobic core. We combine advanced solid-state nuclear magnetic resonance spectroscopy with high-fidelity molecular dynamics simulations to demonstrate how the rotational dynamics of choline headgroups remain nearly unchanged (slightly faster) with incorporation of cholesterol into a phospholipid membrane, contrasting the well-known extreme slowdown of the other phospholipid segments. Notably, our results suggest a new paradigm in which phospholipid dipole headgroups interact as quasi-freely rotating flexible dipoles at the interface, independent of the properties in the hydrophobic region.

Sterols are required for the coordinated assembly of lipid droplets in developing seeds

Lipid droplets (LDs) are intracellular organelles critical for energy storage and lipid metabolism. They are typically composed of an oil core coated by a monolayer of phospholipids and proteins such as oleosins. The mechanistic details of LD biogenesis remain poorly defined. However, emerging evidence suggest that their formation is a spatiotemporally regulated process, occurring at specific sites of the endoplasmic reticulum defined by a specific set of lipids and proteins. Here, we show that sterols are required for formation of oleosin-coated LDs in Arabidopsis. Analysis of sterol pathway mutants revealed that deficiency in several ∆⁵-sterols accounts for the phenotype. Importantly, mutants deficient in these sterols also display reduced LD number, increased LD size and reduced oil content in seeds. Collectively, our data reveal a role of sterols in coordinating the synthesis of oil and oleosins and their assembly into LDs, highlighting the importance of membrane lipids in regulating LD biogenesis.

Lipid droplet biogenesis originates at the endoplasmic reticulum and is defined by a specific set of lipids and proteins. Here, the authors show that sterols play an important role in coordinating oil and oleosin biosynthesis for the formation of lipid droplets in plant leaves and seeds.

Functional Group Distributions, Partition Coefficients, and Resistance Factors in Lipid Bilayers Using Site Identification by Ligand Competitive Saturation

Small molecules such as metabolites and drugs must pass through the membrane of the cell, a barrier primarily comprising phospholipid bilayers and embedded proteins. To better understand the process of passive diffusion, knowledge of the ability of various functional groups to partition across bilayers and the associated energetics would be of utility. In the present study, the site identification by ligand competitive saturation (SILCS) methodology has been applied to sample the distributions of a diverse set of chemical solutes representing the functional groups of small molecules across phospholipid bilayers composed of 0.9:0.1 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphocholine/cholesterol and a mixture of 0.52:0.18:0.3 1,2-dioleoyl-sn-glycero-3-phospho-l-serine/1,2-dioleoyl-sn-glycero-3-phosphocholine/cholesterol used in parallel artificial membrane permeability assay experiments. A combination of oscillating chemical potential grand canonical Monte Carlo and molecular dynamics in the SILCS simulations was applied to achieve solute sampling through the bilayers and surrounding aqueous environment from which the distribution of solutes and the functional groups they represent were obtained. Results show differential distribution of aliphatic versus aromatic groups with the former having increased sampling in the center of the bilayers versus in the region of the glycerol linker for the latter. Variations in the distribution of different polar groups are evident, with large differences between negative acetate and positive methylammonium with accumulation of the polar-neutral and acetate solutes above the bilayer head groups. Conversion of the distributions to absolute free energies allows for a detailed understanding of energetics of functional groups in different regions of the bilayers and for calculation of absolute free-energy profiles of multifunctional drug-like molecules across the bilayers from which partition coefficients and resistance factors suitable for insertion into the homogenous solubility-diffusion equation for calculation of permeability were obtained. Comparisons of the calculated bilayer/solution partition coefficients with 1-octanol/water experimental data for both drug-like molecules and the solutes show overall good agreement, validating the calculated distributions and associated absolute free-energy profiles.

A metabolic reaction-diffusion model for PKCα translocation via PIP2 hydrolysis in an endothelial cell

Hydrolysis of the phospholipid phosphatidylinositol 4,5-bisphosphate (PIP2) at the cell membrane induces the release of inositol 1,4,5-trisphosphate (IP3) into the cytoplasm and diffusion of diacylglycerol (DAG) through the membrane, respectively. Release of IP3 subsequently increases Ca2+ levels in the cytoplasm, which results in activation of protein kinase C α (PKCα) by Ca2+ and DAG, and finally the translocation of PKCα from the cytoplasm to the membrane. In this study, we developed a metabolic reaction-diffusion framework to simulate PKCα translocation via PIP2 hydrolysis in an endothelial cell. A three-dimensional cell model, divided into membrane and cytoplasm domains, was reconstructed from confocal microscopy images. The associated metabolic reactions were divided into their corresponding domain; PIP2 hydrolysis at the membrane domain resulted in DAG diffusion at the membrane domain and IP3 release into the cytoplasm domain. In the cytoplasm domain, Ca2+ was released from the endoplasmic reticulum, and IP3, Ca2+, and PKCα diffused through the cytoplasm. PKCα bound Ca2+ at, and diffused through, the cytoplasm, and was finally activated by binding with DAG at the membrane. Using our model, we analyzed IP3 and DAG dynamics, Ca2+ waves, and PKCα translocation in response to a microscopic stimulus. We found a qualitative agreement between our simulation results and our experimental results obtained by live-cell imaging. Interestingly, our results suggest that PKCα translocation is dominated by DAG dynamics. This three-dimensional reaction-diffusion mathematical framework could be used to investigate the link between PKCα activation in a cell and cell function.

Structural Transformation in Vesicles upon Hydrolysis of Phosphatidylethanolamine and Phosphatidylcholine with Phospholipase C

This study provides insights into dynamic nanostructural changes in phospholipid systems during hydrolysis with phospholipase C, the fate of the hydrolysis products, and the kinetics of lipolysis. The effect of lipid restructuring of the vesicle was investigated using small-angle X-ray scattering and cryogenic scanning electron microscopy. The rate and extent of phospholipid hydrolysis were quantified using nuclear magnetic resonance. Hydrolysis of two phospholipids, phosphatidylethanolamine (PE) and phosphatidylcholine (PC), results in the cleavage of the molecular headgroup, causing two strikingly different changes in lipid self-assembly. The diacylglycerol product of PC escapes the lipid bilayer, whereas the diacylglycerol product adopts a different configuration within the lipid bilayer of the PE vesicles. These results are then discussed concerning the change of the lipid configuration upon the lipid membrane and its potential implications in vivo, which is of significant importance for the detailed understanding of the fate of lipidic particles and the rational design of enzyme-responsive lipid-based drug delivery systems.

Alteration of Membrane Physicochemical Properties by Two Factors for Membrane Protein Integration

After a nascent chain of a membrane protein emerges from the ribosomal tunnel, the protein is integrated into the cell membrane. This process is controlled by a series of proteinaceous molecular devices, such as signal recognition particles and Sec translocons. In addition to these proteins, we discovered two endogenous components regulating membrane protein integration in the inner membrane of Escherichia coli. The integration is blocked by diacylglycerol (DAG), whereas the blocking is relieved by a glycolipid named membrane protein integrase (MPIase). Here, we investigated the influence of these integration-blocking and integration-promoting factors on the physicochemical properties of membrane lipids via solid-state NMR and fluorescence measurements. These factors did not have destructive effects on membrane morphology because the membrane maintained its lamellar structure and did not fuse in the presence of DAG and/or MPIase at their effective concentrations. We next focused on membrane flexibility. DAG did not affect the mobility of the membrane surface, whereas the sugar chain in MPIase was highly mobile and enhanced the flexibility of membrane lipid headgroups. Comparison with a synthetic MPIase analog revealed the effects of the long sugar chain on membrane properties. The acyl chain order inside the membrane was increased by DAG, whereas the increase was cancelled by the addition of MPIase. MPIase also loosened the membrane lipid packing. Focusing on the transbilayer movement, MPIase reduced the rapid flip-flop motion of DAG. On the other hand, MPIase could not compensate for the diminished lateral diffusion by DAG. These results suggest that by manipulating the membrane lipids dynamics, DAG inhibits the protein from contacting the inner membrane, whereas the flexible long sugar chain of MPIase increases the opportunity for interaction between the membrane and the protein, leading to membrane integration of the newly formed protein.

Phospholipids: Pulling Back the Actin Curtain for Granule Delivery to the Immune Synapse

Phosphoinositides, together with the phospholipids phosphatidylserine and phosphatidic acid, are important components of the plasma membrane acting as second messengers that, with diacylglycerol, regulate a diverse range of signaling events converting extracellular changes into cellular responses. Local changes in their distribution and membrane charge on the inner leaflet of the plasma membrane play important roles in immune cell function. Here we discuss their distribution and regulators highlighting the importance of membrane changes across the immune synapse on the cytoskeleton and the impact on the function of cytotoxic T lymphocytes.

Fate of Liposomes in the Presence of Phospholipase C and D: From Atomic to Supramolecular Lipid Arrangement

Understanding the origins of lipid membrane bilayer rearrangement in response to external stimuli is an essential component of cell biology and the bottom-up design of liposomes for biomedical applications. The enzymes phospholipase C and D (PLC and PLD) both cleave the phosphorus-oxygen bonds of phosphate esters in phosphatidylcholine (PC) lipids. The atomic position of this hydrolysis reaction has huge implications for the stability of PC-containing self-assembled structures, such as the cell wall and lipid-based vesicle drug delivery vectors. While PLC converts PC to diacylglycerol (DAG), the interaction of PC with PLD produces phosphatidic acid (PA). Here we present a combination of small-angle scattering data and all-atom molecular dynamics simulations, providing insights into the effects of atomic-scale reorganization on the supramolecular assembly of PC membrane bilayers upon enzyme-mediated incorporation of DAG or PA. We observed that PC liposomes completely disintegrate in the presence of PLC, as conversion of PC to DAG progresses. At lower concentrations, DAG molecules within fluid PC bilayers form hydrogen bonds with backbone carbonyl oxygens in neighboring PC molecules and burrow into the hydrophobic region. This leads initially to membrane thinning followed by a swelling of the lamellar phase with increased DAG. At higher DAG concentrations, localized membrane tension causes a change in lipid phase from lamellar to the hexagonal and micellar cubic phases. Molecular dynamics simulations show that this destabilization is also caused in part by the decreased ability of DAG-containing PC membranes to coordinate sodium ions. Conversely, PLD-treated PC liposomes remain stable up to extremely high conversions to PA. Here, the negatively charged PA headgroup attracts significant amounts of sodium ions from the bulk solution to the membrane surface, leading to a swelling of the coordinated water layer. These findings are a vital step toward a fundamental understanding of the degradation behavior of PC lipid membranes in the presence of these clinically relevant enzymes, and toward the rational design of diagnostic and drug delivery technologies for phospholipase-dysregulation-based diseases.

Construction of Epidermal Growth Factor Receptor Peptide Magnetic Nanovesicles with Lipid Bilayers for Enhanced Capture of Liver Cancer Circulating Tumor Cells

Highly effective targeted tumor recognition via vectors is crucial for cancer detection. In contrast to antibodies and proteins, peptides are direct targeting ligands with a low molecular weight. In the present study, a peptide magnetic nanovector platform containing a lipid bilayer was designed using a peptide amphiphile (PA) as a skeleton material in a controlled manner without surface modification. Fluorescein isothiocyanate-labeled epidermal growth factor receptor (EGFR) peptide nanoparticles (NPs) could specifically bind to EGFR-positive liver tumor cells. EGFR peptide magnetic vesicles (EPMVs) could efficiently recognize and separate hepatoma carcinoma cells from cell solutions and treated blood samples (ratio of magnetic EPMVs versus anti-EpCAM NPs: 3.5 ± 0.29). Analysis of the circulating tumor cell (CTC) count in blood samples from 32 patients with liver cancer showed that EPMVs could be effectively applied for CTC capture. Thus, this nanoscale, targeted cargo-packaging technology may be useful for designing cancer diagnostic systems.

Distribution and dynamics of quinones in the lipid bilayer mimicking the inner membrane of mitochondria

Quinone and its analogues (Q) constitute an important class of compounds that perform key electron transfer reactions in oxidative- and photo-phosphorylation. In the inner membrane of mitochondria, ubiquinone molecules undergo continuous redox transitions enabling electron transfer between the respiratory complexes. In such a dynamic system undergoing continuous turnover for ATP synthesis, an uninterrupted supply of substrate molecules is absolutely necessary. In the current work, we have performed atomistic molecular dynamics simulations and free energy calculations to assess the structure, dynamics, and localization of quinone and its analogues in a lipid bilayer, whose composition mimics the one in the inner mitochondrial membrane. The results show that there is a strong tendency of both quinone and quinol molecules to localize in the vicinity of the lipids' acyl groups, right under the lipid head group region. Additionally, we observe a second location in the middle of the bilayer where quinone molecules tend to stabilize. Translocation of quinone through a lipid bilayer is very fast and occurs in 10-100ns time scale, whereas the translocation of quinol is at least an order of magnitude slower. We suggest that this has important mechanistic implications given that the localization of Q ensures maximal occupancy of the Q-binding sites or Q-entry points in electron transport chain complexes, thereby maintaining an optimal turnover rate for ATP synthesis.

Toward Atomistic Resolution Structure of Phosphatidylcholine Headgroup and Glycerol Backbone at Different Ambient Conditions

Phospholipids are essential building blocks of biological membranes. Despite a vast amount of very accurate experimental data, the atomistic resolution structures sampled by the glycerol backbone and choline headgroup in phoshatidylcholine bilayers are not known. Atomistic resolution molecular dynamics simulations have the potential to resolve the structures, and to give an arrestingly intuitive interpretation of the experimental data, but only if the simulations reproduce the data within experimental accuracy. In the present work, we simulated phosphatidylcholine (PC) lipid bilayers with 13 different atomistic models, and compared simulations with NMR experiments in terms of the highly structurally sensitive C-H bond vector order parameters. Focusing on the glycerol backbone and choline headgroups, we showed that the order parameter comparison can be used to judge the atomistic resolution structural accuracy of the models. Accurate models, in turn, allow molecular dynamics simulations to be used as an interpretation tool that translates these NMR data into a dynamic three-dimensional representation of biomolecules in biologically relevant conditions. In addition to lipid bilayers in fully hydrated conditions, we reviewed previous experimental data for dehydrated bilayers and cholesterol-containing bilayers, and interpreted them with simulations. Although none of the existing models reached experimental accuracy, by critically comparing them we were able to distill relevant chemical information: (1) increase of choline order parameters indicates the P-N vector tilting more parallel to the membrane, and (2) cholesterol induces only minor changes to the PC (glycerol backbone) structure. This work has been done as a fully open collaboration, using nmrlipids.blogspot.fi as a communication platform; all the scientific contributions were made publicly on this blog. During the open research process, the repository holding our simulation trajectories and files ( https://zenodo.org/collection/user-nmrlipids ) has become the most extensive publicly available collection of molecular dynamics simulation trajectories of lipid bilayers.

Alteration of lipid membrane structure and dynamics by diacylglycerols with unsaturated chains

Diacylglycerols (DAGs) with unsaturated acyl chains play many important roles in biomembranes, such as a second messenger and activator for protein kinase C. In this study, three DAGs of distinctly different chain unsaturations (i.e. di16:0DAG (DPG), 16:0-18:1DAG (POG), and di18:1DAG (DOG)) are studied using atomistic MD simulation to compare their roles in the structure and dynamics of 16:0-18:1phosphatidylcholine (POPC) membranes. All three DAGs are able to produce the so-called 'condensing effect' in POPC membranes: decreasing area-per-lipid, and increasing acyl chain order and bilayer thickness. Our visual and quantitative analyses clearly show that DAG with unsaturated chains induce larger spacing between POPC headgroups, compared with DAG with saturated chains; this particular effect has long been hypothesized to be crucial for activating enzymes and receptors in cell membranes. DAGs with unsaturated chains are also located closer to the bilayer/aqueous interface than DPG and are more effective in slowing down lateral diffusion of molecules. We show that DAG molecules seek the "umbrella coverage" from neighboring phospholipid headgroups - similar to cholesterol. Unlike cholesterol, DAGs also hide their chains from water by laterally inserting their chains into the surrounding. Thus, acyl chains of DAG are more spread and disordered than those of PC due to the insertion. By calculating the potential of mean force (PMF) for POPC in POPC/DAG bilayers, we found that all three DAGs can significantly increase the free energy barrier for POPC to flip-flop, but only DAGs with unsaturated chains can additionally increase the free energy of POPC desorption.

Aluminum ions alter the function of non-specific phospholipase C through the changes in plasma membrane physical properties

The first indication of the aluminum (Al) toxicity in plants growing in acidic soils is the cessation of root growth, but the detailed mechanism of Al effect is unknown. Here we examined the impact of Al stress on the activity of non-specific phospholipase C (NPC) in the connection with the processes related to the plasma membrane using fluorescently labeled phosphatidylcholine. We observed a rapid and significant decrease of labeled diacylglycerol (DAG), product of NPC activity, in Arabidopsis seedlings treated with AlCl₃. Interestingly, an application of the membrane fluidizer, benzyl alcohol, restored the level of DAG during Al treatment. Our observations suggest that the activity of NPC is affected by Al-induced changes in plasma membrane physical properties.

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.

Cholesterol modulates function of connexin 43 gap junction channel via PKC pathway in H9c2 cells

It has been shown that cholesterol modulates activity of protein kinase C (PKC), and PKC phosphorylates connexin 43 (Cx43) to regulate its function, respectively. However, it is not known whether cholesterol modulates function of Cx43 through regulating activity of PKC. In the present study, we demonstrated that cholesterol enrichment reduced the dye transfer ability of Cx43 in cultured H9c2 cells. Western blot analysis indicated that cholesterol enrichment enhanced the phosphorylated state of Cx43. Immunofluorescent images showed that cholesterol enrichment made the Cx43 distribution from condensed to diffused manner in the interface between the cells. In cholesterol enriched cells, PKC antagonists partially restored the dye transfer ability among the cells, downregulated the phosphorylation of Cx43 and redistributed Cx43 from the diffused manner to the condensed manner in the cell interface. In addition, reduction of cholesterol level suppressed PKC activity to phosphorylate Cx43 and restored Cx43 function in PKC agonist-treated cells. Furthermore, we demonstrated that cholesterol enrichment upregulated the phosphorylated state of Cx43 at Ser368, while PKC antagonists reversed the effect. Taken together, cholesterol level in the cells plays important roles in regulating Cx43 function through activation of the PKC signaling pathway.

Complex roles of hybrid lipids in the composition, order, and size of lipid membrane domains

Hybrid lipids (HL) are phospholipids with one saturated chain and one unsaturated chain. HL are hypothesized to act as linactants (i.e., 2D surfactants) in cell membranes, reducing line tension and creating nanoscopic lipid domains. Here we compare three hybrid lipids of different chain unsaturation (16:0-18:1PC (POPC), 16:0-18:2PC (PLPC), and 16:0-20:4PC (PAPC)) in their abilities to alter the composition, line tension, order, and compactness of lipid domains. We found that the liquid-ordered (Lo) and liquid-disordered (Ld) lipid domains in PAPC/di18:0PC(DSPC)/cholesterol and PLPC/DSPC/cholesterol mixtures are micrometer-sized, and only the POPC/DSPC/cholesterol system has nanoscopic domains. The results indicate that some HLs with polyunsaturated chains are not linactants, and the monounsaturated POPC displays both properties of weak linactants and "Ld-phase" lipids such as di18:1PC (DOPC). The obtained phase boundaries from giant unilamellar vesicles (GUV) show that both POPC and PLPC partition well in the Lo phases. Our MD simulations reveal that these hybrid lipids decrease the order and compactness of Lo domains. Thus, hybrid lipids distinguish themselves from other lipid groups in this combined "partitioning and loosening" ability, which could explain why the Lo domains of GUVs, which often do not contain HL, are more compact than the raft domains in cell membranes. Our line tension measurement and Monte Carlo simulation both show that even the monounsaturated POPC is a weak linactant with only modest ability to occupy domain boundaries and reduce line tension. A more important property of HLs is that they can reduce physical property differences of Lo and Ld bulk domains, which also reduces line tension at domain boundaries.

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.

Conical lipids in flat bilayers induce packing defects similar to that induced by positive curvature

In biological membranes, changes in lipid composition or mechanical deformations produce defects in the geometrical arrangement of lipids, thus allowing the adsorption of certain peripheral proteins. Here, we perform molecular dynamics simulations on bilayers containing a cylindrical lipid (PC) and a conical lipid (DOG). Profiles of atomic density and lateral pressure across the bilayer show differences in the acyl chain region due to deeper partitioning of DOG compared to PC. However, such analyses are less informative for the interfacial region where peripheral proteins adsorb. To circumvent this limitation, we develop, to our knowledge, a new method of membrane surface analysis. This method allows the identification of chemical defects, where hydrocarbon chains are accessible to the solvent, and geometrical defects, i.e., voids deeper than the glycerol backbone. The size and number of both types of defects increase with the number of monounsaturated acyl chains in PC and with the introduction of DOG, although the defects do not colocalize with the conical lipid. Interestingly, the size and probability of the defects promoted by DOG resemble those induced by positive curvature, thus explaining why conical lipids and positive curvature can both drive the adsorption of peripheral proteins that use hydrophobic residues as membrane anchors.

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

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

Molecular properties of various structurally defined sphingomyelins -- correlation of structure with function

Sphingomyelins are important phospholipids in plasma membranes of most cells. Because of their dominantly saturated nature, they affect the lateral structure of membranes, and contribute to the regulation of cholesterol distribution within membranes, and in cells. However, the abundance of molecular species present in cells also implies that sphingomyelins have other, more specific functions. Many of these functions are currently unknown, but are under extensive study. Mostly model membrane studies have shown that sphingomyelins (and other sphingolipids), in contrast to glycerophospholipids, have important hydrogen bonding properties which in several important ways confer specific functional properties to this abundant class of membrane phospholipids. The often very asymmetric nature of sphingomyelins, arising from mismatch in length between the long chain base and N-acyl chains, also impose specific properties (e.g., interdigitation) to sphingomyelins not seen with glycerophospholipids. In this review, the latest sphingomyelin literature will be scrutinized, and an effort will be made to correlate the molecular structure of sphingomyelin with functional properties. In particular, the effects of head group properties, interfacial hydrogen bonding, long chain base hydroxylation, N-acyl chain hydroxylation, and N-acyl chain methyl-branching will be discussed.

Cholesterol and POPC segmental order parameters in lipid membranes: solid state 1H-13C NMR and MD simulation studies

The concentration of cholesterol in cell membranes affects membrane fluidity and thickness, and might regulate different processes such as the formation of lipid rafts. Since interpreting experimental data from biological membranes is rather intricate, investigations on simple models with biological relevance are necessary to understand the natural systems. We study the effect of cholesterol on the molecular structure of multi-lamellar vesicles (MLVs) composed of 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphocholine (POPC), a phospholipid ubiquitous in cell membranes, with compositions in the range 0-60 mol% cholesterol. Order parameters, |S(CH)|, are experimentally determined by using (1)H-(13)C solid-state nuclear magnetic resonance (NMR) spectroscopy with segmental detail for all parts of both the cholesterol and POPC molecules, namely the ring system and alkyl chain of the sterol, as well as the glycerol backbone, choline headgroup and the sn-1 and sn-2 acyl chains of POPC. With increasing cholesterol concentration the acyl chains gradually adopt a more extended conformation while the orientation and dynamics of the polar groups are rather unaffected. Additionally, we perform classical molecular dynamics simulations on virtual bilayers mimicking the POPC-cholesterol MLVs investigated by NMR. Good agreement between experiments and simulations is found for the cholesterol alignment in the bilayer and for the |S(CH)| profiles of acyl chains below 15 mol% cholesterol. Deviations occur for the choline headgroup and glycerol backbone parts of POPC, as well as for the phospholipid and cholesterol alkyl chains at higher cholesterol concentrations. The unprecedented detail of the NMR data enables a more complete comparison between simulations and experiments on POPC-cholesterol bilayers and may aid in developing more realistic model descriptions of biological membranes.