Anisotropic spontaneous curvatures in lipid membranes

Glycocalyx-induced formation of membrane tubes

Tubular membrane structures are ubiquitous in cells and in the membranes of intracellular organelles such as the Golgi complex and the endoplasmic reticulum. Tubulation plays essential roles in numerous biological processes, including filopodia growth, trafficking, ion transport, and cellular motility. Understanding the fundamental mechanism of the formation of membrane tubes is thus an important problem in the fields of biology and biophysics. Although extensive studies have shown that tubes can be formed due to localized forces acting on the membrane or by the spontaneous curvature induced by membrane-bound proteins, little is known about how membrane tubes are induced by glycocalyx, a sugar-rich layer at the cell surface. In this work, we develop a biophysical model that combines polymer physics theory and the Canham-Helfrich membrane theory to investigate how the glycocalyx generates cylindrical tubular protrusions on the cell membrane. Our results show that the glycocalyx alone can induce the formation of tubular membrane structures. This tube formation involves a first-order shape transition without any externally applied force or other curvature-inducing mechanisms. We also find there exist critical values of glycocalyx grafting density and glycopolymer length needed to induce the formation of tubular structures. The presence of a vertical actin force, line tension, and spontaneous curvature reduce this critical grafting density and length of polymer that triggers the formation of membrane tube, which suggests that the glycocalyx makes tube formation energetically more favorable when combined with an actin force, line tension, and spontaneous curvature.

Biophysical modeling of membrane curvature generation and curvature sensing by the glycocalyx

Generation of membrane curvature is fundamental to cellular function. Recent studies have established that the glycocalyx, a sugar-rich polymer layer at the cell surface, can generate membrane curvature. While there have been some theoretical efforts to understand the interplay between the glycocalyx and membrane bending, there remain open questions about how the properties of the glycocalyx affect membrane bending. For example, the relationship between membrane curvature and the density of glycosylated proteins on its surface remains unclear. In this work, we use polymer brush theory to develop a detailed biophysical model of the energetic interactions of the glycocalyx with the membrane. Using this model, we identify the conditions under which the glycocalyx can both generate and sense curvature. Our model predicts that the extent of membrane curvature generated depends on the grafting density of the glycocalyx and the backbone length of the polymers constituting the glycocalyx. Furthermore, when coupled with the intrinsic membrane properties such as spontaneous curvature and a line tension along the membrane, the curvature generation properties of the glycocalyx are enhanced. These predictions were tested experimentally by examining the propensity of glycosylated transmembrane proteins to drive the assembly of highly curved filopodial protrusions at the plasma membrane of adherent mammalian cells. Our model also predicts that the glycocalyx has curvature-sensing capabilities, in agreement with the results of our experiments. Thus, our study develops a quantitative framework for mapping the properties of the glycocalyx to the curvature generation capability of the membrane.

Glycocalyx-induced formation of membrane tubes

Tubular membrane structures are ubiquitous in cells and in the membranes of intracellular organelles such as the Golgi complex and the endoplasmic reticulum. Tubulation plays essential roles in numerous biological processes, including filopodia growth, trafficking, ion transport, and cellular motility. Understanding the fundamental mechanism of the formation of membrane tubes is thus an important problem in the fields of biology and biophysics. Though extensive studies have shown that tubes can be formed due to localized forces acting on the membrane or by the curvature induced by membrane-bound proteins, little is known about how membrane tubes are induced by glycocalyx, a sugar-rich layer at the cell surface. In this work, we develop a biophysical model that combines polymer physics theory and the Canham-Helfrich membrane theory to investigate how the glycocalyx generates cylindrical tubular protrusions on the cell membrane. Our results show that the glycocalyx alone can induce the formation of tubular membrane structures. This tube formation involves a first-order shape transition without any externally applied force or other curvature-inducing mechanisms. We also find that critical values of glycocalyx grafting density and glycopolymer length are needed to induce the formation of tubular structures. The presence of vertical actin force, line tension, and spontaneous curvature reduces the critical grafting density and length of polymer that triggers the formation of membrane tube, which suggests that the glycocalyx makes tube formation energetically more favorable when combined with an actin force, line tension, and spontaneous curvature.

Significance Statement

