Multiscale modeling of lipids and lipid bilayers

A multiscale modeling approach is applied for simulations of lipids and lipid assemblies on mesoscale. First, molecular dynamics simulation of initially disordered system of lipid molecules in water within all-atomic model was carried out. On the next stage, structural data obtained from the molecular dynamics (MD) simulation were used to build a coarse-grained (ten sites) lipid model, with effective interaction potentials computed by the inverse Monte Carlo method. Finally, several simulations of the coarse-grained model on longer length- and time-scale were performed, both within Monte Carlo and molecular dynamics simulations: a periodical sample of lipid molecules ordered in bilayer, a free sheet of such bilayer without periodic boundary conditions, formation of vesicle from a plain membrane, process of self-assembly of lipids randomly dispersed in volume. It was shown that the coarse-grained model, developed exclusively from all-atomic simulation data, reproduces well all the basic features of lipids in water solution.     

Mutations and Lethality in Simulated Prebiotic Networks

The Graded Autocatalysis Replication Domain (GARD) model describes an origin of life scenario which involves non-covalent compositional assemblies, made of monomeric mutually catalytic molecules. GARD constitutes an alternative to informational biopolymers as a mechanism of primordial inheritance. In the present work, we examined the effect of mutations, one of the most fundamental mechanisms for evolution, in the context of the networks of mutual interaction within GARD prebiotic assemblies. We performed a systematic analysis analogous to single and double gene deletions within GARD. While most deletions have only a small effect on both growth rate and molecular composition of the assemblies, ~10% of the deletions caused lethality, or sometimes showed enhanced fitness. Analysis of 14 different network properties on 2,000 different GARD networks indicated that lethality usually takes place when the deleted node has a high molecular count, or when it is a catalyst for such node. A correlation was also found between lethality and node degree centrality, similar to what is seen in real biological networks. Addressing double knockout mutations, our results demonstrate the occurrence of both synthetic lethality and extragenic suppression within GARD networks, and convey an attempt to correlate synthetic lethality to network node-pair properties. The analyses presented help establish GARD as a workable alternative prebiotic scenario, suggesting that life may have begun with large molecular networks of low fidelity, that later underwent evolutionary compaction and fidelity augmentation.     

Modulation of entry of enveloped viruses by cholesterol and sphingolipids (Review)

Enveloped animal viruses infect host cells by fusion of viral and target membranes. This crucial fusion event occurs either with the plasma membrane of the host cells at the physiological pH or with the endosomal membranes at low pH and is triggered by specific glycoproteins in the virus envelope. Both lipids and proteins play critical and co-operative roles in the fusion process. Interactions of viral proteins with their receptors direct which membranes fuse and viral fusion proteins then drive the process. These fusion proteins operate on lipid assemblies, whose physical and mechanical properties are equally important to the proper functioning of the process. Lipids contribute to the viral fusion process by virtue of their distinct chemical structure, composition and/or their preferred partitioning into specific microdomains in the plasma membrane called 'rafts'. An involvement of lipid rafts in viral entry and membrane fusion has been examined recently. However, the mechanism(s) by which lipids as dynamic raft components control viral envelope-glycoprotein-triggered fusion is not clear. This paper will review literature findings on the contribution of the two raft-associated lipids, cholesterol and sphingolipids in viral entry.     

The molecular roots of compositional inheritance.

Non-covalent compositional assemblies, made of monomeric mutually catalytic molecules, constitute an alternative to alphabet-based informational biopolymers as a mechanism of primordial inheritance. Such assemblies appear implicitly in many "Metabolism First" origin of life scenarios, and more explicitly in the Graded Autocatalysis Replication Domain (GARD) model [Segréet al. (2000). Proc. Natl Acad. Sci. U.S.A.97, 4112-4117]. In the present work, we provide a detailed analysis of the quantitative molecular roots of such behavior. It is demonstrated that the fidelity of reproduction provided by a newly defined heritability measure eta(*)(s), strongly depends on the values of molecular recognition parameters and on assembly size. We find that if the catalytic rate acceleration coefficients are distributed normally, transfer of compositional information becomes impossible, due to frequent "compositional error catastrophes". In contrast, if the catalytic acceleration rates obey a lognormal distribution, as actually predicted by a statistical formalism for molecular repertoires, high reproduction fidelity is obtained. There is also a clear dependence on assembly size N, whereby maximal eta is seen in a narrow range around N approximately 3.5 N(G)/lambda, where N(G)is the size of the primordial molecular repertoire and lambda is a molecular interaction statistical parameter. Such relationships help define the physicochemical conditions that could underlie the early steps in pre-biotic evolution.     

