Relaxation dynamics of fluid membranes

OsSRF8 interacts with OsINP1 and OsDAF1 to regulate pollen aperture formation in rice

In higher plants, mature male gametophytes have distinct apertures. After pollination, pollen grains germinate, and a pollen tube grows from the aperture to deliver sperm cells to the embryo sac, completing fertilization. In rice, the pollen aperture has a single-pore structure with a collar-like annulus and a plug-like operculum. A crucial step in aperture development is the formation of aperture plasma membrane protrusion (APMP) at the distal polar region of the microspore during the late tetrad stage. Previous studies identified OsINP1 and OsDAF1 as essential regulators of APMP and pollen aperture formation in rice, but their precise molecular mechanisms remain unclear. We demonstrate that the Poaceae-specific OsSRF8 gene, encoding a STRUBBELIG-receptor family 8 protein, is essential for pollen aperture formation in Oryza sativa. Mutants lacking functional OsSRF8 exhibit defects in APMP and pollen aperture formation, like loss-of-function OsINP1 mutants. OsSRF8 is specifically expressed during early anther development and initially diffusely distributed in the microsporocytes. At the tetrad stage, OsSRF8 is recruited by OsINP1 to the pre-aperture region through direct protein-protein interaction, promoting APMP formation. The OsSRF8-OsINP1 complex then recruits OsDAF1 to the APMP site to co-regulate annulus formation. Our findings provide insights into the mechanisms controlling pollen aperture formation in cereal species.

(2+δ)-dimensional theory of the electromechanics of lipid membranes: Electrostatics

The coupling of electric fields to the mechanics of lipid membranes gives rise to intriguing electromechanical behavior, as, for example, evidenced by the deformation of lipid vesicles in external electric fields. Electromechanical effects are relevant for many biological processes, such as the propagation of action potentials in axons and the activation of mechanically gated ion channels. Currently, a theoretical framework describing the electromechanical behavior of arbitrarily curved and deforming lipid membranes does not exist. Purely mechanical models commonly treat lipid membranes as two-dimensional surfaces, ignoring their finite thickness. While holding analytical and numerical merit, this approach cannot describe the coupling of lipid membranes to electric fields and is thus unsuitable for electromechanical models. In a sequence of articles, we derive an effective surface theory of the electromechanics of lipid membranes, called the (2+δ)-dimensional theory, which has the advantages of surface descriptions while accounting for finite thickness effects. The present article proposes a generic dimension reduction procedure relying on low-order spectral expansions. This procedure is applied to the electrostatics of lipid membranes to obtain the (2+δ)-dimensional theory that captures potential differences across and electric fields within lipid membranes. This model is tested on different geometries relevant for lipid membranes, showing good agreement with the corresponding three-dimensional electrostatics theory.

Dynamic framework for large-scale modeling of membranes and peripheral proteins

In this chapter, we present a novel computational framework to study the dynamic behavior of extensive membrane systems, potentially in interaction with peripheral proteins, as an alternative to conventional simulation methods. The framework effectively describes the complex dynamics in protein-membrane systems in a mesoscopic particle-based setup. Furthermore, leveraging the hydrodynamic coupling between the membrane and its surrounding solvent, the coarse-grained model grounds its dynamics in macroscopic kinetic properties such as viscosity and diffusion coefficients, marrying the advantages of continuum- and particle-based approaches. We introduce the theoretical background and the parameter-space optimization method in a step-by-step fashion, present the hydrodynamic coupling method in detail, and demonstrate the application of the model at each stage through illuminating examples. We believe this modeling framework to hold great potential for simulating membrane and protein systems at biological spatiotemporal scales, and offer substantial flexibility for further development and parametrization.

Kinetic pathway and micromechanics of fusion/fission for polyelectrolyte vesicles

