Quantitative model of Ras-phosphoinositide 3-kinase signalling cross-talk based on co-operative molecular assembly

Insight from the maximal activation of the signal transduction excitable network in Dictyostelium discoideum

Cell migration requires the coordination of an excitable signal transduction network involving Ras and PI3K pathways with cytoskeletal activity. We show that expressing activated Ras GTPase-family proteins in cells lacking PTEN or other mutations which increase cellular protrusiveness transforms cells into a persistently activated state. Leading- and trailing-edge markers were found exclusively at the cell perimeter and the cytosol, respectively, of the dramatically flattened cells. In addition, the lifetimes of dynamic actin puncta were increased where they overlapped with actin waves, suggesting a mechanism for the coupling between these two networks. All of these phenotypes could be reversed by inhibiting signal transduction. Strikingly, maintaining cells in this state of constant activation led to a form of cell death by catastrophic fragmentation. These findings provide insight into the feedback loops that control excitability of the signal transduction network, which drives migration.

Excitable Signal Transduction Networks in Directed Cell Migration

Although directed migration of eukaryotic cells may have evolved to escape nutrient depletion, it has been adopted for an extensive range of physiological events during development and in the adult organism. The subversion of these movements results in disease, such as cancer. Mechanisms of propulsion and sensing are extremely diverse, but most eukaryotic cells move by extending actin-filled protrusions termed macropinosomes, pseudopodia, or lamellipodia or by extension of blebs. In addition to motility, directed migration involves polarity and directional sensing. The hundreds of gene products involved in these processes are organized into networks of parallel and interconnected pathways. Many of these components are activated or inhibited coordinately with stimulation and on each spontaneously extended protrusion. Moreover, these networks display hallmarks of excitability, including all-or-nothing responsiveness and wave propagation. Cellular protrusions result from signal transduction waves that propagate outwardly from an origin and drive cytoskeletal activity. The range of the propagating waves and hence the size of the protrusions can be altered by lowering or raising the threshold for network activation, with larger and wider protrusions favoring gliding or oscillatory behavior over amoeboid migration. Here, we evaluate the variety of models of excitable networks controlling directed migration and outline critical tests. We also discuss the utility of this emerging view in producing cell migration and in integrating the various extrinsic cues that direct migration.

Exosomal miRs in Lung Cancer: A Mathematical Model

Lung cancer, primarily non-small-cell lung cancer (NSCLC), is the leading cause of cancer deaths in the United States and worldwide. While early detection significantly improves five-year survival, there are no reliable diagnostic tools for early detection. Several exosomal microRNAs (miRs) are overexpressed in NSCLC, and have been suggested as potential biomarkers for early detection. The present paper develops a mathematical model for early stage of NSCLC with emphasis on the role of the three highest overexpressed miRs, namely miR-21, miR-205 and miR-155. Simulations of the model provide quantitative relationships between the tumor volume and the total mass of each of the above miRs in the tumor. Because of the positive correlation between these miRs in the tumor tissue and in the blood, the results of the paper may be viewed as a first step toward establishing a combination of miRs 21, 205, 155 and possibly other miRs as serum biomarkers for early detection of NSCLC.

Dimerization-based control of cooperativity

Cooperativity of ligand-receptor binding influences the input-output behavior of a biochemical system and thus is an important determinant of its physiological function. Canonically, such cooperativity is understood in terms of ligand-receptor binding affinity, where an initial binding event changes the affinity for subsequent binding events. Here, we demonstrate that dimerization-a simple yet pervasive signaling motif across biology-can have significant control over cooperativity and even dominate over the canonical mechanism. Through an exhaustive parameter sensitivity analysis of a general kinetic model for signal-mediated dimerization, we show that quantitative modulation of dimerization processes can reinforce, eliminate, and even reverse cooperativity imposed by the canonical allosteric ligand-receptor binding affinity mechanism. The favored accumulation of stoichiometrically asymmetric dimers (those with ligand-receptor stoichiometry of 1 : 2) is a major determinant of dimerization-based cooperativity control. However, simulations demonstrate that favoring accumulation of such stoichiometrically asymmetric dimers can either increase or decrease cooperativity, and thus the quantitative relationship between stoichiometrically asymmetric dimers and cooperativity is highly dependent on the parameter values of the particular system of interest. These results suggest that the dimerization motif provides a novel mechanism for both generating and quantitatively tuning cooperativity that, due to the ubiquity of dimerization motifs in biochemical systems, may play a major role in a host of biological functions. Thus, the canonical, allosteric view of cooperativity is incomplete without considering dimerization effects, which is of particular importance as dimerization is often a necessary feature of the allosteric mechanism.

