Dynamics of Cdc42 network embodies a Turing-type mechanism of yeast cell polarity

Versatile system cores as a conceptual basis for generality in cell and developmental biology

The discovery of general principles underlying the complexity and diversity of cellular and developmental systems is a central and long-standing aim of biology. While new technologies collect data at an ever-accelerating rate, there is growing concern that conceptual progress is not keeping pace. We contend that this is due to a paucity of conceptual frameworks that support meaningful generalizations. This led us to develop the core and periphery (C&P) hypothesis, which posits that many biological systems can be decomposed into a highly versatile core with a large behavioral repertoire and a specific periphery that configures said core to perform one particular function. Versatile cores tend to be widely reused across biology, which confers generality to theories describing them. Here, we introduce this concept and describe examples at multiple scales, including Turing patterning, actomyosin dynamics, multi-cellular morphogenesis, and vertebrate gastrulation. We also sketch its evolutionary basis and discuss key implications and open questions. We propose that the C&P hypothesis could unlock new avenues of conceptual progress in mesoscale biology.

Multiscale Modeling of Bistability in the Yeast Polarity Circuit

Cell polarity refers to the asymmetric distribution of proteins and other molecules along a specified axis within a cell. Polarity establishment is the first step in many cellular processes. For example, directed growth or migration requires the formation of a cell front and back. In many cases, polarity occurs in the absence of spatial cues. That is, the cell undergoes symmetry breaking. Understanding the molecular mechanisms that allow cells to break symmetry and polarize requires computational models that span multiple spatial and temporal scales. Here, we apply a multiscale modeling approach to examine the polarity circuit of yeast. In addition to symmetry breaking, experiments revealed two key features of the yeast polarity circuit: bistability and rapid dismantling of the polarity site following a loss of signal. We used modeling based on ordinary differential equations (ODEs) to investigate mechanisms that generate these behaviors. Our analysis revealed that a model involving positive and negative feedback acting on different time scales captured both features. We then extend our ODE model into a coarse-grained reaction-diffusion equation (RDE) model to capture the spatial profiles of polarity factors. After establishing that the coarse-grained RDE model qualitatively captures key features of the polarity circuit, we expand it to more accurately capture the biochemical reactions involved in the system. We convert the expanded model to a particle-based model that resolves individual molecules and captures fluctuations that arise from the stochastic nature of biochemical reactions. Our models assume that negative regulation results from negative feedback. However, experimental observations do not rule out the possibility that negative regulation occurs through an incoherent feedforward loop. Therefore, we conclude by using our RDE model to suggest how negative feedback might be distinguished from incoherent feedforward regulation.

Testing models of cell cortex wave generation by Rho GTPases

The Rho GTPases pattern the cell cortex in a variety of fundamental cell-morphogenetic processes including division, wound repair, and locomotion. It has recently become apparent that this patterning arises from the ability of the Rho GTPases to self-organize into static and migrating spots, contractile pulses, and propagating waves in cells from yeasts to mammals 1 . These self-organizing Rho GTPase patterns have been explained by a variety of theoretical models which require multiple interacting positive and negative feedback loops. However, it is often difficult, if not impossible, to discriminate between different models simply because the available experimental data do not simultaneously capture the dynamics of multiple molecular concentrations and biomechanical variables at fine spatial and temporal resolution. Specifically, most studies typically provide either the total Rho GTPase signal or the Rho GTPase activity as reported by various sensors, but not both. Therefore, it remains largely unknown how membrane accumulation of Rho GTPases (i.e., Rho membrane enrichment) is related to Rho activity. Here we dissect the dynamics of RhoA by simultaneously imaging both total RhoA and active RhoA in the regime of acute cortical excitability 2 , characterized by pronounced waves of Rho activity and F-actin polymerization 3-5 . We find that within nascent waves, accumulation of active RhoA precedes that of total RhoA, and we exploit this finding to distinguish between two popular theoretical models previously used to explain propagating cortical Rho waves.

Patterning of the cell cortex by Rho GTPases

The Rho GTPases - RHOA, RAC1 and CDC42 - are small GTP binding proteins that regulate basic biological processes such as cell locomotion, cell division and morphogenesis by promoting cytoskeleton-based changes in the cell cortex. This regulation results from active (GTP-bound) Rho GTPases stimulating target proteins that, in turn, promote actin assembly and myosin 2-based contraction to organize the cortex. This basic regulatory scheme, well supported by in vitro studies, led to the natural assumption that Rho GTPases function in vivo in an essentially linear matter, with a given process being initiated by GTPase activation and terminated by GTPase inactivation. However, a growing body of evidence based on live cell imaging, modelling and experimental manipulation indicates that Rho GTPase activation and inactivation are often tightly coupled in space and time via signalling circuits and networks based on positive and negative feedback. In this Review, we present and discuss this evidence, and we address one of the fundamental consequences of coupled activation and inactivation: the ability of the Rho GTPases to self-organize, that is, direct their own transition from states of low order to states of high order. We discuss how Rho GTPase self-organization results in the formation of diverse spatiotemporal cortical patterns such as static clusters, oscillatory pulses, travelling wave trains and ring-like waves. Finally, we discuss the advantages of Rho GTPase self-organization and pattern formation for cell function.

Rho of Plants patterning: linking mathematical models and molecular diversity

ROPs (Rho of Plants) are plant specific small GTPases involved in many membrane patterning processes and play important roles in the establishment and communication of cell polarity. These small GTPases can produce a wide variety of patterns, ranging from a single cluster in tip-growing root hairs and pollen tubes to an oriented stripe pattern controlling protoxylem cell wall deposition. For an understanding of what controls these various patterns, models are indispensable. Consequently, many modelling studies on small GTPase patterning exist, often focusing on yeast or animal cells. Multiple patterns occurring in plants, however, require the stable co-existence of multiple active ROP clusters, which does not occur with the most common yeast/animal models. The possibility of such patterns critically depends on the precise model formulation. Additionally, different small GTPases are usually treated interchangeably in models, even though plants possess two types of ROPs with distinct molecular properties, one of which is unique to plants. Furthermore, the shape and even the type of ROP patterns may be affected by the cortical cytoskeleton, and cortex composition and anisotropy differ dramatically between plants and animals. Here, we review insights into ROP patterning from modelling efforts across kingdoms, as well as some outstanding questions arising from these models and recent experimental findings.

Turing pattern formation on the sphere is robust to the removal of a hole

