Self-organized shape dynamics of active surfaces

Cell-intrinsic mechanical regulation of plasma membrane accumulation at the cytokinetic furrow

Cytokinesis is the process where the mother cell's cytoplasm separates into daughter cells. This is driven by an actomyosin contractile ring that produces cortical contractility and drives cleavage furrow ingression, resulting in the formation of a thin intercellular bridge. While cytoskeletal reorganization during cytokinesis has been extensively studied, less is known about the spatiotemporal dynamics of the plasma membrane. Here, we image and model plasma membrane lipid and protein dynamics on the cell surface during leukemia cell cytokinesis. We reveal an extensive accumulation and folding of the plasma membrane at the cleavage furrow and the intercellular bridge, accompanied by a depletion and unfolding of the plasma membrane at the cell poles. These membrane dynamics are caused by two actomyosin-driven biophysical mechanisms: the radial constriction of the cleavage furrow causes local compression of the apparent cell surface area and accumulation of the plasma membrane at the furrow, while actomyosin cortical flows drag the plasma membrane toward the cell division plane as the furrow ingresses. The magnitude of these effects depends on the plasma membrane fluidity, cortex adhesion, and cortical contractility. Overall, our work reveals cell-intrinsic mechanical regulation of plasma membrane accumulation at the cleavage furrow that is likely to generate localized differences in membrane tension across the cytokinetic cell. This may locally alter endocytosis, exocytosis, and mechanotransduction, while also serving as a self-protecting mechanism against cytokinesis failures that arise from high membrane tension at the intercellular bridge.

Quantitative comparison of cell-cell detachment force in different experimental setups

We compare three different setups for measuring cell-cell adhesion. We show that the measured strength depends on the type of setup that is used. For identical cells different assays measure different detachment forces. This can be understood from the fact that cell-cell detachment is a global property of the system. We also analyse the role of external force and line tension on contact angle and cell-cell detachment. Comparison with the experiments suggest that viscous forces play an important role in the process. We dedicate this article to Fyl Pincus who for many of us is an example to be followed not only for outstanding science but also for a marvelous human behavior.

Cell intrinsic mechanical regulation of plasma membrane accumulation at the cytokinetic furrow

Cytokinesis is the process where the mother cell's cytoplasm separates into daughter cells. This is driven by an actomyosin contractile ring that produces cortical contractility and drives cleavage furrow ingression, resulting in the formation of a thin intercellular bridge. While cytoskeletal reorganization during cytokinesis has been extensively studied, little is known about the spatiotemporal dynamics of the plasma membrane. Here, we image and model plasma membrane lipid and protein dynamics on the cell surface during leukemia cell cytokinesis. We reveal an extensive accumulation and folding of plasma membrane at the cleavage furrow and the intercellular bridge, accompanied by a depletion and unfolding of plasma membrane at the cell poles. These membrane dynamics are caused by two actomyosin-driven biophysical mechanisms: the radial constriction of the cleavage furrow causes local compression of the apparent cell surface area and accumulation of the plasma membrane at the furrow, while actomyosin cortical flows drag the plasma membrane towards the cell division plane as the furrow ingresses. The magnitude of these effects depends on the plasma membrane fluidity, cortex adhesion and cortical contractility. Overall, our work reveals cell intrinsic mechanical regulation of plasma membrane accumulation at the cleavage furrow that is likely to generate localized differences in membrane tension across the cytokinetic cell. This may locally alter endocytosis, exocytosis and mechanotransduction, while also serving as a self-protecting mechanism against cytokinesis failures that arise from high membrane tension at the intercellular bridge.

Shaping epithelial lumina under pressure

The formation of fluid- or gas-filled lumina surrounded by epithelial cells pervades development and disease. We review the balance between lumen pressure and mechanical forces from the surrounding cells that governs lumen formation. We illustrate the mechanical side of this balance in several examples of increasing complexity, and discuss how recent work is beginning to elucidate how nonlinear and active mechanics and anisotropic biomechanical structures must conspire to overcome the isotropy of pressure to form complex, non-spherical lumina.

Collective chemomechanical oscillations in active hydrogels

We report on the collective response of an assembly of chemomechanical Belousov-Zhabotinsky (BZ) hydrogel beads. We first demonstrate that a single isolated spherical BZ hydrogel bead with a radius below a critical value does not oscillate, whereas an assembly of the same BZ hydrogel beads presents chemical oscillation. A BZ chemical model with an additional flux of chemicals out of the BZ hydrogel captures the experimentally observed transition from oxidized nonoscillating to oscillating BZ hydrogels and shows this transition is due to a flux of inhibitors out of the BZ hydrogel. The model also captures the role of neighboring BZ hydrogel beads in decreasing the critical size for an assembly of BZ hydrogel beads to oscillate. We finally leverage the quorum sensing behavior of the collective to trigger their chemomechanical oscillation and discuss how this collective effect can be used to enhance the oscillatory strain of these active BZ hydrogels. These findings could help guide the eventual fabrication of a swarm of autonomous, communicating, and motile hydrogels.

Odd elasticity and topological waves in active surfaces

