Stressed Lipid Droplets: How Neutral Lipids Relieve Surface Tension and Membrane Expansion Drives Protein Association

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.

The lipid droplet lipidome

Lipid droplets (LDs) are intracellular organelles with a hydrophobic core formed by neutral lipids surrounded by a phospholipid monolayer harboring a variety of regulatory and enzymatically active proteins. Over the last few decades, our understanding of LD biology has evolved significantly. Nowadays, LDs are appreciated not just as passive energy storage units, but rather as active players in the regulation of lipid metabolism and quality control machineries. To fulfill their functions in controlling cellular metabolic states, LDs need to be highly dynamic and responsive organelles. A large body of evidence supports a dynamic nature of the LD proteome and its contact sites with other organelles. However, much less is known about the lipidome of LDs. Numerous examples clearly indicate the intrinsic link between LD lipids and proteins, calling for a deeper characterization of the LD lipidome in various physiological and pathological settings. Here, we reviewed the current state of knowledge in the field of the LD lipidome, providing a brief overview of the lipid classes and their molecular species present within the neutral core and phospholipid monolayer.

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.

Tld1 is a regulator of triglyceride lipolysis that demarcates a lipid droplet subpopulation

Cells store lipids in the form of triglyceride (TG) and sterol ester (SE) in lipid droplets (LDs). Distinct pools of LDs exist, but a pervasive question is how proteins localize to and convey functions to LD subsets. Here, we show that the yeast protein YDR275W/Tld1 (for TG-associated LD protein 1) localizes to a subset of TG-containing LDs and reveal it negatively regulates lipolysis. Mechanistically, Tld1 LD targeting requires TG, and it is mediated by two distinct hydrophobic regions (HRs). Molecular dynamics simulations reveal that Tld1's HRs interact with TG on LDs and adopt specific conformations on TG-rich LDs versus SE-rich LDs in yeast and human cells. Tld1-deficient yeast display no defect in LD biogenesis but exhibit elevated TG lipolysis dependent on lipase Tgl3. Remarkably, overexpression of Tld1, but not LD protein Pln1/Pet10, promotes TG accumulation without altering SE pools. Finally, we find that Tld1-deficient cells display altered LD mobilization during extended yeast starvation. We propose that Tld1 senses TG-rich LDs and regulates lipolysis on LD subpopulations.

Deletion of ORM2 Causes Oleic Acid-Induced Growth Defects in Saccharomyces cerevisiae

The endoplasmic reticulum (ER) resident proteins of the Orm family (Orm1p and Orm2p) play an essential regulatory role in sphingolipid metabolism and proteostasis of Saccharomyces cerevisiae. Sphingolipid metabolism and its relationship with yeast ORM1 and ORM2 have been studied widely, but its position in phospholipids and neutral lipids requires further studies. We found that the deletion of ORM2 reduced phospholipid levels, but orm1Δ had shown no significant alteration of phospholipids. On the contrary, neutral lipid levels and lipid droplet (LD) numbers were increased in both orm1∆ and orm2∆ cells. Unlike orm1Δ, free fatty acid (FFA) levels were steeply accumulated in orm2∆ cells, and deletion of ORM2 made the cells more sensitive towards oleic acid toxicity. Misregulation of fatty acids has been implicated in the causation of several lipid metabolic disorders. It is imminent to comprehend the control mechanisms of free fatty acid homeostasis and its pathophysiology. Our study has provided experimental evidence of ORM2 role in the lipid and fatty acid metabolism of yeast.

Tld1 is a novel regulator of triglyceride lipolysis that demarcates a lipid droplet subpopulation

