Domain formation in the plasma membrane: roles of nonequilibrium lipid transport and membrane proteins

Heterogeneous biological membranes regulate protein partitioning via fluctuating diffusivity

Cell membranes phase separate into ordered L o and disordered L d domains depending on their compositions. This membrane compartmentalization is heterogeneous and regulates the localization of specific proteins related to cell signaling and trafficking. However, it is unclear how the heterogeneity of the membranes affects the diffusion and localization of proteins in L o and L d domains. Here, using Langevin dynamics simulations coupled with the phase-field (LDPF) method, we investigate several tens of milliseconds-scale diffusion and localization of proteins in heterogeneous biological membrane models showing phase separation into L o and L d domains. The diffusivity of proteins exhibits temporal fluctuations depending on the field composition. Increases in molecular concentrations and domain preference of the molecule induce subdiffusive behavior due to molecular collisions by crowding and confinement effects, respectively. Moreover, we quantitatively demonstrate that the protein partitioning into the L o domain is determined by the difference in molecular diffusivity between domains, molecular preference of domain, and molecular concentration. These results pave the way for understanding how biological reactions caused by molecular partitioning may be controlled in heterogeneous media. Moreover, the methodology proposed here is applicable not only to biological membrane systems but also to the study of diffusion and localization phenomena of molecules in various heterogeneous systems.

Phase separation on surfaces in the presence of matter exchange

We present a field theory to describe the composition of a surface spontaneously exchanging matter with its bulk environment. By only assuming matter conservation in the system, we show with extensive numerical simulations that, depending on the matter exchange rates, a complex patterned composition distribution emerges on the surface. For one-dimensional systems we show analytically and numerically that coarsening is arrested and as a consequence domains have a characteristic length scale. Our results show that the causes of heterogeneous lipid composition in cellular membranes may be justified in simple physical terms.

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.

Effects of Passive Phospholipid Flip-Flop and Asymmetric External Fields on Bilayer Phase Equilibria

Compositional asymmetry between the leaflets of bilayer membranes modifies their phase behavior and is thought to influence other important features such as mechanical properties and protein activity. We address here how phase behavior is affected by passive phospholipid flip-flop, such that the compositional asymmetry is not fixed. We predict transitions from "pre-flip-flop" behavior to a restricted set of phase equilibria that can persist in the presence of passive flip-flop. Surprisingly, such states are not necessarily symmetric. We further account for external symmetry breaking, such as a preferential substrate interaction, and show how this can stabilize strongly asymmetric equilibrium states. Our theory explains several experimental observations of flip-flop-mediated changes in phase behavior and shows how domain formation and compositional asymmetry can be controlled in concert, by manipulating passive flip-flop rates and applying external fields.

A Rationale for Mesoscopic Domain Formation in Biomembranes

Cell plasma membranes display a dramatically rich structural complexity characterized by functional sub-wavelength domains with specific lipid and protein composition. Under favorable experimental conditions, patterned morphologies can also be observed in vitro on model systems such as supported membranes or lipid vesicles. Lipid mixtures separating in liquid-ordered and liquid-disordered phases below a demixing temperature play a pivotal role in this context. Protein-protein and protein-lipid interactions also contribute to membrane shaping by promoting small domains or clusters. Such phase separations displaying characteristic length-scales falling in-between the nanoscopic, molecular scale on the one hand and the macroscopic scale on the other hand, are named mesophases in soft condensed matter physics. In this review, we propose a classification of the diverse mechanisms leading to mesophase separation in biomembranes. We distinguish between mechanisms relying upon equilibrium thermodynamics and those involving out-of-equilibrium mechanisms, notably active membrane recycling. In equilibrium, we especially focus on the many mechanisms that dwell on an up-down symmetry breaking between the upper and lower bilayer leaflets. Symmetry breaking is an ubiquitous mechanism in condensed matter physics at the heart of several important phenomena. In the present case, it can be either spontaneous (domain buckling) or explicit, i.e., due to an external cause (global or local vesicle bending properties). Whenever possible, theoretical predictions and simulation results are confronted to experiments on model systems or living cells, which enables us to identify the most realistic mechanisms from a biological perspective.

Dynamic Scaling of Exosome Sizes