Odd elasticity describes active elastic systems whose stress-strain relationship is not compatible with a potential energy. As the requirement of energy conservation is lifted from linear elasticity, new antisymmetric (odd) components appear in the elastic tensor. In this work we study the odd elasticity and non-Hermitian wave dynamics of active surfaces, specifically plates of moderate thickness. These odd moduli can endow the vibrational modes of the plate with a nonzero topological invariant known as the first Chern number. Within continuum elastic theory, we show that the Chern number is related to the presence of unidirectional shearing waves that are hosted at the plate's boundary. We show that the existence of these chiral edge waves hinges on a distinctive two-step mechanism. Unlike electronic Chern insulators where the magnetic field at the same time gaps the spectrum and imparts chirality, here the finite thickness of the sample gaps the shear modes, and the odd elasticity makes them chiral.

Forceful patterning: theoretical principles of mechanochemical pattern formation

Biological pattern formation is essential for generating and maintaining spatial structures from the scale of a single cell to tissues and even collections of organisms. Besides biochemical interactions, there is an important role for mechanical and geometrical features in the generation of patterns. We review the theoretical principles underlying different types of mechanochemical pattern formation across spatial scales and levels of biological organization.

TubULAR: tracking in toto deformations of dynamic tissues via constrained maps

A common motif in biology is the arrangement of cells into tubes, which further transform into complex shapes. Traditionally, analysis of dynamic tissues has relied on inspecting static snapshots, live imaging of cross-sections or tracking isolated cells in three dimensions. However, capturing the interplay between in-plane and out-of-plane behaviors requires following the full surface as it deforms and integrating cell-scale motions into collective, tissue-scale deformations. Here, we present an analysis framework that builds in toto maps of tissue deformations by following tissue parcels in a static material frame of reference. Our approach then relates in-plane and out-of-plane behaviors and decomposes complex deformation maps into elementary contributions. The tube-like surface Lagrangian analysis resource (TubULAR) provides an open-source implementation accessible either as a standalone toolkit or as an extension of the ImSAnE package used in the developmental biology community. We demonstrate our approach by analyzing shape change in the embryonic Drosophila midgut and beating zebrafish heart. The method naturally generalizes to in vitro and synthetic systems and provides ready access to the mechanical mechanisms relating genetic patterning to organ shape change.

Shape transitions in a network model of active elastic shells

Morphogenesis involves the transformation of initially simple shapes, such as multicellular spheroids, into more complex 3D shapes. These shape changes are governed by mechanical forces including molecular motor-generated forces as well as hydrostatic fluid pressure, both of which are actively regulated in living matter through mechano-chemical feedback. Inspired by autonomous, biophysical shape change, such as occurring in the model organism hydra, we introduce a minimal, active, elastic model featuring a network of springs in a globe-like spherical shell geometry. In this model there is coupling between activity and the shape of the shell: if the local curvature of a filament represented by a spring falls below a critical value, its elastic constant is actively changed. This results in deformation of the springs that changes the shape of the shell. By combining excitation of springs and pressure regulation, we show that the shell undergoes a transition from spheroidal to either elongated ellipsoidal or a different spheroidal shape, depending on pressure. There exists a critical pressure at which there is an abrupt change from ellipsoids to spheroids, showing that pressure is potentially a sensitive switch for material shape. We thus offer biologically inspired design principles for autonomous shape transitions in active elastic shells.

Diffusion coefficients and MSD measurements on curved membranes and porous media

We study some geometric aspects that influence the transport properties of particles that diffuse on curved surfaces. We compare different approaches to surface diffusion based on the Laplace-Beltrami operator adapted to predict concentration along entire membranes, confined subdomains along surfaces, or within porous media. Our goal is to summarize, firstly, how diffusion in these systems results in different types of diffusion coefficients and mean square displacement measurements, and secondly, how these two factors are affected by the concavity of the surface, the shape of the possible barriers or obstacles that form the available domains, the sinuosity, tortuosity, and constrictions of the trajectories and even how the observation plane affects the measurements of the diffusion. In addition to presenting a critical and organized comparison between different notions of MSD, in this review, we test the correspondence between theoretical predictions and numerical simulations by performing finite element simulations and illustrate some situations where diffusion theory can be applied. We briefly reviewed computational schemes for understanding surface diffusion and finally, discussed how this work contributes to understanding the role of surface diffusion transport properties in porous media and their relationship to other transport processes.

Wildebeest herds on rolling hills: Flocking on arbitrary curved surfaces