Despite the wide existence of vesicles in living cells as well as their important applications like drug delivery, the underlying mechanism of vesicle fusion/fission remains under debate. Classical models cannot fully explain recent observations in experiments and simulations. Here, we develop a constrained self-consistent field theory that allows tracking the shape evolution and free energy as a function of center-of-mass separation distance. Fusion and fission are described in a unified framework. Both the kinetic pathway and the mechanical response can be simultaneously captured. By taking vesicles formed by polyelectrolytes as a model system, we predict discontinuous transitions between the three morphologies: parent vesicle with a single cavity, hemifission/hemifusion, and two separated child vesicles, as a result of breaking topological isomorphism. With the increase in inter-vesicle repulsion, we observe a great reduction in the cleavage energy, indicating that vesicle fission can be achieved without hemifission, in good agreement with simulation results. The force-extension relationship elucidates typical plasticity for separating two vesicles. The super extensibility in the mechanical response of vesicle is in stark contrast to soft particles with other morphologies, such as cylinder and sphere. Our work elucidates the fundamental physical chemistry based on intrinsic topological features of vesicle fusion/fission, which provides insights into various phenomena observed in experiments and simulations.

Axis convergence in C. elegans embryos

Embryos develop in a surrounding that guides key aspects of their development. For example, the anteroposterior (AP) body axis is always aligned with the geometric long axis of the surrounding eggshell in fruit flies and worms. The mechanisms that ensure convergence of the AP axis with the long axis of the eggshell remain unresolved. We investigate axis convergence in early C. elegans development, where the nascent AP axis, when misaligned, actively re-aligns to converge with the long axis of the egg. We identify two physical mechanisms that underlie axis convergence. First, bulk cytoplasmic flows, driven by actomyosin cortical flows, can directly reposition the AP axis. Second, active forces generated within the pseudocleavage furrow, a transient actomyosin structure similar to a contractile ring, can drive a mechanical re-orientation such that it becomes positioned perpendicular to the long axis of the egg. This in turn ensures AP axis convergence. Numerical simulations, together with experiments that either abolish the pseudocleavage furrow or change the shape of the egg, demonstrate that the pseudocleavage-furrow-dependent mechanism is a major driver of axis convergence. We conclude that active force generation within the actomyosin cortical layer drives axis convergence in the early nematode.

Geometry-structure models for liquid crystal interfaces, drops and membranes: wrinkling, shape selection and dissipative shape evolution

We review our recent contributions to anisotropic soft matter models for liquid crystal interfaces, drops and membranes, emphasizing validations with experimental and biological data, and with related theory and simulation literature. The presentation aims to illustrate and characterize the rich output and future opportunities of using a methodology based on the liquid crystal-membrane shape equation applied to static and dynamic pattern formation phenomena. The geometry of static and kinetic shapes is usually described with dimensional curvatures that co-mingle shape and curvedness. In this review, we systematically show how the application of a novel decoupled shape-curvedness framework to practical and ubiquitous soft matter phenomena, such as the shape of drops and tactoids and bending of evolving membranes, leads to deeper quantitative insights than when using traditional dimensional mean and Gaussian curvatures. The review focuses only on (1) statics of wrinkling and shape selection in liquid crystal interfaces and membranes; (2) kinetics and dissipative dynamics of shape evolution in membranes; and (3) computational methods for shape selection and shape evolution; due to various limitations other important topics are excluded. Finally, the outlook follows a similar structure. The main results include: (1) single and multiple wavelength corrugations in liquid crystal interfaces appear naturally in the presence of surface splay and bend orientation distortions with scaling laws governed by ratios of anchoring-to-isotropic tension energy; adding membrane elasticity to liquid crystal anchoring generates multiple scales wrinkling as in tulips; drops of liquid crystals encapsulates in membranes can adopt, according to the ratios of anchoring/tension/bending, families of shapes as multilobal, tactoidal, and serrated as observed in biological cells. (2) Mapping the liquid crystal director to a membrane unit normal. The dissipative shape evolution model with irreversible thermodynamics for flows dominated by bending rates, yields new insights. The model explains the kinetic stability of cylinders, while spheres and saddles are attractors. The model also adds to the evolving understanding of outer hair cells in the inner ear. (3) Computational soft matter geometry includes solving shape equations, trajectories on energy and orientation landscapes, and shape-curvedness evolutions on entropy production landscape with efficient numerical methods and adaptive approaches.

