Shape dynamics, lipid hydrodynamics, and the complex viscoelasticity of bilayer membranes [corrected]

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.

A Review of Continuum Mechanics for Mechanical Deformation of Lipid Membranes

Mechanical deformation of lipid membranes plays important roles in various cellular tasks. Curvature deformation and lateral stretching are two major energy contributions to the mechanical deformation of lipid membranes. In this paper, continuum theories for these two major membrane deformation events were reviewed. Theories based on curvature elasticity and lateral surface tension were introduced. Numerical methods as well as biological applications of the theories were discussed.

A Lagrangian Thin-Shell Finite Element Method for Interacting Particles on Fluid Membranes

A recurring motif in soft matter and biophysics is modeling the mechanics of interacting particles on fluid membranes. One of the main outstanding challenges in these applications is the need to model the strong coupling between the substrate deformation and the particles' positions as the latter freely move on the former. This work presents a thin-shell finite element formulation based on subdivision surfaces to compute equilibrium configurations of a thin fluid shell with embedded particles. We use a variational Lagrangian framework to couple the mechanics of the particles and the substrate without having to resort to ad hoc constraints to anchor the particles to the surface. Unlike established methods for such systems, the particles are allowed to move between elements of the finite element mesh. This is achieved by parametrizing the particle locations on the reference configuration. Using the Helfrich-Canham energy as a model for fluid shells, we present the finite element method's implementation and an efficient search algorithm required to locate particles on the reference mesh. Several analyses with varying numbers of particles are finally presented reproducing symmetries observed in the classic Thomson problem and showcasing the coupling between interacting particles and deformable membranes.

The many faces of membrane tension: Challenges across systems and scales

Our understanding of the role of membrane tension in the field of membrane biophysics is rapidly evolving from a passive construct to an active player in a variety of cellular phenomena. Membrane tension has been shown to be a key regulator of many cellular processes ranging including trafficking, ion channel activation, and the invasion of red blood cells by malaria parasites. Recent experimental advances in cells, including the development of a fluorescent tension reporter, have shown that membrane tension is heterogeneous. In this mini-review, I summarize the recent advances in membrane tension measurements and discuss the contributions from different cellular constituents such as the cortical cytoskeleton. Then, I will explore how these different complexities can be considered in biophysical models of different scales. Finally, I will elaborate on the need for iterations between models and experiments as technologies in both fields advance to enable us to obtain critical insights into the physiological role of membrane tension as a critical component of mechanotransduction.

Modeling Receptor Motility along Advecting Lipid Membranes

This work aims to overview multiphysics mechanobiological computational models for receptor dynamics along advecting cell membranes. Continuum and statistical models of receptor motility are the two main modeling methodologies identified in reviewing the state of the art. Within the former modeling class, a further subdivision based on different biological purposes and processes of proteins’ motion is recognized; cell adhesion, cell contractility, endocytosis, and receptor relocations on advecting membranes are the most relevant biological processes identified in which receptor motility is pivotal. Numerical and/or experimental methods and approaches are highlighted in the exposure of the reviewed works provided by the literature, pertinent to the topic of the present manuscript. With a main focus on the continuum models of receptor motility, we discuss appropriate multiphyisics laws to model the mass flux of receptor proteins in the reproduction of receptor relocation and recruitment along cell membranes to describe receptor–ligand chemical interactions, and the cell’s structural response. The mass flux of receptor modeling is further supported by a discussion on the methodology utilized to evaluate the protein diffusion coefficient developed over the years.

Absolute vs Convective Instabilities and Front Propagation in Lipid Membrane Tubes

We analyze the stability of biological membrane tubes, with and without a base flow of lipids. Membrane dynamics are completely specified by two dimensionless numbers: the well-known Föppl-von Kármán number Γ and the recently introduced Scriven-Love number SL, respectively quantifying the base tension and base flow speed. For unstable tubes, the growth rate of a local perturbation depends only on Γ, whereas SL governs the absolute versus convective nature of the instability. Furthermore, nonlinear simulations of unstable tubes reveal an initially localized disturbance result in propagating fronts, which leave a thin atrophied tube in their wake. Depending on the value of Γ, the thin tube is connected to the unperturbed regions via oscillatory or monotonic shape transitions-reminiscent of recent experimental observations on the retraction and atrophy of axons. We elucidate our findings through a weakly nonlinear analysis, which shows membrane dynamics may be approximated by a model of the class of extended Fisher-Kolmogorov equations. Our study sheds light on the pattern selection mechanism in axonal shapes by recognizing the existence of two Lifshitz points, at which the front dynamics undergo steady-to-oscillatory bifurcations.