The collective behavior of active agents, whether herds of wildebeest or microscopic actin filaments propelled by molecular motors, is an exciting frontier in biological and soft matter physics. Almost three decades ago, Toner and Tu developed a continuum theory of the collective action of flocks, or herds, that helped launch the modern field of active matter. One challenge faced when applying continuum active matter theories to living phenomena is the complex geometric structure of biological environments. Both macroscopic and microscopic herds move on asymmetric curved surfaces, like undulating grass plains or the surface layers of cells or embryos, which can render problems analytically intractable. In this paper, we present a formulation of the Toner-Tu flocking theory that uses the finite element method to solve the governing equations on arbitrary curved surfaces. First, we test the developed formalism and its numerical implementation in channel flow with scattering obstacles and on cylindrical and spherical surfaces, comparing our results to analytical solutions. We then progress to surfaces with arbitrary curvature, moving beyond previously accessible problems to explore herding behavior on a variety of landscapes. This approach allows the investigation of transients and dynamic solutions not revealed by analytic methods. It also enables versatile incorporation of new geometries and boundary conditions and efficient sweeps of parameter space. Looking forward, the paper presented here lays the groundwork for a dialogue between Toner-Tu theory and data on collective motion in biologically relevant geometries, from drone footage of migrating animal herds to movies of microscopic cytoskeletal flows within cells.

Modelling contractile ring formation and division to daughter cells for simulating proliferative multicellular dynamics

Cell proliferation is a fundamental process underlying embryogenesis, homeostasis, wound healing, and cancer. The process involves multiple events during each cell cycle, such as cell growth, contractile ring formation, and division to daughter cells, which affect the surrounding cell population geometrically and mechanically. However, existing methods do not comprehensively describe the dynamics of multicellular structures involving cell proliferation at a subcellular resolution. In this study, we present a novel model for proliferative multicellular dynamics at the subcellular level by building upon the nonconservative fluid membrane (NCF) model that we developed in earlier research. The NCF model utilizes a dynamically-rearranging closed triangular mesh to depict the shape of each cell, enabling us to analyze cell dynamics over extended periods beyond each cell cycle, during which cell surface components undergo dynamic turnover. The proposed model represents the process of cell proliferation by incorporating cell volume growth and contractile ring formation through an energy function and topologically dividing each cell at the cleavage furrow formed by the ring. Numerical simulations demonstrated that the model recapitulated the process of cell proliferation at subcellular resolution, including cell volume growth, cleavage furrow formation, and division to daughter cells. Further analyses suggested that the orientation of actomyosin stress in the contractile ring plays a crucial role in the cleavage furrow formation, i.e., circumferential orientation can form a cleavage furrow but isotropic orientation cannot. Furthermore, the model replicated tissue-scale multicellular dynamics, where the successive proliferation of adhesive cells led to the formation of a cell sheet and stratification on the substrate. Overall, the proposed model provides a basis for analyzing proliferative multicellular dynamics at subcellular resolution.

Geometry-induced patterns through mechanochemical coupling

Intracellular protein patterns regulate a variety of vital cellular processes such as cell division and motility, which often involve dynamic cell-shape changes. These changes in cell shape may in turn affect the dynamics of pattern-forming proteins, hence leading to an intricate feedback loop between cell shape and chemical dynamics. While several computational studies have examined the rich resulting dynamics, the underlying mechanisms are not yet fully understood. To elucidate some of these mechanisms, we explore a conceptual model for cell polarity on a dynamic one-dimensional manifold. Using concepts from differential geometry, we derive the equations governing mass-conserving reaction-diffusion systems on time-evolving manifolds. Analyzing these equations mathematically, we show that dynamic shape changes of the membrane can induce pattern-forming instabilities in parts of the membrane, which we refer to as regional instabilities. Deformations of the local membrane geometry can also (regionally) suppress pattern formation and spatially shift already existing patterns. We explain our findings by applying and generalizing the local equilibria theory of mass-conserving reaction-diffusion systems. This allows us to determine a simple onset criterion for geometry-induced pattern-forming instabilities, which is linked to the phase-space structure of the reaction-diffusion system. The feedback loop between membrane shape deformations and reaction-diffusion dynamics then leads to a surprisingly rich phenomenology of patterns, including oscillations, traveling waves, and standing waves, even if these patterns do not occur in systems with a fixed membrane shape. Our paper reveals that the local conformation of the membrane geometry acts as an important dynamical control parameter for pattern formation in mass-conserving reaction-diffusion systems.

Model of membrane deformations driven by a surface pH gradient

Many cellular organelles are membrane-bound structures with complex membrane composition and shape. Their shapes have been observed to depend on the metabolic state of the organelle and the mechanisms that couple biochemical pathways and membrane shape are still actively investigated. Here, we study a model coupling inhomogeneities in the lipid composition and membrane geometry via a generalized Helfrich free energy. We derive the resulting stress tensor, the Green's function for a tubular membrane, and compute the phase diagram of the induced deformations. We then apply this model to study the deformation of mitochondria cristae described as membrane tubes supporting a pH gradient at its surface. This gradient in turn controls the lipid composition of the membrane via the protonation or deprotonation of cardiolipins, which are acid-based lipids known to be crucial for mitochondria shape and functioning. Our model predicts the appearance of tube deformations resembling the observed shape changes of cristea when submitted to a proton gradient.

A cell membrane model that reproduces cortical flow-driven cell migration and collective movement

