Dynamic phase separation of fluid membranes with rigid inclusions

Investigating the entropic nature of membrane-mediated interactions driving the aggregation of peripheral proteins

Peripheral membrane-associated proteins are known to accumulate on the surface of biomembranes as a result of membrane-mediated interactions. For a pair of rotationally-symmetric curvature-inducing proteins, membrane mechanics at the low-temperature limit predicts pure repulsion. On the other hand, temperature-dependent entropic forces arise between pairs of stiff-binding proteins suppressing membrane fluctuations. These Casimir-like interactions have thus been suggested as candidates for attractive forces leading to aggregation. With dense assemblies of peripheral proteins on the membrane, both these abstractions encounter short-range and multi-body complications. Here, we make use of a particle-based membrane model augmented with flexible peripheral proteins to quantify purely membrane-mediated interactions and investigate their underlying nature. We introduce a continuous reaction coordinate corresponding to the progression of protein aggregation. We obtain free energy and entropy landscapes for different surface concentrations along this reaction coordinate. In parallel, we investigate time-dependent estimates of membrane entropy corresponding to membrane undulations and coarse-grained director field and how they change dynamically with protein aggregation. Congruent outcomes of the two approaches point to the conclusion that for low surface concentrations, interactions with an entropic nature may drive the aggregation. But at high concentrations, enthalpic contributions due to concerted membrane deformation by protein clusters are dominant.

Thermodynamics and Kinetics of Aggregation of Flexible Peripheral Membrane Proteins

Biomembrane remodeling is essential for cellular trafficking, with membrane-binding peripheral proteins playing a key role in it. Significant membrane remodeling as in endo- and exocytosis is often due to aggregates of many proteins with direct or membrane-mediated interactions. Understanding this process via computer simulations is extremely challenging: protein-membrane systems involve time and length scales that make atomistic simulations impractical, while most coarse-grained models fall short in resolving dynamics and physical effects of protein and membrane flexibility. Here, we develop a coarse-grained model of the bilayer membrane bestrewed with rotationally symmetric flexible proteins, parametrized to reflect local curvatures and lateral dynamics of proteins. We investigate the kinetics, equilibrium distributions, and the free energy landscape governing the formation and breakup of protein clusters on the surface of the membrane. We demonstrate how the flexibility of the proteins as well as their surface concentration play deciding roles in highly selective macroscopic aggregation behavior.

Fluctuations and conformational stability of a membrane patch with curvature inducing inclusions

Membranes with curvature inducing inclusions display a range of cooperative phenomena, which can be linked to biomembrane function, e.g. membrane tubulation, vesiculation, softening and spontaneous tension. We investigate how these phenomena are related for a fluctuating, framed membrane through analysis of a descretized membrane model by Monte Carlo simulation techniques. The membrane model is based on a dynamically triangulated surface equipped with non-interacting, up-down symmetry breaking inclusions where only terms coupled linearly to mean-curvature are maintained. We show that the lateral configurational entropy plays a key role for the mechanical properties of the semi-flexible membrane, e.g. a pronounced softening at intermediate inclusion coverages of the membrane and generation of membrane tension. Tensionless framed membranes will remain quasi-flat up to some threshold coverage, where a shape instability occurs with formation of pearling or tubular membranes, which below full coverage is associated with segregation of inclusions between the curved and flat membrane geometries. For inclusions with preference for highly curved membranes the instability appears at dilute inclusion coverages and is accompanied by strong configurational fluctuations.

Membrane-Mediated Cooperativity of Proteins

Besides direct protein-protein interactions, indirect interactions mediated by membranes play an important role for the assembly and cooperative function of proteins in membrane shaping and adhesion. The intricate shapes of biological membranes are generated by proteins that locally induce membrane curvature. Indirect curvature-mediated interactions between these proteins arise because the proteins jointly affect the bending energy of the membranes. These curvature-mediated interactions are attractive for crescent-shaped proteins and are a driving force in the assembly of the proteins during membrane tubulation. Membrane adhesion results from the binding of receptor and ligand proteins that are anchored in the apposing membranes. The binding of these proteins strongly depends on nanoscale shape fluctuations of the membranes, leading to a fluctuation-mediated binding cooperativity. A length mismatch between receptor-ligand complexes in membrane adhesion zones causes repulsive curvature-mediated interactions that are a driving force for the length-based segregation of proteins during membrane adhesion.

Mechanism of Shiga Toxin Clustering on Membranes

The bacterial Shiga toxin interacts with its cellular receptor, the glycosphingolipid globotriaosylceramide (Gb3 or CD77), as a first step to entering target cells. Previous studies have shown that toxin molecules cluster on the plasma membrane, despite the apparent lack of direct interactions between them. The precise mechanism by which this clustering occurs remains poorly defined. Here, we used vesicle and cell systems and computer simulations to show that line tension due to curvature, height, or compositional mismatch, and lipid or solvent depletion cannot drive the clustering of Shiga toxin molecules. By contrast, in coarse-grained computer simulations, a correlation was found between clustering and toxin nanoparticle-driven suppression of membrane fluctuations, and experimentally we observed that clustering required the toxin molecules to be tightly bound to the membrane surface. The most likely interpretation of these findings is that a membrane fluctuation-induced force generates an effective attraction between toxin molecules. Such force would be of similar strength to the electrostatic force at separations around 1 nm, remain strong at distances up to the size of toxin molecules (several nanometers), and persist even beyond. This force is predicted to operate between manufactured nanoparticles providing they are sufficiently rigid and tightly bound to the plasma membrane, thereby suggesting a route for the targeting of nanoparticles to cells for biomedical applications.

