Budding of crystalline domains in fluid membranes

Curvature and confinement effects on chiral liquid crystal morphologies

Chiral liquid crystals (ChLCs) exhibit an inherent twist that originates at the molecular scale and can extend over multiple length scales when unconstrained. Under confinement, the twist is thwarted, leading to formation of defects in the molecular order that offer distinct optical responses and opportunities for colloidal driven assembly. Past studies have explored spheroidal confinement down to the nanoscopic regime, where curved boundaries produce surface defects to accommodate topological constraints and restrict the propagation of cuboidal defect networks. Similarly, strict confinement in channels and shells has been shown to give rise to escaped configurations and skyrmions. However, little is known about the role of extrinsic curvature in the development of cholesteric textures and Blue Phases (BP). In this paper, we examine the palette of morphologies that arises when ChLCs are confined in toroidal and cylindrical cavities. The equilibrium morphologies are obtained following an annealing strategy of a Landau-de Gennes free energy functional. Three dimensionless groups are identified to build phase diagrams: the natural twist, the ratio of elastic energies, and the circumscription of a BP cell. Curvature is shown to introduce helical features that are first observed as a Double Twist, and progress to Chiral Ribbons and, ultimately, Helical BP and BP. Chiral ribbons are examined as useful candidates for driven assembly given their tunability and robustness.

Dynamics of fluid bilayer vesicles: Soft meshes and robust curvature energy discretization

Continuum models like the Helfrich Hamiltonian are widely used to describe fluid bilayer vesicles. Here we study the molecular dynamics compatible dynamics of the vertices of two-dimensional meshes representing the bilayer, whose in-plane motion is only weakly constrained. We show (i) that Jülicher's discretization of the curvature energy offers vastly superior robustness for soft meshes compared to the commonly employed expression by Gommper and Kroll and (ii) that for sufficiently soft meshes, the typical behavior of fluid bilayer vesicles can emerge even if the mesh connectivity remains fixed throughout the simulations. In particular, soft meshes can accommodate large shape transformations, and the model can generate the typical ℓ^{-4} signal for the amplitude of surface undulation modes of nearly spherical vesicles all the way up to the longest wavelength modes. Furthermore, we compare results for Newtonian, Langevin, and Brownian dynamics simulations of the mesh vertices to demonstrate that the internal friction of the membrane model is negligible, making it suitable for studying the internal dynamics of vesicles via coupling to hydrodynamic solvers or particle-based solvent models.

Effect of Mechanical Properties on Multicomponent Shell Patterning

Self-organized shells are fundamental in biological compartmentalization. They protect genomic material or enclose enzymes to aid the metabolic process. Studies of crystalline shells have shown the importance of the mechanical properties of building units in the shell morphology. However, the mechanism underlying the morphology of multicomponent assemblies is still poorly understood. Here, we analyze multicomponent closed shells that have different mechanical properties. By minimizing elastic energy, we show that heterogeneous bending rigidities regulate the surface pattern into circular, spikes, and ridge shapes. Interestingly, our continuum elasticity model recovers the patterns that have been proposed in bacterial microcompartments (BMCs), which are self-organized protein shells that aid the breakdown of complex molecules and allow bacteria to survive in hostile environments. In addition, our work elucidates the principles of pattern formation that can be used to design and engineer multicomponent microcompartments with a specific surface distribution of the components.

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.

Chemically controlled pattern formation in self-oscillating elastic shells

Patterns and morphology develop in living systems such as embryos in response to chemical signals. To understand and exploit the interplay of chemical reactions with mechanical transformations, chemomechanical polymer systems have been synthesized by attaching chemicals into hydrogels. In this work, we design autonomous responsive elastic shells that undergo morphological changes induced by chemical reactions. We couple the local mechanical response of the gel with the chemical processes on the shell. This causes swelling and deswelling of the gel, generating diverse morphological changes, including periodic oscillations. We further introduce a mechanical instability and observe buckling-unbuckling dynamics with a response time delay. Moreover, we investigate the mechanical feedback on the chemical reaction and demonstrate the dynamic patterns triggered by an initial deformation. We show the chemical characteristics that account for the shell morphology and discuss the future designs for autonomous responsive materials.

More than just a barrier: using physical models to couple membrane shape to cell function

The correct execution of many cellular processes, such as division and motility, requires the cell to adopt a specific shape. Physically, these shapes are determined by the interplay of the plasma membrane and internal cellular driving factors. While the plasma membrane defines the boundary of the cell, processes inside the cell can result in the generation of forces that deform the membrane. These processes include protein binding, the assembly of protein superstructures, and the growth and contraction of cytoskeletal networks. Due to the complexity of the cell, relating observed membrane deformations back to internal processes is a challenging problem. Here, we review cell shape changes in endocytosis, cell adhesion, cell migration and cell division and discuss how by modeling membrane deformations we can investigate the inner working principles of the cell.

Buckling transitions and soft-phase invasion of two-component icosahedral shells

