Lipid simulations: a perspective on lipids in action

Structural and thermodynamic properties of bulk triglycerides and triglyceride/water mixtures reproduced using a polarizable coarse-grained model

Triglycerides (TGs) play important roles in renewable energies, food production, medicine, and metabolism in organisms. Here, we developed a novel coarse-grained (CG) force field (FF) for triglycerides to reproduce both the structural and thermodynamic properties of bulk TGs, TG/air interfaces, and TG/water mixtures using molecular dynamics (MD) simulations. We rigorously optimized the bonded and nonbonded force parameters between the CG beads of TGs and nonbonded force parameters between TG beads and polarizable CG water beads by employing an efficient meta-multilinear interpolation parameterization algorithm recently developed by us. This CG FF performs very well in reproducing the percolating network of the TG bulk phase self-assembled in water and a variety of molecular conformations predicted by all-atom MD simulations. More importantly, it also correctly reproduces multiple experimentally measurable macroscopic thermodynamic properties, including the density and surface tensions of both the TG/air and TG/water interfaces. This paves the way for studying more complicated systems involving TGs on a large scale.

Study of trioleoylglycerol two-layer and adiposome cross-section mimicking four-layer systems through atomic-level simulations

Adiposomes are artificially prepared lipid droplet (LD)-mimetic structures, which, unlike LDs, do not harbor proteins. The dynamics of interaction between triacylglycerols (TAGs), drug molecule, and phospholipids in adiposomes is currently not well-established. Trioleoylglycerol (TOG) molecule was divided into three parts: two oleoyl tails and one 2-monooleoylglycerol (MOG). Forcefield parameters for two oleoyl tails were adopted from the AMBER18 repository while that of the MOG forcefield was taken from the literature. Charge correction was performed on the MOG forcefield before its utilization. After charge correction, the resulting TOG molecule had zero charge. TOG bilayer (2L) and tetralayer (4L) systems were prepared and simulated. TOG bilayer (2L) systems-modeled from two different initial conformations, the TOG3 conformation and the TOG2:1 conformation-showed that TOG2:1 conformation was more prevailing irrespective of the starting conformation and was subsequently used in further simulations. The hydrated TOG 2L system showed TOG-water solution solubility of 0.051 mol L-1 which is near experimental values. This validated the correct parameterization of the TOG molecule. The simulations of 4L systems showed stable membrane behaviors toward the end of simulations. It was also observed that in the 4L system, the TOG molecules showed the formation of micelles with the drug molecule. Almost six TOGs remained continuously in contact with the drug molecule throughout the simulation. The availability of charge-corrected TOG parameterization is expected to equip future studies with a framework for molecular dynamics simulations of adiposomes and/or LDs at the atomic level.

The beginning and the end of SNARE-induced membrane fusion

Membrane fusion is not a spontaneous process. Physiologically, the formation of coiled-coil protein complexes, the SNAREpins, bridges the membrane of a vesicle and a target membrane, brings them in close contact, and provides the energy necessary for their fusion. In this review, we utilize results from in vitro experiments and simple physics and chemistry models to dissect the kinetics and energetics of the fusion process from the encounter of the two membranes to the full expansion of a fusion pore. We find three main energy barriers that oppose the fusion process: SNAREpin initiation, fusion pore opening, and expansion. SNAREpin initiation is inherent to the proteins and makes in vitro fusion kinetic experiments rather slow. The kinetics are physiologically accelerated by effectors. The energy barriers that precede pore opening and pore expansion can be overcome by several SNAREpins acting in concert.

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.

Understanding the interactions between sebum triglycerides and water: a molecular dynamics simulation study

Sebum triglycerides actively contribute to the water transport across the sebum layerviathe percolating network of triglyceride head groups.

The biophysical properties of plasmalogens originating from their unique molecular architecture