In many cells, the existence of glycocalyx, a thick layer of polymer meshwork comprising proteins and complex sugar chains coating the outside of the cell membrane, regulates the formation of membrane tubes. Here, we propose a theoretical model that combines polymer physics theory and the Canham-Helfrich membrane theory to study the formation of cylindrical tubular protrusions induced by the glycocalyx. Our findings indicate that glycocalyx plays an important role in the formation of membrane tubes. We find that there exists critical grafting density and length of polymer that triggers the formation of membrane tubes, and the glycocalyx-induced tube formation is facilitated when combined with actin forces, line tension, and spontaneous curvature. Our theoretical model has implications for understanding how biological membranes may form tubular structures.

Dimensionality reduction in bulk-boundary reaction-diffusion systems

Intracellular protein patterns regulate many vital cellular functions, such as the processing of spatiotemporal information or the control of shape deformations. To do so, pattern-forming systems can be sensitive to the cell geometry by means of coupling the protein dynamics on the cell membrane to dynamics in the cytosol. Recent studies demonstrated that modeling the cytosolic dynamics in terms of an averaged protein pool disregards possibly crucial aspects of the pattern formation, most importantly concentration gradients normal to the membrane. At the same time, the coupling of two domains (surface and volume) with different dimensions renders many standard tools for the numerical analysis of self-organizing systems inefficient. Here, we present a generic framework for projecting the cytosolic dynamics onto the lower-dimensional surface that respects the influence of cytosolic concentration gradients in static and evolving geometries. This method uses a priori physical information about the system to approximate the cytosolic dynamics by a small number of dominant characteristic concentration profiles (basis), akin to basis transformations of finite element methods. As a proof of concept, we apply our framework to a toy model for volume-dependent interrupted coarsening, evaluate the accuracy of the results for various basis choices, and discuss the optimal basis choice for biologically relevant systems. Our analysis presents an efficient yet accurate method for analyzing pattern formation with surface-volume coupling in evolving geometries.

Dynamins combine mechano-constriction and membrane remodeling to enable two-step mitochondrial fission via a 'snap-through' instability

Mitochondrial fission is controlled by dynamin proteins, the dysregulation of which is correlated with diverse diseases. Fission dynamins are GTP hydrolysis-driven mechanoenzymes that self-oligomerize into helical structures that constrict membrane to achieve fission, but details are not well understood. However, dynamins can also remodel membranes by inducing negative Gaussian curvature, the type of curvature required for completion of fission. Here, we examine how these drastically different mechanisms synergistically exert their influences on a membrane, via a mechanical model calibrated with small-angle X-ray scattering structural data. We find that free dynamin can trigger a "snap-through instability" that enforces a shape transition from an oligomer-confined cylindrical membrane to a drastically narrower catenoid-shaped neck within the spontaneous hemi-fission regime, in a manner that depends critically on the length of the confined tube. These results indicate how the combination of dynamin assembly, and paradoxically disassembly, can lead to diverse pathways to scission.

Comparing Multifunctional Viral and Eukaryotic Proteins for Generating Scission Necks in Membranes

Deterministic formation of membrane scission necks by protein machinery with multiplexed functions is critical in biology. A microbial example is M2 viroporin, a proton pump from the influenza A virus that is multiplexed with membrane remodeling activity to induce budding and scission in the host membrane during viral maturation. In comparison, the dynamin family constitutes a class of eukaryotic proteins implicated in mitochondrial fission, as well as various budding and endocytosis pathways. In the case of Dnm1, the mitochondrial fission protein in yeast, the membrane remodeling activity is multiplexed with mechanoenzyme activity to create fission necks. It is not clear why these functions are combined in these scission processes, which occur in drastically different compositions and solution conditions. In general, direct experimental access to changing neck sizes induced by individual proteins or peptide fragments is challenging due to the nanoscale dimensions and influence of thermal fluctuations. Here, we use a mechanical model to estimate the size of scission necks by leveraging small-angle X-ray scattering structural data of protein-lipid systems under different conditions. The influence of interfacial tension, lipid composition, and membrane budding morphology on the size of the induced scission necks is systematically investigated using our data and molecular dynamic simulations. We find that the M2 budding protein from the influenza A virus has robust pH-dependent membrane activity that induces nanoscopic necks within the range of spontaneous hemifission for a broad range of lipid compositions. In contrast, the sizes of scission necks generated by mitochondrial fission proteins strongly depend on lipid composition, which suggests a role for mechanical constriction.