What is the optimal distribution of two types of crystalline phases on the surface of icosahedral shells, such as of many viral capsids? We here investigate the distribution of a thin layer of soft material on a crystalline convex icosahedral shell. We demonstrate how the shapes of spherical viruses can be understood from the perspective of elasticity theory of thin two-component shells. We develop a theory of shape transformations of an icosahedral shell upon addition of a softer, but still crystalline, material onto its surface. We show how the soft component "invades" the regions with the highest elastic energy and stress imposed by the 12 topological defects on the surface. We explore the phase diagram as a function of the surface fraction of the soft material, the shell size, and the incommensurability of the elastic moduli of the rigid and soft phases. We find that, as expected, progressive filling of the rigid shell by the soft phase starts from the most deformed regions of the icosahedron. With a progressively increasing soft-phase coverage, the spherical segments of domes are filled first (12 vertices of the shell), then the cylindrical segments connecting the domes (30 edges) are invaded, and, ultimately, the 20 flat faces of the icosahedral shell tend to be occupied by the soft material. We present a detailed theoretical investigation of the first two stages of this invasion process and develop a model of morphological changes of the cone structure that permits noncircular cross sections. In conclusion, we discuss the biological relevance of some structures predicted from our calculations, in particular for the shape of viral capsids.

Competing pathways for the invagination of clathrin-coated membranes

Clathrin-mediated endocytosis is the major pathway by which eukaryotic cells take up extracellular material, but it is still elusive which physical pathways are being taken during membrane invagination. From a continuum point of view, it can be driven by increases in coat stiffness, preferred curvature or line tension. Here we develop a comprehensive theoretical framework that can be solved analytically and that predicts the consequences of these different scenarios. We find that for the case of increasing stiffness or preferred curvature, curvature will be acquired gradually with growth, while for increasing line tension, the lattice must have grown to a certain size before a flat-to-curved transition can occur. At low membrane tension, the critical value for coat stiffness is 30 kBT, for preferred curvature it is 200 nm, and for line tension it is 6 pN. For high membrane tension, critical coat stiffness is 150 kBT and critical preferred curvature is 70 nm. In the mixed case when a coat with finite rigidity but increasing line tension is considered, a cup-to-sphere transition can occur for a line tension of 6 pN. The flat-to-curved and the cup-to-sphere transitions driven by line tension are both suppressed by high membrane tension.

Role of combined cell membrane and wall mechanical properties regulated by polarity signals in cell budding

Budding yeast, Saccharomyces cerevisiae, serves as a prime biological model to study mechanisms underlying asymmetric growth. Previous studies have shown that prior to bud emergence, polarization of a conserved small GTPase Cdc42 must be established on the cell membrane of a budding yeast. Additionally, such polarization contributes to the delivery of cell wall remodeling enzymes and hydrolase from cytosol through the membrane, to change the mechanical properties of the cell wall. This leads to the hypothesis that Cdc42 and its associated proteins at least indirectly regulate cell surface mechanical properties. However, how the surface mechanical properties in the emerging bud are changed and whether such change is important are not well understood. To test several hypothesised mechanisms, a novel three-dimensional coarse-grained particle-based model has been developed which describes inhomogeneous mechanical properties of the cell surface. Model simulations predict alternation of the levels of stretching and bending stiffness of the cell surface in the bud region by the polarized Cdc42 signals is essential for initiating bud formation. Model simulations also suggest that bud shape depends strongly on the distribution of the polarized signaling molecules while the neck width of the emerging bud is strongly impacted by the mechanical properties of the chitin and septin rings. Moreover, the temporal change of the bud mechanical properties is shown to affect the symmetry of the bud shape. The 3D model of asymmetric cell growth can also be used for studying viral budding and other vegetative reproduction processes performed via budding, as well as detailed studies of cell growth.

Shapes of fluid membranes with chiral edges

We carry out Monte Carlo simulations of a colloidal fluid membrane with a free edge and composed of chiral rodlike viruses. The membrane is modeled by a triangular mesh of beads connected by bonds in which the bonds and beads are free to move at each Monte Carlo step. Since the constituent viruses are experimentally observed to twist only near the membrane edge, we use an effective energy that favors a particular sign of the geodesic torsion of the edge. The effective energy also includes the membrane bending stiffness, edge bending stiffness, and edge tension. We find three classes of membrane shapes resulting from the competition of the various terms in the free energy: branched shapes, chiral disks, and vesicles. Increasing the edge bending stiffness smooths the membrane edge, leading to correlations among the membrane normals at different points along the edge. The normalized power spectrum for edge displacements shows a peak with increasing preferred geodesic torsion. We also consider membrane shapes under an external force by fixing the distance between two ends of the membrane and finding the shape for increasing values of the distance between the two ends. As the distance increases, the membrane twists into a ribbon, with the force eventually reaching a plateau.

Shape selection and mis-assembly in viral capsid formation by elastic frustration

The successful assembly of a closed protein shell (or capsid) is a key step in the replication of viruses and in the production of artificial viral cages for bio/nanotechnological applications. During self-assembly, the favorable binding energy competes with the energetic cost of the growing edge and the elastic stresses generated due to the curvature of the capsid. As a result, incomplete structures such as open caps, cylindrical or ribbon-shaped shells may emerge, preventing the successful replication of viruses. Using elasticity theory and coarse-grained simulations, we analyze the conditions required for these processes to occur and their significance for empty virus self-assembly. We find that the outcome of the assembly can be recast into a universal phase diagram showing that viruses with high mechanical resistance cannot be self-assembled directly as spherical structures. The results of our study justify the need of a maturation step and suggest promising routes to hinder viral infections by inducing mis-assembly.

Mechanical conditions for stable symmetric cell constriction