Subdiffusive-Brownian crossover in membrane proteins: a generalized Langevin equation-based approach

In this work, we propose a generalized Langevin equation-based model to describe the lateral diffusion of a protein in a lipid bilayer. The memory kernel is represented in terms of a viscous (instantaneous) and an elastic (noninstantaneous) component modeled through a Dirac δ function and a three-parameter Mittag-Leffler type function, respectively. By imposing a specific relationship between the parameters of the three-parameter Mittag-Leffler function, the different dynamical regimes-namely ballistic, subdiffusive, and Brownian, as well as the crossover from one regime to another-are retrieved. Within this approach, the transition time from the ballistic to the subdiffusive regime and the spectrum of relaxation times underlying the transition from the subdiffusive to the Brownian regime are given. The reliability of the model is tested by comparing the mean-square displacement derived in the framework of this model and the mean-square displacement of a protein diffusing in a membrane calculated through molecular dynamics simulations.

Biomembranes undergo complex, non-axisymmetric deformations governed by Kirchhoff-Love kinematicsand revealed by a three-dimensional computational framework

Biomembranes play a central role in various phenomena like locomotion of cells, cell-cell interactions, packaging and transport of nutrients, transmission of nerve impulses, and in maintaining organelle morphology and functionality. During these processes, the membranes undergo significant morphological changes through deformation, scission, and fusion. Modelling the underlying mechanics of such morphological changes has traditionally relied on reduced order axisymmetric representations of membrane geometry and deformation. Axisymmetric representations, while robust and extensively deployed, suffer from their inability to model-symmetry breaking deformations and structural bifurcations. To address this limitation, a three-dimensional computational mechanics framework for high fidelity modelling of biomembrane deformation is presented. The proposed framework brings together Kirchhoff–Love thin-shell kinematics, Helfrich-energy-based mechanics, and state-of-the-art numerical techniques for modelling deformation of surface geometries. Lipid bilayers are represented as spline-based surface discretizations immersed in a three-dimensional space; this enables modelling of a wide spectrum of membrane geometries, boundary conditions, and deformations that are physically admissible in a three-dimensional space. The mathematical basis of the framework and its numerical machinery are presented, and their utility is demonstrated by modelling three classical, yet non-trivial, membrane deformation problems: formation of tubular shapes and their lateral constriction, Piezo1-induced membrane footprint generation and gating response, and the budding of membranes by protein coats during endocytosis. For each problem, the full three-dimensional membrane deformation is captured, potential symmetry-breaking deformation paths identified, and various case studies of boundary and load conditions are presented. Using the endocytic vesicle budding as a case study, we also present a ‘phase diagram’ for its symmetric and broken-symmetry states.

Mitochondrial cristae modeled as an out-of-equilibrium membrane driven by a proton field

As the places where most of the fuel of the cell, namely, ATP, is synthesized, mitochondria are crucial organelles in eukaryotic cells. The shape of the invaginations of the mitochondria inner membrane, known as a crista, has been identified as a signature of the energetic state of the organelle. However, the interplay between the rate of ATP synthesis and the crista shape remains unclear. In this work, we investigate the crista membrane deformations using a pH-dependent Helfrich model, maintained out of equilibrium by a diffusive flux of protons. This model gives rise to shape changes of a cylindrical invagination, in particular to the formation of necks between wider zones under variable, and especially oscillating, proton flux.

A physical mechanism of TANGO1-mediated bulky cargo export

The endoplasmic reticulum (ER)-resident protein TANGO1 assembles into a ring around ER exit sites (ERES), and links procollagens in the ER lumen to COPII machinery, tethers, and ER-Golgi intermediate compartment (ERGIC) in the cytoplasm (Raote et al., 2018). Here, we present a theoretical approach to investigate the physical mechanisms of TANGO1 ring assembly and how COPII polymerization, membrane tension, and force facilitate the formation of a transport intermediate for procollagen export. Our results indicate that a TANGO1 ring, by acting as a linactant, stabilizes the open neck of a nascent COPII bud. Elongation of such a bud into a transport intermediate commensurate with bulky procollagens is then facilitated by two complementary mechanisms: (i) by relieving membrane tension, possibly by TANGO1-mediated fusion of retrograde ERGIC membranes and (ii) by force application. Altogether, our theoretical approach identifies key biophysical events in TANGO1-driven procollagen export.