The formation of buds on the cell membrane of budding yeast cells is thought to be driven by reactions and diffusion involving the protein Cdc42. These processes can be described by a coupled system of partial differential equations known as the Schnakenberg system. The Schnakenberg system is known to exhibit diffusion-driven pattern formation, thus providing a mechanism for bud formation. However, it is not known how the accumulation of bud scars on the cell membrane affect the ability of the Schnakenberg system to form patterns. We have approached this problem by modelling a bud scar on the cell membrane with a hole on the sphere. We have studied how the spectrum of the Laplace-Beltrami operator, which determines the resulting pattern, is affected by the size of the hole, and by numerically solving the Schnakenberg system on a sphere with a hole using the finite element method. Both theoretical predictions and numerical solutions show that pattern formation is robust to the introduction of a bud scar of considerable size, which lends credence to the hypothesis that bud formation is driven by diffusion-driven instability.

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.

A simple reaction-diffusion system as a possible model for the origin of chemotaxis

Chemotaxis is a directed cell movement in response to external chemical stimuli. In this paper, we propose a simple model for the origin of chemotaxis - namely how a directed movement in response to an external chemical signal may occur based on purely reaction-diffusion equations reflecting inner working of the cells. The model is inspired by the well-studied role of the rho-GTPase Cdc42 regulator of cell polarity, in particular in yeast cells. We analyse several versions of the model to better understand its analytic properties and prove global regularity in one and two dimensions. Using computer simulations, we demonstrate that in the framework of this model, at least in certain parameter regimes, the speed of the directed movement appears to be proportional to the size of the gradient of signalling chemical. This coincides with the form of the chemical drift in the most studied mean field model of chemotaxis, the Keller-Segel equation.

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.

Particle-based simulations reveal two positive feedback loops allow relocation and stabilization of the polarity site during yeast mating

Many cells adjust the direction of polarized growth or migration in response to external directional cues. The yeast Saccharomyces cerevisiae orient their cell fronts (also called polarity sites) up pheromone gradients in the course of mating. However, the initial polarity site is often not oriented towards the eventual mating partner, and cells relocate the polarity site in an indecisive manner before developing a stable orientation. During this reorientation phase, the polarity site displays erratic assembly-disassembly behavior and moves around the cell cortex. The mechanisms underlying this dynamic behavior remain poorly understood. Particle-based simulations of the core polarity circuit revealed that molecular-level fluctuations are unlikely to overcome the strong positive feedback required for polarization and generate relocating polarity sites. Surprisingly, inclusion of a second pathway that promotes polarity site orientation generated relocating polarity sites with properties similar to those observed experimentally. This pathway forms a second positive feedback loop involving the recruitment of receptors to the cell membrane and couples polarity establishment to gradient sensing. This second positive feedback loop also allows cells to stabilize their polarity site once the site is aligned with the pheromone gradient.

Membrane nanodomains: Dynamic nanobuilding blocks of polarized cell growth

Cell polarity is intimately linked to numerous biological processes, such as oriented plant cell division, particular asymmetric division, cell differentiation, cell and tissue morphogenesis, and transport of hormones and nutrients. Cell polarity is typically initiated by a polarizing cue that regulates the spatiotemporal dynamic of polarity molecules, leading to the establishment and maintenance of polar domains at the plasma membrane. Despite considerable progress in identifying key polarity regulators in plants, the molecular and cellular mechanisms underlying cell polarity formation have yet to be fully elucidated. Recent work suggests a critical role for membrane protein/lipid nanodomains in polarized morphogenesis in plants. One outstanding question is how the spatiotemporal dynamics of signaling nanodomains are controlled to achieve robust cell polarization. In this review, we first summarize the current state of knowledge on potential regulatory mechanisms of nanodomain dynamics, with a special focus on Rho-like GTPases from plants. We then discuss the pavement cell system as an example of how cells may integrate multiple signals and nanodomain-involved feedback mechanisms to achieve robust polarity. A mechanistic understanding of nanodomains' roles in plant cell polarity is still in the early stages and will remain an exciting area for future investigations.

Non-asymptotic transients away from steady states determine cellular responsiveness to dynamic spatial-temporal signals

Majority of the theory on cell polarization and the understanding of cellular sensing and responsiveness to localized chemical cues has been based on the idea that non-polarized and polarized cell states can be represented by stable asymptotic switching between them. The existing model classes that describe the dynamics of signaling networks underlying polarization are formulated within the framework of autonomous systems. However these models do not simultaneously capture both, robust maintenance of polarized state longer than the signal duration, and retained responsiveness to signals with complex spatial-temporal distribution. Based on recent experimental evidence for criticality organization of biochemical networks, we challenge the current concepts and demonstrate that non-asymptotic signaling dynamics arising at criticality uniquely ensures optimal responsiveness to changing chemoattractant fields. We provide a framework to characterize non-asymptotic dynamics of system's state trajectories through a non-autonomous treatment of the system, further emphasizing the importance of (long) transient dynamics, as well as the necessity to change the mathematical formalism when describing biological systems that operate in changing environments.

Cdc42 mobility and membrane flows regulate fission yeast cell shape and survival

Local Cdc42 GTPase activation promotes polarized exocytosis, resulting in membrane flows that deplete low-mobility membrane-associated proteins from the growth region. To investigate the self-organizing properties of the Cdc42 secretion-polarization system under membrane flow, we developed a reaction-diffusion particle model. The model includes positive feedback activation of Cdc42, hydrolysis by GTPase-activating proteins (GAPs), and flow-induced displacement by exo/endocytosis. Simulations show how polarization relies on flow-induced depletion of low mobility GAPs. To probe the role of Cdc42 mobility in the fission yeast Schizosaccharomyces pombe, we changed its membrane binding properties by replacing its prenylation site with 1, 2 or 3 repeats of the Rit1 C terminal membrane binding domain (ritC), yielding alleles with progressively lower unbinding and diffusion rates. Concordant modelling predictions and experimental observations show that lower Cdc42 mobility results in lower Cdc42 activation level and wider patches. Indeed, while Cdc42-1ritC cells are viable and polarized, Cdc42-2ritC polarize poorly and Cdc42-3ritC is inviable. The model further predicts that GAP depletion increases Cdc42 activity at the expense of loss of polarization. Experiments confirm this prediction, as deletion of Cdc42 GAPs restores viability to Cdc42-3ritC cells. Our combined experimental and modelling studies demonstrate how membrane flows are an integral part of Cdc42-driven pattern formation.

Spatiotemporal Coordination of Rac1 and Cdc42 at the Whole Cell Level during Cell Ruffling