ABC lipid transporters: extruders, flippases, or flopless activators?

Many mammalian ABC transporters move membrane lipids to acceptor lipid assemblies in the extracellular aqueous milieu. Because the desorption from the membrane costs more energy than provided by two ATPs, the transporter probably only translocates the lipid to a partially hydrophilic site on its extracellular face. From this high-energy site, the lipid may efficiently move to the acceptor, which ideally is bound to the transporter, or, in the absence of an acceptor, fall back into the membrane. If the lipid originated from the cytosolic membrane surface, this represents lipid flop and is probably a side activity of the transporters.     

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.     

Lipids and lipid domains in the peroxisomal membrane of the yeast Yarrowia lipolytica

Biological membranes have unique and highly diverse compositions of their lipid constituents. At present, we have only partial understanding of how membrane lipids and lipid domains regulate the structural integrity and functionality of cellular organelles, maintain the unique molecular composition of each organellar membrane by orchestrating the intracellular trafficking of membrane-bound proteins and lipids, and control the steady-state levels of numerous signaling molecules generated in biological membranes. Similar to other organellar membranes, a single lipid bilayer enclosing the peroxisome, an organelle known for its essential role in lipid metabolism, has a unique lipid composition and organizes some of its lipid and protein components into distinctive assemblies. This review highlights recent advances in our knowledge of how lipids and lipid domains of the peroxisomal membrane regulate the processes of peroxisome assembly and maintenance in the yeast Yarrowia lipolytica. We critically evaluate the molecular mechanisms through which lipid constituents of the peroxisomal membrane control these multistep processes and outline directions for future research in this field.     

Adaptive Lipid Packing and Bioactivity in Membrane Domains

Lateral compositional and physicochemical heterogeneity is a ubiquitous feature of cellular membranes on various length scales, from molecular assemblies to micrometric domains. Segregated lipid domains of increased local order, referred to as rafts, are believed to be prominent features in eukaryotic plasma membranes; however, their exact nature (i.e. size, lifetime, composition, homogeneity) in live cells remains difficult to define. Here we present evidence that both synthetic and natural plasma membranes assume a wide range of lipid packing states with varying levels of molecular order. These states may be adapted and specifically tuned by cells during active cellular processes, as we show for stimulated insulin secretion. Most importantly, these states regulate both the partitioning of molecules between coexisting domains and the bioactivity of their constituent molecules, which we demonstrate for the ligand binding activity of the glycosphingolipid receptor GM1. These results confirm the complexity and flexibility of lipid-mediated membrane organization and reveal mechanisms by which this flexibility could be functionalized by cells.     

Metabolic changes of membrane lipid composition in Acholeplasma laidlawii by hydrocarbons, alcohols, and detergents: arguments for effects on lipid packing

The packing of lipids into different aggregates, such as spheres, rods, or bilayers, is dependent on the hydrophobic volume, the hydrocarbon-water interfacial area, and the hydrocarbon chain length of the participating molecules, according to the self-assembly theory [Israelachvili, J. N., Marcelja, S., & Horn, R. G. (1980) Q. Rev. Biophys. 13, 121-200]. The origin of the participating molecules should be of no importance with respect to their abilities to affect the above-mentioned parameters. In this investigation, Acholeplasma laidlawii, with a defined acyl chain composition of the membrane lipids, has been grown in the presence of three different classes of foreign molecules, known to partition into model and biological membranes. This results in an extensive metabolic alteration in the lipid polar head group composition, which is expressed as changes in the molar ratio between the lipids monoglucosyldiglyceride (MGDG) and diglucosyldiglyceride (DGDG), forming reversed hexagonal and lamellar phases in excess water, respectively. The formation of nonlamellar phases by A. laidlawii lipids depends critically upon the MGDG concentration [Lindblom, G., Brentel, I., Sj?lund, M., Wikander, G., & Wieslander, A. (1986) Biochemistry (preceding paper in this issue)]. The foreign molecules tested belong to the following groups: nonpolar organic solvents, alcohols, and detergents. Their effects on the gel to liquid crystalline phase transition temperature (Tm), on the order parameter of the acyl chains, and on the phase equilibria between lamellar and nonlamellar liquid crystalline phases in lipid-water model systems are known in several instances.(ABSTRACT TRUNCATED AT 250 WORDS)     