Cells store lipids in the form of triglyceride (TG) and sterol-ester (SE) in lipid droplets (LDs). Distinct pools of LDs exist, but a pervasive question is how proteins localize to and convey functions to LD subsets. Here, we show the yeast protein YDR275W/Tld1 (for TG-associated LD protein 1) localizes to a subset of TG-containing LDs, and reveal it negatively regulates lipolysis. Mechanistically, Tld1 LD targeting requires TG, and is mediated by two distinct hydrophobic regions (HRs). Molecular dynamics simulations reveal Tld1 HRs interact with TG on LDs and adopt specific conformations on TG-rich LDs versus SE-rich LDs in yeast and human cells. Tld1-deficient yeast display no defect in LD biogenesis, but exhibit elevated TG lipolysis dependent on lipase Tgl3. Remarkably, over-expression of Tld1, but not LD protein Pln1/Pet10, promotes TG accumulation without altering SE pools. Finally, we find Tld1-deficient cells display altered LD mobilization during extended yeast starvation. We propose Tld1 senses TG-rich LDs and regulates lipolysis on LD subpopulations.

Capturing the Liquid-Crystalline Phase Transformation: Implications for Protein Targeting to Sterol Ester-Rich Lipid Droplets

Lipid droplets are essential organelles that store and traffic neutral lipids. The phospholipid monolayer surrounding their neutral lipid core engages with a highly dynamic proteome that changes according to cellular and metabolic conditions. Recent work has demonstrated that when the abundance of sterol esters increases above a critical concentration, such as under conditions of starvation or high LDL exposure, the lipid droplet core can undergo an amorphous to liquid-crystalline phase transformation. Herein, we study the consequences of this transformation on the physical properties of lipid droplets that are thought to regulate protein association. Using simulations of different sterol-ester concentrations, we have captured the liquid-crystalline phase transformation at the molecular level, highlighting the alignment of sterol esters in alternating orientations to form concentric layers. We demonstrate how ordering in the core permeates into the neutral lipid/phospholipid interface, changing the magnitude and nature of neutral lipid intercalation and inducing ordering in the phospholipid monolayer. Increased phospholipid packing is concomitant with altered surface properties, including smaller area per phospholipid and substantially reduced packing defects. Additionally, the ordering of sterol esters in the core causes less hydration in more ordered regions. We discuss these findings in the context of their expected consequences for preferential protein recruitment to lipid droplets under different metabolic conditions.

Seipin concentrates distinct neutral lipids via interactions with their acyl chain carboxyl esters

Lipid droplets (LDs) are essential for cellular lipid homeostasis by storing diverse neutral lipids (NLs), such as triacylglycerol (TAG), steryl esters (SE), and retinyl esters (RE). A proper assembly of TAG-containing LDs at the ER requires Seipin, a conserved protein often mutated in lipodystrophies. Here, we show that the yeast Seipin Sei1 and its partner Ldb16 also promote the storage of other NL in LDs. Importantly, this role of Sei1/Ldb16 is evolutionarily conserved as expression of human-Seipin restored normal SE-containing LDs in yeast Seipin mutants. As in the case of TAG, the formation of SE-containing LDs requires interactions between hydroxyl-residues in human Seipin or yeast Ldb16 with NL carboxyl esters. These findings provide a universal mechanism for Seipin-mediated LD formation and suggest a model for how Seipin distinguishes NLs from aliphatic phospholipid acyl chains in the center of the membrane bilayer.

Chain flexibility of medicinal lipids determines their selective partitioning into lipid droplets

In guiding lipid droplets (LDs) to serve as storage vessels that insulate high-value lipophilic compounds in cells, we demonstrate that chain flexibility of lipids determines their selective migration in intracellular LDs. Focusing on commercially important medicinal lipids with biogenetic similarity but structural dissimilarity, we computationally and experimentally validate that LD remodeling should be differentiated between overproduction of structurally flexible squalene and that of rigid zeaxanthin and β-carotene. In molecular dynamics simulations, worm-like flexible squalene is readily deformed to move through intertwined chains of triacylglycerols in the LD core, whereas rod-like rigid zeaxanthin is trapped on the LD surface due to a high free energy barrier in diffusion. By designing yeast cells with either much larger LDs or with a greater number of LDs, we observe that intracellular storage of squalene significantly increases with LD volume expansion, but that of zeaxanthin and β-carotene is enhanced through LD surface broadening; as visually evidenced, the outcomes represent internal penetration of squalene and surface localization of zeaxanthin and β-carotene. Our study shows the computational and experimental validation of selective lipid migration into a phase-separated organelle and reveals LD dynamics and functionalization.