Rho-GTPases are central regulators within a complex signaling network that controls cytoskeletal organization and cell movement. The network includes multiple GTPases, such as the most studied Rac1, Cdc42, and RhoA, along with their numerous effectors that provide mutual regulation through feedback loops. Here we investigate the temporal and spatial relationship between Rac1 and Cdc42 during membrane ruffling, using a simulation model that couples GTPase signaling with cell morphodynamics and captures the GTPase behavior observed with FRET-based biosensors. We show that membrane velocity is regulated by the kinetic rate of GTPase activation rather than the concentration of active GTPase. Our model captures both uniform and polarized ruffling. We also show that cell-type specific time delays between Rac1 and Cdc42 activation can be reproduced with a single signaling motif, in which the delay is controlled by feedback from Cdc42 to Rac1. The resolution of our simulation output matches those of time-lapsed recordings of cell dynamics and GTPase activity. Our data-driven modeling approach allows us to validate simulation results with quantitative precision using the same pipeline for the analysis of simulated and experimental data.

A tractable physical model for the yeast polarity predicts epistasis and fitness

Accurate phenotype prediction based on genetic information has numerous societal applications, such as crop design or cellular factories. Epistasis, when biological components interact, complicates modelling phenotypes from genotypes. Here we show an approach to mitigate this complication for polarity establishment in budding yeast, where mechanistic information is abundant. We coarse-grain molecular interactions into a so-called mesotype, which we combine with gene expression noise into a physical cell cycle model. First, we show with computer simulations that the mesotype allows validation of the most current biochemical polarity models by quantitatively matching doubling times. Second, the mesotype elucidates epistasis emergence as exemplified by evaluating the predicted mutational effect of key polarity protein Bem1p when combined with known interactors or under different growth conditions. This example also illustrates how unlikely evolutionary trajectories can become more accessible. The tractability of our biophysically justifiable approach inspires a road-map towards bottom-up modelling complementary to statistical inferences.

This article is part of the theme issue ‘Interdisciplinary approaches to predicting evolutionary biology’.

Spatiotemporal coordination of Rac1 and Cdc42 at the whole cell level during cell ruffling

Rho-GTPases are central regulators within a complex signaling network that controls the cytoskeletal organization and cell movement. This network includes multiple GTPases, such as the most studied Rac1, Cdc42, and RhoA, and their numerous effectors that provide mutual regulation and feedback loops. Here we investigate the temporal and spatial relationship between Rac1 and Cdc42 during membrane ruffling using a simulation model which couples GTPase signaling with cell morphodynamics to capture the GTPase behavior observed with FRET-based biosensors. We show that membrane velocity is regulated by the kinetic rate of GTPase activation rather than the concentration of active GTPase. Our model captures both uniform and polarized ruffling. We also show that cell-type specific time delays between Rac1 and Cdc42 activation can be reproduced with a single signaling motif, in which the delay is controlled by feedback from Cdc42 to Rac1. The resolution of our simulation output matches those of the time-lapsed recordings of cell dynamics and GTPase activity. This approach allows us to validate simulation results with quantitative precision using the same pipeline for the analysis of simulated and experimental data.

Preferential localization of a single spot in reaction-diffusion systems on non-spherical surfaces

The present work systematically examines the effect of breaking the rotational symmetry of a surface on the spot positioning in reaction-diffusion (RD) systems. In particular, we study analytically and numerically the steady-state positioning of a single spot in RD systems on a prolate and an oblate ellipsoid. We adapt perturbative techniques to perform a linear stability analysis of the RD system on both ellipsoids. Furthermore, the spot positionings in the steady states of non-linear RD equations are obtained numerically on both ellipsoids. Our analysis suggests that preferential spot positioning can be observed on non-spherical surfaces. The present work may provide useful insights into the role of cell geometry on various symmetry-breaking mechanisms in cellular processes.

Complexity and self-organization in the evolution of cell polarization

Cellular life exhibits order and complexity, which typically increase over the course of evolution. Cell polarization is a well-studied example of an ordering process that breaks the internal symmetry of a cell by establishing a preferential axis. Like many cellular processes, polarization is driven by self-organization, meaning that the macroscopic pattern emerges as a consequence of microscopic molecular interactions at the biophysical level. However, the role of self-organization in the evolution of complex protein networks remains obscure. In this Review, we provide an overview of the evolution of polarization as a self-organizing process, focusing on the model species Saccharomyces cerevisiae and its fungal relatives. Moreover, we use this model system to discuss how self-organization might relate to evolutionary change, offering a shift in perspective on evolution at the microscopic scale.

Oligomerization of peripheral membrane proteins provides tunable control of cell surface polarity

Asymmetric distributions of peripheral membrane proteins define cell polarity across all kingdoms of life. Non-linear positive feedback on membrane binding is essential to amplify and stabilize these asymmetries, but how specific molecular sources of non-linearity shape polarization dynamics remains poorly understood. Here we show that the ability to oligomerize, which is common to many peripheral membrane proteins, can play a profound role in shaping polarization dynamics in simple feedback circuits. We show that size-dependent binding avidity and mobility of membrane-bound oligomers endow polarity circuits with several key properties. Size-dependent membrane binding avidity confers a form of positive feedback on the accumulation of oligomer subunits. Although insufficient by itself, this sharply reduces the amount of additional feedback required for spontaneous emergence and stable maintenance of polarized states. Size-dependent oligomer mobility makes symmetry breaking and stable polarity more robust with respect to variation in subunit diffusivities and cell sizes, and slows the approach to a final stable spatial distribution, allowing cells to "remember" polarity boundaries imposed by transient external cues. Together, these findings reveal how oligomerization of peripheral membrane proteins can provide powerful and highly tunable sources of non-linear feedback in biochemical circuits that govern cell surface polarity. Given its prevalence and widespread involvement in cell polarity, we speculate that self-oligomerization may have provided an accessible path to evolving simple polarity circuits.

The Cdc42 GAP Rga6 promotes monopolar outgrowth of spores

The molecular mechanisms underlying the establishment of the monopolar growth of fission yeast spores have been less characterized. Here, we report that the Cdc42 GTPase-activating protein (GAP) Rga6 is required for promoting monopolar growth during spore germination. The absence of Rga6 increases the number of spores that grow in a bipolar fashion. Rga6 decorates the non-growing cortical region, binds phosphatidylinositol 4,5-bisphosphate, and colocalizes with the phosphatidylinositol 4,5-bisphosphate-binding protein Opy1. Overexpression of Opy1 diminishes the cortical localization of Rga6. The characteristic localization of Rga6 on the cell cortex depends on the C-terminal PBR region of Rga6. Moreover, engineered chimera composed of the Rga6 C-terminal PBR region fused to the GAP domain of Rga3 or Rga4 are sufficient to rescue the spore growth phenotype caused by the absence of Rga6. Hence, our work establishes a paradigm in which the lipid composition of the plasma membrane directs polarized cell growth by specifying the cortical localization of a GAP protein.

