Membrane tubule formation by banana-shaped proteins with or without transient network structure

Spatiotemporal pattern formation of membranes induced by surface molecular binding/unbinding

Nonequilibrium membrane pattern formation is studied using meshless membrane simulation. We consider that molecules bind to either surface of a bilayer membrane and move to the opposite leaflet by flip-flop. When binding does not modify the membrane properties and the transfer rates among the three states are cyclically symmetric, the membrane exhibits spiral-wave and homogeneous-cycling modes at high and low binding rates, respectively, as in an off-lattice cyclic Potts model. When binding changes the membrane spontaneous curvature, these spatiotemporal dynamics are coupled with microphase separation. When two symmetric membrane surfaces are in thermal equilibrium, the membrane domains form 4.8.8 tiling patterns in addition to stripe and spot patterns. In nonequilibrium conditions, moving biphasic domains and time-irreversible fluctuating patterns appear. The domains move ballistically or diffusively depending on the conditions.

Particle Deformability Enables Control of Interactions between Membrane-Anchored Nanoparticles

Nanoparticles adsorbed on a membrane can induce deformations of the membrane that give rise to effective interactions between the particles. Previous studies have focused primarily on rigid nanoparticles with fixed shapes. However, DNA origami technology has enabled the creation of deformable nanostructures with controllable shapes and mechanical properties, presenting new opportunities to modulate interactions between particles adsorbed on deformable surfaces. Here we use coarse-grained molecular dynamics simulations to investigate deformable, hinge-like nanostructures anchored to lipid membranes via cholesterol anchors. We characterize deformations of the particles and membrane as a function of the hinge stiffness. Flexible particles adopt open configurations to conform to a flat membrane, whereas stiffer particles induce deformations of the membrane. We further show that particles spontaneously aggregate and that cooperative effects lead to changes in their shape when they are close together. Using umbrella sampling methods, we quantify the effective interaction between two particles and show that stiffer hinge-like particles experience stronger and longer-ranged attraction. Our results demonstrate that interactions between deformable, membrane-anchored nanoparticles can be controlled by modifying mechanical properties of the particles, suggesting new ways to modulate the self-assembly of particles on deformable surfaces.

Curvature sensing of curvature-inducing proteins with internal structure

Many types of peripheral and transmembrane proteins can sense and generate membrane curvature. Laterally isotropic proteins and crescent proteins with twofold rotational symmetry, such as Bin/Amphiphysin/Rvs superfamily proteins, have been studied theoretically. However, proteins often have an asymmetric structure or a higher rotational symmetry. We studied theoretically the curvature sensing of proteins with asymmetric structures and structural deformations. First, we examined proteins consisting of two rodlike segments. When proteins have mirror symmetry, their sensing ability is similar to that of single-rod proteins; hence, with increasing protein density on a cylindrical membrane tube, a second- or first-order transition occurs at a middle or small tube radius, respectively. As asymmetry is introduced, this transition becomes a continuous change and metastable states appear at high protein densities. Protein with threefold, fivefold, or higher rotational symmetry has laterally isotropic bending energy. However, when a structural deformation is allowed, the protein can have a preferred orientation and stronger curvature sensing.

Membrane domain formation induced by binding/unbinding of curvature-inducing molecules on both membrane surfaces

The domain formation of curvature-inducing molecules, such as peripheral or transmembrane proteins and conical surfactants, is studied in thermal equilibrium and nonequilibrium steady states using meshless membrane simulations. These molecules can bind to both surfaces of a bilayer membrane and also move to the opposite leaflet by a flip-flop. Under symmetric conditions for the two leaflets, the membrane domains form checkerboard patterns in addition to striped and spot patterns. The unbound membrane stabilizes the vertices of the checkerboard. Under asymmetric conditions, the domains form kagome-lattice and thread-like patterns. In the nonequilibrium steady states, a flow of the binding molecules between the upper and lower solutions can occur via flip-flop.

Modulation of self-organizing circuits at deforming membranes by intracellular and extracellular factors

Mechanical forces exerted to the plasma membrane induce cell shape changes. These transient shape changes trigger, among others, enrichment of curvature-sensitive molecules at deforming membrane sites. Strikingly, some curvature-sensing molecules not only detect membrane deformation but can also alter the amplitude of forces that caused to shape changes in the first place. This dual ability of sensing and inducing membrane deformation leads to the formation of curvature-dependent self-organizing signaling circuits. How these cell-autonomous circuits are affected by auxiliary parameters from inside and outside of the cell has remained largely elusive. Here, we explore how such factors modulate self-organization at the micro-scale and its emerging properties at the macroscale.

Membrane curvature as a signal to ensure robustness of diverse cellular processes

An increasing corpus of research has demonstrated that membrane shape, generated either by the external environment of the cell or by intrinsic mechanisms such as cytokinesis and vesicle or organelle formation, is an important parameter in the control of diverse cellular processes. In this review we discuss recent findings that demonstrate how membrane curvature (from nanometer to micron length-scales) alters protein function. We describe an expanding toolkit for experimentally modulating membrane curvature to reveal effects on protein function, and discuss how membrane curvature - far from being a passive consequence of the physical environment and the internal protein activity of a cell - is an important signal that controls protein affinity and enzymatic activity to ensure robust forward progression of key processes within the cell.

Membrane shape deformation induced by curvature-inducing proteins consisting of chiral crescent binding and intrinsically disordered domains

Curvature-inducing proteins containing a bin/amphiphysin/Rvs domain often have intrinsically disordered domains. Recent experiments have shown that these disordered chains enhance curvature sensing and generation. Here, we report on the modification of protein–membrane interactions by disordered chains using meshless membrane simulations. The protein and bound membrane are modeled together as a chiral crescent protein rod with two excluded-volume chains. As the chain length increases, the repulsion between them reduces the cluster size of the proteins. It induces spindle-shaped vesicles and a transition between arc-shaped and circular protein assemblies in a disk-shaped vesicle. For flat membranes, an intermediate chain length induces many tubules owing to the repulsion between the protein assemblies, whereas longer chains promote perpendicular elongation of tubules. Moreover, protein rods with zero rod curvature and sufficiently long chains stabilize the spherical buds. For proteins with a negative rod curvature, an intermediate chain length induces a rugged membrane with branched protein assemblies, whereas longer chains induce the formation of tubules with periodic concave-ring structures.