Nonaxisymmetric Shapes of Biological Membranes from Locally Induced Curvature

In various biological processes such as endocytosis and caveolae formation, the cell membrane is locally deformed into curved morphologies. Previous models to study membrane morphologies resulting from locally induced curvature often only consider the possibility of axisymmetric shapes-an indeed unphysical constraint. Past studies predict that the cell membrane buds at low resting tensions and stalls at a flat pit at high resting tensions. In this work, we lift the restriction to axisymmetry to study all possible membrane morphologies. Only if the resting tension of the membrane is low, we reproduce axisymmetric membrane morphologies. When the resting tension is moderate to high, we show that 1) axisymmetric membrane pits are unstable and 2) nonaxisymmetric ridge-shaped structures are energetically favorable. Furthermore, we find the interplay between intramembrane viscous flow and the rate of induced curvature affects the membrane's ability to transition into nonaxisymmetric ridges and axisymmetric buds. In particular, we show that axisymmetric buds are favored when the induced curvature is rapidly increased, whereas nonaxisymmetric ridges are favored when the curvature is slowly increased. Our results hold relevant implications for biological processes such as endocytosis and physical phenomena like phase separation in lipid bilayers.

A mechanical model reveals that non-axisymmetric buckling lowers the energy barrier associated with membrane neck constriction

Membrane neck formation is essential for scission, which, as recent experiments on tubules have demonstrated, can be location dependent. The diversity of biological machinery that can constrict a neck such as dynamin, actin, ESCRTs and BAR proteins, and the range of forces and deflection over which they operate, suggest that the constriction process is functionally mechanical and robust to changes in biological environment. In this study, we used a mechanical model of the lipid bilayer to systematically investigate the influence of location, symmetry constraints, and helical forces on membrane neck constriction. Simulations from our model demonstrated that the energy barriers associated with constriction of a membrane neck are location-dependent. Importantly, if symmetry restrictions are relaxed, then the energy barrier for constriction is dramatically lowered and the membrane buckles at lower values of forcing parameters. Our simulations also show that constriction due to helical proteins further reduces the energy barrier for neck formation when compared to cylindrical proteins. These studies establish that despite different molecular mechanisms of neck formation in cells, the mechanics of constriction naturally leads to a loss of symmetry that can lower the energy barrier to constriction.

Geometric coupling of helicoidal ramps and curvature-inducing proteins in organelle membranes

Cellular membranes display an incredibly diverse range of shapes, both in the plasma membrane and at membrane bound organelles. These morphologies are intricately related to cellular functions, enabling and regulating fundamental membrane processes. However, the biophysical mechanisms at the origin of these complex geometries are not fully understood from the standpoint of membrane-protein coupling. In this study, we focused on a minimal model of helicoidal ramps representative of specialized endoplasmic reticulum compartments. Given a helicoidal membrane geometry, we asked what is the distribution of spontaneous curvature required to maintain this shape at mechanical equilibrium? Based on the Helfrich energy of elastic membranes with spontaneous curvature, we derived the shape equation for minimal surfaces, and applied it to helicoids. We showed the existence of switches in the sign of the spontaneous curvature associated with geometric variations of the membrane structures. Furthermore, for a prescribed gradient of spontaneous curvature along the exterior boundaries, we identified configurations of the helicoidal ramps that are confined between two infinitely large energy barriers. Overall our results suggest possible mechanisms for geometric control of helicoidal ramps in membrane organelles based on curvature-inducing proteins.

Interplay of elastic instabilities and viscoelasticity in the finite deformation of thin membranes