Many fundamental biological processes are dependent on cellular migration. Although the mechanical mechanisms of single-cell migration are relatively well understood, those underlying migration of multiple cells adhered to each other in a cluster, referred to as cluster migration, are poorly understood. A key reason for this knowledge gap is that many forces-including contraction forces from actomyosin networks, hydrostatic pressure from the cytosol, frictional forces from the substrate, and forces from adjacent cells-contribute to cell cluster movement, making it challenging to model, and ultimately elucidate, the final result of these forces. This paper describes a two-dimensional cell membrane model that represents cells on a substrate with polygons and expresses various mechanical forces on the cell surface, keeping these forces balanced at all times by neglecting cell inertia. The model is discrete but equivalent to a continuous model if appropriate replacement rules for cell surface segments are chosen. When cells are given a polarity, expressed by a direction-dependent surface tension reflecting the location dependence of contraction and adhesion on a cell boundary, the cell surface begins to flow from front to rear as a result of force balance. This flow produces unidirectional cell movement, not only for a single cell but also for multiple cells in a cluster, with migration speeds that coincide with analytical results from a continuous model. Further, if the direction of cell polarity is tilted with respect to the cluster center, surface flow induces cell cluster rotation. The reason why this model moves while keeping force balance on cell surface (i.e., under no net forces from outside) is because of the implicit inflow and outflow of cell surface components through the inside of the cell. An analytical formula connecting cell migration speed and turnover rate of cell surface components is presented.

Patterning of morphogenetic anisotropy fields

Orientational order, encoded in anisotropic fields, plays an important role during the development of an organism. A striking example of this is the freshwater polyp Hydra, where topological defects in the muscle fiber orientation have been shown to localize to key features of the body plan. This body plan is organized by morphogen concentration gradients, raising the question how muscle fiber orientation, morphogen gradients and body shape interact. Here, we introduce a minimal model that couples nematic orientational order to the gradient of a morphogen field. We show that on a planar surface, alignment to a radial concentration gradient can induce unbinding of topological defects, as observed during budding and tentacle formation in Hydra, and stabilize aster/vortex-like defects, as observed at a Hydra's mouth. On curved surfaces mimicking the morphologies of Hydra in various stages of development-from spheroid to adult-our model reproduces the experimentally observed reorganization of orientational order. Our results suggest how gradient alignment and curvature effects may work together to control orientational order during development and lay the foundations for future modeling efforts that will include the tissue mechanics that drive shape deformations.

Curvature in Biological Systems: Its Quantification, Emergence, and Implications across the Scales

Surface curvature both emerges from, and influences the behavior of, living objects at length scales ranging from cell membranes to single cells to tissues and organs. The relevance of surface curvature in biology is supported by numerous experimental and theoretical investigations in recent years. In this review, first, a brief introduction to the key ideas of surface curvature in the context of biological systems is given and the challenges that arise when measuring surface curvature are discussed. Giving an overview of the emergence of curvature in biological systems, its significance at different length scales becomes apparent. On the other hand, summarizing current findings also shows that both single cells and entire cell sheets, tissues or organisms respond to curvature by modulating their shape and their migration behavior. Finally, the interplay between the distribution of morphogens or micro-organisms and the emergence of curvature across length scales is addressed with examples demonstrating these key mechanistic principles of morphogenesis. Overall, this review highlights that curved interfaces are not merely a passive by-product of the chemical, biological, and mechanical processes but that curvature acts also as a signal that co-determines these processes.

Homeotic compartment curvature and tension control spatiotemporal folding dynamics

Shape is a conspicuous and fundamental property of biological systems entailing the function of organs and tissues. While much emphasis has been put on how tissue tension and mechanical properties drive shape changes, whether and how a given tissue geometry influences subsequent morphogenesis remains poorly characterized. Here, we explored how curvature, a key descriptor of tissue geometry, impinges on the dynamics of epithelial tissue invagination. We found that the morphogenesis of the fold separating the adult Drosophila head and thorax segments is driven by the invagination of the Deformed (Dfd) homeotic compartment. Dfd controls invagination by modulating actomyosin organization and in-plane epithelial tension via the Tollo and Dystroglycan receptors. By experimentally introducing curvature heterogeneity within the homeotic compartment, we established that a curved tissue geometry converts the Dfd-dependent in-plane tension into an inward force driving folding. Accordingly, the interplay between in-plane tension and tissue curvature quantitatively explains the spatiotemporal folding dynamics. Collectively, our work highlights how genetic patterning and tissue geometry provide a simple design principle driving folding morphogenesis during development.

Active morphogenesis of patterned epithelial shells

Shape transformations of epithelial tissues in three dimensions, which are crucial for embryonic development or in vitro organoid growth, can result from active forces generated within the cytoskeleton of the epithelial cells. How the interplay of local differential tensions with tissue geometry and with external forces results in tissue-scale morphogenesis remains an open question. Here, we describe epithelial sheets as active viscoelastic surfaces and study their deformation under patterned internal tensions and bending moments. In addition to isotropic effects, we take into account nematic alignment in the plane of the tissue, which gives rise to shape-dependent, anisotropic active tensions and bending moments. We present phase diagrams of the mechanical equilibrium shapes of pre-patterned closed shells and explore their dynamical deformations. Our results show that a combination of nematic alignment and gradients in internal tensions and bending moments is sufficient to reproduce basic building blocks of epithelial morphogenesis, including fold formation, budding, neck formation, flattening, and tubulation.