Prognostic biomarkers in phase II trial of cetuximab-containing induction and chemoradiation in resectable HNSCC: Eastern cooperative oncology group E2303

PURPOSE

We sought to evaluate the correlation between tissue biomarker expression (using standardized, quantitative immunofluorescence) and clinical outcome in the E2303 trial.

EXPERIMENTAL DESIGN

Sixty-three eligible patients with operable stage III/IV head and neck squamous cell cancer (HNSCC) participated in the Eastern Cooperative Oncology Group (ECOG) 2303 phase II trial of induction chemotherapy with weekly cetuximab, paclitaxel, and carboplatin followed by chemoradiation with the same regimen. A tissue microarray (TMA) was constructed and EGF receptor (EGFR), ERK1/2, Met, Akt, STAT3, β-catenin, E-cadherin, EGFR Variant III, insulin-like growth factor-1 receptor, NF-κB, p53, PI3Kp85, PI3Kp110a, PTEN, NRAS, and pRb protein expression levels were assessed using automated quantitative protein analysis (AQUA). For each dichotomized biomarker, overall survival (OS), progression-free survival (PFS), and event-free survival (EFS) were estimated by the Kaplan-Meier method and compared using log-rank tests. Multivariable Cox proportional hazards models were used to estimate HRs and test for significance.

RESULTS

Forty-two of 63 patients with TMA data on at least one biomarker were included in the biomarker analysis. Tumor extracellular signal-regulated kinase (ERK)1/2 levels were significantly associated with PFS [HR (low/high), 3.29; P = 0.026] and OS [HR (low/high), 4.34; P = 0.008]. On multivariable Cox regression analysis, ERK1/2 remained significantly associated with OS (P = 0.024) and PFS (P = 0.022) after controlling for primary site (oropharynx vs. non-oropharynx) and disease stage (III vs. IV), respectively. Clustering analysis revealed that clusters indicative of activated RAS/MAPK/ERK and/or PI3K/Akt pathways were associated with inferior OS and/or PFS and maintained significance in multivariable analysis.

CONCLUSIONS

These results implicate PI3K/Akt and RAS/MAPK/ERK pathways in resistance to cetuximab-containing chemoradiation in HNSCC. Large prospective studies are required to validate these results.

Computational model of VEGFR2 pathway to ERK activation and modulation through receptor trafficking

Vascular Endothelial Growth Factor (VEGF) signal transduction is central to angiogenesis in development and in pathological conditions such as cancer, retinopathy and ischemic diseases. We constructed and validated a computational model of VEGFR2 trafficking and signaling, to study the role of receptor trafficking kinetics in modulating ERK phosphorylation in VEGF-stimulated endothelial cells. Trafficking parameters were optimized and validated against four previously published in vitro experiments. Based on these parameters, model simulations demonstrated interesting behaviors that may be highly relevant to understanding VEGF signaling in endothelial cells. First, at moderate VEGF doses, VEGFR2 phosphorylation and ERK phosphorylation are related in a log-linear fashion, with a stable duration of ERK activation; but with higher VEGF stimulation, phosphoERK becomes saturated, and its duration increases. Second, a large endosomal fraction of VEGFR2 makes the ERK activation reaction network less sensitive to perturbations in VEGF dosage. Third, extracellular-matrix-bound VEGF binds and activates VEGFR2, but by internalizing at a slower rate, matrix-bound VEGF-induced intracellular ERK phosphorylation is predicted to be greater in magnitude and more sustained, in agreement with experimental evidence. Fourth, different endothelial cell types appear to have different trafficking rates, which result in different levels of endosomal receptor localization and different ERK response profiles.

Computational Model of Gab1/2-Dependent VEGFR2 Pathway to Akt Activation

