Atomistic simulations of phosphatidylcholines and cholesteryl esters in high-density lipoprotein-sized lipid droplet and trilayer: clues to cholesteryl ester transport and storage

CHARMM-GUI Membrane Builder for Lipid Droplet Modeling and Simulation

Lipid droplets (LDs) are organelles that are necessary for eukaryotic and prokaryotic metabolism and energy storage. They have a unique structure consisting of a spherical phospholipid monolayer encasing neutral lipids such as triacylglycerol (TAG). LDs have garnered increased interest for their implications in disease and for drug delivery applications. Consequently, there is an increased need for tools to study their structure, composition, and dynamics in biological contexts. In this work, we utilize CHARMM-GUI Membrane Builder to simulate and analyze LDs with and without a plant LD protein, oleosin. The results show that Membrane Builder can generate biologically relevant all-atom LD systems with relatively short equilibration times using a new TAG library having optimized headgroup parameters. TAG molecules originally inserted into a lipid bilayer aggregate in the membrane center, forming a TAG-only core flanked by two monolayers. The TAG-only core thickness stably grows with increasing TAG mole fraction. A 70 % TAG system has a core that is thick enough to house oleosin without its interactions with the distal leaflet or disruption of its secondary structure. We hope that Membrane Builder can aid in the future study of LD systems, including their structure and dynamics with and without proteins.

Hepatic cholesterol biosynthesis and dioxin-induced dysregulation: A multiscale computational approach

Humans are constantly exposed to lipophilic persistent organic pollutants (POPs) that accumulate in fatty foods. Among the numerous POPs, dioxins, in particular 2,3,7,8-tetrachlorodibenzo-p-dioxin (TCDD), can impact several organ systems. While the hazard is clearly recognized, it is still difficult to develop a comprehensive understanding of the overall health impacts of dioxins. As chemical toxicity testing is steadily adopting new approach methodologies (NAMs), it becomes imperative to develop computational models that can bridge the data gaps between in vitro testing and in vivo outcomes. As an effort to address this challenge, we propose a multiscale computational approach using a "template-and-anchor" (T&A) structure. A template is a high-level umbrella model that permits the integration of information from various, detailed anchor models. In the present study, we use this T&A approach to describe the effect of TCDD on cholesterol dynamics. Specifically, we represent hepatic cholesterol biosynthesis as an anchor model that is perturbed by TCDD, leading to steatosis, along with alterations of plasma cholesterol. In the future, incorporating pertinent information from all anchor models into the template model will allow the characterization of the global effects of dioxin, which can subsequently be translated into overall - and ultimately personalized - human health risk assessment.

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.

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.

Profiling of non-polar lipids in tears of contact lens wearers during the day

PURPOSE

This study explored whether the non-polar lipids in the human tear fluid lipidome show diurnal variation with and without contact lens wear. It also addressed the relationship between changes in ocular comfort during the day with the level of non-polar lipids.

METHODS

Tear samples were collected in the morning and evening with and without contact lenses using fine glass capillary tubes and were analysed by chip-based nano-electrospray ionization tandem mass spectrometric techniques. Tear levels of cholesteryl esters (CE), wax esters (WE) and triacylglycerides (TAG) were quantified.

RESULTS

TAG 48:0, 52:0 and WE 26:0/16:0, and 27:0/17:0 increased from morning to evening. TAG 52:2, WE 21:0/16:0, 21:0/18:1 and 28:0/18:1 decreased during the day when no lenses were worn. CE 21:0 was the only non-polar lipid that increased from morning to evening in contact lens wear. WE 21:0/16:0 and 27:0/17:0 were lower in the morning in contact lens wear compared to no lens wear (p ≤ 0.05). The level of non-polar lipids did not correlate with ocular comfort at the end of the day.

CONCLUSION

Even though the level of some of non-polar lipid species changed from morning to evening the total level of major tear non-polar lipids remained unchanged during the day with and without contact lens wear. The effect of change in the quantity and structure of lipid species on tear stability and ocular comfort warrants more investigation.

Structure and intermolecular interactions in spheroidal high-density lipoprotein subpopulations

•

High-density lipoprotein subpopulations have unique surface profiles and dynamics.