A model is proposed for characterizing exosome size distributions based on dynamic scaling of domain growth on the limiting membrane of multivesicular bodies in the established exosome biogenesis pathway. The scaling exponent in this model captures the asymmetry of exosome size distributions, which are notably right-skewed to larger vesicles, independent of the minimum detectable vesicle size. Analyses of exosome size distributions obtained by cryogenic transmission electron microscopy imaging and nanoparticle tracking show, respectively, that the scaling exponent is sensitive to the state of the cell source for exosomes in cell culture supernatants and can distinguish exosome size distributions in serum samples taken from cancer patients relative to those from healthy donors. Finally, we comment on mechanistic differences between our dynamic scaling model and random fragmentation models used to describe size distributions of synthetic vesicles.

Impact of lipid rafts on the T-cell-receptor and peptide-major-histocompatibility-complex interactions under different measurement conditions

The interactions between T-cell receptor (TCR) and peptide-major-histocompatibility complex (pMHC), which enable T-cell development and initiate adaptive immune responses, have been intensively studied. However, a central issue of how lipid rafts affect the TCR-pMHC interactions remains unclear. Here, by using a statistical-mechanical membrane model, we show that the binding affinity of TCR and pMHC anchored on two apposing cell membranes is significantly enhanced because of the lipid raft-induced signaling protein aggregation. This finding may provide an alternative insight into the mechanism of T-cell activation triggered by very low densities of pMHC. In the case of cell-substrate adhesion, our results indicate that the loss of lateral mobility of the proteins on the solid substrate leads to the inhibitory effect of lipid rafts on TCR-pMHC interactions. Our findings help to understand why different experimental methods for measuring the impact of lipid rafts on the receptor-ligand interactions have led to contradictory conclusions.

Physical mechanisms of micro- and nanodomain formation in multicomponent lipid membranes

This article summarizes a variety of physical mechanisms proposed in the literature, which can generate micro- and nanodomains in multicomponent lipid bilayers and biomembranes. It mainly focusses on lipid-driven mechanisms that do not involve direct protein-protein interactions. Specifically, it considers (i) equilibrium mechanisms based on lipid-lipid phase separation such as critical cluster formation close to critical points, and multiple domain formation in curved geometries, (ii) equilibrium mechanisms that stabilize two-dimensional microemulsions, such as the effect of linactants and the effect of curvature-composition coupling in bilayers and monolayers, and (iii) non-equilibrium mechanisms induced by the interaction of a biomembrane with the cellular environment, such as membrane recycling and the pinning effects of the cytoplasm. Theoretical predictions are discussed together with simulations and experiments. The presentation is guided by the theory of phase transitions and critical phenomena, and the appendix summarizes the mathematical background in a concise way within the framework of the Ginzburg-Landau theory. This article is part of a Special Issue entitled: Lipid order/lipid defects and lipid-control of protein activity edited by Dirk Schneider.

Nucleation of symmetric domains in the coupled leaflets of a bilayer

We study the kinetics governing the attainment of inter-leaflet domain symmetry in a phase-separating amphiphilic bilayer. "Indirect" inter-leaflet coupling via hydrophobic mismatch can induce an instability towards a metastable pattern of locally asymmetric domains upon quenching from high temperature. This necessitates a nucleation step to form the conventional symmetric domains, which are favoured by a "direct" inter-leaflet coupling. We model the energetics for a symmetric domain to nucleate from the metastable state, and find that an interplay between hydrophobic mismatch and thickness stretching/compression causes the effective hydrophobic mismatch, and thus line tension, to depend on domain size. This leads to strong departure from classical nucleation theory. We speculate on implications for cell membrane rafts or clusters, whose size may be of similar magnitude to estimated critical radii for domain symmetry.

Live cell plasma membranes do not exhibit a miscibility phase transition over a wide range of temperatures

Lipid/cholesterol mixtures derived from cell membranes as well as their synthetic reconstitutions exhibit well-defined miscibility phase transitions and critical phenomena near physiological temperatures. This suggests that lipid/cholesterol-mediated phase separation plays a role in the organization of live cell membranes. However, macroscopic lipid-phase separation is not generally observed in cell membranes, and the degree to which properties of isolated lipid mixtures are preserved in the cell membrane remain unknown. A fundamental property of phase transitions is that the variation of tagged particle diffusion with temperature exhibits an abrupt change as the system passes through the transition, even when the two phases are distributed in a nanometer-scale emulsion. We support this using a variety of Monte Carlo and atomistic simulations on model lipid membrane systems. However, temperature-dependent fluorescence correlation spectroscopy of labeled lipids and membrane-anchored proteins in live cell membranes shows a consistently smooth increase in the diffusion coefficient as a function of temperature. We find no evidence of a discrete miscibility phase transition throughout a wide range of temperatures: 14-37 °C. This contrasts the behavior of giant plasma membrane vesicles (GPMVs) blebbed from the same cells, which do exhibit phase transitions and macroscopic phase separation. Fluorescence lifetime analysis of a DiI probe in both cases reveals a significant environmental difference between the live cell and the GPMV. Taken together, these data suggest the live cell membrane may avoid the miscibility phase transition inherent to its lipid constituents by actively regulating physical parameters, such as tension, in the membrane.

