A Geometric Tension Dynamics Model of Epithelial Convergent Extension.

Epithelial tissue elongation by convergent extension is a key motif of animal morphogenesis. On a coarse scale, cell motion resembles laminar fluid flow; yet in contrast to a fluid, epithelial cells adhere to each other and maintain the tissue layer under actively generated internal tension. To resolve this apparent paradox, we formulate a model in which tissue flow occurs through adiabatic remodelling of the cellular force balance causing local cell rearrangement. We propose that the gradual shifting of the force balance is caused by positive feedback on myosin-generated cytoskeletal tension. Shifting force balance within a tension network causes active T1s oriented by the global anisotropy of tension. Rigidity of cells against shape changes converts the oriented internal rearrangements into net tissue deformation. Strikingly, we find that the total amount of tissue extension depends on the initial magnitude of anisotropy and on cellular packing order. T1s degrade this order so that tissue flow is self-limiting. We explain these findings by showing that coordination of T1s depends on coherence in local tension configurations, quantified by a certain order parameter in tension space. Our model reproduces the salient tissue- and cell-scale features of germ band elongation during Drosophila gastrulation, in particular the slowdown of tissue flow after approximately twofold extension concomitant with a loss of order in tension configurations. This suggests local cell geometry contains morphogenetic information and yields predictions testable in future experiments. Furthermore, our focus on defining biologically controlled active tension dynamics on the manifold of force-balanced states may provide a general approach to the description of morphogenetic flow.

Redundancy and the role of protein copy numbers in the cell polarization machinery of budding yeast.

How can a self-organized cellular function evolve, adapt to perturbations, and acquire new sub-functions? To make progress in answering these basic questions of evolutionary cell biology, we analyze, as a concrete example, the cell polarity machinery of Saccharomyces cerevisiae. This cellular module exhibits an intriguing resilience: it remains operational under genetic perturbations and recovers quickly and reproducibly from the deletion of one of its key components. Using a combination of modeling, conceptual theory, and experiments, we propose that multiple, redundant self-organization mechanisms coexist within the protein network underlying cell polarization and are responsible for the module's resilience and adaptability. Based on our mechanistic understanding of polarity establishment, we hypothesize that scaffold proteins, by introducing new connections in the existing network, can increase the redundancy of mechanisms and thus increase the evolvability of other network components. Moreover, our work gives a perspective on how a complex, redundant cellular module might have evolved from a more rudimental ancestral form.

The Geometric Basis of Epithelial Convergent Extension.

Shape changes of epithelia during animal development, such as convergent extension, are achieved through concerted mechanical activity of individual cells. While much is known about the corresponding large scale tissue flow and its genetic drivers, key open questions regard the cell-scale mechanics, e.g. internal vs external driving forces, and coordination, e.g. bottom-up self-organization vs top-down genetic instruction. To address these questions, we develop a quantitative, model-based analysis framework to relate cell geometry to local tension in recently obtained timelapse imaging data of gastrulating Drosophila embryos. This analysis provides a systematic decomposition of cell shape changes and T1-rearrangements into internally driven, active, and externally driven, passive, contributions. Specifically, we find evidence that germ band extension is driven by active T1 processes that self-organize through positive feedback acting on tensions. More generally, our findings suggest that epithelial convergent extension results from controlled transformation of internal force balance geometry which we quantify with a novel quantification tool for local tension configurations.

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.

Directing Min protein patterns with advective bulk flow.

The Min proteins constitute the best-studied model system for pattern formation in cell biology. We theoretically predict and experimentally show that the propagation direction of in vitro Min protein patterns can be controlled by a hydrodynamic flow of the bulk solution. We find downstream propagation of Min wave patterns for low MinE:MinD concentration ratios, upstream propagation for large ratios, but multistability of both propagation directions in between. Whereas downstream propagation can be described by a minimal model that disregards MinE conformational switching, upstream propagation can be reproduced by a reduced switch model, where increased MinD bulk concentrations on the upstream side promote protein attachment. Our study demonstrates that a differential flow, where bulk flow advects protein concentrations in the bulk, but not on the surface, can control surface-pattern propagation. This suggests that flow can be used to probe molecular features and to constrain mathematical models for pattern-forming systems.

Learning a conserved mechanism for early neuroectoderm morphogenesis.

Morphogenesis is the process whereby the body of an organism develops its target shape. The morphogen BMP is known to play a conserved role across bilaterian organisms in determining the dorsoventral (DV) axis. Yet, how BMP governs the spatio-temporal dynamics of cytoskeletal proteins driving morphogenetic flow remains an open question. Here, we use machine learning to mine a morphodynamic atlas of Drosophila development, and construct a mathematical model capable of predicting the coupled dynamics of myosin, E-cadherin, and morphogenetic flow. Mutant analysis shows that BMP sets the initial condition of this dynamical system according to the following signaling cascade: BMP establishes DV pair-rule-gene patterns that set-up an E-cadherin gradient which in turn creates a myosin gradient in the opposite direction through mechanochemical feedbacks. Using neural tube organoids, we argue that BMP, and the signaling cascade it triggers, prime the conserved dynamics of neuroectoderm morphogenesis from fly to humans.