Binding of anisotropic curvature-inducing proteins onto membrane tubes

Bin/Amphiphysin/Rvs superfamily proteins and other curvature-inducing proteins have anisotropic shapes and anisotropically bend biomembranes. Here, we report how the anisotropic proteins bind the membrane tube and are orientationally ordered using mean-field theory including an orientation-dependent excluded volume. The proteins exhibit a second-order or first-order nematic transition with increasing protein density depending on the radius of the membrane tube. The tube curvatures for the maximum protein binding and orientational order are different and varied by the protein density and rigidity. As the external force along the tube axis increases, a first-order transition from a large tube radius with low protein density to a small radius with high density occurs once, and subsequently, the protein orientation tilts to the tube-axis direction. When an isotropic bending energy is used for the proteins with an elliptic shape, the force-dependence curves become symmetric and the first-order transition occurs twice. This theory quantitatively reproduces the results of meshless membrane simulation for short proteins, whereas deviations are seen for long proteins owing to the formation of protein clusters.

Membrane-Mediated Interactions Between Protein Inclusions

Integral or peripheral membrane proteins, or protein oligomers often get close to each other on cell membranes and carry out biological tasks in a collective manner. In addition to electrostatic and van der Waals interactions, those proteins also experience membrane-mediated interactions, which may be necessary for their functionality. The membrane-mediated interactions originate from perturbation of lipid membranes by the presence of protein inclusions, and have been the subject of intensive research in membrane biophysics. Here we review both theoretical and numerical studies of such interactions for membrane proteins and for nanoparticles bound to lipid membranes.

Binding of curvature-inducing proteins onto tethered vesicles

A tethered vesicle, which consists of a cylindrical membrane tube and a spherical vesicle, is produced by a mechanical force that is experimentally imposed by optical tweezers and a micropipette. This tethered vesicle is employed for examining the curvature sensing of curvature-inducing proteins. In this study, we clarify how the binding of proteins with a laterally isotropic spontaneous curvature senses and generates the membrane curvatures of the tethered vesicle using mean-field theory and meshless membrane simulation. The force-dependence curves of the protein density in the membrane tube and the tube curvature are reflection symmetric and point symmetric, respectively, from the force point, in which the tube has a sensing curvature. The bending rigidity and spontaneous curvature of the bound proteins can be estimated from these force-dependence curves. First-order transitions can occur between low and high protein densities in the tube at both low and high force amplitudes. The simulation results of the homogeneous phases agree very well with the theoretical predictions. In addition, beaded-necklace-like tubes with microphase separation are found in the simulation.

Dynamic mechanochemical feedback between curved membranes and BAR protein self-organization

In many physiological situations, BAR proteins reshape membranes with pre-existing curvature (templates), contributing to essential cellular processes. However, the mechanism and the biological implications of this reshaping process remain unclear. Here we show, both experimentally and through modelling, that BAR proteins reshape low curvature membrane templates through a mechanochemical phase transition. This phenomenon depends on initial template shape and involves the co-existence and progressive transition between distinct local states in terms of molecular organization (protein arrangement and density) and membrane shape (template size and spherical versus cylindrical curvature). Further, we demonstrate in cells that this phenomenon enables a mechanotransduction mode, in which cellular stretch leads to the mechanical formation of membrane templates, which are then reshaped into tubules by BAR proteins. Our results demonstrate the interplay between membrane mechanics and BAR protein molecular organization, integrating curvature sensing and generation in a comprehensive framework with implications for cell mechanical responses.

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.

Three-dimensional morphodynamic simulations of macropinocytic cups

Macropinocytosis refers to the non-specific uptake of extracellular fluid, which plays ubiquitous roles in cell growth, immune surveillance, and virus entry. Despite its widespread occurrence, it remains unclear how its initial cup-shaped plasma membrane extensions form without any external solid support, as opposed to the process of particle uptake during phagocytosis. Here, by developing a computational framework that describes the coupling between the bistable reaction-diffusion processes of active signaling patches and membrane deformation, we demonstrated that the protrusive force localized to the edge of the patches can give rise to a self-enclosing cup structure, without further assumptions of local bending or contraction. Efficient uptake requires a balance among the patch size, magnitude of protrusive force, and cortical tension. Furthermore, our model exhibits cyclic cup formation, coexistence of multiple cups, and cup-splitting, indicating that these complex morphologies self-organize via a common mutually-dependent process of reaction-diffusion and membrane deformation.

A bacterial membrane sculpting protein with BAR domain-like activity

Bin/Amphiphysin/RVS (BAR) domain proteins belong to a superfamily of coiled-coil proteins influencing membrane curvature in eukaryotes and are associated with vesicle biogenesis, vesicle-mediated protein trafficking, and intracellular signaling. Here, we report a bacterial protein with BAR domain-like activity, BdpA, from Shewanella oneidensis MR-1, known to produce redox-active membrane vesicles and micrometer-scale outer membrane extensions (OMEs). BdpA is required for uniform size distribution of membrane vesicles and influences scaffolding of OMEs into a consistent diameter and curvature. Cryo-TEM reveals that a strain lacking BdpA produces lobed, disordered OMEs rather than membrane tubules or narrow chains produced by the wild-type strain. Overexpression of BdpA promotes OME formation during planktonic growth of S. oneidensis where they are not typically observed. Heterologous expression results in OME production in Marinobacter atlanticus and Escherichia coli. Based on the ability of BdpA to alter membrane architecture in vivo, we propose that BdpA and its homologs comprise a newly identified class of bacterial BAR domain-like proteins.

Vesicle budding induced by binding of curvature-inducing proteins