Modeling the interplay between protein and lipid aggregation in supported membranes

We present a theoretical model that deals with the complex interplay between lipid segregation and the self-aggregation of lipid-attached proteins. The model, in contrast to previous ones that consider proteins only as passive elements affecting the lipid distribution, describes the system including three terms: the dynamic interactions between protein monomers, the interactions between lipid components, and a mixed term considering protein-lipid interactions. It is used to explain experimental results performed on a well-defined system in which a self-aggregating soluble bacterial cytoskeletal protein polymerizes on a lipid bilayer containing two lipid components. All the elements considered in a previously described protein model, including torsion of the monomers within the filament, are needed to account for the observed filament shapes. The model also points out that lipid segregation can affect the length and curvature of the filaments and that the dynamic behavior of the lipids and proteins can have different time scales, giving rise to memory effects. This simple model that considers a dynamic protein assembly on a fluid and active lipid surface can be easily extended to other biologically relevant situations in which the interplay between protein and lipid aggregation is needed to fully describe the system.

The molecular structure of the liquid-ordered phase of lipid bilayers

Molecular dynamics simulations reveal substructures within the liquid-ordered phase of lipid bilayers. These substructures, identified in a 10 μs all-atom trajectory of liquid-ordered/liquid-disordered coexistence (L(o)/L(d)) are composed of saturated hydrocarbon chains packed with local hexagonal order and separated by interstitial regions enriched in cholesterol and unsaturated chains. Lipid hydrocarbon chain order parameters calculated from the L(o) phase are in excellent agreement with (2)H NMR measurements; the local hexagonal packing is also consistent with (1)H-MAS NMR spectra of the L(o) phase, NMR diffusion experiments, and small-angle X-ray and neutron scattering. The balance of cholesterol-rich to local hexagonal order is proposed to control the partitioning of membrane components into the L(o) regions. The latter have been frequently associated with formation of so-called rafts, platforms in the plasma membranes of cells that facilitate interaction between components of signaling pathways.

Sphingolipid domains in the plasma membranes of fibroblasts are not enriched with cholesterol

The plasma membranes of mammalian cells are widely expected to contain domains that are enriched with cholesterol and sphingolipids. In this work, we have used high-resolution secondary ion mass spectrometry to directly map the distributions of isotope-labeled cholesterol and sphingolipids in the plasma membranes of intact fibroblast cells. Although acute cholesterol depletion reduced sphingolipid domain abundance, cholesterol was evenly distributed throughout the plasma membrane and was not enriched within the sphingolipid domains. Thus, we rule out favorable cholesterol-sphingolipid interactions as dictating plasma membrane organization in fibroblast cells. Because the sphingolipid domains are disrupted by drugs that depolymerize the cells actin cytoskeleton, cholesterol must instead affect the sphingolipid organization via an indirect mechanism that involves the cytoskeleton.

Interaction of virions with membrane glycolipids

Cellular membranes contain various lipids including glycolipids (GLs). The hydrophilic head groups of GLs extend from the membrane into the aqueous environment outside the cell where they act as recognition sites for specific interactions. The first steps of interaction of virions with cells often include contacts with GLs. To clarify the details of such contacts, we have used the total internal reflection fluorescence microscopy to explore the interaction of individual unlabelled virus-like particles (or, more specifically, norovirus protein capsids), which are firmly bound to a lipid bilayer, and fluorescent vesicles containing glycosphingolipids (these lipids form a subclass of GLs). The corresponding binding kinetics were earlier found to be kinetically limited, while the detachment kinetics were logarithmic over a wide range of time. Here, the detachment rate is observed to dramatically decrease with increasing concentration of glycosphingolipids from 1% to 8%. This effect has been analytically explained by using a generic model describing the statistics of bonds in the contact area between a virion and a lipid membrane. Among other factors, the model takes the formation of GL domains into account. Our analysis indicates that in the system under consideration, such domains, if present, have a characteristic size smaller than the contact area between the vesicle and the virus-like particle.