Pneumatic structures and actuators are found in a variety of natural and engineered systems such as dielectric actuators, soft robots, plants and fungi cells, or even the vocal sac of frogs. These structures are often subjected to mechanical instabilities arising from the thinning of their cross section and that may be harvested to perform mechanical work at a low energetic cost. While most of our understanding of this unstable behavior is for purely elastic membranes, real materials including lipid bilayers, elastomers, and connective tissues typically display a time-dependent viscoelastic response. This paper thus explores the role of viscous effects on the nature of this elastic instability when such membranes are dynamically inflated. For this, we first introduce an extension of the transient network theory to describe the finite strain viscoelastic response of membranes, enabling an elegant formulation while keeping a close connection with the dynamics of the underlying polymer network. We then combine experiments and simulations to analyze the viscoelastic behavior of an inflated blister made of a commercial adhesive tape (VHB 4905). Our results show that the viscous component induces a rich spectrum of behaviors bounded by two well-known elastic solutions corresponding to very high and very low inflation rates. We also show that membrane relaxation may induce unwanted buckling when it is subjected to cyclic inflations at certain frequencies. These results have clear implications for the inflation and mechanical work performed by time-dependent pneumatic structures and instability-based actuators.

pH sensing by lipids in membranes: The fundamentals of pH-driven migration, polarization and deformations of lipid bilayer assemblies

Most biological molecules contain acido-basic groups that modulate their structure and interactions. A consequence is that pH gradients, local heterogeneities and dynamic variations are used by cells and organisms to drive or regulate specific biological functions including energetic metabolism, vesicular traffic, migration and spatial patterning of tissues in development. While the direct or regulatory role of pH in protein function is well documented, the role of hydrogen and hydroxyl ions in modulating the properties of lipid assemblies such as bilayer membranes is only beginning to be understood. Here, we review approaches using artificial lipid vesicles that have been instrumental in providing an understanding of the influence of pH gradients and local variations on membrane vectorial motional processes: migration, membrane curvature effects promoting global or local deformations, crowding generation by segregative polarization processes. In the case of pH induced local deformations, an extensive theoretical framework is given and an application to a specific biological issue, namely the structure and stability of mitochondrial cristae, is described. This article is part of a Special Issue entitled: Emergence of Complex Behavior in Biomembranes edited by Marjorie Longo.

Size-dependent formation of membrane nanotubes: continuum modeling and molecular dynamics simulations

With continuum theory and molecular dynamics simulations we demonstrated that the lipid membrane upon extraction exhibits size- and tension-dependent mechanical behaviors, and different structural lipid rearrangements in different leaflets.

Multi-physics interactions drive VEGFR2 relocation on endothelial cells

Vascular Endothelial Growth Factor Receptor-2 (VEGFR2) is a pro-angiogenic receptor, expressed on endothelial cells (ECs). Although biochemical pathways that follow the VEGFR2 activation are well established, knowledge about the dynamics of receptors on the plasma membrane remains limited. Ligand stimulation induces the polarization of ECs and the relocation of VEGFR2, either in cell protrusions or in the basal aspect in cells plated on ligand-enriched extracellular matrix (ECM). We develop a mathematical model in order to simulate the relocation of VEGFR2 on the cell membrane during the mechanical adhesion of cells onto a ligand-enriched substrate. Co-designing the in vitro experiments with the simulations allows identifying three phases of the receptor dynamics, which are controlled respectively by the high chemical reaction rate, by the mechanical deformation rate, and by the diffusion of free receptors on the membrane. The identification of the laws that regulate receptor polarization opens new perspectives toward developing innovative anti-angiogenic strategies through the modulation of EC activation.

Irreversible thermodynamics of curved lipid membranes

The theory of irreversible thermodynamics for arbitrarily curved lipid membranes is presented here. The coupling between elastic bending and irreversible processes such as intramembrane lipid flow, intramembrane phase transitions, and protein binding and diffusion is studied. The forms of the entropy production for the irreversible processes are obtained, and the corresponding thermodynamic forces and fluxes are identified. Employing the linear irreversible thermodynamic framework, the governing equations of motion along with appropriate boundary conditions are provided.

Probing Elastic and Viscous Properties of Phospholipid Bilayers Using Neutron Spin Echo Spectroscopy

The elastic and viscous properties of self-assembled amphiphilic membranes dictate the intricate hierarchy of their structure and dynamics ranging from the diffusion of individual molecules to the large-scale deformation of the membrane. We previously demonstrated that neutron spin echo spectroscopy measurements of model amphiphilic membranes can access the naturally occurring submicrosecond membrane motions, such as bending and thickness fluctuations. Here we show how the experimentally measured fluctuation parameters can be used to determine the inherent membrane properties and demonstrate how membrane viscosity and compressibility modulus are influenced by lipid composition in a series of simple phosphatidylcholine bilayers with different tail lengths as a function of temperature. This approach highlights the interdependence of the bilayer elastic and viscous properties and the collective membrane dynamics and opens new avenues to investigating the mechanical properties of more complex and biologically inspired systems.