Cell constriction is a decisive step for division in many cells. However, its physical pathway remains poorly understood, calling for a quantitative analysis of the forces required in different cytokinetic scenarios. Using a model cell composed by a flexible membrane (actin cortex and cell membrane) that encloses the cytoplasm, we study the mechanical conditions necessary for stable symmetric constriction under radial equatorial forces using analytical and numerical methods. We deduce that stable symmetric constriction requires positive effective spontaneous curvature, while spontaneous constriction requires a spontaneous curvature higher than the characteristic inverse cell size. Surface tension reduction (for example by actin cortex growth and membrane trafficking) increases the stability and spontaneity of cellular constriction. A reduction of external pressure also increases stability and spontaneity. Cells with prolate lobes (elongated cells) require lower stabilization forces than oblate-shaped cells (discocytes). We also show that the stability and spontaneity of symmetric constriction increase as constriction progresses. Our quantitative results settle the physical requirements for stable cytokinesis, defining a quantitative framework to analyze the mechanical role of the different constriction machinery and cytokinetic pathways found in real cells, so contributing to a deeper quantitative understanding of the physical mechanism of the cell division process.

Entry modes of ellipsoidal nanoparticles on a membrane during clathrin-mediated endocytosis

The membrane wrapping and internalization of nanoparticles, such as viruses and drug nanocarriers, through clathrin-mediated endocytosis (CME) are vitally important for intracellular transport. During CME, the shape of the particle plays crucial roles in the determination of particle-membrane interactions, but much of the previous work has been focused on spherical particles. In this work, we develop a stochastic model to study the CME of ellipsoidal particles. In our model, the deformation of the membrane and wrapping of the nanoparticles are driven by the accumulation of clathrin lattices, which is stimulated by the ligand-receptor interactions. Using our model, we systematically investigate the effect of particle shape (ellipsoids with different aspect ratios) on the CME. Our results show three entry modes: tip-first, tilted, and laying-down modes, used by ellipsoidal nanoparticles for internalization depending on the aspect ratio. Certain ellipsoids are able to take multiple entry modes for internalization. Interestingly, the prolate ellipsoid with an aspect ratio of 0.42 can be internalized with a significantly reduced number of ligand-receptor bonds. Particles which can be internalized with fewer bonds are excellent candidates for transcellular drug delivery. Moreover, our results demonstrate that internalization of ellipsoids with intermediate aspect ratios is easier than that of particles with low and high aspect ratios. Our model and simulations provide critical mechanistic insights into CME of ellipsoidal particles, and represent a viable platform for optimal design of nanoparticles for targeted drug delivery applications.

The equilibrium structure of self-assembled protein nano-cages

Understanding how highly symmetric, robust, monodisperse protein nano-cages self-assemble can have major applications in various areas of bio-nanotechnology, such as drug delivery, biomedical imaging and gene therapy. We develop a model to investigate the assembly of protein subunits into the structures with different sizes and symmetries. Using Monte Carlo simulation, we obtain global minimum energy structures. Our results suggest that the physical properties including the spontaneous curvature, flexibility and bending rigidity of coat proteins are sufficient to predict the size, symmetry and shape selectivity of the assembly products. Further, on a thermodynamic basis, we discuss the polymorphism of nano-cages observed in assembly experiments.

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.

Structural templating of J-aggregates: Visualizing bis(monoacylglycero)phosphate domains in live cells

Identifying the key structural and dynamical determinants that drive the association of biomolecules, whether in solution, or perhaps more importantly in a membrane environment, has critical implications for our understanding of cellular dynamics, processes, and signaling. With recent advances in high-resolution imaging techniques, from the development of new molecular labels to technical advances in imaging methodologies and platforms, researchers are now reaping the benefits of being able to directly characterize and quantify local dynamics, structures, and conformations in live cells and tissues. These capabilities are providing unique insights into association stoichiometries, interactions, and structures on sub-micron length scales. We previously examined the role of lipid headgroup chemistry and phase state in guiding the formation of pseudoisocyanine (PIC) dye J-aggregates on supported planar bilayers [Langmuir, 25, 10719]. We describe here how these same J-aggregates can report on the in situ formation of organellar membrane domains in live cells. Live cell hyperspectral confocal microscopy using GFP-conjugated GTPase markers of early (Rab5) and late (Rab7) endosomes revealed that the PIC J-aggregates were confined to domains on either the limiting membrane or intralumenal vesicles (ILV) of late endosomes, known to be enriched in the anionic lipid bis(monoacylglycero)phosphate (BMP). Correlated confocal fluorescence - atomic force microscopy performed on endosomal membrane-mimetic supported planar lipid bilayers confirmed BMP-specific templating of the PIC J-aggregates. These data provide strong evidence for the formation of BMP-rich lipid domains during multivesicular body formation and portend the application of structured dye aggregates as markers of cellular membrane domain structure, size, and formation.

Shape Morphologies of Icosahedral Two-Component Vesicles

What are the features of partitioning of crystalline materials on the surface of a two-component icosahedral vesicle? We model the response of the rigid hardly stretchable crystalline icosahedra upon addition of a softer component on its surface. We demonstrate how the soft phase "invades" the shell regions with the highest elastic energy density around 12 5-fold topological defects. We explore the phase diagram of these inhomogeneous shells as a function of the soft material fraction, shell radius, and elastic moduli of the two phases. The findings are compared with the recent computer simulation findings, and their biological relevance, for example, for the structure of icosahedral viruses, is also discussed.