•

Relative hydrophobic surface area decreases with increasing lipoprotein size.

•

Core lipid exposure at the lipoprotein surface decreases with increasing size.

•

Cholesterol molecules localise near apolipoprotein A-I central helices.

•

Lipid and protein interactions stabilise multifoil models of apolipoprotein A-I.

Human serum high-density lipoproteins (HDLs) are a population of small, dense protein-lipid aggregates that are crucial for intravascular lipid trafficking and are protective against cardiovascular disease. The spheroidal HDL subfraction can be separated by size and density into five major subpopulations with distinct molecular compositions and unique biological functionalities: HDL_3c, HDL_3b, HDL_3a, HDL_2a and HDL_2b. Representative molecular models of these five subpopulations were developed and characterised for the first time in the presence of multiple copies of its primary protein component apolipoprotein A-I (apoA-I) using coarse-grained molecular dynamics simulations. Each HDL model exhibited size, morphological and compositional profiles consistent with experimental observables. With increasing particle size the separation of core and surface molecules became progressively more defined, resulting in enhanced core lipid mixing, reduced core lipid exposure at the surface, and the formation of an interstitial region between core and surface molecules in HDL_2b. Cholesterol molecules tended to localise around the central helix-5 of apoA-I, whilst triglyceride molecules predominantly interacted with aromatic, hydrophobic residues located within the terminal helix-10 across all subpopulation models. The three intermediate HDL models exhibited similar surface profiles despite having distinct molecular compositions. ApoA-I in trefoil, quatrefoil and pentafoil arrangements across the surface of HDL particles exhibited significant warping and twisting, but largely retained intermolecular contacts between adjacent apoA-I chains. Representative HDL subpopulations differed in particle size, morphology, intermolecular interaction profiles and lipid and protein dynamics. These findings reveal how different HDL subpopulations might exhibit distinct functional associations depending on particle size, form and composition.

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.

Time-Resolved Fluorescence Anisotropy of a Molecular Rotor Resolves Microscopic Viscosity Parameters in Complex Environments

Understanding viscosity in complex environments remains a largely unanswered question despite its importance in determining reaction rates in vivo. Here, time-resolved fluorescence anisotropy imaging (TR-FAIM) is combined with fluorescent molecular rotors (FMRs) to simultaneously determine two non-equivalent viscosity-related parameters in complex heterogeneous environments. The parameters, FMR rotational correlation time and lifetime, are extracted from fluorescence anisotropy decays, which in heterogeneous environments show dip-and-rise behavior due to multiple dye populations. Decays of this kind are found both in artificially constructed adiposomes and in live cell lipid droplet organelles. Molecular dynamics simulations are used to assign each population to nano-environments within the lipid systems. The less viscous population corresponds to the state showing an average 25° tilt to the lipid membrane normal, and the more viscous population to the state showing an average 55° tilt. This combined experimental and simulation approach enables a comprehensive description of the FMR probe behavior within viscous nano-environments in complex, biological systems.

Multiscale Simulations of Biological Membranes: The Challenge To Understand Biological Phenomena in a Living Substance

Biological membranes are tricky to investigate. They are complex in terms of molecular composition and structure, functional over a wide range of time scales, and characterized by nonequilibrium conditions. Because of all of these features, simulations are a great technique to study biomembrane behavior. A significant part of the functional processes in biological membranes takes place at the molecular level; thus computer simulations are the method of choice to explore how their properties emerge from specific molecular features and how the interplay among the numerous molecules gives rise to function over spatial and time scales larger than the molecular ones. In this review, we focus on this broad theme. We discuss the current state-of-the-art of biomembrane simulations that, until now, have largely focused on a rather narrow picture of the complexity of the membranes. Given this, we also discuss the challenges that we should unravel in the foreseeable future. Numerous features such as the actin-cytoskeleton network, the glycocalyx network, and nonequilibrium transport under ATP-driven conditions have so far received very little attention; however, the potential of simulations to solve them would be exceptionally high. A major milestone for this research would be that one day we could say that computer simulations genuinely research biological membranes, not just lipid bilayers.

Molecular dynamics simulations of lipid nanodiscs

