Membrane shape remodeling by protein crowding

Stability of a biomembrane tube covered with proteins

Membrane tubes are essential structural features in cells that facilitate biomaterial transport and inter- and intracellular signaling. The shape of these tubes can be regulated by the proteins that surround and adhere to them. We study the stability of a biomembrane tube coated with proteins by combining linear stability analysis, out-of-equilibrium hydrodynamic calculations, and numerical solutions of a Helfrich-like membrane model. Our analysis demonstrates that both long- and short-wavelength perturbations can destabilize the tubes. Numerical simulations confirm the derived linear stability criteria and yield the nonlinearly perturbed vesicle shapes. Our study highlights the interplay between membrane shape and protein density, where the shape instability concurs with a redistribution of proteins into a banded pattern.

Sequence Sensitivity in Membrane Remodeling by Polyampholyte Condensates

Intrinsically disordered peptides (IDPs) have been found to undergo liquid-liquid phase separation (LLPS) and produce complex coacervates that play numerous regulatory roles in the cell. Recent experimental studies have discovered that LLPS at or near the membrane surface helps in the biomolecular organization during signaling events and can significantly alter the membrane morphology. However, the molecular mechanism and microscopic details of such processes still remain unclear. Here we study the effect of polyampholyte and polyelectrolyte condensation on two different anionic membranes, as they represent a majority of naturally occurring IDPs. The polyampholytes are fifty-residue polymers, made of glutamate(E) and lysine(K) with different charge patterns. The polyelectrolytes are separate chains of E25 and K25. We first calibrate the MARTINI v3.0 force field and then perform long-time-scale coarse-grained molecular dynamics simulations. We find that condensates formed by all the polyampholytes get adsorbed on the membrane. However, the strong polyampholytes (i.e., blocky sequences) can remodel the membranes more prominently than the weaker ones (i.e., scrambled sequences). Condensates formed by the blocky sequences induce a significant negative curvature (∼0.1 nm-1) and local demixing of lipids, whereas those by the scrambled sequences tend to wet the membrane to a greater extent without generating significant curvature or demixing. We perform several microscopic analyses to characterize the nature of the interaction between membranes and these condensates. Our analyses of interaction energetics reveal that membrane remodeling and/or wetting are favored by enhanced interactions between polyampholytes with lipids and the counterions.

Membranes get in shape: Biophysics of curving bilayers

Membrane curvature is ubiquitous and essential in cell biology. Curved membranes have several distinct features, including specific protein and lipid sorting, distinct lipid ordering, and changes in transbilayer stress. Curvature also interplays with membrane tension to generate forces that change membrane shape. This research highlight summarizes recent contributions to this topic published in Biophysical Journal.

Defining permeability of curved membranes in molecular dynamics simulations

Many phospholipid membranes in the cell have a high curvature; for instance, in caveolae, mitochondrial crystae, nanotubes, membrane pearls, small liposomes, or exosomes. Molecular dynamics (MD) simulations are a computational tool to gain insight in the transport behavior at the atomic scale. Membrane permeability is a key kinetic property that might be affected in these highly curved membranes. Unfortunately, the geometry of highly curved membranes creates ambiguity in the permeability value, even with an arbitrarily large factor purely based on geometry, caused by the radial flux not being a constant value in steady state. In this contribution, the ambiguity in permeability for liposomes is countered by providing a new permeability definition. First, the inhomogeneous solubility diffusion model based on the Smoluchowski equation is solved analytically under radial symmetry, from which the entrance and escape permeabilities are defined. Next, the liposome permeability is defined guided by the criterion that a flat and curved membrane should have equal permeability, in case these were to be carved out from an imaginary homogeneous medium. With this criterion, our new definition allows for a fair comparison of flat and curved membranes. The definition is then transferred to the counting method, which is a practical computational approach to derive permeability by counting complete membrane crossings. Finally, the usability of the approach is illustrated with MD simulations of diphosphatidylcholine (DPPC) bilayers, without or with some cholesterol content. Our new liposome permeability definition allows us to compare a spherically shaped membrane with its flat counterpart, thus showcasing how the curvature effect on membrane transport may be assessed.