Modeling the Mechanics of Cell Division: Influence of Spontaneous Membrane Curvature, Surface Tension, and Osmotic Pressure

Many cell division processes have been conserved throughout evolution and are being revealed by studies on model organisms such as bacteria, yeasts, and protozoa. Cellular membrane constriction is one of these processes, observed almost universally during cell division. It happens similarly in all organisms through a mechanical pathway synchronized with the sequence of cytokinetic events in the cell interior. Arguably, such a mechanical process is mastered by the coordinated action of a constriction machinery fueled by biochemical energy in conjunction with the passive mechanics of the cellular membrane. Independently of the details of the constriction engine, the membrane component responds against deformation by minimizing the elastic energy at every constriction state following a pathway still unknown. In this paper, we address a theoretical study of the mechanics of membrane constriction in a simplified model that describes a homogeneous membrane vesicle in the regime where mechanical work due to osmotic pressure, surface tension, and bending energy are comparable. We develop a general method to find approximate analytical expressions for the main descriptors of a symmetrically constricted vesicle. Analytical solutions are obtained by combining a perturbative expansion for small deformations with a variational approach that was previously demonstrated valid at the reference state of an initially spherical vesicle at isotonic conditions. The analytic approximate results are compared with the exact solution obtained from numerical computations, getting a good agreement for all the computed quantities (energy, area, volume, constriction force). We analyze the effects of the spontaneous curvature, the surface tension and the osmotic pressure in these quantities, focusing especially on the constriction force. The more favorable conditions for vesicle constriction are determined, obtaining that smaller constriction forces are required for positive spontaneous curvatures, low or negative membrane tension and hypertonic media. Conditions for spontaneous constriction at a given constriction force are also determined. The implications of these results for biological cell division are discussed. This work contributes to a better quantitative understanding of the mechanical pathway of cellular division, and could assist the design of artificial divisomes in vesicle-based self-actuated microsystems obtained from synthetic biology approaches.

Thermodynamic free energy methods to investigate shape transitions in bilayer membranes

The conformational free energy landscape of a system is a fundamental thermodynamic quantity of importance particularly in the study of soft matter and biological systems, in which the entropic contributions play a dominant role. While computational methods to delineate the free energy landscape are routinely used to analyze the relative stability of conformational states, to determine phase boundaries, and to compute ligand-receptor binding energies its use in problems involving the cell membrane is limited. Here, we present an overview of four different free energy methods to study morphological transitions in bilayer membranes, induced either by the action of curvature remodeling proteins or due to the application of external forces. Using a triangulated surface as a model for the cell membrane and using the framework of dynamical triangulation Monte Carlo, we have focused on the methods of Widom insertion, thermodynamic integration, Bennett acceptance scheme, and umbrella sampling and weighted histogram analysis. We have demonstrated how these methods can be employed in a variety of problems involving the cell membrane. Specifically, we have shown that the chemical potential, computed using Widom insertion, and the relative free energies, computed using thermodynamic integration and Bennett acceptance method, are excellent measures to study the transition from curvature sensing to curvature inducing behavior of membrane associated proteins. The umbrella sampling and WHAM analysis has been used to study the thermodynamics of tether formation in cell membranes and the quantitative predictions of the computational model are in excellent agreement with experimental measurements. Furthermore, we also present a method based on WHAM and thermodynamic integration to handle problems related to end-point-catastrophe that are common in most free energy methods.

On the modeling of endocytosis in yeast

The cell membrane deforms during endocytosis to surround extracellular material and draw it into the cell. Results of experiments on endocytosis in yeast show general agreement that 1) actin polymerizes into a network of filaments exerting active forces on the membrane to deform it, and 2) the large-scale membrane deformation is tubular in shape. In contrast, there are three competing proposals for precisely how the actin filament network organizes itself to drive the deformation. We use variational approaches and numerical simulations to address this competition by analyzing a meso-scale model of actin-mediated endocytosis in yeast. The meso-scale model breaks up the invagination process into three stages: 1) initiation, where clathrin interacts with the membrane via adaptor proteins; 2) elongation, where the membrane is then further deformed by polymerizing actin filaments; and 3) pinch-off. Our results suggest that the pinch-off mechanism may be assisted by a pearling-like instability. We rule out two of the three competing proposals for the organization of the actin filament network during the elongation stage. These two proposals could be important in the pinch-off stage, however, where additional actin polymerization helps break off the vesicle. Implications and comparisons with earlier modeling of endocytosis in yeast are discussed.

The effect of spontaneous curvature on a two-phase vesicle

Vesicles are membrane-bound structures commonly known for their roles in cellular transport and the shape of a vesicle is determined by its surrounding membrane (lipid bilayer). When the membrane is composed of different lipids, it is natural for the lipids of similar molecular structure to migrate towards one another (via spinodal decomposition), creating a multi-phase vesicle. In this article, we consider a two-phase vesicle model which is driven by nature's propensity to maintain a minimal state of elastic energy. The model assumes a continuum limit, thereby treating the membrane as a closed three-dimensional surface. The main purpose of this study is to reveal the complexity of the Helfrich two-phase vesicle model with non-zero spontaneous curvature and provide further evidence to support the relevance of spontaneous curvature as a modelling parameter. In this paper, we illustrate the complexity of the Helfrich two-phase model by providing multiple examples of undocumented solutions and energy hysteresis. We also investigate the influence of spontaneous curvature on morphological effects and membrane phenomena such as budding and fusion.