Vesicle budding induced by protein binding that generates an isotropic spontaneous curvature is studied using a mean-field theory. Many spherical buds are formed via protein binding. As the binding chemical potential increases, the proteins first bind to the buds and then to the remainder of the vesicle. For a high spontaneous curvature and/or high bending rigidity of the bound membrane, it is found that a first-order transition occurs between a small number of large buds and a large number of small buds. These two states coexist around the transition point. The proposed scheme is simple and easily applicable to many interaction types, so we investigate the effects of interprotein interactions, the protein-insertion-induced changes in area, the variation of the saddle-splay modulus, and the area-difference-elasticity energy. The differences in the preferred curvatures for curvature sensing and generation are also clarified.

Recent developments in membrane curvature sensing and induction by proteins

BACKGROUND

Membrane-bound intracellular organelles have characteristic shapes attributed to different local membrane curvatures, and these attributes are conserved across species. Over the past decade, it has been confirmed that specific proteins control the large curvatures of the membrane, whereas many others due to their specific structural features can sense the curvatures and bind to the specific geometrical cues. Elucidating the interplay between sensing and induction is indispensable to understand the mechanisms behind various biological processes such as vesicular trafficking and budding.

SCOPE OF REVIEW

We provide an overview of major classes of membrane proteins and the mechanisms of curvature sensing and induction. We then discuss the importance of membrane elastic characteristics to induce the membrane shapes similar to intracellular organelles. Finally, we survey recently available assays developed for studying the curvature sensing and induction by many proteins.

MAJOR CONCLUSIONS

Recent theoretical/computational modeling along with experimental studies have uncovered fascinating connections between lipid membrane and protein interactions. However, the phenomena of protein localization and synchronization to generate spatiotemporal dynamics in membrane morphology are yet to be fully understood.

GENERAL SIGNIFICANCE

The understanding of protein-membrane interactions is essential to shed light on various biological processes. This further enables the technological applications of many natural proteins/peptides in therapeutic treatments. The studies of membrane dynamic shapes help to understand the fundamental functions of membranes, while the medicinal roles of various macromolecules (such as proteins, peptides, etc.) are being increasingly investigated.

A theory of ordering of elongated and curved proteins on membranes driven by density and curvature

Cell membranes interact with a myriad of curvature-active proteins that control membrane morphology and are responsible for mechanosensation and mechanotransduction. Some of these proteins, such as those containing BAR domains, are curved and elongated, and hence may adopt different states of orientational order, from isotropic to maximize entropy to nematic as a result of crowding or to adapt to the curvature of the underlying membrane. Here, extending the classical work of Onsager for ordering in hard particle systems and that of [E. S. Nascimento et al., Phys. Rev. E, 2017, 96, 022704], we develop a mean-field density functional theory to predict the orientational order and evaluate the free energy of ensembles of elongated and curved objects on curved membranes. This theory depends on the microscopic properties of the particles and explains how a density-dependent isotropic-to-nematic transition is modified by anisotropic curvature. We also examine the coexistence of isotropic and nematic phases. This theory predicts how ordering depends on geometry but we assume here that the geometry is fixed. It also lays the ground to understand the interplay between membrane reshaping by BAR proteins and molecular order, examined by [Le Roux et al., submitted, 2020].

On the Role of Curved Membrane Nanodomains, and Passive and Active Skeleton Forces in the Determination of Cell Shape and Membrane Budding

Biological membranes are composed of isotropic and anisotropic curved nanodomains. Anisotropic membrane components, such as Bin/Amphiphysin/Rvs (BAR) superfamily protein domains, could trigger/facilitate the growth of membrane tubular protrusions, while isotropic curved nanodomains may induce undulated (necklace-like) membrane protrusions. We review the role of isotropic and anisotropic membrane nanodomains in stability of tubular and undulated membrane structures generated or stabilized by cyto- or membrane-skeleton. We also describe the theory of spontaneous self-assembly of isotropic curved membrane nanodomains and derive the critical concentration above which the spontaneous necklace-like membrane protrusion growth is favorable. We show that the actin cytoskeleton growth inside the vesicle or cell can change its equilibrium shape, induce higher degree of segregation of membrane nanodomains or even alter the average orientation angle of anisotropic nanodomains such as BAR domains. These effects may indicate whether the actin cytoskeleton role is only to stabilize membrane protrusions or to generate them by stretching the vesicle membrane. Furthermore, we demonstrate that by taking into account the in-plane orientational ordering of anisotropic membrane nanodomains, direct interactions between them and the extrinsic (deviatoric) curvature elasticity, it is possible to explain the experimentally observed stability of oblate (discocyte) shapes of red blood cells in a broad interval of cell reduced volume. Finally, we present results of numerical calculations and Monte-Carlo simulations which indicate that the active forces of membrane skeleton and cytoskeleton applied to plasma membrane may considerably influence cell shape and membrane budding.

Binding, unbinding and aggregation of crescent-shaped nanoparticles on nanoscale tubular membranes

Using molecular dynamics simulations of a coarse-grained implicit solvent model, we investigate the binding of crescent-shaped nanoparticles (NPs) on tubular lipid membranes. The NPs adhere to the membrane through their concave side. We found that the binding/unbinding transition is first-order, with the threshold binding energy being higher than the unbinding threshold, and the energy barrier between the bound and unbound states at the transition that increases with increasing the NP's arclength L_np or curvature mismatch μ = R_c/R_np, where R_c and R_np are the radii of curvature of the tubular membrane and the NP, respectively. Furthermore, we found that the threshold binding energy increases with increasing either L_np or μ. NPs with curvature larger than that of the tubule (μ > 1) lie perpendicularly to the tubule's axis. However, for μ smaller than a specific arclength-dependent mismatch μ*, the NPs are tilted with respect to the tubule's axis, with the tilt angle that increases with decreasing μ. We also investigated the self-assembly of the NPs on the tubule at relatively weak adhesion strength and found that for μ > 1 and high values of L_np, the NPs self-assemble into linear chains, and lie side-by-side. For μ < μ* and high L_np, the NPs also self-assemble into chains, while being tilted with respect to the tubule's axis.

Pattern formation in reaction-diffusion system on membrane with mechanochemical feedback

