The molecular face of lipid rafts in model membranes

Lateral Organization and Dynamics of the Realistic Plasma Membrane

Large-scale simulations of realistic crowded cell membranes can bridge the gap between the simulations and experiments. However, the compositional complexity and structural asymmetry of cell membranes continue to pose significant challenges in computational biology. Recent advances in understanding native membranes, including their composition and protein structures, enable us to construct a highly realistic model of the mammalian plasma membrane. Using this model, we explore the organization and dynamics of biological cell membranes at the molecular level. We found that the interaction preferences of protein-lipid mediate the formation of dynamic clusters of nonrandomly distributed proteins, accompanied by heterogeneous structural properties and anomalous diffusion. These evolving dynamic clusters intertwine to form a highly complex and continuously changing protein network. Our study provides significant insights into the intricate lateral dynamic organization of cell membranes.

Interaction of lipid domains originating from differential domain-monolayer contact energy

We consider a flat membrane containing pure lipid domains located in the membrane monolayers and separated in the membrane plane. We assume the energy of contact along the membrane mid-surface between a domain and the underlying monolayer to be different from that between the two monolayers. We theoretically analyse the effect of the differential contact energy on the elastic deformations of tilt and splay in the membrane monolayers and the resulting interaction between two domains situated in the apposed monolayers. We demonstrate that the character of this interaction depends on the ratio, η, between the domain rigidity and that of a regular membrane monolayer. For the rigidity ratio smaller than a critical value, η < η* ≈ 3, the domain interaction is predicted to be attractive for all inter-monolayer distances. For the super-critical values of the rigidity ratio, η > η*, the interaction is repulsive for small distances and attractive for large distances with a certain equilibrium inter-domain separation corresponding to a vanishing interaction force. The predicted attractive interaction is proposed to favor the registration in the membrane plane of apposed domains as observed in most domain-containing membranes.

Signatures of glassy dynamics in highly ordered lipid bilayers with emergence of soft dynamic channels

Over the last few decades, extensive investigations on spatial and dynamic heterogeneity have been performed on carefully reconstituted biological lipid membranes. Characterizing the molecular features in heterogeneous membranes is extremely challenging due to the experimentally inaccessible time- and length-scales of these emergent systems. In this context, simulations can provide important insights into molecular-level interactions leading to membrane heterogeneity and associated functions. To that end, we use the non-affine displacement (NAD) framework (a concept borrowed from the physics of granular materials) to faithfully capture molecular-scale local membrane order in simulated heterogeneous bilayers. In our latest application of NAD, we investigate the temperature-dependent spatial and temporal organization on microsecond trajectories of liquid-ordered bilayer systems at all-atom resolution (DPPC/DOPC/CHOL: 0.55:0.15:0.30; 40 × 40 nm2 with a total of 5600 lipids and 2 × 106 atoms). Lateral organization in these large bilayer patches shows noticeable dynamic heterogeneity despite their liquid-ordered nature. Moreover, our NAD analyses reveal soft fluid channels within the tightly packed membrane reminiscent of the classical two-component Kob-Andersen glass-forming binary mixture. Hence, we characterized these systems using classical glass physics markers for dynamic heterogeneities such as overlap, four-point susceptibility, Van Hove, and intermediate scattering functions to quantify the multiple time scales underlying the lipid dynamics. Our analyses reveal that highly ordered membrane systems can have glass-like dynamics with distinct soft fluid channels inside them. Biologically, these dynamic channels could act as conduits for facilitating molecular encounters for biological functions even in highly ordered phases such as lipid nanodomains and rafts.

Mechanistic Insight on Ethanol Driven Swelling and Disruption of Cholesterol Containing Biomimetic Vesicles From Coarse-Grained Molecular Dynamics

We have performed coarse-grained (CG) molecular dynamics (MD) simulations to delineate the impact of ethanol (EtOH) on cholesterol (CHOL) containing biomimetic bilayer and vesicle composed of 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphocholine (POPC) lipids. We have first deduced the missing interaction parameters for the POPC-CHOL-EtOH-water system within the SPICA/SDK CG force-field (CG-FF). By monitoring the electron density profiles, the orientational order parameter, and reproducing the all-atom MD-derived free energy for the insertion of ethanol from the bulk aqueous phase to the hydrophobic core of the POPC-CHOL lipid bilayer, we successfully determined all the missing non-bonding interaction parameters for the POPC-CHOL-EtOH-water system. The proposed force field was applied to investigate the effect of ethanol at various concentrations on unilamellar vesicles composed of POPC and cholesterol. It was found that 40 mol% or more concentration of ethanol is required to disintegrate or rupture the POPC-CHOL vesicle membranes. While cholesterol offers some resilience against the detrimental effects of ethanol, we still observe an increase in vesicle size (swelling) and a contraction in the bilayer thickness (thinning) as ethanol concentration rises from 0 to 30 mol%. At ethanol concentrations exceeding 30 mol%, the vesicles become increasingly susceptible to disintegration due to enhanced penetration of ethanol and water molecules into the hydrophobic core of the membranes.

The Importance of Bilayer Asymmetry in Biological Membranes: Insights from Model Membranes

This mini-review intends to highlight the importance of bilayer asymmetry. Biological membranes are complex structures that are a physical barrier separating the external environment from the cellular content. This complex bilayer comprises an extensive lipid repertory, suggesting that the different lipid structures might play a role in the membrane. Interestingly, this vast repertory of lipids is asymmetrically distributed between leaflets that form the lipid bilayer. Here, we discuss the properties of the plasma membrane from the perspective of experimental model membranes, consisting of simplified and controlled in vitro systems. We summarize some crucial features of the exoplasmic (outer) and cytoplasmic (inner) leaflets observed through investigations using symmetric and asymmetric membranes. Symmetric model membranes for the exoplasmic leaflet have a unique lipid composition that might form a coexistence of phases, namely the liquid disordered and liquid order phases. These phase domains may appear in different sizes and shapes depending on lipid composition and lipid-lipid interactions. In contrast, symmetric model membranes for the cytoplasmic leaflet form a fluid phase. We discuss the outcomes reported in the literature for asymmetric bilayers, which vary according to lipid compositions and, consequently, reflect different intra- and inter-leaflet interactions. Interestingly, the asymmetric bilayer could show induced domains in the inner leaflet, or it could decrease the tendency of the outer leaflet to phase separation. If cells regulate the lipid composition of the plasma membrane, they can adjust the existence and sizes of the domains by tuning the lipid composition.

Understanding the coarse-grained free energy landscape of phospholipids and their phase separation

The cell membrane exhibits lateral heterogeneity due to the preferential association among the large number of lipid species that constitute the membrane. In particular, the preferential association of cholesterol (CHOL) with saturated lipids into ordered domains has been an area of intense investigation. The large spatiotemporal scales that comprise spontaneous domain formation largely precludes computational investigation via conventional all-atom molecular dynamics. We demonstrate here that molecular coarse-grained (CG) models, obtained from the bottom-up, i.e., via statistical mechanical renormalization of atomistic models, are capable of spontaneous assembly and phase separation for two model raft-like systems, DLiPC/DPPC/CHOL and DOPC/DPPC/CHOL. The resulting bottom-up CG models exhibit spontaneous self-assembly and phase separation and recapitulate the structural correlations of the underlying atomistic models. The accuracy and fast dynamics of these CG models constitute an effective means of bypassing the limited spatiotemporal scales of atomistic simulations. As the first bottom-up CG models of lipid phase separation, the CG models in this work provide an informative analysis for further construction of bottom-up CG models transferable across a range of lipid compositions.

