Molecular Dynamics Simulations Reveal Leaflet Coupling in Compositionally Asymmetric Phase-Separated Lipid Membranes

Modeling asymmetric cell membranes at all-atom resolution

Molecular dynamics (MD) simulations are a useful tool when studying the properties of membranes as they allow for a molecular view of lipid interactions with proteins, nucleic acids, or small molecules. While model membranes are usually symmetric in their lipid composition between leaflets and include a small number of lipid components, physiological membranes are highly complex and vary in the level of asymmetry. Simulation studies have shown that changes in leaflet asymmetry can alter the properties of a membrane. It is therefore necessary to carefully build asymmetric membranes to accurately simulate membranes. This chapter carefully describes the different methods for building asymmetric membranes and the advantages/disadvantages of each method. The simplest methods involve building a membrane with either an equal number of lipids per leaflet or an equal initial surface area (SA) estimated by the area per lipid. More detailed methods include combining two symmetric membranes of equal SA or altering an asymmetric membrane and adjusting the number of lipids after equilibration to minimize an observable such as differential stress (0-DS). More complex methods that require specific simulation software are also briefly described. The challenges and assumptions are listed for each method which should help guide the researcher to choose the best method for their unique MD simulation of an asymmetric membrane.

Atomistic Insight into the Lipid Nanodomains of Synaptic Vesicles

Membrane curvature, once regarded as a passive consequence of membrane composition and cellular architecture, has been shown to actively modulate various properties of the cellular membrane. These changes could also lead to segregation of the constituents of the membrane, generating nanodomains with precise biological properties. Proteins often linked with neurodegeneration (e.g., tau, alpha-synuclein) exhibit an unintuitive affinity for synaptic vesicles in neurons, which are reported to lack distinct, ordered nanodomains based on their composition. In this study, all-atom molecular dynamics simulations are used to study a full-scale synaptic vesicle of realistic Gaussian curvature and its effect on the membrane dynamics and lipid nanodomain organization. Compelling indicators of nanodomain formation, from the perspective of composition, surface areas per lipid, order parameter, and domain lifetime, are identified in the vesicle membrane, which are absent in a flat bilayer of the same lipid composition. Therefore, our study supports the idea that curvature may induce phase separation in an otherwise fluid, disordered membrane.

Building Asymmetric Lipid Bilayers for Molecular Dynamics Simulations: What Methods Exist and How to Choose One?

The compositional asymmetry of biological membranes has attracted significant attention over the last decade. Harboring more differences from symmetric membranes than previously appreciated, asymmetric bilayers have proven quite challenging to study with familiar concepts and techniques, leaving many unanswered questions about the reach of the asymmetry effects. One particular area of active research is the computational investigation of composition- and number-asymmetric lipid bilayers with molecular dynamics (MD) simulations. Offering a high level of detail into the organization and properties of the simulated systems, MD has emerged as an indispensable tool in the study of membrane asymmetry. However, the realization that results depend heavily on the protocol used for constructing the asymmetric bilayer models has sparked an ongoing debate about how to choose the most appropriate approach. Here we discuss the underlying source of the discrepant results and review the existing methods for creating asymmetric bilayers for MD simulations. Considering the available data, we argue that each method is well suited for specific applications and hence there is no single best approach. Instead, the choice of a construction protocol-and consequently, its perceived accuracy-must be based primarily on the scientific question that the simulations are designed to address.

Full scale structural, mechanical and dynamical properties of HIV-1 liposomes

Enveloped viruses are enclosed by a lipid membrane inside of which are all of the components necessary for the virus life cycle; viral proteins, the viral genome and metabolites. Viral envelopes are lipid bilayers that adopt morphologies ranging from spheres to tubes. The envelope is derived from the host cell during viral replication. Thus, the composition of the bilayer depends on the complex constitution of lipids from the host-cell's organelle(s) where assembly and/or budding of the viral particle occurs. Here, molecular dynamics (MD) simulations of authentic, asymmetric HIV-1 liposomes are used to derive a unique level of resolution of its full-scale structure, mechanics and dynamics. Analysis of the structural properties reveal the distribution of thicknesses of the bilayers over the entire liposome as well as its global fluctuations. Moreover, full-scale mechanical analyses are employed to derive the global bending rigidity of HIV-1 liposomes. Finally, dynamical properties of the lipid molecules reveal important relationships between their 3D diffusion, the location of lipid-rafts and the asymmetrical composition of the envelope. Overall, our simulations reveal complex relationships between the rich lipid composition of the HIV-1 liposome and its structural, mechanical and dynamical properties with critical consequences to different stages of HIV-1's life cycle.

Deciphering lipid transfer between and within membranes with time-resolved small-angle neutron scattering

