Cryo-EM images of phase-separated lipid bilayer vesicles analyzed with a machine-learning approach.

Lateral lipid heterogeneity (i.e., raft formation) in biomembranes plays a functional role in living cells. Three-component mixtures of low- and high-melting lipids plus cholesterol offer a simplified experimental model for raft domains in which a liquid-disordered (Ld) phase coexists with a liquid-ordered (Lo) phase. Using such models, we recently showed that cryogenic electron microscopy (cryo-EM) can detect phase separation in lipid vesicles based on differences in bilayer thickness. However, the considerable noise within cryo-EM data poses a significant challenge for accurately determining the membrane phase state at high spatial resolution. To this end, we have developed an image-processing pipeline that utilizes machine learning (ML) to predict the bilayer phase in projection images of lipid vesicles. Importantly, the ML method exploits differences in both the thickness and molecular density of Lo compared to Ld, which leads to improved phase identification. To assess accuracy, we used artificial images of phase-separated lipid vesicles generated from all-atom molecular dynamics simulations of Lo and Ld phases. Synthetic ground-truth data sets mimicking a series of compositions along a tieline of Ld + Lo coexistence were created and then analyzed with various ML models. For all tieline compositions, we find that the ML approach can correctly identify the bilayer phase with >90% accuracy, thus providing a means to isolate the intensity profiles of coexisting Ld and Lo phases, as well as accurately determine domain-size distributions, number of domains, and phase-area fractions. The method described here provides a framework for characterizing nanoscopic lateral heterogeneities in membranes and paves the way for a more detailed understanding of raft properties in biological contexts.

Molecular simulations and NMR reveal how lipid fluctuations affect membrane mechanics.

Lipid bilayers form the main matrix of functional cell membranes, and their dynamics underlie a host of physical and biological processes. Here we show that elastic membrane properties and collective molecular dynamics (MD) are related by the mean-square amplitudes (order parameters) and relaxation rates (correlation times) of lipid acyl chain motions. We performed all-atom MD simulations of liquid-crystalline bilayers that allow direct comparison with carbon-hydrogen (CH) bond relaxations measured with NMR spectroscopy. Previous computational and theoretical approaches have assumed isotropic relaxation, which yields inaccurate description of lipid chain dynamics and incorrect data interpretation. Instead, the new framework includes a fixed bilayer normal (director axis) and restricted anisotropic motion of the CH bonds in accord with their segmental order parameters, enabling robust validation of lipid force fields. Simulated spectral densities of thermally excited CH bond fluctuations exhibited well-defined spin-lattice (Zeeman) relaxations analogous to those in NMR measurements. Their frequency signature could be fit to a simple power-law function, indicative of nematic-like collective dynamics. Moreover, calculated relaxation rates scaled as the squared order parameters yielding an apparent κC modulus for bilayer bending. Our results show a strong correlation with κC values obtained from solid-state NMR studies of bilayers without and with cholesterol as validated by neutron spin-echo measurements of membrane elasticity. The simulations uncover a critical role of interleaflet coupling in membrane mechanics and thus provide important insights into molecular sites of emerging elastic properties within lipid bilayers.

Direct label-free imaging of nanodomains in biomimetic and biological membranes by cryogenic electron microscopy.

The nanoscale organization of biological membranes into structurally and compositionally distinct lateral domains is believed to be central to membrane function. The nature of this organization has remained elusive due to a lack of methods to directly probe nanoscopic membrane features. We show here that cryogenic electron microscopy (cryo-EM) can be used to directly image coexisting nanoscopic domains in synthetic and bioderived membranes without extrinsic probes. Analyzing a series of single-component liposomes composed of synthetic lipids of varying chain lengths, we demonstrate that cryo-EM can distinguish bilayer thickness differences as small as 0.5 Å, comparable to the resolution of small-angle scattering methods. Simulated images from computational models reveal that features in cryo-EM images result from a complex interplay between the atomic distribution normal to the plane of the bilayer and imaging parameters. Simulations of phase-separated bilayers were used to predict two sources of contrast between coexisting ordered and disordered phases within a single liposome, namely differences in membrane thickness and molecular density. We observe both sources of contrast in biomimetic membranes composed of saturated lipids, unsaturated lipids, and cholesterol. When extended to isolated mammalian plasma membranes, cryo-EM reveals similar nanoscale lateral heterogeneities. The methods reported here for direct, probe-free imaging of nanodomains in unperturbed membranes open new avenues for investigation of nanoscopic membrane organization.