Regulate PD-L1's membrane orientation thermodynamics with hydrophobic nanoparticles

Tumor cells can escape from immune killing by binding their programmed death ligand-1 (PD-L1) to the programmed cell death protein 1 (PD-1) of T cells. These immune checkpoint proteins (PD-L1/PD-1) have become very important drug targets, since blocking PD-L1 or PD-1 can recover the killing capability of T cells against tumor cells. Instead of targeting the binding interface between PD-L1 and PD-1, we explored the possibility of regulating the membrane orientation thermodynamics of PD-L1 with ligand-modified ultra-small hydrophobic nanoparticles (NPs) using μs-scale coarse-grained molecular dynamics (MD) simulations in this work. Our MD results indicate that embedded hydrophobic NPs can significantly change the membrane orientation thermodynamics of the extracellular domain of PD-L1, enhancing the probability in the "stand up" state for better binding to PD-1. Meanwhile, embedded hydrophobic NPs promote the tilt of the transmembrane domain of PD-L1. Besides, effects on both extracellular and transmembrane domains are determined by the ligand length and NP concentration. Our study may provide an alternative strategy to achieve PD-L1-related immunotherapy with nanomedicine.

Enhancing the Assembly Properties of Bottom-Up Coarse-Grained Phospholipids

A plethora of key biological events occur at the cellular membrane where the large spatiotemporal scales necessitate dimensionality reduction or coarse-graining approaches over conventional all-atom molecular dynamics simulation. Constructing coarse-grained descriptions of membranes systematically from statistical mechanical principles has largely remained challenging due to the necessity of capturing amphipathic self-assembling behavior in coarse-grained models. We show that bottom-up coarse-grained lipid models can possess metastable morphological behavior and that this potential metastability has ramifications for accurate development and training. We in turn develop a training algorithm which evades metastability issues by linking model training to self-assembling behavior, and demonstrate its robustness via construction of solvent-free coarse-grained models of various phospholipid membranes, including lipid species such as phosphatidylcholines, phosphatidylserines, sphingolipids, and cholesterol. The resulting coarse-grained lipid models are orders of magnitude faster than their atomistic counterparts while retaining structural fidelity and constitute a promising direction for the development of coarse-grained models of realistic cell membranes.

All-Atom Membrane Builder via Multiscale Simulation

I present an automated and flexible tool designed for constructing bilayer membranes at all-atom (AA) resolution. The builder initiates the construction and equilibration of bilayer membranes at Martini coarse-grained (CG) resolution, followed by resolution enhancement to the atomic level using the accompanying backmapping tool. Notably, this tool enables users to create bilayer membranes with user-defined lipid compositions and protein structures, while also offering the flexibility to accommodate new lipid types. To assess the simplicity and robustness of the tool, I demonstrate the construction of several membranes incorporating protein structures. The tool is freely available at github.com/ksy141/mstool.

Transient pores in hemifusion diaphragms

Exchange of material across two membranes, as in the case of synaptic neurotransmitter release from a vesicle, involves the formation and poration of a hemifusion diaphragm (HD). The nontrivial geometry of the HD leads to environment-dependent control, regarding the stability and dynamics of the pores required for this kind of exocytosis. This work combines particle simulations, field-based calculations, and phenomenological modeling to explore the factors influencing the stability, dynamics, and possible control mechanisms of pores in HDs. We find that pores preferentially form at the HD rim, and that their stability is sensitive to a number of factors, including the three line tensions, membrane tension, HD size, and the ability of lipids to "flip-flop" across leaflets. Along with a detailed analysis of these factors, we discuss ways that vesicles or cells may use them to open and close pores and thereby quickly and efficiently transport material.

Investigating the Influence of Photoswitchable Lipids on the Structure and Dynamics of Lipid Membranes: Fundamentals and Potential Applications

In this work, we delve into the impact of photoisomerization of photoswitchable lipids (PSLs) on the membrane structure and dynamics at a molecular level. Through all-atom molecular dynamics simulations, we explore how UV irradiation-induced trans-to-cis isomerization of these lipids, particularly the azobenzene-derivatized phosphatidylcholine (AzoPC) lipid, influences the structure and dynamics of a simplified lipid membrane, mimicking those of E. coli bacteria across different temperatures. Our findings align with previous experimental observations regarding membrane properties and offer insights into localized effects and microscopic heterogeneity. Additionally, we estimate the relaxation time scale of the lipid membrane following AzoPC photoisomerization. Moreover, we demonstrate the feasibility of photoactivated drug release, exemplified by the controlled liberation of doxorubicin, an anticancer agent, through the membrane, suggesting the potential of PSLs in engineering photoactivated liposomes, coined as photoazosomes, for precise targeted drug delivery applications.

There and back again: bridging meso- and nano-scales to understand lipid vesicle patterning

We describe a complete methodology to bridge the scales between nanoscale molecular dynamics and (micrometer) mesoscale Monte Carlo simulations in lipid membranes and vesicles undergoing phase separation, in which curving molecular species are furthermore embedded. To go from the molecular to the mesoscale, we notably appeal to physical renormalization arguments enabling us to rigorously infer the mesoscale interaction parameters from its molecular counterpart. We also explain how to deal with the physical timescales at stake at the mesoscale. Simulating the as-obtained mesoscale system enables us to equilibrate the long wavelengths of the vesicles of interest, up to the vesicle size. Conversely, we then backmap from the meso- to the nano-scale, which enables us to equilibrate in turn the short wavelengths down to the molecular length-scales. By applying our approach to the specific situation of patterning a vesicle membrane, we show that macroscopic membranes can thus be equilibrated at all length-scales in achievable computational time offering an original strategy to address the fundamental challenge of timescale in simulations of large bio-membrane systems.

Measuring the Transmembrane Registration of Lipid Domains in Droplet Interface Bilayers through Tensiometry

Diverse collections of lipids self-assemble into domains within biological membranes, and these domains are typically organized in both the transverse and lateral directions of the membrane. The ability of the membrane to link these domains across the membrane's interior grants cells control over features on the external cellular surface. Numerous hypothesized factors drive the cross-membrane (or transverse) coupling of lipid domains. In this work we seek to isolate these transverse lipid-lipid influences in a simple model system using droplet interface bilayers (DIBs) to better understand the associated mechanics. DIBs enable symmetric and asymmetric combinations of domain-forming lipid mixtures within a model bilayer, and the evolving energetics of the membrane may be tracked using drop-shape analysis. We find that symmetric distributions of domain-forming lipids produce long-lasting, gradual shifts in the DIB membrane energetics that are not observed in asymmetric distributions of the lipids where the domain-forming lipids are only within one leaflet. The approach selected for this work provides experimental measurement of the mismatch penalty associated with antiregistered lipid domains as well as measurements of the influence of rafts on DIB behaviors with suggestions for their future use as a model platform.

Building complex membranes with Martini 3