Affinity purification of lipid vesicles

We present a novel column chromatography technique for recovery and purification of lipid vesicles, which can be extended to other macromolecular assemblies. This technique is based on reversible binding of biotinylated lipids to monomeric avidin. Unlike the very strong binding of biotin and biotin-functionalized molecules to streptavidin, the interaction between biotin-functionalized molecules and monomeric avidin can be disrupted effectively by ligand competition from free biotin. In this work, biotin-functionalized lipids (biotin-PEG-PE) were incorporated into synthetic lipid vesicles (DOPC), resulting in unilamellar biotinylated lipid vesicles. The vesicles were bound to immobilized monomeric avidin, washed extensively with buffer, and eluted with a buffer supplemented with free biotin. Increasing the biotinyl lipid molar ratio beyond 0.53% of all lipids did not increase the efficiency of vesicle recovery. A simple adsorption model suggests 1.1 x 10(13) active binding sites/mL of resin with an equilibrium binding constant of K = 1.0 x 10(8) M(-1). We also show that this method is very robust and reproducible and can accommodate vesicles of varying sizes with diverse contents. This method can be scaled up to larger columns and/or high throughput analysis, such as a 96-well plate format.     

NERDSS: A Nonequilibrium Simulator for Multibody Self-Assembly at the Cellular Scale

Currently, a significant barrier to building predictive models of cellular self-assembly processes is that molecular models cannot capture minutes-long dynamics that couple distinct components with active processes, whereas reaction-diffusion models cannot capture structures of molecular assembly. Here, we introduce the nonequilibrium reaction-diffusion self-assembly simulator (NERDSS), which addresses this spatiotemporal resolution gap. NERDSS integrates efficient reaction-diffusion algorithms into generalized software that operates on user-defined molecules through diffusion, binding and orientation, unbinding, chemical transformations, and spatial localization. By connecting the fast processes of binding with the slow timescales of large-scale assembly, NERDSS integrates molecular resolution with reversible formation of ordered, multisubunit complexes. NERDSS encodes models using rule-based formatting languages to facilitate model portability, usability, and reproducibility. Applying NERDSS to steps in clathrin-mediated endocytosis, we design multicomponent systems that can form lattices in solution or on the membrane, and we predict how stochastic but localized dephosphorylation of membrane lipids can drive lattice disassembly. The NERDSS simulations reveal the spatial constraints on lattice growth and the role of membrane localization and cooperativity in nucleating assembly. By modeling viral lattice assembly and recapitulating oscillations in protein expression levels for a circadian clock model, we illustrate the adaptability of NERDSS. NERDSS simulates user-defined assembly models that were previously inaccessible to existing software tools, with broad applications to predicting self-assembly in vivo and designing high-yield assemblies in vitro.     

Domain coupling in asymmetric lipid bilayers

Biological membranes are heterogeneous assemblies of lipids, proteins, and cholesterol that are organized as asymmetric bimolecular leaflets of lipids with embedded proteins. Modulated by the concentration of cholesterol lipids and proteins may segregate into two or more liquid phases with different physical properties that can coexist in the same membrane. In this review, we summarize recent advances on how this situation can be recreated in a supported bilayer format and how this system has been used to demonstrate the induction of ordered lipid domains in lipid compositions that are typical for the inner leaflet by lipid compositions that are typical for the outer leaflet of mammalian plasma membranes. Proteins are shown to differentially target such induced inner leaflet domains.     

Preparation of giant unilamellar vesicles from damp lipid film for better lipid compositional uniformity