How cholesterol stiffens unsaturated lipid membranes.

Cholesterol is an integral component of eukaryotic cell membranes and a key molecule in controlling membrane fluidity, organization, and other physicochemical parameters. It also plays a regulatory function in antibiotic drug resistance and the immune response of cells against viruses, by stabilizing the membrane against structural damage. While it is well understood that, structurally, cholesterol exhibits a densification effect on fluid lipid membranes, its effects on membrane bending rigidity are assumed to be nonuniversal; i.e., cholesterol stiffens saturated lipid membranes, but has no stiffening effect on membranes populated by unsaturated lipids, such as 1,2-dioleoyl-sn-glycero-3-phosphocholine (DOPC). This observation presents a clear challenge to structure-property relationships and to our understanding of cholesterol-mediated biological functions. Here, using a comprehensive approach-combining neutron spin-echo (NSE) spectroscopy, solid-state deuterium NMR (2H NMR) spectroscopy, and molecular dynamics (MD) simulations-we report that cholesterol locally increases the bending rigidity of DOPC membranes, similar to saturated membranes, by increasing the bilayer's packing density. All three techniques, inherently sensitive to mesoscale bending fluctuations, show up to a threefold increase in effective bending rigidity with increasing cholesterol content approaching a mole fraction of 50%. Our observations are in good agreement with the known effects of cholesterol on the area-compressibility modulus and membrane structure, reaffirming membrane structure-property relationships. The current findings point to a scale-dependent manifestation of membrane properties, highlighting the need to reassess cholesterol's role in controlling membrane bending rigidity over mesoscopic length and time scales of important biological functions, such as viral budding and lipid-protein interactions.

Structural and functional consequences of reversible lipid asymmetry in living membranes.

Maintenance of lipid asymmetry across the two leaflets of the plasma membrane (PM) bilayer is a ubiquitous feature of eukaryotic cells. Loss of this asymmetry has been widely associated with cell death. However, increasing evidence points to the physiological importance of non-apoptotic, transient changes in PM asymmetry. Such transient scrambling events are associated with a range of biological functions, including intercellular communication and intracellular signaling. Thus, regulation of interleaflet lipid distribution in the PM is a broadly important but underappreciated cellular process with key physiological and structural consequences. Here, we compile the mounting evidence revealing multifaceted, functional roles of PM asymmetry and transient loss thereof. We discuss the consequences of reversible asymmetry on PM structure, biophysical properties and interleaflet coupling. We argue that despite widespread recognition of broad aspects of membrane asymmetry, its importance in cell biology demands more in-depth investigation of its features, regulation, and physiological and pathological implications.

Probing the Link between Pancratistatin and Mitochondrial Apoptosis through Changes in the Membrane Dynamics on the Nanoscale.

Pancratistatin (PST) is a natural antiviral alkaloid that has demonstrated specificity toward cancerous cells and explicitly targets the mitochondria. PST initiates apoptosis while leaving healthy, noncancerous cells unscathed. However, the manner by which PST induces apoptosis remains elusive and impedes the advancement of PST as a natural anticancer therapeutic agent. Herein, we use neutron spin-echo (NSE) spectroscopy, molecular dynamics (MD) simulations, and supporting small angle scattering techniques to study PST's effect on membrane dynamics using biologically representative model membranes. Our data suggests that PST stiffens the inner mitochondrial membrane (IMM) by being preferentially associated with cardiolipin, which would lead to the relocation and release of cytochrome c. Second, PST has an ordering effect on the lipids and disrupts their distribution within the IMM, which would interfere with the maintenance and functionality of the active forms of proteins in the electron transport chain. These previously unreported findings implicate PST's effect on mitochondrial apoptosis.