Cell-Size Space Regulates the Behavior of Confined Polymers: From Nano- and Micromaterials Science to Biology

Polymer micromaterials in a liquid or gel phase covered with a surfactant membrane are widely used materials in pharmaceuticals, cosmetics, and foods. In particular, cell-sized micromaterials of biopolymer solutions covered with a lipid membrane have been studied as artificial cells to understand cells from a physicochemical perspective. The characteristics and phase transitions of polymers confined to a microscopic space often differ from those in bulk systems. The effect that causes this difference is referred to as the cell-size space effect (CSE), but the specific physicochemical factors remain unclear. This study introduces the analysis of CSE on molecular diffusion, nanostructure transition, and phase separation and presents their main factors, i.e., short- and long-range interactions with the membrane surface and small volume (finite element nature). This serves as a guide for determining the dominant factors of CSE. Furthermore, we also introduce other factors of CSE such as spatial closure and the relationships among space size, the characteristic length of periodicity, the structure size, and many others produced by biomolecular assemblies through the analysis of protein reaction-diffusion systems and biochemical reactions.

Spatial models of pattern formation during phagocytosis

Phagocytosis, the biological process in which cells ingest large particles such as bacteria, is a key component of the innate immune response. Fcγ receptor (FcγR)-mediated phagocytosis is initiated when these receptors are activated after binding immunoglobulin G (IgG). Receptor activation initiates a signaling cascade that leads to the formation of the phagocytic cup and culminates with ingestion of the foreign particle. In the experimental system termed "frustrated phagocytosis", cells attempt to internalize micropatterned disks of IgG. Cells that engage in frustrated phagocytosis form "rosettes" of actin-enriched structures called podosomes around the IgG disk. The mechanism that generates the rosette pattern is unknown. We present data that supports the involvement of Cdc42, a member of the Rho family of GTPases, in pattern formation. Cdc42 acts downstream of receptor activation, upstream of actin polymerization, and is known to play a role in polarity establishment. Reaction-diffusion models for GTPase spatiotemporal dynamics exist. We demonstrate how the addition of negative feedback and minor changes to these models can generate the experimentally observed rosette pattern of podosomes. We show that this pattern formation can occur through two general mechanisms. In the first mechanism, an intermediate species forms a ring of high activity around the IgG disk, which then promotes rosette organization. The second mechanism does not require initial ring formation but relies on spatial gradients of intermediate chemical species that are selectively activated over the IgG patch. Finally, we analyze the models to suggest experiments to test their validity.

Design centering enables robustness screening of pattern formation models

MOTIVATION

Access to unprecedented amounts of quantitative biological data allows us to build and test biochemically accurate reaction-diffusion models of intracellular processes. However, any increase in model complexity increases the number of unknown parameters and, thus, the computational cost of model analysis. To efficiently characterize the behavior and robustness of models with many unknown parameters remains, therefore, a key challenge in systems biology.

RESULTS

We propose a novel computational framework for efficient high-dimensional parameter space characterization of reaction-diffusion models in systems biology. The method leverages the Lp-Adaptation algorithm, an adaptive-proposal statistical method for approximate design centering and robustness estimation. Our approach is based on an oracle function, which predicts for any given point in parameter space whether the model fulfills given specifications. We propose specific oracles to efficiently predict four characteristics of Turing-type reaction-diffusion models: bistability, instability, capability of spontaneous pattern formation and capability of pattern maintenance. We benchmark the method and demonstrate that it enables global exploration of a model's ability to undergo pattern-forming instabilities and to quantify robustness for model selection in polynomial time with dimensionality. We present an application of the framework to pattern formation on the endosomal membrane by the small GTPase Rab5 and its effectors, and we propose molecular mechanisms underlying this system.

AVAILABILITY AND IMPLEMENTATION

Our code is implemented in MATLAB and is available as open source under https://git.mpi-cbg.de/mosaic/software/black-box-optimization/rd-parameter-space-screening.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Cell Repolarization: A Bifurcation Study of Spatio-Temporal Perturbations of Polar Cells

The intrinsic polarity of migrating cells is regulated by spatial distributions of protein activity. Those proteins (Rho-family GTPases, such as Rac and Rho) redistribute in response to stimuli, determining the cell front and back. Reaction-diffusion equations with mass conservation and positive feedback have been used to explain initial polarization of a cell. However, the sensitivity of a polar cell to a reversal stimulus has not yet been fully understood. We carry out a PDE bifurcation analysis of two polarity models to investigate routes to repolarization: (1) a single-GTPase ("wave-pinning") model and (2) a mutually antagonistic Rac-Rho model. We find distinct routes to reversal in (1) vs. (2). We show numerical simulations of full PDE solutions for the RD equations, demonstrating agreement with predictions of the bifurcation results. Finally, we show that simulations of the polarity models in deforming 1D model cells are consistent with biological experiments.

Sensing the shape of a cell with reaction diffusion and energy minimization

Some dividing cells sense their shape by becoming polarized along their long axis. Cell polarity is controlled in part by polarity proteins, like Rho GTPases, cycling between active membrane-bound forms and inactive cytosolic forms, modeled as a "wave-pinning" reaction-diffusion process. Does shape sensing emerge from wave pinning? We show that wave pinning senses the cell's long axis. Simulating wave pinning on a curved surface, we find that high-activity domains migrate to peaks and troughs of the surface. For smooth surfaces, a simple rule of minimizing the domain perimeter while keeping its area fixed predicts the final position of the domain and its shape. However, when we introduce roughness to our surfaces, shape sensing can be disrupted, and high-activity domains can become localized to locations other than the global peaks and valleys of the surface. On rough surfaces, the domains of the wave-pinning model are more robust in finding the peaks and troughs than the minimization rule, although both can become trapped in steady states away from the peaks and valleys. We can control the robustness of shape sensing by altering the Rho GTPase diffusivity and the domain size. We also find that the shape-sensing properties of cell polarity models can explain how domains localize to curved regions of deformed cells. Our results help to understand the factors that allow cells to sense their shape-and the limits that membrane roughness can place on this process.

