Local accumulation of diacylglycerol alters membrane properties nonlinearly due to its transbilayer activity

Elucidating of the metabolic impact of risperidone on brain microvascular endothelial cells using untargeted metabolomics-based LC-MS

Risperidone is useful for the treatment of schizophrenia symptoms; however, it also has side effects, and an overdose can be harmful. The metabolic effects of risperidone at high therapeutic doses and its metabolites have not been elucidated. Endogenous cellular metabolites may be comprehensively analyzed using untargeted metabolomics-based liquid chromatography-mass spectrometry (LC-MS), which can reveal changes in cell regulation and metabolic pathways. By identifying the metabolites and pathway changes using a nontargeted metabolomics-based LC-MS approach, we aimed to shed light on the potential toxicological effects of high-dose risperidone on brain microvascular endothelial cells (MVECs) associated with the human blood brain barrier. A total of 42 metabolites were selected as significant putative metabolites of the toxicological response of high-dose risperidone in MVECs. Six highly correlated pathways were identified, including those involving diacylglycerol, fatty acid, ceramide, glycerophospholipid, amino acid, and tricarboxylic acid metabolism. We demonstrated that methods focused on metabolomics are useful for identifying metabolites that may be used to clarify the mechanism of drug-induced toxicity.

Structure-function relationship of phase-separated liposomes containing diacylglycerol analogues

The composition and morphology of lipid-based nanoparticles can influence their overall in vivo behavior. Previously, we demonstrated that phase separation in liposomes composed of DSPC and a diacylglycerol lipid analogue (DOaG) drives the in vivo biodistribution towards a specific subset of endothelial cells in zebrafish embryos. In the absence of traditional targeting functionalities (e.g., antibodies, ligands), this selectivity is mediated solely by the unique liposome morphology and composition, characterized by a DOaG-rich lipid droplet within the DSPC-rich phospholipid bilayer. The phase separation is induced due to the geometry of DOaG lipid and its ability to create non-bilayer phases in lipid membranes. To investigate the underlying principles of phase separation and to optimize the liposome colloidal stability, we performed a structure-function relationship study by synthesizing a library of DOaG analogues with varying molecular properties, such as the number, length and sn-position of the acyl chains, as well as the degree of saturation or carbonyl substituents. We assessed the ability of these lipid analogues to assemble into phase-separated liposomes and studied their morphology, colloidal stability, and in vivo biodistribution in zebrafish embryos. We found that analogues containing unsaturated, medium length (C16-C18) fatty acids were required to obtain colloidally stable, phase-separated liposomes with cell-specific biodistribution patterns. Moreover, we observed that using the pure DOaG isomer, with acyl chains at the sn-1,3 positions, leads to more colloidally stable liposomes than when a mixture of sn-1,2 and sn-1,3 isomers is used. Similarly, we observed that incorporating a DOaG analogue with fatty tails shorter than DSPC, as well as PEGylation, endows liposomes with long term stability while retaining cell-selective biodistribution. Diacylglycerols are known to promote fusion, lipid polymorphism, signaling and protein recruitment on lipid membranes. In this study, we showed that diacylglycerol derivatives can induce phase separation in liposomes, unlocking the potential for cell-specific targeting in vivo. We believe that these findings can be the foundation for future use of diacylglycerols in lipid-based nanomedicines and could lead to the development of novel targeted delivery strategies.

Molecular dynamics simulations of intracellular lipid droplets: a new tool in the toolbox

Lipid droplets (LDs) are ubiquitous intracellular organelles with a central role in multiple lipid metabolic pathways. However, identifying correlations between their structural properties and their biological activity has proved challenging, owing to their unique physicochemical properties as compared with other cellular membranes. In recent years, molecular dynamics (MD) simulations, a computational methodology allowing the accurate description of molecular assemblies down to their individual components, have been demonstrated to be a useful and powerful approach for studying LD structural and dynamical properties. In this short review, we attempt to highlight, as comprehensively as possible, how MD simulations have contributed to our current understanding of multiple molecular mechanisms involved in LD biology.

Analyzing Lipid Membrane Defects via a Coarse-Grained to Triangulated Surface Map: The Role of Lipid Order and Local Curvature in Molecular Binding