A balance between membrane elasticity and polymerization energy sets the shape of spherical clathrin coats

In endocytosis, scaffolding is one of the mechanisms to create membrane curvature by moulding the membrane into the spherical shape of the clathrin cage. However, the impact of membrane elastic parameters on the assembly and shape of clathrin lattices has never been experimentally evaluated. Here, we show that membrane tension opposes clathrin polymerization. We reconstitute clathrin budding in vitro with giant unilamellar vesicles (GUVs), purified adaptors and clathrin. By changing the osmotic conditions, we find that clathrin coats cause extensive budding of GUVs under low membrane tension while polymerizing into shallow pits under moderate tension. High tension fully inhibits polymerization. Theoretically, we predict the tension values for which transitions between different clathrin coat shapes occur. We measure the changes in membrane tension during clathrin polymerization, and use our theoretical framework to estimate the polymerization energy from these data. Our results show that membrane tension controls clathrin-mediated budding by varying the membrane budding energy.

A relationship between membrane tension and clathrin polymerization during endocytosis has not been experimentally established. Here, the authors show using an in vitro reconstituted system and theoretical modelling that membrane tension regulates clathrin polymerization into spherical cages by varying the membrane budding energy.

Mesoscale computational studies of membrane bilayer remodeling by curvature-inducing proteins

Biological membranes constitute boundaries of cells and cell organelles. These membranes are soft fluid interfaces whose thermodynamic states are dictated by bending moduli, induced curvature fields, and thermal fluctuations. Recently, there has been a flood of experimental evidence highlighting active roles for these structures in many cellular processes ranging from trafficking of cargo to cell motility. It is believed that the local membrane curvature, which is continuously altered due to its interactions with myriad proteins and other macromolecules attached to its surface, holds the key to the emergent functionality in these cellular processes. Mechanisms at the atomic scale are dictated by protein-lipid interaction strength, lipid composition, lipid distribution in the vicinity of the protein, shape and amino acid composition of the protein, and its amino acid contents. The specificity of molecular interactions together with the cooperativity of multiple proteins induce and stabilize complex membrane shapes at the mesoscale. These shapes span a wide spectrum ranging from the spherical plasma membrane to the complex cisternae of the Golgi apparatus. Mapping the relation between the protein-induced deformations at the molecular scale and the resulting mesoscale morphologies is key to bridging cellular experiments across the various length scales. In this review, we focus on the theoretical and computational methods used to understand the phenomenology underlying protein-driven membrane remodeling. Interactions at the molecular scale can be computationally probed by all atom and coarse grained molecular dynamics (MD, CGMD), as well as dissipative particle dynamics (DPD) simulations, which we only describe in passing. We choose to focus on several continuum approaches extending the Canham - Helfrich elastic energy model for membranes to include the effect of curvature-inducing proteins and explore the conformational phase space of such systems. In this description, the protein is expressed in the form of a spontaneous curvature field. The approaches include field theoretical methods limited to the small deformation regime, triangulated surfaces and particle-based computational models to investigate the large-deformation regimes observed in the natural state of many biological membranes. Applications of these methods to understand the properties of biological membranes in homogeneous and inhomogeneous environments of proteins, whose underlying curvature fields are either isotropic or anisotropic, are discussed. The diversity in the curvature fields elicits a rich variety of morphological states, including tubes, discs, branched tubes, and caveola. Mapping the thermodynamic stability of these states as a function of tuning parameters such as concentration and strength of curvature induction of the proteins is discussed. The relative stabilities of these self-organized shapes are examined through free-energy calculations. The suite of methods discussed here can be tailored to applications in specific cellular settings such as endocytosis during cargo trafficking and tubulation of filopodial structures in migrating cells, which makes these methods a powerful complement to experimental studies.

Budding behavior of multi-component vesicles

For multicomponent vesicles, the line tension of domain boundaries and the component-dependent elastic properties (e.g., spontaneous curvatures) are the two most important factors that mediate the budding behaviors of these vesicles. This paper especially focuses on their effects on the budding types and the budding number of a two-component vesicle. We found that the budding number is mainly determined by the component-dependent elastic properties while the budding types mediated by line tensions. A phase diagram is also obtained showing three different types of phase regions: (i) partially budding, (ii) partially bud-off, and (iii) totally bud-off regions. These inspiring results are, however, derived from a very simple spherical-cap model, and have been tested by computer simulations showing good agreement. We emphasize that, besides testing the spherical-cap model, the computer simulation techniques developed in the current work can be easily extended to other systems involving multicomponent vesicles.

Fluctuating shells under pressure