A lipid nanodisc is a discoidal lipid bilayer stabilized by proteins, peptides, or polymers on its edge. Nanodiscs have two important connections to structural biology. The first is associated with high-density lipoprotein (HDL), a particle with a variety of functionalities including lipid transport. Nascent HDL (nHDL) is a nanodisc stabilized by Apolipoprotein A-I (APOA1). Determining the structure of APOA1 and its mimetic peptides in nanodiscs is crucial to understanding pathologies related to HDL maturation and designing effective therapies. Secondly, nanodiscs offer non-detergent membrane-mimicking environments and greatly facilitate structural studies of membrane proteins. Although seemingly similar, natural and synthetic nanodiscs are different in that nHDL is heterogeneous in size, due to APOA1 elasticity, and gradually matures to become spherical. Synthetic nanodiscs, in contrast, should be homogenous, stable, and size-tunable. This report reviews previous molecular dynamics (MD) simulation studies of nanodiscs and illustrates convergence and accuracy issues using results from new multi-microsecond atomistic MD simulations. These new simulations reveal that APOA1 helices take 10-20 μs to rearrange on the nanodisc, while peptides take 2 μs to migrate from the disc surfaces to the edge. These systems can also become kinetically trapped depending on the initial conditions. For example, APOA1 was trapped in a biologically irrelevant conformation for the duration of a 10 μs trajectory; the peptides were similarly trapped for 5 μs. It therefore remains essential to validate MD simulations of these systems with experiments due to convergence and accuracy issues. This article is part of a Special Issue entitled: Emergence of Complex Behavior in Biomembranes edited by Marjorie Longo.

Mechanism and Determinants of Amphipathic Helix-Containing Protein Targeting to Lipid Droplets

Cytosolic lipid droplets (LDs) are the main storage organelles for metabolic energy in most cells. They are unusual organelles that are bounded by a phospholipid monolayer and specific surface proteins, including key enzymes of lipid and energy metabolism. Proteins targeting LDs from the cytoplasm often contain amphipathic helices, but how they bind to LDs is not well understood. Combining computer simulations with experimental studies in vitro and in cells, we uncover a general mechanism for targeting of cytosolic proteins to LDs: large hydrophobic residues of amphipathic helices detect and bind to large, persistent membrane packing defects that are unique to the LD surface. Surprisingly, amphipathic helices with large hydrophobic residues from many different proteins are capable of binding to LDs. This suggests that LD protein composition is additionally determined by mechanisms that selectively prevent proteins from binding LDs, such as macromolecular crowding at the LD surface.

Interdigitation between Triglycerides and Lipids Modulates Surface Properties of Lipid Droplets

Intracellular lipid droplets (LDs) are the main cellular site of metabolic energy storage. Their structure is unique inside the cell, with a core of esterified fatty acids and sterols, mainly triglycerides and sterol esters, surrounded by a single monolayer of phospholipids. Numerous peripheral proteins, including several that were previously associated with intracellular compartments surrounded by a lipid bilayer, have been recently shown to target the surface of LDs, but how they are able to selectively target this organelle remains largely unknown. Here, we use atomistic and coarse-grained molecular dynamics simulations to investigate the molecular properties of the LD surface and to characterize how it differs from that of a lipid bilayer. Our data suggest that although several surface properties are remarkably similar between the two structures, key differences originate from the interdigitation between surface phospholipids and core neutral lipids that occurs in LDs. This property is extremely sensitive to membrane undulations, unlike in lipid bilayers, and it strongly affects both lipid-packing defects and the lateral pressure profile. We observed a marked change in overall surface properties for surface tensions >10 mN/m, indicative of a bimodal behavior. Our simulations provide a comprehensive molecular characterization of the unique surface properties of LDs and suggest how the molecular properties of the surface lipid monolayer can be modulated by the underlying neutral lipids.

Identification of a novel lipid binding motif in apolipoprotein B by the analysis of hydrophobic cluster domains