Inferring biophysical properties of membranes during endocytosis using machine learning

Endocytosis is a fundamental cellular process in eukaryotic cells that facilitates the transport of molecules into the cell. With the help of fluorescence microscopy and electron tomography, researchers have accumulated extensive geometric data of membrane shapes during endocytosis. These data contain rich information about the mechanical properties of membranes, which are hard to access via experiments due to the small dimensions of the endocytic patch. In this study, we propose an approach that combines machine learning with the Helfrich theory of membranes to infer the mechanical properties of membranes during endocytosis from a dataset of membrane shapes extracted from electron tomography. Our results demonstrate that machine learning can output solutions that both match the experimental profile and satisfy the membrane shape equations derived from Helfrich theory. The learning results show that during the early stage of endocytosis, the inferred membrane tension is negative, indicating the presence of strong compressive forces at the boundary of the endocytic invagination. Our method presents a generic framework for extracting membrane information from super-resolution imaging.

Formation of protein-mediated bilayer tubes is governed by a snapthrough transition

Plasma membrane tubes are ubiquitous in cellular membranes and in the membranes of intracellular organelles. They play crucial roles in trafficking, ion transport, and cellular motility. These tubes can be formed due to localized forces acting on the membrane or by the curvature induced by membrane-bound proteins. Here, we present a mathematical framework to model cylindrical tubular protrusions formed by proteins that induce anisotropic spontaneous curvature. Our analysis revealed that the tube radius depends on an effective tension that includes contributions from the bare membrane tension and the protein-induced curvature. We also found that the length of the tube undergoes an abrupt transition from a short, dome-shaped membrane to a long cylinder and this transition is characteristic of a snapthrough instability. Finally, we show that the snapthrough instability depends on the different parameters including coat area, bending modulus, and extent of protein-induced curvature. Our findings have implications for tube formation due to BAR-domain proteins in processes such as endocytosis, t-tubule formation in myocytes, and cristae formation in mitochondria.

Molecular Shape Solution for Mesoscopic Remodeling of Cellular Membranes

Cellular membranes self-assemble from and interact with various molecular species. Each molecule locally shapes the lipid bilayer, the soft elastic core of cellular membranes. The dynamic architecture of intracellular membrane systems is based on elastic transformations and lateral redistribution of these elementary shapes, driven by chemical and curvature stress gradients. The minimization of the total elastic stress by such redistribution composes the most basic, primordial mechanism of membrane curvature-composition coupling (CCC). Although CCC is generally considered in the context of dynamic compositional heterogeneity of cellular membrane systems, in this article we discuss a broader involvement of CCC in controlling membrane deformations. We focus specifically on the mesoscale membrane transformations in open, reservoir-governed systems, such as membrane budding, tubulation, and the emergence of highly curved sites of membrane fusion and fission. We reveal that the reshuffling of molecular shapes constitutes an independent deformation mode with complex rheological properties.This mode controls effective elasticity of local deformations as well as stationary elastic stress, thus emerging as a major regulator of intracellular membrane remodeling.

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.

Mechanical Principles Governing the Shapes of Dendritic Spines

Dendritic spines are small, bulbous protrusions along the dendrites of neurons and are sites of excitatory postsynaptic activity. The morphology of spines has been implicated in their function in synaptic plasticity and their shapes have been well-characterized, but the potential mechanics underlying their shape development and maintenance have not yet been fully understood. In this work, we explore the mechanical principles that could underlie specific shapes using a minimal biophysical model of membrane-actin interactions. Using this model, we first identify the possible force regimes that give rise to the classic spine shapes-stubby, filopodia, thin, and mushroom-shaped spines. We also use this model to investigate how the spine neck might be stabilized using periodic rings of actin or associated proteins. Finally, we use this model to predict that the cooperation between force generation and ring structures can regulate the energy landscape of spine shapes across a wide range of tensions. Thus, our study provides insights into how mechanical aspects of actin-mediated force generation and tension can play critical roles in spine shape maintenance.

The Mechanics and Thermodynamics of Tubule Formation in Biological Membranes