Shapes of biological membranes are dynamically regulated in living cells. Although membrane shape deformation by proteins at thermal equilibrium has been extensively studied, nonequilibrium dynamics have been much less explored. Recently, chemical reaction propagation has been experimentally observed in plasma membranes. Thus, it is important to understand how the reaction-diffusion dynamics are modified on deformable curved membranes. Here, we investigated nonequilibrium pattern formation on vesicles induced by mechanochemical feedback between membrane deformation and chemical reactions, using dynamically triangulated membrane simulations combined with the Brusselator model. We found that membrane deformation changes stable patterns relative to those that occur on a non-deformable curved surface, as determined by linear stability analysis. We further found that budding and multi-spindle shapes are induced by Turing patterns, and we also observed the transition from oscillation patterns to stable spot patterns. Our results demonstrate the importance of mechanochemical feedback in pattern formation on deforming membranes.

Comparative Study of Curvature Sensing Mediated by F-BAR and an Intrinsically Disordered Region of FBP17

Membrane curvature has emerged as an intriguing physical principle underlying biological signaling and membrane trafficking. The CIP4/FBP17/Toca-1 F-BAR subfamily is unique in the BAR family because its structurally folded F-BAR domain does not contain any hydrophobic motifs that insert into membrane. Although widely assumed so, whether the banana-shaped F-BAR domain alone can sense curvature has never been experimentally demonstrated. Using a nanobar-supported lipid bilayer system, we found that the F-BAR domain of FBP17 displayed minimal curvature sensing in vitro. In comparison, an alternatively spliced intrinsically disordered region (IDR) adjacent to the F-BAR domain has the membrane curvature-sensing ability greatly exceeding that of F-BAR domain alone. In living cells, the presence of the IDR delayed the recruitment of FBP17 in curvature-coupled cortical waves. Collectively, we propose that contrary to the common belief, FBP17's curvature-sensing capability largely originates from IDR, and not the F-BAR domain alone.

Detachment of a fluid membrane from a substrate and vesiculation

The detachment dynamics of a fluid membrane with isotropic spontaneous curvature from a flat substrate are studied by using meshless membrane simulations. The membrane is detached from an open edge leading to vesicle formation. With strong adhesion, the competition between the bending and adhesion energies determines the minimum value of the spontaneous curvature for the detachment. In contrast, with weak adhesion, detachment occurs at smaller spontaneous curvatures due to the membrane thermal undulation. When parts of the membrane are pinned on the substrate, the detachment becomes slower and a remaining membrane patch forms straight or concave membrane edges. The edge undulation induces vesiculation of long strips and disk-shaped patches. Therefore, membrane rolling is obtained only for membrane strips shorter than the wavelength for deformation into unduloids. This suggests that the rolling observed for Ca²⁺-dependent membrane-binding proteins annexins A3, A4, A5, and A13 results from the anisotropic spontaneous curvature induced by the proteins.

Cellular Blebs and Membrane Invaginations Are Coupled through Membrane Tension Buffering

Bleb-type cellular protrusions play key roles in a range of biological processes. It was recently found that bleb growth is facilitated by a local supply of membrane from tubular invaginations, but the interplay between the expanding bleb and the membrane tubes remains poorly understood. On the one hand, the membrane area stored in tubes may serve as a reservoir for bleb expansion. On the other hand, the sequestering of excess membrane in stabilized invaginations may effectively increase the cell membrane tension, which suppresses spontaneous protrusions. Here, we investigate this duality through physical modeling and in vivo experiments. In agreement with observations, our model describes the transition into a tube-flattening mode of bleb expansion while also predicting that the blebbing rate is impaired by elevating the concentration of the curved membrane proteins that form the tubes. We show both theoretically and experimentally that the stabilizing effect of tubes could be counterbalanced by the cortical myosin contractility. Our results largely suggest that proteins able to induce membrane tubulation, such as those containing N-BAR domains, can buffer the effective membrane tension-a master regulator of all cell deformations.

Curving Cells Inside and Out: Roles of BAR Domain Proteins in Membrane Shaping and Its Cellular Implications

Many cellular processes rely on precise and timely deformation of the cell membrane. While many proteins participate in membrane reshaping and scission, usually in highly specialized ways, Bin/amphiphysin/Rvs (BAR) domain proteins play a pervasive role, as they not only participate in many aspects of cell trafficking but also are highly versatile membrane remodelers. Subtle changes in the shape and size of the BAR domain can greatly impact the way in which BAR domain proteins interact with the membrane. Furthermore, the activity of BAR domain proteins can be tuned by external physical parameters, and so they behave differently depending on protein surface density, membrane tension, or membrane shape. These proteins can form 3D structures that mold the membrane and alter its liquid properties, even promoting scission under various circumstances.As such, BAR domain proteins have numerous roles within the cell. Endocytosis is among the most highly studied processes in which BAR domain proteins take on important roles. Over the years, a more complete picture has emerged in which BAR domain proteins are tied to almost all intracellular compartments; examples include endosomal sorting and tubular networks in the endoplasmic reticulum and T-tubules. These proteins also have a role in autophagy, and their activity has been linked with cancer. Here, we briefly review the history of BAR domain protein discovery, discuss the mechanisms by which BAR domain proteins induce curvature, and attempt to settle important controversies in the field. Finally, we review BAR domain proteins in the context of a cell, highlighting their emerging roles in cell signaling and organelle shaping.

Shape transition from elliptical to cylindrical membrane tubes induced by chiral crescent-shaped protein rods

Proteins often form chiral assembly structures on a biomembrane. However, the role of the chirality in the interaction with an achiral membrane is poorly understood. Here, we report how chirality of crescent-shaped protein rods changes their assembly and tubulation using meshless membrane simulations. The achiral rods deformed the membrane tube into an elliptical shape by stabilizing the edges of the ellipse. In contrast, the chiral rods formed a helical assembly that generated a cylindrical membrane tube with a constant radius in addition to the elliptical tube. This helical assembly could be further stabilized by the direct side-to-side attraction between the protein rods. The chirality also promotes the tubulation from a flat membrane. These results agree with experimental findings of the constant radius of membrane tubules induced by the Bin/Amphiphysin/Rvs (BAR) superfamily proteins.

Theoretical study of vesicle shapes driven by coupling curved proteins and active cytoskeletal forces