Apolipoprotein B (apoB) is a large amphipathic protein that is the structural scaffold for the formation of several classes of lipoproteins involved in lipid transport throughout the body. The goal of the present study was to identify specific domains in the apoB sequence that contribute to its lipid binding properties. A sequence analysis algorithm was developed to identify stretches of hydrophobic amino acids devoid of charged amino acids, which are referred to as hydrophobic cluster domains (HCDs). This analysis identified 78 HCDs in apoB with hydrophobic stretches ranging from 6 to 26 residues. Each HCD was analyzed in silico for secondary structure and lipid binding properties, and a subset was synthesized for experimental evaluation. One HCD peptide, B38, showed high affinity binding to both isolated HDL and LDL, and could exchange between lipoproteins. All-atom molecular dynamics simulations indicate that B38 inserts 3.7Å below the phosphate plane of the bilayer. B38 forms an unusual α-helix with a broad hydrophobic face and polar serine and threonine residues on the opposite face. Based on this structure, we hypothesized that B38 could efflux cholesterol from cells. B38 showed a 12-fold greater activity than the 5A peptide, a bihelical Class A amphipathic helix (EC₅₀ of 0.2658 vs. 3.188μM; p<0.0001), in promoting cholesterol efflux from ABCA1 expressing BHK-1 cells. In conclusion, we have identified novel domains within apoB that contribute to its lipid biding properties. Additionally, we have discovered a unique amphipathic helix design for efficient ABCA1-specific cholesterol efflux.

Computational studies of plasma lipoprotein lipids

Plasma lipoproteins are macromolecular assemblies of proteins and lipids found in the blood. The lipid components of lipoproteins are amphipathic lipids such as phospholipids (PLs), and unesterified cholesterols (UCs) and hydrophobic lipids such as cholesteryl esters (CEs) and triglycerides (TGs). Since lipoproteins are soft matter supramolecular assemblies easily deformable by thermal fluctuations and they also exist in varying densities and protein/lipid components, a detailed understanding of their structure/function is experimentally difficult. Molecular dynamics (MD) simulation has emerged as a particularly promising way to explore the structure and dynamics of lipoproteins. The purpose of this review is to survey the current status of computational studies of the lipid components of the lipoproteins. Computational studies aim to explore three levels of complexity for the 3-dimensional structural dynamics of lipoproteins at various metabolic stages: (i) lipoprotein particles consist of protein with minimal lipid; (ii) lipoprotein particles consist of PL-rich discoidal bilayer-like lipid particles; (iii) mature circulating lipoprotein particles consist of CE-rich or TG-rich spheroidal lipid-droplet-like particles. Due to energy barriers involved in conversion between these species, other biomolecules also participate in lipoprotein biological assembly. For example: (i) lipid-poor apolipoprotein A-I (apoA-I) interacts with ATP-binding cassette transporter A1 (ABCA1) to produce nascent discoidal high density lipoprotein (dHDL) particles; (ii) lecithin-cholesterol acyltransferase (LCAT) mediates the conversion of UC to CE in dHDL, driving spheroidal HDL (sHDL) formation; (iii) transfer proteins, cholesterol ester transfer protein (CETP) and phospholipid transfer protein (PLTP), transfer both CE and TG and PL, respectively, between lipoprotein particles. Computational studies have the potential to explore different lipoprotein particles at each metabolic stage in atomistic detail. This review discusses the current status of computational methods including all-atom MD (AAMD), coarse-grain MD (CGMD), and MD-simulated annealing (MDSA) and their applications in lipoprotein structural dynamics and biological assemblies. Results from MD simulations are discussed and compared across studies in order to identify key findings, controversies, issues and future directions. This article is part of a Special Issue entitled: Biosimulations edited by Ilpo Vattulainen and Tomasz Róg.

Tear film lipid layer: A molecular level view

Human cornea is covered by an aqueous tear film, and the outermost layer of the tear film is coated by lipids. This so-called tear film lipid layer (TFLL) reduces surface tension of the tear film and helps with the film re-spreading after blinks. Alterations of tear lipids composition and properties are related to dry eye syndrome. Therefore, unveiling structural and functional properties of TFLL is necessary for understanding tear film function under both normal and pathological conditions. Key properties of TFLL, such as resistance against high lateral pressures and ability to spread at the tear film surface, are directly related to the chemical identity of TFLL lipids. Hence, a molecular-level description is required to get better insight into TFLL properties. Molecular dynamics simulations are particularly well suited for this task and they were recently used for investigating TFLL. The present review discusses molecular level organization and properties of TFLL as seen by these simulation studies. This article is part of a Special Issue entitled: Biosimulations edited by Ilpo Vattulainen and Tomasz Róg.