Structural Basis of the Interaction of the G Proteins, Gαi1, Gβ1γ2 and Gαi1β1γ2, with Membrane Microdomains and Their Relationship to Cell Localization and Activity

GPCRs receive signals from diverse messengers and activate G proteins that regulate downstream signaling effectors. Efficient signaling is achieved through the organization of these proteins in membranes. Thus, protein-lipid interactions play a critical role in bringing G proteins together in specific membrane microdomains with signaling partners. Significantly, the molecular basis underlying the membrane distribution of each G protein isoform, fundamental to fully understanding subsequent cell signaling, remains largely unclear. We used model membranes with lipid composition resembling different membrane microdomains, and monomeric, dimeric and trimeric Gi proteins with or without single and multiple mutations to investigate the structural bases of G protein-membrane interactions. We demonstrated that cationic amino acids in the N-terminal region of the Gαi₁ and C-terminal region of the Gγ₂ subunit, as well as their myristoyl, palmitoyl and geranylgeranyl moieties, define the differential G protein form interactions with membranes containing different lipid classes (PC, PS, PE, SM, Cho) and the various microdomains they may form (Lo, Ld, PC bilayer, charged, etc.). These new findings in part explain the molecular basis underlying amphitropic protein translocation to membranes and localization to different membrane microdomains and the role of these interactions in cell signal propagation, pathophysiology and therapies targeted to lipid membranes.

Strength in numbers: effect of protein crowding on the shape of cell membranes

Continuous reshaping of the plasma membrane into pleomorphic shapes is critical for a plethora of cellular functions. How the cell carries out this enigmatic control of membrane remodeling has remained an active research field for decades and several molecular and biophysical mechanisms have shown to be involved in overcoming the energy barrier associated with membrane bending. The reported mechanisms behind membrane bending have been largely concerned with structural protein features, however, in the last decade, reports on the ability of densely packed proteins to bend membranes by protein-protein crowding, have challenged prevailing mechanistic views. Crowding has now been shown to generate spontaneous vesicle formation and tubular morphologies on cell- and model membranes, demonstrating crowding as a relevant player involved in the bending of membranes. Still, current research is largely based on unnatural overexpression of proteins in non-native domains, and together with efforts in modeling, this has led to questioning the in vivo impact of crowding. In this review, we examine this previously overlooked mechanism by summarizing recent advances in the understanding of protein-protein crowding and its prevalence in cellular membrane-shaping processes.

Crowding-induced membrane remodeling: Interplay of membrane tension, polymer density, architecture

The plasma membrane hosts a wide range of biomolecules, mainly proteins and carbohydrates, that mediate cellular interactions with its environment. The crowding of such biomolecules regulates cellular morphologies and cellular trafficking. Recent discoveries have shown that the structure and density of cell surface polymers and hence the signaling machinery change with the state of the cell, especially in cancer progression. The alterations in membrane-attached glycocalyx and glycosylation of proteins and lipids are common features of cancer cells. The overexpression of glycocalyx polymers, such as mucin and hyaluronan, strongly correlates with cancer metastasis. Here, we present a mesoscale biophysics-based model that accounts for the shape regulation of membranes by crowding of membrane-attached biopolymer-glycocalyx and actin networks. Our computational model is based on the dynamically triangulated Monte Carlo model for membranes and coarse-grained representations of polymer chains. The model allows us to investigate the crowding-induced shape transformations in cell membranes in a tension- and graft polymer density-dependent manner. Our results show that the number of membrane protrusions and their shape depend on membrane tension, with higher membrane tension inducing more tubular protrusions than the vesicular shapes formed at low tension at high surface coverage of polymers. The shape transformations occur above the threshold density predicted by the polymer brush theory, but this threshold also depends on the membrane tension. Increasing the size of the polymer, either by changing the length or by adding side chains, is shown to increase the crowding-induced curvature. The effect of crowding is more prominent for flexible polymers than for semiflexible rigid polymers. We also present an extension of the model that incorporates properties of the actin-like filament networks and demonstrate how tubular structures can be generated by biopolymer crowding on the cytosolic side of cell membranes.

