High-resolution structure of coexisting nanoscopic and microscopic lipid domains

Multiscale Modeling of Macromolecular Interactions between Tau-Amylin Oligomers and Asymmetric Lipid Nanodomains That Link Alzheimer's and Diabetic Diseases

The molecular events of protein misfolding and self-aggregation of tau and amylin are associated with the progression of Alzheimer's and diabetes, respectively. Recent studies suggest that tau and amylin can form hetero-tau-amylin oligomers. Those hetero-oligomers are more neurotoxic than homo-tau oligomers. So far, the detailed interactions between the hetero-oligomers and the neuronal membrane are unknown. Using multiscale MD simulations, the lipid binding and protein folding behaviors of hetero-oligomers on asymmetric lipid nanodomains or raft membranes were examined. Our raft membranes contain phase-separated phosphatidylcholine (PC), cholesterol, and anionic phosphatidylserine (PS) or ganglioside (GM1) in one leaflet of the lipid bilayer. The hetero-oligomers bound more strongly to the PS and GM1 than other lipids via the hydrophobic and hydrophilic interactions, respectively, in the raft membranes. The hetero-tetramer disrupted the acyl chain orders of both PC and PS in the PS-containing raft membrane, but only the GM1 in the GM1-containing raft membrane as effectively as the homo-tau-tetramer. We discovered that the alpha-helical content in the heterodimer was greater than the sum of alpha-helical contents from isolated tau and amylin monomers on both raft membranes, indicative of a synergetic effect of tau-amylin interactions in surface-induced protein folding. Our results provide new molecular insights into understanding the cross-talk between Alzheimer's and diabetes.

Quantitative Comparison against Experiments Reveals Imperfections in Force Fields' Descriptions of POPC-Cholesterol Interactions

Cholesterol is a central building block in biomembranes, where it induces orientational order, slows diffusion, renders the membrane stiffer, and drives domain formation. Molecular dynamics (MD) simulations have played a crucial role in resolving these effects at the molecular level; yet, it has recently become evident that different MD force fields predict quantitatively different behavior. Although easily neglected, identifying such limitations is increasingly important as the field rapidly progresses toward simulations of complex membranes mimicking the in vivo conditions: pertinent multicomponent simulations must capture accurately the interactions between their fundamental building blocks, such as phospholipids and cholesterol. Here, we define quantitative quality measures for simulations of binary lipid mixtures in membranes against the C-H bond order parameters and lateral diffusion coefficients from NMR spectroscopy as well as the form factors from X-ray scattering. Based on these measures, we perform a systematic evaluation of the ability of commonly used force fields to describe the structure and dynamics of binary mixtures of palmitoyloleoylphosphatidylcholine (POPC) and cholesterol. None of the tested force fields clearly outperforms the others across the tested properties and conditions. Still, the Slipids parameters provide the best overall performance in our tests, especially when dynamic properties are included in the evaluation. The quality evaluation metrics introduced in this work will, particularly, foster future force field development and refinement for multicomponent membranes using automated approaches.

Lipid domain boundary triggers membrane damage and protein folding of human islet amyloid polypeptide in the early pathogenesis of amyloid diseases

The misfolding and self-aggregation of human Islet Amyloid Polypeptide (hIAPP) are linked to the onset of type 2 diabetes (T2D). However, the mechanism of how the disordered hIAPP aggregates trigger membrane damage leading to the loss of Islet cells in T2D is unknown. Using coarse-grained (CG) and all-atom (AA) molecular dynamics simulations, we have investigated the membrane-disruption behaviors of hIAPP oligomers on the phase-separated lipid nanodomains that mimic the highly heterogeneous lipid raft structures of cell membranes. Our results revealed that hIAPP oligomers preferentially bind to the liquid-ordered and liquid-disordered domain boundary around two hydrophobic residues at L16 and I26, and lipid acyl chain order disruption and beta-sheet formation occur upon hIAPP binding to the membrane surface. We propose that the lipid order disruption and surface-induced beta-sheet formation on the lipid domain boundary represent the early molecular events of membrane damage associated with the early pathogenesis of T2D.

Exploring the binding kinetics and behaviors of self-aggregated beta-amyloid oligomers to phase-separated lipid rafts with or without ganglioside-clusters