Vascular endothelial growth factor (VEGF) signal transduction is central to angiogenesis in development and in pathological conditions such as cancer, retinopathy and ischemic diseases. However, no detailed mass-action models of VEGF receptor signaling have been developed. We constructed and validated the first computational model of VEGFR2 trafficking and signaling, to study the opposing roles of Gab1 and Gab2 in regulation of Akt phosphorylation in VEGF-stimulated endothelial cells. Trafficking parameters were optimized against 5 previously published in vitro experiments, and the model was validated against six independent published datasets. The model showed agreement at several key nodes, involving scaffolding proteins Gab1, Gab2 and their complexes with Shp2. VEGFR2 recruitment of Gab1 is greater in magnitude, slower, and more sustained than that of Gab2. As Gab2 binds VEGFR2 complexes more transiently than Gab1, VEGFR2 complexes can recycle and continue to participate in other signaling pathways. Correspondingly, the simulation results show a log-linear relationship between a decrease in Akt phosphorylation and Gab1 knockdown while a linear relationship was observed between an increase in Akt phosphorylation and Gab2 knockdown. Global sensitivity analysis demonstrated the importance of initial-concentration ratios of antagonistic molecular species (Gab1/Gab2 and PI3K/Shp2) in determining Akt phosphorylation profiles. It also showed that kinetic parameters responsible for transient Gab2 binding affect the system at specific nodes. This model can be expanded to study multiple signaling contexts and receptor crosstalk and can form a basis for investigation of therapeutic approaches, such as tyrosine kinase inhibitors (TKIs), overexpression of key signaling proteins or knockdown experiments.

Coordinate direct input of both KRAS and IGF1 receptor to activation of PI3 kinase in KRAS-mutant lung cancer

Using a panel of non-small cell lung cancer (NSCLC) lines, we show here that MAP-ERK kinase (MEK) and RAF inhibitors are selectively toxic for the KRAS-mutant genotype, whereas phosphoinositide 3-kinase (PI3K), AKT, and mTOR inhibitors are not. IGF1 receptor (IGF1R) tyrosine kinase inhibitors also show selectivity for KRAS-mutant lung cancer lines. Combinations of IGF1R and MEK inhibitors resulted in strengthened inhibition of KRAS-mutant lines and also showed improved effectiveness in autochthonous mouse models of Kras-induced NSCLC. PI3K pathway activity is dependent on basal IGF1R activity in KRAS-mutant, but not wild-type, lung cancer cell lines. KRAS is needed for both MEK and PI3K pathway activity in KRAS-mutant, but not wild-type, lung cancer cells, whereas acute activation of KRAS causes stimulation of PI3K dependent upon IGF1R kinase activity. Coordinate direct input of both KRAS and IGF1R is thus required to activate PI3K in KRAS-mutant lung cancer cells.

SIGNIFICANCE

It has not yet been possible to target RAS proteins directly, so combined targeting of effect or pathways acting downstream of RAS, including RAF/MEK and PI3K/AKT, has been the most favored approach to the treatment of RAS -mutant cancers. This work sheds light on the ability of RASto activate PI3K through direct interaction, indicating that input is also required from a receptor tyrosinekinase, IGF1R in the case of KRAS -mutant lung cancer. This suggests potential novel combination therapeutic strategies for NSCLC.

A mathematical model of CR3/TLR2 crosstalk in the context of Francisella tularensis infection

Complement Receptor 3 (CR3) and Toll-like Receptor 2 (TLR2) are pattern recognition receptors expressed on the surface of human macrophages. Although these receptors are essential components for recognition by the innate immune system, pathogen coordinated crosstalk between them can suppress the production of protective cytokines and promote infection. Recognition of the virulent Schu S4 strain of the intracellular pathogen Francisella tularensis by host macrophages involves CR3/TLR2 crosstalk. Although experimental data provide evidence that Lyn kinase and PI3K are essential components of the CR3 pathway that influences TLR2 activity, additional responsible upstream signaling components remain unknown. In this paper we construct a mathematical model of CR3 and TLR2 signaling in response to F. tularensis. After demonstrating that the model is consistent with experimental results we perform numerical simulations to evaluate the contributions that Akt and Ras-GAP make to ERK inhibition. The model confirms that phagocytosis-associated changes in the composition of the cell membrane can inhibit ERK activity and predicts that Akt and Ras-GAP synergize to inhibit ERK.

Data-driven modelling of receptor tyrosine kinase signalling networks quantifies receptor-specific potencies of PI3K- and Ras-dependent ERK activation

