A coupled bulk-surface model for cell polarisation

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.

Physical confinement and cell proximity increase cell migration rates and invasiveness: A mathematical model of cancer cell invasion through flexible channels

Cancer cell migration between different body parts is the driving force behind cancer metastasis, which is the main cause of mortality of patients. Migration of cancer cells often proceeds by penetration through narrow cavities in locally stiff, yet flexible tissues. In our previous work, we developed a model for cell geometry evolution during invasion, which we extend here to investigate whether leader and follower (cancer) cells that only interact mechanically can benefit from sequential transmigration through narrow micro-channels and cavities. We consider two cases of cells sequentially migrating through a flexible channel: leader and follower cells being closely adjacent or distant. Using Wilcoxon's signed-rank test on the data collected from Monte Carlo simulations, we conclude that the modelled transmigration speed for the follower cell is significantly larger than for the leader cell when cells are distant, i.e. follower cells transmigrate after the leader has completed the crossing. Furthermore, it appears that there exists an optimum with respect to the width of the channel such that cell moves fastest. On the other hand, in the case of closely adjacent cells, effectively performing collective migration, the leader cell moves 12% faster since the follower cell pushes it. This work shows that mechanical interactions between cells can increase the net transmigration speed of cancer cells, resulting in increased invasiveness. In other words, interaction between cancer cells can accelerate metastatic invasion.

The emergence of spatial patterns for compartmental reaction kinetics coupled by two bulk diffusing species with comparable diffusivities

Originating from the pioneering study of Alan Turing, the bifurcation analysis predicting spatial pattern formation from a spatially uniform state for diffusing morphogens or chemical species that interact through nonlinear reactions is a central problem in many chemical and biological systems. From a mathematical viewpoint, one key challenge with this theory for two component systems is that stable spatial patterns can typically only occur from a spatially uniform state when a slowly diffusing 'activator' species reacts with a much faster diffusing 'inhibitor' species. However, from a modelling perspective, this large diffusivity ratio requirement for pattern formation is often unrealistic in biological settings since different molecules tend to diffuse with similar rates in extracellular spaces. As a result, one key long-standing question is how to robustly obtain pattern formation in the biologically realistic case where the time scales for diffusion of the interacting species are comparable. For a coupled one-dimensional bulk-compartment theoretical model, we investigate the emergence of spatial patterns for the scenario where two bulk diffusing species with comparable diffusivities are coupled to nonlinear reactions that occur only in localized 'compartments', such as on the boundaries of a one-dimensional domain. The exchange between the bulk medium and the spatially localized compartments is modelled by a Robin boundary condition with certain binding rates. As regulated by these binding rates, we show for various specific nonlinearities that our one-dimensional coupled PDE-ODE model admits symmetry-breaking bifurcations, leading to linearly stable asymmetric steady-state patterns, even when the bulk diffusing species have equal diffusivities. Depending on the form of the nonlinear kinetics, oscillatory instabilities can also be triggered. Moreover, the analysis is extended to treat a periodic chain of compartments. This article is part of the theme issue 'New trends in pattern formation and nonlinear dynamics of extended systems'.

Front propagation in the shadow wave-pinning model

In this paper we consider a non-local bistable reaction-diffusion equation which is a simplified version of the wave-pinning model of cell polarization. In the small diffusion limit, a typical solution u(x, t) of this model approaches one of the stable states of the bistable nonlinearity in different parts of the spatial domain [Formula: see text], separated by an interface moving at a normal velocity regulated by the integral [Formula: see text]. In what is often referred to as wave-pinning, feedback between mass-conservation and bistablity causes the interface to slow and approach a fixed limit. In the limit of a small diffusivity [Formula: see text], we prove that for any [Formula: see text] the interface can be estimated within [Formula: see text] of the location as predicted using formal asymptotics. We also discuss the sharpness of our result by comparing the formal asymptotic results with numerical simulations.

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.

Numerical investigations of the bulk-surface wave pinning model