Plasmalogens are a unique class of phospholipids that are present in many organisms. Their presence in cell membranes has intrigued researchers for decades due to their unique molecular structure, namely the vinyl-ether bond at the sn-1 position, and their association with brain related disorders. Apparently, based on their amount in the cell membranes, their function is to provide exclusive structural and dynamical properties to these complex molecular assemblies. Yet, many of their physiological roles manifested through their biophysical properties have been challenging to identify. In this review, the biophysical properties of plasmalogens are discussed and compared to other lipid species. The role of plasmalogens is examined in the context of cell membrane function, and some future directions are given.

The mystery of membrane organization: composition, regulation and roles of lipid rafts

Cellular plasma membranes are laterally heterogeneous, featuring a variety of distinct subcompartments that differ in their biophysical properties and composition. A large number of studies have focused on understanding the basis for this heterogeneity and its physiological relevance. The membrane raft hypothesis formalized a physicochemical principle for a subtype of such lateral membrane heterogeneity, in which the preferential associations between cholesterol and saturated lipids drive the formation of relatively packed (or ordered) membrane domains that selectively recruit certain lipids and proteins. Recent studies have yielded new insights into this mechanism and its relevance in vivo, owing primarily to the development of improved biochemical and biophysical technologies.

Hydrocarbons Are Essential for Optimal Cell Size, Division, and Growth of Cyanobacteria

Cyanobacteria are intricately organized, incorporating an array of internal thylakoid membranes, the site of photosynthesis, into cells no larger than other bacteria. They also synthesize C15-C19 alkanes and alkenes, which results in substantial production of hydrocarbons in the environment. All sequenced cyanobacteria encode hydrocarbon biosynthesis pathways, suggesting an important, undefined physiological role for these compounds. Here, we demonstrate that hydrocarbon-deficient mutants of Synechococcus sp. PCC 7002 and Synechocystis sp. PCC 6803 exhibit significant phenotypic differences from wild type, including enlarged cell size, reduced growth, and increased division defects. Photosynthetic rates were similar between strains, although a minor reduction in energy transfer between the soluble light harvesting phycobilisome complex and membrane-bound photosystems was observed. Hydrocarbons were shown to accumulate in thylakoid and cytoplasmic membranes. Modeling of membranes suggests these compounds aggregate in the center of the lipid bilayer, potentially promoting membrane flexibility and facilitating curvature. In vivo measurements confirmed that Synechococcus sp. PCC 7002 mutants lacking hydrocarbons exhibit reduced thylakoid membrane curvature compared to wild type. We propose that hydrocarbons may have a role in inducing the flexibility in membranes required for optimal cell division, size, and growth, and efficient association of soluble and membrane bound proteins. The recent identification of C15-C17 alkanes and alkenes in microalgal species suggests hydrocarbons may serve a similar function in a broad range of photosynthetic organisms.

Data including GROMACS input files for atomistic molecular dynamics simulations of mixed, asymmetric bilayers including molecular topologies, equilibrated structures, and force field for lipids compatible with OPLS-AA parameters

In this Data in Brief article we provide a data package of GROMACS input files for atomistic molecular dynamics simulations of multicomponent, asymmetric lipid bilayers using the OPLS-AA force field. These data include 14 model bilayers composed of 8 different lipid molecules. The lipids present in these models are: cholesterol (CHOL), 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphatidylcholine (POPC), 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphatidylethanolamine (POPE), 1-stearoyl-2-oleoyl-sn-glycero-3-phosphatidyl-ethanolamine (SOPE), 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphatidylserine (POPS), 1-stearoyl-2-oleoyl-sn-glycero-3-phosphatidylserine (SOPS), N-palmitoyl-D-erythro-sphingosyl-phosphatidylcholine (SM16), and N-lignoceroyl-D-erythro-sphingosyl-phosphatidylcholine (SM24). The bilayers׳ compositions are based on lipidomic studies of PC-3 prostate cancer cells and exosomes discussed in Llorente et al. (2013) [1], showing an increase in the section of long-tail lipid species (SOPS, SOPE, and SM24) in the exosomes. Former knowledge about lipid asymmetry in cell membranes was accounted for in the models, meaning that the model of the inner leaflet is composed of a mixture of PC, PS, PE, and cholesterol, while the extracellular leaflet is composed of SM, PC and cholesterol discussed in Van Meer et al. (2008) [2]. The provided data include lipids׳ topologies, equilibrated structures of asymmetric bilayers, all force field parameters, and input files with parameters describing simulation conditions (md.mdp). The data is associated with the research article “Interdigitation of Long-Chain Sphingomyelin Induces Coupling of Membrane Leaflets in a Cholesterol Dependent Manner” (Róg et al., 2016) [3].