The Martini model is a popular force field for coarse-grained simulations. Membranes have always been at the center of its development, with the latest version, Martini 3, showing great promise in capturing more and more realistic behavior. In this chapter we provide a step-by-step tutorial on how to construct starting configurations, run initial simulations and perform dedicated analysis for membrane-based systems of increasing complexity, including leaflet asymmetry, curvature gradients and embedding of membrane proteins.

Characterization of domain formation in complex membranes

A widely known property of lipid membranes is their tendency to undergo a separation into disordered (Ld) and ordered (Lo) domains. This impacts the local structure of the membrane relevant for the physical (e.g., enhanced electroporation) and biological (e.g., protein sorting) significance of these regions. The increase in computing power, advancements in simulation software, and more detailed information about the composition of biological membranes shifts the study of these domains into the focus of classical molecular dynamics simulations. In this chapter, we present a versatile yet robust analysis pipeline that can be easily implemented and adapted for a wide range of lipid compositions. It employs Gaussian-based Hidden Markov Models to predict the hidden order states of individual lipids by describing their structure through the area per lipid and the average SCC order parameters per acyl chain. Regions of the membrane with a high correlation between ordered lipids are identified by employing the Getis-Ord local spatial autocorrelation statistic on a Voronoi tessellation of the lipids. As an example, the approach is applied to two distinct systems at a coarse-grained resolution, demonstrating either a strong tendency towards phase separation (1,2-dipalmitoyl-sn-glycero-3-phosphocholine (DPPC), 1,2-dilinoleoyl-sn-glycero-3-phosphocholine (DIPC), cholesterol) or a weak tendency toward phase separation (1-palmitoyl-2-oleoyl-sn-glycero-3-phosphocholine (POPC), 1-palmitoyl-2-docosahexaenoyl-sn-glycero-3-phosphocholine (PUPC), cholesterol). Explanations of the steps are complemented by coding examples written in Python, providing both a comprehensive understanding and practical guidance for a seamless integration of the workflow into individual projects.

Role of phosphatidic acid lipids on plasma membrane association of the Ebola virus matrix protein VP40

The Ebola virus matrix protein VP40 is responsible for the formation of the viral matrix by localizing at the inner leaflet of the human plasma membrane (PM). Various lipid types, including PI(4,5)P2 (i.e. PIP2) and phosphatidylserine (PS), play active roles in this process. Specifically, the negatively charged headgroups of both PIP2 and PS interact with the basic residues of VP40 and stabilize it at the membrane surface, allowing for eventual egress. Phosphatidic acid (PA), resulting from the enzyme phospholipase D (PLD), is also known to play an active role in viral development. In this work, we performed a biophysical and computational analysis to investigate the effects of the presence of PA on the membrane localization and association of VP40. We used coarse-grained molecular dynamics simulations to quantify VP40 hexamer interactions with the inner leaflet of the PM. Analysis of the local distribution of lipids shows enhanced lipid clustering when PA is abundant in the membrane. We observed that PA lipids have a similar role to that of PS lipids in VP40 association due to the geometry and charge. Complementary experiments performed in cell culture demonstrate competition between VP40 and a canonical PA-binding protein for the PM. Also, inhibition of PA synthesis reduced the detectable budding of virus-like particles. These computational and experimental results provide new insights into the early stages of Ebola virus budding and the role that PA lipids have on the VP40-PM association.

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.

Towards a Unitary Hypothesis of Alzheimer's Disease Pathogenesis

The "amyloid cascade" hypothesis of Alzheimer's disease (AD) pathogenesis invokes the accumulation in the brain of plaques (containing the amyloid-β protein precursor [AβPP] cleavage product amyloid-β [Aβ]) and tangles (containing hyperphosphorylated tau) as drivers of pathogenesis. However, the poor track record of clinical trials based on this hypothesis suggests that the accumulation of these peptides is not the only cause of AD. Here, an alternative hypothesis is proposed in which the AβPP cleavage product C99, not Aβ, is the main culprit, via its role as a regulator of cholesterol metabolism. C99, which is a cholesterol sensor, promotes the formation of mitochondria-associated endoplasmic reticulum (ER) membranes (MAM), a cholesterol-rich lipid raft-like subdomain of the ER that communicates, both physically and biochemically, with mitochondria. We propose that in early-onset AD (EOAD), MAM-localized C99 is elevated above normal levels, resulting in increased transport of cholesterol from the plasma membrane to membranes of intracellular organelles, such as ER/endosomes, thereby upregulating MAM function and driving pathology. By the same token, late-onset AD (LOAD) is triggered by any genetic variant that increases the accumulation of intracellular cholesterol that, in turn, boosts the levels of C99 and again upregulates MAM function. Thus, the functional cause of AD is upregulated MAM function that, in turn, causes the hallmark disease phenotypes, including the plaques and tangles. Accordingly, the MAM hypothesis invokes two key interrelated elements, C99 and cholesterol, that converge at the MAM to drive AD pathogenesis. From this perspective, AD is, at bottom, a lipid disorder.

Insights into the Surface Binding and Structural Interference of Polyphenols with the Membrane Raft Domains in Relation to Their Distinctive Ability to Inhibit Preadipocyte Differentiation in 3T3-L1 Cells

Polyphenols with different structures have shown distinct variations in their ability to inhibit the differentiation of 3T3-L1 preadipocytes. However, the underlying mechanisms for these differences remain unclear. In the present study, the surface binding of polyphenols to different membrane domains was explored using coarse-grained molecular dynamics simulation (CG-MDs). Subsequently, this surface binding was confirmed in the liposome system by microscale thermophoresis. Additionally, the interference of polyphenols on the membrane raft's structure was studied through atomic force microscopy and high-content screening fluorescence microscopy. The results indicated that polyphenols with a differentiation-inhibitory ability, such as epicatechin-3-gallate (ECG) and epicatechin-3-gallate-(4β → 8, 2β → O → 7)-epicatechin-3-gallate (A-type ECG dimer), exhibited strong binding to ordered domains enriched in sphingolipids and cholesterol. This binding led to the structural disruption of membrane rafts by altering their size and shape, with the binding constant of 3.8 μM for ECG and 0.3 μM for A-type ECG dimer, respectively. In contrast, epicatechin (EC) with little differentiation-inhibitory ability had no effects on membrane rafts, and its binding constant with the ordered domain was 380.6 μM. Overall, the surface binding of polyphenols to ordered domains and the resulting disruption of membrane rafts structure might be a fundamental mechanism by which polyphenols inhibited the differentiation of 3T3-L1 preadipocytes.

Seeing the Membrane from Both Sides Now: Lipid Asymmetry and Its Strange Consequences

Almost all biomembranes are constructed as lipid bilayers and, in almost all of these, the two opposing monolayers (leaflets) have distinct lipid compositions. This lipid asymmetry arises through the concerted action of a suite of energy-dependent enzymes that maintain living bilayers in a far-from-equilibrium steady-state. Recent discoveries reveal that lipid compositional asymmetry imparts biophysical asymmetries and that this dualistic organization may have major consequences for cellular physiology. Importantly, while transbilayer asymmetry appears to be an essential, near-ubiquitous characteristic of biological membranes, it has been challenging to reproduce in reconstituted or synthetic systems. Although recent methodological developments have overcome some critical challenges, it remains difficult to extrapolate results from available models to biological systems. Concurrently, there are few experimental approaches for targeted, controlled manipulation of lipid asymmetry in living cells. Thus, the biophysical and functional consequences of membrane asymmetry remain almost wholly unexplored. This perspective summarizes the current state of knowledge and highlights emerging themes that are beginning to make inroads into the fundamental question of why life tends toward asymmetry in its bilayers.