A mechanosensing mechanism controls plasma membrane shape homeostasis at the nanoscale

As cells migrate and experience forces from their surroundings, they constantly undergo mechanical deformations which reshape their plasma membrane (PM). To maintain homeostasis, cells need to detect and restore such changes, not only in terms of overall PM area and tension as previously described, but also in terms of local, nanoscale topography. Here, we describe a novel phenomenon, by which cells sense and restore mechanically induced PM nanoscale deformations. We show that cell stretch and subsequent compression reshape the PM in a way that generates local membrane evaginations in the 100 nm scale. These evaginations are recognized by I-BAR proteins, which triggers a burst of actin polymerization mediated by Rac1 and Arp2/3. The actin polymerization burst subsequently re-flattens the evagination, completing the mechanochemical feedback loop. Our results demonstrate a new mechanosensing mechanism for PM shape homeostasis, with potential applicability in different physiological scenarios.

Inhomogeneous Canham-Helfrich Abscission in Catenoid Necks under Critical Membrane Mosaicity

The mechanical effects of membrane compositional inhomogeneities are analyzed in a process analogous to neck formation in cellular membranes. We cast on the Canham-Helfrich model of fluid membranes with both the spontaneous curvature and the surface tension being non-homogeneous functions along the cell membrane. The inhomogeneous distribution of necking forces is determined by the equilibrium mechanical equations and the boundary conditions as considered in the axisymmetric setting compatible with the necking process. To establish the role played by mechanical inhomogeneity, we focus on the catenoid, a surface of zero mean curvature. Analytic solutions are shown to exist for the spontaneous curvature and the constrictive forces in terms of the border radii. Our theoretical analysis shows that the inhomogeneous distribution of spontaneous curvature in a mosaic-like neck constrictional forces potentially contributes to the membrane scission under minimized work in living cells.

Solving the Incompressible Surface Stokes Equation by Standard Velocity-Correction Projection Methods

In this paper, an effective numerical algorithm for the Stokes equation of a curved surface is presented and analyzed. The velocity field was decoupled from the pressure by the standard velocity correction projection method, and the penalty term was introduced to make the velocity satisfy the tangential condition. The first-order backward Euler scheme and second-order BDF scheme are used to discretize the time separately, and the stability of the two schemes is analyzed. The mixed finite element pair (P2,P1) is applied to discretization of space. Finally, numerical examples are given to verify the accuracy and effectiveness of the proposed method.

Membrane area gain and loss during cytokinesis

In cytokinesis of animal cells, the cell is symmetrically divided into two. Since the cell's volume is conserved, the projected area has to increase to allow for the change of shape. Here we aim to predict how membrane gain and loss adapt during cytokinesis. We work with a kinetic model in which membrane turnover depends on membrane tension and cell shape. We apply this model to a series of calculated vesicle shapes as a proxy for the shape of dividing cells. We find that the ratio of kinetic turnover parameters changes nonmonotonically with cell shape, determined by the dependence of exocytosis and endocytosis on membrane curvature. Our results imply that controlling membrane turnover will be crucial for the successful division of artificial cells.

Visceral organ morphogenesis via calcium-patterned muscle constrictions