Signal transduction networks in mammalian cells, comprising a limited set of interacting biochemical pathways, are accessed by various growth factor and cytokine receptors to elicit distinct cell responses. This raises the question as to how specificity of the stimulus-response relationship is encoded at the molecular level. It has been proposed that specificity arises not only from the activation of unique signalling pathways, but also from quantitative differences in the activation and regulation of shared receptor-proximal signalling proteins. To address such hypotheses, data sets with greater precision and coverage of experimental conditions will need to be acquired, and rigorous frameworks that codify and parameterize the inherently non-linear relationships among signalling activities will need to be developed. In the present study we apply a systematic approach combining quantitative measurements and mathematical modelling to compare the signalling networks accessed by FGF (fibroblast growth factor) and PDGF (platelet-derived growth factor) receptors in mouse fibroblasts, in which the ERK (extracellular-signal-regulated kinase) cascade is activated by Ras- and PI3K (phosphoinositide 3-kinase)-dependent pathways. We show that, whereas the FGF stimulation of PI3K signalling is relatively weak, this deficiency is compensated for by a more potent Ras-dependent activation of ERK. Thus, as the modelling would predict, the ERK pathway is activated to a greater extent in cells co-stimulated with FGF and PDGF, relative to the saturated levels achieved with either ligand alone. It is envisaged that similar approaches will prove valuable in the elucidation of quantitative differences among other closely related receptor signalling networks.

Allosteric modulation of Ras-GTP is linked to signal transduction through RAF kinase

Ras is a key signal transduction protein in the cell. Mutants of Gly(12) and Gln(61) impair GTPase activity and are found prominently in cancers. In wild type Ras-GTP, an allosteric switch promotes disorder to order transition in switch II, placing Gln(61) in the active site. We show that the "on" and "off" conformations of the allosteric switch can also be attained in RasG12V and RasQ61L. Although both mutants have similarly impaired active sites in the on state, RasQ61L stabilizes an anti-catalytic conformation of switch II in the off state of the allosteric switch when bound to Raf. This translates into more potent activation of the MAPK pathway involving Ras, Raf kinase, MEK, and ERK (Ras/Raf/MEK/ERK) in cells transfected with RasQ61L relative to RasG12V. This differential is not observed in the Raf-independent pathway involving Ras, phosphoinositide 3-kinase (PI3K), and Akt (Ras/PI3K/Akt). Using a combination of structural analysis, hydrolysis rates, and experiments in NIH-3T3 cells, we link the allosteric switch to the control of signaling in the Ras/Raf/MEK/ERK pathway, supporting a GTPase-activating protein-independent model for duration of the Ras-Raf complex.

Systematic quantification of negative feedback mechanisms in the extracellular signal-regulated kinase (ERK) signaling network

Cell responses are actuated by tightly controlled signal transduction pathways. Although the concept of an integrated signaling network replete with interpathway cross-talk and feedback regulation is broadly appreciated, kinetic data of the type needed to characterize such interactions in conjunction with mathematical models are lacking. In mammalian cells, the Ras/ERK pathway controls cell proliferation and other responses stimulated by growth factors, and several cross-talk and feedback mechanisms affecting its activation have been identified. In this work, we take a systematic approach to parse the magnitudes of multiple regulatory mechanisms that attenuate ERK activation through canonical (Ras-dependent) and non-canonical (PI3K-dependent) pathways. In addition to regulation of receptor and ligand levels, we consider three layers of ERK-dependent feedback: desensitization of Ras activation, negative regulation of MEK kinase (e.g. Raf) activities, and up-regulation of dual-specificity ERK phosphatases. Our results establish the second of these as the dominant mode of ERK self-regulation in mouse fibroblasts. We further demonstrate that kinetic models of signaling networks, trained on a sufficient diversity of quantitative data, can be reasonably comprehensive, accurate, and predictive in the dynamical sense.

PI3K-dependent cross-talk interactions converge with Ras as quantifiable inputs integrated by Erk

Although it is appreciated that canonical signal-transduction pathways represent dominant modes of regulation embedded in larger interaction networks, relatively little has been done to quantify pathway cross-talk in such networks. Through quantitative measurements that systematically canvas an array of stimulation and molecular perturbation conditions, together with computational modeling and analysis, we have elucidated cross-talk mechanisms in the platelet-derived growth factor (PDGF) receptor signaling network, in which phosphoinositide 3-kinase (PI3K) and Ras/extracellular signal-regulated kinase (Erk) pathways are prominently activated. We show that, while PI3K signaling is insulated from cross-talk, PI3K enhances Erk activation at points both upstream and downstream of Ras. The magnitudes of these effects depend strongly on the stimulation conditions, subject to saturation effects in the respective pathways and negative feedback loops. Motivated by those dynamics, a kinetic model of the network was formulated and used to precisely quantify the relative contributions of PI3K-dependent and -independent modes of Ras/Erk activation.