This review focuses on time-resolved neutron scattering, particularly time-resolved small angle neutron scattering (TR-SANS), as a powerful in situ noninvasive technique to investigate intra- and intermembrane transport and distribution of lipids and sterols in lipid membranes. In contrast to using molecular analogues with potentially large chemical tags that can significantly alter transport properties, small angle neutron scattering relies on the relative amounts of the two most abundant isotope forms of hydrogen: protium and deuterium to detect complex membrane architectures and transport processes unambiguously. This review discusses advances in our understanding of the mechanisms that sustain lipid asymmetry in membranes-a key feature of the plasma membrane of cells-as well as the transport of lipids between membranes, which is an essential metabolic process.

Interdomain exchange and the flip-flop of cholesterol in ternary component lipid membranes and their effects on heterogeneous cholesterol diffusion

Cell membranes are heterogeneous with a variety of lipids, cholesterol, and proteins and are composed of domains of different compositions. Such heterogeneous environments make the transport of cholesterol complicated: cholesterol not only diffuses within a particular domain but also travels between domains. Cholesterol also flip-flops between upper and lower leaflets such that cholesterol may reside both within leaflets and in the central region between two leaflets. How the presence of multiple domains and the interdomain exchange of cholesterol would affect the cholesterol transport, however, remains elusive. In this study, therefore, we perform molecular dynamics simulations up to 100μs for ternary component lipid membranes, which consist of saturated lipids (dipalmitoylphosphatidylcholine, DPPC), unsaturated lipids (dilinoleylphosphatidylcholine, DIPC), and cholesterol. The ternary component membranes in our simulations form two domains readily: DPPC and DIPC domains. We find that the diffusion of cholesterol molecules is much more heterogeneous and non-Gaussian than expected for binary component lipid membranes of lipids and cholesterol. The non-Gaussian parameter of the cholesterol molecules is about four times larger in the ternary component lipid membranes than in the binary component lipid membranes. Such non-Gaussian and heterogeneous transport of cholesterol arises from the interplay among the interdomain kinetics, the different diffusivity of cholesterol in different domains, and the flip-flop of cholesterol. This suggests that in cell membranes that consist of various domains and proteins, the cholesterol transport can be very heterogeneous. We also find that the mechanism of the interdomain exchange differs for different domains: cholesterol tends to exit the DIPC domain along the central region of the membrane for the DIPC-to-DPPC transition, while the cholesterol is likely to exit the DPPC domain within the membrane leaflet for the DPPC-to-DIPC transition. Also, the interdomain exchange kinetics of cholesterol for the DPPC-to-DIPC transition is up to 7.9 times slower than the DIPC-to-DPPC transition.

LiPyphilic: A Python Toolkit for the Analysis of Lipid Membrane Simulations

Molecular dynamics simulations are now widely used to study emergent phenomena in lipid membranes with complex compositions. Here, we present LiPyphilic-a fast, fully tested, and easy-to-install Python package for analyzing such simulations. Analysis tools in LiPyphilic include the identification of cholesterol flip-flop events, the classification of local lipid environments, and the degree of interleaflet registration. LiPyphilic is both force field- and resolution-agnostic, and by using the powerful atom selection language of MDAnalysis, it can handle membranes with highly complex compositions. LiPyphilic also offers two on-the-fly trajectory transformations to (i) fix membranes split across periodic boundaries and (ii) perform nojump coordinate unwrapping. Our implementation of nojump unwrapping accounts for fluctuations in the box volume under the NPT ensemble-an issue that most current implementations have overlooked. The full documentation of LiPyphilic, including installation instructions and links to interactive online tutorials, is available at https://lipyphilic.readthedocs.io/en/latest.

Model Membrane Systems Used to Study Plasma Membrane Lipid Asymmetry

It is well known that the lipid distribution in the bilayer leaflets of mammalian plasma membranes (PMs) is not symmetric. Despite this, model membrane studies have largely relied on chemically symmetric model membranes for the study of lipid–lipid and lipid–protein interactions. This is primarily due to the difficulty in preparing stable, asymmetric model membranes that are amenable to biophysical studies. However, in the last 20 years, efforts have been made in producing more biologically faithful model membranes. Here, we review several recently developed experimental and computational techniques for the robust generation of asymmetric model membranes and highlight a new and particularly promising technique to study membrane asymmetry.

Lipid interactions of an actinoporin pore-forming oligomer

