Local perturbation analysis: a computational tool for biophysical reaction-diffusion models

Turing Instabilities are Not Enough to Ensure Pattern Formation

Symmetry-breaking instabilities play an important role in understanding the mechanisms underlying the diversity of patterns observed in nature, such as in Turing's reaction-diffusion theory, which connects cellular signalling and transport with the development of growth and form. Extensive literature focuses on the linear stability analysis of homogeneous equilibria in these systems, culminating in a set of conditions for transport-driven instabilities that are commonly presumed to initiate self-organisation. We demonstrate that a selection of simple, canonical transport models with only mild multistable non-linearities can satisfy the Turing instability conditions while also robustly exhibiting only transient patterns. Hence, a Turing-like instability is insufficient for the existence of a patterned state. While it is known that linear theory can fail to predict the formation of patterns, we demonstrate that such failures can appear robustly in systems with multiple stable homogeneous equilibria. Given that biological systems such as gene regulatory networks and spatially distributed ecosystems often exhibit a high degree of multistability and nonlinearity, this raises important questions of how to analyse prospective mechanisms for self-organisation.

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

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

From actin waves to mechanism and back: How theory aids biological understanding

Actin dynamics in cell motility, division, and phagocytosis is regulated by complex factors with multiple feedback loops, often leading to emergent dynamic patterns in the form of propagating waves of actin polymerization activity that are poorly understood. Many in the actin wave community have attempted to discern the underlying mechanisms using experiments and/or mathematical models and theory. Here, we survey methods and hypotheses for actin waves based on signaling networks, mechano-chemical effects, and transport characteristics, with examples drawn from Dictyostelium discoideum, human neutrophils, Caenorhabditis elegans, and Xenopus laevis oocytes. While experimentalists focus on the details of molecular components, theorists pose a central question of universality: Are there generic, model-independent, underlying principles, or just boundless cell-specific details? We argue that mathematical methods are equally important for understanding the emergence, evolution, and persistence of actin waves and conclude with a few challenges for future studies.

Identifying patterning behavior in a plant infestation of insect pests

In this study, we developed a mechanistic model formulated as a system of reaction-diffusion equations (RDE) to explore the spatiotemporal dynamics of a theoretical pest with a tillering host plant in a controlled rectangular plant field. Local perturbation analysis, a recently developed method of analysis for wave propagation, was utilized to determine patterning regimes resulting from the local and global behaviors of the slow and fast diffusing components of the RDE system, respectively. Turing analysis was done to show that the RDE system does not exhibit Turing patterns. With bug mortality as the bifurcation parameter, regions with oscillations and stable coexistence of the pest and tillers were identified. Numerical simulations illustrate the patterning regimes in 1D and 2D settings. The oscillations suggest that recurrences in pest infestation is possible. Moreover, simulations showed that patterns produced in the model are strongly influenced by the pests' homogeneous dynamics inside the controlled environment.

Excitable dynamics in a molecularly-explicit model of cell motility: Mixed-mode oscillations and beyond

Mesenchymal cell motility is mainly regulated by two members of the Rho-family of GTPases, called Rac and Rho. The mutual inhibition exerted by these two proteins on each other's activation and the promotion of Rac activation by an adaptor protein called paxillin have been implicated in driving cellular polarization comprised of front (high active Rac) and back (high active Rho) during cell migration. Mathematical modeling of this regulatory network has previously shown that bistability is responsible for generating a spatiotemporal pattern underscoring cellular polarity called wave-pinning when diffusion is included. We previously developed a 6V reaction-diffusion model of this network to decipher the role of Rac, Rho and paxillin (along with other auxiliary proteins) in generating wave-pinning. In this study, we simplify this model through a series of steps into an excitable 3V ODE model comprised of one fast variable (the scaled concentration of active Rac), one slow variable (the maximum paxillin phosphorylation rate - turned into a variable) and a very slow variable (a recovery rate - also turned into a variable). We then explore, through slow-fast analysis, how excitability is manifested by showing that the model can exhibit relaxation oscillations (ROs) as well as mixed-mode oscillations (MMOs) whose underlying dynamics are consistent with a delayed Hopf bifurcation with a canard explosion. By reintroducing diffusion and the scaled concentration of inactive Rac into the model, we obtain a 4V PDE model that generates several unique spatiotemporal patterns that are relevant to cell motility. These patterns are then characterized and their impact on cell motility are explored by employing the cellular potts model (CPM). Our results reveal that wave pinning produces purely very directed motion in CPM, while MMOs allow for meandering and non-motile behaviors to occur. This highlights the role of MMOs as a potential mechanism for mesenchymal cell motility.

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

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