Coarsening and wavelength selection far from equilibrium: A unifying framework based on singular perturbation theory.

Intracellular protein patterns are described by (nearly) mass-conserving reaction-diffusion systems. While these patterns initially form out of a homogeneous steady state due to the well-understood Turing instability, no general theory exists for the dynamics of fully nonlinear patterns. We develop a unifying theory for nonlinear wavelength-selection dynamics in (nearly) mass-conserving two-component reaction-diffusion systems independent of the specific mathematical model chosen. Previous work has shown that these systems support an extremely broad band of stable wavelengths, but the mechanism by which a specific wavelength is selected has remained unclear. We show that an interrupted coarsening process selects the wavelength at the threshold to stability. Based on the physical intuition that coarsening is driven by competition for mass and interrupted by weak source terms that break strict mass conservation, we develop a singular perturbation theory for the stability of stationary patterns. The resulting closed-form analytical expressions enable us to quantitatively predict the coarsening dynamics and the final pattern wavelength. We find excellent agreement with numerical results throughout the diffusion- and reaction-limited regimes of the dynamics, including the crossover region. Further, we show how, in these limits, the two-component reaction-diffusion systems map to generalized Cahn-Hilliard and conserved Allen-Cahn dynamics, therefore providing a link to these two fundamental scalar field theories. The systematic understanding of the length-scale dynamics of fully nonlinear patterns in two-component systems provided here builds the basis to reveal the mechanisms underlying wavelength selection in multicomponent systems with potentially several conservation laws.

Length Regulation Drives Self-Organization in Filament-Motor Mixtures.

Cytoskeletal networks form complex intracellular structures. Here we investigate a minimal model for filament-motor mixtures in which motors act as depolymerases and thereby regulate filament length. Combining agent-based simulations and hydrodynamic equations, we show that resource-limited length regulation drives the formation of filament clusters despite the absence of mechanical interactions between filaments. Even though the orientation of individual remains fixed, collective filament orientation emerges in the clusters, aligned orthogonal to their interfaces.

Bridging scales in a multiscale pattern-forming system.

Self-organized pattern formation is vital for many biological processes. Reaction-diffusion models have advanced our understanding of how biological systems develop spatial structures, starting from homogeneity. However, biological processes inherently involve multiple spatial and temporal scales and transition from one pattern to another over time, rather than progressing from homogeneity to a pattern. To deal with such multiscale systems, coarse-graining methods are needed that allow the dynamics to be reduced to the relevant degrees of freedom at large scales, but without losing information about the patterns at small scales. Here, we present a semiphenomenological approach which exploits mass conservation in pattern formation, and enables reconstruction of information about patterns from the large-scale dynamics. The basic idea is to partition the domain into distinct regions (coarse grain) and determine instantaneous dispersion relations in each region, which ultimately inform about local pattern-forming instabilities. We illustrate our approach by studying the Min system, a paradigmatic model for protein pattern formation. By performing simulations, we first show that the Min system produces multiscale patterns in a spatially heterogeneous geometry. This prediction is confirmed experimentally by in vitro reconstitution of the Min system. Using a recently developed theoretical framework for mass-conserving reaction-diffusion systems, we show that the spatiotemporal evolution of the total protein densities on large scales reliably predicts the pattern-forming dynamics. Our approach provides an alternative and versatile theoretical framework for complex systems where analytical coarse-graining methods are not applicable, and can, in principle, be applied to a wide range of systems with an underlying conservation law.

Bulk-surface coupling identifies the mechanistic connection between Min-protein patterns in vivo and in vitro.

Self-organisation of Min proteins is responsible for the spatial control of cell division in Escherichia coli, and has been studied both in vivo and in vitro. Intriguingly, the protein patterns observed in these settings differ qualitatively and quantitatively. This puzzling dichotomy has not been resolved to date. Using reconstituted proteins in laterally wide microchambers with a well-controlled height, we experimentally show that the Min protein dynamics on the membrane crucially depend on the micro chamber height due to bulk concentration gradients orthogonal to the membrane. A theoretical analysis shows that in vitro patterns at low microchamber height are driven by the same lateral oscillation mode as pole-to-pole oscillations in vivo. At larger microchamber height, additional vertical oscillation modes set in, marking the transition to a qualitatively different in vitro regime. Our work reveals the qualitatively different mechanisms of mass transport that govern Min protein-patterns for different bulk heights and thus shows that Min patterns in cells are governed by a different mechanism than those in vitro.

Wavelength Selection by Interrupted Coarsening in Reaction-Diffusion Systems.