Systems biology of cellular membranes: a convergence with biophysics

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

Mechanics of active surfaces

We derive a fully covariant theory of the mechanics of active surfaces. This theory provides a framework for the study of active biological or chemical processes at surfaces, such as the cell cortex, the mechanics of epithelial tissues, or reconstituted active systems on surfaces. We introduce forces and torques acting on a surface, and derive the associated force balance conditions. We show that surfaces with in-plane rotational symmetry can have broken up-down, chiral, or planar-chiral symmetry. We discuss the rate of entropy production in the surface and write linear constitutive relations that satisfy the Onsager relations. We show that the bending modulus, the spontaneous curvature, and the surface tension of a passive surface are renormalized by active terms. Finally, we identify active terms which are not found in a passive theory and discuss examples of shape instabilities that are related to active processes in the surface.

Theory and algorithms to compute Helfrich bending forces: a review

Cell membranes are vital to shield a cell's interior from the environment. At the same time they determine to a large extent the cell's mechanical resistance to external forces. In recent years there has been considerable interest in the accurate computational modeling of such membranes, driven mainly by the amazing variety of shapes that red blood cells and model systems such as vesicles can assume in external flows. Given that the typical height of a membrane is only a few nanometers while the surface of the cell extends over many micrometers, physical modeling approaches mostly consider the interface as a two-dimensional elastic continuum. Here we review recent modeling efforts focusing on one of the computationally most intricate components, namely the membrane's bending resistance. We start with a short background on the most widely used bending model due to Helfrich. While the Helfrich bending energy by itself is an extremely simple model equation, the computation of the resulting forces is far from trivial. At the heart of these difficulties lies the fact that the forces involve second order derivatives of the local surface curvature which by itself is the second derivative of the membrane geometry. We systematically derive and compare the different routes to obtain bending forces from the Helfrich energy, namely the variational approach and the thin-shell theory. While both routes lead to mathematically identical expressions, so-called linear bending models are shown to reproduce only the leading order term while higher orders differ. The main part of the review contains a description of various computational strategies which we classify into three categories: the force, the strong and the weak formulation. We finally give some examples for the application of these strategies in actual simulations.

Analysis of diffusion in curved surfaces and its application to tubular membranes

Diffusion of particles in curved surfaces is inherently complex compared with diffusion in a flat membrane, owing to the nonplanarity of the surface. The consequence of such nonplanar geometry on diffusion is poorly understood but is highly relevant in the case of cell membranes, which often adopt complex geometries. To address this question, we developed a new finite element approach to model diffusion on curved membrane surfaces based on solutions to Fick's law of diffusion and used this to study the effects of geometry on the entry of surface-bound particles into tubules by diffusion. We show that variations in tubule radius and length can distinctly alter diffusion gradients in tubules over biologically relevant timescales. In addition, we show that tubular structures tend to retain concentration gradients for a longer time compared with a comparable flat surface. These findings indicate that sorting of particles along the surfaces of tubules can arise simply as a geometric consequence of the curvature without any specific contribution from the membrane environment. Our studies provide a framework for modeling diffusion in curved surfaces and suggest that biological regulation can emerge purely from membrane geometry.

Two-phase vesicles: a study on evolutionary and stationary models

In the current article, the dynamic evolution of two-phase vesicles is presented as an extension to a previous stationary model and based on an equilibrium of local forces. In the simplified model, ignoring the effects of membrane inertia, a dynamic equilibrium between the membrane bending potential and local fluid friction is considered in each phase. The equilibrium equations at the domain borders are completed by extended introduction of membrane section reactions. We show that in some cases, the results of stationary and evolutionary models are in agreement with each other and also with experimental observations, while in others the two models differ markedly. The value of our approach is that we can account for unresponsive points of uncertainty using our equations with the local velocity of the lipid membranes and calculating the intermediate states (shapes) in the consequent evolutionary, or response, path.

Relaxation dynamics of two-component fluid bilayer membranes