The bulk-surface wave pinning model is a reaction-diffusion system for studying cell polarisation. It is constituted by a surface reaction-diffusion equation, coupled to a bulk diffusion equation with a non-linear boundary condition. Cell polarisation arises as the surface component develops specific patterns. Since proteins diffuse much faster in the cell interior than on the membrane, in the literature, the bulk component is often assumed to be spatially homogeneous. Therefore, the model can be reduced to a single surface equation. However, in real applications a spatially non-uniform bulk component might be an important player to take into account. In this paper, we study, through numerical computations, the role of the bulk component and, more specifically, how different bulk diffusion rates might affect the polarisation response. We find that the bulk component is indeed a key factor in determining the surface polarisation response. Moreover, for certain geometries, it is the spatial heterogeneity of the bulk component that triggers the polarisation response, which might not be possible in a reduced model. Understanding how polarisation depends on bulk diffusivity might be crucial when studying models of migrating cells, which are naturally subject to domain deformation.

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.

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.

Computational framework for single-cell spatiotemporal dynamics of optogenetic membrane recruitment

We describe a modular computational framework for analyzing cell-wide spatiotemporal signaling dynamics in single-cell microscopy experiments that accounts for the experiment-specific geometric and diffractive complexities that arise from heterogeneous cell morphologies and optical instrumentation. Inputs are unique cell geometries and protein concentrations derived from confocal stacks and spatiotemporally varying environmental stimuli. After simulating the system with a model of choice, the output is convolved with the microscope point-spread function for direct comparison with the observable image. We experimentally validate this approach in single cells with BcLOV4, an optogenetic membrane recruitment system for versatile control over cell signaling, using a three-dimensional non-linear finite element model with all parameters experimentally derived. The simulations recapitulate observed subcellular and cell-to-cell variability in BcLOV4 signaling, allowing for inter-experimental differences of cellular and instrumentation origins to be elucidated and resolved for improved interpretive robustness. This single-cell approach will enhance optogenetics and spatiotemporally resolved signaling studies.

Pattern forming systems coupling linear bulk diffusion to dynamically active membranes or cells

Some analytical and numerical results are presented for pattern formation properties associated with novel types of reaction-diffusion (RD) systems that involve the coupling of bulk diffusion in the interior of a multi-dimensional spatial domain to nonlinear processes that occur either on the domain boundary or within localized compartments that are confined within the domain. The class of bulk-membrane system considered herein is derived from an asymptotic analysis in the limit of small thickness of a thin domain that surrounds the bulk medium. When the bulk domain is a two-dimensional disk, a weakly nonlinear analysis is used to characterize Turing and Hopf bifurcations that can arise from the linearization around a radially symmetric, but spatially non-uniform, steady-state of the bulk-membrane system. In a singularly perturbed limit, the existence and linear stability of localized membrane-bound spike patterns is analysed for a Gierer-Meinhardt activator-inhibitor model that includes bulk coupling. Finally, the emergence of collective intracellular oscillations is studied for a class of PDE-ODE bulk-cell model in a bounded two-dimensional domain that contains spatially localized, but dynamically active, circular cells that are coupled through a linear bulk diffusion field. Applications of such coupled bulk-membrane or bulk-cell systems to some biological systems are outlined, and some open problems in this area are discussed. This article is part of the theme issue 'Recent progress and open frontiers in Turing's theory of morphogenesis'.

Modern perspectives on near-equilibrium analysis of Turing systems