Ca2+-dependent mechanism of membrane insertion and destabilization by the SARS-CoV-2 fusion peptide.

Cell penetration after recognition of the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) virus by the ACE2 receptor and the fusion of its viral envelope membrane with cellular membranes are the early steps of infectivity. A region of the Spike protein of the virus, identified as the "fusion peptide" (FP), is liberated at its N-terminal site by a specific cleavage occurring in concert with the interaction of the receptor-binding domain of the Spike. Studies have shown that penetration is enhanced by the required binding of Ca2+ ions to the FPs of coronaviruses, but the mechanisms of membrane insertion and destabilization remain unclear. We have predicted the preferred positions of Ca2+ binding to the SARS-CoV-2-FP, the role of Ca2+ ions in mediating peptide-membrane interactions, the preferred mode of insertion of the Ca2+-bound SARS-CoV-2-FP, and consequent effects on the lipid bilayer from extensive atomistic molecular dynamics simulations and trajectory analyses. In a systematic sampling of the interactions of the Ca2+-bound peptide models with lipid membranes, SARS-CoV-2-FP penetrated the bilayer and disrupted its organization only in two modes involving different structural domains. In one, the hydrophobic residues F833/I834 from the middle region of the peptide are inserted. In the other, more prevalent mode, the penetration involves residues L822/F823 from the LLF motif, which is conserved in CoV-2-like viruses, and is achieved by the binding of Ca2+ ions to the D830/D839 and E819/D820 residue pairs. FP penetration is shown to modify the molecular organization in specific areas of the bilayer, and the extent of membrane binding of the SARS-CoV-2 FP is significantly reduced in the absence of Ca2+ ions. These findings provide novel mechanistic insights regarding the role of Ca2+ in mediating SARS-CoV-2 fusion and provide a detailed structural platform to aid the ongoing efforts in rational design of compounds to inhibit SARS-CoV-2 cell entry.

Vesicle Viewer: Online visualization and analysis of small-angle scattering from lipid vesicles.

Small-angle X-ray and neutron scattering are among the most powerful experimental techniques for investigating the structure of biological membranes. Much of the critical information contained in small-angle scattering (SAS) data is not easily accessible to researchers who have limited time to analyze results by hand or to nonexperts who may lack the necessary scientific background to process such data. Easy-to-use data visualization software can allow them to take full advantage of their SAS data and maximize the use of limited resources. To this end, we developed an internet-based application called Vesicle Viewer to visualize and analyze SAS data from unilamellar lipid bilayer vesicles. Vesicle Viewer utilizes a modified scattering density profile (SDP) analysis called EZ-SDP in which key bilayer structural parameters, such as area per lipid and bilayer thickness, are easily and robustly determined. Notably, we introduce a bilayer model that is able to describe an asymmetric bilayer, whether it be chemically or isotopically asymmetric. The application primarily uses Django, a Python package specialized for the development of robust web applications. In addition, several other libraries are used to support the more technical aspects of the project; notable examples are Matplotlib (for graphs) and NumPy (for calculations). By eliminating the barrier of downloading and installing software, this web-based application will allow scientists to analyze their own vesicle scattering data using their preferred operating system. The web-based application can be found at https://vesicleviewer.dmarquardt.ca/.

Membrane lipids are both the substrates and a mechanistically responsive environment of TMEM16 scramblase proteins.