Characterizing spontaneous Ca2+ local transients in OPCs using computational modeling

Spontaneous Ca2+ local transients (SCaLTs) in isolated oligodendrocyte precursor cells are largely regulated by the following fluxes: store-operated Ca2+ entry (SOCE), Na+/Ca2+ exchange, Ca2+ pumping through Ca2+-ATPases, and Ca2+-induced Ca2+-release through ryanodine receptors and inositol-trisphosphate receptors. However, the relative contributions of these fluxes in mediating fast spiking and the slow baseline oscillations seen in SCaLTs remain incompletely understood. Here, we developed a stochastic spatiotemporal computational model to simulate SCaLTs in a homogeneous medium with ionic flow between the extracellular, cytoplasmic, and endoplasmic-reticulum compartments. By simulating the model and plotting both the histograms of SCaLTs obtained experimentally and from the model as well as the standard deviation of inter-SCaLT intervals against inter-SCaLT interval averages of multiple model and experimental realizations, we revealed the following: (1) SCaLTs exhibit very similar characteristics between the two data sets, (2) they are mostly random, (3) they encode information in their frequency, and (4) their slow baseline oscillations could be due to the stochastic slow clustering of inositol-trisphosphate receptors (modeled as an Ornstein-Uhlenbeck noise process). Bifurcation analysis of a deterministic temporal version of the model showed that the contribution of fluxes to SCaLTs depends on the parameter regime and that the combination of excitability, stochasticity, and mixed-mode oscillations are responsible for irregular spiking and doublets in SCaLTs. Additionally, our results demonstrated that blocking each flux reduces SCaLTs' frequency and that the reverse (forward) mode of Na+/Ca2+ exchange decreases (increases) SCaLTs. Taken together, these results provide a quantitative framework for SCaLT formation in oligodendrocyte precursor cells.

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.

Cells use molecular working memory to navigate inchanging chemoattractant fields

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

Cellular Tango: how extracellular matrix adhesion choreographs Rac-Rho signaling and cell movement

The small GTPases Rac and Rho are known to regulate eukaryotic cell shape, promoting front protrusion (Rac) or rear retraction (Rho) of the cell edge. Such cell deformation changes the contact and adhesion of cell to the extracellular matrix (ECM), while ECM signaling through integrin receptors also affects GTPase activity. We develop and investigate a model for this three-way feedback loop in 1D and 2D spatial domains, as well as in a fully deforming 2D cell shapes with detailed adhesion-bond biophysics. The model consists of reaction-diffusion equations solved numerically with open-source software, Morpheus, and with custom-built cellular Potts model simulations. We find a variety of patterns and cell behaviors, including persistent polarity, flipped front-back cell polarity oscillations, spiral waves, and random protrusion-retraction. We show that the observed spatial patterns depend on the cell shape, and vice versa.

Spots, stripes, and spiral waves in models for static and motile cells : GTPase patterns in cells

The polarization and motility of eukaryotic cells depends on assembly and contraction of the actin cytoskeleton and its regulation by proteins called GTPases. The activity of GTPases causes assembly of filamentous actin (by GTPases Cdc42, Rac), resulting in protrusion of the cell edge. Mathematical models for GTPase dynamics address the spontaneous formation of patterns and nonuniform spatial distributions of such proteins in the cell. Here we revisit the wave-pinning model for GTPase-induced cell polarization, together with a number of extensions proposed in the literature. These include introduction of sources and sinks of active and inactive GTPase (by the group of A. Champneys), and negative feedback from F-actin to GTPase activity. We discuss these extensions singly and in combination, in 1D, and 2D static domains. We then show how the patterns that form (spots, waves, and spirals) interact with cell boundaries to create a variety of interesting and dynamic cell shapes and motion.

Membrane Tension Can Enhance Adaptation to Maintain Polarity of Migrating Cells

