GDI-mediated cell polarization in yeast provides precise spatial and temporal control of Cdc42 signaling

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

Polarized exocytosis induced by local Cdc42 GTPase activity results in membrane flows that deplete low-mobility membrane-associated proteins. A reaction-diffusion particle model comprising Cdc42 positive feedback activation, hydrolysis by GTPase-activating proteins (GAPs), and flow-induced displacement by exo/endocytosis shows that flow-induced depletion of low mobility GAPs promotes polarization. We modified Cdc42 mobility in Schizosaccharomyces pombe by replacing its prenylation site with 1, 2 or 3 repeats of the Rit C-terminal membrane-binding domain (ritC), yielding alleles with progressively lower mobility and increased flow-coupling. While Cdc42-1ritC cells are viable and polarized, Cdc42-2ritC polarize poorly and Cdc42-3ritC are inviable, in agreement with model's predictions. Deletion of Cdc42 GAPs restores viability to Cdc42-3ritC cells, verifying the model's prediction that GAP deletion increases Cdc42 activity at the expense of polarization. Our work demonstrates how membrane flows are an integral part of Cdc42-driven pattern formation and require Cdc42-GTP to turn over faster than the surface on which it forms.

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.

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.

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.

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.

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.

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’.

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.

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.

Orientation of Cell Polarity by Chemical Gradients

Accurate decoding of spatial chemical landscapes is critical for many cell functions. Eukaryotic cells decode local chemical gradients to orient growth or movement in productive directions. Recent work on yeast model systems, whose gradient sensing pathways display much less complexity than those in animal cells, has suggested new paradigms for how these very small cells successfully exploit information in noisy and dynamic pheromone gradients to identify their mates. Pheromone receptors regulate a polarity circuit centered on the conserved Rho-family GTPase, Cdc42. The polarity circuit contains both positive and negative feedback pathways, allowing spontaneous symmetry breaking and also polarity site disassembly and relocation. Cdc42 orients the actin cytoskeleton, leading to focused vesicle traffic that promotes movement of the polarity site and also reshapes the cortical distribution of receptors at the cell surface. In this article, we review the advances from work on yeasts and compare them with the excitable signaling pathways that have been revealed in chemotactic animal cells. Expected final online publication date for the Annual Review of Biophysics, Volume 51 is May 2022. Please see http://www.annualreviews.org/page/journal/pubdates for revised estimates.

Formation of self-organizing functionally distinct Rho of plants domains involves a reduced mobile population

Rho family proteins are central to the regulation of cell polarity in eukaryotes. Rho of Plants-Guanyl nucleotide Exchange Factor (ROPGEF) can form self-organizing polar domains following co-expression with an Rho of Plants (ROP) and an ROP GTPase-Activating Protein (ROPGAP). Localization of ROPs in these domains has not been demonstrated, and the mechanisms underlying domain formation and function are not well understood. Here we show that six different ROPs form self-organizing domains when co-expressed with ROPGEF3 and GAP1 in Nicotiana benthamiana or Arabidopsis (Arabidopsis thaliana). Domain formation was associated with ROP-ROPGEF3 association, reduced ROP mobility, as revealed by time-lapse imaging and Fluorescence Recovery After Photobleaching beam size analysis, and was independent of Rho GTP Dissociation Inhibitor mediated recycling. The domain formation depended on the ROPs' activation/inactivation cycles and interaction with anionic lipids via a C-terminal polybasic domain. Coexpression with the microtubule-associated protein ROP effector INTERACTOR OF CONSTITUTIVELY ACTIVE ROP 1 (ICR1) revealed differential function of the ROP domains in the ability to recruit ICR1. Taken together, the results reveal mechanisms underlying self-organizing ROP domain formation and function.

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.

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.

Turing Patterning in Stratified Domains

Reaction-diffusion processes across layered media arise in several scientific domains such as pattern-forming E. coli on agar substrates, epidermal-mesenchymal coupling in development, and symmetry-breaking in cell polarization. We develop a modeling framework for bilayer reaction-diffusion systems and relate it to a range of existing models. We derive conditions for diffusion-driven instability of a spatially homogeneous equilibrium analogous to the classical conditions for a Turing instability in the simplest nontrivial setting where one domain has a standard reaction-diffusion system, and the other permits only diffusion. Due to the transverse coupling between these two regions, standard techniques for computing eigenfunctions of the Laplacian cannot be applied, and so we propose an alternative method to compute the dispersion relation directly. We compare instability conditions with full numerical simulations to demonstrate impacts of the geometry and coupling parameters on patterning, and explore various experimentally relevant asymptotic regimes. In the regime where the first domain is suitably thin, we recover a simple modulation of the standard Turing conditions, and find that often the broad impact of the diffusion-only domain is to reduce the ability of the system to form patterns. We also demonstrate complex impacts of this coupling on pattern formation. For instance, we exhibit non-monotonicity of pattern-forming instabilities with respect to geometric and coupling parameters, and highlight an instability from a nontrivial interaction between kinetics in one domain and diffusion in the other. These results are valuable for informing design choices in applications such as synthetic engineering of Turing patterns, but also for understanding the role of stratified media in modulating pattern-forming processes in developmental biology and beyond.

Targeting and signaling of Rho of plants guanosine triphosphatases require synergistic interaction between guanine nucleotide inhibitor and vesicular trafficking