Structure and Dynamics of Phospholipid Nanodiscs from All-Atom and Coarse-Grained Simulations

We investigated structural and dynamical properties of nanodiscs comprising dimyristoylphosphatidylcholine (DMPC) lipids and major scaffold protein MSP1Δ(1-22) from human apolipoprotein A-1 using combined all-atom and coarse-grained (CG) molecular dynamics (MD) simulations. The computational efficiency of the Martini-CG force field enables the spontaneous self-assembly of lipids and scaffold proteins into stable nanodisc structures on time scales up to tens of microseconds. Subsequent all-atom and CG-MD simulations reveal that the lipids in the nanodisc have lower configurational entropy and higher acyl tail order than in a lamellar bilayer phase. These altered average properties arise from rather differential behavior of lipids, depending on their location in the nanodisc. Since the scaffold proteins exert constrictive forces from the outer rim of the disc toward its center, lipids at the center of the nanodisc are highly ordered, whereas annular lipids that are in contact with the MSP proteins are remarkably disordered due to perturbed packing. Although specific differences between all-atom and CG simulations are also evident, the results obtained at both levels of resolution are in overall good agreement with each other and provide atomic level interpretations of recent experiments. Thus, the present study highlights the applicability of multiscale simulation approaches for nanodisc systems and opens the way for future applications, including the study of nanodisc-embedded membrane proteins.

Atomistic MD simulation reveals the mechanism by which CETP penetrates into HDL enabling lipid transfer from HDL to CETP

Inhibition of cholesterol ester transfer protein (CETP), a protein mediating transfer of neutral lipids between lipoproteins, has been proposed as a means to elevate atheroprotective HDL subpopulations and thereby reduce atherosclerosis. However, off-target and adverse effects of the inhibition have raised doubts about the molecular mechanism of CETP-HDL interaction. Recent experimental findings have demonstrated the penetration of CETP into HDL. However, atomic level resolution of CETP penetration into HDL, a prerequisite for a better understanding of CETP functionality and HDL atheroprotection, is missing. We constructed an HDL particle that mimics the actual human HDL mass composition and investigated for the first time, by large-scale atomistic molecular dynamics, the interaction of an upright CETP with a human HDL-mimicking model. The results demonstrated how CETP can penetrate the HDL particle surface, with the formation of an opening in the N barrel domain end of CETP, put in evidence the major anchoring role of a tryptophan-rich region of this domain, and unveiled the presence of a phenylalanine barrier controlling further access of HDL-derived lipids to the tunnel of CETP. The findings reveal novel atomistic details of the CETP-HDL interaction mechanism and can provide new insight into therapeutic strategies.

Distribution of neutral lipids in the lipid droplet core

Cholesteryl esters (CEs) are a form of cholesterol (CHOL) storage in the living cells, as opposed to free CHOL. CEs are major constituents of low density lipoprotein particles. Therefore, CEs are implicated in provoking atherosclerosis. Arranged into cytoplasmic lipid droplets (LDs), CEs are stored intracellularly. They can also be transported extracellularly by means of lipoproteins. In this work, large-scale molecular dynamics (MD) simulations are used to characterize the molecular structure of LDs containing various fractions (10-50 mol %) of cholesteryl oleate (CO) with respect to triolein (TO) fraction. The simulated LDs were covered by a phospholipid monolayer formed by a mixture of 1-palmitoyl-2-oleoylphosphatidylcholine, POPC (75 mol %), and 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphoethanolamine, POPE (25 mol %), molecules. We report that most CO molecules are located within the hydrophobic core of LDs, whereas a small fraction (0.3-1.9 mol %) penetrates the monolayer. The solubility of CO in the phospholipid monolayer is relatively small. Due to a good miscibility with TO molecules, CO forms a liquid phase inside LD at 333 K. There is long-range order in the liquid TO-CO droplet core up to 8 nm from the phospholipid interface, resulting from the structuring of hydrophilic groups. This structuring slowly decays in the direction toward the LD center of mass. No sorting of TO and CO is detected, irrespective of the molar fractions simulated. The distribution of CO within the LDs is significant in determining the rate of their hydrolysis by surface-bound enzyme lipases, and thus has a subsequent impact on the levels of CO in plasma and LDLs.