Thermal fluctuations strongly modify the large length-scale elastic behavior of cross-linked membranes, giving rise to scale-dependent elastic moduli. Whereas thermal effects in flat membranes are well understood, many natural and artificial microstructures are modeled as thin elastic shells. Shells are distinguished from flat membranes by their nonzero curvature, which provides a size-dependent coupling between the in-plane stretching modes and the out-of-plane undulations. In addition, a shell can support a pressure difference between its interior and its exterior. Little is known about the effect of thermal fluctuations on the elastic properties of shells. Here, we study the statistical mechanics of shape fluctuations in a pressurized spherical shell, using perturbation theory and Monte Carlo computer simulations, explicitly including the effects of curvature and an inward pressure. We predict novel properties of fluctuating thin shells under point indentations and pressure-induced deformations. The contribution due to thermal fluctuations increases with increasing ratio of shell radius to thickness and dominates the response when the product of this ratio and the thermal energy becomes large compared with the bending rigidity of the shell. Thermal effects are enhanced when a large uniform inward pressure acts on the shell and diverge as this pressure approaches the classical buckling transition of the shell. Our results are relevant for the elasticity and osmotic collapse of microcapsules.

Curvature-driven effective attraction in multicomponent membranes

We study closed liquid membranes that segregate into three phases due to differences in the chemical and physical properties of its components. The shape and in-plane membrane arrangement of the phases are coupled through phase-specific bending energies and line tensions. We use simulated annealing Monte Carlo simulations to find low-energy structures, allowing both phase arrangement and membrane shape to relax. The three-phase system is the simplest one in which there are multiple interface pairs, allowing us to analyze interfacial preferences and pairwise distinct line tensions. We observe the system's preference for interface pairs that maximize differences in spontaneous curvature. From a pattern selection perspective, this acts as an effective attraction between phases of most disparate spontaneous curvature. We show that this effective attraction is robust enough to persist even when the interface between these phases is the most penalized by line tension. This effect is driven by geometry and not by any explicit component-component interaction.

Minimal mesoscale model for protein-mediated vesiculation in clathrin-dependent endocytosis

In eukaryotic cells, the internalization of extracellular cargo via the endocytic machinery is an important regulatory process required for many essential cellular functions. The role of cooperative protein-protein and protein-membrane interactions in the ubiquitous endocytic pathway in mammalian cells, namely the clathrin-dependent endocytosis, remains unresolved. We employ the Helfrich membrane Hamiltonian together with surface evolution methodology to address how the shapes and energetics of vesicular-bud formation in a planar membrane are stabilized by presence of the clathrin-coat assembly. Our results identify a unique dual role for the tubulating protein epsin: multiple epsins localized spatially and orientationally collectively play the role of a curvature inducing capsid; in addition, epsin serves the role of an adapter in binding the clathrin coat to the membrane. Our results also suggest an important role for the clathrin lattice, namely in the spatial- and orientational-templating of epsins. We suggest that there exists a critical size of the coat above which a vesicular bud with a constricted neck resembling a mature vesicle is stabilized. Based on the observed strong dependence of the vesicle diameter on the bending rigidity, we suggest that the variability in bending stiffness due to variations in membrane composition with cell type can explain the experimentally observed variability on the size of clathrin-coated vesicles, which typically range 50-100 nm. Our model also provides estimates for the number of epsins involved in stabilizing a coated vesicle, and without any direct fitting reproduces the experimentally observed shapes of vesicular intermediates as well as their probability distributions quantitatively, in wildtype as well as CLAP IgG injected neuronal cell experiments. We have presented a minimal mesoscale model which quantitatively explains several experimental observations on the process of vesicle nucleation induced by the clathrin-coated assembly prior to vesicle scission in clathrin dependent endocytosis.

Domain-driven morphogenesis of cellular membranes

Cellular membrane systems delimit and organize the intracellular space. Most of the morphological rearrangements in cells involve the coordinated remodeling of the lipid bilayer, the core of the membranes. This process is generally thought to be initiated and coordinated by specialized protein machineries. Nevertheless, it has become increasingly evident that the most essential part of the geometric information and energy required for membrane remodeling is supplied via the cooperative and synergistic action of proteins and lipids, as cellular shapes are constructed using the intrinsic dynamics, plasticity and self-organizing capabilities provided by the lipid bilayer. Here, we analyze the essential role of proteo-lipid membrane domains in conducting and coordinating morphological remodeling in cells.

Budding and vesiculation induced by conical membrane inclusions

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

Vesicle formation by self-assembly of membrane-bound matrix proteins into a fluidlike budding domain

The shape of enveloped viruses depends critically on an internal protein matrix, yet it remains unclear how the matrix proteins control the geometry of the envelope membrane. We found that matrix proteins purified from Newcastle disease virus adsorb on a phospholipid bilayer and condense into fluidlike domains that cause membrane deformation and budding of spherical vesicles, as seen by fluorescent and electron microscopy. Measurements of the electrical admittance of the membrane resolved the gradual growth and rapid closure of a bud followed by its separation to form a free vesicle. The vesicle size distribution, confined by intrinsic curvature of budding domains, but broadened by their merger, matched the virus size distribution. Thus, matrix proteins implement domain-driven mechanism of budding, which suffices to control the shape of these proteolipid vesicles.

Monte Carlo study of crystalline order and defects on weakly curved surfaces

We numerically study the ground states of particles interacting via a repulsive Yukawa potential on two rigid substrates shaped as isolated and periodically arranged bumps characterized by a spatially varying Gaussian curvature. Below a critical aspect ratio that describes the substrate deformation, the lattice is frustrated, but defect free. A further increase of the aspect ratio triggers defect unbinding transitions that lower the total potential energy by introducing dislocations either in isolation or within grain boundaries. In the presence of very strong deformations, isolated disclinations are nucleated. We show that the character and spatial distribution of defects observed in the ground state reflect the symmetries and periodicity of the two model surfaces investigated in this study.

