Holmes WR , Liao L , Bement W , Edelstein-Keshet L .

	FIGURE 1:. Schematic diagram showing the Rho-Abr-Cdc42 cross-talk as modeled in Simon et al. (2013). Here we consider the possible effects of PKCs (blue arrows) on either the basal activation rates of RhoA and Cdc42 or the magnitude of positive feedback from Cdc42 to itself and from Rho via Abr to itself.
	FIGURE 2:. Model geometry, control simulation, and failure of the simplest model. (a) Geometry of the wound, showing the Rho and Cdc42 zones (green and red, here and throughout all figures) and our one-dimensional model (x is distance from the wound edge). (b) A control simulation showing the distributions of Rho, Cdc42, and Abr for a “wild-type” model cell, with parameter values in Supplemental Table S1. Initial conditions are shown as dashed lines, final profiles 20 s later are shown as solid lines. (c–f) Model results conflict with experimental observations when we assume that PKCs affect only the basal activation or inactivation rates of Rho and Cdc42. (c) PKCβ DN (full suppression of basal activation rates bR, bC) fails to extinguish the Cdc42 zone. (d) PKCβ OE (bR, bC increased by a factor of 1.33) hardly changes the zone intensities, failing to make them significantly brighter. (e–f) PKCη is assumed to influence basal inactivation rates (dR and dC) by the indicated fold change. (e) PKCη DN (no major disagreement). (f) PKCη OE suppresses Rho but allows Cdc42 to invade the region next to the wound edge, contrary to what is observed experimentally. In all cases, simulations are run for 20 s to determine how zones localize before wound closure.
	FIGURE 3:. (a) Reduced Rho–Abr phase plane showing the Abr nullcline (black) and sample Rho nullclines (control as well as DN, OE2, and OE3) for multiple values of γAR (0.5, 1, 2, and threefold multiples of the value stated in Supplemental Table S1; all other parameters are as in that table). Intersections (shown in purple) are steady-state Rho intensity levels in the Rho zone, adjacent to the wound, and the black dot is background Rho activity level. (b) Bifurcation diagram showing steady-state Rho activity levels as a function of the Rho–Abr positive feedback strength γAR. The lower solid line represents the background Rho level; the upper branch represents the Rho zone intensity. Note that the two steady states can coexist over a range of the parameter values and that the intensity of the zone (but not of the background) increases sharply as the positive feedback strength γAR increases. (c) As in a, but for several values of the basal activation rate (bR). (d) As in b, but for the parameter bR. Note that both background and Rho-zone intensities increase with bR but that the zone intensity increases less strongly than in b. (e) Simulation of the model “control” (wild-type cells) using parameters in Supplemental Table S1. (Initial conditions, dashed lines; final profiles, solid lines.) (f) Similar to e, but with both feedback strengths γAR and γC multiplied by a factor of 2 to mimic PKCβ OE (equivalent to 2× in a). Note that zones get significantly “brighter.” (g) Similar to f, but with a threefold increase in the feedback strengths. (h) Similar to e–g, but with a factor of 0.5 reduction in feedback strengths to mimic PKCβ DN. Here all zones are abolished.
	FIGURE 4:. (a) Diagram showing the Cdc42 rates of activation (red curve) and three sample inactivation (black) rates vs. Cdc42 activity levels. The three black lines correspond to PKCη OE, WT, and DN (dηC = −0.15, 0, 0.1, respectively). Intersections of red and black curves are Cdc42 steady-state levels, with red dots representing the elevated Cdc42 intensities in the Cdc42 zone and the blue dots (practically all the same) indicating background Cdc42 activity levels. All other parameters are as in Supplemental Table S1. (b) Cdc42 bifurcation plot with respect to the strength of PKCη-mediated Cdc42 inactivation. (c) Control (WT) simulation for comparison purposes. (d) Spatial behavior. Formation of a single Cdc42 zone for spatially uniform PKCη OE (dηC and dηR = 0.1 with η = 1 uniformly on the domain.) (e) Simulation of reduced inactivation for PKCη DN: the parameters dηC and dηR are reduced, but zone intensity hardly changes.
	FIGURE 5:. Model predictions for graded PKCη activity and predicted zone inversion. (a) Bifurcation plot for the case where PKCη jointly influences Rho and Cdc42. For PKCη activity levels higher than the value indicated by green (respectively red) line, the Rho (respectively Cdc42) zone cannot form due to excessive inactivation (bistability is required for each zone formation and only exists below those PKCη activity levels). (b) The assumed PKCη spatial distribution (linearly graded, with value η = 3 at the wound edge and η = 0 at 20 μm away). Red and green lines have the same meaning as in a. (c) Control simulation for comparison. (d) Graded PKCη overexpression (dηR = dηC = 0.1) suppresses the Rho zone. The Cdc42 zone persists and stays in the same location. (e) Same as in d, but with an extended Rho initial condition that represents some fraction (20%) of PKCη OE cells. The initial condition for Rho is the same as in previous simulations but is extended to a distance of 20 μm from the wound edge (mimicking a possible extension of the zone of initial Rho activity). Zone inversion occurs.