Organ architecture is often composed of multiple laminar tissues arranged in concentric layers. During morphogenesis, the initial geometry of visceral organs undergoes a sequence of folding, adopting a complex shape that is vital for function. Genetic signals are known to impact form, yet the dynamic and mechanical interplay of tissue layers giving rise to organs' complex shapes remains elusive. Here, we trace the dynamics and mechanical interactions of a developing visceral organ across tissue layers, from subcellular to organ scale in vivo. Combining deep tissue light-sheet microscopy for in toto live visualization with a novel computational framework for multilayer analysis of evolving complex shapes, we find a dynamic mechanism for organ folding using the embryonic midgut of Drosophila as a model visceral organ. Hox genes, known regulators of organ shape, control the emergence of high-frequency calcium pulses. Spatiotemporally patterned calcium pulses trigger muscle contractions via myosin light chain kinase. Muscle contractions, in turn, induce cell shape change in the adjacent tissue layer. This cell shape change collectively drives a convergent extension pattern. Through tissue incompressibility and initial organ geometry, this in-plane shape change is linked to out-of-plane organ folding. Our analysis follows tissue dynamics during organ shape change in vivo, tracing organ-scale folding to a high-frequency molecular mechanism. These findings offer a mechanical route for gene expression to induce organ shape change: genetic patterning in one layer triggers a physical process in the adjacent layer – revealing post-translational mechanisms that govern shape change.

The morphological role of ligand inhibitors in blocking receptor- and clathrin-mediated endocytosis

Cells often internalize particles through endocytic pathways that involve the binding between cell receptors and particle ligands, which drives the cell membrane to wrap the particle into a delivery vesicle. Previous findings showed that receptor-mediated endocytosis is impossible for spherical particles smaller than a minimum size because of the energy barrier created by membrane bending. In this study, we investigate the morphological role of ligand inhibitors in blocking endocytosis, inspired by antibodies that inhibit virus ligands to prevent infection. While ligand inhibitors have the obvious effect of reducing the driving force due to adhesion, they also have a nontrivial (morphological) impact on the entropic and elastic energy of the system. We determine the necessary conditions for endocytosis by considering the additional energy barrier due to the membrane bending to wrap the inhibiting protrusions. We find that inhibitors increase the minimum radius previously reported, depending on their density and size. In addition, we extend this result to the case of clathrin-mediated endocytosis, which is the most common pathway for virus entry. The assembly of a clathrin coat with a spontaneous curvature increases the energy barrier and sets a maximum particle size (in agreement with experimental observations on spherical particles). Our investigation suggests that morphological considerations can inform the optimal design of neutralizing viral antibodies and new strategies for targeted nanomedicine.

Microfluidics Approach to the Mechanical Properties of Red Blood Cell Membrane and Their Effect on Blood Rheology

In this article, we describe the general features of red blood cell membranes and their effect on blood flow and blood rheology. We first present a basic description of membranes and move forward to red blood cell membranes’ characteristics and modeling. We later review the specific properties of red blood cells, presenting recent numerical and experimental microfluidics studies that elucidate the effect of the elastic properties of the red blood cell membrane on blood flow and hemorheology. Finally, we describe specific hemorheological pathologies directly related to the mechanical properties of red blood cells and their effect on microcirculation, reviewing microfluidic applications for the diagnosis and treatment of these diseases.

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

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

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.

The Mechanics and Thermodynamics of Tubule Formation in Biological Membranes

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

Non-Equilibrium Large-Scale Membrane Transformations Driven by MinDE Biochemical Reaction Cycles

The MinDE proteins from E. coli have received great attention as a paradigmatic biological pattern-forming system. Recently, it has surfaced that these proteins do not only generate oscillating concentration gradients driven by ATP hydrolysis, but that they can reversibly deform giant vesicles. In order to explore the potential of Min proteins to actually perform mechanical work, we introduce a new model membrane system, flat vesicle stacks on top of a supported lipid bilayer. MinDE oscillations can repeatedly deform these flat vesicles into tubules and promote progressive membrane spreading through membrane adhesion. Dependent on membrane and buffer compositions, Min oscillations further induce robust bud formation. Altogether, we demonstrate that under specific conditions, MinDE self-organization can result in work performed on biomimetic systems and achieve a straightforward mechanochemical coupling between the MinDE biochemical reaction cycle and membrane transformation.

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.

Dynamics of passive and active membrane tubes