Recent discoveries about functional mechanisms of proteins in the TMEM16 family of phospholipid scramblases have illuminated the dual role of the membrane as both the substrate and a mechanistically responsive environment in the wide range of physiological processes and genetic disorders in which they are implicated. This is highlighted in the review of recent findings from our collaborative investigations of molecular mechanisms of TMEM16 scramblases that emerged from iterative functional, structural, and computational experimentation. In the context of this review, we present new MD simulations and trajectory analyses motivated by the fact that new structural information about the TMEM16 scramblases is emerging from cryo-EM determinations in lipid nanodiscs. Because the functional environment of these proteins in in vivo and in in vitro is closer to flat membranes, we studied comparatively the responses of the membrane to the TMEM16 proteins in flat membranes and nanodiscs. We find that bilayer shapes in the nanodiscs are very different from those observed in the flat membrane systems, but the function-related slanting of the membrane observed at the nhTMEM16 boundary with the protein is similar in the nanodiscs and in the flat bilayers. This changes, however, in the bilayer composed of longer-tail lipids, which is thicker near the phospholipid translocation pathway, which may reflect an enhanced tendency of the long tails to penetrate the pathway and create, as shown previously, a nonconductive environment. These findings support the correspondence between the mechanistic involvement of the lipid environment in the flat membranes, and the nanodiscs. © 2019 Wiley Periodicals, Inc.

A New Computational Method for Membrane Compressibility: Bilayer Mechanical Thickness Revisited.

Because lipid bilayers can bend and stretch in ways similar to thin elastic sheets, physical models of bilayer deformation have utilized mechanical constants such as the moduli for bending rigidity (κC) and area compressibility (KA). However, the use of these models to quantify the energetics of membrane deformation associated with protein-membrane interactions, and the membrane response to stress is often hampered by the shortage of experimental data suitable for the estimation of the mechanical constants of various lipid mixtures. Although computational tools such as molecular dynamics simulations can provide alternative means to estimate KA values, current approaches suffer significant technical limitations. Here, we present a novel, to our knowledge, computational framework that allows for a direct estimation of KA values for individual bilayer leaflets. The theory is based on the concept of elasticity and derives KA from real-space analysis of local thickness fluctuations sampled in molecular dynamics simulations. We explore and validate the model on a large set of single and multicomponent bilayers of different lipid compositions and sizes, simulated at different temperatures. The calculated bilayer compressibility moduli agree with values estimated previously from experiments and those obtained from a standard computational method based on a series of constrained tension simulations. We further validate our framework in a comparison with an existing polymer brush model and confirm the polymer brush model's predicted linear relationship with proportionality coefficient of 24, using elastic parameters calculated from the simulation trajectories. The robustness of the results that emerge from the method allows us to revisit the origins of the bilayer mechanical (compressible) thickness and in particular its dependence on acyl-chain unsaturation and the presence of cholesterol.

Gramicidin Increases Lipid Flip-Flop in Symmetric and Asymmetric Lipid Vesicles.

Unlike most transmembrane proteins, phospholipids can migrate from one leaflet of the membrane to the other. Because this spontaneous lipid translocation (flip-flop) tends to be very slow, cells facilitate the process with enzymes that catalyze the transmembrane movement and thereby regulate the transbilayer lipid distribution. Nonenzymatic membrane-spanning proteins with unrelated primary functions have also been found to accelerate lipid flip-flop in a nonspecific manner and by various hypothesized mechanisms. Using deuterated phospholipids, we examined the acceleration of flip-flop by gramicidin channels, which have well-defined structures and known functions, features that make them ideal candidates for probing the protein-membrane interactions underlying lipid flip-flop. To study compositionally and isotopically asymmetric proteoliposomes containing gramicidin, we expanded a recently developed protocol for the preparation and characterization of lipid-only asymmetric vesicles. Channel incorporation, conformation, and function were examined with small angle x-ray scattering, circular dichroism, and a stopped-flow spectrofluorometric assay, respectively. As a measure of lipid scrambling, we used differential scanning calorimetry to monitor the effect of gramicidin on the melting transition temperatures of the two bilayer leaflets. The two calorimetric peaks of the individual leaflets merged into a single peak over time, suggestive of scrambling, and the effect of the channel on the transbilayer lipid distribution in both symmetric 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphocholine and asymmetric 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphocholine/1,2-dimyristoyl-sn-glycero-3-phosphocholine vesicles was quantified from proton NMR measurements. Our results show that gramicidin increases lipid flip-flop in a complex, concentration-dependent manner. To determine the molecular mechanism of the process, we used molecular dynamics simulations and further computational analysis of the trajectories to estimate the extent of membrane deformation. Together, the experimental and computational approaches were found to constitute an effective means for studying the effects of transmembrane proteins on lipid distribution in both symmetric and asymmetric model membranes.