Lipid packing defects are known to serve as quantitative indicators for protein binding to lipid membranes. This paper presents a protocol for mapping molecular lipid detail onto a triangulated continuum leaflet representation. Besides establishing the desired forward counterpart to the existing inverse TS2CG map, this coarse-grained to triangulated surface (CG2TS) map enables straightforward extraction of the defect characteristics for any membrane geometry found in nature. We have applied our protocol to investigate the role of local curvature and varying lipid packing on the defect constant π. We find that the defect size is greatly influenced by both factors, arguing strongly against the usual assignment of a single defect constant in the case of more realistic membrane conditions. An important discovery is that lipids in the gel phase produce larger defects, or a higher π, in domains of high (local) curvature than the same lipid in a liquid phase of any curvature. This finding suggests that membranes featuring very ordered lipid packing can bind proteins via large defects in curved regions. Finally, we propose a route for estimating defect constants directly from the standard membrane properties. Identifying the precise role of composition, lipid (tail) order, and (local) curvature in defects for the irregular lipid structures that are (temporally) present in many biological processes is instrumental for obtaining fundamental insight as well as for a rational design of membrane binding targets.

Triglycerides Stabilize Water/Organic Interfaces of Changing Area via Conformational Flexibility

The role of triglycerides (TGs) in both natural and synthetic biological membranes has long been the subject of study, involving metabolism, disease, and colloidal synthesis. TGs have been found to be critical components for successful liposomal encapsulation via a water/oil/water double emulsion, which this work endeavors to explain. TGs can occupy multiple positions in biological membranes. The glycerol backbone can reside at the water/organic interface, adjacent to phospholipid headgroups ("m" conformation), typically with relatively low (<3%) solubility. The glycerol backbone can also occupy hydrophobic regions, where it is isolated from water ("h" or "oil" conformation). This can occur in either midmembrane positions or phospholipid-coated lipid droplets (LDs). These conformations can be distinguished using 13C-nuclear magnetic resonance spectroscopy (NMR), which determines the degree of hydration of the TG backbone. Using this method, it was revealed that TGs transition from "m" to "h" conformation as the organic solvent is removed via evaporation. A new transitional TG backbone position has been identified with a level of hydration between "m" and "h". These results suggest that TGs can temporarily coat and stabilize the large water/organic interfaces present after emulsification. As the organic solvent is removed and interfaces shrink, the TGs recede into midmembrane spaces or bud off into LDs, which are confirmed via transmission electron microscopy (TEM) and can be removed via centrifugation. Encapsulation efficiency is found to be inversely related to both the saturation and length of the TG acyl chains, indicating that membrane fluidization is a key property arising from the presence of TGs. Beyond clarification of a mechanism for high-efficiency liposomal encapsulation, these results implicate TGs as components that are able to stabilize biological membrane transitions involving a changing interfacial area and curvature. This role for TGs may be of use in the formulation of drug delivery systems as well as in the investigation of membrane transitions in life sciences.

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.

Catalytic core function of yeast Pah1 phosphatidate phosphatase reveals structural insight into its membrane localization and activity control

The PAH1-encoded phosphatidate (PA) phosphatase is a major source of diacylglycerol for the production of the storage lipid triacylglycerol and a key regulator for the de novo phospholipid synthesis in Saccharomyces cerevisiae. The catalytic function of Pah1 depends on its membrane localization which is mediated through its phosphorylation by multiple protein kinases and dephosphorylation by the Nem1-Spo7 protein phosphatase complex. The full-length Pah1 is composed of a catalytic core (N-LIP and HAD-like domains, amphipathic helix, and the WRDPLVDID domain) and non-catalytic regulatory sequences (intrinsically disordered regions, RP domain, and acidic tail) for phosphorylation and interaction with Nem1-Spo7. How the catalytic core regulates Pah1 localization and cellular function is not clear. In this work, we analyzed a variant of Pah1 (i.e., Pah1-CC (catalytic core)) that is composed only of the catalytic core. Pah1-CC expressed on a low-copy plasmid complemented the pah1Δ mutant phenotypes (e.g., nuclear/ER membrane expansion, reduced levels of triacylglycerol, and lipid droplet formation) without requiring Nem1-Spo7. The cellular function of Pah1-CC was supported by its PA phosphatase activity mostly associated with the membrane fraction. Although functional, Pah1-CC was distinct from Pah1 in the protein and enzymological properties, which include overexpression toxicity, association with heat shock proteins, and significant reduction of the Vmax value. These findings on the Pah1 catalytic core enhance the understanding of its structural requirements for membrane localization and activity control.