Most eukaryotic cells are polarized. Common toolbox regulating cell polarization includes Rho guanosine triphosphatases (GTPases), in which spatiotemporal activation is regulated by a plethora of regulators. Rho of plants (ROPs) are the only Rho GTPases in plants. Although vesicular trafficking was hinted in the regulation of ROPs, it was unclear where vesicle-carried ROP starts, whether it is dynamically regulated, and which components participate in vesicle-mediated ROP targeting. In addition, although vesicle trafficking and guanine nucleotide inhibitor (GDI) pathways in Rho signaling have been extensively studied in yeast, it is unknown whether the two pathways interplay. Unclear are also cellular and developmental consequences of their interaction in multicellular organisms. Here, we show that the dynamic targeting of ROP through vesicles requires coat protein complex II and ADP-ribosylation factor 1-mediated post-Golgi trafficking. Trafficking of vesicle-carried ROPs between the plasma membrane and the trans-Golgi network is mediated through adaptor protein 1 and sterol-mediated endocytosis. Finally, we show that GDI and vesicle trafficking synergistically regulate cell polarization and ROP targeting, suggesting that the establishment and maintenance of cell polarity is regulated by an evolutionarily conserved mechanism.

Fission Yeast Polarization: Modeling Cdc42 Oscillations, Symmetry Breaking, and Zones of Activation and Inhibition

Cells polarize for growth, motion, or mating through regulation of membrane-bound small GTPases between active GTP-bound and inactive GDP-bound forms. Activators (GEFs, GTP exchange factors) and inhibitors (GAPs, GTPase activating proteins) provide positive and negative feedbacks. We show that a reaction-diffusion model on a curved surface accounts for key features of polarization of model organism fission yeast. The model implements Cdc42 membrane diffusion using measured values for diffusion coefficients and dissociation rates and assumes a limiting GEF pool (proteins Gef1 and Scd1), as in prior models for budding yeast. The model includes two types of GAPs, one representing tip-localized GAPs, such as Rga3; and one representing side-localized GAPs, such as Rga4 and Rga6, that we assume switch between fast and slow diffusing states. After adjustment of unknown rate constants, the model reproduces active Cdc42 zones at cell tips and the pattern of GEF and GAP localization at cell tips and sides. The model reproduces observed tip-to-tip oscillations with periods of the order of several minutes, as well as asymmetric to symmetric oscillations transitions (corresponding to NETO "new end take off"), assuming the limiting GEF amount increases with cell size.

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.

Cell size sets the diameter of the budding yeast contractile ring

The formation and maintenance of subcellular structures and organelles with a well-defined size is a key requirement for cell function, yet our understanding of the underlying size control mechanisms is limited. While budding yeast cell polarization and subsequent assembly of a septin ring at the site of bud formation has been successfully used as a model for biological self-assembly processes, the mechanisms that set the size of the septin ring at the bud neck are unknown. Here, we use live-cell imaging and genetic manipulation of cell volume to show that the septin ring diameter increases with cell volume. This cell-volume-dependence largely accounts for modulations of ring size due to changes in ploidy and genetic manipulation of cell polarization. Our findings suggest that the ring diameter is set through the dynamic interplay of septin recruitment and Cdc42 polarization, establishing it as a model for size homeostasis of self-assembling organelles.

Budding yeast cell polarization is known to self-assemble, but it is still not clear what controls the size of the resulting septin ring. Here the authors show that the septin ring diameter is set by cell volume, ensuring that larger cells have larger rings.

Quasi-periodic migration of single cells on short microlanes

Cell migration on microlanes represents a suitable and simple platform for the exploration of the molecular mechanisms underlying cell cytoskeleton dynamics. Here, we report on the quasi-periodic movement of cells confined in stripe-shaped microlanes. We observe persistent polarized cell shapes and directed pole-to-pole motion within the microlanes. Cells depolarize at one end of a given microlane, followed by delayed repolarization towards the opposite end. We analyze cell motility via the spatial velocity distribution, the velocity frequency spectrum and the reversal time as a measure for depolarization and spontaneous repolarization of cells at the microlane ends. The frequent encounters of a boundary in the stripe geometry provides a robust framework for quantitative investigations of the cytoskeleton protrusion and repolarization dynamics. In a first advance to rigorously test physical models of cell migration, we find that the statistics of the cell migration is recapitulated by a Cellular Potts model with a minimal description of cytoskeleton dynamics. Using LifeAct-GFP transfected cells and microlanes with differently shaped ends, we show that the local deformation of the leading cell edge in response to the tip geometry can locally either amplify or quench actin polymerization, while leaving the average reversal times unaffected.

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.

Intracellular organization in cell polarity - placing organelles into the polarity loop

Many studies have investigated the processes that support polarity establishment and maintenance in cells. On the one hand, polarity complexes at the cell cortex and their downstream signaling pathways have been assigned as major regulators of polarity. On the other hand, intracellular organelles and their polarized trafficking routes have emerged as important components of polarity. In this Review, we argue that rather than trying to identify the prime 'culprit', now it is time to consider all these players as a collective. We highlight that understanding the intimate coordination between the polarized cell cortex and the intracellular compass that is defined by organelle positioning is essential to capture the concept of polarity. After briefly reviewing how polarity emerges from a dynamic maintenance of cellular asymmetries, we highlight how intracellular organelles and their associated trafficking routes provide diverse feedback for dynamic cell polarity maintenance. We argue that the asymmetric organelle compass is an indispensable element of the polarity network.