Core-based lipid nanoparticles as a nanoplatform for delivery of near-infrared fluorescent imaging agents

Pyropheophorbide a (Pyro) is a near-infrared (NIR) fluorescent dye and photosensitizer with high quantum yield that makes the dye suitable for tumor treatment both as an imaging and therapy agent. We have designed and synthesized a series of a Pyro-based NIR probes, based on the conjugation of Pyro with lipids. The nature of our probes requires the use of a lipophilic carrier to deliver the probes to cancer cell membranes. To address this, we have utilized lipid-based nanoparticles (LNPs) consisting of PEGylated lipids, which form the nanoparticle shell, and a lipid core. To endow the LNPs with targeting properties, nitrilotriacetic acid (NTA) lipids were included in the composition that enables the non-covalent attachment of His-tag targeting proteins preserving their functional activity. We found that the nature of the core molecules influence the nanoparticle size, shelf-life and stability at physiological temperature. Two different Pyro-lipid conjugates were loaded either into the core or shell of the LNPs. The conjugates revealed differential ability to be accumulated in the cell membrane of the target cells with time. Thus, the modular organization of the core-shell LNPs allows facile adjustment of their composition with goal to fine tuning the nanoparticle properties for in vivo application.

Interfacial properties of high-density lipoprotein-like lipid droplets with different lipid and apolipoprotein A-I compositions

The surface properties of high-density lipoproteins (HDLs) are important because different enzymes bind and carry out their functions at the surface of HDL particles during metabolic processes. However, the surface properties of HDL and other lipoproteins are poorly known because they cannot be directly measured for nanoscale particles with contemporary experimental methods. In this work, we carried out coarse-grained molecular dynamics simulations to study the concentration of core lipids in the surface monolayer and the interfacial tension of droplets resembling HDL particles. We simulated lipid droplets composed of different amounts of phospholipids, cholesterol esters (CEs), triglycerides (TGs), and apolipoprotein A-Is. Our results reveal that the amount of TGs in the vicinity of water molecules in the phospholipid monolayer is 25-50% higher compared to the amount of CEs in a lipid droplet with a mixed core of an equal amount of TG and CE. In addition, the correlation time for the exchange of molecules between the core and the monolayer is significantly longer for TGs compared to CEs. This suggests that the chemical potential of TG is lower in the vicinity of aqueous phase but the free-energy barrier for the translocation between the monolayer and the core is higher compared to CEs. From the point of view of enzymatic modification, this indicates that TG molecules are more accessible from the aqueous phase. Further, our results point out that CE molecules decrease the interfacial tension of HDL-like lipid droplets whereas TG keeps it constant while the amount of phospholipids varies.

Interfacial tension and surface pressure of high density lipoprotein, low density lipoprotein, and related lipid droplets

Lipid droplets play a central role in energy storage and metabolism on a cellular scale. Their core is comprised of hydrophobic lipids covered by a surface region consisting of amphiphilic lipids and proteins. For example, high and low density lipoproteins (HDL and LDL, respectively) are essentially lipid droplets surrounded by specific proteins, their main function being to transport cholesterol. Interfacial tension and surface pressure of these particles are of great interest because they are related to the shape and the stability of the droplets and to protein adsorption at the interface. Here we use coarse-grained molecular-dynamics simulations to consider a number of related issues by calculating the interfacial tension in protein-free lipid droplets, and in HDL and LDL particles mimicking physiological conditions. First, our results suggest that the curvature dependence of interfacial tension becomes significant for particles with a radius of ∼5 nm, when the area per molecule in the surface region is <1.4 nm(2). Further, interfacial tensions in the used HDL and LDL models are essentially unaffected by single apo-proteins at the surface. Finally, interfacial tensions of lipoproteins are higher than in thermodynamically stable droplets, suggesting that HDL and LDL are kinetically trapped into a metastable state.