Membrane-Active Thermoresponsive Block Copolymers Containing a Diacylglycerol-Based Segment: RAFT Synthesis, Doxorubicin Encapsulation, and Evaluation of Cytotoxicity against Breast Cancer Cells

Herein, we report the formation of drug delivery systems from original thermoresponsive block copolymers containing lipid-based segments. Two acrylate monomers derived from palmitic- or oleic-acid-based diacylglycerols (DAGs) were synthesized and polymerized by the reversible addition-fragmentation chain transfer (RAFT) method. Well-defined DAG-based polymers with targeted molar masses and narrow molar mass distributions were next used as macro-chain transfer agents (macro-CTAs) for the polymerization of N-isopropylacrylamide (NIPAAm) or N-vinylcaprolactam (NVCL). The obtained amphiphilic block copolymers were formed into polymeric nanoparticles (PNPs) with and without encapsulated doxorubicin and characterized. Their biological assessment indicated appropriate cytocompatibility with the representatives of normal cells. Furthermore, compared to the free drug, increased cytotoxicity and apoptosis or necrosis induction in breast cancer cells was documented, including a highly aggressive and invasive triple-negative MDA-MB-231 cell line.

Diacylglycerol-dependent hexamers of the SNARE-assembling chaperone Munc13-1 cooperatively bind vesicles

Munc13-1 is essential for vesicle docking and fusion at the active zone of synapses. Here, we report that Munc13-1 self-assembles into molecular clusters within diacylglycerol-rich microdomains present in phospholipid bilayers. Although the copy number of Munc13-1 molecules in these clusters has a broad distribution, a systematic Poisson analysis shows that this is most likely the result of two molecular species: monomers and mainly hexameric oligomers. Each oligomer is able to capture one vesicle independently. Hexamers have also been observed in crystals of Munc13-1 that form between opposed phospholipid bilayers [K. Grushin, R. V. Kalyana Sundaram, C. V. Sindelar, J. E. Rothman, Proc. Natl. Acad. Sci. U.S.A. 119, e2121259119 (2022)]. Mutations targeting the contacts stabilizing the crystallographic hexagons also disrupt the isolated hexamers, suggesting they are identical. Additionally, these mutations also convert vesicle binding from a cooperative to progressive mode. Our study provides an independent approach showing that Munc13-1 can form mainly hexamers on lipid bilayers each capable of vesicle capture.

Pragmatic Coarse-Graining of Proteins: Models and Applications

The molecular details involved in the folding, dynamics, organization, and interaction of proteins with other molecules are often difficult to assess by experimental techniques. Consequently, computational models play an ever-increasing role in the field. However, biological processes involving large-scale protein assemblies or long time scale dynamics are still computationally expensive to study in atomistic detail. For these applications, employing coarse-grained (CG) modeling approaches has become a key strategy. In this Review, we provide an overview of what we call pragmatic CG protein models, which are strategies combining, at least in part, a physics-based implementation and a top-down experimental approach to their parametrization. In particular, we focus on CG models in which most protein residues are represented by at least two beads, allowing these models to retain some degree of chemical specificity. A description of the main modern pragmatic protein CG models is provided, including a review of the most recent applications and an outlook on future perspectives in the field.

Untargeted Lipidomic Approach for Studying Different Nervous System Tissues of the Murine Model of Krabbe Disease

Krabbe disease is a rare neurodegenerative disease with an autosomal recessive character caused by a mutation in the GALC gene. The mutation leads to an accumulation of psychosine and a subsequent degeneration of oligodendrocytes and Schwann cells. Psychosine is the main biomarker of the disease. The Twitcher mouse is the most commonly used animal model to study Krabbe disease. Although there are many references to this model in the literature, the lipidomic study of nervous system tissues in the Twitcher model has received little attention. This study focuses on the comparison of the lipid profiles of four nervous system tissues (brain, cerebellum, spinal cord, and sciatic nerve) in the Twitcher mouse compared to the wild-type mouse. Altogether, approximately 230 molecular species belonging to 19 lipid classes were annotated and quantified. A comparison at the levels of class, molecular species, and lipid building blocks showed significant differences between the two groups, particularly in the sciatic nerve. The in-depth study of the lipid phenotype made it possible to hypothesize the genes and enzymes involved in the changes. The integration of metabolic data with genetic data may be useful from a systems biology perspective to gain a better understanding of the molecular basis of the disease.

Novel Roles for Diacylglycerol in Synaptic Vesicle Priming and Release Revealed by Complete Reconstitution of Core Protein Machinery