Lipid binding kinetics and energetics of self-aggregated and disordered beta-amyloid oligomers of various sizes, from solution to lipid raft surfaces, were investigated using MD simulations. Our systems include small (monomers to tetramers) and larger (octamers and dodecamers) oligomers. Our lipid rafts contain saturated and unsaturated phosphatidylcholine (PC), cholesterol, and with or without asymmetrically distributed monosialotetrahexosylganglioside (GM1). All rafts exhibited dynamic and structurally diversified domains including liquid-ordered (Lo), liquid-disordered (Ld), and interfacial Lod domains. For rafts without GM1, all oligomers bound to the Lod domain. For GM1-containing rafts, all small oligomers and most larger oligomers bound specifically to the GM1-clusters embedded in the Lo domain. Lipid-protein binding energies followed an order of GM1 >> unsaturated PC > saturated PC > cholesterol for all rafts. In addition, protein-induced membrane structural disruption increased progressively with the size of the oligomer for the annular lipids surrounding the membrane-bound protein in non-GM1-containing rafts. We propose that the tight binding of beta-amyloid oligomers to the GM1-clusters and the structural perturbation of lipids surrounding the membrane-bound proteins at the Lod domain are early molecular events of the beta-amyloid aggregation process on neuronal membrane surfaces that trigger the onset of Alzheimer's.

Macrocyclic Peptide-Conjugated Tip for Fast and Selective Molecular Recognition Imaging by High-Speed Atomic Force Microscopy

Fast and selective recognition of molecules at the nanometer scale without labeling is a much desired but still challenging goal to achieve. Here, we show the use of high-speed atomic force microscopy (HS-AFM) for real-time and real-space recognition of unlabeled membrane receptors using tips conjugated with small synthetic macrocyclic peptides. The single-molecule recognition method is validated by experiments on the human hepatocyte growth factor receptor (hMET), which selectively binds to the macrocyclic peptide aMD4. By testing and comparing aMD4 synthesized with linkers of different lengths and rigidities, we maximize the interaction between the functionalized tip and hMET added to both a mica surface and supported lipid bilayers. Phase contrast imaging by HS-AFM enables us to discriminate nonlabeled hMET against the murine MET homologue, which does not bind to aMD4. Moreover, using ligands and linkers of small size, we achieve minimal deterioration of the spatial resolution in simultaneous topographic imaging. The versatility of macrocyclic peptides in detecting unlimited types of membrane receptors with high selectivity and the fast imaging by HS-AFM broaden the range of future applications of this method for molecular recognition without labeling.

Increasing complexity in small-angle X-ray and neutron scattering experiments: from biological membrane mimics to live cells

Small-angle X-ray and neutron scattering are well-established, non-invasive experimental techniques to interrogate global structural properties of biological membrane mimicking systems under physiologically relevant conditions. Recent developments, both in bottom-up sample preparation techniques for increasingly complex model systems, and in data analysis techniques have opened the path toward addressing long standing issues of biological membrane remodelling processes. These efforts also include emerging quantitative scattering studies on live cells, thus enabling a bridging of molecular to cellular length scales. Here, we review recent progress in devising compositional models for joint small-angle X-ray and neutron scattering studies on diverse membrane mimics - with a specific focus on membrane structural coupling to amphiphatic peptides and integral proteins - and live Escherichia coli. In particular, we outline the present state-of-the-art in small-angle scattering methods applied to complex membrane systems, highlighting how increasing system complexity must be followed by an advance in compositional modelling and data-analysis tools.

Structure and Thermotropic Behavior of Bovine- and Porcine-Derived Exogenous Lung Surfactants

Two commercial exogenous pulmonary surfactants, Curosurf and Survanta, are investigated. Their thermotropic behavior and associated structural changes for the samples in bulk are characterized and described. For Survanta, the obtained results of differential scanning calorimetry showed a thermogram with three peaks on heating and only a single peak on cooling. Curosurf on the other hand, presents calorimetric thermograms with only one peak in both the heating and cooling scans. This distinct thermotropic behavior between the two pulmonary surfactants, a consequence of their particular compositions, is associated with structural changes that were evaluated by simultaneous small- and wide-angle X-ray scattering experiments with in situ temperature variation. Interestingly, for temperatures below ∼35 °C for Curosurf and ∼53 °C for Survanta, the scattering data indicated the coexistence of two lamellar phases with different carbon chain organizations. For temperatures above these limits, the coexistence of phases disappears, giving rise to a fluid phase in both pulmonary surfactants, with multilamelar vesicles for Curosurf and unilamellar vesicles for Survanta. This process is quasi-reversible under cooling, and advanced data analysis for the scattering data indicated differences in the structural and elastic properties of the pulmonary surfactants. The detailed and systematic investigation shown in this work expands on the knowledge of the structure and thermodynamic behavior of Curosurf and Survanta, being relevant from both physiological and biophysical perspectives and also providing a basis for further studies on other types of pulmonary surfactants.