Eukaryote cells have a flexible shape, which dynamically changes according to the function performed by the cell. One mechanism for deforming the cell membrane into the desired shape is through the expression of curved membrane proteins. Furthermore, these curved membrane proteins are often associated with the recruitment of the cytoskeleton, which then applies active forces that deform the membrane. This coupling between curvature and activity was previously explored theoretically in the linear limit of small deformations, and low dimensionality. Here we explore the unrestricted shapes of vesicles that contain active curved membrane proteins, in three-dimensions, using Monte-Carlo numerical simulations. The activity of the proteins is in the form of protrusive forces that push the membrane outwards, as may arise from the cytoskeleton of the cell due to actin or microtubule polymerization occurring near the membrane. For proteins that have an isotropic convex shape, the additional protrusive force enhances their tendency to aggregate and form membrane protrusions (buds). In addition, we find another transition from deformed spheres with necklace type aggregates, to flat pancake-shaped vesicles, where the curved proteins line the outer rim. This second transition is driven by the active forces, coupled to the spontaneous curvature, and the resulting configurations may shed light on the formation of sheet-like protrusions and lamellipodia of adhered and motile cells.

ACAP1 assembles into an unusual protein lattice for membrane deformation through multiple stages

Studies on the Bin-Amphiphysin-Rvs (BAR) domain have advanced a fundamental understanding of how proteins deform membrane. We previously showed that a BAR domain in tandem with a Pleckstrin Homology (PH domain) underlies the assembly of ACAP1 (Arfgap with Coil-coil, Ankryin repeat, and PH domain I) into an unusual lattice structure that also uncovers a new paradigm for how a BAR protein deforms membrane. Here, we initially pursued computation-based refinement of the ACAP1 lattice to identify its critical protein contacts. Simulation studies then revealed how ACAP1, which dimerizes into a symmetrical structure in solution, is recruited asymmetrically to the membrane through dynamic behavior. We also pursued electron microscopy (EM)-based structural studies, which shed further insight into the dynamic nature of the ACAP1 lattice assembly. As ACAP1 is an unconventional BAR protein, our findings broaden the understanding of the mechanistic spectrum by which proteins assemble into higher-ordered structures to achieve membrane deformation.

Unusual Organization of I-BAR Proteins on Tubular and Vesicular Membranes

Protein-mediated membrane remodeling is a ubiquitous and critical process for proper cellular function. Inverse Bin/Amphiphysin/Rvs (I-BAR) domains drive local membrane deformation as a precursor to large-scale membrane remodeling. We employ a multiscale approach to provide the molecular mechanism of unusual I-BAR domain-driven membrane remodeling at a low protein surface concentration with near-atomistic detail. We generate a bottom-up coarse-grained model that demonstrates similar membrane-bound I-BAR domain aggregation behavior as our recent Mesoscopic Membrane with Explicit Proteins model. Together, these models bridge several length scales and reveal an aggregation behavior of I-BAR domains. We find that at low surface coverage (i.e., low bound protein density), I-BAR domains form transient, tip-to-tip strings on periodic flat membrane sheets. Inside of lipid bilayer tubules, we find linear aggregates parallel to the axis of the tubule. Finally, we find that I-BAR domains form tip-to-tip aggregates around the edges of membrane domes. These results are supported by in vitro experiments showing low curvature bulges surrounded by I-BAR domains on giant unilamellar vesicles. Overall, our models reveal new I-BAR domain aggregation behavior in membrane tubules and on the surface of vesicles at low surface concentration that add insight into how I-BAR domain proteins may contribute to certain aspects of membrane remodeling in cells.

Minimal coarse-grained models for molecular self-organisation in biology

The molecular machinery of life is largely created via self-organisation of individual molecules into functional assemblies. Minimal coarse-grained models, in which a whole macromolecule is represented by a small number of particles, can be of great value in identifying the main driving forces behind self-organisation in cell biology. Such models can incorporate data from both molecular and continuum scales, and their results can be directly compared to experiments. Here we review the state of the art of models for studying the formation and biological function of macromolecular assemblies in living organisms. We outline the key ingredients of each model and their main findings. We illustrate the contribution of this class of simulations to identifying the physical mechanisms behind life and diseases, and discuss their future developments.

Membrane Morphologies Induced by Arc-Shaped Scaffolds Are Determined by Arc Angle and Coverage

The intricate shapes of biological membranes such as tubules and membrane stacks are induced by proteins. In this article, we systematically investigate the membrane shapes induced by arc-shaped scaffolds such as proteins and protein complexes with coarse-grained modeling and simulations. We find that arc-shaped scaffolds induce membrane tubules at membrane coverages larger than a threshold of ∼40%, irrespective of their arc angle. The membrane morphologies at intermediate coverages below this tubulation threshold, in contrast, strongly depend on the arc angle. Scaffolds with arc angles of about 60°, akin to N-BAR domains, do not change the membrane shape at coverages below the tubulation threshold, whereas scaffolds with arc angles larger than about 120° induce double-membrane stacks at intermediate coverages. The scaffolds stabilize the curved membrane edges that connect the membrane stacks, as suggested for complexes of reticulon proteins. Our results provide general insights on the determinants of membrane shaping by arc-shaped scaffolds.

Modeling Membrane Curvature Generation due to Membrane⁻Protein Interactions

To alter and adjust the shape of the plasma membrane, cells harness various mechanisms of curvature generation. Many of these curvature generation mechanisms rely on the interactions between peripheral membrane proteins, integral membrane proteins, and lipids in the bilayer membrane. Mathematical and computational modeling of membrane curvature generation has provided great insights into the physics underlying these processes. However, one of the challenges in modeling these processes is identifying the suitable constitutive relationships that describe the membrane free energy including protein distribution and curvature generation capability. Here, we review some of the commonly used continuum elastic membrane models that have been developed for this purpose and discuss their applications. Finally, we address some fundamental challenges that future theoretical methods need to overcome to push the boundaries of current model applications.

Organizing membrane-curving proteins: the emerging dynamical picture