The role of intracellular interactions in the collective polarization of tissues and its interplay with cellular geometry

Planar cell polarity (PCP), the long-range in-plane polarization of epithelial tissues, provides directional information that guides a multitude of developmental processes at cellular and tissue levels. While it is manifest that cells utilize both intracellular and intercellular interactions, the coupling between the two modules, essential to the coordination of collective polarization, remains an active area of investigation. We propose a generalized reaction-diffusion model to study the role of intracellular interactions in the emergence of long-range polarization, and show that the nonlocality of cytoplasmic interactions, i.e. coupling of membrane proteins localized on different cell-cell junctions, is of vital importance to the faithful detection of weak directional signals, and becomes increasingly more crucial to the stability of polarization against the deleterious effects of large geometric irregularities. We demonstrate that nonlocal interactions are necessary for geometric information to become accessible to the PCP components. The prediction of the model regarding polarization in elongated tissues, is shown to be in agreement with experimental observations, where the polarity emerges perpendicular to the axis of elongation. Core PCP is adopted as a model pathway, in term of which we interpret the model parameters. To this end, we introduce three distinct classes of mutations, (I) in membrane proteins, (II) in cytoplasmic proteins, and (III) local enhancement of geometric disorder. Comparing the in silico and in vivo phenotypes, we show that our model successfully recapitulates the salient phenotypic features of these mutations. Exploring the parameter space helps us shed light on the role of cytoplasmic proteins in cell-cell communications, and make falsifiable predictions regarding the cooperation of cytoplasmic and membrane proteins in the establishment of long-range polarization.

Planar cell polarity (PCP) is an indispensable and conserved pathway in morphogenesis. In spite of the advances in understanding the different modules of PCP, a comprehensive picture of the intracellular protein-protein interactions necessary for the emergence of long-range tissue polarity is still lacking. In order to address this question, we devised a generalized reaction-diffusion model, through which we investigated the role of cytoplasmic interactions in PCP pathways. The length scale of intracellular interactions is demonstrated to be crucial to the stability of the cytoplasmic segregation of membrane proteins in disordered tissues, as well as the capacity of polarization field for detecting the gradient and geometrical cues. Finally, three classes of mutants are investigated within the context of our model. Comparison with the in vivo observations allows us to infer the major contributions of cytoplasmic proteins to the emergence of tissue polarity, and make testable predictions regarding the cooperation of cytoplasmic and membrane proteins in the coordination of collective polarization.

Protein Recruitment through Indirect Mechanochemical Interactions

Some of the key proteins essential for important cellular processes are capable of recruiting other proteins from the cytosol to phospholipid membranes. The physical basis for this cooperativity of binding is, surprisingly, still unclear. Here, we suggest a general feedback mechanism that explains cooperativity through mechanochemical coupling mediated by the mechanical properties of phospholipid membranes. Our theory predicts that protein recruitment, and therefore also protein pattern formation, involves membrane deformation and is strongly affected by membrane composition.

Regulation of Cdc42 and its effectors in epithelial morphogenesis

Cdc42 - a member of the small Rho GTPase family - regulates cell polarity across organisms from yeast to humans. It is an essential regulator of polarized morphogenesis in epithelial cells, through coordination of apical membrane morphogenesis, lumen formation and junction maturation. In parallel, work in yeast and Caenorhabditis elegans has provided important clues as to how this molecular switch can generate and regulate polarity through localized activation or inhibition, and cytoskeleton regulation. Recent studies have revealed how important and complex these regulations can be during epithelial morphogenesis. This complexity is mirrored by the fact that Cdc42 can exert its function through many effector proteins. In epithelial cells, these include atypical PKC (aPKC, also known as PKC-3), the P21-activated kinase (PAK) family, myotonic dystrophy-related Cdc42 binding kinase beta (MRCKβ, also known as CDC42BPB) and neural Wiskott-Aldrich syndrome protein (N-WASp, also known as WASL). Here, we review how the spatial regulation of Cdc42 promotes polarity and polarized morphogenesis of the plasma membrane, with a focus on the epithelial cell type.

Self-organization principles of intracellular pattern formation

Dynamic patterning of specific proteins is essential for the spatio-temporal regulation of many important intracellular processes in prokaryotes, eukaryotes and multicellular organisms. The emergence of patterns generated by interactions of diffusing proteins is a paradigmatic example for self-organization. In this article, we review quantitative models for intracellular Min protein patterns in Escherichia coli, Cdc42 polarization in Saccharomyces cerevisiae and the bipolar PAR protein patterns found in Caenorhabditis elegans. By analysing the molecular processes driving these systems we derive a theoretical perspective on general principles underlying self-organized pattern formation. We argue that intracellular pattern formation is not captured by concepts such as ‘activators’, ‘inhibitors’ or ‘substrate depletion’. Instead, intracellular pattern formation is based on the redistribution of proteins by cytosolic diffusion, and the cycling of proteins between distinct conformational states. Therefore, mass-conserving reaction–diffusion equations provide the most appropriate framework to study intracellular pattern formation. We conclude that directed transport, e.g. cytosolic diffusion along an actively maintained cytosolic gradient, is the key process underlying pattern formation. Thus the basic principle of self-organization is the establishment and maintenance of directed transport by intracellular protein dynamics.

This article is part of the theme issue ‘Self-organization in cell biology’.

Comparison of Deterministic and Stochastic Regime in a Model for Cdc42 Oscillations in Fission Yeast