CGCompiler: Automated Coarse-Grained Molecule Parametrization via Noise-Resistant Mixed-Variable Optimization

Coarse-grained force fields (CG FFs) such as the Martini model entail a predefined, fixed set of Lennard-Jones parameters (building blocks) to model virtually all possible nonbonded interactions between chemically relevant molecules. Owing to its universality and transferability, the building-block coarse-grained approach has gained tremendous popularity over the past decade. The parametrization of molecules can be highly complex and often involves the selection and fine-tuning of a large number of parameters (e.g., bead types and bond lengths) to optimally match multiple relevant targets simultaneously. The parametrization of a molecule within the building-block CG approach is a mixed-variable optimization problem: the nonbonded interactions are discrete variables, whereas the bonded interactions are continuous variables. Here, we pioneer the utility of mixed-variable particle swarm optimization in automatically parametrizing molecules within the Martini 3 coarse-grained force field by matching both structural (e.g., RDFs) as well as thermodynamic data (phase-transition temperatures). For the sake of demonstration, we parametrize the linker of the lipid sphingomyelin. The important advantage of our approach is that both bonded and nonbonded interactions are simultaneously optimized while conserving the search efficiency of vector guided particle swarm optimization (PSO) methods over other metaheuristic search methods such as genetic algorithms. In addition, we explore noise-mitigation strategies in matching the phase-transition temperatures of lipid membranes, where nucleation and concomitant hysteresis introduce a dominant noise term within the objective function. We propose that noise-resistant mixed-variable PSO methods can both improve and automate parametrization of molecules within building-block CG FFs, such as Martini.

Anomalous lateral diffusion of lipids during the fluid/gel phase transition of a lipid membrane

A lipid membrane undergoes a phase transition from fluid to gel phase upon changing external thermodynamic conditions, such as decreasing temperature and increasing pressure. Extremophilic organisms face the challenge of preventing this deleterious phase transition. The main focus of their adaptive strategy is to facilitate effective temperature sensing through sensor proteins, relying on the drastic changes in packing density and membrane fluidity during the phase transition. Although the changes in packing density parameters due to the fluid/gel phase transition are studied in detail, the impact on membrane fluidity is less explored in the literature. Understanding the lateral diffusive dynamics of lipids in response to temperature, particularly during the fluid/gel phase transition, is albeit crucial. Here we have simulated the phase transition of a single component lipid membrane composed of dipalmitoylphosphatidylcholine (DPPC) lipids using a coarse-grained (CG) model and studied the changes of the structural and dynamical properties. It is observed that near the phase transition point, both fluid and gel phase domains coexist together. The dynamics remains highly non-Gaussian for a long time even when the mean square displacement reaches the Fickian regime at a much earlier time. This Fickian yet non-Gaussian diffusion (FnGD) is a characteristic of a highly heterogeneous system, previously observed for the lateral diffusion of lipids in raft mimetic membranes having liquid-ordered and liquid-disordered phases co-existing together. We have analyzed the molecular trajectories and calculated the jump-diffusion of the lipids, stemming from sudden jump translations, using a translational jump-diffusion (TJD) approach. An overwhelming contribution of the jump-diffusion of the lipids is observed suggesting anomalous diffusion of lipids during fluid/gel phase transition of the membrane. These results are important in unravelling the intricate nature of lipid diffusion during the phase transition of the membrane and open up a new possibility of investigating the most significant change of membrane properties during phase transition, which can be effectively sensed by proteins.

Experimentally determined leaflet-leaflet phase diagram of an asymmetric lipid bilayer

We have determined the partial leaflet-leaflet phase diagram of an asymmetric lipid bilayer at ambient temperature using asymmetric giant unilamellar vesicles (aGUVs). Symmetric GUVs with varying amounts of 1,2-dipalmitoyl-sn-glycero-3-phosphocholine and DOPC (1,2-dioleoyl-sn-glycero-3-phosphocholine) were hemifused to a supported lipid bilayer (SLB) composed of DOPC, resulting in lipid exchange between their outer leaflets. The GUVs and SLB contained a red and green lipid fluorophore, respectively, thus enabling the use of confocal fluorescence imaging to determine both the extent of lipid exchange (quantified for individual vesicles by the loss of red intensity and gain of green intensity) and the presence or absence of phase separation in aGUVs. Consistent with previous reports, we found that hemifusion results in large variation in outer leaflet exchange for individual GUVs, which allowed us to interrogate the phase behavior at multiple points within the asymmetric composition space of the binary mixture. When initially symmetric GUVs showed coexisting gel and fluid domains, aGUVs with less than ~50% outer leaflet exchange were also phase-separated. In contrast, aGUVs with greater than 50% outer leaflet exchange were uniform and fluid. In some cases, we also observed three coexisting bilayer-spanning phases: two registered phases and an anti-registered phase. These results suggest that a relatively large unfavorable midplane interaction between ordered and disordered phases in opposing leaflets (i.e., a midplane surface tension) can overwhelm the driving force for lateral phase separation within one of the leaflets, resulting in an asymmetric bilayer with two uniformly mixed leaflets that is poised to phase-separate upon leaflet scrambling.

Understanding the free-energy landscape of phase separation in lipid bilayers using molecular dynamics

Liquid-liquid phase separation inside the cell often results in biological condensates that can critically affect cell homeostasis. Such phase separation events occur in multiple parts of cells, including the cell membranes, where the "lipid raft" hypothesis posits the formation of ordered domains floating in a sea of disordered lipids. The resulting lipid domains often have functional roles. However, the thermodynamics of lipid phase separation and their resulting mechanistic effects on cell function and dysfunction are poorly understood. Understanding such complex phenomena in cell membranes, with their diverse lipid compositions, is exceptionally difficult. For these reasons, simple model systems that can recapitulate similar behavior are widely used to study this phenomenon. Despite these simplifications, the timescale and length scales of domain formation pose a challenge for molecular dynamics (MD) simulations. Thus, most MD studies focus on spontaneous lipid phase separation-essentially measuring the sign (but not the amplitude) of the free-energy change upon separation-rather than directly interrogating the thermodynamics. Here, we propose a proof-of-concept pipeline that can directly measure this free energy by combining coarse-grained MD with enhanced sampling protocols using a novel collective variable. This approach will be a useful tool to help connect the thermodynamics of phase separation with the mechanistic insights already available from MD simulations.

Investigation of nano- and microdomains formed by ceramide 1 phosphate in lipid bilayers