Here we introduce the full functional reconstitution of genetically-validated core protein machinery (SNAREs, Munc13, Munc18, Synaptotagmin, Complexin) for synaptic vesicle priming and release in a geometry that enables detailed characterization of the fate of docked vesicles both before and after release is triggered with Ca 2+ . Using this novel setup, we discover new roles for diacylglycerol (DAG) in regulating vesicle priming and Ca 2+- triggered release involving the SNARE assembly chaperone Munc13. We find that low concentrations of DAG profoundly accelerate the rate of Ca 2+ -dependent release, and high concentrations reduce clamping and permit extensive spontaneous release. As expected, DAG also increases the number of ready-release vesicles. Dynamic single-molecule imaging of Complexin binding to ready-release vesicles directly establishes that DAG accelerates the rate of SNAREpin assembly mediated by Munc13 and Munc18 chaperones. The selective effects of physiologically validated mutations confirmed that the Munc18-Syntaxin-VAMP2 'template' complex is a functional intermediate in the production of primed, ready-release vesicles, which requires the coordinated action of Munc13 and Munc18.

SIGNIFICANCE STATEMENT

Munc13 and Munc18 are SNARE-associated chaperones that act as "priming" factors, facilitating the formation of a pool of docked, release-ready vesicles and regulating Ca 2+ -evoked neurotransmitter release. Although important insights into Munc18/Munc13 function have been gained, how they assemble and operate together remains enigmatic. To address this, we developed a novel biochemically-defined fusion assay which enabled us to investigate the cooperative action of Munc13 and Munc18 in molecular terms. We find that Munc18 nucleates the SNARE complex, while Munc13 promotes and accelerates the SNARE assembly in a DAG-dependent manner. The concerted action of Munc13 and Munc18 stages the SNARE assembly process to ensure efficient 'clamping' and formation of stably docked vesicles, which can be triggered to fuse rapidly (∼10 msec) upon Ca 2+ influx.

Comparison between Sickle Cell Disease Patients and Healthy Donors: Untargeted Lipidomic Study of Erythrocytes

Sickle cell disease (SCD) is one of the most common severe monogenic disorders in the world caused by a mutation on HBB gene and characterized by hemoglobin polymerization, erythrocyte rigidity, vaso-occlusion, chronic anemia, hemolysis, and vasculopathy. Recently, the scientific community has focused on the multiple genetic and clinical profiles of SCD. However, the lipid composition of sickle cells has received little attention in the literature. According to recent studies, changes in the lipid profile are strongly linked to several disorders. Therefore, the aim of this study is to dig deeper into lipidomic analysis of erythrocytes in order to highlight any variations between healthy and patient subjects. 241 lipid molecular species divided into 17 classes have been annotated and quantified. Lipidomic profiling of SCD patients showed that over 24% of total lipids were altered most of which are phospholipids. In-depth study of significant changes in lipid metabolism can give an indication of the enzymes and genes involved. In a systems biology scenario, these variations can be useful to improve the understanding of the biochemical basis of SCD and to try to make a score system that could be predictive for the severity of clinical manifestations.

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.

Triglyceride lipolysis triggers liquid crystalline phases in lipid droplets and alters the LD proteome

Lipid droplets (LDs) are reservoirs for triglycerides (TGs) and sterol-esters (SEs), but how these lipids are organized within LDs and influence their proteome remain unclear. Using in situ cryo-electron tomography, we show that glucose restriction triggers lipid phase transitions within LDs generating liquid crystalline lattices inside them. Mechanistically this requires TG lipolysis, which decreases the LD's TG:SE ratio, promoting SE transition to a liquid crystalline phase. Molecular dynamics simulations reveal TG depletion promotes spontaneous TG and SE demixing in LDs, additionally altering the lipid packing of the PL monolayer surface. Fluorescence imaging and proteomics further reveal that liquid crystalline phases are associated with selective remodeling of the LD proteome. Some canonical LD proteins, including Erg6, relocalize to the ER network, whereas others remain LD-associated. Model peptide LiveDrop also redistributes from LDs to the ER, suggesting liquid crystalline phases influence ER-LD interorganelle transport. Our data suggests glucose restriction drives TG mobilization, which alters the phase properties of LD lipids and selectively remodels the LD proteome.

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.

Reprogrammed Lipid Metabolism and the Lipid-Associated Hallmarks of Colorectal Cancer

Simple Summary