Role of neutral lipids in tear fluid lipid layer: coarse-grained simulation study

Tear fluid lipid layer (TFLL) residing at the air-water interface of tears has been recognized to play an important role in the development of dry eye syndrome. Yet, the composition, structure, and mechanical properties of TFLL are only partly known. Here, we report results of coarse-grained simulations of a lipid layer comprising phospholipids, free fatty acids, cholesteryl esters, and triglycerides at the air-water interface to shed light on the properties of TFLL. We consider structural as well as dynamical properties of the lipid layer as a function of surface pressure. Simulations revealed that neutral lipids reside heterogeneously between phospholipids at relatively low pressures but form a separate hydrophobic phase with increasing surface pressure, transforming the initial lipid monolayer to a two-layered structure. When the model of TFLL was compared to a one-component phospholipid monolayer system, we found drastic differences in both structural and dynamical properties that explain the prominent role of neutral lipids as stabilizers of the TFLL. Based on our results, we suggest that neutral lipids are able to increase the stability of the TFLL by modulating its dynamical and structural behavior, which is important for the proper function of tear film.

Lipid exchange mechanism of the cholesteryl ester transfer protein clarified by atomistic and coarse-grained simulations

Cholesteryl ester transfer protein (CETP) transports cholesteryl esters, triglycerides, and phospholipids between different lipoprotein fractions in blood plasma. The inhibition of CETP has been shown to be a sound strategy to prevent and treat the development of coronary heart disease. We employed molecular dynamics simulations to unravel the mechanisms associated with the CETP-mediated lipid exchange. To this end we used both atomistic and coarse-grained models whose results were consistent with each other. We found CETP to bind to the surface of high density lipoprotein (HDL) -like lipid droplets through its charged and tryptophan residues. Upon binding, CETP rapidly (in about 10 ns) induced the formation of a small hydrophobic patch to the phospholipid surface of the droplet, opening a route from the core of the lipid droplet to the binding pocket of CETP. This was followed by a conformational change of helix X of CETP to an open state, in which we found the accessibility of cholesteryl esters to the C-terminal tunnel opening of CETP to increase. Furthermore, in the absence of helix X, cholesteryl esters rapidly diffused into CETP through the C-terminal opening. The results provide compelling evidence that helix X acts as a lid which conducts lipid exchange by alternating the open and closed states. The findings have potential for the design of novel molecular agents to inhibit the activity of CETP.

Coronary heart disease is a major cause of death in the Western societies. One of the most promising interventions to prevent and slow down the progress of coronary heart disease is the elevation of high density lipoprotein (HDL) levels in circulation. Animal models together with early clinical studies have shown that the inhibition of cholesteryl ester transfer protein (CETP) is a promising strategy to achieve higher HDL levels. However, drugs with acceptable side-effects for CETP-inhibition do not yet exist, although the next generation CETP inhibitor (anacetrapib) has great potential in this regard. In this study, our objective is to gain more detailed information regarding the interactions of CETP with lipoprotein particles. We show how the CETP-lipoprotein complex is formed and how lipid exchange between CETP and lipoprotein particles takes place. Our findings help to understand in a mechanistic way how CETP-mediated lipid exchange occurs and how it could be exploited in the design of new and more efficient molecular agents against coronary heart disease.

Lipid simulations: a perspective on lipids in action

In this article, we provide an overview of lipid simulations, describing how a computer can be used as a laboratory for lipid research. We briefly discuss the methodology of lipid simulations followed by a number of topical applications that show the benefit of computer modeling for complementing experiments. In particular, we show examples of cases in which simulations have made predictions of novel phenomena that have later been confirmed by experimental studies. Overall, the applications discussed in this article focus on the most recent state of the art and aim to provide a perspective of where the field of lipid simulations stands at the moment.

"Sticky" and "promiscuous", the yin and yang of apolipoprotein A-I termini in discoidal high-density lipoproteins: a combined computational-experimental approach