Migratory cells are known to adapt to environments that contain wide-ranging levels of chemoattractant. Although biochemical models of adaptation have been previously proposed, here, we discuss a different mechanism based on mechanosensing, in which the interaction between biochemical signaling and cell tension facilitates adaptation. We describe and analyze a model of mechanochemical-based adaptation coupling a mechanics-based physical model of cell tension coupled with the wave-pinning reaction-diffusion model for Rac GTPase activity. The mathematical analysis of this model, simulations of a simplified one-dimensional cell geometry, and two-dimensional finite element simulations of deforming cells reveal that as a cell protrudes under the influence of high stimulation levels, tension-mediated inhibition of Rac signaling causes the cell to polarize even when initially overstimulated. Specifically, tension-mediated inhibition of Rac activation, which has been experimentally observed in recent years, facilitates this adaptation by countering the high levels of environmental stimulation. These results demonstrate how tension-related mechanosensing may provide an alternative (and potentially complementary) mechanism for cell adaptation.

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.

Rac activation is key to cell motility and directionality: An experimental and modelling investigation

Cell migration is a tightly-regulated process that involves protein gradients formed by the Rho family of GTPases, including Rho and Rac. The front (rear) of cells is generally characterized by higher active Rac (Rho) and lower active Rho (Rac) concentrations. Protein clusters, called adhesions, that anchor cells to their external environment have been shown to be dynamic and small (stable and large) at the cell front (rear), forming the force-transmission points necessary for persistent movement. Differences in adhesion sizes and dynamics have been linked to gradients in Rac and Rho activity. Here, we study the effects of Rac activation and gradients in Rac and Rho concentrations and activities on cellular polarity and adhesion size using mathematical and experimental approaches. The former is accomplished by expanding an existing reaction-diffusion model to a 2D domain utilizing stochastic dynamics. The model revealed that a hysteresis between the induced/uninduced states (corresponding to higher/lower Rac concentrations, respectively) along with Rac and Rho activation gradients, generated by chemical cues, were vital for forming polarity. Experimentally, the induced state was generated by increasing the cellular βPIX (a Rac-GEF) level and/or decreasing ROCK (a Rac-GAP effector protein) activity with Y-27632 (a ROCK-inhibitor). In agreement with the simulations, our results showed that cells with elevated RacGTP migrated faster, indicating more robust cellular polarization. However, the directionality of cells was not changed significantly, suggesting that external and/or internal physical or chemical cues were needed. Complementing the faster migration observed, adhesions were smaller, generating the phenotype expected with the induced state.

Small GTPase patterning: How to stabilise cluster coexistence

Many biological processes have to occur at specific locations on the cell membrane. These locations are often specified by the localised activity of small GTPase proteins. Some processes require the formation of a single cluster of active GTPase, also called unipolar polarisation (here “polarisation”), whereas others need multiple coexisting clusters. Moreover, sometimes the pattern of GTPase clusters is dynamically regulated after its formation. This raises the question how the same interacting protein components can produce such a rich variety of naturally occurring patterns. Most currently used models for GTPase-based patterning inherently yield polarisation. Such models may at best yield transient coexistence of at most a few clusters, and hence fail to explain several important biological phenomena. These existing models are all based on mass conservation of total GTPase and some form of direct or indirect positive feedback. Here, we show that either of two biologically plausible modifications can yield stable coexistence: including explicit GTPase turnover, i.e., breaking mass conservation, or negative feedback by activation of an inhibitor like a GAP. Since we start from two different polarising models our findings seem independent of the precise self-activation mechanism. By studying the net GTPase flows among clusters, we provide insight into how these mechanisms operate. Our coexistence models also allow for dynamical regulation of the final pattern, which we illustrate with examples of pollen tube growth and the branching of fungal hyphae. Together, these results provide a better understanding of how cells can tune a single system to generate a wide variety of biologically relevant patterns.

Mathematical Models for Cell Migration with Real-Time Cell Cycle Dynamics

The fluorescent ubiquitination-based cell cycle indicator, also known as FUCCI, allows the visualization of the G1 and S/G2/M cell cycle phases of individual cells. FUCCI consists of two fluorescent probes, so that cells in the G1 phase fluoresce red and cells in the S/G2/M phase fluoresce green. FUCCI reveals real-time information about cell cycle dynamics of individual cells, and can be used to explore how the cell cycle relates to the location of individual cells, local cell density, and different cellular microenvironments. In particular, FUCCI is used in experimental studies examining cell migration, such as malignant invasion and wound healing. Here we present, to our knowledge, new mathematical models that can describe cell migration and cell cycle dynamics as indicated by FUCCI. The fundamental model describes the two cell cycle phases, G1 and S/G2/M, which FUCCI directly labels. The extended model includes a third phase, early S, which FUCCI indirectly labels. We present experimental data from scratch assays using FUCCI-transduced melanoma cells, and show that the predictions of spatial and temporal patterns of cell density in the experiments can be described by the fundamental model. We obtain numerical solutions of both the fundamental and extended models, which can take the form of traveling waves. These solutions are mathematically interesting because they are a combination of moving wavefronts and moving pulses. We derive and confirm a simple analytical expression for the minimum wave speed, as well as exploring how the wave speed depends on the spatial decay rate of the initial condition.