Lipid membranes play key roles in cells, such as in trafficking, division, infection, remodeling of organelles, among others. The key step in all these processes is creating membrane curvature, typically under the control of many anchored, adhered or included proteins. However, it has become clear that the membrane itself can mediate the interactions among proteins to produce highly ordered assemblies. Computer simulations are ideally suited to investigate protein organization and the dynamics of membrane remodeling at near-micron scales, something that is extremely challenging to tackle experimentally. We review recent computational efforts in modeling protein-caused membrane deformation mechanisms, specifically focusing on coarse-grained simulations. We highlight work that exposed the membrane-mediated ordering of proteins into lines, meshwork, spirals and other assemblies, in what seems to be a very generic mechanism driven by a combination of short and long-ranged forces. Modulating the mechanical properties of membranes is an underexplored signaling mechanism in various processes deserving of more attention in the near future.

The 2018 biomembrane curvature and remodeling roadmap

The importance of curvature as a structural feature of biological membranes has been recognized for many years and has fascinated scientists from a wide range of different backgrounds. On the one hand, changes in membrane morphology are involved in a plethora of phenomena involving the plasma membrane of eukaryotic cells, including endo- and exocytosis, phagocytosis and filopodia formation. On the other hand, a multitude of intracellular processes at the level of organelles rely on generation, modulation, and maintenance of membrane curvature to maintain the organelle shape and functionality. The contribution of biophysicists and biologists is essential for shedding light on the mechanistic understanding and quantification of these processes. Given the vast complexity of phenomena and mechanisms involved in the coupling between membrane shape and function, it is not always clear in what direction to advance to eventually arrive at an exhaustive understanding of this important research area. The 2018 Biomembrane Curvature and Remodeling Roadmap of Journal of Physics D: Applied Physics addresses this need for clarity and is intended to provide guidance both for students who have just entered the field as well as established scientists who would like to improve their orientation within this fascinating area.

Biophysics of membrane curvature remodeling at molecular and mesoscopic lengthscales

At the micron scale, where cell organelles display an amazing complexity in their shape and organization, the physical properties of a biological membrane can be better-understood using continuum models subject to thermal (stochastic) undulations. Yet, the chief orchestrators of these complex and intriguing shapes are a specialized class of membrane associating often peripheral proteins called curvature remodeling proteins (CRPs) that operate at the molecular level through specific protein-lipid interactions. We review multiscale methodologies to model these systems at the molecular as well as at the mesoscopic and cellular scales, and also present a free energy perspective of membrane remodeling through the organization and assembly of CRPs. We discuss the morphological space of nearly planar to highly curved membranes, methods to include thermal fluctuations, and review studies that model such proteins as curvature fields to describe the emergent curved morphologies. We also discuss several mesoscale models applied to a variety of cellular processes, where the phenomenological parameters (such as curvature field strength) are often mapped to models of real systems based on molecular simulations. Much insight can be gained from the calculation of free energies of membranes states with protein fields, which enable accurate mapping of the state and parameter values at which the membrane undergoes morphological transformations such as vesiculation or tubulation. By tuning the strength, anisotropy, and spatial organization of the curvature-field, one can generate a rich array of membrane morphologies that are highly relevant to shapes of several cellular organelles. We review applications of these models to budding of vesicles commonly seen in cellular signaling and trafficking processes such as clathrin mediated endocytosis, sorting by the ESCRT protein complexes, and cellular exocytosis regulated by the exocyst complex. We discuss future prospects where such models can be combined with other models for cytoskeletal assembly, and discuss their role in understanding the effects of cell membrane tension and the mechanics of the extracellular microenvironment on cellular processes.

Membrane re-modelling by BAR domain superfamily proteins via molecular and non-molecular factors

Lipid membranes are structural components of cell surfaces and intracellular organelles. Alterations in lipid membrane shape are accompanied by numerous cellular functions, including endocytosis, intracellular transport, and cell migration. Proteins containing Bin-Amphiphysin-Rvs (BAR) domains (BAR proteins) are unique, because their structures correspond to the membrane curvature, that is, the shape of the lipid membrane. BAR proteins present at high concentration determine the shape of the membrane, because BAR domain oligomers function as scaffolds that mould the membrane. BAR proteins co-operate with various molecular and non-molecular factors. The molecular factors include cytoskeletal proteins such as the regulators of actin filaments and the membrane scission protein dynamin. Lipid composition, including saturated or unsaturated fatty acid tails of phospholipids, also affects the ability of BAR proteins to mould the membrane. Non-molecular factors include the external physical forces applied to the membrane, such as tension and friction. In this mini-review, we will discuss how the BAR proteins orchestrate membrane dynamics together with various molecular and non-molecular factors.

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.

Curvature-Mediated Assembly of Janus Nanoparticles on Membrane Vesicles

Besides direct particle-particle interactions, nanoparticles adsorbed to biomembranes experience indirect interactions that are mediated by the membrane curvature arising from particle adsorption. In this Letter, we show that the curvature-mediated interactions of adsorbed Janus particles depend on the initial curvature of the membrane prior to adsorption, that is, on whether the membrane initially bulges toward or away from the particles in our simulations. The curvature-mediated interaction can be strongly attractive for Janus particles adsorbed to the outside of a membrane vesicle, which initially bulges away from the particles. For Janus particles adsorbed to the vesicle inside, in contrast, the curvature-mediated interactions are repulsive. We find that the area fraction of the adhesive Janus particle surface is an important control parameter for the curvature-mediated interaction and assembly of the particles, besides the initial membrane curvature.

Long-Range Organization of Membrane-Curving Proteins

Biological membranes have a central role in mediating the organization of membrane-curving proteins, a dynamic process that has proven to be challenging to probe experimentally. Using atomic force microscopy, we capture the hierarchically organized assemblies of Bin/amphiphysin/Rvs (BAR) proteins on supported lipid membranes. Their structure reveals distinct long linear aggregates of proteins, regularly spaced by up to 300 nm. Employing accurate free-energy calculations from large-scale coarse-grained computer simulations, we found that the membrane mediates the interaction among protein filaments as a combination of short- and long-ranged interactions. The long-ranged component acts at strikingly long distances, giving rise to a variety of micron-sized ordered patterns. This mechanism may contribute to the long-ranged spatiotemporal control of membrane remodeling by proteins in the cell.

