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.
Historically, intracellular signal transduction has been characterized in terms of distinct pathways, comprised of linear, sequential activation processes. This concept is exemplified by the canonical mitogen‐activated protein kinase (MAPK) cascades, such as the Ras→Raf→MEK→extracellular signal‐regulated kinase (Erk) pathway in mammals. Our current understanding of signal‐transduction networks includes more complex interactions, including those between the classically defined pathways (cross‐talk) and those responsible for feedback regulation/reinforcement. Although there is a great deal of qualitative information available as to how signaling networks are wired, an as yet unmet challenge is their systematic quantification; to understand cell regulation at the molecular level, we need to know the relative magnitudes of classical and cross‐talk interactions that converge and collaborate to activate key nodes in the network. In most signaling networks, Erk isoforms are both master integrators of upstream inputs and master controllers of transcription factors and other effectors (Kolch, 2000).
In this study, we systematically combine quantitative biochemical measurements and computational modeling to quantify the magnitudes of cross‐talk and negative feedback interactions in a signaling network. We investigated signaling mediated by platelet‐derived growth factor (PDGF) receptors in fibroblasts, leading to activation of the Ras/Erk pathway and exceptionally robust activation of phosphoinositide 3‐kinases (PI3Ks), which control responses such as directed cell migration, survival, and proliferation through the production of specific lipid second messengers (Vanhaesebroeck et al, 2001; Engelman et al, 2006; Hawkins et al, 2006). Although cross‐talk between the PI3K and Ras/Erk pathways has been studied extensively, a definitive pattern of regulation has not yet emerged, because the two pathways apparently affect each other in various ways and in a context‐dependent manner.
We measured Erk phosphorylation and PI3K‐dependent Akt phosphorylation in PDGF‐stimulated NIH 3T3 fibroblasts by quantitative immunoblotting for an array of 126 experimental conditions, sampling different combinations of ligand dose, stimulation time, and molecular manipulation; the Erk data set is shown in Figure 1. Considering biological replicates and parallel determination of total Erk and Akt levels, this data set comprises 2772 total measurements. As shown in Figure 1B and C, blocking the activity of either Ras or PI3K, by the expression of dominant‐negative (S17N) H‐Ras or incubation with a pharmacological inhibitor, respectively, only partially reduces PDGF‐stimulated Erk phosphorylation; the degree of inhibition depends, in a complex manner, on both stimulation dose and time. By comparison, the Akt phosphorylation results showed that the PI3K pathway is not significantly affected by perturbations affecting Ras and Erk; cross‐talk is apparently unidirectional, from PI3K to Ras/Erk, in this network.
We went on to show that simultaneous inhibition of Ras and PI3K almost completely abolished PDGF‐stimulated Erk phosphorylation, indicating that Ras and PI3K are responsible for all of the major pathways from PDGF receptors to Erk, and at least one mode of PI3K‐dependent cross‐talk to Erk is independent of Ras. We further refined this conceptual model by measuring Ras‐GTP, the active form of Ras, for selected experimental conditions using a quantitative enzymatic assay (Scheele et al, 1995; Kaur et al, 2006). These experiments showed that PI3K‐dependent cross‐talk affects the Erk pathway both downstream and upstream of Ras, and they also served to characterize a known negative feedback mechanism affecting Ras as a consequence of signaling through the Erk cascade. A second regulatory loop that we characterized for this pathway is the increased expression of MAPK phosphatase‐1 (MKP‐1), which responds to and antagonizes Erk activation.
Motivated by the dynamics revealed in our unique data set, 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. In the context of the model, the magnitudes of the Ras‐ and PI3K‐dependent inputs converging on MEK/Erk determines the saturability of Erk phosphorylation with respect to PDGF dose and also the degree to which the response adapts, and conversely, the observed dynamics can be used to quantify the magnitudes of the inputs. Unspecified parameter values were estimated using a Monte Carlo‐based algorithm that directly and globally compares the model output to the experimental data; thus, an ensemble of 10 000 parameter sets that fit the data almost equally well was assembled (Brown and Sethna, 2003; Gutenkunst et al, 2007) and compared with the corresponding data as shown in Figure 7. Taken together, the data force the model to reconcile time‐ and PDGF dose‐dependent features of the network observed under the various experimental conditions tested.
Analysis of the parameter sets chosen by the algorithm revealed a consistent ratio of PI3K‐ and Ras‐dependent contributions to the dual phosphorylation of MEK, the kinase activity directly upstream of Erk. We formulated a single number, the MEK activation comparator (MAC), which compares the capacities of the two pathways to generate dually phosphorylated MEK. Importantly, the MAC quantifies these inputs in a way that uncouples them from negative feedback effects. This analysis revealed that, although the PI3K‐dependent MEK activation pathway is predicted to be intrinsically much less potent than the Ras‐dependent pathway under maximal PDGF stimulation conditions, feedback regulation of Ras renders the PI3K‐dependent pathway somewhat more important.
We also used our computational approach to generate hypothetical predictions with an eye toward future experiments. Although inhibition of PI3K affects cross‐talk interactions both upstream and downstream of Ras, the model ensemble predicts unique kinetic signatures that might be expected if either mechanism were silenced selectively, which could help validate the point of action of a particular pathway from PI3K to Erk. In turn, characterization of cross‐talk mechanisms along these lines would allow both conceptual and quantitative models to be refined.
The integration of high‐quality, high‐resolution measurements, directly compared with data‐driven kinetic models, is demonstrated as a means to quantify cross‐talk interactions in signaling networks. For the field to move forward, an issue that will need to be confronted is the generality of signaling networks across receptor and cell types. Certainly, major disparities in the regulation of signaling in different contexts, even among receptors of the same class, should be expected, and more quantitative approaches will be needed to unravel the mechanisms that distinguish particular signaling systems.
To meet the challenge of quantifying crosstalk interactions in signal transduction networks, we have assembled a unique data set comprised of biochemical measurements that systematically canvas an array of stimulation and molecular perturbation conditions.
The experimental results show that, in the platelet‐derived growth factor (PDGF) receptor signaling network, activation of the mitogen‐activated protein kinases (MAPKs) Erk‐1/2 is mediated by two major pathways: one that is strictly Ras‐dependent, and another that is strictly dependent on crosstalk from phosphoinositide 3‐kinase (PI3K); PI3K‐dependent crosstalk and negative feedback affecting Ras, and a second negative feedback loop involving upregulation of MAPK phosphatase activity, were also characterized.
Further quantification of the constituent processes in the network was achieved through data‐driven kinetic modeling and a Monte Carlo, ensemble‐based parameter estimation scheme, and metrics were devised that compare the relative magnitudes of PI3K‐dependent and ‐independent signaling modes.
Another demonstrated use of the model is the generation of ensemble predictions, which might help validate the point(s) of action of a specific pathway mediating crosstalk in a signaling network.
Signal transduction is traditionally characterized in terms of intracellular pathways, which govern outcomes such as cell proliferation, survival, migration, and differentiation. Cross‐talk interactions, which couple distinct pathways, are recognized for their importance in cell regulation yet remain poorly defined because of their complex nature (Hunter, 2000). They allow information flow and regulation to be distributed across multiple pathways, and perturbations targeting specific signaling molecules might therefore have unexpected effects (Bray, 1990). Here, we describe the systematic and quantitative characterization of cross‐talk between two major signaling pathways, phosphoinositide 3‐kinase (PI3K) and Ras/extracellular signal‐regulated kinase (Erk), which play prominent roles in signaling networks accessed by many cell surface receptors. Signaling mediated by platelet‐derived growth factor (PDGF) receptors, members of the receptor tyrosine kinase class of signal transducers (Claesson‐Welsh, 1994; Schlessinger, 2000), is relevant to wound healing, development, and cancer (Heldin and Westermark, 1999) and is distinguished by potent activation of type IA PI3Ks, which control a host of cellular responses through the production of specific lipid second messengers (Vanhaesebroeck et al, 2001; Engelman et al, 2006; Hawkins et al, 2006). PDGF receptors also enhance signaling through Ras proteins by promoting Ras guanine nucleotide exchange factor (Ras‐GEF) activity and hence an increase in GTP‐bound Ras; in this state, Ras activates the canonical Raf/mitogen‐activated protein kinase (MAPK) or Erk kinase (MEK)/Erk cascade. Erk isoforms are both master integrators of upstream signals and master controllers of transcription factors and other effectors (Kolch, 2000).
Although cross‐talk between the PI3K and Ras/Erk pathways has been studied extensively, a definitive pattern of regulation has not yet emerged, because the two pathways apparently affect each other in various ways and in a context‐dependent manner. Despite striking evidence that PI3Ks bind to and are activated by Ras‐GTP (Rodriguez‐Viciana et al, 1994), the effects of blocking Ras signaling on growth factor‐stimulated PI3K activation range from partial inhibition (Rodriguez‐Viciana et al, 1994; Klinghoffer et al, 1996) to very little or no inhibition (Klippel et al, 1996; Kaur et al, 2006). Numerous reports have also implicated PI3K in the activation or regulation of Ras/Erk signaling, but again with disparate conclusions being drawn (discussed in Duckworth and Cantley, 1997). Although PI3K inhibitors apparently abrogate Erk activation but not Ras‐GTP loading in certain contexts (Cross et al, 1994; Karnitz et al, 1995; Duckworth and Cantley, 1997; King et al, 1997), other reports indicate a role for PI3K upstream of Ras (Hu et al, 1995; DePaolo et al, 1996; Hawes et al, 1996; Wennström and Downward, 1999).
Through a kinetic analysis of PI3K and Ras/Erk signaling in the PDGF receptor system, complemented by quantitative modeling (Tyson et al, 2003; Ma'ayan et al, 2005; Kholodenko, 2006), we find that PI3K signaling affects Erk activation both upstream and downstream of Ras, and that Ras and PI3K account for most, if not all, of the pathways from PDGF receptors to Erk. The magnitudes of the Ras‐ and PI3K‐dependent effects are shown to depend strongly on the stimulation conditions, in ways that are readily explained based on the saturability of the respective pathways and negative feedback mechanisms. Interestingly, cross‐talk from the Ras/Erk pathway to PI3K is shown to be far less significant, an observation that we speculate is tied to the ability of PDGF receptors to directly and potently activate PI3K signaling. Through modeling and computation, the magnitudes of PI3K‐dependent cross‐talk mechanisms are quantified, and predictions are generated to inform further refinement of those mechanisms.
Ras and PI3K make independent contributions to Erk signaling, whereas PI3K is effectively insulated from cross‐talk
To systematically evaluate cross‐talk between the PI3K and Ras/Erk pathways, Erk phosphorylation and PI3K‐dependent Akt phosphorylation in PDGF‐stimulated NIH 3T3 fibroblasts were measured by quantitative immunoblotting for an array of 126 experimental conditions, sampling different combinations of ligand dose, stimulation time, and molecular manipulation (Figures 1 and 2). Considering biological replicates and parallel determination of total Erk and Akt levels, this data set comprises 2772 total measurements. We have confirmed on multiple occasions that the immunoblot imaging system used produces a linear response as a function of epitope loaded, over a sufficiently broad dynamic range (unpublished results).
The results reveal that when either Ras or PI3K is inhibited, by the expression of dominant‐negative (S17N) H‐Ras or incubation with LY294002 compound, respectively, PDGF‐stimulated Erk phosphorylation is partially inhibited (Figure 1A–C). Expression of S17N H‐Ras sequesters Ras‐GEFs and thus affects activation of all Ras isoforms, and we confirmed previously that it prevents PDGF‐stimulated Ras‐GTP loading in our cells (Kaur et al, 2006); we also confirmed that the control compound LY303511, which is pharmacologically similar to LY294002 except that it does not inhibit PI3K (Knight et al, 2006; Gharbi et al, 2007), does not affect Erk phosphorylation significantly (Supplementary Figure S1). The degree of Erk inhibition depends strongly on both stimulation dose and time; no single combination of these conditions is indicative of the relative roles of Ras and PI3K in activating Erk. At short stimulation times, PI3K inhibition significantly affects Erk activation at low but not high PDGF concentrations, in accord with a previous study (Duckworth and Cantley, 1997); however, even at saturating PDGF concentrations, Erk phosphorylation in PI3K‐inhibited cells is less sustained, with significantly lower values than the control at later stimulation times. The effects of Ras inhibition are similar, except that Erk phosphorylation is less sensitive to PDGF concentration and shows a lower initial rate at higher PDGF doses (compare Figure 1B and C). As we will show, these subtle differences in the kinetics between PI3K‐ and Ras‐inhibited cells are indicative of distinct pathways integrated by Erk.
Expression of constitutively active (G12V) H‐Ras yields not only an elevated level of Erk phosphorylation prior to stimulation but also a somewhat reduced responsiveness to PDGF (Figure 1D), consistent with high constitutive expression of MAPK phosphatase (MKP) as reported previously (Kaur et al, 2006); indeed, the basal expression of MKP‐1 is roughly 17‐fold higher in G12V Ras‐expressing cells (results not shown). Thus, Erk signaling in these cells is the net result of both higher phosphorylation and higher dephosphorylation rates. Erk phosphorylation is mediated by both G12V H‐Ras and PDGF‐responsive signaling, which presumably includes activation of endogenous Ras isoforms and pathways that depend on PI3K.
By comparison, the PI3K pathway is not significantly activated through cross‐talk from Ras. Neither expression of dominant‐negative Ras nor application of the MEK inhibitor PD098059 has a significant effect on Akt phosphorylation (Figure 2A–C); the latter treatment is expected to enhance Ras‐GTP loading by relieving feedback desensitization of the Ras‐GEF mSos (Klarlund et al, 1995; Langlois et al, 1995; Waters et al, 1995), but the only effect of the drug is a slight reduction of the peak Akt phosphorylation. Expression of constitutively active Ras results in a higher basal Akt phosphorylation level in unstimulated cells but does not grossly affect PDGF‐stimulated Akt phosphorylation. When Akt phosphorylation is far from saturation, the PDGF‐stimulated response is incremented by the basal level, suggesting an approximately additive relationship; there is also evidence that the saturated Akt phosphorylation response is more sustained in G12V Ras‐expressing cells (Figure 2D). Therefore, whereas inhibiting activation of endogenous Ras did not significantly affect PDGF‐stimulated PI3K/Akt signaling in this study, high levels of Ras‐GTP cooperate with PDGF receptors in a subtle way to activate this pathway (Jiménez et al, 2002; Kaur et al, 2006).
These results suggest two possible Erk activation mechanisms, or combinations thereof: Ras and PI3K promote Erk signaling independently, or PI3K lies upstream of Ras in the same pathway, in which case at least one additional pathway would be responsible for the residual Erk phosphorylation seen when either Ras or PI3K is inhibited. In each of three independent experiments, simultaneous inhibition of Ras and PI3K almost completely abolished PDGF‐stimulated Erk phosphorylation, indicating that Ras and PI3K are responsible for all of the major pathways from PDGF receptors to Erk, and at least one mode of PI3K‐dependent cross‐talk to Erk operates in parallel with Ras (Figure 3). Incidentally, this result also provides further confirmation that Ras function is effectively blocked in our S17N Ras‐expressing cells.
PI3K enhances Erk activation at points both upstream and downstream of Ras
Having established a role for PI3K‐dependent signaling to Erk in parallel with Ras, we sought to determine whether or not PI3K also partially contributes to the PDGF receptor‐mediated activation of Ras. Using a coupled enzymatic assay (Scheele et al, 1995; Kaur et al, 2006), we measured the kinetics of PDGF‐stimulated Ras‐GTP loading for selected conditions (Figure 4). In each of three independent experiments, incubation with LY294002 reduced the initial rate of Ras‐GTP loading at a low PDGF concentration (50 pM), whereas Ras‐GTP loading kinetics were minimally affected at 1 nM PDGF. Interestingly, the plateau level of Ras‐GTP in PI3K‐inhibited cells was similar to that of control cells, even at the lower dose of PDGF.
As alluded to in the previous section, Ras‐GTP loading is affected by a known negative feedback loop that destabilizes Ras‐GEF recruitment, and accordingly, incubation with PD098059 yields a higher Ras‐GTP level (Figure 4). The feedback loop might also explain the effect of PI3K inhibition seen at lower PDGF concentrations. Erk phosphorylation is dramatically reduced under those conditions (Figure 1C), and thus the feedback regulation might be relieved. The net effect of reduced but stable Ras‐GEF recruitment would be a lower initial rate and more sustained level of Ras‐GTP loading.
Multiple negative feedback loops shape the integrated Erk phosphorylation response
In addition to feedback regulation of Ras‐GEF, we also characterized the upregulation of MKP expression induced by Erk and other MAPKs (Brondello et al, 1997; Bhalla et al, 2002; Plows et al, 2002) (Figure 5A and B). MKP‐1 expression is upregulated ∼3‐fold in response to 1 nM PDGF, after a lag time of approximately 15 min. At 30 pM PDGF, which elicits roughly half‐maximal Erk phosphorylation (Figure 1B and C), MKP‐1 expression increases only slightly, suggestive of a signaling threshold for triggering the feedback loop. A similar MKP‐1 response is observed when either Ras or PI3K is inhibited, with a partial reduction in the fold induction (Figure 5B); this outcome can be attributed at least in part to the diminished Erk phosphorylation response, integrated over time, when either of the two major pathways is disrupted. As expected, inhibition of MEK/Erk or c‐Jun N‐terminal kinase (JNK) (but not of p38; results not shown) significantly reduces the induction of MKP‐1 expression, and simultaneous inhibition of MEK and JNK blocks MKP‐1 upregulation completely (Figure 5C).
Parsing the magnitudes of PI3K‐dependent cross‐talk interactions in the PDGF receptor signaling network
Synthesizing the data assembled, a conceptual model of the PDGF receptor signaling network is now clear (Figure 6). Yet, the unique dose–response and kinetic information in the data set allow for a more quantitative description through mathematical modeling and analysis. The relatively coarse granularity of our model, described in detail in the Supplementary information, reflects a careful balance between the level of molecular detail included and both the uncertainty of those details and the ability to specify model parameters. In the context of the model, the magnitudes of the Ras‐ and PI3K‐dependent inputs converging on MEK/Erk determine the saturability of Erk phosphorylation with respect to PDGF dose and also the degree to which the response adapts, and conversely, the observed dynamics can be used to quantify the magnitudes of the inputs.
A total of 34 unspecified parameter values were estimated using a Monte Carlo‐based algorithm that directly and globally compares the model output to the experimental data (Supplementary information); thus, an ensemble of 10 000 parameter sets that fit the data almost equally well was assembled. This approach embraces the inherent ‘sloppiness’ of kinetic models, which prohibits the reliable identification of precise parameter values yet allows for predictions based on collective fits to data (Gutenkunst et al, 2007). The model output is expressed as an ensemble average (Brown and Sethna, 2003; Violin et al, 2008) and compared with the corresponding data (Figure 7). The remaining data used to constrain the algorithm indicate the degree of saturation in the MEK/Erk pathway (Supplementary Figure S2). Taken together, the data force the model to reconcile time‐ and PDGF dose‐dependent features of the network observed under the various experimental conditions tested, including subtle differences between the two control experiments (DMSO vehicle only versus empty pBM‐puro vector). Considering the global constraints imposed by the data and the necessarily simplified nature of the kinetic model, it is our assessment that the model performs quite well; however, one feature that it was not able to reconcile well is the sensitivity of the Erk phosphorylation response in S17N Ras‐expressing cells stimulated with the lower doses of PDGF (Figure 7; bottom row, second plot from left). We speculate that the apparent sensitivity of the response, relative to the model output, is diminished by cell variability. It is known that variability in the cells’ sensitivity to stimulation (the EC50 of the dose–response curve) produces a flatter dose–response curve for the population (Altan‐Bonnet and Germain, 2005). The model also does not capture the transient nature of the S17N Ras/30 pM PDGF curve, suggesting that ligand depletion or other negative feedback effects might need to be characterized and accounted for in future refinements of the model. Indeed, the transient nature of the Akt phosphorylation response at low PDGF concentrations (Figure 2) provides additional evidence for such effects affecting PI3K‐dependent signaling.
Parameter statistics (Supplementary information) show that many parameters tend to be pegged at arbitrarily high values, suggesting ways to simplify the model without a significant change in output. For example, the algorithm consistently picked parameter sets in which the enzymes operate far from saturation. This observation facilitates our analysis of Ras‐ and PI3K‐dependent MEK phosphorylation, in which we examine the catalytic efficiency of each phosphorylation step relative to that of the opposing dephosphorylation reaction (Figure 8A). The results reveal a consistent ratio of PI3K‐ and Ras‐dependent contributions, as indicated by points clustered around a diagonal line, with the Ras‐dependent pathway tending to be more potent under maximal PDGF stimulation. To further quantify these pathways, we formulated a single number, the MEK activation comparator (MAC), which compares their capacities to generate dually phosphorylated MEK (ratio of PI3K‐dependent/Ras‐dependent; see Supplementary information). Over the parameter set ensemble, the median MAC value is 0.22, with quartile values of 0.17 and 0.30; that is, the PI3K‐dependent MEK activation pathway is predicted to be intrinsically less potent than the Ras‐dependent pathway under maximal PDGF stimulation. Importantly, the MAC quantifies these inputs in a way that uncouples them from negative feedback effects. As shown in Figure 7, feedback desensitization attenuates the maximum stimulation of Ras‐GTP by roughly two‐thirds. Adjusting the MAC to account for this reduction, which presumably affects both Ras‐dependent MEK phosphorylation steps equally, yields a median ratio of 1.6 (quartile values of 1.2 and 2.1); that is, Ras‐GEF desensitization renders the PI3K‐dependent pathway more important. This is consistent with the observation that PI3K inhibition affects Erk phosphorylation somewhat more dramatically than inhibition of Ras‐GTP loading.
Upstream of Ras, a similar analysis was performed for PI3K‐dependent and ‐independent modes of Ras‐GEF recruitment, which are characterized by specific model parameters (Figure 8B). Taken together, these modes determine the saturability of Ras‐GTP loading with respect to PDGF dose. Although most of the parameter sets in the ensemble included a significant PI3K‐dependent contribution, a small fraction of them ignored the PI3K‐dependent mode, which we attribute to the subtlety of the data constraining it (the effect of PI3K inhibition on Ras‐GTP accumulation). Defining the GEF recruitment comparator (GRC) as the ratio of PI3K‐dependent/PI3K‐independent modes at low PDGF concentrations (see Supplementary information), the median GRC value for the ensemble is 1.4, with quartile values of 0.7 and 2.0. This analysis suggests that the two Ras activation modes are almost equally important.
Model predictions and potential for model refinement
We have shown how data‐driven analysis of a kinetic model can be used to quantify inputs to signaling pathways. Another use of this approach is to generate hypothetical predictions with an eye toward future experiments. Here, we focus on the PI3K‐dependent cross‐talk impinging upstream and downstream of Ras (Figure 9). Although both mechanisms are blocked by PI3K inhibition, the model ensemble predicts unique kinetic signatures that might be expected if either mechanism were silenced selectively. At lower PDGF concentrations, it is predicted that blocking either mechanism would yield Erk phosphorylation levels that are intermediate between control and PI3K‐inhibited conditions, but the key discriminator is the Erk phosphorylation kinetics observed with a high PDGF dose; there, if the cross‐talk mechanism upstream of Ras is blocked, the kinetics are predicted to be similar to control conditions, whereas if the mechanism downstream of Ras is blocked, the kinetics are predicted to be similar to the PI3K‐inhibited case (Figure 9). The basis for this result stems from the saturability of Ras‐GEF recruitment, wherein high PDGF concentrations stimulate maximal Ras‐GTP loading even in the absence of PI3K signaling.
A unique data set was collected to systematically characterize multiple cross‐talk interactions between the PI3K and Ras/Erk pathways and regulatory feedback loops in the PDGF receptor signaling network. PI3K is strongly activated and plays a prominent role in this system, and Ras‐ and PI3K‐dependent pathways converge as distinct inputs to Erk activation that we have quantified experimentally and through computational modeling. At lower PDGF concentrations, Erk activation qualitatively follows AND logic, with both inputs required for a robust response, whereas it follows OR logic at higher PDGF concentrations, wherein either pathway is sufficient. The basis for these observations lies in the saturability of the pathways. Half‐maximal PDGF receptor phosphorylation is elicited by ∼0.5–1 nM PDGF‐BB (Park et al, 2003), concentrations that yield maximal activation of PI3K and Ras in our cells; downstream of Ras, we also find partial saturation of MEK or Erk. The magnitudes of the inputs to the Erk pathway also affect adaptation of the response imposed by negative feedback regulation.
In contrast, cross‐talk from the Ras/Erk pathway to PI3K is far less significant, which might be considered surprising in light of the interaction between Ras‐GTP and PI3K catalytic subunits. Apparently, Ras does not appreciably complement the already tight binding between PDGF receptors and PI3K (Kazlauskas and Cooper, 1990), whereas Ras might play a more significant role in conjunction with receptors/adaptors that bind PI3K regulatory subunit with lower affinity (Kaur et al, 2006).
To advance this line of inquiry, the PI3K‐dependent intermediates that execute the cross‐talk to Erk will need to be better and more quantitatively characterized. One candidate pathway is the phosphorylation of Raf‐1 on Ser338, which has been implicated in Raf‐1 activation. Phosphorylation of Raf‐1 on Ser338 is catalyzed by isoforms of p21‐activated kinase (Pak), which can be activated downstream of PI3K and either Rac or Cdc42 (Chaudhary et al, 2000; Sun et al, 2000). Raf‐1 mutants that cannot bind Ras‐GTP can be activated in cells by co‐expression of Pak3 and constitutively active Cdc42, suggesting a Ras‐independent pathway (Sun et al, 2000), whereas other evidence suggests a role for Ras in the activation of Raf‐1 by Rac/Cdc42 and Pak (Li et al, 2001; Zang et al, 2002). Consistent with this literature, we found that Raf‐1 phosphorylation on Ser338 strongly depends on both PI3K and Ras in our cells (results not shown); from the standpoint of cross‐talk from PI3K to Erk, this finding is inconclusive, because it does not account for the Ras‐independent pathway in our cells. Other aspects of Raf‐1 regulation might account for the pathway, as there is ample evidence that Ser338 phosphorylation is not a strong correlate of Raf‐1 activity (Sun et al, 2000; Chiloeches et al, 2001; Oehrl et al, 2003; Beeser et al, 2005); regulation of other Raf isoforms and other putative MEK kinases, such as Pak (Beeser et al, 2005; Park et al, 2007) and 3‐phosphoinositide‐dependent protein kinase 1 (Sato et al, 2004), might also be considered. It is also worth noting that such PI3K‐dependent modes of Erk activation could mask a subtle negative regulation of the Erk pathway, through Akt for example (Zimmermann and Moelling, 1999); however, overexpression of Akt‐1 does not affect PDGF‐stimulated Erk phosphorylation in our cells, despite a two‐fold increase in phosphorylated Akt (results not shown). Finally, although at least one important mode of PI3K‐dependent cross‐talk to Erk is Ras independent, PI3K lipid products apparently contribute upstream of Ras as well, perhaps through direct recruitment of Ras‐specific GEF activity or through mediators such as centaurin‐α1 (Hayashi et al, 2006).
If such mechanisms are to be tested, mathematical model predictions provide kinetic signatures that could help validate the point of action of a particular pathway from PI3K to Erk. For example, if a particular kinase contributing to the Ras‐independent Erk activation pathway were blocked, the model predicts that the time courses of Erk phosphorylation at both low and high PDGF concentrations would be perturbed relative to control (Figure 9). Furthermore, the extent of perturbation relative to complete PI3K inhibition, especially at high PDGF concentrations, reveals the degree to which the kinase is playing a major role in the pathway. Such an analysis would require quantitative measurements of the type presented here. Importantly, the extent of Erk inhibition in the dominant‐negative Ras background would lend further evidence and opportunity for model validation. In turn, characterization of cross‐talk mechanisms along these lines allows both conceptual and quantitative models to be refined. The data set generated here, together with new data, will be instrumental in iteratively testing more detailed models of the PDGF receptor signaling network.
Another issue that needs to be confronted is the generality of signal‐transduction networks across receptor and cell types. Certainly, major disparities in the regulation of signaling in different contexts, even among receptors of the same class, should be expected. Indeed, it is instructive to compare the PDGF receptor network to the proposed regulatory structure ascribed to epidermal growth factor (EGF) receptor, which has been the subject of numerous quantitative modeling efforts (Wiley et al, 2003). Recent experimental and modeling work has provided compelling evidence that EGF receptor signaling through PI3K and Erk are coupled primarily through regulation of Gab1 (Kiyatkin et al, 2006). According to that model, phosphorylated Gab1 recruits PI3K, and Gab1 responds to 3′ phosphoinositides, constituting a positive feedback loop; Gab1 also aids in the activation of Ras. Further, Gab1 and therefore PI3K are subject to a negative feedback loop mediated by Erk. The evidence suggests that the inability of EGF receptors to interact directly with PI3K demands the utilization of Gab1 and thus dictates the important features of that network. In contrast, PDGF receptors activate PI3K directly and strongly, perhaps explaining why PI3K is not subject to regulation by Erk in our cells; other evidence suggests that Gab1 is prominently utilized to recruit PI3K in response to EGF, but not PDGF, in the same cell background (Zhang et al, 2002). Clearly, control at the network level is highly context dependent, and more quantitative approaches will be needed to unravel the mechanisms that distinguish particular signaling systems.
Materials and methods
All tissue culture reagents were from Invitrogen (Carlsbad, CA). Human recombinant PDGF‐BB was from Peprotech (Rocky Hill, NJ). Antibodies against Ras (Y13‐259, agarose conjugated), Akt‐1/2 N terminus, Erk‐1 C terminus, and MKP‐1 were from Santa Cruz Biotechnology (Santa Cruz, CA). Phospho‐specific antibodies against Akt pSer473, Erk pThr202/pTyr204, and Raf‐1 pSer338 were from Cell Signaling Technology (Beverly, MA). Pharmacological inhibitors were from Calbiochem (San Diego, CA); where applicable, cells were pre‐incubated with the inhibitor for 30–60 min prior to PDGF stimulation. Unless otherwise noted, all other reagents were from Sigma‐Aldrich (St Louis, MO).
Plasmids, retroviral infection, and cell culture
The retroviral vectors bearing S17N and G12V H‐Ras, cloned into pBM‐IRES‐Puro, and their introduction into the ecotropic ϕNX packaging cell line by the calcium phosphate precipitation method, were described previously (Kaur et al, 2006). NIH 3T3 fibroblasts (American Type Culture Collection, Manassas, VA) were cultured at 37°C, 5% CO2 in Dulbecco's modified Eagle's medium supplemented with 10% fetal bovine serum, 2 mM l‐glutamine, and the antibiotics penicillin and streptomycin. Where applicable, NIH 3T3 cells were serially infected with retrovirus bearing empty vector, S17N H‐Ras, or G12V H‐Ras and selected using puromycin prior to each experiment, as described previously (Kaur et al, 2006).
Lysate preparation and biochemical assays
Cells were serum‐starved for 4 h prior to stimulation. Detergent lysates were prepared for quantitative immunoblotting, and immunoblots were performed using enhanced chemiluminescence, as described previously (Park et al, 2003). Blots comparing lysates prepared on the same day, representing either different inhibitor treatments or different cell variants and respective control conditions, were performed in parallel and exposed at the same time. The Bio‐Rad Fluor S‐Max system, which gives a linear response with respect to light output, was used, and band intensity was quantified using local background subtraction. Total cellular protein concentrations were determined by Micro BCA assay (Pierce, Rockford, IL).
To determine the amounts of Ras‐bound GTP (from all Ras isoforms), eluted from anti‐Ras immunoprecipitates (Y13‐259, agarose conjugated), the coupled nucleoside‐5′‐diphosphate kinase/luciferase assay (Scheele et al, 1995) was used; the determination of the sum of GDP+GTP was achieved by first converting GDP to GTP. These assays were performed as described previously (Kaur et al, 2006).
Data normalization and statistical analysis
Biochemical data were first normalized by an appropriate loading control as stated in the figure legends. To evaluate the consistency of relative trends across independent experiments, the data were further normalized as follows (see also Supplementary Figure S3). Suppose there are n experiments for which data are collected at m time points. During each of the n experiments, the same control is run. Let Xij denote the experimental measurement for the control in the ith experiment at the jth time point. Often the quantity of interest Yij (for example, the fraction of Erk or Akt in the phosphorylated state) is related to Xij by an unknown scale factor, αi. Under ideal conditions, the control would not vary from experiment to experiment. Therefore, we sought the set of αi values that minimize the suitable error metric F, defined as follows.
where j is the mean value that results for time point j. The function F is proportional to the mean coefficient of variation for the control time series. The minimization is subject to a constraint that eliminates the trivial solution, αi=0 for all i. The time course of cells stimulated with 1 nM PDGF under control conditions (DMSO or empty vector) was used to determine the normalization factors, and each factor was applied to the remaining data collected on the same day (control conditions with lower PDGF concentrations and either inhibitor or Ras variant expression conditions, measured in parallel). Where applicable, normalized means for control and perturbed conditions were compared by Student's t‐test, using the most conservative model (two‐tails, unequal variance).
Kinetic modeling and analysis
A mathematical model was constructed based on known or plausible signaling mechanisms, with appropriate simplifications, as described in detail in the Supplementary information. Briefly, PDGF receptor binding, dimerization, and endocytosis, and the production of 3′ PI lipids by receptor‐recruited PI3K, are modeled as previously described (Park et al, 2003; Schneider and Haugh, 2005). Activated receptors and 3′ PI lipids mediate the recruitment of Ras‐GEF activity and Ras‐GTP loading; for simplicity, the roles of adaptors such as Grb2 and Shc are implicit. Ras‐GTP and 3′ PI lipids independently mediate the activation of one or more MEK kinases, and MEK and Erk are sequentially activated by dual phosphorylations. Activated Erk governs negative feedback through desensitization of Ras‐GEF and through upregulation of MKP activity.
The model has 18 time‐dependent variables and a total of 44 rate constants or other adjustable parameters, of which 10 were assigned fixed values based on previous work. Remaining parameters were specified using a computational scheme based on the Metropolis algorithm, which is well suited for kinetic parameter estimation constrained by biochemical data (Brown and Sethna, 2003; Violin et al, 2008). Our variation of the method, implemented in MATLAB (The MathWorks; codes available upon request), works as follows (see Supplementary information for details). Starting with initial parameter values, the model output is computed and aligned with the data, and the quality of fit is evaluated. Next, new values are chosen according to normal distributions centered on the previous values; the width of the distribution is a parameter of the algorithm. If specified error criteria are met, the new parameter set is accepted and used as the nexus for choosing the next parameter values; otherwise, the procedure is repeated with the previous parameter set. If the starting parameter set (determined based on a separate error minimization subroutine) is nearly optimal, then the parameter sets generated in this manner constitute an ensemble that fit the data almost equally well. After compiling the ensemble, the model output is recalculated for each parameter set, and at each time point, an ensemble mean and standard deviation are calculated.
We thank Timothy Elston (Department of Pharmacology, University of North Carolina‐Chapel Hill) for helpful discussions on data normalization and parameter estimation practices. This study was supported by National Institutes of Health grant R01‐GM067739.
Conflict of Interest
The author declares that he has no conflict of interest.
Supplementary text S1 and supplementary Figures S1, S2 and S3 [msb20094-sup-0001.pdf]
PDGF receptor network
SBML (level 2, version 4) model [msb20094-sup-0002.zip]
- Copyright © 2009 EMBO and Nature Publishing Group