Apolipoprotein (apo) A-I-containing lipoproteins in the form of high-density lipoproteins (HDL) are inversely correlated with atherosclerosis. Because HDL is a soft form of condensed matter easily deformable by thermal fluctuations, the molecular mechanisms for HDL remodeling are not well understood. A promising approach to understanding HDL structure and dynamics is molecular dynamics (MD). In the present study, two computational strategies, MD simulated annealing (MDSA) and MD temperature jump, were combined with experimental particle reconstitution to explore molecular mechanisms for phospholipid- (PL-) rich HDL particle remodeling. The N-terminal domains of full-length apoA-I were shown to be "sticky", acting as a molecular latch largely driven by salt bridges, until, at a critical threshold of particle size, the associated domains released to expose extensive hydrocarbon regions of the PL to solvent. The "sticky" N-termini also associate with other apoA-I domains, perhaps being involved in N-terminal loops suggested by other laboratories. Alternatively, the overlapping helix 10 C-terminal domains of apoA-I were observed to be extremely mobile or "promiscuous", transiently exposing limited hydrocarbon regions of PL. Based upon these models and reconstitution studies, we propose that separation of the N-terminal domains, as particles exceed a critical size, triggers fusion between particles or between particles and membranes, while the C-terminal domains of apoA-I drive the exchange of polar lipids down concentration gradients between particles. This hypothesis has significant biological relevance since lipid exchange and particle remodeling are critically important processes during metabolism of HDL particles at every step in the antiatherogenic process of reverse cholesterol transport.

Molecular organization of the tear fluid lipid layer

The tear fluid protects the corneal epithelium from drying out as well as from invasion by pathogens. It also provides cell nutrients. Similarly to lung surfactant, it is composed of an aqueous phase covered by a lipid layer. Here we describe the molecular organization of the anterior lipid layer of the tear film. Artificial tear fluid lipid layers (ATFLLs) composed of egg yolk phosphatidylcholine (60 mol %), free fatty acids (20 mol %), cholesteryl oleate (10 mol %), and triglycerides (10 mol %) were deposited on the air-water interface and their physico-chemical behavior was compared to egg-yolk phosphatidylcholine monolayers by using Langmuir-film balance techniques, x-ray diffraction, and imaging techniques as well as in silico molecular level simulations. At low surface pressures, ATFLLs were organized at the air-water interface as heterogeneous monomolecular films. Upon compression the ATFLLs collapsed toward the air phase and formed hemispherelike lipid aggregates. This transition was reversible upon relaxation. These results were confirmed by molecular-level simulations of ATFLL, which further provided molecular-scale insight into the molecular distributions inside and dynamics of the tear film. Similar type of behavior is observed in lung surfactant but the folding takes place toward the aqueous phase. The results provide novel information of the function of lipids in the tear fluid.

Role of lipids in spheroidal high density lipoproteins

We study the structure and dynamics of spherical high density lipoprotein (HDL) particles through coarse-grained multi-microsecond molecular dynamics simulations. We simulate both a lipid droplet without the apolipoprotein A-I (apoA-I) and the full HDL particle including two apoA-I molecules surrounding the lipid compartment. The present models are the first ones among computational studies where the size and lipid composition of HDL are realistic, corresponding to human serum HDL. We focus on the role of lipids in HDL structure and dynamics. Particular attention is paid to the assembly of lipids and the influence of lipid-protein interactions on HDL properties. We find that the properties of lipids depend significantly on their location in the particle (core, intermediate region, surface). Unlike the hydrophobic core, the intermediate and surface regions are characterized by prominent conformational lipid order. Yet, not only the conformations but also the dynamics of lipids are found to be distinctly different in the different regions of HDL, highlighting the importance of dynamics in considering the functionalization of HDL. The structure of the lipid droplet close to the HDL-water interface is altered by the presence of apoA-Is, with most prominent changes being observed for cholesterol and polar lipids. For cholesterol, slow trafficking between the surface layer and the regimes underneath is observed. The lipid-protein interactions are strongest for cholesterol, in particular its interaction with hydrophobic residues of apoA-I. Our results reveal that not only hydrophobicity but also conformational entropy of the molecules are the driving forces in the formation of HDL structure. The results provide the first detailed structural model for HDL and its dynamics with and without apoA-I, and indicate how the interplay and competition between entropy and detailed interactions may be used in nanoparticle and drug design through self-assembly.