Two Coexisting Membrane Structures Are Defined by Lateral and Transbilayer Interactions between Sphingomyelin and Cholesterol

The structure of fully hydrated bilayers composed of equimolar proportions of palmitoylsphingomyelin (PSM) and cholesterol has been examined by synchrotron X-ray powder diffraction and atomistic molecular dynamics (MD) simulations. Two coexisting bilayer structures, which are distinguished by the transbilayer phosphate-phosphate distance of coupled PSM molecules, are observed by diffraction at 37 °C. The MD simulations reveal that PSM molecules in the thicker membrane are characterized by more ordered, more extended, and less interdigitated hydrocarbon tails compared to those in the thinner membrane. Intermolecular hydrogen bonds further distinguish the two bilayer structures, and we observe the disruption of a sphingomyelin intermolecular hydrogen bond network induced by the proximity of cholesterol. Through an unsupervised clustering of interatomic distances, we show for the first time that the asymmetry of phospholipids is important in driving their interactions with cholesterol. We identify four distinct modes of interaction, two of which lead to the dehydration of cholesterol. These two modes of interaction provide the first description of precise physical mechanisms underlying the umbrella model, which itself explains how phospholipids may shield cholesterol from water. The most dehydrating mode of interaction is particular to the N-acylated fatty acid moiety of PSM and thus may explain the long-held observation that cholesterol preferentially mixes with sphingomyelins over glycerophospholipids.

Membrane Thinning Induces Sorting of Lipids and the Amphipathic Lipid Packing Sensor (ALPS) Protein Motif

Heterogeneities (e.g., membrane proteins and lipid domains) and deformations (e.g., highly curved membrane regions) in biological lipid membranes cause lipid packing defects that may trigger functional sorting of lipids and membrane-associated proteins. To study these phenomena in a controlled and efficient way within molecular simulations, we developed an external field protocol that artificially enhances packing defects in lipid membranes by enforcing local thinning of a flat membrane region. For varying lipid compositions, we observed strong thinning-induced depletion or enrichment, depending on the lipid's intrinsic shape and its effect on a membrane's elastic modulus. In particular, polyunsaturated and lysolipids are strongly attracted to regions high in packing defects, whereas phosphatidylethanolamine (PE) lipids and cholesterol are strongly repelled from it. Our results indicate that externally imposed changes in membrane thickness, area, and curvature are underpinned by shared membrane elastic principles. The observed sorting toward the thinner membrane region is in line with the sorting expected for a positively curved membrane region. Furthermore, we have demonstrated that the amphipathic lipid packing sensor (ALPS) protein motif, a known curvature and packing defect sensor, is effectively attracted to thinner membrane regions. By extracting the force that drives amphipathic molecules toward the thinner region, our thinning protocol can directly quantify and score the lipid packing sensing of different amphipathic molecules. In this way, our protocol paves the way toward high-throughput exploration of potential defect- and curvature-sensing motifs, making it a valuable addition to the molecular simulation toolbox.

Magainin 2 and PGLa in Bacterial Membrane Mimics I: Peptide-Peptide and Lipid-Peptide Interactions

We addressed the onset of synergistic activity of the two well-studied antimicrobial peptides magainin 2 (MG2a) and PGLa using lipid-only mimics of Gram-negative cytoplasmic membranes. Specifically, we coupled a joint analysis of small-angle x-ray and neutron scattering experiments on fully hydrated lipid vesicles in the presence of MG2a and L18W-PGLa to all-atom and coarse-grained molecular dynamics simulations. In agreement with previous studies, both peptides, as well as their equimolar mixture, were found to remain upon adsorption in a surface-aligned topology and to induce significant membrane perturbation, as evidenced by membrane thinning and hydrocarbon order parameter changes in the vicinity of the inserted peptide. These effects were particularly pronounced for the so-called synergistic mixture of 1:1 (mol/mol) L18W-PGLa/MG2a and cannot be accounted for by a linear combination of the membrane perturbations of two peptides individually. Our data are consistent with the formation of parallel heterodimers at concentrations below a synergistic increase of dye leakage from vesicles. Our simulations further show that the heterodimers interact via salt bridges and hydrophobic forces, which apparently makes them more stable than putatively formed antiparallel L18W-PGLa and MG2a homodimers. Moreover, dimerization of L18W-PGLa and MG2a leads to a relocation of the peptides within the lipid headgroup region as compared to the individual peptides. The early onset of dimerization of L18W-PGLa and MG2a at low peptide concentrations consequently appears to be key to their synergistic dye-releasing activity from lipid vesicles at high concentrations.