In the nearly seven decades since the publication of Alan Turing's work on morphogenesis, enormous progress has been made in understanding both the mathematical and biological aspects of his proposed reaction-diffusion theory. Some of these developments were nascent in Turing's paper, and others have been due to new insights from modern mathematical techniques, advances in numerical simulations and extensive biological experiments. Despite such progress, there are still important gaps between theory and experiment, with many examples of biological patterning where the underlying mechanisms are still unclear. Here, we review modern developments in the mathematical theory pioneered by Turing, showing how his approach has been generalized to a range of settings beyond the classical two-species reaction-diffusion framework, including evolving and complex manifolds, systems heterogeneous in space and time, and more general reaction-transport equations. While substantial progress has been made in understanding these more complicated models, there are many remaining challenges that we highlight throughout. We focus on the mathematical theory, and in particular linear stability analysis of 'trivial' base states. We emphasize important open questions in developing this theory further, and discuss obstacles in using these techniques to understand biological reality. This article is part of the theme issue 'Recent progress and open frontiers in Turing's theory of morphogenesis'.

Narrow escape problem in two-shell spherical domains

Intracellular transport in living cells is often spatially inhomogeneous with an accelerated effective diffusion close to the cell membrane and a ballistic motion away from the centrosome due to active transport along actin filaments and microtubules, respectively. Recently it was reported that the mean first passage time (MFPT) for transport to a specific area on the cell membrane is minimal for an optimal actin cortex width. In this paper, we ask whether this optimization in a two-compartment domain can also be achieved by passive Brownian particles. We consider a Brownian motion with different diffusion constants in the two shells and a potential barrier between the two, and we investigate the narrow escape problem by calculating the MFPT for Brownian particles to reach a small window on the external boundary. In two and three dimensions, we derive asymptotic expressions for the MFPT in the thin cortex and small escape region limits confirmed by numerical calculations of the MFPT using the finite-element method and stochastic simulations. From this analytical and numeric analysis, we finally extract the dependence of the MFPT on the ratio of diffusion constants, the potential barrier height, and the width of the outer shell. The first two are monotonous, whereas the last one may have a minimum for a sufficiently attractive cortex, for which we propose an analytical expression of the potential barrier height matching very well the numerical predictions.

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.

Signaling gradients in surface dynamics as basis for planarian regeneration

Based on experimental data, we introduce and analyze a system of reaction-diffusion equations for the regeneration of planarian flatworms. We model dynamics of head and tail cells expressing positional control genes that translate into localized signals which in turn guide stem cell differentiation. Tissue orientation and positional information are encoded in a long range wnt-related signaling gradient. Our system correctly reproduces typical cut and graft experiments, and improves on previous models by preserving polarity in regeneration over orders of magnitude in body size during growth phases. Key to polarity preservation in our model flatworm is the sensitivity of cell differentiation to gradients of wnt-related signals relative to the tissue surface. This process is particularly relevant in small tissue layers close to cuts during their healing, and modeled in a robust fashion through dynamic boundary conditions.

Long time behavior and stable patterns in high-dimensional polarity models of asymmetric cell division

Asymmetric cell division is one of the fundamental processes to create cell diversity in the early stage of embryonic development. During this process, the polarity formation in the cell membrane has been considered as a key process by which the entire polarity formation in the cytosol is controlled, and it has been extensively studied in both experiments and mathematical models. Nonetheless, a mathematically rigorous analysis of the polarity formation in the asymmetric cell division has been little explored, particularly for bulk-surface models. In this article, we deal with polarity models proposed for describing the PAR polarity formation in the asymmetric cell division of a C. elegans embryo. Using a simpler but mathematically consistent model, we exhibit the long time behavior of the polarity formation of a bulk-surface cell. Moreover, we mathematically prove the existence of stable polarity solutions of the model equation in an arbitrary high-dimensional domain and analyse how the boundary position of polarity domain is determined. Our results propose that the existence and dynamics of the polarity in the asymmetric cell division can be understood universally in terms of basic mathematical structures.

On the influence of cell shape on dynamic reaction-diffusion polarization patterns