Lipid droplet (LD) is a highly dynamic organelle capable of regulating lipid metabolism, storage and transportation. Here, by combining molecular dynamics simulations and microbial LD engineering, the authors demonstrate that the structural flexibility of lipids is one of decisive factors in selective partitioning into LDs.

Efficient Quantification of Lipid Packing Defect Sensing by Amphipathic Peptides: Comparing Martini 2 and 3 with CHARMM36

In biological systems, proteins can be attracted to curved or stretched regions of lipid bilayers by sensing hydrophobic defects in the lipid packing on the membrane surface. Here, we present an efficient end-state free energy calculation method to quantify such sensing in molecular dynamics simulations. We illustrate that lipid packing defect sensing can be defined as the difference in mechanical work required to stretch a membrane with and without a peptide bound to the surface. We also demonstrate that a peptide's ability to concurrently induce excess leaflet area (tension) and elastic softening─a property we call the "characteristic area of sensing" (CHAOS)─and lipid packing sensing behavior are in fact two sides of the same coin. In essence, defect sensing displays a peptide's propensity to generate tension. The here-proposed mechanical pathway is equally accurate yet, computationally, about 40 times less costly than the commonly used alchemical pathway (thermodynamic integration), allowing for more feasible free energy calculations in atomistic simulations. This enabled us to directly compare the Martini 2 and 3 coarse-grained and the CHARMM36 atomistic force fields in terms of relative binding free energies for six representative peptides including the curvature sensor ALPS and two antiviral amphipathic helices (AH). We observed that Martini 3 qualitatively reproduces experimental trends while producing substantially lower (relative) binding free energies and shallower membrane insertion depths compared to atomistic simulations. In contrast, Martini 2 tends to overestimate (relative) binding free energies. Finally, we offer a glimpse into how our end-state-based free energy method can enable the inverse design of optimal lipid packing defect sensing peptides when used in conjunction with our recently developed evolutionary molecular dynamics (Evo-MD) method. We argue that these optimized defect sensors─aside from their biomedical and biophysical relevance─can provide valuable targets for the development of lipid force fields.

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 CYTOLD and ERTOLD pathways for lipid droplet-protein targeting

Lipid droplets (LDs) are the main organelles for lipid storage, and their surfaces contain unique proteins with diverse functions, including those that facilitate the deposition and mobilization of LD lipids. Among organelles, LDs have an unusual structure with an organic, hydrophobic oil phase covered by a phospholipid monolayer. The unique properties of LD monolayer surfaces require proteins to localize to LDs by distinct mechanisms. Here we review the two pathways known to mediate direct LD protein localization: the CYTOLD pathway mediates protein targeting from the cytosol toLDs, and the ERTOLD pathway functions in protein targeting from the endoplasmic reticulum toLDs. We describe the emerging principles for each targeting pathway in animal cells and highlight open questions in the field.

Origin of gradients in lipid density and surface tension between connected lipid droplet and bilayer

We combined theory and experiments to depict physical parameters modulating the phospholipid (PL) density and tension equilibrium between a bilayer and an oil droplet in contiguity. This situation is encountered during a neutral lipid (NL) droplet formation in the endoplasmic reticulum. We set up macroscopic and microscopic models to uncover free parameters and the origin of molecular interactions controlling the PL densities of the droplet monolayer and the bilayer. The established physical laws and predictions agreed with experiments performed with droplet-embedded vesicles. We found that the droplet monolayer is always by a few percent (∼10%) less packed with PLs than the bilayer. Such a density gradient arises from PL-NL interactions on the droplet, which are lower than PL-PL trans interactions in the bilayer, i.e., interactions between PLs belonging to different leaflets of the bilayer. Finally, despite the pseudo-surface tension for the water/PL acyl chains in the bilayer being higher than the water/NL surface tension, the droplet monolayer always has a higher surface tension than the bilayer because of its lower PL density. Thus, a PL density gradient is mandatory to maintain the mechanical and thermodynamic equilibrium of the droplet-bilayer continuity. Our study sheds light on the origin of the molecular interactions responsible for the unique surface properties of lipid droplets compared with cellular bilayer membranes.

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.

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.