A versatile cortical pattern-forming circuit based on Rho, F-actin, Ect2, and RGA-3/4

Many cells can generate complementary traveling waves of actin filaments (F-actin) and cytoskeletal regulators. This phenomenon, termed cortical excitability, results from coupled positive and negative feedback loops of cytoskeletal regulators. The nature of these feedback loops, however, remains poorly understood. We assessed the role of the Rho GAP RGA-3/4 in the cortical excitability that accompanies cytokinesis in both frog and starfish. RGA-3/4 localizes to the cytokinetic apparatus, "chases" Rho waves in an F-actin-dependent manner, and when coexpressed with the Rho GEF Ect2, is sufficient to convert the normally quiescent, immature Xenopus oocyte cortex into a dramatically excited state. Experiments and modeling show that changing the ratio of RGA-3/4 to Ect2 produces cortical behaviors ranging from pulses to complex waves of Rho activity. We conclude that RGA-3/4, Ect2, Rho, and F-actin form the core of a versatile circuit that drives a diverse range of cortical behaviors, and we demonstrate that the immature oocyte is a powerful model for characterizing these dynamics.

Forced and spontaneous symmetry breaking in cell polarization

How does breaking the symmetry of an equation alter the symmetry of its solutions? Here, we systematically examine how reducing underlying symmetries from spherical to axisymmetric influences the dynamics of an archetypal model of cell polarization, a key process of biological spatial self-organization. Cell polarization is characterized by nonlinear and non-local dynamics, but we overcome the theory challenges these traits pose by introducing a broadly applicable numerical scheme allowing us to efficiently study continuum models in a wide range of geometries. Guided by numerical results, we discover a dynamical hierarchy of timescales that allows us to reduce relaxation to a purely geometric problem of area-preserving geodesic curvature flow. Through application of variational results, we analytically construct steady states on a number of biologically relevant shapes. In doing so, we reveal non-trivial solutions for symmetry breaking.

Mode selection mechanism in traveling and standing waves revealed by Min wave reconstituted in artificial cells

Reaction-diffusion coupling (RDc) generates spatiotemporal patterns, including two dynamic wave modes: traveling and standing waves. Although mode selection plays a substantial role in the spatiotemporal organization of living cell molecules, the mechanism for selecting each wave mode remains elusive. Here, we investigated a wave mode selection mechanism using Min waves reconstituted in artificial cells, emerged by the RDc of MinD and MinE. Our experiments and theoretical analysis revealed that the balance of membrane binding and dissociation from the membrane of MinD determines the mode selection of the Min wave. We successfully demonstrated that the transition of the wave modes can be regulated by controlling this balance and found hysteresis characteristics in the wave mode transition. These findings highlight a previously unidentified role of the balance between activators and inhibitors as a determinant of the mode selection of waves by RDc and depict an unexplored mechanism in intracellular spatiotemporal pattern formations.

Cells use molecular working memory to navigate inchanging chemoattractant fields

In order to migrate over large distances, cells within tissues and organisms rely on sensing local gradient cues which are irregular, conflicting, and changing over time and space. The mechanism how they generate persistent directional migration when signals are disrupted, while still remaining adaptive to signal's localization changes remain unknown. Here we find that single cells utilize a molecular mechanism akin to a working memory to satisfy these two opposing demands. We derive theoretically that this is characteristic for receptor networks maintained away from steady states. Time-resolved live-cell imaging of Epidermal growth factor receptor (EGFR) phosphorylation dynamics shows that cells transiently memorize position of encountered signals via slow-escaping remnant of the polarized signaling state, a dynamical 'ghost', driving memory-guided persistent directional migration. The metastability of this state further enables migrational adaptation when encountering new signals. We thus identify basic mechanism of real-time computations underlying cellular navigation in changing chemoattractant fields.

Chemotropism and Cell-Cell Fusion in Fungi

Fungi exhibit an enormous variety of morphologies, including yeast colonies, hyphal mycelia, and elaborate fruiting bodies. This diversity arises through a combination of polar growth, cell division, and cell fusion. Because fungal cells are nonmotile and surrounded by a protective cell wall that is essential for cell integrity, potential fusion partners must grow toward each other until they touch and then degrade the intervening cell walls without impacting cell integrity. Here, we review recent progress on understanding how fungi overcome these challenges. Extracellular chemoattractants, including small peptide pheromones, mediate communication between potential fusion partners, promoting the local activation of core cell polarity regulators to orient polar growth and cell wall degradation. However, in crowded environments, pheromone gradients can be complex and potentially confusing, raising the question of how cells can effectively find their partners. Recent findings suggest that the cell polarity circuit exhibits searching behavior that can respond to pheromone cues through a remarkably flexible and effective strategy called exploratory polarization.

A positive feedback loop involving the Spa2 SHD domain contributes to focal polarization

Focal polarization is necessary for finely arranged cell-cell interactions. The yeast mating projection, with its punctate polarisome, is a good model system for this process. We explored the critical role of the polarisome scaffold protein Spa2 during yeast mating with a hypothesis motivated by mathematical modeling and tested by in vivo experiments. Our simulations predicted that two positive feedback loops generate focal polarization, including a novel feedback pathway involving the N-terminal domain of Spa2. We characterized the latter using loss-of-function and gain-of-function mutants. The N-terminal region contains a Spa2 Homology Domain (SHD) which is conserved from yeast to humans, and when mutated largely reproduced the spa2Δ phenotype. Our work together with published data show that the SHD domain recruits Msb3/4 that stimulates Sec4-mediated transport of Bud6 to the polarisome. There, Bud6 activates Bni1-catalyzed actin cable formation, recruiting more Spa2 and completing the positive feedback loop. We demonstrate that disrupting this loop at any point results in morphological defects. Gain-of-function perturbations partially restored focal polarization in a spa2 loss-of-function mutant without restoring localization of upstream components, thus supporting the pathway order. Thus, we have collected data consistent with a novel positive feedback loop that contributes to focal polarization during pheromone-induced polarization in yeast.

Spatiotemporal development of coexisting wave domains of Rho activity in the cell cortex