The distribution of signaling molecules following mechanical or chemical stimulation of a cell defines cell polarization, with regions of high active Cdc42 at the front and low active Cdc42 at the rear. As reaction-diffusion phenomena between signaling molecules, such as Rho GTPases, define the gradient dynamics, we hypothesize that the cell shape influences the maintenance of the “front-to-back” cell polarization patterns. We investigated the influence of cell shape on the Cdc42 patterns using an established computational polarization model. Our simulation results showed that not only cell shape but also Cdc42 and Rho-related (in)activation parameter values affected the distribution of active Cdc42. Despite an initial Cdc42 gradient, the in silico results showed that the maximal Cdc42 concentration shifts in the opposite direction, a phenomenon we propose to call “reverse polarization”. Additional in silico analyses indicated that “reverse polarization” only occurred in a particular parameter value space that resulted in a balance between inactivation and activation of Rho GTPases. Future work should focus on a mathematical description of the underpinnings of reverse polarization, in combination with experimental validation using, for example, dedicated FRET-probes to spatiotemporally track Rho GTPase patterns in migrating cells. In summary, the findings of this study enhance our understanding of the role of cell shape in intracellular signaling.

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.

Mathematical modelling in cell migration: tackling biochemistry in changing geometries

Directed cell migration poses a rich set of theoretical challenges. Broadly, these are concerned with (1) how cells sense external signal gradients and adapt; (2) how actin polymerisation is localised to drive the leading cell edge and Myosin-II molecular motors retract the cell rear; and (3) how the combined action of cellular forces and cell adhesion results in cell shape changes and net migration. Reaction-diffusion models for biological pattern formation going back to Turing have long been used to explain generic principles of gradient sensing and cell polarisation in simple, static geometries like a circle. In this minireview, we focus on recent research which aims at coupling the biochemistry with cellular mechanics and modelling cell shape changes. In particular, we want to contrast two principal modelling approaches: (1) interface tracking where the cell membrane, interfacing cell interior and exterior, is explicitly represented by a set of moving points in 2D or 3D space and (2) interface capturing. In interface capturing, the membrane is implicitly modelled analogously to a level line in a hilly landscape whose topology changes according to forces acting on the membrane. With the increased availability of high-quality 3D microscopy data of complex cell shapes, such methods will become increasingly important in data-driven, image-based modelling to better understand the mechanochemistry underpinning cell motion.

Bridging from single to collective cell migration: A review of models and links to experiments

Mathematical and computational models can assist in gaining an understanding of cell behavior at many levels of organization. Here, we review models in the literature that focus on eukaryotic cell motility at 3 size scales: intracellular signaling that regulates cell shape and movement, single cell motility, and collective cell behavior from a few cells to tissues. We survey recent literature to summarize distinct computational methods (phase-field, polygonal, Cellular Potts, and spherical cells). We discuss models that bridge between levels of organization, and describe levels of detail, both biochemical and geometric, included in the models. We also highlight links between models and experiments. We find that models that span the 3 levels are still in the minority.

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.

CDC-42 Interactions with Par Proteins Are Critical for Proper Patterning in Polarization

Many cells rearrange proteins and other components into spatially distinct domains in a process called polarization. This asymmetric patterning is required for a number of biological processes including asymmetric division, cell migration, and embryonic development. Proteins involved in polarization are highly conserved and include members of the Par and Rho protein families. Despite the importance of these proteins in polarization, it is not yet known how they interact and regulate each other to produce the protein localization patterns associated with polarization. In this study, we develop and analyse a biologically based mathematical model of polarization that incorporates interactions between Par and Rho proteins that are consistent with experimental observations of CDC-42. Using minimal network and eFAST sensitivity analyses, we demonstrate that CDC-42 is predicted to reinforce maintenance of anterior PAR protein polarity which in turn feedbacks to maintain CDC-42 polarization, as well as supporting posterior PAR protein polarization maintenance. The mechanisms for polarity maintenance identified by these methods are not sufficient for the generation of polarization in the absence of cortical flow. Additional inhibitory interactions mediated by the posterior Par proteins are predicted to play a role in the generation of Par protein polarity. More generally, these results provide new insights into the role of CDC-42 in polarization and the mutual regulation of key polarity determinants, in addition to providing a foundation for further investigations.

Influence of myosin activity and mechanical impact on keratocyte polarization