Smectic monolayer confined on a sphere: topology at the particle scale

The impact of topology on the structure of a smectic monolayer confined to a sphere is explored by particle-resolved computer simulations of hard rods. The orientations of the particles are tangential to the sphere and either free or restricted to a prescribed director field with a latitude or longitude orderings. Depending on the imprinted topology, a wealth of different states are found including equatorial smectic with isotropic poles, equatorial smectic with empty poles, a broken egg-shell like modulated smectic, a capped nematic with equatorial bald patches, equatorial nematic with empty poles, and a situation with 4 or 8 half-strength topological defects. Potentially these states could be verified in experiments with Pickering emulsions of droplets with colloidal rods. The unique nature of dipolar structures consisting of positive and negative half-strength disclinations is revealed. These structures, classified by their density and interaction with other defects in the system, relieve the strain of the poles by separating closely positioned half-strength defects. The proximity of these structures to the half-strength defects might enhance the structural diffusion of the defects across the system.

Acceleration and suppression of banana-shaped-protein-induced tubulation by addition of small membrane inclusions of isotropic spontaneous curvatures

Membrane tubulation induced by banana-shaped protein rods is investigated by using coarse-grained meshless membrane simulations. It is found that tubulation is promoted by laterally isotropic membrane inclusions that generate the same sign of spontaneous curvature as the adsorbed protein rods. The inclusions are concentrated in the tubules and reduce the bending energy of the tip of the tubules. On the other hand, inclusions with an opposite curvature suppress tubulation by percolated-network formation at a high protein-rod density while they induce the formation of a spherical membrane bud at a low rod density. When equal amounts of the two types of inclusions (with positive and negative curvatures) are added, their effects cancel each other for the first short period but later the tubulation is slowly accelerated. Positive surface tension suppresses tubulation. Our results suggest that the cooperation of scaffolding of BAR (Bin/Amphiphysin/Rvs) domains and isotropic membrane inclusions is important for tubulation.

Computational and theoretical approaches for studies of a lipid recognition protein on biological membranes

Many cellular functions, including cell signaling and related events, are regulated by the association of peripheral membrane proteins (PMPs) with biological membranes containing anionic lipids, e.g., phosphatidylinositol phosphate (PIP). This association is often mediated by lipid recognition modules present in many PMPs. Here, I summarize computational and theoretical approaches to investigate the molecular details of the interactions and dynamics of a lipid recognition module, the pleckstrin homology (PH) domain, on biological membranes. Multiscale molecular dynamics simulations using combinations of atomistic and coarse-grained models yielded results comparable to those of actual experiments and could be used to elucidate the molecular mechanisms of the formation of protein/lipid complexes on membrane surfaces, which are often difficult to obtain using experimental techniques. Simulations revealed some modes of membrane localization and interactions of PH domains with membranes in addition to the canonical binding mode. In the last part of this review, I address the dynamics of PH domains on the membrane surface. Local PIP clusters formed around the proteins exhibit anomalous fluctuations. This dynamic change in protein-lipid interactions cause temporally fluctuating diffusivity of proteins, i.e., the short-term diffusivity of the bound protein changes substantially with time, and may in turn contribute to the formation/dissolution of protein complexes in membranes.

Systems biology of cellular membranes: a convergence with biophysics

Systems biology and systems medicine have played an important role in the last two decades in shaping our understanding of biological processes. While systems biology is synonymous with network maps and '-omics' approaches, it is not often associated with mechanical processes. Here, we make the case for considering the mechanical and geometrical aspects of biological membranes as a key step in pushing the frontiers of systems biology of cellular membranes forward. We begin by introducing the basic components of cellular membranes, and highlight their dynamical aspects. We then survey the functions of the plasma membrane and the endomembrane system in signaling, and discuss the role and origin of membrane curvature in these diverse cellular processes. We further give an overview of the experimental and modeling approaches to study membrane phenomena. We close with a perspective on the converging futures of systems biology and membrane biophysics, invoking the need to include physical variables such as location and geometry in the study of cellular membranes. WIREs Syst Biol Med 2017, 9:e1386. doi: 10.1002/wsbm.1386 For further resources related to this article, please visit the WIREs website.

The mesoscopic membrane with proteins (MesM-P) model

We present the Mesoscopic Membrane with Proteins (MesM-P) model, an extension of a previously developed elastic membrane model for mesoscale simulations of lipid membranes. MesM-P employs a discrete mesoscopic quasi-particle approach to model protein-facilitated shape and topology changes of the lipid membrane on length and time scales inaccessible to all-atom and quasimolecular coarse-grained molecular dynamics simulations. We investigate the ability of MesM-P to model the behavior of large lipid vesicles as a function of bound protein density. We find four distinct mechanisms for protein aggregation on the surface of the membrane, depending on membrane stiffness and protein spontaneous curvature. We also establish a connection between MesM-P and the results of higher resolution coarse-grained molecular dynamics simulations.

Construction of Nuclear Envelope Shape by a High-Genus Vesicle with Pore-Size Constraint

Nuclear pores have an approximately uniform distribution in the nuclear envelope of most living cells. Hence, the morphology of the nuclear envelope is a spherical stomatocyte with a high genus. We have investigated the morphology of high-genus vesicles under pore-size constraint using dynamically triangulated membrane simulations. Bending-energy minimization without volume or other constraints produces a circular-cage stomatocyte, where the pores are aligned in a circular line on an oblate bud. As the pore radius is reduced, the circular-pore alignment is more stabilized than a random pore distribution on a spherical bud. However, we have clarified the conditions for the formation of a spherical stomatocyte: a small perinuclear volume, osmotic pressure within nucleoplasm, and repulsion between the pores. When area-difference elasticity is taken into account, the formation of cylindrical or budded tubules from the stomatocyte and discoidal stomatocyte is found.

Membrane structure formation induced by two types of banana-shaped proteins