Biological membranes are renowned for their intricate complexity, with the formation of membrane domains being pivotal to the successful execution of numerous cellular processes. However, due to their nanoscale characteristics, these domains are often understudied, as the experimental techniques required for quantitative investigation present significant challenges. In this study we employ spot-variation z-scan fluorescence correlation spectroscopy (svzFCS) tailored for artificial lipid vesicles of varying composition and combine this approach with high-resolution imaging. This method has been harnessed to examine the lipid-segregation behavior of distinct types of ceramide-1-phosphate (C1P), a crucial class of signaling molecules, within these membranes. Moreover, we provide a quantitative portrayal of the lipid membranes studied and the domains induced by C1P at both nano and microscales. Given the lack of definitive conclusions from the experimental data obtained, it was supplemented with comprehensive in silico studies-including the analysis of diffusion coefficient via molecular dynamics and domain populations via Monte Carlo simulations. This approach enhanced our insight into the dynamic behavior of these molecules within model lipid membranes, confirming that nano- and microdomains can co-exist in lipid vesicles.

The Membrane-Mediated Interaction of Liquid-Ordered Lipid Domains in the Presence of Amphipathic Peptides

The lipid membranes of living cells are composed of a large number of lipid types and can undergo phase separation with the formation of nanometer-scale liquid-ordered lipid domains, also called rafts. Raft coalescence, i.e., the fusion of lipid domains, is involved in important cell processes, such as signaling and trafficking. In this work, within the framework of the theory of elasticity of lipid membranes, we explore how amphipathic peptides adsorbed on lipid membranes may affect the domain-domain fusion processes. We show that the elastic deformations of lipid membranes drive amphipathic peptides to the boundary of lipid domains, which leads to an increase in the average energy barrier of the domain-domain fusion, even if the surface concentration of amphipathic peptides is low and the domain boundaries are only partially occupied by the peptides. This inhibition of the fusion of lipid domains may lead to negative side effects of using amphipathic peptides as antimicrobial agents.

Understanding the Free Energy Landscape of Phase Separation in Lipid Bilayers using Molecular Dynamics

Liquid-liquid phase separation (LLPS) inside the cell often results in biological condensates that can critically impact cell homeostasis. Such phase separation events occur in multiple parts of cells, including the cell membranes, where the so-called "lipid raft" hypothesis posits the formation of ordered domains floating in a sea of disordered lipids. The resulting lipid domains often have functional roles. However, the thermodynamics of lipid phase separation and their resulting mechanistic effects on cell function and dysfunction are poorly understood. Understanding such complex phenomena in cell membranes, with their diverse lipid compositions, is exceptionally difficult. For this reasons, simple model systems that can recapitulate similar behavior are widely used to study this phenomenon. Despite these simplifications, the timescale and and length scales of domain formation pose a challenge for molecular dynamics (MD) simulations. Thus, most MD studies focus on spontaneous lipid phase separation - essentially measuring the sign (but not the amplitude) of the free energy change upon separation - rather than directly interrogating the thermodynamics. Here, we propose a proof-of-concept pipeline that can directly measure this free energy by combining coarse-grained MD with enhanced sampling protocols using a novel collective variable. This approach will be a useful tool to help connect the thermodynamics of phase separation with the mechanistic insights already available from molecular dynamics simulations.

SIGNIFICANCE

Standard molecular dynamics simulations can determine the sign the free energy change upon phase separation, but not the amplitude. We present a new method to determine the phase separation free energy for lipid membranes, based on a enhanced sampling using the weighted ensemble method combined with a novel collective variable, validated using coarse-grained simulations applied to several simple systems. The new method will be valuable as a way to develop models that connect molecular-level structural features to the thermodynamics of phase separation.

Membrane Domain Anti-Registration Induces an Intrinsic Transmembrane Potential

Plasma membrane segregation into various nanoscale membrane domains is driven by distinct interactions between diverse lipids and proteins. Among them, liquid-ordered (Lo) membrane domains are defined as "lipid rafts" and liquid-disordered (Ld) ones as "lipid non-rafts". Using model membrane systems, both intra-leaflet and inter-leaflet dynamics of these membrane domains are widely studied. Nevertheless, the biological impact of the latter, which is accompanied by membrane domain registration/anti-registration, is far from clear. Hence, in this work, we studied the biological relevance of the membrane domain anti-registration using both all-atom molecular dynamics (MD) simulations and confocal fluorescence microscopy. All-atom MD simulations suggested an intrinsic transmembrane potential for the case of the membrane anti-registration (Lo/Ld). Meanwhile, confocal fluorescence microscopy experiments of HeLa and 293T cell lines indicated that membrane cholesterol depletion could significantly alter the transmembrane potential of cells. Considering differences in the cholesterol content between Lo and Ld membrane domains, our confocal fluorescence microscopy experiments are consistent with our all-atom MD simulations. In short, membrane domain anti-registration induces local membrane asymmetry and, thus, an intrinsic transmembrane potential.

Energy Transfer Characteristics of Lipid Bilayer Membranes of Liposomes Examined with Picosecond Time-Resolved Raman Spectroscopy

A number of biochemical reactions proceed inside biomembranes. Since the rate of a chemical reaction is influenced by chemical properties of the surrounding environment, it is important to examine the chemical environment inside the biomembranes. Although the energy transfer characteristics are a basic and important property of a reaction medium, experimental investigation of the thermal conducting capabilities of the biomembranes is a challenging task. We have examined the energy transfer characteristics of lipid bilayer membranes of liposomes, a good model system for the biomembrane, with picosecond time-resolved Raman spectroscopy. The cooling kinetics of the first excited singlet (S1) state of trans-stilbene solubilized within the lipid bilayer membranes is observed as a peak shift of the 1570 cm-1 Raman band of S1 trans-stilbene. The cooling rate constant of S1 trans-stilbene is obtained in six lipid bilayer membranes formed by phospholipids with different hydrocarbon chains, DSPC, DPPC, DMPC, DLPC, DOPC, and egg-PC. We estimate the thermal diffusivity of the lipid bilayer membranes with a known correlation between the cooling rate constant and the thermal diffusivity of the solvent. The thermal diffusivity estimated for the liquid-crystal-phase lipid bilayer membranes is 8.9 × 10-8 to 9.4 × 10-8 m2 s-1, while that for the gel-phase lipid bilayer membranes is 8.4 × 10-8 to 8.5 × 10-8 m2 s-1. The difference in thermal diffusivity between the two phases is explained by a one-dimensional diffusion equation of heat.

Utilizing Machine Learning to Greatly Expand the Range and Accuracy of Bottom-Up Coarse-Grained Models through Virtual Particles

Coarse-grained (CG) models parametrized using atomistic reference data, i.e., "bottom up" CG models, have proven useful in the study of biomolecules and other soft matter. However, the construction of highly accurate, low resolution CG models of biomolecules remains challenging. We demonstrate in this work how virtual particles, CG sites with no atomistic correspondence, can be incorporated into CG models within the context of relative entropy minimization (REM) as latent variables. The methodology presented, variational derivative relative entropy minimization (VD-REM), enables optimization of virtual particle interactions through a gradient descent algorithm aided by machine learning. We apply this methodology to the challenging case of a solvent-free CG model of a 1,2-dioleoyl-sn-glycero-3-phosphocholine (DOPC) lipid bilayer and demonstrate that introduction of virtual particles captures solvent-mediated behavior and higher-order correlations which REM alone cannot capture in a more standard CG model based only on the mapping of collections of atoms to the CG sites.