Utilising Onsager's variational formulation, we derive dynamical equations for the relaxation of a fluid membrane tube in the limit of small deformation, allowing for a contrast of solvent viscosity across the membrane and variations in surface tension due to membrane incompressibility. We compute the relaxation rates, recovering known results in the case of purely axis-symmetric perturbations and making new predictions for higher order (azimuthal) m-modes. We analyse the long and short wavelength limits of these modes by making use of various asymptotic arguments. We incorporate stochastic terms to our dynamical equations suitable to describe both passive thermal forces and non-equilibrium active forces. We derive expressions for the fluctuation amplitudes, an effective temperature associated with active fluctuations, and the power spectral density for both the thermal and active fluctuations. We discuss an experimental assay that might enable measurement of these fluctuations to infer the properties of the active noise. Finally we discuss our results in the context of active membranes more generally and give an overview of some open questions in the field.

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.

Shear-Driven Instabilities of Membrane Tubes and Dynamin-Induced Scission

Motivated by the mechanics of dynamin-mediated membrane tube fission, we analyze the stability of fluid membrane tubes subjected to shear flow in azimuthal direction. We find a novel helical instability driven by the membrane shear flow which results in a nonequilibrium steady state for the tube fluctuations. This instability has its onset at shear rates that may be physiologically accessible under the action of dynamin and could also be probed using in vitro experiments on membrane nanotubes, e.g., using magnetic tweezers. We discuss how such an instability may play a role in the mechanism for dynamin-mediated membrane tube fission.

Large-scale simulation of biomembranes incorporating realistic kinetics into coarse-grained models

Biomembranes are two-dimensional assemblies of phospholipids that are only a few nanometres thick, but form micrometre-sized structures vital to cellular function. Explicit molecular modelling of biologically relevant membrane systems is computationally expensive due to the large number of solvent particles and slow membrane kinetics. Coarse-grained solvent-free membrane models offer efficient sampling but sacrifice realistic kinetics, thereby limiting the ability to predict pathways and mechanisms of membrane processes. Here, we present a framework for integrating coarse-grained membrane models with continuum-based hydrodynamics. This framework facilitates efficient simulation of large biomembrane systems with large timesteps, while achieving realistic equilibrium and non-equilibrium kinetics. It helps to bridge between the nanometer/nanosecond spatiotemporal resolutions of coarse-grained models and biologically relevant time- and lengthscales. As a demonstration, we investigate fluctuations of red blood cells, with varying cytoplasmic viscosities, in 150-milliseconds-long trajectories, and compare kinetic properties against single-cell experimental observations.

Geometry and dynamics of lipid membranes: The Scriven-Love number

The equations governing lipid membrane dynamics in planar, spherical, and cylindrical geometries are presented here. Unperturbed and first-order perturbed equations are determined and nondimensionalized. In membrane systems with a nonzero base flow, perturbed in-plane and out-of-plane quantities are found to vary over different length scales. A new dimensionless number, named the Scriven-Love number, and the well-known Föppl-von Kármán number result from a scaling analysis. The Scriven-Love number compares out-of-plane forces arising from the in-plane, intramembrane viscous stresses to the familiar elastic bending forces, while the Föppl-von Kármán number compares tension to bending forces. Both numbers are calculated in past experimental works, and span a wide range of values in various biological processes across different geometries. In situations with large Scriven-Love and Föppl-von Kármán numbers, the dynamical response of a perturbed membrane is dominated by out-of-plane viscous and surface tension forces-with bending forces playing a negligible role. Calculations of non-negligible Scriven-Love numbers in various biological processes and in vitro experiments show in-plane intramembrane viscous flows cannot generally be ignored when analyzing lipid membrane behavior.

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.

Mechanics of nuclear membranes