The Rho family GTPases are molecular switches that regulate cytoskeletal dynamics and cell movement through a complex spatiotemporal organization of their activity. In Patiria miniata (starfish) oocytes under in vitro experimental conditions (with overexpressed Ect2, induced expression of Δ90 cyclin B, and roscovitine treatment), such activity generates multiple co-existing regions of coherent propagation of actin waves. Here we use computational modeling to investigate the development and properties of such wave domains. The model reveals that the formation of wave domains requires a balance between the activation and inhibition in the Rho signaling motif. Intriguingly, the development of the wave domains is preceded by a stage of low-activity quasi-static patterns, which may not be readily observed in experiments. Spatiotemporal patterns of this stage and the different paths of their destabilization define the behavior of the system in the later high-activity (observable) stage. Accounting for a strong intrinsic noise allowed us to achieve good quantitative agreement between simulated dynamics in different parameter regimes of the model and different wave dynamics in Patiria miniata and wild type Xenopus laevis (frog) data. For quantitative comparison of simulated and experimental results, we developed an automated method of wave domain detection, which revealed a sharp reversal in the process of pattern formation in starfish oocytes. Overall, our findings provide an insight into spatiotemporal regulation of complex and diverse but still computationally reproducible cell-level actin dynamics.

Endocytosis in the context-dependent regulation of individual and collective cell properties

Endocytosis allows cells to transport particles and molecules across the plasma membrane. In addition, it is involved in the termination of signalling through receptor downmodulation and degradation. This traditional outlook has been substantially modified in recent years by discoveries that endocytosis and subsequent trafficking routes have a profound impact on the positive regulation and propagation of signals, being key for the spatiotemporal regulation of signal transmission in cells. Accordingly, endocytosis and membrane trafficking regulate virtually every aspect of cell physiology and are frequently subverted in pathological conditions. Two key aspects of endocytic control over signalling are coming into focus: context-dependency and long-range effects. First, endocytic-regulated outputs are not stereotyped but heavily dependent on the cell-specific regulation of endocytic networks. Second, endocytic regulation has an impact not only on individual cells but also on the behaviour of cellular collectives. Herein, we will discuss recent advancements in these areas, highlighting how endocytic trafficking impacts complex cell properties, including cell polarity and collective cell migration, and the relevance of these mechanisms to disease, in particular cancer.

Cdc42 GTPase Activating Proteins (GAPs) Regulate Generational Inheritance of Cell Polarity and Cell Shape in Fission Yeast

The highly conserved small GTPase Cdc42 regulates polarized cell growth and morphogenesis from yeast to humans. We previously reported that Cdc42 activation exhibits oscillatory dynamics at cell tips of Schizosaccharomyces pombe cells. Mathematical modeling suggests that this dynamic behavior enables a variety of symmetric and asymmetric Cdc42 activation distributions to coexist in cell populations. For individual wild-type cells, however, Cdc42 distribution is initially asymmetrical and becomes more symmetrical as cell volume increases, enabling bipolar growth activation. To explore whether different patterns of Cdc42 activation are possible in vivo, we examined S. pombe rga4∆ mutant cells, lacking the Cdc42 GTPase-activating protein (GAP) Rga4. We found that monopolar rga4∆ mother cells divide asymmetrically leading to the emergence of both symmetric and asymmetric Cdc42 distributions in rga4∆ daughter cells. Motivated by different hypotheses that can mathematically reproduce the unequal fate of daughter cells, we used genetic screening to identify mutants that alter the rga4∆ phenotype. We found that the unequal distribution of active Cdc42 GTPase is consistent with an unequal inheritance of another Cdc42 GAP, Rga6, in the two daughter cells. Our findings highlight the crucial role of Cdc42 GAP localization in maintaining consistent Cdc42 activation and growth patterns across generations.

A novel stochastic simulation approach enables exploration of mechanisms for regulating polarity site movement

Cells polarize their movement or growth toward external directional cues in many different contexts. For example, budding yeast cells grow toward potential mating partners in response to pheromone gradients. Directed growth is controlled by polarity factors that assemble into clusters at the cell membrane. The clusters assemble, disassemble, and move between different regions of the membrane before eventually forming a stable polarity site directed toward the pheromone source. Pathways that regulate clustering have been identified but the molecular mechanisms that regulate cluster mobility are not well understood. To gain insight into the contribution of chemical noise to cluster behavior we simulated clustering using the reaction-diffusion master equation (RDME) framework to account for molecular-level fluctuations. RDME simulations are a computationally efficient approximation, but their results can diverge from the underlying microscopic dynamics. We implemented novel concentration-dependent rate constants that improved the accuracy of RDME-based simulations, allowing us to efficiently investigate how cluster dynamics might be regulated. Molecular noise was effective in relocating clusters when the clusters contained low numbers of limiting polarity factors, and when Cdc42, the central polarity regulator, exhibited short dwell times at the polarity site. Cluster stabilization occurred when abundances or binding rates were altered to either lengthen dwell times or increase the number of polarity molecules in the cluster. We validated key results using full 3D particle-based simulations. Understanding the mechanisms cells use to regulate the dynamics of polarity clusters should provide insights into how cells dynamically track external directional cues.

Cells localize polarity molecules in a small region of the plasma membrane forming a polarity cluster that directs functions such as migration, reproduction, and growth. Guided by external signals, these clusters move across the membrane allowing cells to reorient growth or motion. The polarity molecules continuously and randomly shuttle between the cluster and the cell cytosol and, as a result, the number and distribution of molecules at the cluster constantly changes. Here we present an improved stochastic simulation algorithm to investigate how such molecular-scale fluctuations induce cluster movement across the cell membrane. Unexpectedly, cluster mobility does not correlate with variations in total molecule abundance within the cluster, but rather with changes in the spatial distribution of molecules that form the cluster. Cluster motion is faster when polarity molecules are scarce and when they shuttle rapidly between the cluster and the cytosol. Our results suggest that cells control cluster mobility by regulating the abundance of polarity molecules and biochemical reactions that affect the time molecules spend at the cluster. We provide insights into how cells harness random molecular behavior to perform functions important for survival, such as detecting the direction of external signals.

Chemotactic movement of a polarity site enables yeast cells to find their mates

How small eukaryotic cells can interpret dynamic, noisy, and spatially complex chemical gradients to orient growth or movement is poorly understood. We address this question using Saccharomyces cerevisiae, where cells orient polarity up pheromone gradients during mating. Initial orientation is often incorrect, but polarity sites then move around the cortex in a search for partners. We find that this movement is biased by local pheromone gradients across the polarity site: that is, movement of the polarity site is chemotactic. A bottom-up computational model recapitulates this biased movement. The model reveals how even though pheromone-bound receptors do not mimic the shape of external pheromone gradients, nonlinear and stochastic effects combine to generate effective gradient tracking. This mechanism for gradient tracking may be applicable to any cell that searches for a target in a complex chemical landscape.

Type V myosin focuses the polarisome and shapes the tip of yeast cells