A molecular insight into the electro-transfer of small molecules through electropores driven by electric fields

The transport of chemical compounds across the plasma membrane into the cell is relevant for several biological and medical applications. One of the most efficient techniques to enhance this uptake is reversible electroporation. Nevertheless, the detailed molecular mechanism of transport of chemical species (dyes, drugs, genetic materials, …) following the application of electric pulses is not yet fully elucidated. In the past decade, molecular dynamics (MD) simulations have been conducted to model the effect of pulsed electric fields on membranes, describing several aspects of this phenomenon. Here, we first present a comprehensive review of the results obtained so far modeling the electroporation of lipid membranes, then we extend these findings to study the electrotransfer across lipid bilayers subject to microsecond pulsed electric fields of Tat11, a small hydrophilic charged peptide, and of siRNA. We use in particular a MD simulation protocol that allows to characterize the transport of charged species through stable pores. Unexpectedly, our results show that for an electroporated bilayer subject to transmembrane voltages in the order of 500mV, i.e. consistent with experimental conditions, both Tat11 and siRNA can translocate through nanoelectropores within tens of ns. We discuss these results in comparison to experiments in order to rationalize the mechanism of drug uptake by cells. This article is part of a Special Issue entitled: Biosimulations edited by Ilpo Vattulainen and Tomasz Róg.

Role of charged lipids in membrane structures - Insight given by simulations

Lipids and proteins are the main components of cell membranes. It is becoming increasingly clear that lipids, in addition to providing an environment for proteins to work in, are in many cases also able to modulate the structure and function of those proteins. Particularly charged lipids such as phosphatidylinositols and phosphatidylserines are involved in several examples of such effects. Molecular dynamics simulations have proved an invaluable tool in exploring these aspects. This so-called computational microscope can provide both complementing explanations for the experimental results and guide experiments to fruitful directions. In this paper, we review studies that have utilized molecular dynamics simulations to unravel the roles of charged lipids in membrane structures. We focus on lipids as active constituents of the membranes, affecting both general membrane properties as well as non-lipid membrane components, mainly proteins. This article is part of a Special Issue entitled: Biosimulations edited by Ilpo Vattulainen and Tomasz Róg.

Functionalized lipids and surfactants for specific applications

Synthetic lipids and surfactants that do not exist in biological systems have been used for the last few decades in both basic and applied science. The most notable applications for synthetic lipids and surfactants are drug delivery, gene transfection, as reporting molecules, and as support for structural lipid biology. In this review, we describe the potential of the synergistic combination of computational and experimental methodologies to study the behavior of synthetic lipids and surfactants embedded in lipid membranes and liposomes. We focused on select cases in which molecular dynamics simulations were used to complement experimental studies aiming to understand the structure and properties of new compounds at the atomistic level. We also describe cases in which molecular dynamics simulations were used to design new synthetic lipids and surfactants, as well as emerging fields for the application of these compounds. This article is part of a Special Issue entitled: Biosimulations edited by Ilpo Vattulainen and Tomasz Róg.

Computer modelling studies of the bilayer/water interface

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

Characterization of thylakoid lipid membranes from cyanobacteria and higher plants by molecular dynamics simulations