Fluid Phase Coexistence in Biological Membrane: Insights from Local Nonaffine Deformation of Lipids

Lateral heterogeneities in biomembranes play a crucial role in various physiological functions of the cell. Such heterogeneities lead to demixing of lipid constituents and formation of distinct liquid domains in the membrane. We study lateral heterogeneities in terms of topological rearrangements of lipids to identify the liquid-liquid phase coexistence in model membranes. Using ideas from the physics of amorphous systems and glasses, we calculate the degree of nonaffine deformation associated with individual lipids to characterize the liquid-ordered (L_o) and liquid-disordered (L_d) regions in model lipid bilayers. We explore the usage of this method on all-atom and coarse-grained lipid bilayer trajectories. This method is helpful in defining the instantaneous L_o-L_d domain boundaries in complex multicomponent bilayer systems. The characterization is also used to highlight the effect of line-active molecules on the phase boundaries and domain mixing. Overall, we propose a framework to explore the molecular origin of spatial and dynamical heterogeneity in biomembrane systems, which can be exploited not only in computer simulations but also in experiments.

Exploring the impact of proteins on the line tension of a phase-separating ternary lipid mixture

The separation of lipid mixtures into thermodynamically stable phase-separated domains is dependent on lipid composition, temperature, and system size. Using molecular dynamics simulations, the line tension between thermodynamically stable lipid domains formed from ternary mixtures of di-C16:0 PC:di-C18:2 PC:cholesterol at 40:40:20 mol. % ratio was investigated via two theoretical approaches. The line tension was found to be 3.1 ± 0.2 pN by capillary wave theory and 4.7 ± 3.7 pN by pressure tensor anisotropy approaches for coarse-grained models based on the Martini force field. Using an all-atom model of the lipid membrane based on the CHARMM36 force field, the line tension was found to be 3.6 ± 0.9 pN using capillary wave theory and 1.8 ± 2.2 pN using pressure anisotropy approaches. The discrepancy between estimates of the line tension based on capillary wave theory and pressure tensor anisotropy methods is discussed. Inclusion of protein in Martini membrane lipid mixtures was found to reduce the line tension by 25%-35% as calculated by the capillary wave theory approach. To further understand and predict the behavior of proteins in phase-separated membranes, we have formulated an analytical Flory-Huggins model and parameterized it against the simulation results. Taken together these results suggest a general role for proteins in reducing the thermodynamic cost associated with domain formation in lipid mixtures and quantifies the thermodynamic driving force promoting the association of proteins to domain interfaces.

Combined effects of osmotic and hydrostatic pressure on multilamellar lipid membranes in the presence of PEG and trehalose

We studied the interaction of lipid membranes with the disaccharide trehalose (TRH), which is known to stabilize biomembranes against various environmental stress factors. Generally, stress factors include low/high temperature, shear, osmotic and hydrostatic pressure. Small-angle X-ray-scattering was applied in combination with fluorescence spectroscopy and calorimetric measurements to get insights into the influence of trehalose on the supramolecular structure, hydration level, and elastic and thermodynamic properties as well as phase behavior of the model biomembrane DMPC, covering a large region of the temperature, osmotic and hydrostatic pressure phase space. We observed distinct effects of trehalose on the topology of the lipid's supramolecular structure. Trehalose, unlike osmotic pressure induced by polyethylene glycol, leads to a decrease of lamellar order and a swelling of multilamellar vesicles, which is attributable to direct interactions between the membrane and trehalose. Our results revealed a distinct biphasic concentration dependence of the observed effects of trehalose. While trehalose intercalates between the polar head groups at low concentrations, the effects after saturation are dominated by the exclusion of trehalose from the membrane surface.

Membrane Lipid Nanodomains