Membrane tubulation is a ubiquitous process that occurs both at the plasma membrane and on the membranes of intracellular organelles. These tubulation events are known to be mediated by forces applied on the membrane either due to motor proteins, by polymerization of the cytoskeleton, or due to the interactions between membrane proteins binding onto the membrane. The numerous experimental observations of tube formation have been amply supported by mathematical modeling of the associated membrane mechanics and have provided insights into the force-displacement relationships of membrane tubes. Recent advances in quantitative biophysical measurements of membrane-protein interactions and tubule formation have necessitated the need for advances in modeling that will account for the interplay of multiple aspects of physics that occur simultaneously. Here, we present a comprehensive review of experimental observations of tubule formation and provide context from the framework of continuum modeling. Finally, we explore the scope for future research in this area with an emphasis on iterative modeling and experimental measurements that will enable us to expand our mechanistic understanding of tubulation processes in cells.

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.

Non-uniform distribution of myosin-mediated forces governs red blood cell membrane curvature through tension modulation

The biconcave disk shape of the mammalian red blood cell (RBC) is unique to the RBC and is vital for its circulatory function. Due to the absence of a transcellular cytoskeleton, RBC shape is determined by the membrane skeleton, a network of actin filaments cross-linked by spectrin and attached to membrane proteins. While the physical properties of a uniformly distributed actin network interacting with the lipid bilayer membrane have been assumed to control RBC shape, recent experiments reveal that RBC biconcave shape also depends on the contractile activity of nonmuscle myosin IIA (NMIIA) motor proteins. Here, we use the classical Helfrich-Canham model for the RBC membrane to test the role of heterogeneous force distributions along the membrane and mimic the contractile activity of sparsely distributed NMIIA filaments. By incorporating this additional contribution to the Helfrich-Canham energy, we find that the RBC biconcave shape depends on the ratio of forces per unit volume in the dimple and rim regions of the RBC. Experimental measurements of NMIIA densities at the dimple and rim validate our prediction that (a) membrane forces must be non-uniform along the RBC membrane and (b) the force density must be larger in the dimple than the rim to produce the observed membrane curvatures. Furthermore, we predict that RBC membrane tension and the orientation of the applied forces play important roles in regulating this force-shape landscape. Our findings of heterogeneous force distributions on the plasma membrane for RBC shape maintenance may also have implications for shape maintenance in different cell types.

The spectrin-actin network of the membrane skeleton plays an important role in controlling specialized cell membrane morphology. In the paradigmatic red blood cell (RBC), where actin filaments are present exclusively in the membrane skeleton, recent experiments reveal that nonmuscle myosin IIA (NMIIA) motor contractility maintains the RBC biconcave disk shape. In this study, we have identified criteria for micron-scale distributions of NMIIA forces at the membrane required to maintain the biconcave disk shape of an RBC in the resting condition. Supported by experimental measurements of RBC NMIIA distribution, we showed that a heterogeneous force distribution with a larger force density at the dimple is able to capture the experimentally observed biconcave morphology of an RBC with better accuracy compared to previous models that did not consider the heterogeneity in the force distribution. Furthermore, we showed that the biconcave geometry of the RBC is closely regulated by the effective membrane tension and the direction of applied forces on the membrane. These findings can be generalized to any force-mediated membrane shape, providing insight into the role of actomyosin forces in prescribing and maintaining the morphology of different cell types.

Numerical Model for the Determination of Erythrocyte Mechanical Properties and Wall Shear Stress in vivo From Intravital Microscopy