Oscillations occur in a wide variety of essential cellular processes, such as cell cycle progression, circadian clocks and calcium signaling in response to stimuli. It remains unclear how intrinsic stochasticity can influence these oscillatory systems. Here, we focus on oscillations of Cdc42 GTPase in fission yeast. We extend our previous deterministic model by Xu and Jilkine to construct a stochastic model, focusing on the fast diffusion case. We use SSA (Gillespie's algorithm) to numerically explore the low copy number regime in this model, and use analytical techniques to study the long-time behavior of the stochastic model and compare it to the equilibria of its deterministic counterpart. Numerical solutions suggest noisy limit cycles exist in the parameter regime in which the deterministic system converges to a stable limit cycle, and quasi-cycles exist in the parameter regime where the deterministic model has a damped oscillation. Near an infinite period bifurcation point, the deterministic model has a sustained oscillation, while stochastic trajectories start with an oscillatory mode and tend to approach deterministic steady states. In the low copy number regime, metastable transitions from oscillatory to steady behavior occur in the stochastic model. Our work contributes to the understanding of how stochastic chemical kinetics can affect a finite-dimensional dynamical system, and destabilize a deterministic steady state leading to oscillations.

A theoretical framework for planar polarity establishment through interpretation of graded cues by molecular bridges

Planar polarity is a widespread phenomenon found in many tissues, allowing cells to coordinate morphogenetic movements and function. A common feature of animal planar polarity systems is the formation of molecular bridges between cells, which become polarised along a tissue axis. We propose that these bridges provide a general mechanism by which cells interpret different forms of tissue gradients to coordinate directional information. We illustrate this using a generalised and consistent modelling framework, providing a conceptual basis for understanding how different mechanisms of gradient function can generate planar polarity. We make testable predictions of how different gradient mechanisms can influence polarity direction.

Summary: This Hypothesis uses a theoretical framework to explore how molecular bridges provide a general mechanism to interpret different forms of tissue gradients to establish planar polarity.

Minimal in vitro systems shed light on cell polarity

Cell polarity - the morphological and functional differentiation of cellular compartments in a directional manner - is required for processes such as orientation of cell division, directed cellular growth and motility. How the interplay of components within the complexity of a cell leads to cell polarity is still heavily debated. In this Review, we focus on one specific aspect of cell polarity: the non-uniform accumulation of proteins on the cell membrane. In cells, this is achieved through reaction-diffusion and/or cytoskeleton-based mechanisms. In reaction-diffusion systems, components are transformed into each other by chemical reactions and are moving through space by diffusion. In cytoskeleton-based processes, cellular components (i.e. proteins) are actively transported by microtubules (MTs) and actin filaments to specific locations in the cell. We examine how minimal systems - in vitro reconstitutions of a particular cellular function with a minimal number of components - are designed, how they contribute to our understanding of cell polarity (i.e. protein accumulation), and how they complement in vivo investigations. We start by discussing the Min protein system from Escherichia coli, which represents a reaction-diffusion system with a well-established minimal system. This is followed by a discussion of MT-based directed transport for cell polarity markers as an example of a cytoskeleton-based mechanism. To conclude, we discuss, as an example, the interplay of reaction-diffusion and cytoskeleton-based mechanisms during polarity establishment in budding yeast.

The shared role of the Rsr1 GTPase and Gic1/Gic2 in Cdc42 polarization

The Cdc42 GTPase plays a central role in polarity development in many species. In budding yeast, Cdc42 is essential for polarized growth at the proper site and also for spontaneous cell polarization in the absence of spatial cues. Cdc42 polarization is critical for multiple events in the G1 phase prior to bud emergence, including bud-site assembly, polarization of the actin cytoskeleton, and septin filament assembly to form a ring at the new bud site. Yet the mechanism by which Cdc42 polarizes is not fully understood. Here we report that biphasic Cdc42 polarization in the G1 phase is coupled to stepwise assembly of the septin ring for bud emergence. We show that the Rsr1 GTPase shares a partially redundant role with Gic1 and Gic2, two related Cdc42 effectors, in the first phase of Cdc42 polarization in haploid cells. We propose that the first phase of Cdc42 polarization is mediated by positive feedback loops that function in parallel-one involving Rsr1 via local activation of Cdc42 in response to spatial cues and another involving Gic1 or Gic2 via reduction of diffusion of active Cdc42.

Spherical Caps in Cell Polarization

Intracellular symmetry breaking plays a key role in wide range of biological processes, both in single cells and in multicellular organisms. An important class of symmetry-breaking mechanisms relies on the cytoplasm/membrane redistribution of proteins that can autocatalytically promote their own recruitment to the plasma membrane. We present an analytical construction and a comprehensive parametric analysis of stable localized patterns in a reaction-diffusion model of such a mechanism in a spherical cell. The constructed patterns take the form of high-concentration patches localized into spherical caps, similar to the patterns observed in the studies of symmetry breaking in single cells and early embryos.

The effect of cell geometry on polarization in budding yeast

The localization (or polarization) of proteins on the membrane during the mating of budding yeast (Saccharomyces cerevisiae) is an important model system for understanding simple pattern formation within cells. While there are many existing mathematical models of polarization, for both budding and mating, there are still many aspects of this process that are not well understood. In this paper we set out to elucidate the effect that the geometry of the cell can have on the dynamics of certain models of polarization. Specifically, we look at several spatial stochastic models of Cdc42 polarization that have been adapted from published models, on a variety of tip-shaped geometries, to replicate the shape change that occurs during the growth of the mating projection. We show here that there is a complex interplay between the dynamics of polarization and the shape of the cell. Our results show that while models of polarization can generate a stable polarization cap, its localization at the tip of mating projections is unstable, with the polarization cap drifting away from the tip of the projection in a geometry dependent manner. We also compare predictions from our computational results to experiments that observe cells with projections of varying lengths, and track the stability of the polarization cap. Lastly, we examine one model of actin polarization and show that it is unlikely, at least for the models studied here, that actin dynamics and vesicle traffic are able to overcome this effect of geometry.

MinE conformational switching confers robustness on self-organized Min protein patterns

Many fundamental cellular processes are spatially regulated by self-organized protein patterns, which are often based on nucleotide-binding proteins that switch their nucleotide state upon interaction with a second, activating protein. For reliable function, these protein patterns must be robust against parameter changes, although the basis for such robustness is generally elusive. Here we take a combined theoretical and experimental approach to the Escherichia coli Min system, a paradigmatic system for protein self-organization. By mathematical modeling and in vitro reconstitution of mutant proteins, we demonstrate that the robustness of pattern formation is dramatically enhanced by an interlinked functional switching of both proteins, rather than one. Such interlinked functional switching could be a generic means of obtaining robustness in biological pattern-forming systems.

Protein patterning is vital for many fundamental cellular processes. This raises two intriguing questions: Can such intrinsically complex processes be reduced to certain core principles and, if so, what roles do the molecular details play in individual systems? A prototypical example for protein patterning is the bacterial Min system, in which self-organized pole-to-pole oscillations of MinCDE proteins guide the cell division machinery to midcell. These oscillations are based on cycling of the ATPase MinD and its activating protein MinE between the membrane and the cytoplasm. Recent biochemical evidence suggests that MinE undergoes a reversible, MinD-dependent conformational switch from a latent to a reactive state. However, the functional relevance of this switch for the Min network and pattern formation remains unclear. By combining mathematical modeling and in vitro reconstitution of mutant proteins, we dissect the two aspects of MinE’s switch, persistent membrane binding and a change in MinE’s affinity for MinD. Our study shows that the MinD-dependent change in MinE’s binding affinity for MinD is essential for patterns to emerge over a broad and physiological range of protein concentrations. Mechanistically, our results suggest that conformational switching of an ATPase-activating protein can lead to the spatial separation of its distinct functional states and thereby confer robustness on an intracellular protein network with vital roles in bacterial cell division.

Spatial modeling of the membrane-cytosolic interface in protein kinase signal transduction

The spatial architecture of signaling pathways and the interaction with cell size and morphology are complex, but little understood. With the advances of single cell imaging and single cell biology, it becomes crucial to understand intracellular processes in time and space. Activation of cell surface receptors often triggers a signaling cascade including the activation of membrane-attached and cytosolic signaling components, which eventually transmit the signal to the cell nucleus. Signaling proteins can form steep gradients in the cytosol, which cause strong cell size dependence. We show that the kinetics at the membrane-cytosolic interface and the ratio of cell membrane area to the enclosed cytosolic volume change the behavior of signaling cascades significantly. We suggest an estimate of average concentration for arbitrary cell shapes depending on the cell volume and cell surface area. The normalized variance, known from image analysis, is suggested as an alternative measure to quantify the deviation from the average concentration. A mathematical analysis of signal transduction in time and space is presented, providing analytical solutions for different spatial arrangements of linear signaling cascades. Quantification of signaling time scales reveals that signal propagation is faster at the membrane than at the nucleus, while this time difference decreases with the number of signaling components in the cytosol. Our investigations are complemented by numerical simulations of non-linear cascades with feedback and asymmetric cell shapes. We conclude that intracellular signal propagation is highly dependent on cell geometry and, thereby, conveys information on cell size and shape to the nucleus.

Particle-based simulations of polarity establishment reveal stochastic promotion of Turing pattern formation

Polarity establishment, the spontaneous generation of asymmetric molecular distributions, is a crucial component of many cellular functions. Saccharomyces cerevisiae (yeast) undergoes directed growth during budding and mating, and is an ideal model organism for studying polarization. In yeast and many other cell types, the Rho GTPase Cdc42 is the key molecular player in polarity establishment. During yeast polarization, multiple patches of Cdc42 initially form, then resolve into a single front. Because polarization relies on strong positive feedback, it is likely that the amplification of molecular-level fluctuations underlies the generation of multiple nascent patches. In the absence of spatial cues, these fluctuations may be key to driving polarization. Here we used particle-based simulations to investigate the role of stochastic effects in a Turing-type model of yeast polarity establishment. In the model, reactions take place either between two molecules on the membrane, or between a cytosolic and a membrane-bound molecule. Thus, we developed a computational platform that explicitly simulates molecules at and near the cell membrane, and implicitly handles molecules away from the membrane. To evaluate stochastic effects, we compared particle simulations to deterministic reaction-diffusion equation simulations. Defining macroscopic rate constants that are consistent with the microscopic parameters for this system is challenging, because diffusion occurs in two dimensions and particles exchange between the membrane and cytoplasm. We address this problem by empirically estimating macroscopic rate constants from appropriately designed particle-based simulations. Ultimately, we find that stochastic fluctuations speed polarity establishment and permit polarization in parameter regions predicted to be Turing stable. These effects can operate at Cdc42 abundances expected of yeast cells, and promote polarization on timescales consistent with experimental results. To our knowledge, our work represents the first particle-based simulations of a model for yeast polarization that is based on a Turing mechanism.

Many cells need to generate and maintain biochemical signals in specific subcellular regions. This phenomenon is broadly called polarity establishment, and is important in fundamental processes such as cell migration and differentiation. A key polarity factor found in diverse organisms, including yeast and humans, is the protein Cdc42. In yeast, Cdc42-dependent polarization occurs through a self-reinforcing biochemical signaling loop. Directional cues can guide polarity establishment, but interestingly, yeast can polarize in the absence of such a cue. The mechanism thought to underlie this symmetry breaking involves the amplification of inhomogeneities in molecular distributions that arise from molecular-level fluctuations. We investigated the effects of random fluctuations on polarization by performing particle-based simulations of the Cdc42 signaling network. We found that fluctuations can facilitate polarization, allowing faster polarization, and polarization over a broader range of concentrations. Our observations may help understand how polarity works in other systems.

Experimental evolution and proximate mechanisms in biology

Biological functions – studied by molecular, systems and behavioral biology – are referred to as proximate mechanisms. Why and how they have emerged from the course of evolution are referred to as ultimate mechanisms. Despite the conceptual and technical schism between the disciplines that focus on each, studies from one side can benefit the other. Experimental evolution is an emerging field at the crossroads of functional and evolutionary biology. Herein microorganisms and mammalian cell lines evolve in well-controlled laboratory environments over multiple generations. Phenotypic changes arising from the process are then characterized in genetics and function to understand the evolutionary process. While providing empirical tests to evolutionary questions, such studies also offer opportunities of new insights into proximate mechanisms. Experimental evolution optimizes biological systems by means of adaptation; the adapted systems with their mutations present unique perturbed states of the systems that generate new and often unexpected output/performance. Hence, learning about these states not only adds to but also might deepen knowledge on the proximate processes. To demonstrate this point, five examples in experimental evolution are introduced, and their relevance to functional biology explicated. In some examples, from evolution experiments, updates were made to known proximate processes – gene regulation and cell polarization. In some examples, new contexts were found for known proximate processes – cell division and drug resistance of cancer. In one example, a new cellular mechanism was discovered. These cases identify ways the approach of experimental evolution can be used to ask questions in functional biology.

Many roads to symmetry breaking: molecular mechanisms and theoretical models of yeast cell polarity

Mathematical modeling has been instrumental in identifying common principles of cell polarity across diverse systems. These principles include positive feedback loops that are required to destabilize a spatially uniform state of the cell. The conserved small G-protein Cdc42 is a master regulator of eukaryotic cellular polarization. Here we discuss recent developments in studies of Cdc42 polarization in budding and fission yeasts and demonstrate that models describing symmetry-breaking polarization can be classified into six minimal classes based on the structure of positive feedback loops that activate and localize Cdc42. Owing to their generic system-independent nature, these model classes are also likely to be relevant for the G-protein–based symmetry-breaking systems of higher eukaryotes. We review experimental evidence pro et contra different theoretically plausible models and conclude that several parallel and non–mutually exclusive mechanisms are likely involved in cellular polarization of yeasts. This potential redundancy needs to be taken into consideration when interpreting the results of recent cell-rewiring studies.

A microtubule-based minimal model for spontaneous and persistent spherical cell polarity

We propose a minimal model for the spontaneous and persistent generation of polarity in a spherical cell based on dynamic microtubules and a single mobile molecular component. This component, dubbed the polarity factor, binds to microtubules nucleated from a centrosome located in the center of the cell, is subsequently delivered to the cell membrane, where it diffuses until it unbinds. The only feedback mechanism we impose is that the residence time of the microtubules at the membrane increases with the local density of the polarity factor. We show analytically that this system supports a stable unipolar symmetry-broken state for a wide range of parameters. We validate the predictions of the model by 2D particle-based simulations. Our model provides a route towards the creation of polarity in a minimal cell-like environment using a biochemical reconstitution approach.

Evolutionary dynamics in the fungal polarization network, a mechanistic perspective

Polarity establishment underlies proper cell cycle completion across virtually all organisms. Much progress has been made in generating an understanding of the structural and functional components of this process, especially in model species. Here we focus on the evolutionary dynamics of the fungal polarization protein network in order to determine general components and mechanistic principles, species- or lineage-specific adaptations and the evolvability of the network. The currently available genomic and proteomic screens in a variety of fungal species have shown three main characteristics: (1) certain proteins, processes and functions are conserved throughout the fungal clade; (2) orthologous functions can never be assumed, as various cases have been observed of homologous loci with dissimilar functions; (3) species have, typically, various species- or lineage-specific proteins incorporated in their polarization network. Further large-scale comparative and experimental studies, including those on non-model species representing the great fungal diversity, are needed to gain a better understanding of the evolutionary dynamics and generalities of the polarization network in fungi.

Local sampling paints a global picture: Local concentration measurements sense direction in complex chemical gradients

Detecting and interpreting extracellular spatial signals is essential for cellular orientation within complex environments, such as during directed cell migration or growth in multicellular development. Although the molecular understanding of how cells read spatial signals like chemical gradients is still lacking, recent work has revealed that stochastic processes at different temporal and spatial scales are at the core of this gradient sensing process in a wide range of eukaryotes. Fast biochemical reactions like those underlying GTPase activity dynamics form a functional module together with slower cell morphological changes driven by membrane remodelling. This biochemical-morphological module explores the environment by stochastic local concentration sampling to determine the source of the gradient signal, enabling efficient signal detection and interpretation before polarised growth or migration towards the gradient source is initiated. Here we review recent data describing local sampling and propose a model of local fast and slow feedback counteracted by gradient-dependent substrate limitation to be at the core of gradient sensing by local sampling.

A Two-Tiered Mechanism Enables Localized Cdc42 Signaling during Enterocyte Polarization

Signaling by the small GTPase Cdc42 governs a diverse set of cellular processes that contribute to tissue morphogenesis. Since these processes often require highly localized signaling, Cdc42 activity must be clustered in order to prevent ectopic signaling. During cell polarization, apical Cdc42 signaling directs the positioning of the nascent apical membrane. However, the molecular mechanisms that drive Cdc42 clustering during polarity establishment are largely unknown. Here, we demonstrate that during cell polarization localized Cdc42 signaling is enabled via activity-dependent control of Cdc42 mobility. By performing photoconversion experiments, we show that inactive Cdc42-GDP is 30-fold more mobile than active Cdc42-GTP. This switch in apical mobility originates from a dual mechanism involving RhoGDI-mediated membrane dissociation of Cdc42-GDP and Tuba-mediated immobilization of Cdc42-GTP. Interference with either mechanism affects Cdc42 clustering and as a consequence impairs Cdc42-mediated apical membrane clustering. We therefore identify a molecular network, comprised of Cdc42, the guanine nucleotide exchange factor (GEF) Tuba, and RhoGDI, that enables differential diffusion of inactive and active Cdc42 and is required to establish localized Cdc42 signaling during enterocyte polarization.

SParSE++: improved event-based stochastic parameter search

Background

Despite the increasing availability of high performance computing capabilities, analysis and characterization of stochastic biochemical systems remain a computational challenge. To address this challenge, the Stochastic Parameter Search for Events (SParSE) was developed to automatically identify reaction rates that yield a probabilistic user-specified event. SParSE consists of three main components: the multi-level cross-entropy method, which identifies biasing parameters to push the system toward the event of interest, the related inverse biasing method, and an optional interpolation of identified parameters. While effective for many examples, SParSE depends on the existence of a sufficient amount of intrinsic stochasticity in the system of interest. In the absence of this stochasticity, SParSE can either converge slowly or not at all.

Results

We have developed SParSE++, a substantially improved algorithm for characterizing target events in terms of system parameters. SParSE++ makes use of a series of novel parameter leaping methods that accelerate the convergence rate to the target event, particularly in low stochasticity cases. In addition, the interpolation stage is modified to compute multiple interpolants and to choose the optimal one in a statistically rigorous manner. We demonstrate the performance of SParSE++ on four example systems: a birth-death process, a reversible isomerization model, SIRS disease dynamics, and a yeast polarization model. In all four cases, SParSE++ shows significantly improved computational efficiency over SParSE, with the largest improvements resulting from analyses with the strictest error tolerances.

Conclusions

As researchers continue to model realistic biochemical systems, the need for efficient methods to characterize target events will grow. The algorithmic advancements provided by SParSE++ fulfill this need, enabling characterization of computationally intensive biochemical events that are currently resistant to analysis.

Electronic supplementary material

The online version of this article (doi:10.1186/s12918-016-0367-z) contains supplementary material, which is available to authorized users.

A framework for discrete stochastic simulation on 3D moving boundary domains

We have developed a method for modeling spatial stochastic biochemical reactions in complex, three-dimensional, and time-dependent domains using the reaction-diffusion master equation formalism. In particular, we look to address the fully coupled problems that arise in systems biology where the shape and mechanical properties of a cell are determined by the state of the biochemistry and vice versa. To validate our method and characterize the error involved, we compare our results for a carefully constructed test problem to those of a microscale implementation. We demonstrate the effectiveness of our method by simulating a model of polarization and shmoo formation during the mating of yeast. The method is generally applicable to problems in systems biology where biochemistry and mechanics are coupled, and spatial stochastic effects are critical.

Multistability and dynamic transitions of intracellular Min protein patterns

Cells owe their internal organization to self-organized protein patterns, which originate and adapt to growth and external stimuli via a process that is as complex as it is little understood. Here, we study the emergence, stability, and state transitions of multistable Min protein oscillation patterns in live Escherichia coli bacteria during growth up to defined large dimensions. De novo formation of patterns from homogenous starting conditions is observed and studied both experimentally and in simulations. A new theoretical approach is developed for probing pattern stability under perturbations. Quantitative experiments and simulations show that, once established, Min oscillations tolerate a large degree of intracellular heterogeneity, allowing distinctly different patterns to persist in different cells with the same geometry. Min patterns maintain their axes for hours in experiments, despite imperfections, expansion, and changes in cell shape during continuous cell growth. Transitions between multistable Min patterns are found to be rare events induced by strong intracellular perturbations. The instances of multistability studied here are the combined outcome of boundary growth and strongly nonlinear kinetics, which are characteristic of the reaction-diffusion patterns that pervade biology at many scales.

Geometry-induced protein pattern formation

Protein patterns are known to adapt to cell shape and serve as spatial templates that choreograph downstream processes like cell polarity or cell division. However, how can pattern-forming proteins sense and respond to the geometry of a cell, and what mechanistic principles underlie pattern formation? Current models invoke mechanisms based on dynamic instabilities arising from nonlinear interactions between proteins but neglect the influence of the spatial geometry itself. Here, we show that patterns can emerge as a direct result of adaptation to cell geometry, in the absence of dynamical instability. We present a generic reaction module that allows protein densities robustly to adapt to the symmetry of the spatial geometry. The key component is an NTPase protein that cycles between nucleotide-dependent membrane-bound and cytosolic states. For elongated cells, we find that the protein dynamics generically leads to a bipolar pattern, which vanishes as the geometry becomes spherically symmetrical. We show that such a reaction module facilitates universal adaptation to cell geometry by sensing the local ratio of membrane area to cytosolic volume. This sensing mechanism is controlled by the membrane affinities of the different states. We apply the theory to explain AtMinD bipolar patterns in [Formula: see text] EcMinDE Escherichia coli. Due to its generic nature, the mechanism could also serve as a hitherto-unrecognized spatial template in many other bacterial systems. Moreover, the robustness of the mechanism enables self-organized optimization of protein patterns by evolutionary processes. Finally, the proposed module can be used to establish geometry-sensitive protein gradients in synthetic biological systems.

ATP8B1-mediated spatial organization of Cdc42 signaling maintains singularity during enterocyte polarization

The disease-associated phospholipid flippase ATP8B1 decreases Cdc42 mobility at the apical membrane to ensure the formation of a single apical domain and to maintain healthy lumen architecture.

During yeast cell polarization localization of the small GTPase, cell division control protein 42 homologue (Cdc42) is clustered to ensure the formation of a single bud. Here we show that the disease-associated flippase ATPase class I type 8b member 1 (ATP8B1) enables Cdc42 clustering during enterocyte polarization. Loss of this regulation results in increased apical membrane size with scattered apical recycling endosomes and permits the formation of more than one apical domain, resembling the singularity defect observed in yeast. Mechanistically, we show that to become apically clustered, Cdc42 requires the interaction between its polybasic region and negatively charged membrane lipids provided by ATP8B1. Disturbing this interaction, either by ATP8B1 depletion or by introduction of a Cdc42 mutant defective in lipid binding, increases Cdc42 mobility and results in apical membrane enlargement. Re-establishing Cdc42 clustering, by tethering it to the apical membrane or lowering its diffusion, restores normal apical membrane size in ATP8B1-depleted cells. We therefore conclude that singularity regulation by Cdc42 is conserved between yeast and human and that this regulation is required to maintain healthy tissue architecture.

A Cellular System for Spatial Signal Decoding in Chemical Gradients

Directional cell growth requires that cells read and interpret shallow chemical gradients, but how the gradient directional information is identified remains elusive. We use single-cell analysis and mathematical modeling to define the cellular gradient decoding network in yeast. Our results demonstrate that the spatial information of the gradient signal is read locally within the polarity site complex using double-positive feedback between the GTPase Cdc42 and trafficking of the receptor Ste2. Spatial decoding critically depends on low Cdc42 activity, which is maintained by the MAPK Fus3 through sequestration of the Cdc42 activator Cdc24. Deregulated Cdc42 or Ste2 trafficking prevents gradient decoding and leads to mis-oriented growth. Our work discovers how a conserved set of components assembles a network integrating signal intensity and directionality to decode the spatial information contained in chemical gradients.

Dendritic spine geometry can localize GTPase signaling in neurons

Dendritic spines are the postsynaptic terminals of most excitatory synapses in the mammalian brain. Learning and memory are associated with long-lasting structural remodeling of dendritic spines through an actin-mediated process regulated by the Rho-family GTPases RhoA, Rac, and Cdc42. These GTPases undergo sustained activation after synaptic stimulation, but whereas Rho activity can spread from the stimulated spine, Cdc42 activity remains localized to the stimulated spine. Because Cdc42 itself diffuses rapidly in and out of the spine, the basis for the retention of Cdc42 activity in the stimulated spine long after synaptic stimulation has ceased is unclear. Here we model the spread of Cdc42 activation at dendritic spines by means of reaction-diffusion equations solved on spine-like geometries. Excitable behavior arising from positive feedback in Cdc42 activation leads to spreading waves of Cdc42 activity. However, because of the very narrow neck of the dendritic spine, wave propagation is halted through a phenomenon we term geometrical wave-pinning. We show that this can account for the localization of Cdc42 activity in the stimulated spine, and, of interest, retention is enhanced by high diffusivity of Cdc42. Our findings are broadly applicable to other instances of signaling in extreme geometries, including filopodia and primary cilia.

Role of competition between polarity sites in establishing a unique front

Polarity establishment in many cells is thought to occur via positive feedback that reinforces even tiny asymmetries in polarity protein distribution. Cdc42 and related GTPases are activated and accumulate in a patch of the cortex that defines the front of the cell. Positive feedback enables spontaneous polarization triggered by stochastic fluctuations, but as such fluctuations can occur at multiple locations, how do cells ensure that they make only one front? In polarizing cells of the model yeast Saccharomyces cerevisiae, positive feedback can trigger growth of several Cdc42 clusters at the same time, but this multi-cluster stage rapidly evolves to a single-cluster state, which then promotes bud emergence. By manipulating polarity protein dynamics, we show that resolution of multi-cluster intermediates occurs through a greedy competition between clusters to recruit and retain polarity proteins from a shared intracellular pool.