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

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.

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.

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.

Cancer immune therapy using engineered ‛tail-flipping' nanoliposomes targeting alternatively activated macrophages

Alternatively-activated, M2-like tumor-associated macrophages (TAM) strongly contribute to tumor growth, invasiveness and metastasis. Technologies to disable the pro-tumorigenic function of these TAMs are of high interest to immunotherapy research. Here we show that by designing engineered nanoliposomes bio-mimicking peroxidated phospholipids that are recognised and internalised by scavenger receptors, TAMs can be targeted. Incorporation of phospholipids possessing a terminal carboxylate group at the sn-2 position into nanoliposome bilayers drives their uptake by M2 macrophages with high specificity. Molecular dynamics simulation of the lipid bilayer predicts flipping of the sn-2 tail towards the aqueous phase, while molecular docking data indicates interaction of the tail with Scavenger Receptor Class B type 1 (SR-B1). In vivo, the engineered nanoliposomes are distributed specifically to M2-like macrophages and, upon delivery of the STAT6 inhibitor (AS1517499), zoledronic acid or muramyl tripeptide, these cells promote reduction of the premetastatic niche and/or tumor growth. Altogether, we demonstrate the efficiency and versatility of our engineered “tail-flipping” nanoliposomes in a pre-clinical model, which paves the way to their development as cancer immunotherapeutics in humans.

Tumor-associated macrophages are mostly pro-tumorigenic, due to their re-programming by the tumor microenvironment. Here authors show that nanoliposomes, incorporating phospholipids with a flipping-tail chain, are engulfed specifically by intratumoral, alternatively activated macrophages, while delivering a cargo that converts these cells into anti-tumor macrophages.

Seipin transmembrane segments critically function in triglyceride nucleation and lipid droplet budding from the membrane

Lipid droplets (LDs) are organelles formed in the endoplasmic reticulum (ER) to store triacylglycerol (TG) and sterol esters. The ER protein seipin is key for LD biogenesis. Seipin forms a cage-like structure, with each seipin monomer containing a conserved hydrophobic helix and two transmembrane (TM) segments. How the different parts of seipin function in TG nucleation and LD budding is poorly understood. Here, we utilized molecular dynamics simulations of human seipin, along with cell-based experiments, to study seipin’s functions in protein–lipid interactions, lipid diffusion, and LD maturation. An all-atom simulation indicates that seipin TM segment residues and hydrophobic helices residues located in the phospholipid tail region of the bilayer attract TG. Simulating larger, growing LDs with coarse-grained models, we find that the seipin TM segments form a constricted neck structure to facilitate conversion of a flat oil lens into a budding LD. Using cell experiments and simulations, we also show that conserved, positively charged residues at the end of seipin’s TM segments affect LD maturation. We propose a model in which seipin TM segments critically function in TG nucleation and LD growth.

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.

Key Factors Governing Initial Stages of Lipid Droplet Formation

Lipid droplets (LDs) are neutral lipid storage organelles surrounded by a phospholipid (PL) monolayer. LD biogenesis from the endoplasmic reticulum is driven by phase separation of neutral lipids, overcoming surface tension and membrane deformation. However, the core biophysics of the initial steps of LD formation remains relatively poorly understood. Here, we use a tunable, phenomenological coarse-grained model to study triacylglycerol (TG) nucleation in a bilayer membrane. We show that PL rigidity has a strong influence on TG lensing and membrane remodeling: when membrane rigidity increases, TG clusters remain more planar with high anisotropy but a minor degree of phase nucleation. This finding is confirmed by advanced sampling simulations that calculate nucleation free energy as a function of the degree of nucleation and anisotropy. We also show that asymmetric tension, controlled by the number of PL molecules on each membrane leaflet, determines the budding direction. A TG lens buds in the direction of the monolayer containing excess PL molecules to allow for better PL coverage of TG, consistent with the reported experiments. Finally, two governing mechanisms of the LD growth, Ostwald ripening and merging, are observed. Taken together, this study characterizes the interplay between two thermodynamic quantities during the initial LD phases, the TG bulk free energy and membrane remodeling energy.

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.

Mechanistic Understanding from Molecular Dynamics in Pharmaceutical Research 2: Lipid Membrane in Drug Design

We review the use of molecular dynamics (MD) simulation as a drug design tool in the context of the role that the lipid membrane can play in drug action, i.e., the interaction between candidate drug molecules and lipid membranes. In the standard "lock and key" paradigm, only the interaction between the drug and a specific active site of a specific protein is considered; the environment in which the drug acts is, from a biophysical perspective, far more complex than this. The possible mechanisms though which a drug can be designed to tinker with physiological processes are significantly broader than merely fitting to a single active site of a single protein. In this paper, we focus on the role of the lipid membrane, arguably the most important element outside the proteins themselves, as a case study. We discuss work that has been carried out, using MD simulation, concerning the transfection of drugs through membranes that act as biological barriers in the path of the drugs, the behavior of drug molecules within membranes, how their collective behavior can affect the structure and properties of the membrane and, finally, the role lipid membranes, to which the vast majority of drug target proteins are associated, can play in mediating the interaction between drug and target protein. This review paper is the second in a two-part series covering MD simulation as a tool in pharmaceutical research; both are designed as pedagogical review papers aimed at both pharmaceutical scientists interested in exploring how the tool of MD simulation can be applied to their research and computational scientists interested in exploring the possibility of a pharmaceutical context for their research.