The mechanical properties and deformability of Red Blood Cells (RBCs) are important determinants of blood rheology and microvascular hemodynamics. The objective of this study is to quantify the mechanical properties and wall shear stress experienced by the RBC membrane during capillary plug flow in vivo utilizing high speed video recording from intravital microscopy, biomechanical modeling, and computational methods. Capillaries were imaged in the rat cremaster muscle pre- and post-RBC transfusion of stored RBCs for 2-weeks. RBC membrane contours were extracted utilizing image processing and parametrized. RBC parameterizations were used to determine updated deformation gradient and Lagrangian Green strain tensors for each point along the parametrization and for each frame during plug flow. The updated Lagrangian Green strain and Displacement Gradient tensors were numerically fit to the Navier-Lame equations along the parameterized boundary to determined Lame's constants. Mechanical properties and wall shear stress were determined before and transfusion, were grouped in three populations of erythrocytes: native cells (NC) or circulating cells before transfusion, and two distinct population of cells after transfusion with stored cells (SC1 and SC2). The distinction, between the heterogeneous populations of cells present after the transfusion, SC1 and SC2, was obtained through principle component analysis (PCA) of the mechanical properties along the membrane. Cells with the first two principle components within 3 standard deviations of the mean, were labeled as SC1, and those with the first two principle components greater than 3 standard deviations from the mean were labeled as SC2. The calculated shear modulus average was 1.1±0.2, 0.90±0.15, and 12 ± 8 MPa for NC, SC1, and SC2, respectively. The calculated young's modulus average was 3.3±0.6, 2.6±0.4, and 32±20 MPa for NC, SC1, and SC2, respectively. o our knowledge, the methods presented here are the first estimation of the erythrocyte mechanical properties and shear stress in vivo during capillary plug flow. In summary, the methods introduced in this study may provide a new avenue of investigation of erythrocyte mechanics in the context of hematologic conditions that adversely affect erythrocyte mechanical properties.

Principles of self-organization and load adaptation by the actin cytoskeleton during clathrin-mediated endocytosis

Force generation by actin assembly shapes cellular membranes. An experimentally constrained multiscale model shows that a minimal branched actin network is sufficient to internalize endocytic pits against membrane tension. Around 200 activated Arp2/3 complexes are required for robust internalization. A newly developed molecule-counting method determined that ~200 Arp2/3 complexes assemble at sites of clathrin-mediated endocytosis in human cells. Simulations predict that actin self-organizes into a radial branched array with growing ends oriented toward the base of the pit. Long actin filaments bend between attachment sites in the coat and the base of the pit. Elastic energy stored in bent filaments, whose presence was confirmed by cryo-electron tomography, contributes to endocytic internalization. Elevated membrane tension directs more growing filaments toward the base of the pit, increasing actin nucleation and bending for increased force production. Thus, spatially constrained actin filament assembly utilizes an adaptive mechanism enabling endocytosis under varying physical constraints.

The outer membrane of a cell is a tight but elastic barrier that controls what enters or leaves the cell. Large molecules typically cannot cross this membrane unaided. Instead, to enter the cell, they must be packaged into a pocket of the membrane that is then pulled inside. This process, called endocytosis, shuttles material into a cell hundreds of times a minute.

Endocytosis relies on molecular machines that assemble and disassemble at the membrane as required. One component, a protein called actin, self-assembles near the membrane into long filaments with many repeated subunits. These filaments grow against the membrane, pulling it inwards. But it was not clear how actin filaments organize in such a way that allows them to pull on the membrane with enough force – and without a template to follow.

Akamatsu et al. set about identifying how actin operates during endocytosis by using computer simulations that were informed by measurements made in living cells. The simulations included information about the location of actin and other essential molecules, along with the details of how these molecules work individually and together. Akamatsu et al. also developed a method to count the numbers of molecules of a key protein at individual sites of endocytosis. High-resolution imaging was then used to create 3D pictures of actin and endocytosis in action in human cells grown in the laboratory.

The analysis showed the way actin filaments arrange themselves depends on the starting positions of a few key molecules that connect to actin. Imaging confirmed that, like a pole-vaulting pole, the flexible actin filaments bend to store energy and then release it to pull the membrane inwards during endocytosis. Finally, the simulations predicted that the collection of filaments adapts its shape and size in response to the resistance of the elastic membrane. This makes the system opportunistic and adaptable to the unpredictable environment within cells.

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.

Guided by curvature: shaping cells by coupling curved membrane proteins and cytoskeletal forces

Eukaryote cells have flexible membranes that allow them to have a variety of dynamical shapes. The shapes of the cells serve important biological functions, both for cells within an intact tissue, and during embryogenesis and cellular motility. How cells control their shapes and the structures that they form on their surface has been a subject of intensive biological research, exposing the building blocks that cells use to deform their membranes. These processes have also drawn the interest of theoretical physicists, aiming to develop models based on physics, chemistry and nonlinear dynamics. Such models explore quantitatively different possible mechanisms that the cells can employ to initiate the spontaneous formation of shapes and patterns on their membranes. We review here theoretical work where one such class of mechanisms was investigated: the coupling between curved membrane proteins, and the cytoskeletal forces that they recruit. Theory indicates that this coupling gives rise to a rich variety of membrane shapes and dynamics, while experiments indicate that this mechanism appears to drive many cellular shape changes.This article is part of the theme issue 'Self-organization in cell biology'.