Budding and vesiculation induced by conical membrane inclusions

Conical inclusions in a lipid bilayer generate an overall spontaneous curvature of the membrane that depends on concentration and geometry of the inclusions. Examples are integral and attached membrane proteins, viruses, and lipid domains. We propose an analytical model to study budding and vesiculation of the lipid bilayer membrane, which is based on the membrane bending energy and the translational entropy of the inclusions. If the inclusions are placed on a membrane with similar curvature radius, their repulsive membrane-mediated interaction is screened. Therefore, for high inclusion density the inclusions aggregate, induce bud formation, and finally vesiculation. Already with the bending energy alone our model allows the prediction of bud radii. However, in case the inclusions induce a single large vesicle to split into two smaller vesicles, bending energy alone predicts that the smaller vesicles have different sizes whereas the translational entropy favors the formation of equal-sized vesicles. Our results agree well with those of recent computer simulations.

Cylindrical inclusions in a copolymer membrane

The membrane-mediated interaction between two parallel, cylindrical inclusions is investigated by using the self-consistent field theory (SCFT). The rodlike inclusions are located within the interior of the bilayer and are enveloped by two monolayers. They may exhibit one of the two basic types of behaviors involving pinching two monolayers together and swelling them outward. For different parameters, we calculate the density profile of the deformation membrane, the associated interaction free energy, as well as the conformational entropy of polymer chains. The similarity of the two types of interaction potentials is the qualitative characteristics. An energy barrier separates an attractive from a repulsive region; the repulsive region is preceded by a weak attraction at a large distance. The difference between them, which is due to the different contact environments around the rods, lies in the appearance of a small barrier at a short distance in the pinching structure. Particular emphasis is put on the closely energetic and entropic analyses of the interaction potential. We show that the chemical potential energy has provided a qualitative trend and roughly dominated the basic shape of the interaction potential; the amphiphile entropy in the swelling structure and the solvent entropy in the pinching structure, combined with the corresponding chemical potential energy, are responsible for the repulsive barrier at an intermediate distance and for the weak attraction at a large distance, respectively. The influence of inclusion hydrophobicity on the interaction potential is taken into account. In particular, the pinching and swelling structures can appear and can transform into each other in a system at intermediate hydrophobicity.

Receptor aggregation by intermembrane interactions: a Monte Carlo study

The lateral organization of receptors on cell surfaces is critically important to their function; many receptors transmit transmembrane signals when redistributed into clusters, while the response of others is potentiated by their aggregation. Cell-cell contact can play a crucial role in receptor aggregation, even when the bonds between receptors on one cell and ligands on the other are monovalent. Monte Carlo simulations on a two-membrane model were carried out to determine whether weak enthalpic interactions among receptors in one membrane, and among ligands in another, can work synergistically to give large-scale clustering when the two membranes are brought into contact. The simulations give support to such a clustering mechanism. In addition, because clustering is a cooperative process akin to a phase separation, individual receptors and ligands may undergo repeated binding and unbinding while in a clustered "phase," and a single ligand could interact with multiple different receptor partners. The results suggest a resolution of the dichotomy between serial triggering and aggregation models of T cell activation.

Dynamic simulations of membranes with cytoskeletal interactions

We describe a simulation algorithm for the dynamics of elastic membrane sheets over long length and time scales. Our model includes implicit hydrodynamic coupling between membrane and surrounding solvent and allows for arbitrary external forces acting on the membrane surface. In particular, the methodology is well suited to studying membranes in interaction with cytoskeletal filaments. We present results for the thermal undulations of a lipid bilayer attached to a regular network of spectrin filaments as a model for the red blood cell membrane. The dynamic fluctuations of the bilayer over the spectrin network are quantified and used to predict the macroscopic diffusion constant of band 3 on the surface of the red blood cell. We find that thermal undulations likely play a role in the mobility of band 3 in the plane of the erythrocyte membrane.

Effect of membrane characteristics on phase separation and domain formation in cholesterol-lipid mixtures

We examine, using an analytical mean-field model, the distribution of cholesterol in a lipid bilayer. The model accounts for the perturbation of lipid packing induced by the embedded cholesterol, in a manner similar to that of transmembrane proteins. We find that the membrane-induced interactions between embedded cholesterol molecules vary as a function of the cholesterol content. Thus, the effective lipid-cholesterol interaction is concentration-dependent. Moreover, it transitions from repulsive to attractive to repulsive as the cholesterol content increases. As the concentration of cholesterol in the bilayer exceeds a critical value, phase separation occurs. The coexistence between cholesterol-rich and cholesterol-poor domains is universal for any bilayer parameters, although the composition of the cholesterol-rich phase varies as a function of the lipid properties. Although we do not assume specific cholesterol-lipid interactions or the formation of a lipid-cholesterol cluster, we find that the composition of the cholesterol-rich domains is constant, independent of the cholesterol content in the bilayer.

Dynamics of pinned membranes with application to protein diffusion on the surface of red blood cells

We present a theoretical treatment and simulation algorithm for the dynamics of Helfrich elastic membrane surfaces in the presence of general harmonic perturbations and hydrodynamic coupling to the surrounding solvent. In the limit of localized and strong interactions, this harmonic model can be used to pin the membrane to intracellular/intercellular structures. We consider the case of pinning to the cytoskeleton and use such a model to estimate the macroscopic diffusion constant for band 3 protein on the surface of human erythrocytes. Comparison to experimental results suggests that thermal undulations of the membrane surface should play a significant role in protein mobility on the red blood cell.