Spontaneous phosphoinositide 3-kinase signaling dynamics drive spreading and random migration of fibroblasts

During directed cell migration (chemotaxis), cytoskeletal dynamics are stimulated and spatially biased by phosphoinositide 3-kinase (PI3K) and other signal transduction pathways. Live-cell imaging using total internal reflection fluorescence (TIRF) microscopy revealed that, in the absence of soluble cues, 3'-phosphoinositides are enriched in a localized and dynamic fashion during active spreading and random migration of mouse fibroblasts on adhesive surfaces. Surprisingly, we found that PI3K activation is uncoupled from classical integrin-mediated pathways and feedback from the actin cytoskeleton. Inhibiting PI3K significantly impairs cell motility, both in the context of normal spreading and when microtubules are dissociated, which induces a dynamic protrusion phenotype as seen by TIRF in our cells. Accordingly, during random migration, 3'-phosphoinositides are frequently localized to regions of membrane protrusion and correlate quantitatively with the direction and persistence of cell movement. These results underscore the importance of localized PI3K signaling not only in chemotaxis but also in basal motility/migration of fibroblasts.

Signal transduction at point-blank range: analysis of a spatial coupling mechanism for pathway crosstalk

The plasma membrane provides a physical platform for the orchestration of molecular interactions and biochemical conversions involved in the early stages of receptor-mediated signal transduction in living cells. In that context, we introduce here the concept of spatial coupling, wherein simultaneous recruitment of different enzymes to the same receptor scaffold facilitates crosstalk between different signaling pathways through the local release and capture of activated signaling molecules. To study the spatiotemporal dynamics of this mechanism, we have developed a Brownian dynamics modeling approach and applied it to the receptor-mediated activation of Ras and the cooperative recruitment of phosphoinositide 3-kinase (PI3K) by activated receptors and Ras. Various analyses of the model simulations show that cooperative assembly of multimolecular complexes nucleated by activated receptors is facilitated by the local release and capture of membrane-anchored signaling molecules (such as active Ras) from/by receptor-bound signaling proteins. In the case of Ras/PI3K crosstalk, the model predicts that PI3K is more likely to be recruited by activated receptors bound or recently visited by the enzyme that activates Ras. By this mechanism, receptor-bound PI3K is stabilized through short-range, diffusion-controlled capture of active Ras and Ras/PI3K complexes released from the receptor complex. We contend that this mechanism is a means by which signaling pathways are propagated and spatially coordinated for efficient crosstalk between them.

Sensing and integration of Erk and PI3K signals by Myc

The transcription factor Myc plays a central role in regulating cell-fate decisions, including proliferation, growth, and apoptosis. To maintain a normal cell physiology, it is critical that the control of Myc dynamics is precisely orchestrated. Recent studies suggest that such control of Myc can be achieved at the post-translational level via protein stability modulation. Myc is regulated by two Ras effector pathways: the extracellular signal-regulated kinase (Erk) and phosphatidylinositol 3-kinase (PI3K) pathways. To gain quantitative insight into Myc dynamics, we have developed a mathematical model to analyze post-translational regulation of Myc via sequential phosphorylation by Erk and PI3K. Our results suggest that Myc integrates Erk and PI3K signals to result in various cellular responses by differential stability control of Myc protein isoforms. Such signal integration confers a flexible dynamic range for the system output, governed by stability change. In addition, signal integration may require saturation of the input signals, leading to sensitive signal integration to the temporal features of the input signals, insensitive response to their amplitudes, and resistance to input fluctuations. We further propose that these characteristics of the protein stability control module in Myc may be commonly utilized in various cell types and classes of proteins.

Computational models of tandem SRC homology 2 domain interactions and application to phosphoinositide 3-kinase