Peptide-Induced Lipid Flip-Flop in Asymmetric Liposomes Measured by Small Angle Neutron Scattering.

Despite the prevalence of lipid transbilayer asymmetry in natural plasma membranes, most biomimetic model membranes studied are symmetric. Recent advances have helped to overcome the difficulties in preparing asymmetric liposomes in vitro, allowing for the examination of a larger set of relevant biophysical questions. Here, we investigate the stability of asymmetric bilayers by measuring lipid flip-flop with time-resolved small-angle neutron scattering (SANS). Asymmetric large unilamellar vesicles with inner bilayer leaflets containing predominantly 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphocholine (POPC) and outer leaflets composed mainly of 1,2-dimyristoyl-sn-glycero-3-phosphocholine (DMPC) displayed slow spontaneous flip-flop at 37 ◦C (half-time, t1/2 = 140 h). However, inclusion of peptides, namely, gramicidin, alamethicin, melittin, or pHLIP (i.e., pH-low insertion peptide), accelerated lipid flip-flop. For three of these peptides (i.e., pHLIP, alamethicin, and melittin), each of which was added externally to preformed asymmetric vesicles, we observed a completely scrambled bilayer in less than 2 h. Gramicidin, on the other hand, was preincorporated during the formation of the asymmetric liposomes and showed a time resolvable 8-fold increase in the rate of lipid asymmetry loss. These results point to a membrane surface-related (e.g., adsorption/insertion) event as the primary driver of lipid scrambling in the asymmetric model membranes of this study. We discuss the implications of membrane peptide binding, conformation, and insertion on lipid asymmetry.

Accurate In Silico Modeling of Asymmetric Bilayers Based on Biophysical Principles.

Technological advances in the last decade have enabled the study of ever more complex and physiologically relevant model membranes to help dispel the mystery surrounding the role of plasma membrane asymmetry in various cellular processes. The slowly accumulating body of experimental data is fueling renewed interest in and the need for computational methods to support interpretations and address a wide range of problems that are still not amenable to direct experimental study. The specific appeal of molecular dynamics simulations for this purpose lies in their ability to access information at atomic resolution, which is useful for the formulation of testable mechanistic hypotheses. But, the range of questions that can be addressed reliably with such simulations is determined by the appropriate construction and simulation of asymmetric bilayer models. One essential way to achieve this goal is to follow rigorous biophysical criteria and principles. In this context, we show that the requirement for a robust comparison between the properties of simulated asymmetric and symmetric model membranes is for the tension in each bilayer leaflet to be zero. Commonly used methods for constructing asymmetric bilayers, including matching the average areas of the leaflets from the corresponding symmetric systems, do not ensure zero leaflet tension, thus precluding physically realistic changes in the areas of the two leaflets. We present, to our knowledge, a new method for identifying the ideal lipid packing in bilayers with different leaflet compositions that achieves the zero-tension goal, and discuss the basic principles underlying the biophysically correct computational study of asymmetric membranes.

Cholesterol Promotes Protein Binding by Affecting Membrane Electrostatics and Solvation Properties.

Binding of the retroviral structural protein Gag to the cellular plasma membrane is mediated by the protein's matrix (MA) domain. Prominent among MA-PM interactions is electrostatic attraction between the positively charged MA domain and the negatively charged plasma membrane inner leaflet. Previously, we reported that membrane association of HIV-1 Gag, as well as purified Rous sarcoma virus (RSV) MA and Gag, depends strongly on the presence of acidic lipids and is enhanced by cholesterol (Chol). The mechanism underlying this enhancement was unclear. Here, using a broad set of in vitro and in silico techniques we addressed molecular mechanisms of association between RSV MA and model membranes, and investigated how Chol enhances this association. In neutron scattering experiments with liposomes in the presence or absence of Chol, MA preferentially interacted with preexisting POPS-rich clusters formed by nonideal lipid mixing, binding peripherally to the lipid headgroups with minimal perturbation to the bilayer structure. Molecular dynamics simulations showed a stronger MA-bilayer interaction in the presence of Chol, and a large Chol-driven increase in lipid packing and membrane surface charge density. Although in vitro MA-liposome association is influenced by disparate variables, including ionic strength and concentrations of Chol and charged lipids, continuum electrostatic theory revealed an underlying dependence on membrane surface potential. Together, these results conclusively show that Chol affects RSV MA-membrane association by making the electrostatic potential at the membrane surface more negative, while decreasing the penalty for lipid headgroup desolvation. The presented approach can be applied to other viral and nonviral proteins.