In cell migration, polarization is the process by which a stationary cell breaks symmetry and initiates motion. Although a lot is known about the mechanisms involved in cell polarization, the role played by myosin contraction remains unclear. In addition, cell polarization by mechanical impact has received little attention. Here, we study the influence of myosin activity on cell polarization and the initiation of motion induced by mechanical cues using a computational model for keratocytes. The model accounts for cell deformation, the dynamics of myosin and the signaling protein RhoA (a member of the Rho GTPases family), as well as the forces acting on the actomyosin network. Our results show that the attainment of a steady polarized state depends on the strength of myosin down- or up-regulation and that myosin upregulation favors cell polarization. Our results also confirm the existence of a threshold level for cell polarization, which is determined by the level of polarization of the Rho GTPases at the time the external stimuli vanish. In all, this paper shows that capturing the interactions between the signaling proteins (Rho GTPases for keratocytes) and the compounds of the motile machinery in a moving cell is crucial to study cell polarization.

Simple Rho GTPase Dynamics Generate a Complex Regulatory Landscape Associated with Cell Shape

Migratory cells exhibit a variety of morphologically distinct responses to their environments that manifest in their cell shape. Some protrude uniformly to increase substrate contacts, others are broadly contractile, some polarize to facilitate migration, and yet others exhibit mixtures of these responses. Prior studies have identified a discrete collection of shapes that the majority of cells display and demonstrated that activity levels of the cytoskeletal regulators Rac1 and RhoA GTPase regulate those shapes. Here, we use computational modeling to assess whether known GTPase dynamics can give rise to a sufficient diversity of spatial signaling states to explain the observed shapes. Results show that the combination of autoactivation and mutually antagonistic cross talk between GTPases, along with the conservative membrane binding, generates a wide array of distinct homogeneous and polarized regulatory phenotypes that arise for fixed model parameters. From a theoretical perspective, these results demonstrate that simple GTPase dynamics can generate complex multistability in which six distinct stable steady states (three homogeneous and three polarized) coexist for a fixed set of parameters, each of which naturally maps to an observed morphology. From a biological perspective, although we do not explicitly model the cytoskeleton or resulting cell morphologies, these results, along with prior literature linking GTPase activity to cell morphology, support the hypothesis that GTPase signaling dynamics can generate the broad morphological characteristics observed in many migratory cell populations. Further, the observed diversity may be the result of cells populating a complex morphological landscape generated by GTPase regulation rather than being the result of intrinsic cell-cell variation. These results demonstrate that Rho GTPases may have a central role in regulating the broad characteristics of cell shape (e.g., expansive, contractile, polarized, etc.) and that shape heterogeneity may be (at least partly) a reflection of the rich signaling dynamics regulating the cytoskeleton rather than intrinsic cell heterogeneity.

Cell Size, Mechanical Tension, and GTPase Signaling in the Single Cell

Cell polarization requires redistribution of specific proteins to the nascent front and back of a eukaryotic cell. Among these proteins are Rac and Rho, members of the small GTPase family that regulate the actin cytoskeleton. Rac promotes actin assembly and protrusion of the front edge, whereas Rho activates myosin-driven contraction at the back. Mathematical models of cell polarization at many levels of detail have appeared. One of the simplest based on "wave-pinning" consists of a pair of reaction-diffusion equations for a single GTPase. Mathematical analysis of wave-pinning so far is largely restricted to static domains in one spatial dimension. Here, we extend the analysis to cells that change in size, showing that both shrinking and growing cells can lose polarity. We further consider the feedback between mechanical tension, GTPase activation, and cell deformation in both static, growing, shrinking, and moving cells. Special cases (spatially uniform cell chemistry, the absence or presence of mechanical feedback) are analyzed, and the full model is explored by simulations in 1D. We find a variety of novel behaviors, including "dilution-induced" oscillations of Rac activity and cell size, as well as gain or loss of polarization and motility in the model cell.

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.