Effect of line tension on the lateral organization of lipid membranes

The principles of organization and functioning of cellular membranes are currently not well understood. The raft hypothesis suggests the existence of domains or rafts in cell membranes, which behave as protein and lipid platforms. They have a functional role in important cellular processes, like protein sorting or cell signaling, among others. Theoretical work suggests that the interfacial energy at the domain edge, also known as line tension, is a key parameter determining the distribution of domain sizes, but there is little evidence of how line tension affects membrane organization. We have investigated the effects of the line tension on the formation and stability of liquid ordered domains in model lipid bilayers with raft-like composition by means of time-lapse confocal microscopy coupled to atomic force microscopy. We varied the hydrophobic mismatch between the two phases, and consequently the line tension, by modifying the thickness of the disordered phase with phosphatidylcholines of different acyl chain length. The temperature of domain formation, the dynamics of domain growth, and the distribution of domain sizes depend strongly on the thickness difference between the domains and the surrounding membrane, which is related to line tension. When considering line tension calculated from a theoretical model, our results revealed a linear increase of the temperature of domain formation and domain growth rate with line tension. Domain budding was also shown to depend on height mismatch. Our experiments contribute significantly to our knowledge of the physical-chemical parameters that control membrane organization. Importantly, the general trends observed can be extended to cellular membranes.

Domain-induced budding in buckling membranes

We present a phase field model on buckling membranes to analyze phase separation and budding on soft membranes. By numerically integrating dynamic equations, it turns out that the formation of caps is greatly influenced by the presence of a little excess area due to the surface area constraint. When cap-shaped domains are created, domain coalescence is mainly observed not between domains with same budding directions, but between domains with opposite budding directions, because the bending energy between two domains is larger in the former case. Although we do not introduce spontaneous curvature like Helfrich model, we obtain some suggestions related to the slow dynamics of the phase separation on vesicles.

Defect scars on flexible surfaces with crystalline order

The crystallography of two-dimensional particle packings on flexible surfaces of spherical topology is investigated. Examples are viral capsids and crystalline vesicles. Computer simulations of dynamically triangulated surfaces are employed to study the shape and structure of lattice defects as a function of the Föppl-von Kármán number gamma. We find that grain-boundary scars become much more fuzzy with increasing temperature, that the size of grain-boundary scars saturates with increasing vesicle radius, and that the buckling transition shifts to higher values of gamma due to the presence of scars.

Shape transitions in lipid membranes and protein mediated vesicle fusion and fission

In the cell, the plasma membrane is often densely decorated by transmembrane proteins. The morphology and dynamics of the membrane are strongly influenced by the presence of proteins. In this paper, we use a coarse-grained model to explore the composite membrane-protein system and develop a simulation methodology based on thermodynamic integration to examine free energy changes during membrane shape transitions. The authors show that a critical concentration of conical membrane proteins or proteins with nonzero spontaneous curvature can drive the formation of small vesicles. The driving force of vesicle budding stems from the preference of proteins to gather in regions of high curvature. A sufficiently high concentration of proteins therefore can influence the topology of the membrane. The biological significance of our results is discussed.

Self-consistent field theory simulations of block copolymer assembly on a sphere

Recently there has been a strong interest in the area of defect formation in ordered structures on curved surfaces. Here we explore the closely related topic of self-assembly in thin block copolymer melt films confined to the surface of a sphere. Our study is based on a self-consistent field theory (SCFT) model of block copolymers that is numerically simulated by spectral collocation with a spherical harmonic basis and an extension of the Rasmussen-Kalosakas operator splitting algorithm [J. Polym. Sci. Part B: Polym. Phys. 40, 1777 (2002)]. In this model, we assume that the composition of the thin block copolymer film varies only in longitude and colatitude and is constant in the radial direction. Using this approach we are able to study the formation of defects in the lamellar and cylindrical phases, and their dependence on sphere radius. Specifically, we compute ground-state (i.e., lowest-energy) configurations on the sphere for both the cylindrical and lamellar phases. Grain boundary scars are also observed in our simulations of the cylindrical phase when the sphere radius surpasses a threshold value R_{c} approximately 5d , where d is the natural lattice spacing of the cylindrical phase, which is consistent with theoretical predictions [Bowick, Phys. Rev. B 62, 8738 (2000); Bausch, Science 299, 1716 (2003)]. A strong segregation limit approximate free energy is also presented, along with simple microdomain packing arguments, to shed light on the observed SCFT simulation results.

Crystallography on curved surfaces

We study static and dynamical properties that distinguish 2D crystals constrained to lie on a curved substrate from their flat-space counterparts. A generic mechanism of dislocation unbinding in the presence of varying Gaussian curvature is presented in the context of a model surface amenable to full analytical treatment. We find that glide diffusion of isolated dislocations is suppressed by a binding potential of purely geometrical origin. Finally, the energetics and biased diffusion dynamics of point defects such as vacancies and interstitials are explained in terms of their geometric potential.

Budding and fission of a multiphase vesicle