Abreu-Blanco, Coordination of Rho family GTPase activities to orchestrate cytoskeleton responses during cell wound repair. 2014, Pubmed

Abreu-Blanco, Coordination of Rho family GTPase activities to orchestrate cytoskeleton responses during cell wound repair. 2014, Pubmed
Amano, Regulation and functions of Rho-associated kinase. 2000, Pubmed
Bement, Rho GTPase activity zones and transient contractile arrays. 2006, Pubmed
Benink, Concentric zones of active RhoA and Cdc42 around single cell wounds. 2005, Pubmed , Xenbase
Bishop, Rho GTPases and their effector proteins. 2000, Pubmed
Burkel, A Rho GTPase signal treadmill backs a contractile array. 2012, Pubmed , Xenbase
Chuang, Abr and Bcr are multifunctional regulators of the Rho GTP-binding protein family. 1995, Pubmed
Clark, Integration of single and multicellular wound responses. 2009, Pubmed , Xenbase
Dovas, Serine 34 phosphorylation of rho guanine dissociation inhibitor (RhoGDIalpha) links signaling from conventional protein kinase C to RhoGTPase in cell adhesion. 2010, Pubmed
Dovas, RhoGDI: multiple functions in the regulation of Rho family GTPase activities. 2005, Pubmed
Edelstein-Keshet, From simple to detailed models for cell polarization. 2013, Pubmed
Ferrell, Bistability in cell signaling: How to make continuous processes discontinuous, and reversible processes irreversible. 2001, Pubmed , Xenbase
Holmes, Modelling cell polarization driven by synthetic spatially graded Rac activation. 2012, Pubmed
Holmes, An efficient, nonlinear stability analysis for detecting pattern formation in reaction diffusion systems. 2014, Pubmed
Holmes, Local perturbation analysis: a computational tool for biophysical reaction-diffusion models. 2015, Pubmed
Jilkine, Mathematical model for spatial segregation of the Rho-family GTPases based on inhibitory crosstalk. 2007, Pubmed
Kono, Proteasomal degradation resolves competition between cell polarization and cellular wound healing. 2012, Pubmed
Lin, Synthetic spatially graded Rac activation drives cell polarization and movement. 2012, Pubmed
Mandato, Contraction and polymerization cooperate to assemble and close actomyosin rings around Xenopus oocyte wounds. 2001, Pubmed , Xenbase
Marée, Polarization and movement of keratocytes: a multiscale modelling approach. 2006, Pubmed
Moon, Rho GTPase-activating proteins in cell regulation. 2003, Pubmed
Mori, ASYMPTOTIC AND BIFURCATION ANALYSIS OF WAVE-PINNING IN A REACTION-DIFFUSION MODEL FOR CELL POLARIZATION. 2011, Pubmed
Nobes, Rho, rac and cdc42 GTPases: regulators of actin structures, cell adhesion and motility. 1995, Pubmed
Olofsson, Rho guanine dissociation inhibitors: pivotal molecules in cellular signalling. 1999, Pubmed
Rossman, GEF means go: turning on RHO GTPases with guanine nucleotide-exchange factors. 2005, Pubmed
Simon, Pattern formation of Rho GTPases in single cell wound healing. 2013, Pubmed , Xenbase
Sonnemann, Wound repair: toward understanding and integration of single-cell and multicellular wound responses. 2011, Pubmed
Tyson, Sniffers, buzzers, toggles and blinkers: dynamics of regulatory and signaling pathways in the cell. 2003, Pubmed
Vaughan, Control of local Rho GTPase crosstalk by Abr. 2011, Pubmed , Xenbase
Vaughan, Lipid domain-dependent regulation of single-cell wound repair. 2014, Pubmed , Xenbase