Geometric instability catalyzes mitochondrial fission

The mitochondrial membrane undergoes extreme remodeling during fission. While a few membrane-squeezing proteins are recognized as the key drivers of fission, there is a growing body of evidence that strongly suggests that conical lipids play a critical role in regulating mitochondrial morphology and fission. However, the mechanisms by which proteins and lipids cooperate to execute fission have not been quantitatively investigated. Here, we computationally model the squeezing of the largely tubular mitochondrion and show that proteins and conical lipids can act synergistically to trigger buckling instability and achieve extreme constriction. More remarkably, the study reveals that the conical lipids can act with different fission proteins to induce hierarchical instabilities and create increasingly narrow and stable constrictions. We reason that this geometric plasticity imparts significant robustness to the fission reaction by arresting the elastic tendency of the membrane to rebound during protein polymerization and depolymerization cycles. Our in vitro study validates protein–lipid cooperativity in constricting membrane tubules. Overall, our work presents a general mechanism for achieving drastic topological remodeling in cellular membranes.

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.

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.

Curvature-Induced Spatial Ordering of Composition in Lipid Membranes

Phase segregation of membranal components, such as proteins, lipids, and cholesterols, leads to the formation of aggregates or domains that are rich in specific constituents. This process is important in the interaction of the cell with its surroundings and in determining the cell's behavior and fate. Motivated by published experiments on curvature-modulated phase separation in lipid membranes, we formulate a mathematical model aiming at studying the spatial ordering of composition in a two-component biomembrane that is subjected to a prescribed (imposed) geometry. Based on this model, we identified key nondimensional quantities that govern the biomembrane response and performed numerical simulations to quantitatively explore their influence. We reproduce published experimental observations and extend them to surfaces with geometric features (imposed geometry) and lipid phases beyond those used in the experiments. In addition, we demonstrate the possibility for curvature-modulated phase separation above the critical temperature and propose a systematic procedure to determine which mechanism, the difference in bending stiffness or difference in spontaneous curvatures of the two phases, dominates the coupling between shape and composition.

Anisotropic Membrane Curvature Sensing by Amphipathic Peptides

Many proteins and peptides have an intrinsic capacity to sense and induce membrane curvature, and play crucial roles for organizing and remodeling cell membranes. However, the molecular driving forces behind these processes are not well understood. Here, we describe an approach to study curvature sensing by simulating the interactions of single molecules with a buckled lipid bilayer. We analyze three amphipathic antimicrobial peptides, a class of membrane-associated molecules that specifically target and destabilize bacterial membranes, and find qualitatively different sensing characteristics that would be difficult to resolve with other methods. Our findings provide evidence for direction-dependent curvature sensing mechanisms in amphipathic peptides and challenge existing theories of hydrophobic insertion. The buckling approach is generally applicable to a wide range of curvature-sensing molecules, and our results provide strong motivation to develop new experimental methods to track position and orientation of membrane proteins.

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'.

Thermal fluctuations of vesicles and nonlinear curvature elasticity--implications for size-dependent renormalized bending rigidity and vesicle size distribution

Both closed and open biological membranes noticeably undulate at physiological temperatures. These thermal fluctuations influence a broad range of biophysical phenomena, ranging from self-assembly to adhesion. In particular, the experimentally measured thermal fluctuation spectra also provide a facile route to the assessment of mechanical and certain other physical properties of biological membranes. The theoretical assessment of thermal fluctuations, be it for closed vesicles or the simpler case of flat open lipid bilayers, is predicated upon assuming that the elastic curvature energy is a quadratic functional of the curvature tensor. However, a qualitatively correct description of several phenomena such as binding-unbinding transition, vesicle-to-bicelle transition, appearance of hats and saddles among others, appears to require consideration of constitutively nonlinear elasticity that includes fourth order curvature contributions rather than just quadratic. In particular, such nonlinear considerations are relevant in the context of large-curvature or small-sized vesicles. In this work we discuss the statistical mechanics of closed membranes (vesicles) incorporating both constitutive and geometrical nonlinearities. We derive results for the renormalized bending rigidity of small vesicles and show that significant stiffening may occur for sub-20 nm vesicle sizes. Our closed-form results may also be used to determine nonlinear curvature elasticity properties from either experimentally measured fluctuation spectra or microscopic calculations such as molecular dynamics. Finally, in the context of our results on thermal fluctuations of vesicles and nonlinear curvature elasticity, we reexamine the problem of determining the size distribution of vesicles and obtain results that reconcile well with experimental observations. However, our results are somewhat paradoxical. Specifically, the molecular dynamics predictions for the thermo-mechanical behavior of small vesicles of prior studies appear to be inconsistent with the nonlinear elastic properties that we estimate by fitting to the experimentally determined vesicle size-distribution trends and data.