The actinoporins are cytolytic toxins produced by sea anemones. Upon encountering a membrane, preferably containing sphingomyelin, they oligomerize and insert their N-terminal helix into the membrane, forming a pore. Whether sphingomyelin is specifically recognized by the protein or simply induces phase coexistence in the membrane has been debated. Here, we perform multi-microsecond molecular dynamics simulations of an octamer of fragaceatoxin C, a member of the actinoporin family, in lipid bilayers containing either pure 1,2-dioleoyl-sn-glycero-3-phosphocholine (DOPC) or a 1:1 mixture of DOPC and palmitoyl sphingomyelin (PSM). The complex is highly stable in both environments, with only slight fraying of the inserted helices near their N-termini. Analyzing the structural parameters of the mixed membrane in the course of the simulation, we see signs of a phase transition for PSM in the inner leaflet of the bilayer. In both leaflets, cross-interactions between lipids of different type decrease over time. Surprisingly, the aromatic loop thought to be responsible for sphingomyelin recognition interacts more with DOPC than PSM by the end of the simulation. These results support the notion that the key membrane property that actinoporins recognize is lipid phase coexistence.

Bilayer compositional asymmetry influences the nanoscopic to macroscopic phase domain size transition

The eukaryotic plasma membrane (PM) exhibits lipid mixing heterogeneities known as lipid rafts. These lipid rafts, the result of liquid-liquid phase separation, can be modeled by coexisting liquid ordered (Lo) and liquid disordered (Ld) domains. Four-lipid component systems with a high-melting lipid, a nanodomain-inducing low-melting lipid, a macrodomain-inducing low-melting lipid, and cholesterol (chol) can give rise to domains of different sizes. These four-component systems have been characterized in experiments, yet there are few studies that model the asymmetric distribution of lipids actually found in the PM. We used molecular dynamics (MD) simulations to analyze the transition from nanoscopic to macroscopic domains in symmetric and in asymmetric model membranes. Using coarse-grained MD simulations, we found that asymmetry promotes macroscopic domain growth in a case where symmetric systems exhibit nanoscopic domains. Also, macroscopic domain formation in symmetric systems is highly dependent on registration of like phases in the cytoplasmic and exoplasmic leaflets. Using united-atom MD simulations, we found that symmetric Lo domains are only slightly more ordered than asymmetric Lo domains. We also found that large Lo domains in our asymmetric systems induce a slight chain ordering in the apposed cytoplasmic regions. The chol fractions of phase-separated Lo and Ld domains of the exoplasmic leaflet were unchanged whether the system was symmetric or asymmetric.

Creating Asymmetric Phospholipid Vesicles via Exchange With Lipid-Coated Silica Nanoparticles

Recently, effort has been placed into fabricating model free-floating asymmetric lipid membranes, such as asymmetric vesicles. Here, we report on the use of lipid-coated silica nanoparticles to exchange lipids with initially symmetric vesicles to generate composition-controlled asymmetric vesicles. Our method relies on the simple and natural exchange of lipids between membranes through an aqueous medium. Using a selected temperature, time, and ratio of lipid-coated silica nanoparticles to vesicles, we produced a desired highly asymmetric leaflet composition. At this point, the silica nanoparticles were removed by centrifugation, leaving the asymmetric vesicles in solution. In the present work, the asymmetric vesicles were composed of isotopically distinct dipalmitoylphosphatidylcholine lipids. Lipid asymmetry was detected by both small-angle neutron scattering (SANS) and proton nuclear magnetic resonance (¹H NMR). The rate at which the membrane homogenizes at 75 °C was also assessed.

Model Plasma Membrane Exhibits a Microemulsion in Both Leaves Providing a Foundation for "Rafts"

We consider a model lipid plasma membrane, one that describes the outer leaf as consisting of sphingomyelin, phosphatidylcholine, and cholesterol and the inner leaf of phosphatidylethanolamine, phosphatidylserine, phosphatidylcholine, and cholesterol. Their relative compositions are taken from experiment; the cholesterol freely interchanges between leaves. Fluctuations in local composition are coupled to fluctuations in the local membrane curvature, as in the Leibler-Andelman mechanism. Structure factors of components in both leaves display a peak at nonzero wavevector. This indicates that the disordered fluid membrane is characterized by structure of the corresponding wavelength. The scale is given by membrane properties: its bending modulus and its surface tension, which arises from the membrane's connections to the cytoskeleton. From measurements on the plasma membrane, this scale is on the order of 100 nm. We find that the membrane can be divided into two different kinds of domains that differ not only in their composition but also in their curvature. The first domain in the outer, exoplasmic leaf is rich in cholesterol and sphingomyelin, whereas the inner, cytoplasmic leaf is rich in phosphatidylserine and phosphatidylcholine. The second kind of domain is rich in phosphatidylcholine in the outer leaf and in cholesterol and phosphatidylethanolamine in the inner leaf. The theory provides a tenable basis for the origin of structure in the plasma membrane and an illuminating picture of the organization of lipids therein.