Intracellular signal transduction proteins typically utilize multiple interaction domains for proper targeting, and thus a broad diversity of distinct signaling complexes may be assembled. Considering the coordination of only two such domains, as in tandem Src homology 2 (SH2) domain constructs, gives rise to a kinetic scheme that is not adequately described by simple models used routinely to interpret in vitro binding measurements. To analyze the interactions between tandem SH2 domains and bisphosphorylated peptides, we formulated detailed kinetic models and applied them to the phosphoinositide 3-kinase p85 regulatory subunit/platelet-derived growth factor beta-receptor system. Data for this system from different in vitro assay platforms, including surface plasmon resonance, competition binding, and isothermal titration calorimetry, were reconciled to estimate the magnitude of the cooperativity characterizing the sequential binding of the high and low affinity SH2 domains (C-SH2 and N-SH2, respectively). Compared with values based on an effective volume approximation, the estimated cooperativity is 3 orders of magnitude lower, indicative of significant structural constraints. Homodimerization of full-length p85 was found to be an alternative mechanism for high avidity binding to phosphorylated platelet-derived growth factor receptors, which would render the N-SH2 domain dispensable for receptor binding.

Antagonism of platelet-derived growth factor receptor in non small cell lung cancer: rationale and investigations

Molecules that target growth and survival pathways in cancer cells have revolutionized the treatment of cancer. Imatinib mesylate is one such agent inhibiting the tyrosine kinase that results from the Bcr-Abl translocation. Imatinib is also a potent inhibitor of the platelet-derived growth factor receptor. The platelet-derived growth factor receptor is crucial in the regulation of interstitial fluid pressure, as well as in the function of pericytes. Increased interstitial fluid pressure is a common feature of solid tumors and is thought to impede transcapillary transport of chemotherapy. Preclinical data show that platelet-derived growth factor receptor antagonism decreases interstitial fluid pressure, augments intratumoral concentration of chemotherapy, and impairs tumor growth. Pericytes are important cells in the vascular support structure of tumors regulating endothelial cell survival and directing capillary growth. Preclinical data suggest that dual targeting of pericytes and endothelial cells is a more effective antiangiogenic strategy than antiendothelial monotherapy. Two phase II studies in advanced non-small cell lung cancer are currently under way with imatinib. The first trial evaluates the use of intermittent imatinib and weekly paclitaxel in elderly patients. The second trial evaluates a novel maintenance strategy of imatinib and the antivascular endothelial growth factor antibody bevacizumab after first-line chemotherapy with bevacizumab. These trials should indicate whether encouraging preclinical data can be translated into clinical benefit in non-small cell lung cancer.

A mathematical model for neutrophil gradient sensing and polarization

Directed cell migration in response to chemical cues, also known as chemotaxis, is an important physiological process involved in wound healing, foraging, and the immune response. Cell migration requires the simultaneous formation of actin polymers at the leading edge and actomyosin complexes at the sides and back of the cell. An unresolved question in eukaryotic chemotaxis is how the same chemoattractant signal determines both the cell's front and back. Recent experimental studies have begun to reveal the biochemical mechanisms necessary for this polarized cellular response. We propose a mathematical model of neutrophil gradient sensing and polarization based on experimentally characterized biochemical mechanisms. The model demonstrates that the known dynamics for Rho GTPase and phosphatidylinositol-3-kinase (PI3K) activation are sufficient for both gradient sensing and polarization. In particular, the model demonstrates that these mechanisms can correctly localize the “front” and “rear” pathways in response to both uniform concentrations and gradients of chemical attractants, including in actin-inhibited cells. Furthermore, the model predictions are robust to the values of many parameters. A key result of the model is the proposed coincidence circuit involving PI3K and Ras that obviates the need for the “global inhibitors” proposed, though never experimentally verified, in many previous mathematical models of eukaryotic chemotaxis. Finally, experiments are proposed to (in)validate this model and further our understanding of neutrophil chemotaxis.

Neutrophils target sites of infection and inflammation by sensing chemical signals produced by damaged tissue and infecting microbes and then move in the direction where their concentration is greatest. An open question is how neutrophils integrate this information to determine the direction of motility. We present a mathematical model for the intracellular signaling network regulating polarization and chemotaxis in neutrophils. We demonstrate how the activation of two antagonizing pathways robustly establishes the front and back of the migrating cell. The model is able to reproduce a number of experimental studies, and new experiments are proposed to test different aspects of the model. A key result is the characterization of a coincidence circuit involving phosphatidylinositol-3-kinase (PI3K) and Ras. We demonstrate that this circuit plays a critical role in selectively localizing F-actin to the front of the cell and actomyosin complexes to the rear. As directed motility in response to chemical cues is critical in a number of processes including wound healing and tumor metastasis, the results and insights gained from the model may be applicable to other cell types and organisms.