Interacting active surfaces: A model for three-dimensional cell aggregates

We introduce a modelling and simulation framework for cell aggregates in three dimensions based on interacting active surfaces. Cell mechanics is captured by a physical description of the acto-myosin cortex that includes cortical flows, viscous forces, active tensions, and bending moments. Cells interact with each other via short-range forces capturing the effect of adhesion molecules. We discretise the model equations using a finite element method, and provide a parallel implementation in C++. We discuss examples of application of this framework to small and medium-sized aggregates: we consider the shape and dynamics of a cell doublet, a planar cell sheet, and a growing cell aggregate. This framework opens the door to the systematic exploration of the cell to tissue-scale mechanics of cell aggregates, which plays a key role in the morphogenesis of embryos and organoids.

Three-dimensional chiral morphodynamics of chemomechanical active shells

Morphogenesis of active shells such as cells is a fundamental chemomechanical process that often exhibits three-dimensional (3D) large deformations and chemical pattern dynamics simultaneously. Here, we establish a chemomechanical active shell theory accounting for mechanical feedback and biochemical regulation to investigate the symmetry-breaking and 3D chiral morphodynamics emerging in the cell cortex. The active bending and stretching of the elastic shells are regulated by biochemical signals like actomyosin and RhoA, which, in turn, exert mechanical feedback on the biochemical events via deformation-dependent diffusion and inhibition. We show that active deformations can trigger chemomechanical bifurcations, yielding pulse spiral waves and global oscillations, which, with increasing mechanical feedback, give way to traveling or standing waves subsequently. Mechanical feedback is also found to contribute to stabilizing the polarity of emerging patterns, thus ensuring robust morphogenesis. Our results reproduce and unravel the experimentally observed solitary and multiple spiral patterns, which initiate asymmetric cleavage in Xenopus and starfish embryogenesis. This study underscores the crucial roles of mechanical feedback in cell development and also suggests a chemomechanical framework allowing for 3D large deformation and chemical signaling to explore complex morphogenesis in living shell-like structures.

Self-organization in amoeboid motility

Amoeboid motility has come to refer to a spectrum of cell migration modes enabling a cell to move in the absence of strong, specific adhesion. To do so, cells have evolved a range of motile surface movements whose physical principles are now coming into view. In response to external cues, many cells—and some single-celled-organisms—have the capacity to turn off their default migration mode. and switch to an amoeboid mode. This implies a restructuring of the migration machinery at the cell scale and suggests a close link between cell polarization and migration mediated by self-organizing mechanisms. Here, I review recent theoretical models with the aim of providing an integrative, physical picture of amoeboid migration.

Active Nematic Defects and Epithelial Morphogenesis

Inspired by recent experiments that highlight the role of nematic defects in the morphogenesis of epithelial tissues, we develop a minimal framework to study the dynamics of an active curved surface driven by its nematic texture. Allowing the surface to evolve via relaxational dynamics leads to a theory linking nematic defect dynamics, cellular division rates, and Gaussian curvature. Regions of large positive (negative) curvature and positive (negative) growth are colocalized with the presence of positive (negative) defects. In an ex-vivo setting of cultured murine neural progenitor cells, we show that our framework is consistent with the observed cell accumulation at positive defects and depletion at negative defects. In an in-vivo setting, we show that the defect configuration consisting of a bound +1 defect state, which is stabilized by activity, surrounded by two -1/2 defects can create a stationary ring configuration of tentacles, consistent with observations of a basal marine invertebrate Hydra.

Theory of defect-mediated morphogenesis

Growing experimental evidence indicates that topological defects could serve as organizing centers in the morphogenesis of tissues. Here, we provide a quantitative explanation for this phenomenon, rooted in the buckling theory of deformable active polar liquid crystals. Using a combination of linear stability analysis and computational fluid dynamics, we demonstrate that active layers, such as confined cell monolayers, are unstable to the formation of protrusions in the presence of disclinations. The instability originates from an interplay between the focusing of the elastic forces, mediated by defects, and the renormalization of the system's surface tension by the active flow. The posttransitional regime is also characterized by several complex morphodynamical processes, such as oscillatory deformations, droplet nucleation, and active turbulence. Our findings offer an explanation of recent observations on tissue morphogenesis and shed light on the dynamics of active surfaces in general.

Morphogenesis in Synthetic Chemical Cells

We present an efficient model for describing morphogenesis and the emergence of spatiotemporal structures in synthetic chemical cells. This work is motivated by an experimental setup used for testing Turing's theory of morphogenesis. The model developed is based on the general theory of chemically active droplets, which combines the classical theory of phase separation with reaction-diffusion systems. Through the 2D calculations, we find the six spatiotemporal structures predicted by Turing in 1952 and experimentally observed, in a 1D array of droplets. Moreover, under Turing instability, with a determined chemical wavelength, the system undergoes morphogenesis. This theoretical approach provides a useful tool for understanding the physical differentiation through the direct calculation of the osmotic pressure in each cell as the chemical reaction occurs.