Colorectal cancer (CRC) is the third-most diagnosed cancer and the second-leading cause of cancer-related deaths worldwide. Limitations in early and accurate diagnosis of CRC gives rise to poor patient survival. Advancements in analytical techniques have improved our understanding of the cellular and metabolic changes occurring in CRC and potentiate avenues for improved diagnostic and therapeutic strategies. Lipids are metabolites with important biological functions; however, their role in CRC is poorly understood. Here, we provide an in-depth review of the recent literature concerning lipid alterations in CRC and propose eight lipid metabolism-associated hallmarks of CRC.

Abstract

Lipids have diverse structures, with multifarious regulatory functions in membrane homeostasis and bioenergetic metabolism, in mediating functional protein–lipid and protein–protein interactions, as in cell signalling and proliferation. An increasing body of evidence supports the notion that aberrant lipid metabolism involving remodelling of cellular membrane structure and changes in energy homeostasis and signalling within cancer-associated pathways play a pivotal role in the onset, progression, and maintenance of colorectal cancer (CRC) and their tumorigenic properties. Recent advances in analytical lipidome analysis technologies have enabled the comprehensive identification and structural characterization of lipids and, consequently, our understanding of the role they play in tumour progression. However, despite progress in our understanding of cancer cell metabolism and lipidomics, the key lipid-associated changes in CRC have yet not been explicitly associated with the well-established ‘hallmarks of cancer’ defined by Hanahan and Weinberg. In this review, we summarize recent findings that highlight the role of reprogrammed lipid metabolism in CRC and use this growing body of evidence to propose eight lipid metabolism-associated hallmarks of colorectal cancer, and to emphasize their importance and linkages to the established cancer hallmarks.

Acyl-CoA:wax alcohol acyltransferase 2 modulates the cone visual cycle in mouse retina

The daylight and color vision of diurnal vertebrates depends on cone photoreceptors. The capability of cones to operate and respond to changes in light brightness even under high illumination is attributed to their fast rate of recovery to the ground photosensitive state. This process requires the rapid replenishing of photoisomerized visual chromophore (11-cis-retinal) to regenerate cone visual pigments. Recently, several gene candidates have been proposed to contribute to the cone-specific retinoid metabolism, including acyl-CoA wax alcohol acyltransferase 2 (AWAT2, aka MFAT). Here, we evaluated the role of AWAT2 in the regeneration of visual chromophore by the phenotypic characterization of Awat2^-/- mice. The global absence of AWAT2 enzymatic activity did not affect gross retinal morphology or the rate of visual chromophore regeneration by the canonical RPE65-dependent visual cycle. Analysis of Awat2 expression indicated the presence of the enzyme throughout the murine retina, including the retinal pigment epithelium (RPE) and Müller cells. Electrophysiological recordings revealed reduced maximal rod and cone dark-adapted responses in AWAT2-deficient mice compared to control mice. While rod dark adaptation was not affected by the lack of AWAT2, M-cone dark adaptation both in isolated retina and in vivo was significantly suppressed. Altogether, these results indicate that while AWAT2 is not required for the normal operation of the canonical visual cycle, it is a functional component of the cone-specific visual chromophore regenerative pathway.

Computational Studies of Lipid Droplets

Lipid droplets (LDs) are intracellular organelles whose primary function is energy storage. Known to emerge from the endoplasmic reticulum (ER) bilayer, LDs have a unique structure with a core consisting of neutral lipids, triacylglycerol (TG) or sterol esters (SE), surrounded by a phospholipid (PL) monolayer and decorated by proteins that come and go throughout their complex lifecycle. In this Feature Article, we review recent developments in computational studies of LDs, a rapidly growing area of research. We highlight how molecular dynamics (MD) simulations have provided valuable molecular-level insight into LD targeting and LD biogenesis. Additionally, we review the physical properties of TG from different force fields compared with experimental data. Possible future directions and challenges are discussed.

The Bak core dimer focuses triacylglycerides in the membrane