Cellular nuclei are bound by two uniformly separated lipid membranes that are fused with each other at numerous donut-shaped pores. These membranes are structurally supported by an array of distinct proteins with distinct mechanical functions. As a result, the nuclear envelope possesses unique mechanical properties, which enables it to resist cytoskeletal forces. Here, we review studies that are beginning to provide quantitative insights into nuclear membrane mechanics. We discuss how the mechanical properties of the fused nuclear membranes mediate their response to mechanical forces exerted on the nucleus and how structural reinforcement by different nuclear proteins protects the nuclear membranes against rupture. We also highlight some open questions in nuclear envelope mechanics, and discuss their relevance in the context of health and disease.

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.

The Effects of Intra-membrane Viscosity on Lipid Membrane Morphology: Complete Analytical Solution

We present a linear theory of lipid membranes which accommodates the effects of intra-membrane viscosity into the model of deformations. Within the Monge parameterization, a linearized version of the shape equation describing membrane morphology is derived. Admissible boundary conditions are taken from the existing non-linear model but reformulated and adopted to the present framework. We obtain a complete analytical expression illustrating the deformations of lipid membrane subjected to the influences of intra-membrane viscosity. The result predicts wrinkle phenomena in the event of membrane-substrate interactions. Finally, we mention that the obtained solutions reduce to those from the classical shape equation when the viscosity effects are removed.

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.

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.

Generalized Boussinesq-Scriven surface fluid model with curvature dissipation for liquid surfaces and membranes

Curvature dissipation is relevant in synthetic and biological processes, from fluctuations in semi-flexible polymer solutions, to buckling of liquid columns, tomembrane cell wall functioning. We present a micromechanical model of curvature dissipation relevant to fluid membranes and liquid surfaces based on a parallel surface parameterization and a stress constitutive equation appropriate for anisotropic fluids and fluid membranes.The derived model, aimed at high curvature and high rate of change of curvature in liquid surfaces and membranes, introduces additional viscous modes not included in the widely used 2D Boussinesq-Scriven rheological constitutive equation for surface fluids.The kinematic tensors that emerge from theparallel surface parameterization are the interfacial rate of deformation and the surface co-rotational Zaremba-Jaumann derivative of the curvature, which are used to classify all possibledissipative planar and non-planar modes. The curvature dissipation function that accounts for bending, torsion and twist rates is derived and analyzed under several constraints, including the important inextensional bending mode.A representative application of the curvature dissipation model to the periodic oscillation in nano-wrinkled outer hair cells show how and why curvature dissipation decreases with frequency, and why the 100kHz frequency range is selected. These results contribute to characterize curvature dissipation in membranes and liquid surfaces.

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.

Coupling of lipid membrane elasticity and in-plane dynamics

Biomembranes exhibit liquid and solid features concomitantly with their in-plane fluidity and elasticity tightly regulated by cells. Here, we present experimental evidence supporting the existence of the dynamics-elasticity correlations for lipid membranes and propose a mechanism involving molecular packing densities to explain them. This paper thereby unifies, at the molecular level, the aspects of the continuum mechanics long used to model the two membrane features. This ultimately may elucidate the universal physical principles governing the cellular phenomena involving biomembranes.

Signatures of Mechanosensitive Gating

The question of how mechanically gated membrane channels open and close is notoriously difficult to address, especially if the protein structure is not available. This perspective highlights the relevance of micropipette-aspirated single-particle tracking-used to obtain a channel's diffusion coefficient, D, as a function of applied membrane tension, σ-as an indirect assay for determining functional behavior in mechanosensitive channels. While ensuring that the protein remains integral to the membrane, such methods can be used to identify not only the gating mechanism of a protein, but also associated physical moduli, such as torsional and dilational rigidity, which correspond to the protein's effective shape change. As an example, three distinct D-versus-σ "signatures" are calculated, corresponding to gating by dilation, gating by tilt, and gating by a combination of both dilation and tilt. Both advantages and disadvantages of the approach are discussed.

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.

Energy dissipation in flows through curved spaces