Endocytic proteins drive vesicle growth via instability in high membrane tension environment

Clathrin-mediated endocytosis (CME) is a key pathway for transporting cargo into cells via membrane vesicles; it plays an integral role in nutrient import, signal transduction, neurotransmission, and cellular entry of pathogens and drug-carrying nanoparticles. Because CME entails substantial local remodeling of the plasma membrane, the presence of membrane tension offers resistance to bending and hence, vesicle formation. Experiments show that in such high-tension conditions, actin dynamics is required to carry out CME successfully. In this study, we build on these pioneering experimental studies to provide fundamental mechanistic insights into the roles of two key endocytic proteins-namely, actin and BAR proteins-in driving vesicle formation in high membrane tension environment. Our study reveals an actin force-induced "snap-through instability" that triggers a rapid shape transition from a shallow invagination to a highly invaginated tubular structure. We show that the association of BAR proteins stabilizes vesicles and induces a milder instability. In addition, we present a rather counterintuitive role of BAR depolymerization in regulating the shape evolution of vesicles. We show that the dissociation of BAR proteins, supported by actin-BAR synergy, leads to considerable elongation and squeezing of vesicles. Going beyond the membrane geometry, we put forth a stress-based perspective for the onset of vesicle scission and predict the shapes and composition of detached vesicles. We present the snap-through transition and the high in-plane stress as possible explanations for the intriguing direct transformation of broad and shallow invaginations into detached vesicles in BAR mutant yeast cells.

Membrane tension and peripheral protein density mediate membrane shape transitions

Endocytosis is a ubiquitous eukaryotic membrane budding, vesiculation and internalization process fulfilling numerous roles including compensation of membrane area increase after bursts of exocytosis. The mechanism of the coupling between these two processes to enable homeostasis is not well understood. Recently, an ultrafast endocytosis (UFE) pathway was revealed with a speed significantly exceeding classical clathrin-mediated endocytosis (CME). Membrane tension reduction is a potential mechanism by which endocytosis can be rapidly activated at remote sites. Here, we provide experimental evidence for a mechanism whereby membrane tension reduction initiates membrane budding and tubulation mediated by endocytic proteins, such as endophilin A1. We find that shape instabilities occur at well-defined membrane tensions and surface densities of endophilin. From our data, we obtain a membrane shape stability diagram that shows remarkable consistency with a quantitative model. This model applies to all laterally diffusive curvature-coupling proteins and therefore a wide range of endocytic proteins.

Phenomenology based multiscale models as tools to understand cell membrane and organelle morphologies

An intriguing question in cell biology is "how do cells regulate their shape?" It is commonly believed that the observed cellular morphologies are a result of the complex interaction among the lipid molecules (constituting the cell membrane), and with a number of other macromolecules, such as proteins. It is also believed that the common biophysical processes essential for the functioning of a cell also play an important role in cellular morphogenesis. At the cellular scale-where typical dimensions are in the order of micrometers-the effects arising from the molecular scale can either be modeled as equilibrium or non-equilibrium processes. In this chapter, we discuss the dynamically triangulated Monte Carlo technique to model and simulate membrane morphologies at the cellular scale, which in turn can be used to investigate several questions related to shape regulation in cells. In particular, we focus on two specific problems within the framework of isotropic and anisotropic elasticity theories: namely, (i) the origin of complex, physiologically relevant, membrane shapes due to the interaction of the membrane with curvature remodeling proteins, and (ii) the genesis of steady state cellular shapes due to the action of non-equilibrium forces that are generated by the fission and fusion of transport vesicles and by the binding and unbinding of proteins from the parent membrane.