Giant unilamellar vesicles (GUVs) containing cholesterol often have a wide distribution in lipid composition. In this study, GUVs of 1,2-dioleoyl-sn-glycero-3-phosphocholine(DOPC)/1,2-distearoyl-sn-glycero-3-phosphocholine(DSPC)/cholesterol and 1,2-diphytanoyl-sn-glycero-3-phosphocholine(diPhyPC)/1,2-dipalmitoyl-sn-glycero-3-phosphocholine(DPPC)/cholesterol were prepared from dry lipid films using the standard electroformation method as well as a modified method from damp lipid films, which are made from compositional uniform liposomes prepared using the Rapid Solvent Exchange (RSE) method. We quantified the lipid compositional distributions of GUV by measuring the miscibility transition temperature of GUVs using fluorescence microscopy, since a narrower distribution in the transition temperature should correspond to a more uniform distribution in GUV lipid composition. Cholesterol molecules can demix from other lipids in dry state and form cholesterol crystals. Using optical microscopy, micron-sized crystals were observed in some dry lipid films. Thus, a major cause of GUV lipid compositional heterogeneity is the demixing of lipids in the dry film state. By avoiding the dry film state, GUVs prepared from damp lipid films have a better uniformity in lipid composition, and the standard deviations of miscibility transition temperature are about 2.5 times smaller than that of GUVs prepared from dry lipid films. Comparing the two ternary systems, diPhyPC/DPPC/cholesterol GUVs has a larger cholesterol compositional heterogeneity, which directly correlates with the low maximum solubility of cholesterol in diPhyPC lipid bilayers (40.2±0.5mol%) measured by light scattering. Our data indicate that cholesterol interacts far less favorably with diPhyPC than it does with other PCs. The damp lipid film method also has a potential of preparing GUVs from cell membranes containing native proteins without going through a dry state.     

Mesoscale organization of domains in the plasma membrane – beyond the lipid raft

The plasma membrane is compartmentalized into several distinct regions or domains, which show a broad diversity in both size and lifetime. The segregation of lipids and membrane proteins is thought to be driven by the lipid composition itself, lipid–protein interactions and diffusional barriers. With regards to the lipid composition, the immiscibility of certain classes of lipids underlies the “lipid raft” concept of plasmalemmal compartmentalization. Historically, lipid rafts have been described as cholesterol and (glyco)sphingolipid-rich regions of the plasma membrane that exist as a liquid-ordered phase that are resistant to extraction with non-ionic detergents. Over the years the interest in lipid rafts grew as did the challenges with studying these nanodomains. The term lipid raft has fallen out of favor with many scientists and instead the terms “membrane raft” or “membrane nanodomain” are preferred as they connote the heterogeneity and dynamic nature of the lipid-protein assemblies. In this article, we will discuss the classical lipid raft hypothesis and its limitations. This review will also discuss alternative models of lipid-protein interactions, annular lipid shells, and larger membrane clusters. We will also discuss the mesoscale organization of plasmalemmal domains including visible structures such as clathrin-coated pits and caveolae.     

Yeast Lipids Can Phase-separate into Micrometer-scale Membrane Domains

The lipid raft concept proposes that biological membranes have the potential to form functional domains based on a selective interaction between sphingolipids and sterols. These domains seem to be involved in signal transduction and vesicular sorting of proteins and lipids. Although there is biochemical evidence for lipid raft-dependent protein and lipid sorting in the yeast Saccharomyces cerevisiae, direct evidence for an interaction between yeast sphingolipids and the yeast sterol ergosterol, resulting in membrane domain formation, is lacking. Here we show that model membranes formed from yeast total lipid extracts possess an inherent self-organization potential resulting in liquid-disordered-liquid-ordered phase coexistence at physiologically relevant temperature. Analyses of lipid extracts from mutants defective in sphingolipid metabolism as well as reconstitution of purified yeast lipids in model membranes of defined composition suggest that membrane domain formation depends on specific interactions between yeast sphingolipids and ergosterol. Taken together, these results provide a mechanistic explanation for lipid raft-dependent lipid and protein sorting in yeast.     

Chiral DNA packaging in DNA-cationic liposome assemblies.