Aggregation on a membrane of particles undergoing active exchange with a reservoir

We investigate the dynamics of clusters made of aggregating particles on a membrane which exchanges particles with a reservoir. Exchanges are driven by chemical reactions which supply energy to the system, leading to the establishment of a non-equilibrium steady state. We predict the distribution of cluster size at steady state. We show in particular that in a regime, that cannot exist at equilibrium, the distribution is bimodal: the membrane is mainly populated of single particles and finite-size clusters. This work is motivated by the observations that have revealed the existence of submicrometric clusters of proteins in biological membranes.

Minimal model of plasma membrane heterogeneity requires coupling cortical actin to criticality

We present a minimal model of plasma membrane heterogeneity that combines criticality with connectivity to cortical cytoskeleton. The development of this model was motivated by recent observations of micron-sized critical fluctuations in plasma membrane vesicles that are detached from their cortical cytoskeleton. We incorporate criticality using a conserved order parameter Ising model coupled to a simple actin cytoskeleton interacting through point-like pinning sites. Using this minimal model, we recapitulate several experimental observations of plasma membrane raft heterogeneity. Small (r ∼ 20 nm) and dynamic fluctuations at physiological temperatures arise from criticality. Including connectivity to the cortical cytoskeleton disrupts large fluctuations, prevents macroscopic phase separation at low temperatures (T ≤ 22°C), and provides a template for long-lived fluctuations at physiological temperature (T = 37°C). Cytoskeleton-stabilized fluctuations produce significant barriers to the diffusion of some membrane components in a manner that is weakly dependent on the number of pinning sites and strongly dependent on criticality. More generally, we demonstrate that critical fluctuations provide a physical mechanism for organizing and spatially segregating membrane components by providing channels for interaction over large distances.

Monodisperse domains by proteolytic control of the coarsening instability

The coarsening instability typically disrupts steady-state cluster-size distributions. We show that degradation coupled to the cluster size, such as arising from biological proteolysis, leads to a fixed-point cluster size. Stochastic evaporative and condensative fluxes determine the width of the fixed-point size distribution. At the fixed point, we show how the peak size and width depend on number, interactions, and proteolytic rate. This proteolytic size-control mechanism is consistent with the phenomenology of pseudopilus length control in the general secretion pathway of bacteria.

Nanodomain stabilization dynamics in plasma membranes of biological cells

We discover that a synergistically amplifying role of stabilizing membrane proteins and continuous lipid recycling can explain the physics governing the stability, polydispersity, and dynamics of lipid raft domains in plasma membranes of biological cells. We establish the conjecture using a generalized order parameter based on theoretical formalism, endorsed by detailed scaling arguments and domain mapping. Quantitative agreements with morphological distributions of raft complexes, as obtained from Förster resonance energy transfer based visualization, support the present theoretical conjecture.

Membrane composition-mediated protein-protein interactions

The authors investigate membrane composition-mediated interactions between proteins adsorbed onto a two-component lipid bilayer close to critical demixing using coarse-grained molecular dynamics simulations and a phenomenological Ginzburg-Landau theory. The simulations consist of three-bead lipids and platelike proteins, which adsorb onto the membrane by binding preferentially to one of the two lipid species. The composition profile around one protein and the pair correlation function between two proteins are measured and compared to the analytical predictions. The theoretical framework is applicable to any scalar field embedded in the membrane, and although in this work the authors treat flat membranes, the methodology extends readily to curved geometries. Neglecting fluctuations, both lipid composition profile and induced protein pair potential are predicted to follow a zeroth order modified Bessel function of the second kind with the same characteristic decay length. These predictions are consistent with our molecular dynamics simulations, except that the interaction range is found to be larger than the single profile correlation length.

Lipid microdomains: structural correlations, fluctuations, and formation mechanisms

Compositional lipid microdomains ("lipid rafts") in mammalian plasma membranes are believed to facilitate many important cellular processes. While several physically distinct scenarios predicting the presence of finite-sized microdomains in vivo have been proposed in the past, direct experimental verification or falsification of model predictions has remained elusive. Herein, we demonstrate that the combination of the spatial correlation and temporal fluctuation spectra of the lipid domains can be employed to unambiguously differentiate between the existing theoretical scenarios. Furthermore, the differentiation of the raft formation mechanisms using this methodology can be achieved by collecting data at physiologically relevant conditions without the need to tune control parameters.