We present a model of bi-phasic vesicles in the limit of large surface tension. In this regime, the vesicle is completely stretched and well described by two spherical caps with a fold, which concentrates the membrane stress. The conservation laws and geometric constraints restrict the space of possible shapes to a pair of solutions labeled by a parameter sigma given by line tension/pressure. For a given value of sigma, the two solutions differ by the length of the interface between domains. For a critical value, sigmac, the two vesicle shapes become identical and no connected solution exists above this critical value. This model sheds new light on two proposed mechanisms (osmotic shocks and molecule absorption) to explain the budding and the fission in recent experiments.

Self-assembled pattern formation of block copolymers on the surface of the sphere using self-consistent field theory

The spherical surface is spatially discretized with triangular lattices to numerically calculate the Laplace-Beltrami operator contained in the self-consistent field theory (SCFT) equations using a finite volume method. Based on this method we have developed a spherical alternating-direction implicit (ADI) scheme for the first time to help extend real-space implementation of SCFT in 2D flat space to the surface of the sphere. By using this method, we simulate the equilibrium microphase separation morphology of block copolymers including AB diblocks, ABC linear triblocks and ABC star triblock copolymers occurred on the spherical surface. In general, two classes of microphase separation morphologies such as striped patterns for compositionally symmetric block copolymers and spotted patterns for asymmetric compositions have been found. In contrast to microphase separation morphology in 2D flat space, the geometrical characteristics of a sphere has a large influence on the self-assembled morphology. For striped patterns, several of spiral-form and ring-form patterns are found by changing the ratio of the radius of a sphere to the averaging width of the stripes. The specific pattern such as the striped and spotted pattern with intrinsic dislocations or defects stems from formed periodic patterns due to microphase separation of block copolymers arranged on the curved surface.

Pressure tuning of the morphology of heterogeneous lipid vesicles: a two-photon-excitation fluorescence microscopy study

We used a technique that allows us to visualize local and morphological changes of the membrane of more component giant unilamellar vesicles due to high pressure perturbation. Under these conditions, thermally induced processes are largely suppressed, and the bending rigidity and line tension are influenced by pressure-induced changes in lipid molecular packing and shape only. We studied the effect of pressure on the lateral organization and morphology of the model raft system DOPC (1,2-dioleoyl-sn-glycero-3-phosphocholine)/sphingomyelin/cholesterol as well as of the fluid mixture POPC (1-palmitoyl-2-oleoyl-sn-glycero-3-phosphocholine)/DLPC (1,2-dilauroyl-sn-glycero-3-phosphocholine) by two-photon excitation fluorescence microscopy. The pressure-dependent experiments were carried out using a sample cell made from a thin fused silica capillary. The use of Laurdan as fluorescence label allowed us to also follow the lipid phase state by calculating the generalized polarization (GP) values of the vesicles and extracting their average value. During the compression cycle, a reduction in the volume of the vesicles is observed, accompanied by an increase of the average GP value, indicating an increasingly tighter packing of the lipids. Interestingly, the two systems studied show phenomena of budding and fission, and these at surprisingly low pressures of 200-300 bar. Moreover, these budding processes are not directly related to phase transitions to an overall ordered conformational state of the lipid membrane, which occur at much higher pressures. The topological changes of the lipid vesicles are irreversible and exhibit a different behavior depending on whether the pressure is increased or decreased. The results are discussed in light of the various contributions to the free energy functional of lipid vesicles. Finally, the biological relevance of these studies is highlighted.

Measuring the elasticity of clathrin-coated vesicles via atomic force microscopy

Using a new scheme based on atomic force microscopy (AFM), we investigate mechanical properties of clathrin-coated vesicles (CCVs). CCVs are multicomponent protein and lipid complexes of approximately 100 nm diameter that are implicated in many essential cell-trafficking processes. Our AFM imaging resolves clathrin lattice polygons and provides height deformation in quantitative response to AFM-substrate compression force. We model CCVs as multilayered elastic spherical shells and, from AFM measurements, estimate their bending rigidity to be 285 +/- 30 k(B)T, i.e., approximately 20 times that of either the outer clathrin cage or inner vesicle membrane. Further analysis reveals a flexible coupling between the clathrin coat and the membrane, a structural property whose modulation may affect vesicle biogenesis and cellular function.

Implicit solvent simulation models for biomembranes

Fully atomic simulation strategies are infeasible for the study of many processes of interest to membrane biology, biophysics and biochemistry. We review various coarse-grained simulation methodologies with special emphasis on methods and models that do not require the explicit simulation of water. Examples from our own research demonstrate that such models have potential for simulating a variety of biologically relevant phenomena at the membrane surface.

Phase separation patterns for diblock copolymers on spherical surfaces: a finite volume method

We explore phase separation on spherical surfaces by solving the Cahn-Hilliard equation modified for diblock copolymers using a finite volume method. The spherical surface is discretized into almost uniform triangles by employing successive dyadic refinements of the spherical icosahedron, a methodology that avoids potential mathematical and numerical problems related to the poles in spherical coordinates. The finite volume method is based on averaging Voronoi cells built from triangular meshes to calculate the Laplace-Beltrami operator on the curved surface, which greatly improves both the accuracy and speed of calculation as compared to the conventional finite difference method. By using this method we simulate the phase separation of diblock copolymers on a spherical surface. It is found that stable and intrinsic defects, which would not occur in a flat space after sufficient annealing, appear in the periodic arrangement of the domains on the curved surface due to the distinct Euler characteristic of the surface.