Membrane-Active Peptides Attack Cell Membranes in a Lipid-Regulated Curvature-Generating Mode

Membrane-active peptides (MAPs) exhibit great potential in biomedical applications due to their unique ability to overcome the cell membrane barrier. However, the interactions between MAPs and membranes are complex, and little is known about the possibility of MAP action being specific to certain types of membranes. In this study, a combination of molecular dynamics simulations and theoretical analysis was utilized to investigate the interactions between typical MAPs and realistic cell membrane systems. Remarkably, the simulations revealed that MAPs can attack membranes by generating and sensing positive mean curvature, which is dependent on lipid composition. Furthermore, theoretical calculations demonstrated that this lipid-regulated curvature-based membrane attack mechanism is an integrated result of multiple effects, including peptide-induced membrane wedge and softening effects, the lipid shape effect, the area-difference elastic effect, and the boundary edge effect of formed peptide-lipid nanodomains. This study enhances our comprehension of MAP-membrane interactions and highlights the potential for developing membrane-specific MAP-based agents.

Efficient calculation of the free energy for protein partitioning using restraining potentials

An approach for the efficient simulation of phase-separated lipid bilayers, for use in the calculation of equilibrium free energies of partitioning between lipid domains, is proposed. The methodology exploits restraint potentials and rectangular aspect ratios that enforce lipid phase separation, allowing for the simulation of smaller systems that approximately reproduce bulk behavior. The utility of this approach is demonstrated through the calculation of potentials of mean force for the translation of a transmembrane protein between lipid domains. The impact of the imposed restraints on lipid tail ordering and lipid packing are explored, providing insight into how restraints can best be employed to compute accurate free-energy surfaces. This approach should be useful in the accurate calculation of equilibrium partition coefficients for transmembrane protein partitioning in heterogeneous membranes, providing insight into the thermodynamic driving forces that control this fundamental biophysical phenomenon.

Exploring the Role of Anionic Lipid Nanodomains in the Membrane Disruption and Protein Folding of Human Islet Amyloid Polypeptide Oligomers on Lipid Membrane Surfaces Using Multiscale Molecular Dynamics Simulations

The aggregation of human Islet Amyloid Polypeptide (hIAPP) on cell membranes is linked to amyloid diseases. However, the physio-chemical mechanisms of how these hIAPP aggregates trigger membrane damage are unclear. Using coarse-grained and all-atom molecular dynamics simulations, we investigated the role of lipid nanodomains in the presence or absence of anionic lipids, phosphatidylserine (PS), and a ganglioside (GM1), in the membrane disruption and protein folding behaviors of hIAPP aggregates on phase-separated raft membranes. Our raft membranes contain liquid-ordered (Lo), liquid-disordered (Ld), mixed Lo/Ld (Lod), PS-cluster, and GM1-cluster nanosized domains. We observed that hIAPP aggregates bound to the Lod domain in the absence of anionic lipids, but also to the GM1-cluster- and PS-cluster-containing domains, with stronger affinity in the presence of anionic lipids. We discovered that L16 and I26 are the lipid anchoring residues of hIAPP binding to the Lod and PS-cluster domains. Finally, significant lipid acyl chain order disruption in the annular lipid shells surrounding the membrane-bound hIAPP aggregates and protein folding, particularly beta-sheet formation, in larger protein aggregates were evident. We propose that the interactions of hIAPP and both non-anionic and anionic lipid nanodomains represent key molecular events of membrane damage associated with the pathogenesis of amyloid diseases.

Phase separation in a ternary DPPC/DOPC/POPC system with reducing hydration

The maintenance of plasma membrane structure is vital for the viability of cells. Disruption of this structure can lead to cell death. One important example is the macroscopic phase separation observed during dehydration associated with desiccation and freezing, often leading to loss of permeability and cell death. It has previously been shown that the hybrid lipid 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphocholine (POPC) can act as a line-active component in ternary lipid systems, inhibiting macroscopic phase separation and stabilising membrane microdomains in lipid vesicles [1]. The domain size is found to decrease with increasing POPC concentration until complete mixing is observed. However, no such studies have been carried out at reduced hydration. To examine if this phase separation is unique to vesicles in excess water, we have conducted studies on several binary and ternary model membrane systems at both reduced hydration ("powder" type samples and oriented membrane stacks) and in excess water (supported lipid bilayers) at 0.2 mol fraction POPC, in the range where microdomain stabilisation is reported. Differential scanning calorimetry (DSC) and Fourier transform infrared spectroscopy (FTIR) are used to map phase transition temperatures, with X-ray and neutron scattering providing details of the changes in lipid packing and phase information within these boundaries. Atomic force microscopy (AFM) is used to image bilayers on a substrate in excess water. In all cases, macroscopic phase separation was observed rather than microdomain formation at this molar ratio. Thus POPC does not stabilise microdomains under these conditions, regardless of the type of model membrane, hydration or temperature. Thus we conclude that the driving force for separation under these conditions overcomes any linactant effects of the hybrid lipid.

Phosphatidylinositol 4,5-Bisphosphate Sensing Lipid Raft via Inter-Leaflet Coupling Regulated by Acyl Chain Length of Sphingomyelin

Phosphatidylinositol 4,5-bisphosphate (PIP₂) is an important molecule located at the inner leaflet of cell membrane, where it serves as anchoring sites for a cohort of membrane-associated molecules and as a broad-reaching signaling intermediate. The lipid raft is thought as the major platform recruiting proteins for signal transduction and also known to mediate PIP₂ accumulation across the membrane. While the significance of this cross-membrane coupling is increasingly appreciated, it remains unclear whether and how PIP₂ senses the dynamic change of the ordered lipid domains over the packed hydrophobic core of the bilayer. Herein, by means of molecular dynamic simulation, we reveal that inner PIP₂ molecules can sense the outer lipid domain via inter-leaflet coupling, and the coupling manner is dictated by the acyl chain length of sphingomyelin (SM) partitioned to the lipid domain. Shorter SM promotes membrane domain registration, whereby PIP₂ accumulates beneath the domain across the membrane. In contrast, the anti-registration is thermodynamically preferred if the lipid domain has longer SM due to the hydrophobic mismatch between the corresponding acyl chains in SM and PIP₂. In this case, PIP₂ is expelled by the domain with a higher diffusivity. These results provide molecular insights into the regulatory mechanism of correlation between the outer lipid domain and inner PIP₂, both of which are critical components for cell signal transduction.

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.

CHARMM-GUI Membrane Builder: Past, Current, and Future Developments and Applications

Molecular dynamics simulations of membranes and membrane proteins serve as computational microscopes, revealing coordinated events at the membrane interface. As G protein-coupled receptors, ion channels, transporters, and membrane-bound enzymes are important drug targets, understanding their drug binding and action mechanisms in a realistic membrane becomes critical. Advances in materials science and physical chemistry further demand an atomistic understanding of lipid domains and interactions between materials and membranes. Despite a wide range of membrane simulation studies, generating a complex membrane assembly remains challenging. Here, we review the capability of CHARMM-GUI Membrane Builder in the context of emerging research demands, as well as the application examples from the CHARMM-GUI user community, including membrane biophysics, membrane protein drug-binding and dynamics, protein-lipid interactions, and nano-bio interface. We also provide our perspective on future Membrane Builder development.