Lipid membranes can spontaneously organize their components into domains of different sizes and properties. The organization of membrane lipids into nanodomains might potentially play a role in vital functions of cells and organisms. Model membranes represent attractive systems to study lipid nanodomains, which cannot be directly addressed in living cells with the currently available methods. This review summarizes the knowledge on lipid nanodomains in model membranes and exposes how their specific character contrasts with large-scale phase separation. The overview on lipid nanodomains in membranes composed of diverse lipids (e.g., zwitterionic and anionic glycerophospholipids, ceramides, glycosphingolipids) and cholesterol aims to evidence the impact of chemical, electrostatic, and geometric properties of lipids on nanodomain formation. Furthermore, the effects of curvature, asymmetry, and ions on membrane nanodomains are shown to be highly relevant aspects that may also modulate lipid nanodomains in cellular membranes. Potential mechanisms responsible for the formation and dynamics of nanodomains are discussed with support from available theories and computational studies. A brief description of current fluorescence techniques and analytical tools that enabled progress in lipid nanodomain studies is also included. Further directions are proposed to successfully extend this research to cells.

Regimes of Complex Lipid Bilayer Phases Induced by Cholesterol Concentration in MD Simulation

Cholesterol is essential to the formation of phase-separated lipid domains in membranes. Lipid domains can exist in different thermodynamic phases depending on the molecular composition and play significant roles in determining structure and function of membrane proteins. We investigate the role of cholesterol in the structure and dynamics of ternary lipid mixtures displaying phase separation using molecular dynamics simulations, employing a physiologically relevant span of cholesterol concentration. We find that cholesterol can induce formation of three regimes of phase behavior: 1) miscible liquid-disordered bulk, 2) phase-separated, domain-registered coexistence of liquid-disordered and liquid-ordered domains, and 3) phase-separated, domain-antiregistered coexistence of liquid-disordered and newly identified nanoscopic gel domains composed of cholesterol threads we name "cholesterolic gel" domains. These findings are validated and discussed in the context of current experimental knowledge, models of cholesterol spatial distributions, and models of ternary lipid-mixture phase separation.

Cholesterol Flip-Flop Impacts Domain Registration in Plasma Membrane Models

The plasma membrane is a highly complex multicomponent system that is central to the functioning of cells. Cholesterol, a key lipid component of the plasma membrane, promotes the formation of nanodomains. These nanodomains are often correlated across the two membrane leaflets, but the underlying physical mechanism remains unclear. Using coarse-grained molecular dynamics simulations, we investigate the influence of cholesterol flip-flop on membrane properties, in particular, the interleaflet coupling of cholesterol-enriched domains. We show that the cholesterol density correlation between the leaflets of an average mammalian plasma membrane is significantly reduced by suppressing interleaflet cholesterol population. Our results suggest an amplifying role of cholesterol in signal transduction across the leaflets.

Complex dynamics at the nanoscale in simple biomembranes

Nature is known to engineer complex compositional and dynamical platforms in biological membranes. Understanding this complex landscape requires techniques to simultaneously detect membrane re-organization and dynamics at the nanoscale. Using super-resolution stimulated emission depletion (STED) microscopy coupled with fluorescence correlation spectroscopy (FCS), we reveal direct experimental evidence of dynamic heterogeneity at the nanoscale in binary phospholipid-cholesterol bilayers. Domain formation on the length scale of ~200–600 nm due to local cholesterol compositional heterogeneity is found to be more prominent at high cholesterol content giving rise to distinct intra-domain lipid dynamics. STED-FCS reveals unique dynamical crossover phenomena at length scales of ~100–150 nm within each of these macroscopic regions. The extent of dynamic heterogeneity due to intra-domain hindered lipid diffusion as reflected from the crossover length scale, is driven by cholesterol packing and organization, uniquely influenced by phospholipid type. These results on simple binary model bilayer systems provide novel insights into pathways leading to the emergence of complex nanodomain substructures with implications for a wide variety of membrane mediated cellular events.

Complex biomembrane mimetics on the sub-nanometer scale

Biomimetic lipid vesicles are indispensable tools for gaining insight into the biophysics of cell physiology on the molecular level. The level of complexity of these model systems has steadily increased, and now spans from domain-forming lipid mixtures to asymmetric lipid bilayers. Here, we review recent progress in the development and application of elastic neutron and X-ray scattering techniques for studying these systems in situ and under physiologically relevant conditions on the nanometer to sub-nanometer length scales. In particular, we focus on: (1) structural details of coexisting liquid-ordered and liquid-disordered domains, including their thickness and lipid packing mismatch as a function of a size transition from nanoscopic to macroscopic domains; (2) membrane-mediated protein partitioning into lipid domains; (3) the role of the aqueous medium in tuning interactions between membranes and domains; and (4) leaflet-specific structure in asymmetric bilayers and passive lipid flip-flop.