Pattern formation, localized and running pulsation on active spherical membranes

Active force generation by an actin-myosin cortex coupled to a cell membrane allows the cell to deform, respond to the environment, and mediate cell motility and division. Several membrane-bound activator proteins move along it and couple to the membrane curvature. Besides, they can act as nucleating sites for the growth of filamentous actin. Actin polymerization can generate a local outward push on the membrane. Inward pull from the contractile actomyosin cortex can propagate along the membrane via actin filaments. We use coupled evolution of fields to perform linear stability analysis and numerical calculations. As activity overcomes the stabilizing factors such as surface tension and bending rigidity, the spherical membrane shows instability towards pattern formation, localized pulsation, and running pulsation between poles. We present our results in terms of phase diagrams and evolutions of the coupled fields. They have relevance for living cells and can be verified in experiments on artificial cell-like constructs.

Forced into shape: Mechanical forces in Drosophila development and homeostasis

Mechanical forces play a central role in shaping tissues during development and maintaining epithelial integrity in homeostasis. In this review, we discuss the roles of mechanical forces in Drosophila development and homeostasis, starting from the interplay of mechanics with cell growth and division. We then discuss several examples of morphogenetic processes where complex 3D structures are shaped by mechanical forces, followed by a closer look at patterning processes. We also review the role of forces in homeostatic processes, including cell elimination and wound healing. Finally, we look at the interplay of mechanics and developmental robustness and discuss open questions in the field, as well as novel approaches that will help tackle them in the future.

A C++ expression system for partial differential equations enables generic simulations of biological hydrodynamics

We present a user-friendly and intuitive C++ expression system to implement numerical simulations of continuum biological hydrodynamics. The expression system allows writing simulation programs in near-mathematical notation and makes codes more readable, more compact, and less error-prone. It also cleanly separates the implementation of the partial differential equation model from the implementation of the numerical methods used to discretize it. This allows changing either of them with minimal changes to the source code. The presented expression system is implemented in the high-performance computing platform OpenFPM, supporting simulations that transparently parallelize on multi-processor computer systems. We demonstrate that our expression system makes it easier to write scalable codes for simulating biological hydrodynamics in space and time. We showcase the present framework in numerical simulations of active polar fluids, as well as in classic simulations of fluid dynamics from the incompressible Navier-Stokes equations to Stokes flow in a ball. The presented expression system accelerates scalable simulations of spatio-temporal models that encode the physics and material properties of tissues in order to algorithmically study morphogenesis.

Chemically controlled pattern formation in self-oscillating elastic shells

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

Active flows and deformable surfaces in development

We review progress in active hydrodynamic descriptions of flowing media on curved and deformable manifolds: the state-of-the-art in continuum descriptions of single-layers of epithelial and/or other tissues during development. First, after a brief overview of activity, flows and hydrodynamic descriptions, we highlight the generic challenge of identifying the dependence on dynamical variables of so-called active kinetic coefficients- active counterparts to dissipative Onsager coefficients. We go on to describe some of the subtleties concerning how curvature and active flows interact, and the issues that arise when surfaces are deformable. We finish with a broad discussion around the utility of such theories in developmental biology. This includes limitations to analytical techniques, challenges associated with numerical integration, fitting-to-data and inference, and potential tools for the future, such as discrete differential geometry.

Engineering life in synthetic systems

The second EMBO-EMBL Symposium 'Synthetic Morphogenesis: From Gene Circuits to Tissue Architecture' was held virtually in March 2021, with participants from all over the world joining from the comfort of their sofas to discuss synthetic morphogenesis at large. Leading scientists from a range of disciplines, including developmental biology, physics, chemistry and computer science, covered a gamut of topics from the principles of cell and tissue organization, patterning and gene regulatory networks, to synthetic approaches for exploring evolutionary and developmental biology principles. Here, we describe some of the high points.

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.

Modeling mechanochemical pattern formation in elastic sheets of biological matter

Inspired by active shape morphing in developing tissues and biomaterials, we investigate two generic mechanochemical models where the deformations of a thin elastic sheet are driven by, and in turn affect, the concentration gradients of a chemical signal. We develop numerical methods to study the coupled elastic deformations and chemical concentration kinetics, and illustrate with computations the formation of different patterns depending on shell thickness, strength of mechanochemical coupling and diffusivity. In the first model, the sheet curvature governs the production of a contractility inhibitor and depending on the threshold in the coupling, qualitatively different patterns occur. The second model is based on the stress-dependent activity of myosin motors and demonstrates how the concentration distribution patterns of molecular motors are affected by the long-range deformations generated by them. Since the propagation of mechanical deformations is typically faster than chemical kinetics (of molecular motors or signaling agents that affect motors), we describe in detail and implement a numerical method based on separation of timescales to effectively simulate such systems. We show that mechanochemical coupling leads to long-range propagation of patterns in disparate systems through elastic instabilities even without the diffusive or advective transport of the chemicals.