The thylakoid membrane is mainly composed of non-common lipids, so called galactolipids. Despite the importance of these lipids for the function of the photosynthetic reaction centers, the molecular organization of these membranes is largely unexplored. Here we use multiscale molecular dynamics simulations to characterize the thylakoid membrane of both cyanobacteria and higher plants. We consider mixtures of up to five different galactolipids plus phosphatidylglycerol to represent these complex membranes. We find that the different lipids generally mix well, although nanoscale heterogeneities are observed especially in case of the plant membrane. The fluidity of the cyanobacterial membrane is markedly reduced compared to the plant membrane, even considering elevated temperatures at which thermophilic cyanobacteria are found. We also find that the plant membrane more readily undergoes a phase transformation to an inverted hexagonal phase. We furthermore characterized the conformation and dynamics of the cofactors plastoquinone and plastoquinol, revealing of the fast flip-flop rates for the non-reduced form. Together, our results provide a molecular view on the dynamical organization of the thylakoid membrane.

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.

Microsecond molecular dynamics simulations of lipid mixing

Molecular dynamics (MD) simulations of membranes are often hindered by the slow lateral diffusion of lipids and the limited time scale of MD. In order to study the dynamics of mixing and characterize the lateral distribution of lipids in converged mixtures, we report microsecond-long all-atom MD simulations performed on the special-purpose machine Anton. Two types of mixed bilayers, POPE:POPG (3:1) and POPC:cholesterol (2:1), as well as a pure POPC bilayer, were each simulated for up to 2 μs. These simulations show that POPE:POPG and POPC:cholesterol are each fully miscible at the simulated conditions, with the final states of the mixed bilayers similar to a random mixture. By simulating three POPE:POPG bilayers at different NaCl concentrations (0, 0.15, and 1 M), we also examined the effect of salt concentration on lipid mixing. While an increase in NaCl concentration is shown to affect the area per lipid, tail order, and lipid lateral diffusion, the final states of mixing remain unaltered, which is explained by the largely uniform increase in Na(+) ions around POPE and POPG. Direct measurement of water permeation reveals that the POPE:POPG bilayer with 1 M NaCl has reduced water permeability compared with those at zero or low salt concentration. Our calculations provide a benchmark to estimate the convergence time scale of all-atom MD simulations of lipid mixing. Additionally, equilibrated structures of POPE:POPG and POPC:cholesterol, which are frequently used to mimic bacterial and mammalian membranes, respectively, can be used as starting points of simulations involving these membranes.

Molecular recognition of gangliosides and their potential for cancer immunotherapies

Gangliosides are sialic-acid-containing glycosphingolipids expressed on all vertebrate cells. They are primarily positioned in the plasma membrane with the ceramide part anchored in the membrane and the glycan part exposed on the surface of the cell. These lipids have highly diverse structures, not the least with respect to their carbohydrate chains, with N-acetylneuraminic acid (NeuAc) and N-glycolylneuraminic acid (NeuGc) being the two most common sialic-acid residues in mammalian cells. Generally, human healthy tissue is deficient in NeuGc, but this molecule is expressed in tumors and in human fetal tissues, and was hence classified as an onco-fetal antigen. Gangliosides perform important functions through carbohydrate-specific interactions with proteins, for example, as receptors in cell–cell recognition, which can be exploited by viruses and other pathogens, and also by regulating signaling proteins, such as the epidermal growth factor receptor (EGFR) and the vascular endothelial growth factor receptor (VEGFR), through lateral interaction in the membrane. Through both mechanisms, tumor-associated gangliosides may affect malignant progression, which makes them attractive targets for cancer immunotherapies. In this review, we describe how proteins recognize gangliosides, focusing on the molecular recognition of gangliosides associated with cancer immunotherapy, and discuss the importance of these molecules in cancer research.

Methodologies for the analysis of instantaneous lipid diffusion in MD simulations of large membrane systems