We theoretically investigate the relaxation dynamics of a nearly flat binary lipid bilayer membrane by taking into account the membrane tension, hydrodynamics of the surrounding fluid, inter-monolayer friction and mutual diffusion. Mutual diffusion is the collective irreversible process that leads to homogenization of the density difference between the two lipid species. We find that two relaxation modes associated with the mutual diffusion appear in addition to the three previously discussed relaxation modes reflecting the bending and compression of the membrane. Because of the symmetry, only one of the two diffusive modes is coupled to the bending mode. The two diffusive modes are much slower than the bending and compression modes in the entire realistic wave number range. This means that the long time relaxation behavior is dominated by the mutual diffusion in binary membranes. The two diffusive modes become even slower in the vicinity of the unstable region towards phase separation, while the other modes are almost unchanged. In short time scales, on the other hand, the lipid composition heterogeneity induces in-plane compression and bending of the bilayer.

A stable numerical method for the dynamics of fluidic membranes

We develop a finite element scheme to approximate the dynamics of two and three dimensional fluidic membranes in Navier-Stokes flow. Local inextensibility of the membrane is ensured by solving a tangential Navier-Stokes equation, taking surface viscosity effects of Boussinesq-Scriven type into account. In our approach the bulk and surface degrees of freedom are discretized independently, which leads to an unfitted finite element approximation of the underlying free boundary problem. Bending elastic forces resulting from an elastic membrane energy are discretized using an approximation introduced by Dziuk (Numer Math 111:55-80, 2008). The obtained numerical scheme can be shown to be stable and to have good mesh properties. Finally, the evolution of membrane shapes is studied numerically in different flow situations in two and three space dimensions. The numerical results demonstrate the robustness of the method, and it is observed that the conservation properties are fulfilled to a high precision.

Hydrodynamics of bilayer membranes with diffusing transmembrane proteins

We consider the hydrodynamics of lipid bilayers containing transmembrane proteins of arbitrary shape. This biologically-motivated problem is relevant to the cell membrane, whose fluctuating dynamics play a key role in phenomena ranging from cell migration, intercellular transport, and cell communication. Using Onsager's variational principle, we derive the equations that govern the relaxation dynamics of the membrane shape, of the mass densities of the bilayer leaflets, and of the diffusing proteins' concentration. With our generic formalism, we obtain several results on membrane dynamics. We find that proteins that span the bilayer increase the intermonolayer friction coefficient. The renormalization, which can be significant, is in inverse proportion to the protein's mobility. Second, we find that asymmetric proteins couple to the membrane curvature and to the difference in monolayer densities. For practically all accessible membrane tensions (σ > 10(-8) N m(-1)) we show that the protein density is the slowest relaxing variable. Furthermore, its relaxation rate decreases at small wavelengths due to the coupling to curvature. We apply our formalism to the large-scale diffusion of a concentrated protein patch. We find that the diffusion profile is not self-similar, owing to the wavevector dependence of the effective diffusion coefficient.

Numerical computations of the dynamics of fluidic membranes and vesicles

Vesicles and many biological membranes are made of two monolayers of lipid molecules and form closed lipid bilayers. The dynamical behavior of vesicles is very complex and a variety of forms and shapes appear. Lipid bilayers can be considered as a surface fluid and hence the governing equations for the evolution include the surface (Navier-)Stokes equations, which in particular take the membrane viscosity into account. The evolution is driven by forces stemming from the curvature elasticity of the membrane. In addition, the surface fluid equations are coupled to bulk (Navier-)Stokes equations. We introduce a parametric finite-element method to solve this complex free boundary problem and present the first three-dimensional numerical computations based on the full (Navier-)Stokes system for several different scenarios. For example, the effects of the membrane viscosity, spontaneous curvature, and area difference elasticity (ADE) are studied. In particular, it turns out, that even in the case of no viscosity contrast between the bulk fluids, the tank treading to tumbling transition can be obtained by increasing the membrane viscosity. Besides the classical tank treading and tumbling motions, another mode (called the transition mode in this paper, but originally called the vacillating-breathing mode and subsequently also called trembling, transition, and swinging mode) separating these classical modes appears and is studied by us numerically. We also study how features of equilibrium shapes in the ADE and spontaneous curvature models, like budding behavior or starfish forms, behave in a shear flow.

Protein-induced membrane curvature alters local membrane tension