Apoptosis, the intrinsic programmed cell death process, is mediated by the Bcl-2 family members Bak and Bax. Activation via formation of symmetric core dimers and oligomerization on the mitochondrial outer membrane (MOM) leads to permeabilization and cell death. Although this process is linked to the MOM, the role of the membrane in facilitating such pores is poorly understood. We recently described Bak core domain dimers, revealing lipid binding sites and an initial role of lipids in oligomerization. Here we describe simulations that identified localized clustering and interaction of triacylglycerides (TAGs) with a minimized Bak dimer construct. Coalescence of TAGs occurred beneath this Bak dimer, mitigating dimer-induced local membrane thinning and curvature in representative coarse-grain MOM and model membrane systems. Furthermore, the effects observed as a result of coarse-grain TAG cluster formation was concentration dependent, scaling from low physiological MOM concentrations to those found in other organelles. We find that increasing the TAG concentration in liposomes mimicking the MOM decreased the ability of activated Bak to permeabilize these liposomes. These results suggest that the presence of TAGs within a Bak-lipid membrane preserves membrane integrity and is associated with reduced membrane stress, suggesting a possible role of TAGs in Bak-mediated apoptosis.

Recharging your fats: CHARMM36 parameters for neutral lipids triacylglycerol and diacylglycerol

Neutral lipids (NLs) are an abundant class of cellular lipids. They are characterized by the total lack of charged chemical groups in their structure, and, as a consequence, they play a major role in intracellular lipid storage. NLs that carry a glycerol backbone, such as triacylglycerols (TGs) and diacylglycerols (DGs), are also involved in the biosynthetic pathway of cellular phospholipids, and they have recently been the subject of numerous structural investigations by means of atomistic molecular dynamics simulations. However, conflicting results on the physicochemical behavior of NLs were observed depending on the nature of the atomistic force field used. Here, we show that current phospholipid-derived CHARMM36 parameters for DGs and TGs cannot adequately reproduce interfacial properties of these NLs because of excessive hydrophilicity at the glycerol-ester region. By following a CHARMM36-consistent parameterization strategy, we develop improved parameters for both TGs and DGs that are compatible with both cutoff-based and particle mesh Ewald schemes for the treatment of Lennard-Jones interactions. We show that our improved parameters can reproduce interfacial properties of NLs and their behavior in more complex lipid assemblies. We discuss the implications of our findings in the context of intracellular lipid storage and NLs’ cellular activity.

Molecular-dynamics-simulation-guided membrane engineering allows the increase of membrane fatty acid chain length in Saccharomyces cerevisiae

The use of lignocellulosic-based fermentation media will be a necessary part of the transition to a circular bio-economy. These media contain many inhibitors to microbial growth, including acetic acid. Under industrially relevant conditions, acetic acid enters the cell predominantly through passive diffusion across the plasma membrane. The lipid composition of the membrane determines the rate of uptake of acetic acid, and thicker, more rigid membranes impede passive diffusion. We hypothesized that the elongation of glycerophospholipid fatty acids would lead to thicker and more rigid membranes, reducing the influx of acetic acid. Molecular dynamics simulations were used to predict the changes in membrane properties. Heterologous expression of Arabidopsis thaliana genes fatty acid elongase 1 (FAE1) and glycerol-3-phosphate acyltransferase 5 (GPAT5) increased the average fatty acid chain length. However, this did not lead to a reduction in the net uptake rate of acetic acid. Despite successful strain engineering, the net uptake rate of acetic acid did not decrease. We suggest that changes in the relative abundance of certain membrane lipid headgroups could mitigate the effect of longer fatty acid chains, resulting in a higher net uptake rate of acetic acid.

Synaptotagmins at the endoplasmic reticulum-plasma membrane contact sites maintain diacylglycerol homeostasis during abiotic stress

Endoplasmic reticulum-plasma membrane contact sites (ER-PM CS) play fundamental roles in all eukaryotic cells. Arabidopsis thaliana mutants lacking the ER-PM protein tether synaptotagmin1 (SYT1) exhibit decreased PM integrity under multiple abiotic stresses, such as freezing, high salt, osmotic stress, and mechanical damage. Here, we show that, together with SYT1, the stress-induced SYT3 is an ER-PM tether that also functions in maintaining PM integrity. The ER-PM CS localization of SYT1 and SYT3 is dependent on PM phosphatidylinositol-4-phosphate and is regulated by abiotic stress. Lipidomic analysis revealed that cold stress increased the accumulation of diacylglycerol at the PM in a syt1/3 double mutant relative to wild-type while the levels of most glycerolipid species remain unchanged. In addition, the SYT1-green fluorescent protein fusion preferentially binds diacylglycerol in vivo with little affinity for polar glycerolipids. Our work uncovers a SYT-dependent mechanism of stress adaptation counteracting the detrimental accumulation of diacylglycerol at the PM produced during episodes of abiotic stress.