The assembly of banana-shaped rodlike proteins on membranes and the associated membrane shape transformations are investigated by analytical theory and coarse-grained simulations. The membrane-mediated interactions between two banana-shaped inclusions are derived theoretically using a point-like formalism based on fixed anisotropic curvatures, both for zero surface tension and for finite surface tension. On a larger scale, the interactions between the assemblies of such rodlike inclusions are determined analytically. Meshless membrane simulations are performed in the presence of a large number of inclusions of two types, corresponding to the curved rods of opposite curvatures, both for flat membranes and vesicles. Rods of the same type aggregate into linear assemblies perpendicular to the rod axis, leading to membrane tubulation. However, rods of the other type, those of opposite curvature, are attracted to the lateral sides of these assemblies, and stabilize a straight bump structure that prevents tubulation. When the two types of rods have almost opposite curvatures, the bumps attract one another, forming a striped structure. Positive surface tension is found to stabilize stripe formation. The simulation results agree well with the theoretical predictions provided the point-like curvatures of the model are scaled-down to account for the effective flexibility of the simulated rods.

On the role of external force of actin filaments in the formation of tubular protrusions of closed membrane shapes with anisotropic membrane components

Biological membranes are composed of different components and there is no a priori reason to assume that all components are isotropic. It was previously shown that the anisotropic properties of membrane components may explain the stability of membrane tubular protrusions even without the application of external force. Our theoretical study focuses on the role of anisotropic membrane components in the stability of membrane tubular structures generated or stabilized by actin filaments. We show that the growth of the actin cytoskeleton inside the vesicle can induce the partial lateral segregation of different membrane components. The entropy of mixing of membrane components hinders the total lateral segregation of the anisotropic and isotropic membrane components. Self-assembled aggregates formed by anisotropic membrane components facilitate the growth of long membrane tubular protrusions. Protrusive force generated by actin filaments favors strong segregation of membrane components by diminishing the opposing effect of mixing entropy.

The Matrix protein M1 from influenza C virus induces tubular membrane invaginations in an in vitro cell membrane model

Matrix proteins from enveloped viruses play an important role in budding and stabilizing virus particles. In order to assess the role of the matrix protein M1 from influenza C virus (M1-C) in plasma membrane deformation, we have combined structural and in vitro reconstitution experiments with model membranes. We present the crystal structure of the N-terminal domain of M1-C and show by Small Angle X-Ray Scattering analysis that full-length M1-C folds into an elongated structure that associates laterally into ring-like or filamentous polymers. Using negatively charged giant unilamellar vesicles (GUVs), we demonstrate that M1-C full-length binds to and induces inward budding of membrane tubules with diameters that resemble the diameter of viruses. Membrane tubule formation requires the C-terminal domain of M1-C, corroborating its essential role for M1-C polymerization. Our results indicate that M1-C assembly on membranes constitutes the driving force for budding and suggest that M1-C plays a key role in facilitating viral egress.

How curvature-generating proteins build scaffolds on membrane nanotubes

Bin/Amphiphysin/Rvs (BAR) domain proteins control the curvature of lipid membranes in endocytosis, trafficking, cell motility, the formation of complex subcellular structures, and many other cellular phenomena. They form 3D assemblies that act as molecular scaffolds to reshape the membrane and alter its mechanical properties. It is unknown, however, how a protein scaffold forms and how BAR domains interact in these assemblies at protein densities relevant for a cell. In this work, we use various experimental, theoretical, and simulation approaches to explore how BAR proteins organize to form a scaffold on a membrane nanotube. By combining quantitative microscopy with analytical modeling, we demonstrate that a highly curving BAR protein endophilin nucleates its scaffolds at the ends of a membrane tube, contrary to a weaker curving protein centaurin, which binds evenly along the tube's length. Our work implies that the nature of local protein-membrane interactions can affect the specific localization of proteins on membrane-remodeling sites. Furthermore, we show that amphipathic helices are dispensable in forming protein scaffolds. Finally, we explore a possible molecular structure of a BAR-domain scaffold using coarse-grained molecular dynamics simulations. Together with fluorescence microscopy, the simulations show that proteins need only to cover 30-40% of a tube's surface to form a rigid assembly. Our work provides mechanical and structural insights into the way BAR proteins may sculpt the membrane as a high-order cooperative assembly in important biological processes.

Physical basis of some membrane shaping mechanisms

In vesicular transport pathways, membrane proteins and lipids are internalized, externalized or transported within cells, not by bulk diffusion of single molecules, but embedded in the membrane of small vesicles or thin tubules. The formation of these 'transport carriers' follows sequential events: membrane bending, fission from the donor compartment, transport and eventually fusion with the acceptor membrane. A similar sequence is involved during the internalization of drug or gene carriers inside cells. These membrane-shaping events are generally mediated by proteins binding to membranes. The mechanisms behind these biological processes are actively studied both in the context of cell biology and biophysics. Bin/amphiphysin/Rvs (BAR) domain proteins are ideally suited for illustrating how simple soft matter principles can account for membrane deformation by proteins. We review here some experimental methods and corresponding theoretical models to measure how these proteins affect the mechanics and the shape of membranes. In more detail, we show how an experimental method employing optical tweezers to pull a tube from a giant vesicle may give important quantitative insights into the mechanism by which proteins sense and generate membrane curvature and the mechanism of membrane scission.This article is part of the themed issue 'Soft interfacial materials: from fundamentals to formulation'.

Multiscale simulations of protein-facilitated membrane remodeling

Protein-facilitated shape and topology changes of cell membranes are crucial for many biological processes, such as cell division, protein trafficking, and cell signaling. However, the inherently multiscale nature of membrane remodeling presents a considerable challenge for understanding the mechanisms and physics that drive this process. To address this problem, a multiscale approach that makes use of a diverse set of computational and experimental techniques is required. The atomistic simulations provide high-resolution information on protein-membrane interactions. Experimental techniques, like electron microscopy, on the other hand, resolve high-order organization of proteins on the membrane. Coarse-grained (CG) and mesoscale computational techniques provide the intermediate link between the two scales and can give new insights into the underlying mechanisms. In this Review, we present the recent advances in multiscale computational approaches established in our group. We discuss various CG and mesoscale approaches in studying the protein-mediated large-scale membrane remodeling.