Phase Transition and Phase Separation in Realistic Thylakoid Lipid Membrane of Marine Algae in All-Atom Simulations

Thylakoid membranes are specialized membranes predominantly composed of uncommon galacto- and sulfolipids, having distinct roles in photosynthesis. Large acyl chain variety and richness in polyunsaturated fatty acid (PUFA) content of thylakoid lipids further add to the compositional complexity. The function of these membrane systems is intimately dependent on the fluidity of its lipid matrix, which is strongly modulated by the lipid composition and temperature. The present work, employing extensive atomistic simulations, provides the first atomistic view of the phase transition and domain coexistence in a model membrane composed of thylakoid lipids of a commercially important red alga Gracilaria corticata between 10 and 40 °C. The growth and photosynthetic activity of marine algae are greatly influenced by the seawater temperature. So far, little is known about the molecular organization of lipids in thylakoid membranes, in particular their adaptive arrangements under temperature stress. Our simulations show that the algal thylakoid membrane undergoes a transition from a gel-like phase at a low temperature, 10-15 °C, to a homogeneous liquid-crystalline phase at a high temperature, 40 °C. Clear evidence of spontaneous phase separation into coexisting nanoscale domains is detected at intermediate temperatures nearing the optimal growth temperature range. Particularly, at 25-30 °C, we identified the formation of a stable ripple phase, where the gel-like domains rich in saturated and nearly hexagonally packed lipids were separated from fluid-like domains enriched in lipids containing PUFA chains. The phase separation is driven by the spontaneous and preferential segregation of lipids into differentially ordered domains, mainly depending on the acyl chain types. Cholesterol impairs the phase transition and the emergence of domains and induces a fairly uniform liquid-ordered phase in the membrane over the temperatures studied. This work improves the understanding of the properties and reorganization of lipids in the thylakoid membrane in response to temperature variation.

Cyclic lipopeptides as membrane fusion inhibitors against SARS-CoV-2: New tricks for old dogs

With the resurgence of the coronavirus pandemic, the repositioning of FDA-approved drugs against coronovirus and finding alternative strategies for antiviral therapy are both important. We previously identified the viral lipid envelope as a potential target for the prevention and treatment of SARS-CoV-2 infection with plant alkaloids (Shekunov et al., 2021). Here, we investigated the effects of eleven cyclic lipopeptides (CLPs), including well-known antifungal and antibacterial compounds, on the liposome fusion triggered by calcium, polyethylene glycol 8000, and a fragment of SARS-CoV-2 fusion peptide (816-827) by calcein release assays. Differential scanning microcalorimetry of the gel-to-liquid-crystalline and lamellar-to-inverted hexagonal phase transitions and confocal fluorescence microscopy demonstrated the relation of the fusion inhibitory effects of CLPs to alterations in lipid packing, membrane curvature stress and domain organization. The antiviral effects of CLPs were evaluated in an in vitro Vero-based cell model, and aculeacin A, anidulafugin, iturin A, and mycosubtilin attenuated the cytopathogenicity of SARS-CoV-2 without specific toxicity.

Coarse-Grained Molecular Dynamics of pH-Sensitive Lipids

pH-sensitive lipids represent a class of lipids that can be protonated and destabilized in acidic environments, as they become positively charged in response to low-pH conditions. They can be incorporated into lipidic nanoparticles such as liposomes, which are able to change their properties and allow specific drug delivery at the acidic conditions encountered in some pathological microenvironments. In this work, we used coarse-grained molecular-dynamic simulations to study the stability of neutral and charged lipid bilayers containing POPC (1-palmitoyl-2-oleoyl-sn-glycero-3-phosphocholine) and various kinds of ISUCA ((F)2-(imidazol-1-yl)succinic acid)-derived lipids, which can act as pH-sensitive molecules. In order to explore such systems, we used a MARTINI-derived forcefield, previously parameterized using all-atom simulation results. We calculated the average area per lipid, the second-rank order parameter and the lipid diffusion coefficient of both lipid bilayers made of pure components and mixtures of lipids in different proportions, under neutral or acidic conditions. The results show that the use of ISUCA-derived lipids disturbs the lipid bilayer structure, with the effect being particularly marked under acidic conditions. Although more-in depth studies on these systems must be carried out, these initial results are encouraging and the lipids designed in this research could be a good basis for developing new pH-sensitive liposomes.

Cofilin-Membrane Interactions: Electrostatic Effects in Phosphoinositide Lipid Binding

The actin cytoskeleton interacts with the cell membrane primarily through the indirect interactions of actin-binding proteins such as cofilin-1. The molecular mechanisms underlying the specific interactions of cofilin-1 with membrane lipids are still unclear. Here, we performed coarse-grain molecular dynamics simulations of cofilin-1 with complex lipid bilayers to analyze the specificity of protein-lipid interactions. We observed the maximal interactions with phosphoinositide (PIP) lipids, especially PIP₂ and PIP₃ lipids. A good match was observed between the residues predicted to interact and previous experimental studies. The clustering of PIP lipids around the membrane bound protein leads to an overall lipid demixing and gives rise to persistent membrane curvature. Further, through a series of control simulations, we observe that both electrostatics and geometry are critical for specificity of lipid binding. Our current study is a step towards understanding the physico-chemical basis of cofilin-PIP lipid interactions.

In Vitro Membrane Platform for the Visualization of Water Impermeability across the Liquid-Ordered Phase under Hypertonic Conditions

Passive water penetration across the cell membrane by osmotic diffusion is essential for the homeostasis of cell volume, in addition to the protein-assisted active transportation of water. Since membrane components can regulate water permeability, controlling compositional variation during the volume regulatory process is a prerequisite for investigating the underlying mechanisms of water permeation and related membrane dynamics. However, the lack of a viable in vitro membrane platform in hypertonic solutions impedes advanced knowledge of cell volume regulation processes, especially cholesterol-enriched lipid domains called lipid rafts. By reconstituting the liquid-ordered (L_o) domain as a likeness of lipid rafts, we verified suppressed water permeation across the L_o domains, which had yet to be confirmed with experimental demonstrations despite a simulation approach. With the help of direct transfer of the L_o domains from vesicles to supported lipid membranes, the biological roles of lipid composition in suppressed water translocation were experimentally confirmed. Additionally, the improvement in membrane stability under hypertonic conditions was demonstrated based on molecular dynamics simulations.

Exploring Membrane Binding Targets of Disordered Human Tau Aggregates on Lipid Rafts Using Multiscale Molecular Dynamics Simulations

The self-aggregation of tau, a microtubule-binding protein, has been linked to the onset of Alzheimer's Disease. Recent studies indicate that the disordered tau aggregates, or oligomers, are more toxic than the ordered fibrils found in the intracellular neurofibrillary tangles of tau. At present, details of tau oligomer interactions with lipid rafts, a model of neuronal membranes, are not known. Using molecular dynamics simulations, the lipid-binding events, membrane-damage, and protein folding of tau oligomers on various lipid raft surfaces were investigated. Tau oligomers preferred to bind to the boundary domains (Lod) created by the coexisting liquid-ordered (Lo) and liquid-disordered (Ld) domains in the lipid rafts. Additionally, stronger binding of tau oligomers to the ganglioside (GM1) and phosphatidylserine (PS) domains, and subsequent protein-induced lipid chain order disruption and beta-sheet formation were detected. Our results suggest that GM1 and PS domains, located exclusively in the outer and inner leaflets, respectively, of the neuronal membranes, are specific membrane domain targets, whereas the Lod domains are non-specific targets, of tau oligomers binding to neurons. The molecular details of these specific and non-specific tau bindings to lipid rafts may provide new insights into understanding membrane-associated tauopathies leading to Alzheimer's Disease.