Subnanometer Structure of an Asymmetric Model Membrane: Interleaflet Coupling Influences Domain Properties.

Cell membranes possess a complex three-dimensional architecture, including nonrandom lipid lateral organization within the plane of a bilayer leaflet, and compositional asymmetry between the two leaflets. As a result, delineating the membrane structure-function relationship has been a highly challenging task. Even in simplified model systems, the interactions between bilayer leaflets are poorly understood, due in part to the difficulty of preparing asymmetric model membranes that are free from the effects of residual organic solvent or osmotic stress. To address these problems, we have modified a technique for preparing asymmetric large unilamellar vesicles (aLUVs) via cyclodextrin-mediated lipid exchange in order to produce tensionless, solvent-free aLUVs suitable for a range of biophysical studies. Leaflet composition and structure were characterized using isotopic labeling strategies, which allowed us to avoid the use of bulky labels. NMR and gas chromatography provided precise quantification of the extent of lipid exchange and bilayer asymmetry, while small-angle neutron scattering (SANS) was used to resolve bilayer structural features with subnanometer resolution. Isotopically asymmetric POPC vesicles were found to have the same bilayer thickness and area per lipid as symmetric POPC vesicles, demonstrating that the modified exchange protocol preserves native bilayer structure. Partial exchange of DPPC into the outer leaflet of POPC vesicles produced chemically asymmetric vesicles with a gel/fluid phase-separated outer leaflet and a uniform, POPC-rich inner leaflet. SANS was able to separately resolve the thicknesses and areas per lipid of coexisting domains, revealing reduced lipid packing density of the outer leaflet DPPC-rich phase compared to typical gel phases. Our finding that a disordered inner leaflet can partially fluidize ordered outer leaflet domains indicates some degree of interleaflet coupling, and invites speculation on a role for bilayer asymmetry in modulating membrane lateral organization.

Computational modeling of the N-terminus of the human dopamine transporter and its interaction with PIP2 -containing membranes.

The dopamine transporter (DAT) is a transmembrane protein belonging to the family of neurotransmitter:sodium symporters (NSS). Members of the NSS are responsible for the clearance of neurotransmitters from the synaptic cleft, and for their translocation back into the presynaptic nerve terminal. The DAT contains long intracellular N- and C-terminal domains that are strongly implicated in the transporter function. The N-terminus (N-term), in particular, regulates the reverse transport (efflux) of the substrate through DAT. Currently, the molecular mechanisms of the efflux remain elusive in large part due to lack of structural information on the N-terminal segment. Here we report a computational model of the N-term of the human DAT (hDAT), obtained through an ab initio structure prediction, in combination with extensive atomistic molecular dynamics (MD) simulations in the context of a lipid membrane. Our analysis reveals that whereas the N-term is a highly dynamic domain, it contains secondary structure elements that remain stable in the long MD trajectories of interactions with the bilayer (totaling >2.2 µs). Combining MD simulations with continuum mean-field modeling we found that the N-term engages with lipid membranes through electrostatic interactions with the charged lipids PIP2 (phosphatidylinositol 4,5-Biphosphate) or PS (phosphatidylserine) that are present in these bilayers. We identify specific motifs along the N-term implicated in such interactions and show that differential modes of N-term/membrane association result in differential positioning of the structured segments on the membrane surface. These results will inform future structure-based studies that will elucidate the mechanistic role of the N-term in DAT function.