Mem3DG: Modeling membrane mechanochemical dynamics in 3D using discrete differential geometry

Biomembranes adopt varying morphologies that are vital to cellular functions. Many studies use computational modeling to understand how various mechanochemical factors contribute to membrane shape transformations. Compared with approximation-based methods (e.g., finite element method [FEM]), the class of discrete mesh models offers greater flexibility to simulate complex physics and shapes in three dimensions; its formulation produces an efficient algorithm while maintaining coordinate-free geometric descriptions. However, ambiguities in geometric definitions in the discrete context have led to a lack of consensus on which discrete mesh model is theoretically and numerically optimal; a bijective relationship between the terms contributing to both the energy and forces from the discrete and smooth geometric theories remains to be established. We address this and present an extensible framework, Mem3DG, for modeling 3D mechanochemical dynamics of membranes based on discrete differential geometry (DDG) on triangulated meshes. The formalism of DDG resolves the inconsistency and provides a unifying perspective on how to relate the smooth and discrete energy and forces. To demonstrate, Mem3DG is used to model a sequence of examples with increasing mechanochemical complexity: recovering classical shape transformations such as 1) biconcave disk, dumbbell, and unduloid; and 2) spherical bud on spherical, flat-patch membrane; investigating how the coupling of membrane mechanics with protein mobility jointly affects phase and shape transformation. As high-resolution 3D imaging of membrane ultrastructure becomes more readily available, we envision Mem3DG to be applied as an end-to-end tool to simulate realistic cell geometry under user-specified mechanochemical conditions.

Coacervation of poly-electrolytes in the presence of lipid bilayers: mutual alteration of structure and morphology

Many intrinsically disordered peptides have been shown to undergo liquid-liquid phase separation and form complex coacervates, which play various regulatory roles in the cell. Recent experimental studies found that such phase separation processes may also occur at the lipid membrane surface and help organize biomolecules during signaling events; in some cases, phase separation of proteins at the membrane surface was also observed to lead to significant remodeling of the membrane morphology. The molecular mechanisms that govern the interactions between complex coacervates and lipid membranes and the impacts of such interactions on their structure and morphology, however, remain unclear. Here we study the coacervation of poly-glutamate (E30) and poly-lysine (K30) in the presence of lipid bilayers of different compositions. We carry out explicit-solvent coarse-grained molecular dynamics simulations by using the MARTINI (v3.0) force-field. We find that more than 20% anionic lipids are required for the coacervate to form stable contact with the bilayer. Upon wetting, the coacervate induces negative curvature to the bilayer and facilitates local lipid demixing, without any peptide insertion. The magnitude of negative curvature, extent of lipid demixing, and asphericity of the coacervate increase with the concentration of anionic lipids. Overall, we observe a decrease in the number of contacts among the polyelectrolytes as the droplet spreads over the bilayer. Therefore, unlike previous suggestions, interactions among polyelectrolytes do not constitute a driving force for the membrane bending upon wetting by the coacervate. Rather, analysis of interaction energy components suggests that bending of the membrane is favored by enhanced interactions between polyelectrolytes with lipids as well as with counterions. Kinetic studies reveal that, at the studied polyelectrolyte concentrations, the coacervate formation precedes bilayer wetting.

On the Role of Electrostatic Repulsion in Topological Defect-Driven Membrane Fission

Within a modified Langevin Poisson–Boltzmann model of electric double layers, we derived an analytical expression for osmotic pressure between two charged surfaces. The orientational ordering of the water dipoles as well as the space dependencies of electric potentials, electric fields, and osmotic pressure between two charged spheres were taken into account in the model. Thus, we were able to capture the interaction between the parent cell and connected daughter vesicle or the interactions between neighbouring beads in necklace-like membrane protrusions. The predicted repulsion between them can facilitate the topological antidefect-driven fission of membrane daughter vesicles and the fission of beads of undulated membrane protrusions.