Wavelength selection in reaction-diffusion systems can be understood as a coarsening process that is interrupted by counteracting processes at certain wavelengths. We first show that coarsening in mass-conserving systems is driven by self-amplifying mass transport between neighboring high-density domains. We derive a general coarsening criterion and show that coarsening is generically uninterrupted in two-component systems that conserve mass. The theory is then generalized to study interrupted coarsening and anticoarsening due to weakly broken mass conservation, providing a general path to analyze wavelength selection in pattern formation far from equilibrium.

Flow Induced Symmetry Breaking in a Conceptual Polarity Model.

Important cellular processes, such as cell motility and cell division, are coordinated by cell polarity, which is determined by the non-uniform distribution of certain proteins. Such protein patterns form via an interplay of protein reactions and protein transport. Since Turing's seminal work, the formation of protein patterns resulting from the interplay between reactions and diffusive transport has been widely studied. Over the last few years, increasing evidence shows that also advective transport, resulting from cytosolic and cortical flows, is present in many cells. However, it remains unclear how and whether these flows contribute to protein-pattern formation. To address this question, we use a minimal model that conserves the total protein mass to characterize the effects of cytosolic flow on pattern formation. Combining a linear stability analysis with numerical simulations, we find that membrane-bound protein patterns propagate against the direction of cytoplasmic flow with a speed that is maximal for intermediate flow speed. We show that the mechanism underlying this pattern propagation relies on a higher protein influx on the upstream side of the pattern compared to the downstream side. Furthermore, we find that cytosolic flow can change the membrane pattern qualitatively from a peak pattern to a mesa pattern. Finally, our study shows that a non-uniform flow profile can induce pattern formation by triggering a regional lateral instability.

Pattern localization to a domain edge.

The formation of protein patterns inside cells is generically described by reaction-diffusion models. The study of such systems goes back to Turing, who showed how patterns can emerge from a homogenous steady state when two reactive components have different diffusivities (e.g., membrane-bound and cytosolic states). However, in nature, systems typically develop in a heterogeneous environment, where upstream protein patterns affect the formation of protein patterns downstream. Examples for this are the polarization of Cdc42 adjacent to the previous bud site in budding yeast and the formation of an actin-recruiter ring that forms around a PIP3 domain in macropinocytosis. This suggests that previously established protein patterns can serve as a template for downstream proteins and that these downstream proteins can "sense" the edge of the template. A mechanism for how this edge sensing may work remains elusive. Here we demonstrate and analyze a generic and robust edge-sensing mechanism, based on a two-component mass-conserving reaction-diffusion (McRD) model. Our analysis is rooted in a recently developed theoretical framework for McRD systems, termed local equilibria theory. We extend this framework to capture the spatially heterogeneous reaction kinetics due to the template. This enables us to graphically construct the stationary patterns in the phase space of the reaction kinetics. Furthermore, we show that the protein template can trigger a regional mass-redistribution instability near the template edge, leading to the accumulation of protein mass, which eventually results in a stationary peak at the template edge. We show that simple geometric criteria on the reactive nullcline's shape predict when this edge-sensing mechanism is operational. Thus, our results provide guidance for future studies of biological systems and for the design of synthetic pattern forming systems.

Design of biochemical pattern forming systems from minimal motifs.

Although molecular self-organization and pattern formation are key features of life, only very few pattern-forming biochemical systems have been identified that can be reconstituted and studied in vitro under defined conditions. A systematic understanding of the underlying mechanisms is often hampered by multiple interactions, conformational flexibility and other complex features of the pattern forming proteins. Because of its compositional simplicity of only two proteins and a membrane, the MinDE system from Escherichia coli has in the past years been invaluable for deciphering the mechanisms of spatiotemporal self-organization in cells. Here, we explored the potential of reducing the complexity of this system even further, by identifying key functional motifs in the effector MinE that could be used to design pattern formation from scratch. In a combined approach of experiment and quantitative modeling, we show that starting from a minimal MinE-MinD interaction motif, pattern formation can be obtained by adding either dimerization or membrane-binding motifs. Moreover, we show that the pathways underlying pattern formation are recruitment-driven cytosolic cycling of MinE and recombination of membrane-bound MinE, and that these differ in their in vivo phenomenology.

Capping protein-controlled actin polymerization shapes lipid membranes.

Arp2/3 complex-mediated actin assembly at cell membranes drives the formation of protrusions or endocytic vesicles. To identify the mechanism by which different membrane deformations can be achieved, we reconstitute the basic membrane deformation modes of inward and outward bending in a confined geometry by encapsulating a minimal set of cytoskeletal proteins into giant unilamellar vesicles. Formation of membrane protrusions is favoured at low capping protein (CP) concentrations, whereas the formation of negatively bent domains is promoted at high CP concentrations. Addition of non-muscle myosin II results in full fission events in the vesicle system. The different deformation modes are rationalized by simulations of the underlying transient nature of the reaction kinetics. The relevance of the regulatory mechanism is supported by CP overexpression in mouse melanoma B16-F1 cells and therefore demonstrates the importance of the quantitative understanding of microscopic kinetic balances to address the diverse functionality of the cytoskeleton.