Adsorption of proteins onto membranes can alter the local membrane curvature. This phenomenon has been observed in biological processes such as endocytosis, tubulation, and vesiculation. However, it is not clear how the local surface properties of the membrane, such as membrane tension, change in response to protein adsorption. In this article, we show that the partial differential equations arising from classical elastic model of lipid membranes, which account for simultaneous changes in shape and membrane tension due to protein adsorption in a local region, cannot be solved for nonaxisymmetric geometries using straightforward numerical techniques; instead, a viscous-elastic formulation is necessary to fully describe the system. Therefore, we develop a viscous-elastic model for inhomogeneous membranes of the Helfrich type. Using the newly available viscous-elastic model, we find that the lipids flow to accommodate changes in membrane curvature during protein adsorption. We show that, at the end of protein adsorption process, the system sustains a residual local tension to balance the difference between the actual mean curvature and the imposed spontaneous curvature. We also show that this change in membrane tension can have a functional impact such as altered response to pulling forces in the presence of proteins.

Interplay of packing and flip-flop in local bilayer deformation. How phosphatidylglycerol could rescue mitochondrial function in a cardiolipin-deficient yeast mutant

In a previous work, we have shown that a spatially localized transmembrane pH gradient, produced by acid micro-injection near the external side of cardiolipin-containing giant unilamellar vesicles, leads to the formation of tubules that retract after the dissipation of this gradient. These tubules have morphologies similar to mitochondrial cristae. The tubulation effect is attributable to direct phospholipid packing modification in the outer leaflet, that is promoted by protonation of cardiolipin headgroups. In this study, we compare the case of cardiolipin-containing giant unilamellar vesicles with that of giant unilamellar vesicles that contain phosphatidylglycerol (PG). Local acidification also promotes formation of tubules in the latter. However, compared with cardiolipin-containing giant unilamellar vesicles the tubules are longer, exhibit a visible pearling, and have a much longer lifetime after acid micro-injection is stopped. We attribute these differences to an additional mechanism that increases monolayer surface imbalance, namely inward PG flip-flop promoted by the local transmembrane pH gradient. Simulations using a fully nonlinear membrane model as well as geometrical calculations are in agreement with this hypothesis. Interestingly, among yeast mutants deficient in cardiolipin biosynthesis, only the crd1-null mutant, which accumulates phosphatidylglycerol, displays significant mitochondrial activity. Our work provides a possible explanation of such a property and further emphasizes the salient role of specific lipids in mitochondrial function.

Critical particle sizes for the engulfment of nanoparticles by membranes and vesicles with bilayer asymmetry

The adhesion and engulfment of nanoparticles by biomembranes is essential for many processes such as biomedical imaging, drug delivery, nanotoxicity, and viral infection. Many studies have shown that both surface chemistry, which determines the adhesive strength of the membrane-particle interactions, and particle size represent key parameters for these processes. Here, we show that the asymmetry between the two leaflets of a bilayer membrane provides another key parameter for the engulfment of nanoparticles. The asymmetric membrane prefers to curve in a certain manner as quantitatively described by its spontaneous curvature. We derive two general relationships between particle size, adhesive strength, and spontaneous curvature that determine the instabilities of (i) the nonadhering or free state and (ii) the completely engulfed state of the particle. For model membranes such as lipid or polymer bilayers with a uniform composition, the two relationships lead to two critical particle sizes that determine four distinct engulfment regimes, both for the endocytic and for the exocytic engulfment process. For strong adhesion, the critical particle sizes are on the order of 10 nm, while they are on the order of 1000 nm for weak or ultraweak adhesion. Our theoretical results are therefore accessible to both experimental studies and computer simulations of model membranes. In order to address the more complex process of receptor-mediated endocytosis, we take the adhesion-induced segregation of membrane components into account and consider bound and unbound membrane segments that differ in their spontaneous curvatures. To model protein coats as formed during clathrin-dependent endocytosis, we focus on the case in which the bound membrane segments have a large spontaneous curvature compared to the unbound ones. We derive explicit expressions for the engulfment rate and the uptake of nanoparticles, which both depend on the particle size in a nonmonotonic manner, and provide a quantitative fit to experimental data for clathrin-dependent endocytosis of gold nanoparticles.

Effects of local pH on the formation and regulation of cristae morphologies