Influence of nonequilibrium lipid transport, membrane compartmentalization, and membrane proteins on the lateral organization of the plasma membrane

Compositional lipid domains (lipid rafts) in plasma membranes are believed to be important components of many cellular processes. The mechanisms by which cells regulate the sizes, lifetimes, and spatial localization of these domains are rather poorly understood at the moment. We propose a robust mechanism for the formation of finite-sized lipid raft domains in plasma membranes, the competition between phase separation in an immiscible lipid system and active cellular lipid transport processes naturally leads to the formation of such domains. Simulations of a continuum model reveal that the raft size distribution is broad and the average raft size is strongly dependent on the rates of cellular and interlayer lipid transport processes. We demonstrate that spatiotemporal variations in the recycling may enable the cell to localize larger raft aggregates at specific parts along the membrane. Moreover, we show that membrane compartmentalization may further facilitate spatial localization of the raft domains. Finally, we demonstrate that local interactions with immobile membrane proteins can spatially localize the rafts and lead to further clustering.

Formation and regulation of lipid microdomains in cell membranes: theory, modeling, and speculation

Compositional lipid microdomains ("lipid rafts") in plasma membranes are believed to be important components of many cellular processes. The biophysical mechanisms by which cells regulate the size, lifetime, and spatial localization of these domains are rather poorly understood at the moment. Over the years, experimental studies of raft formation have inspired several phenomenological theories and speculations incorporating a wide variety of thermodynamic assumptions regarding lipid-lipid and lipid-protein interactions, and the potential role of active cellular processes on membrane structure. Here we critically review and discuss these theories, models, and speculations, and present our view on future directions.

Critical dynamics in multicomponent lipid membranes

The formation and dynamics of spatially extended compositional domains in multicomponent lipid membranes both in vivo and in vitro lie at the heart of many important biological and biophysical phenomena. While the thermodynamic basis for domain formation has been explored extensively in the past, the roles of membrane and exterior fluid hydrodynamics on domain formation kinetics have received less attention. A case in point is the impact of hydrodynamics on the dynamics of compositional heterogeneities in lipid membranes in the vicinity of a critical point. In this Rapid Communication it is argued that the asymptotic dynamic behavior of a lipid membrane system in the vicinity of a critical point is strongly influenced by hydrodynamic interactions. More specifically, a mode-coupling argument is developed which predicts a scaling behavior of lipid transport coefficients near the critical point for both symmetric and asymmetric bilayers immersed in a bulk fluid.

Stability of protein-decorated mixed lipid membranes: The interplay of lipid-lipid, lipid-protein, and protein-protein interactions

Membrane-associated proteins are likely to contribute to the regulation of the phase behavior of mixed lipid membranes. To gain insight into the underlying mechanism, we study a thermodynamic model for the stability of a protein-decorated binary lipid layer. Here, proteins interact preferentially with one lipid species and thus locally sequester that species. We aim to specify conditions that lead to an additional macroscopic phase separation of the protein-decorated lipid membrane. Our model is based on a standard mean-field lattice-gas description for both the lipid mixture and the adsorbed protein layer. Besides accounting for the lipid-protein binding strength, we also include attractive lipid-lipid and protein-protein interactions. Our analysis characterizes the decrease in the membrane's critical interaction parameter as a function of the lipid-protein binding strength. For small and large binding strengths we provide analytical expressions; numerical results cover the intermediate range. Our results reiterate the crucial importance of the line tension associated with protein-induced compositional gradients and the presence of attractive lipid-lipid interactions within the membrane. Direct protein-protein attraction effectively increases the line tension and thus tends to further destabilize the membrane.

Multiscale modeling of emergent materials: biological and soft matter

In this review, we focus on four current related issues in multiscale modeling of soft and biological matter. First, we discuss how to use structural information from detailed models (or experiments) to construct coarse-grained ones in a hierarchical and systematic way. This is discussed in the context of the so-called Henderson theorem and the inverse Monte Carlo method of Lyubartsev and Laaksonen. In the second part, we take a different look at coarse graining by analyzing conformations of molecules. This is done by the application of self-organizing maps, i.e., a neural network type approach. Such an approach can be used to guide the selection of the relevant degrees of freedom. Then, we discuss technical issues related to the popular dissipative particle dynamics (DPD) method. Importantly, the potentials derived using the inverse Monte Carlo method can be used together with the DPD thermostat. In the final part we focus on solvent-free modeling which offers a different route to coarse graining by integrating out the degrees of freedom associated with solvent.