Paxillin phosphorylation at serine 273 and its effects on Rac, Rho and adhesion dynamics

Focal adhesions are protein complexes that anchor cells to the extracellular matrix. During migration, the growth and disassembly of these structures are spatiotemporally regulated, with new adhesions forming at the leading edge of the cell and mature adhesions disassembling at the rear. Signalling proteins and structural cytoskeletal components tightly regulate adhesion dynamics. Paxillin, an adaptor protein within adhesions, is one of these proteins. Its phosphorylation at serine 273 (S273) is crucial for maintaining fast adhesion assembly and disassembly. Paxillin is known to bind to a GIT1-βPIX-PAK1 complex, which increases the local activation of the small GTPase Rac. To understand quantitatively the behaviour of this system and how it relates to adhesion assembly/disassembly, we developed a mathematical model describing the dynamics of the small GTPases Rac and Rho as determined by paxillin S273 phosphorylation. Our model revealed that the system possesses bistability, where switching between uninduced (active Rho) and induced (active Rac) states can occur through a change in rate of paxillin phosphorylation or PAK1 activation. The bistable switch is characterized by the presence of memory, minimal change in the levels of active Rac and Rho within the induced and uninduced states, respectively, and the limited regime of monostability associated with the uninduced state. These results were validated experimentally by showing the presence of bimodality in adhesion assembly and disassembly rates, and demonstrating that Rac activity increases after treating Chinese Hamster Ovary cells with okadaic acid (a paxillin phosphatase inhibitor), followed by a modest recovery after 20 min washout. Spatial gradients of phosphorylated paxillin in a reaction-diffusion model gave rise to distinct regions of Rac and Rho activities, resembling polarization of a cell into front and rear. Perturbing several parameters of the model also revealed important insights into how signalling components upstream and downstream of paxillin phosphorylation affect dynamics.

Cellular migration is crucial in both physiological and pathological functions. Maintenance of proper migration and development of aberrant migration are effectuated by cellular machinery involving protein complexes, called adhesions, that anchor the cell to its environment. Over time, these adhesions assemble at the leading edge, as the cell extends forward, anchoring the front of the cells to its substrate, while those at the cell rear disassemble, allowing detachment and forward movement. Their dynamics are controlled by a number of regulatory factors, occurring on both cell-wide and adhesion-level scales. The coordination of these regulatory factors is complex, but insights about their dynamics can be gained from the use of mathematical modeling techniques which integrate many of these components together. Here, we developed several molecularly explicit models to explore how local regulation of paxillin, an adhesion protein, interacts with the activities of Rac and Rho to produce cell-wide polarization associated with motility and directionality. By altering paxillin phosphorylation/dephosphorylation within such models, we have advanced our understanding of how a shift from a non-motile state to a highly motile state occurs. Deciphering these key processes quantitatively thus helped us gain insight into the subcellular factors underlying polarity and movement.

Principles that govern competition or co-existence in Rho-GTPase driven polarization

Rho-GTPases are master regulators of polarity establishment and cell morphology. Positive feedback enables concentration of Rho-GTPases into clusters at the cell cortex, from where they regulate the cytoskeleton. Different cell types reproducibly generate either one (e.g. the front of a migrating cell) or several clusters (e.g. the multiple dendrites of a neuron), but the mechanistic basis for unipolar or multipolar outcomes is unclear. The design principles of Rho-GTPase circuits are captured by two-component reaction-diffusion models based on conserved aspects of Rho-GTPase biochemistry. Some such models display rapid winner-takes-all competition between clusters, yielding a unipolar outcome. Other models allow prolonged co-existence of clusters. We investigate the behavior of a simple class of models and show that while the timescale of competition varies enormously depending on model parameters, a single factor explains a large majority of this variation. The dominant factor concerns the degree to which the maximal active GTPase concentration in a cluster approaches a "saturation point" determined by model parameters. We suggest that both saturation and the effect of saturation on competition reflect fundamental properties of the Rho-GTPase polarity machinery, regardless of the specific feedback mechanism, which predict whether the system will generate unipolar or multipolar outcomes.

Modelling compartmentalization towards elucidation and engineering of spatial organization in biochemical pathways