Interactions between lipids and membrane proteins play a key role in determining the nanoscale dynamic and structural properties of biological membranes. Molecular dynamics (MD) simulations provide a valuable tool for studying membrane models, complementing experimental approaches. It is now possible to simulate large membrane systems, such as simplified models of bacterial and viral envelope membranes. Consequently, there is a pressing need to develop tools to visualize and quantify the dynamics of these immense systems, which typically comprise millions of particles. To tackle this issue, we have developed visual and quantitative analyses of molecular positions and their velocity field using path line, vector field and streamline techniques. This allows us to highlight large, transient flow-like movements of lipids and to better understand crowding within the lipid bilayer. The current study focuses on visualization and analysis of lipid dynamics. However, the methods are flexible and can be readily applied to e.g. proteins and nanoparticles within large complex membranes. The protocols developed here are readily accessible both as a plugin for the molecular visualization program VMD and as a module for the MDAnalysis library.

Lipid nanotechnology

Nanotechnology is a multidisciplinary field that covers a vast and diverse array of devices and machines derived from engineering, physics, materials science, chemistry and biology. These devices have found applications in biomedical sciences, such as targeted drug delivery, bio-imaging, sensing and diagnosis of pathologies at early stages. In these applications, nano-devices typically interface with the plasma membrane of cells. On the other hand, naturally occurring nanostructures in biology have been a source of inspiration for new nanotechnological designs and hybrid nanostructures made of biological and non-biological, organic and inorganic building blocks. Lipids, with their amphiphilicity, diversity of head and tail chemistry, and antifouling properties that block nonspecific binding to lipid-coated surfaces, provide a powerful toolbox for nanotechnology. This review discusses the progress in the emerging field of lipid nanotechnology.

Chemical biology approaches to membrane homeostasis and function

The study of membranes is at a turning point. New theories about membrane structure and function have recently been proposed, however, new technologies, combining chemical, physical, and biochemical approaches are necessary to test these hypotheses. In particular, the NCCR in chemical biology aims to visualize and characterize membrane microdomains and determine their function during hormone signaling.

Receptor-mediated endocytosis of nanoparticles of various shapes

Cellular uptake through endocytosis is crucial for drug delivery and nanomedicine. However, the conditions under which passive endocytosis (i.e., not ATP driven) takes place are not well understood. We report MD simulations of the passive uptake of ligand-coated nanoparticles with varying size, shape, coverage, and membrane-binding strength. We find that the efficiency of passive endocytosis is higher for spherocylindrical particles than for spheres and that endocytosis is suppressed for particles with sharp edges.

Membrane-lipid therapy in operation: the HSP co-inducer BGP-15 activates stress signal transduction pathways by remodeling plasma membrane rafts

Aging and pathophysiological conditions are linked to membrane changes which modulate membrane-controlled molecular switches, causing dysregulated heat shock protein (HSP) expression. HSP co-inducer hydroxylamines such as BGP-15 provide advanced therapeutic candidates for many diseases since they preferentially affect stressed cells and are unlikely have major side effects. In the present study in vitro molecular dynamic simulation, experiments with lipid monolayers and in vivo ultrasensitive fluorescence microscopy showed that BGP-15 alters the organization of cholesterol-rich membrane domains. Imaging of nanoscopic long-lived platforms using the raft marker glycosylphosphatidylinositol-anchored monomeric green fluorescent protein diffusing in the live Chinese hamster ovary (CHO) cell plasma membrane demonstrated that BGP-15 prevents the transient structural disintegration of rafts induced by fever-type heat stress. Moreover, BGP-15 was able to remodel cholesterol-enriched lipid platforms reminiscent of those observed earlier following non-lethal heat priming or membrane stress, and were shown to be obligate for the generation and transmission of stress signals. BGP-15 activation of HSP expression in B16-F10 mouse melanoma cells involves the Rac1 signaling cascade in accordance with the previous observation that cholesterol affects the targeting of Rac1 to membranes. Finally, in a human embryonic kidney cell line we demonstrate that BGP-15 is able to inhibit the rapid heat shock factor 1 (HSF1) acetylation monitored during the early phase of heat stress, thereby promoting a prolonged duration of HSF1 binding to heat shock elements. Taken together, our results indicate that BGP-15 has the potential to become a new class of pharmaceuticals for use in ‘membrane-lipid therapy’ to combat many various protein-misfolding diseases associated with aging.