Physical Characterization of Triolein and Implications for Its Role in Lipid Droplet Biogenesis

Lipid droplets (LDs) are neutral lipid-storing organelles surrounded by a phospholipid (PL) monolayer. At present, how LDs are formed in the endoplasmic reticulum (ER) bilayer is poorly understood. In this study, we present a revised all-atom (AA) triolein (TG) model, the main constituent of the LD core, and characterize its properties in a bilayer membrane to demonstrate the implications of its behavior in LD biogenesis. In bilayer simulations, TG resides at the surface, adopting PL-like conformations (denoted in this work as SURF-TG). Free energy sampling simulation results estimate the barrier for TG relocating from the bilayer surface to the bilayer center to be ∼2 kcal/mol in the absence of an oil lens. SURF-TG is able to modulate membrane properties by increasing PL ordering, decreasing bending modulus, and creating local negative curvature. The other neutral lipid, dioleoyl-glycerol (DAG), also reduces the membrane bending modulus and populates negative curvature regions. A phenomenological coarse-grained (CG) model is also developed to observe larger-scale SURF-TG-mediated membrane deformation. CG simulations confirm that TG nucleates between the bilayer leaflets at a critical concentration when SURF-TG is evenly distributed. However, when one monolayer contains more SURF-TG, the membrane bends toward the other leaflet, followed by TG nucleation if a concentration is higher than the critical threshold. The central conclusion of this study is that SURF-TG is a negative curvature inducer, as well as a membrane modulator. To this end, a model is proposed in which the accumulation of SURF-TG in the luminal leaflet bends the ER bilayer toward the cytosolic side, followed by TG nucleation.

Seipin accumulates and traps diacylglycerols and triglycerides in its ring-like structure

Lipid droplets (LDs) are intracellular organelles responsible for lipid storage, and they emerge from the endoplasmic reticulum (ER) upon the accumulation of neutral lipids, mostly triglycerides (TG), between the two leaflets of the ER membrane. LD biogenesis takes place at ER sites that are marked by the protein seipin, which subsequently recruits additional proteins to catalyze LD formation. Deletion of seipin, however, does not abolish LD biogenesis, and its precise role in controlling LD assembly remains unclear. Here, we use molecular dynamics simulations to investigate the molecular mechanism through which seipin promotes LD formation. We find that seipin clusters TG, as well as its precursor diacylglycerol, inside its unconventional ring-like oligomeric structure and that both its luminal and transmembrane regions contribute to this process. This mechanism is abolished upon mutations of polar residues involved in protein-TG interactions into hydrophobic residues. Our results suggest that seipin remodels the membrane of specific ER sites to prime them for LD biogenesis.

Lipidomic and in-gel analysis of maleic acid co-polymer nanodiscs reveals differences in composition of solubilized membranes

Membrane proteins are key in a large number of physiological and pathological processes. Their study often involves a prior detergent solubilization step, which strips away the membrane and can jeopardize membrane protein integrity. A recent alternative to detergents encompasses maleic acid based copolymers (xMAs), which disrupt the lipid bilayer and form lipid protein nanodiscs (xMALPs) soluble in aqueous buffer. Although xMALPs are often referred to as native nanodiscs, little is known about the resemblance of their lipid and protein content to the native bilayer. Here we have analyzed prokaryotic and eukaryotic xMALPs using lipidomics and in-gel analysis. Our results show that the xMALPs content varies with the chemical properties of the used xMA.

Pre-existing bilayer stresses modulate triglyceride accumulation in the ER versus lipid droplets

Cells store energy in the form of neutral lipids (NLs) packaged into micrometer-sized organelles named lipid droplets (LDs). These structures emerge from the endoplasmic reticulum (ER) at sites marked by the protein seipin, but the mechanisms regulating their biogenesis remain poorly understood. Using a combination of molecular simulations, yeast genetics, and fluorescence microscopy, we show that interactions between lipids' acyl-chains modulate the propensity of NLs to be stored in LDs, in turn preventing or promoting their accumulation in the ER membrane. Our data suggest that diacylglycerol, which is enriched at sites of LD formation, promotes the packaging of NLs into LDs, together with ER-abundant lipids, such as phosphatidylethanolamine. On the opposite end, short and saturated acyl-chains antagonize fat storage in LDs and promote accumulation of NLs in the ER. Our results provide a new conceptual understanding of LD biogenesis in the context of ER homeostasis and function.

The Surface and Hydration Properties of Lipid Droplets