The polarisome is a cortical proteinaceous microcompartment that organizes the growth of actin filaments and the fusion of secretory vesicles in yeasts and filamentous fungi. Polarisomes are compact, spotlike structures at the growing tips of their respective cells. The molecular forces that control the form and size of this microcompartment are not known. Here we identify a complex between the polarisome subunit Pea2 and the type V Myosin Myo2 that anchors Myo2 at the cortex of yeast cells. We discovered a point mutation in the cargo-binding domain of Myo2 that impairs the interaction with Pea2 and consequently the formation and focused localization of the polarisome. Cells carrying this mutation grow round instead of elongated buds. Further experiments and biophysical modeling suggest that the interactions between polarisome-bound Myo2 motors and dynamic actin filaments spatially focus the polarisome and sustain its compact shape.

How cells determine the number of polarity sites

The diversity of cell morphologies arises, in part, through regulation of cell polarity by Rho-family GTPases. A poorly understood but fundamental question concerns the regulatory mechanisms by which different cells generate different numbers of polarity sites. Mass-conserved activator-substrate (MCAS) models that describe polarity circuits develop multiple initial polarity sites, but then those sites engage in competition, leaving a single winner. Theoretical analyses predicted that competition would slow dramatically as GTPase concentrations at different polarity sites increase toward a 'saturation point', allowing polarity sites to coexist. Here, we test this prediction using budding yeast cells, and confirm that increasing the amount of key polarity proteins results in multiple polarity sites and simultaneous budding. Further, we elucidate a novel design principle whereby cells can switch from competition to equalization among polarity sites. These findings provide insight into how cells with diverse morphologies may determine the number of polarity sites.

Function of membrane domains in rho-of-plant signaling

In a crowded environment, establishing interactions between different molecular partners can take a long time. Biological membranes have solved this issue, as they simultaneously are fluid and possess compartmentalized domains. This nanoscale organization of the membrane is often based on weak, local, and multivalent interactions between lipids and proteins. However, from local interactions at the nanoscale, different functional properties emerge at the higher scale, and these are critical to regulate and integrate cellular signaling. Rho of Plant (ROP) proteins are small guanosine triphosphate hydrolase enzymes (GTPases) involved in hormonal, biotic, and abiotic signaling, as well as fundamental cell biological properties such as polarity, vesicular trafficking, and cytoskeleton dynamics. Association with the membrane is essential for ROP function, as well as their precise targeting within micrometer-sized polar domains (i.e. microdomains) and nanometer-sized clusters (i.e. nanodomains). Here, we review our current knowledge about the formation and the maintenance of the ROP domains in membranes. Furthermore, we propose a model for ROP membrane targeting and discuss how the nanoscale organization of ROPs in membranes could determine signaling parameters like signal specificity, amplification, and integration.

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.

De novo design of a reversible phosphorylation-dependent switch for membrane targeting

Modules that switch protein-protein interactions on and off are essential to develop synthetic biology; for example, to construct orthogonal signaling pathways, to control artificial protein structures dynamically, and for protein localization in cells or protocells. In nature, the E. coli MinCDE system couples nucleotide-dependent switching of MinD dimerization to membrane targeting to trigger spatiotemporal pattern formation. Here we present a de novo peptide-based molecular switch that toggles reversibly between monomer and dimer in response to phosphorylation and dephosphorylation. In combination with other modules, we construct fusion proteins that couple switching to lipid-membrane targeting by: (i) tethering a ‘cargo’ molecule reversibly to a permanent membrane ‘anchor’; and (ii) creating a ‘membrane-avidity switch’ that mimics the MinD system but operates by reversible phosphorylation. These minimal, de novo molecular switches have potential applications for introducing dynamic processes into designed and engineered proteins to augment functions in living cells and add functionality to protocells.

The ability to dynamically control protein-protein interactions and localization of proteins is critical in synthetic biological systems. Here the authors develop a peptide-based molecular switch that regulates dimer formation and lipid membrane targeting via reversible phosphorylation.

Multi-spike solutions of a hybrid reaction–transport model

Simulations of classical pattern-forming reaction–diffusion systems indicate that they often operate in the strongly nonlinear regime, with the final steady state consisting of a spatially repeating pattern of localized spikes. In activator–inhibitor systems such as the two-component Gierer–Meinhardt (GM) model, one can consider the singular limit D a  ≪  D h , where D a and D h are the diffusivities of the activator and inhibitor, respectively. Asymptotic analysis can then be used to analyse the existence and linear stability of multi-spike solutions. In this paper, we analyse multi-spike solutions in a hybrid reaction–transport model, consisting of a slowly diffusing activator and an actively transported inhibitor that switches at a rate α between right-moving and left-moving velocity states. Such a model was recently introduced to account for the formation and homeostatic regulation of synaptic puncta during larval development in Caenorhabditis elegans . We exploit the fact that the hybrid model can be mapped onto the classical GM model in the fast switching limit α  → ∞, which establishes the existence of multi-spike solutions. Linearization about the multi-spike solution yields a non-local eigenvalue problem that is used to investigate stability of the multi-spike solution by combining analytical results for α  → ∞ with a graphical construction for finite α .

Cell polarisation in a bulk-surface model can be driven by both classic and non-classic Turing instability

The GTPase Cdc42 is the master regulator of eukaryotic cell polarisation. During this process, the active form of Cdc42 is accumulated at a particular site on the cell membrane called the pole. It is believed that the accumulation of the active Cdc42 resulting in a pole is driven by a combination of activation–inactivation reactions and diffusion. It has been proposed using mathematical modelling that this is the result of diffusion-driven instability, originally proposed by Alan Turing. In this study, we developed, analysed and validated a 3D bulk-surface model of the dynamics of Cdc42. We show that the model can undergo both classic and non-classic Turing instability by deriving necessary conditions for which this occurs and conclude that the non-classic case can be viewed as a limit case of the classic case of diffusion-driven instability. Using three-dimensional Spatio-temporal simulation we predicted pole size and time to polarisation, suggesting that cell polarisation is mainly driven by the reaction strength parameter and that the size of the pole is determined by the relative diffusion.

A mathematical model for bleb regulation in zebrafish primordial germ cells