Change in membrane fluidity induced by polyphenols is highly dependent on the position and number of galloyl groups

The cell membrane fluidity was very important in adipogenesis and galloyl groups on polyphenolic structures could enhance their antiadipogenic activity. However, the effect of polyphenols on membrane fluidity and the role of galloyl groups in fluidity changes remain unclear. Therefore, the present study chose structurally different polyphenols to compare their effects on the membrane morphology and fluidity of 3T3-L1 preadipocytes, and then the reasons behind the changes of membrane fluidity induced by galloylated polyphenols were explored from structural and molecular insights using liposome model and molecular dynamic simulation technology. Our results indicated that galloylated polyphenols could significantly change 3T3-L1 cell membrane morphology and decrease membrane fluidity, while non-galloylated ones could not. The membrane interference effect of polyphenols was enhanced as the number of galloyl groups increased. Morever, the decrease in membrane fluidity induced by galloylated polyphenols was due to the disturbance of polyphenols on lipid alkyl chains in the cell membrane. Galloylated polyphenols could not only locate in the polar head, but also insert into hydrophobic center of lipid bilayer to interfere with the lipid alkyl chains arrangement, thus decreasing the membrane fluidity and showing strong affinity for the membrane. In addition, differences in position of galloyl groups in polyphenols induced distinct effect on cell membranes interactions, thus affecting the binding manner and bioactivity. The results expanded the understanding on the strong antiadipogenic activity of galloylated polyphenols through the aspect of their effects on cell membrane by both experimental and theoretically simulated ways.

Formation of lipid raft nanodomains in homogeneous ternary lipid mixture of POPC/DPSM/cholesterol: Theoretical insights

Lipid rafts, in biological membranes, are cholesterol-rich nanodomains that regulate many protein activities and cellular processes. Understanding the formation of the lipid-raft nanodomains helps us elucidate many complex interactions in the cell. In this study, the formation of lipid-raft nanodomains in a ternary palmitoyl-oleoyl-phosphatidylcholine/stearoyl-sphingomyelin/cholesterol (POPC/DPSM/Chol) lipid mixture, the most realistic surrogate model for biological membranes, has been successfully observed for the first time in-silico using microsecond timescale molecular dynamics simulations. The model reveals the formation of cholesterol-induced nanodomains with raft-like characteristics and their underlying mechanism: the cholesterol molecules segregate themselves into cholesterol nanodomains and then enrich the cholesterol-rich domain with sphingomyelin molecules to form a lipid raft thanks to the weak bonding of cholesterol with sphingomyelin. Besides, it is found that the increase in cholesterol concentration enhances the biophysical properties (e.g., bilayer thickness, area per lipid headgroup, and order parameter) of the lipid raft nanodomains. Such findings suggest that the POPC/DPSM/Chol bilayer is a suitable model to fundamentally extend the nanodomain evolution to investigate their lifetime and kinetics as well as to study protein-lipid interaction, protein-protein interaction, and selection of therapeutic molecules in the presence of lipid rafts.

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.

Lysolipids regulate raft size distribution

The lipid matrix of cellular membranes, directly and indirectly, regulates many vital functions of the cell. The diversity of lipids in membranes leads to the formation of ordered domains called rafts, which play a crucial role in signal transduction, protein sorting and other cellular processes. Rafts are believed to impact the development of different neurodegenerative diseases, such as Alzheimer’s, Parkinson’s, Huntington’s ones, amyotrophic lateral sclerosis, some types of cancer, etc. These diseases correlate with the change in the membrane lipid composition resulting from an oxidative stress, age-related processes, dysfunction of proteins, and many others. In particular, a lot of studies report a significant rise in the level of lysolipids. Physicochemical properties of rafts are determined by membrane composition, in particular, by the content of lysolipids. Lysolipids may thus regulate raft-involving processes. However, the exact mechanism of such regulation is unknown. Although studying rafts in vivo still seems to be rather complicated, liquid-ordered domains are well observed in model systems. In the present study, we used atomic force microscopy (AFM) to examine how lysophospholipids influence the liquid-ordered domains in model ternary membranes. We demonstrated that even a small amount of lysolipids in a membrane significantly impacts domain size depending on the saturation of the lysolipid hydrocarbon tails and the amount of cholesterol. The mixture with the bigger relative fraction of cholesterol was more susceptible to the action of lysolipids. This data helped us to generalize our previous theoretical model of the domain size regulation by lipids with particular molecular shape expanding it to the case of lysolipids and dioleoylglycerol.

Homeoviscous Adaptation of the Lipid Membrane of a Soil Bacterium Surviving under Diurnal Temperature Variation: A Molecular Simulation Perspective

A recent experiment has reported the lipidome remodeling of a soil-based plant-associated bacterium Methylobacterium extorquens due to diurnal temperature variations. The key adaptation strategy is the headgroup-specific remodeling of the acyl chain. To understand the idiosyncratic adaptation at the molecular level, we simulate the model membrane of the same bacterium using the reported lipidome compositions at four different experimental temperatures. We investigate the temperature-dependent packing density and fluidity of the membrane, the constancy of which is key to the homeoviscous adaptation. The results show that complex lipidome remodeling approximately preserves membrane properties under heat and cold stress. The headgroup-specific remodeling of the acyl chain serves to fine-tune the packing density and fluidity of the membrane at different temperatures. While lipids with strongly interacting headgroups are more abundant at higher temperatures, the lipidome is more dominated by lipids with weaker interacting headgroups at lower temperatures. This adaptation alleviates lipid membrane disruption caused by heat and cold stress. This study provides a molecular picture of the homeoviscous adaptation of the realistic lipid membrane of a soil-based bacterium.

Fluid-gel coexistence in lipid membranes under differential stress

A widely conserved property of many biological lipid bilayers is their asymmetry. In addition to having distinct compositions on its two sides, a membrane can also exhibit different tensions in its two leaflets, a state known as differential stress. Here, we examine how this stress can influence the phase behavior of the constituent lipid monolayers of a single-component membrane. For temperatures sufficiently close to, but still above, the main transition, molecular dynamics simulations show the emergence of finite gel domains within the compressed leaflet. We describe the thermodynamics of this phenomenon by adding two empirical single-leaflet free energies for the fluid-gel transition, each evaluated at its respective asymmetry-dependent lipid density. Finite size effects arising in simulation are included in the theory through a geometry-dependent interfacial term. Our model reproduces the phase coexistence observed in simulation. It could therefore be used to connect the "hidden variable" of differential stress to experimentally observable properties of the main phase transition. These ideas could be generalized to any first-order bilayer phase transition in the presence of asymmetry, including liquid-ordered/liquid-disordered phase separation.