Compartmentalization is a fundamental ingredient, central to the functioning of biological systems at multiple levels. At the cellular level, compartmentalization is a key aspect of the functioning of biochemical pathways and an important element used in evolution. It is also being exploited in multiple contexts in synthetic biology. Accurate understanding of the role of compartments and designing compartmentalized systems needs reliable modelling/systems frameworks. We examine a series of building blocks of signalling and metabolic pathways with compartmental organization. We systematically analyze when compartmental ODE models can be used in these contexts, by comparing these models with detailed reaction-transport models, and establishing a correspondence between the two. We build on this to examine additional complexities associated with these pathways, and also examine sample problems in the engineering of these pathways. Our results indicate under which conditions compartmental models can and cannot be used, why this is the case, and what augmentations are needed to make them reliable and predictive. We also uncover other hidden consequences of employing compartmental models in these contexts. Or results contribute a number of insights relevant to the modelling, elucidation, and engineering of biochemical pathways with compartmentalization, at the core of systems and synthetic biology.

A mathematical model coupling polarity signaling to cell adhesion explains diverse cell migration patterns

Protrusion and retraction of lamellipodia are common features of eukaryotic cell motility. As a cell migrates through its extracellular matrix (ECM), lamellipod growth increases cell-ECM contact area and enhances engagement of integrin receptors, locally amplifying ECM input to internal signaling cascades. In contrast, contraction of lamellipodia results in reduced integrin engagement that dampens the level of ECM-induced signaling. These changes in cell shape are both influenced by, and feed back onto ECM signaling. Motivated by experimental observations on melanoma cells lines (1205Lu and SBcl2) migrating on fibronectin (FN) coated topographic substrates (anisotropic post-density arrays), we probe this interplay between intracellular and ECM signaling. Experimentally, cells exhibited one of three lamellipodial dynamics: persistently polarized, random, or oscillatory, with competing lamellipodia oscillating out of phase (Park et al., 2017). Pharmacological treatments, changes in FN density, and substrate topography all affected the fraction of cells exhibiting these behaviours. We use these observations as constraints to test a sequence of hypotheses for how intracellular (GTPase) and ECM signaling jointly regulate lamellipodial dynamics. The models encoding these hypotheses are predicated on mutually antagonistic Rac-Rho signaling, Rac-mediated protrusion (via activation of Arp2/3 actin nucleation) and Rho-mediated contraction (via ROCK phosphorylation of myosin light chain), which are coupled to ECM signaling that is modulated by protrusion/contraction. By testing each model against experimental observations, we identify how the signaling layers interact to generate the diverse range of cell behaviors, and how various molecular perturbations and changes in ECM signaling modulate the fraction of cells exhibiting each. We identify several factors that play distinct but critical roles in generating the observed dynamic: (1) competition between lamellipodia for shared pools of Rac and Rho, (2) activation of RhoA by ECM signaling, and (3) feedback from lamellipodial growth or contraction to cell-ECM contact area and therefore to the ECM signaling level.

Mechanisms of Cell Polarization

Cell polarization is a key step in the migration, development, and organization of eukaryotic cells, both at the single cell and multicellular level. Research on the mechanisms that give rise to polarization of a given cell, and organization of polarity within a tissue has led to new understanding across cellular and developmental biology. In this review, we describe some of the history of theoretical and experimental aspects of the field, as well as some interesting questions and challenges for the future.

Modeling the roles of protein kinase Cβ and η in single-cell wound repair

Wounded cells such as Xenopus oocytes respond to damage by assembly and closure of an array of actin filaments and myosin-2 controlled by Rho GTPases, including Rho and Cdc42. Rho and Cdc42 are patterned around wounds in a characteristic manner, with active Rho concentrating in a ring-like zone inside a larger, ring-like zone of active Cdc42. How this patterning is achieved is unknown, but Rho and Cdc42 at wounds are subject to regulation by other proteins, including the protein kinases C. Specifically, Cdc42 and Rho activity are enhanced by PKCβ and inhibited by PKCη. We adapt a mathematical model of Simon and coworkers to probe the possible roles of these kinases. We show that PKCβ likely affects the magnitude of positive Rho-Abr feedback, whereas PKCη acts on Cdc42 inactivation. The model explains both qualitative and some overall quantitative features of PKC-Rho GTPase regulation. It also accounts for the previous, peculiar observation that ∼ 20% of cells overexpressing PKCη display zone inversions--that is, displacement of active Rho to the outside of the active Cdc42.