Fluid dynamics in intrinsically curved geometries is encountered in many physical systems in nature, ranging from microscopic bio-membranes all the way up to general relativity at cosmological scales. Despite the diversity of applications, all of these systems share a common feature: the free motion of particles is affected by inertial forces originating from the curvature of the embedding space. Here we reveal a fundamental process underlying fluid dynamics in curved spaces: the free motion of fluids, in the complete absence of solid walls or obstacles, exhibits loss of energy due exclusively to the intrinsic curvature of space. We find that local sources of curvature generate viscous stresses as a result of the inertial forces. The curvature- induced viscous forces are shown to cause hitherto unnoticed and yet appreciable energy dissipation, which might play a significant role for a variety of physical systems involving fluid dynamics in curved spaces.

Relaxation dynamics of a compressible bilayer vesicle containing highly viscous fluid

We study the relaxation dynamics of a compressible bilayer vesicle with an asymmetry in the viscosity of the inner and outer fluid medium. First we explore the stability of the vesicle free energy which includes a coupling between the membrane curvature and the local density difference between the two monolayers. Two types of instabilities are identified: a small wavelength instability and a larger wavelength instability. Considering the bulk fluid viscosity and the inter-monolayer friction as the dissipation sources, we next employ Onsager's variational principle to derive the coupled equations both for the membrane and the bulk fluid. The three relaxation modes are coupled to each other due to the bilayer and the spherical structure of the vesicle. Most importantly, a higher fluid viscosity inside the vesicle shifts the crossover mode between the bending and the slipping to a larger value. As the vesicle parameters approach the unstable regions, the relaxation dynamics is dramatically slowed down, and the corresponding mode structure changes significantly. In some limiting cases, our general result reduces to the previously obtained relaxation rates.

Hydrodynamic coupling of particle inclusions embedded in curved lipid bilayer membranes

We develop theory and computational methods to investigate particle inclusions embedded within curved lipid bilayer membranes. We consider the case of spherical lipid vesicles where inclusion particles are coupled through (i) intramembrane hydrodynamics, (ii) traction stresses with the external and trapped solvent fluid, and (iii) intermonolayer slip between the two leaflets of the bilayer. We investigate relative to flat membranes how the membrane curvature and topology augment hydrodynamic responses. We show how both the translational and rotational mobility of protein inclusions are effected by the membrane curvature, ratio of intramembrane viscosity to solvent viscosity, and intermonolayer slip. For general investigations of many-particle dynamics, we also discuss how our approaches can be used to treat the collective diffusion and hydrodynamic coupling within spherical bilayers.

Poiseuille flow in curved spaces

We investigate Poiseuille channel flow through intrinsically curved media, equipped with localized metric perturbations. To this end, we study the flux of a fluid driven through the curved channel in dependence of the spatial deformation, characterized by the parameters of the metric perturbations (amplitude, range, and density). We find that the flux depends only on a specific combination of parameters, which we identify as the average metric perturbation, and derive a universal flux law for the Poiseuille flow. For the purpose of this study, we have improved and validated our recently developed lattice Boltzmann model in curved space by considerably reducing discrete lattice effects.

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.

Mobility Measurements Probe Conformational Changes in Membrane Proteins due to Tension

The function of membrane-embedded proteins such as ion channels depends crucially on their conformation. We demonstrate how conformational changes in asymmetric membrane proteins may be inferred from measurements of their diffusion. Such proteins cause local deformations in the membrane, which induce an extra hydrodynamic drag on the protein. Using membrane tension to control the magnitude of the deformations, and hence the drag, measurements of diffusivity can be used to infer-via an elastic model of the protein-how conformation is changed by tension. Motivated by recent experimental results [Quemeneur et al., Proc. Natl. Acad. Sci. U.S.A. 111, 5083 (2014)], we focus on KvAP, a voltage-gated potassium channel from Aeropyrum pernix. The conformation of KvAP is found to change considerably due to tension, with its "walls," where the protein meets the membrane, undergoing significant angular strains. The torsional stiffness is determined to be 26.8k(B)T per radian at room temperature. This has implications for both the structure and the function of such proteins in the environment of a tension-bearing membrane.

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.