Mechanochemical induction of wrinkling morphogenesis on elastic shells

Morphogenetic dynamics of tissue sheets require coordinated cell shape changes regulated by global patterning of mechanical forces. Inspired by such biological phenomena, we propose a minimal mechanochemical model based on the notion that cell shape changes are induced by diffusible biomolecules that influence tissue contractility in a concentration-dependent manner - and whose concentration is in turn affected by the macroscopic tissue shape. We perform computational simulations of thin shell elastic dynamics to reveal propagating chemical and three-dimensional deformation patterns arising due to a sequence of buckling instabilities. Depending on the concentration threshold that actuates cell shape change, we find qualitatively different patterns. The mechanochemically coupled patterning dynamics are distinct from those driven by purely mechanical or purely chemical factors, and emerge even without diffusion. Using numerical simulations and theoretical arguments, we analyze the elastic instabilities that result from our model and provide simple scaling laws to identify wrinkling morphologies.

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

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

Continuum Theory of Active Phase Separation in Cellular Aggregates

Dense cellular aggregates are common in biology, ranging from bacterial biofilms to organoids, cell spheroids, and tumors. Their dynamics, driven by intercellular forces, is intrinsically out of equilibrium. Motivated by bacterial colonies as a model system, we present a continuum theory to study dense, active, cellular aggregates. We describe the process of aggregate formation as an active phase separation phenomenon, while the merging of aggregates is rationalized as a coalescence of viscoelastic droplets where the key timescales are linked to the turnover of the active force. Our theory provides a general framework for studying the rheology and nonequilibrium dynamics of dense cellular aggregates.

Tools for computational analysis of moving boundary problems in cellular mechanobiology

A cell's ability to change shape is one of the most fundamental biological processes and is essential for maintaining healthy organisms. When the ability to control shape goes awry, it often results in a diseased system. As such, it is important to understand the mechanisms that allow a cell to sense and respond to its environment so as to maintain cellular shape homeostasis. Because of the inherent complexity of the system, computational models that are based on sound theoretical understanding of the biochemistry and biomechanics and that use experimentally measured parameters are an essential tool. These models involve an inherent feedback, whereby shape is determined by the action of regulatory signals whose spatial distribution depends on the shape. To carry out computational simulations of these moving boundary problems requires special computational techniques. A variety of alternative approaches, depending on the type and scale of question being asked, have been used to simulate various biological processes, including cell motility, division, mechanosensation, and cell engulfment. In general, these models consider the forces that act on the system (both internally generated, or externally imposed) and the mechanical properties of the cell that resist these forces. Moving forward, making these techniques more accessible to the non-expert will help improve interdisciplinary research thereby providing new insight into important biological processes that affect human health. This article is categorized under: Cancer > Cancer>Computational Models Cancer > Cancer>Molecular and Cellular Physiology.

Universal Thinning of Liquid Filaments under Dominant Surface Forces

Theory and numerical simulations of the thinning of liquid threads at high superficial concentration of surfactants suggest the existence of an asymptotic regime where surface tension balances surface viscous stresses, leading to an exponential thinning with an e-fold time F(Θ)(3μ_{s}+κ_{s})/σ, where μ_{s} and κ_{s} are the surface shear and dilatational viscosity coefficients, σ is the interfacial tension, Θ=κ_{s}/μ_{s}, and F(Θ) is a universal function. The potential use of this phenomenon to measure the surface viscosity coefficients is discussed.

Cell lineage-dependent chiral actomyosin flows drive cellular rearrangements in early Caenorhabditis elegans development

Proper positioning of cells is essential for many aspects of development. Daughter cell positions can be specified via orienting the cell division axis during cytokinesis. Rotatory actomyosin flows during division have been implied in specifying and reorienting the cell division axis, but how general such reorientation events are, and how they are controlled, remains unclear. We followed the first nine divisions of Caenorhabditis elegans embryo development and demonstrate that chiral counter-rotating flows arise systematically in early AB lineage, but not in early P/EMS lineage cell divisions. Combining our experiments with thin film active chiral fluid theory we identify a mechanism by which chiral counter-rotating actomyosin flows arise in the AB lineage only, and show that they drive lineage-specific spindle skew and cell reorientation events. In conclusion, our work sheds light on the physical processes that underlie chiral morphogenesis in early development.

Shape and size changes of adherent elastic epithelia

Epithelial tissues play a fundamental role in various morphogenetic events during development and early embryogenesis. Although epithelial monolayers are often modeled as two-dimensional (2D) elastic surfaces, they distinguish themselves from conventional thin elastic plates in three important ways- the presence of an apical-basal polarity, spatial variability of cellular thickness, and their nonequilibrium active nature. Here, we develop a minimal continuum model of a planar epithelial tissue as an active elastic material that incorporates all these features. We start from a full three-dimensional (3D) description of the tissue and derive an effective 2D model that captures, through the curvature of the apical surface, both the apical-basal asymmetry and the spatial geometry of the tissue. Crucially, variations of active stresses across the apical-basal axis lead to active torques that can drive curvature transitions. By identifying four distinct sources of activity, we find that bulk active stresses arising from actomyosin contractility and growth compete with boundary active tensions due to localized actomyosin cables and lamellipodial activity to generate the various states spanning the morphospace of a planar epithelium. Our treatment hence unifies 3D shape deformations through the coupled mechanics of apical curvature change and in-plane expansion/contraction of substrate-adhered tissues. Finally, we discuss the implications of our results for some biologically relevant processes such as tissue folding at the onset of lumen formation.