Blebs are cell protrusions generated by local membrane-cortex detachments followed by expansion of the plasma membrane. Blebs are formed by some migrating cells, e.g. primordial germ cells of the zebrafish. While blebs occur randomly at each part of the membrane in unpolarized cells, a polarization process guarantees the occurrence of blebs at a preferential site and thereby facilitates migration toward a specified direction. Little is known about the factors involved in the controlled and directed bleb generation, yet recent studies revealed the influence of an intracellular flow and the stabilizing role of the membrane-cortex linker molecule Ezrin. Based on this information, we develop and analyse a coupled bulk-surface model describing a potential cellular mechanism by which a bleb could be induced at a controlled site. The model rests upon intracellular Darcy flow and a diffusion-advection-reaction system, describing the temporal evolution from a homogeneous to a strongly anisotropic Ezrin distribution. We prove the well-posedness of the mathematical model and show that simulations qualitatively correspond to experimental observations, suggesting that indeed the interaction of an intracellular flow with membrane proteins can be the cause of the Ezrin redistribution accompanying bleb formation.

Mechanisms that ensure monogamous mating in Saccharomyces cerevisiae

Haploid cells of the budding yeast Saccharomyces cerevisiae communicate using secreted pheromones and mate to form diploid zygotes. Mating is monogamous, resulting in the fusion of precisely one cell of each mating type. Monogamous mating in crowded conditions, where cells have access to more than one potential partner, raises the question of how multiple-mating outcomes are prevented. Here we identify mutants capable of mating with multiple partners, revealing the mechanisms that ensure monogamous mating. Before fusion, cells develop polarity foci oriented toward potential partners. Competition between these polarity foci within each cell leads to disassembly of all but one focus, thus favoring a single fusion event. Fusion promotes the formation of heterodimeric complexes between subunits that are uniquely expressed in each mating type. One complex shuts off haploid-specific gene expression, and the other shuts off the ability to respond to pheromone. Zygotes able to form either complex remain monogamous, but zygotes lacking both can re-mate.

Developmental clock and mechanism of de novo polarization of the mouse embryo

Embryo polarization is critical for mouse development; however, neither the regulatory clock nor the molecular trigger that it activates is known. Here, we show that the embryo polarization clock reflects the onset of zygotic genome activation, and we identify three factors required to trigger polarization. Advancing the timing of transcription factor AP-2 gamma (Tfap2c) and TEA domain transcription factor 4 (Tead4) expression in the presence of activated Ras homolog family member A (RhoA) induces precocious polarization as well as subsequent cell fate specification and morphogenesis. Tfap2c and Tead4 induce expression of actin regulators that control the recruitment of apical proteins on the membrane, whereas RhoA regulates their lateral mobility, allowing the emergence of the apical domain. Thus, Tfap2c, Tead4, and RhoA are regulators for the onset of polarization and cell fate segregation in the mouse.

Coordinating cell polarization and morphogenesis through mechanical feedback

Many cellular processes require cell polarization to be maintained as the cell changes shape, grows or moves. Without feedback mechanisms relaying information about cell shape to the polarity molecular machinery, the coordination between cell polarization and morphogenesis, movement or growth would not be possible. Here we theoretically and computationally study the role of a genetically-encoded mechanical feedback (in the Cell Wall Integrity pathway) as a potential coordination mechanism between cell morphogenesis and polarity during budding yeast mating projection growth. We developed a coarse-grained continuum description of the coupled dynamics of cell polarization and morphogenesis as well as 3D stochastic simulations of the molecular polarization machinery in the evolving cell shape. Both theoretical approaches show that in the absence of mechanical feedback (or in the presence of weak feedback), cell polarity cannot be maintained at the projection tip during growth, with the polarization cap wandering off the projection tip, arresting morphogenesis. In contrast, for mechanical feedback strengths above a threshold, cells can robustly maintain cell polarization at the tip and simultaneously sustain mating projection growth. These results indicate that the mechanical feedback encoded in the Cell Wall Integrity pathway can provide important positional information to the molecular machinery in the cell, thereby enabling the coordination of cell polarization and morphogenesis.

The Path towards Predicting Evolution as Illustrated in Yeast Cell Polarity

A bottom-up route towards predicting evolution relies on a deep understanding of the complex network that proteins form inside cells. In a rapidly expanding panorama of experimental possibilities, the most difficult question is how to conceptually approach the disentangling of such complex networks. These can exhibit varying degrees of hierarchy and modularity, which obfuscate certain protein functions that may prove pivotal for adaptation. Using the well-established polarity network in budding yeast as a case study, we first organize current literature to highlight protein entrenchments inside polarity. Following three examples, we see how alternating between experimental novelties and subsequent emerging design strategies can construct a layered understanding, potent enough to reveal evolutionary targets. We show that if you want to understand a cell’s evolutionary capacity, such as possible future evolutionary paths, seemingly unimportant proteins need to be mapped and studied. Finally, we generalize this research structure to be applicable to other systems of interest.

Compete or Coexist? Why the Same Mechanisms of Symmetry Breaking Can Yield Distinct Outcomes

Cellular morphogenesis is governed by the prepattern based on the symmetry-breaking emergence of dense protein clusters. Thus, a cluster of active GTPase Cdc42 marks the site of nascent bud in the baker's yeast. An important biological question is which mechanisms control the number of pattern maxima (spots) and, thus, the number of nascent cellular structures. Distinct flavors of theoretical models seem to suggest different predictions. While the classical Turing scenario leads to an array of stably coexisting multiple structures, mass-conserved models predict formation of a single spot that emerges via the greedy competition between the pattern maxima for the common molecular resources. Both the outcome and the kinetics of this competition are of significant biological importance but remained poorly explored. Recent theoretical analyses largely addressed these questions, but their results have not yet been fully appreciated by the broad biological community. Keeping mathematical apparatus and jargon to the minimum, we review the main conclusions of these analyses with their biological implications in mind. Focusing on the specific example of pattern formation by small GTPases, we speculate on the features of the patterning mechanisms that bypass competition and favor formation of multiple coexisting structures and contrast them with those of the mechanisms that harness competition to form unique cellular structures.

Septin architecture and function in budding yeast

The septins constitute a conserved family of guanosine phosphate-binding and filament-forming proteins widespread across eukaryotic species. Septins appear to have two principal functions. One is to form a cortical diffusion barrier, like the septin collar at the bud neck of Saccharomyces cerevisiae, which prevents movement of membrane-associated proteins between the mother and daughter cells. The second is to serve as a polymeric scaffold for recruiting the proteins required for critical cellular processes to particular subcellular areas. In the last decade, structural information about the different levels of septin organization has appeared, but crucial structural determinants and factors responsible for septin assembly remain largely unknown. This review highlights recent findings on the architecture and function of septins and their remodeling with an emphasis on mitotically dividing budding yeasts.