Lipid droplets (LDs) are energy storage organelles composed of neutral lipids, such as triacylglycerol (TG) and sterol esters, surrounded by a phospholipid (PL) monolayer. Their central role in metabolism, complex life cycle, and unique lipid monolayer surface have garnered great attention over the last decade. In this article, results from the largest and longest all-atom simulations to date suggest that 5-8% of the LD surface is occupied by TG molecules, a number that exceeds the maximal solubility reported for TGs in PL bilayers (2.8%). Two distinct classes of TG molecules that interact with the LD monolayer are found. Those at the monolayer surface (SURF-TG) are ordered like PLs with the glycerol moiety exposed to water, creating a significant amount of chemically unique packing defects, and the acyl chains extended toward the LD center. In contrast, the TGs that intercalate just into the PL tail region (CORE-TG) are disordered and increase the amount of PL packing defects and the PL tail order. The degree of interdigitation caused by CORE-TG is stable and determines the width of the TG-PL overlap, whereas that caused by SURF-TG fluctuates and is highly correlated with the area per PL or the expansion of the monolayer. Thus, when the supply of PLs is limited, SURF-TG may reduce surface tension by behaving as a secondary membrane component. The hydration properties of the simulated LD systems demonstrate ∼10 times more water in the LD core than previously reported. Collectively, the reported surface and hydration properties are expected to play a direct role in the mechanisms by which proteins target and interact with LDs.

Excess diacylglycerol at the endoplasmic reticulum disrupts endomembrane homeostasis and autophagy

Background

When stressed, eukaryotic cells produce triacylglycerol (TAG) to store nutrients and mobilize autophagy to combat internal damage. We and others previously reported that in yeast, elimination of TAG synthesizing enzymes inhibits autophagy under nitrogen starvation, yet the underlying mechanism has remained elusive.

Results

Here, we show that disruption of TAG synthesis led to diacylglycerol (DAG) accumulation and its relocation from the vacuolar membrane to the endoplasmic reticulum (ER). We further show that, beyond autophagy, ER-accumulated DAG caused severe defects in the endomembrane system, including disturbing the balance of ER-Golgi protein trafficking, manifesting in bulging of ER and loss of the Golgi apparatus. Genetic or chemical manipulations that increase consumption or decrease supply of DAG reversed these defects. In contrast, increased amounts of precursors of glycerolipid synthesis, including phosphatidic acid and free fatty acids, did not replicate the effects of excess DAG. We also provide evidence that the observed endomembrane defects do not rely on Golgi-produced DAG, Pkc1 signaling, or the unfolded protein response.

Conclusions

This work identifies DAG as the critical lipid molecule responsible for autophagy inhibition under condition of defective TAG synthesis and demonstrates the disruption of ER and Golgi function by excess DAG as the potential cause of the autophagy defect.

Computer simulations of protein-membrane systems

The interactions between proteins and membranes play critical roles in signal transduction, cell motility, and transport, and they are involved in many types of diseases. Molecular dynamics (MD) simulations have greatly contributed to our understanding of protein-membrane interactions, promoted by a dramatic development of MD-related software, increasingly accurate force fields, and available computer power. In this chapter, we present available methods for studying protein-membrane systems with MD simulations, including an overview about the various all-atom and coarse-grained force fields for lipids, and useful software for membrane simulation setup and analysis. A large set of case studies is discussed.

The conical shape of DIM lipids promotes Mycobacterium tuberculosis infection of macrophages

Phthiocerol dimycocerosate (DIM) is a major virulence factor of the pathogen Mycobacterium tuberculosis (Mtb). While this lipid promotes the entry of Mtb into macrophages, which occurs via phagocytosis, its molecular mechanism of action is unknown. Here, we combined biophysical, cell biology, and modeling approaches to reveal the molecular mechanism of DIM action on macrophage membranes leading to the first step of Mtb infection. Matrix-assisted laser desorption ionization time-of-flight (MALDI-TOF) mass spectrometry showed that DIM molecules are transferred from the Mtb envelope to macrophage membranes during infection. Multiscale molecular modeling and ³¹P-NMR experiments revealed that DIM adopts a conical shape in membranes and aggregates in the stalks formed between 2 opposing lipid bilayers. Infection of macrophages pretreated with lipids of various shapes uncovered a general role for conical lipids in promoting phagocytosis. Taken together, these results reveal how the molecular shape of a mycobacterial lipid can modulate the biological response of macrophages.