Gait-optimized locomotion of wave-driven soft sheets

Inspired by the robust locomotion of limbless animals in a range of environments, the development of soft robots capable of moving by localized swelling, bending, and other forms of differential growth has become a target for soft matter research over the last decade. Engineered soft robots exhibit a wide range of morphologies, but theoretical investigations of soft robot locomotion have largely been limited to slender bodied or one-dimensional examples. Here, we demonstrate design principles regarding the locomotion of two-dimensional soft materials driven by morphoelastic waves along a dry substrate. Focusing on the essential common aspects of many natural and man-made soft actuators, a continuum model is developed which links the deformation of a thin elastic sheet to surface-bound excitation waves. Through a combination of analytic and numerical methods, we investigate the relationship between induced active stress and self-propulsion performance of self-propelling sheets driven by FitzHugh-Nagumo type chemical waves. Examining the role of both sheet geometry and terrain geography on locomotion, our results can provide guidance for the design of more efficient soft crawling devices.

The Poisson Ratio of the Cellular Actin Cortex Is Frequency Dependent

Cell shape changes are vital for many physiological processes such as cell proliferation, cell migration, and morphogenesis. They emerge from an orchestrated interplay of active cellular force generation and passive cellular force response, both crucially influenced by the actin cytoskeleton. To model cellular force response and deformation, cell mechanical models commonly describe the actin cytoskeleton as a contractile isotropic incompressible material. However, in particular at slow frequencies, there is no compelling reason to assume incompressibility because the water content of the cytoskeleton may change. Here, we challenge the assumption of incompressibility by comparing computer simulations of an isotropic actin cortex with tunable Poisson ratio to measured cellular force response. Comparing simulation results and experimental data, we determine the Poisson ratio of the cortex in a frequency-dependent manner. We find that the Poisson ratio of the cortex decreases in the measured frequency regime analogous to trends reported for the Poisson ratio of glassy materials. Our results therefore indicate that actin cortex compression or dilation is possible in response to acting forces at sufficiently fast timescales. This finding has important implications for the parameterization in active gel theories that describe actin cytoskeletal dynamics.

Topology and Morphology of Self-Deforming Active Shells

We present a generic framework for modeling three-dimensional deformable shells of active matter that captures the orientational dynamics of the active particles and hydrodynamic interactions on the shell and with the surrounding environment. We find that the cross talk between the self-induced flows of active particles and dynamic reshaping of the shell can result in conformations that are tunable by varying the form and magnitude of active stresses. We further demonstrate and explain how self-induced topological defects in the active layer can direct the morphodynamics of the shell. These findings are relevant to understanding morphological changes during organ development and the design of bioinspired materials that are capable of self-organization.

Minimal Model of Cellular Symmetry Breaking

The cell cortex, a thin film of active material assembled below the cell membrane, plays a key role in cellular symmetry-breaking processes such as cell polarity establishment and cell division. Here, we present a minimal model of the self-organization of the cell cortex that is based on a hydrodynamic theory of curved active surfaces. Active stresses on this surface are regulated by a diffusing molecular species. We show that coupling of the active surface to a passive bulk fluid enables spontaneous polarization and the formation of a contractile ring on the surface via mechanochemical instabilities. We discuss the role of external fields in guiding such pattern formation. Our work reveals that key features of cellular symmetry breaking and cell division can emerge in a minimal model via general dynamic instabilities.

Probing the Functional Role of Physical Motion in Development

Spatiotemporal organization during development has frequently been proposed to be explainable by reaction-transport models, where biochemical reactions couple to physical motion. However, whereas genetic tools allow causality of molecular players to be dissected via perturbation experiments, the functional role of physical transport processes, such as diffusion and cytoplasmic streaming, frequently remains untestable. This Perspective explores the challenges of validating reaction-transport hypotheses and highlights new opportunities provided by perturbation approaches that specifically target physical transport mechanisms. Using these methods, experimental physics may begin to catch up with molecular biology and find ways to test roles of diffusion and flows in development.

Shape Transformations of Vesicles Induced by Swim Pressure

While the behavior of vesicles in thermodynamic equilibrium has been studied extensively, how active forces control vesicle shape transformations is not understood. Here, we combine theory and simulations to study the shape behavior of vesicles containing active Brownian particles. We show that the combination of active forces, dimensionality, and membrane bending free energy creates a plethora of novel phase transitions. At low swim pressure, the vesicle exhibits a discontinuous transition from a spherical to a prolate shape, which has no counterpart in two dimensions. At high swim pressure it exhibits stochastic spatiotemporal oscillations. Our work helps researchers to understand and control the shape dynamics of membranes in active-matter systems.