Recent studies have indicated that the structural features of DNA-lipid assemblies, dictated by the lipid composition and cationic lipid-to-DNA ratio, critically affect the efficiency of these complexes in acting as vehicles for cellular delivery of genetic material. Using circular dichroism we find that upon binding DNA, positively-charged liposomes induce a secondary conformational transition of the DNA molecules from the native B form to the C motif. Liposomes composed of positively-charged and neutral 'helper' lipids, found to be particularly effective as transfecting agents, induce - in addition to secondary conformational changes - DNA condensation into a left-handed cholesteric-like phase. A structural model is presented according to which two distinct, yet inter-related modes of DNA packaging coexist within such assemblies. The results underline the notion that subtle changes in the components of a supramolecular assembly may substantially modulate the interplay of interactions which dictate its structure and functional properties.     

The lipid network

Natural cell membranes are composed of a remarkable variety of lipids, which provide specific biophysical properties to support membrane protein function. An improved understanding of this complexity of membrane composition may also allow the design of membrane active drugs. Crafting a relevant model of a cell membrane with controlled composition is becoming an art, with the ability to reveal the molecular mechanisms of biological processes and lead to better treatment of pathologies. By matching physiological observations from in vivo experiments to high-resolution information, more easily obtained from in vitro studies, complex interactions at the lipid interface are determined. The role of the lipid network in biological membranes is, therefore, the subject of increasing attention.     

Domain-induced activation of human phospholipase A2 type IIA: local versus global lipid composition

Secretory human phospholipase A2 type IIA (PLA2-IIA) catalyzes the hydrolysis of the sn-2 ester bond in glycerolipids to produce fatty acids and lysolipids. The enzyme is coupled to the inflammatory response, and its specificity toward anionic membrane interfaces suggests a role as a bactericidal agent. PLA2-IIA may also target perturbed native cell membranes that expose anionic lipids to the extracellular face. However, anionic lipid contents in native cells appear lower than the threshold levels necessary for activation. By using phosphatidylcholine/phosphatidylglycerol model systems, we show that local enrichment of anionic lipids into fluid domains triggers PLA2-IIA activity. In addition, the compositional range of enzyme activity is shown to be related to the underlying lipid phase diagram. A comparison is done between PLA2-IIA and snake venom PLA2, which in contrast to PLA2-IIA hydrolyzes both anionic and zwitterionic membranes. In general, this work shows that PLA2-IIA activation can be accomplished through local enrichment of anionic lipids into domains, indicating a mechanism for PLA2-IIA to target perturbed native membranes with low global anionic lipid contents. The results also show that the underlying lipid phase diagram, which determines the lipid composition at a local level, can be used to predict PLA2-IIA activity.     

Lipid Clustering Correlates with Membrane Curvature as Revealed by Molecular Simulations of Complex Lipid Bilayers

Cell membranes are complex multicomponent systems, which are highly heterogeneous in the lipid distribution and composition. To date, most molecular simulations have focussed on relatively simple lipid compositions, helping to inform our understanding of in vitro experimental studies. Here we describe on simulations of complex asymmetric plasma membrane model, which contains seven different lipids species including the glycolipid GM3 in the outer leaflet and the anionic lipid, phosphatidylinositol 4,5-bisphophate (PIP₂), in the inner leaflet. Plasma membrane models consisting of 1500 lipids and resembling the in vivo composition were constructed and simulations were run for 5 µs. In these simulations the most striking feature was the formation of nano-clusters of GM3 within the outer leaflet. In simulations of protein interactions within a plasma membrane model, GM3, PIP₂, and cholesterol all formed favorable interactions with the model α-helical protein. A larger scale simulation of a model plasma membrane containing 6000 lipid molecules revealed correlations between curvature of the bilayer surface and clustering of lipid molecules. In particular, the concave (when viewed from the extracellular side) regions of the bilayer surface were locally enriched in GM3. In summary, these simulations explore the nanoscale dynamics of model bilayers which mimic the in vivo lipid composition of mammalian plasma membranes, revealing emergent nanoscale membrane organization which may be coupled both to fluctuations in local membrane geometry and to interactions with proteins.