Cristae, folded subcompartments of the inner mitochondrial membrane (IMM), have complex and dynamic morphologies. Since cristae are the major site of adenosine triphosphate synthesis, morphological changes of cristae have been studied in relation to functional states of mitochondria. In this sense, investigating the functional and structural significance of cristae may be critical for understanding progressive mitochondrial dysfunction. However, the detailed mechanisms of the formation and regulation of these cristae structures have not been fully elucidated. Among the hypotheses concerning the regulation of cristae morphologies, we exclusively investigate the effects of the local pH gradient on the cristae morphologies by using a numerical model. An area-difference induced curvature of the membrane is modeled as a function of local pH. This curvature is then applied to the finite element model of a closed lipid bilayer in order to find the energetically favorable membrane configuration. From this study, we substantiate the hypothesis that a tubular crista structure can be formed and regulated by the local pH gradient. Through the simulations with various initial conditions, we further demonstrate that the diameter of a crista is mainly determined by the local pH gradient, and the energetically favorable direction of crista growth is perpendicular to the longitudinal axis of a mitochondrion. Finally, the simulation results at the mitochondrial scale suggest that the cristae membrane may have a lower local pH value and/or a higher cardiolipin composition than the other parts of the IMM.

Dynamics and instabilities of lipid bilayer membrane shapes

Biological membranes undergo constant shape remodeling involving the formation of highly curved structures. The lipid bilayer represents the fundamental architecture of the cellular membrane with its shapes determined by the Helfrich curvature bending energy. However, the dynamics of bilayer shape transitions, especially their modulation by membrane proteins, and the resulting shape instabilities, are still not well understood. Here, we review in a unifying manner several theories that describe the fluctuations (i.e. undulations) of bilayer shapes as well as their local coupling with lipid or protein density variation. The coupling between local membrane curvature and lipid density gives rise to a 'slipping mode' in addition to the conventional 'bending mode' for damping the membrane fluctuation. This leads to a number of interesting experimental phenomena regarding bilayer shape dynamics. More importantly, curvature-inducing proteins can couple with membrane shape and eventually render the membrane unstable. A criterion for membrane shape instability is derived from a linear stability analysis. The instability criterion reemphasizes the importance of membrane tension in regulating the stability and dynamics of membrane geometry. Recent progresses in understanding the role of membrane tension in regulating dynamical cellular processes are also reviewed. Protein density is emphasized as a key factor in regulating membrane shape transitions: a threshold density of curvature coupling proteins is required for inducing membrane morphology transitions.

Fluid vesicles in flow

We review the dynamical behavior of giant fluid vesicles in various types of external hydrodynamic flow. The interplay between stresses arising from membrane elasticity, hydrodynamic flows, and the ever present thermal fluctuations leads to a rich phenomenology. In linear flows with both rotational and elongational components, the properties of the tank-treading and tumbling motions are now well described by theoretical and numerical models. At the transition between these two regimes, strong shape deformations and amplification of thermal fluctuations generate a new regime called trembling. In this regime, the vesicle orientation oscillates quasi-periodically around the flow direction while asymmetric deformations occur. For strong enough flows, small-wavelength deformations like wrinkles are observed, similar to what happens in a suddenly reversed elongational flow. In steady elongational flow, vesicles with large excess areas deform into dumbbells at large flow rates and pearling occurs for even stronger flows. In capillary flows with parabolic flow profile, single vesicles migrate towards the center of the channel, where they adopt symmetric shapes, for two reasons. First, walls exert a hydrodynamic lift force which pushes them away. Second, shear stresses are minimal at the tip of the flow. However, symmetry is broken for vesicles with large excess areas, which flow off-center and deform asymmetrically. In suspensions, hydrodynamic interactions between vesicles add up to these two effects, making it challenging to deduce rheological properties from the dynamics of individual vesicles. Further investigations of vesicles and similar objects and their suspensions in steady or time-dependent flow will shed light on phenomena such as blood flow.

Confined bilayers passively regulate shape and stress

Lipid membranes are commonly confined to adjacent subcellular structures or to artificial substrates and particles. We develop an experimental and theoretical framework to investigate the mechanics of confined membranes, including the influence of adhesion, strain, and osmotic pressure. We find that supported lipid bilayers respond to stress by nucleating and evolving spherical and tubular protrusions. In cells, such transformations are generally attributed to proteins. Our results offer insights into the mechanics of cell membranes and can further extend the applications of supported bilayers.