Precisely Ordered Phosphorylation Reactions in the p38 Mitogen-activated Protein (MAP) Kinase Cascade*

Background: MAPK cascades are signaling modules that function as switch generators. Results: In vitro phosphorylation in the p38 MAPK cascade tracked by LC-MS/MS revealed specific phosphorylation intermediates at each level. Conclusion: The p38 MAPK cascade reactions occur through intermediates MEK6/ST* and p38α/TY*. Significance: The precise order of reactions may contribute to the diverse kinetic outputs of the cascades, including those with large Hill coefficients. The MAP kinase cascades, composed of a MAP3K, a MAP2K, and a MAPK, control switch responses to extracellular stimuli and stress in eukaryotes. The most important feature of these modules is thought to be the two double phosphorylation reactions catalyzed by MAP3Ks and MAP2Ks. We addressed whether the reactions are sequential or random in the p38 MAP kinase module. Mass spectrometry was used to track the phosphorylation of the MAP2K MEK6 by two MAP3Ks, TAO2 and ASK1, and the subsequent phosphorylation of p38α by MEK6/S*T* (where S (Ser) and T (Thr) are the two phosphorylation sites and * denotes phosphorylation). Both double phosphorylation reactions are precisely ordered. MEK6 is phosphorylated first on Thr-211 and then on Ser-207 by both MAP3Ks. This is the first demonstration of a precise reaction order for a MAP2K. p38α is phosphorylated first on Tyr-182 and then on Thr-180, the same reaction order observed previously in ERK2. Thus, intermediates were MEK6/ST* and p38α/TY*. Similarly, the phosphorylation of the p38α transcription factor substrate ATF2 occurs in a precise sequence. Progress curves for the appearance of intermediates were fit to kinetic models. The models confirmed the reaction order, revealed processivity in the phosphorylation of MEK6 by ASK1, and suggested that the order of phosphorylation is dictated by both binding and catalysis rates.

MAP 2 kinase cascades mediate decision-making processes throughout eukaryotic life. Changes in cell fate leading to differentiation, proliferation, and apoptosis in response to hormones, as well as metabolic changes in response to stress, involve MAP kinase pathways (1)(2)(3)(4)(5). Each MAP kinase module has a canonical architecture of three kinases, a MAP3K, a MAP2K, and a MAPK. The cascades catalyze a sequence of phosphorylation reactions. MAP3Ks phosphorylate two Ser/ Thr residues in MAP2Ks (6), and MAP2Ks phosphorylate MAPKs on a Tyr and a Thr residue (7), such that MAP2Ks are true dual (Tyr and Ser/Thr) specificity kinases. The doubly phosphorylated MAP kinase is the product of the cascade, a proline-directed Ser/Thr kinase (8). Four distinct MAP kinase cascades have been elucidated, each named for the MAPK-activated, ERK1/2 (extracellular signal-regulated kinase 1 and 2), p38 MAP kinase, JNK (c-Jun kinase), and ERK5 pathways (4). These modules are activated by a variety of inputs, including other kinases, referred to as MAP4Ks, interactions with small G-proteins (such as Ras), or scaffold proteins (such as JNKinteracting protein (JIP)). The modules also have numerous substrates, consisting primarily of transcription factors and other protein kinases. Thus, apparently, the three-tiered modules "plug in" to a variety of cellular contexts. The cascades are further capable of a variety of kinetic outputs as measured by the dose-response activity of the MAPK. In oocytes, the ERK2 MAP kinase module generates highly cooperative or sigmoid responses (9,10). Responses observed in other cell types can be more graded (11)(12)(13). Signal response can further be modified by the strength and duration of inputs (12,14).
Several recent studies have addressed the dose-response behavior of MAP kinase modules in cells (13)(14)(15)(16)(17). Although MAPK phosphorylation occurs through a specific intermediate (13,18), the phosphorylation order in MAP2Ks has not been addressed. The strong sigmoid properties of the cascade occurring in oocytes is thought to be derived from the phosphorylation reactions catalyzed by MAP3Ks and MAP2Ks (19). There are two models (20). One is that both the phosphorylation reactions occur at random but both are required for the activity of the MAP2K (19). A second possibility is that the reactions occur in a precise sequence (as proposed in Fig. 1), and the activity is associated with the second phosphorylation event. In either case, sigmoid behavior should be possible when the reactions are distributive (nonprocessive) (19). Some support for the sequential model is in the literature. It has been established that the MAPK ERK2 is phosphorylated first on Tyr-185 to give ERK2/TY* and then on Thr-183 to give ERK2/T*Y* (18,19,21,22) (where T (Thr) and Y (Tyr) are the two phosphorylation sites and * denotes phosphorylation). Engelberg and colleagues (23) have shown that in the yeast Hog1 pathway, Tyr phosphorylation is not required for activity, but is essential for the salt stress activation of Hog1 in cells (see also Refs. 24 and 25). The phosphorylation of MAP2Ks has not been well studied. Mutational analysis of the role of phosphorylation sites in the activity and activation of MAP2Ks gave mixed results (6,26,27), except that in vitro studies of MEK1 and MEK6 (MAP/ERK kinase 1 and 6) showed that one alanine mutant is phosphorylated more easily than the other (MEK1/AS more than MEK1/SA; MEK6/AT more than MEK6/SA) (27,28). These data are suggestive that phosphorylation of MAP2Ks may be sequential.
Here we phosphorylated MEK6 in vitro, using two different activators, the MAP3K ASK1 (apoptosis signaling kinase-1) (29) and TAO2 (thousand and one kinase-2) (30), tracking the appearance and disappearance of all MEK6 phosphorylation states by LC-MS/MS. We found that MEK6 phosphorylation goes through a precise intermediate, MEK6/ST* (Fig. 1). We also found that p38␣ MAP kinase, the MEK6 substrate, is phosphorylated through p38␣/TY*, as anticipated from studies of ERK2. We further confirmed that p38␣ phosphorylates the transcription factor ATF2 in a specific order (31). Progress curves for the appearance of phosphorylation intermediates and products were fit to kinetic models for the double phosphorylation reactions, which validated the apparent sequence of phosphorylation at each level. The double phosphorylation of MEK6 by TAO2 occurred processively, whereas the rest of the reactions were distributive.

EXPERIMENTAL PROCEDURES
Protein Reagents-MEK6, TAO2, ASK1, p38␣, and ATF2 proteins were expressed and purified following standard procedures (32). A rat cDNA encoding MEK6 was subcloned into a modified pET29b (Novagen) expression vector possessing a C-terminal tobacco etch virus protease recognition sequence and a His tag. This expression vector was modified to encode the kinase-dead MEK6/K82M. Full-length MEK6 and MEK6/ K82M expression and purification were conducted essentially as described in Ref. 33. Viral stock harboring TAO2 (residues 1-451) cDNA was transfected into Sf9 insect cells for protein expression as described previously (28). Human ASK1 (residues 659 -951) was synthesized in Escherichia coli with an N-terminal His tag using the expression vector pLIC (obtained from Stefan Knapp, Structural Genomics Consortium at Oxford (34)). A pET14b plasmid encoding p38␣ described previously (35) was modified to produce kinase-dead p38␣/K53M using the Stratagene QuikChange kit. Wild-type p38␣ and p38␣/ K53M proteins were expressed and purified as described in Ref. 36. A pGEX(KG) plasmid encoding GST-tagged human ATF2 (residues 33-254) was obtained from Melanie Cobb, University of Texas Southwestern (UT Southwestern) Medical Center, expressed in bacteria, and purified on a glutathione-Sepharose column (GE Healthcare) following the manufacturer's protocol.
Phosphorylated Active Kinases-TAO2 and ASK1 are constitutively phosphorylated as expressed in either insect cells (TAO2) or bacteria (ASK1). Doubly phosphorylated p38␣ (p38␣/T*Y*) was obtained by coexpression of p38␣ with MEK6/DD (MEK6/S207D,T211D) in E. coli as described previously (37). MEK6/S*T* was obtained by incubation of ASK1 and MEK6 (in the molar ratio of 1 to 100) at 30°C for 30 min using standard kinase reaction conditions (see below). ASK1 was removed using nickel-Sepharose (GE Healthcare). The phosphorylation states of these proteins were verified by LC-MS/MS (see below).
LC-MS/MS Analysis-HPLC-MS/MS analysis was performed on an LCQ DECA XP ion trap mass spectrometer (ThermoFinnigan) in which the HPLC was coupled in-line to an orthogonal electrospray ionization source. Integration under ion traces corresponding to mass ranges for activation loop peptides were used to acquire raw MS detector responses (see Fig. 2). Time courses for the appearance of peptides were conducted in triplicate, and the different experiments were scaled by the total ion current over an elution window encompassing peptides from the digestion (38).
MS/MS spectra were acquired in a data-dependent mode (recording the masses and fragmentation patterns of the strongest to weakest ions) in each MS scan (ThermoFinnigan). MASCOT software (Matrix Science Ltd., London, UK) (39) was used for identifying peptides from their MS spectrum. Better than 80% of the expected proteolytically generated peptides for MEK6 and p38␣ were identified by the software and included the activation loop-derived peptides.
Cascade Modeling-Kinetic models for the phosphorylation reactions were defined in DynaFit (40,41). The DynaFit software uses the Levenberg-Marquardt algorithm to perform nonlinear least squares regression of progress curve data. Global minimization was achieved by a differential evolution of trial parameter sets (41). Trial parameter values were taken from those obtained for cellular assays of the ERK2 cascade (15). The program outputs standard error for each variable and a correlation matrix for the model parameters. All likely kinetic steps were included in the initial models, including separate binding steps leading to different chemical outcomes (supplemental Fig. S1). However, backward dissociation reactions were eliminated because they were too slow to affect the model, leading to an irreversible mechanism throughout (41,42). Other steps that were too fast to be kinetically relevant were eliminated. Due to correlation between variables, one constant was fixed at each level (based on values from the literature), thus improving the fit and the calculated standard error. Equations used are listed in supplemental Table S2. For ASK1 phosphorylation of MEK6, 224 data points were used to model eight variable parameters, giving a 28:1 data/parameter ratio. For MEK6 phosphorylation of p38␣, 156 data points were used to model five variable parameters, giving a 31:1 data/parameter ratio. In general, all reactions were assumed to be distributive (with dissociation between phosphorylation reactions) (43).

RESULTS
To determine the order of phosphorylation reactions at each level of the p38 MAP kinase cascade, the phosphorylation of MEK6 by two MAP3Ks and the phosphorylation of p38␣ by MEK6 were studied separately. The appearance of phosphorylation at each site of either MEK6 or p38␣ was measured as a function of time of incubation with its activating enzyme. The position and extent of phosphorylation was determined by HPLC-MS/MS of proteolysis-derived peptides. The phosphorylation of ATF2 by p38␣ was measured similarly.
Phosphorylation Intermediates in the Dual Phosphorylation of the MAP2K MEK6 Were Identified by LC-MS/MS-MEK6 was phosphorylated by the MAP3K TAO2. MS/MS spectra of synthetic peptides corresponding to expected chymotrypsinderived activation loop peptides (LVDSVAKTIDAGKPY, supplemental Table S1) in alternative states of phosphorylation were measured. Then chymotrypsin-derived MEK6 activation loop peptides were separated to base line by reverse phase HPLC and identified by tandem mass spectrometry using the MASCOT software (Matrix Science, Inc.) (39) (Fig. 2). After 10 min, substantial quantities of three phosphorylation states were observed (MEK6/ST, MEK6/ST*, MEK6/S*T*), with MS/MS spectra corresponding to those obtained for standard peptides. A small amount of the disfavored intermediate, MEK6/S*T, was also seen. Standard peptides run under similar conditions were used for approximating peptide concentrations (see "Experimental Procedures" and the legend for Fig. 3).
Dual Phosphorylation of MEK6/K82M Occurs through a Single Intermediate, MEK6/K82M/ST*-MEK6/K82M phosphorylation by TAO2 and ASK1 was tracked over 30 min (Fig. 3, A  and B  Dual Phosphorylation of p38␣ by MEK6 Occurs Primarily through the Intermediate p38␣/TY*-MEK6/S*T* was used to phosphorylate kinase-dead p38␣ (p38␣/K53M). HPLC-MS of trypsinized p38␣ was recorded as above to generate progress curves for p38␣ activation loop phosphorylation. The p38␣ activation loop peptides (supplemental Table S1) derived from p38␣/TY, p38␣/TY*, p38␣/T*Y, and p38␣/T*Y* were observed (supplemental Fig. S2). The time course (Fig. 3C) reveals a single dominant intermediate, p38␣/TY*, with small amounts of the alternate intermediate, p38␣/T*Y, also present. This order of phosphorylation, through p38␣/TY*, is the same as that observed in the MAPK ERK2 (13,18). Over the first 2 min, the appearance of p38␣/TY* mirrors the disappearance of p38␣/ TY. Further, the appearance of p38␣/T*Y* occurs with the same apparent rate as the disappearance of p38␣/TY*. These data indicate that the flux of phosphorylation occurs through the intermediate p38␣/TY*. In contrast to the phosphorylation of MEK6, there is a very strong lag in the formation of p38␣/ T*Y*, indicative of a distributive mechanism (little or no processivity). The substrate p38␣ approaches the base line over the time course, in contrast with the MEK6 phosphorylation.
p38␣ Phosphorylates ATF2 through ATF2/TT* in Vitro-We tested whether ordered phosphorylation occurs at the level of MAPK pathway output. Wild-type p38␣ phosphorylates a native substrate, ATF2, on two threonine residues, Thr-69 and Thr-71 (44). Using the same methods as for MEK6 and p38␣ phosphorylation, we tracked tryptic peptides encompassing Thr-69 and Thr-71 (supplemental Table S1) of GST-ATF2 and observed both monophosphorylated intermediates, ATF2/T*T and ATF2/TT* (Fig. 3D). ATF2/TT* is the dominant intermediate observed over the time course, although the alternative intermediate is also present. There is a strong processive component, as indicated by the simultaneous appearance of ATF2/ TT* and ATF2/T*T* as seen previously through in vitro (31) and cellular ATF2 phosphorylation assays (45). The rate of phosphorylation of ATF2 by p38␣ is significantly faster than the rate of phosphorylation of p38␣ by MEK6.
Phosphorylation in the p38␣ MAPK Module Follows a Sequential Mechanism-Given the above data, it appears that the phosphorylation events through the cascade composed of ASK1 (or TAO2), MEK6, and p38␣ follow a sequential mechanism. ASK1 or TAO2 phosphorylates MEK6 through MEK6/ ST* (in the sequence 206 DSVAKTI 212 ); MEK6/S*T* phosphorylates p38␣ through p38␣/TY* (in the sequence 183 TGY 185 ). Finally, p38␣/T*Y* phosphorylates ATF2 through ATF2/TT* Traces represent total ion current normalized selected ion response curves for each proteolytically derived activation loop phosphopeptide. Experiments were performed in triplicate with the error bars representing one standard deviation. Unphosphorylated (black), monophosphorylated (red and green), and doubly phosphorylated (blue) activation loop peptides were tracked for 30 min. In A and B, the sum of peptide signal apparently exceeds the starting value, probably due to errors in the total ion current normalization (38) and matrix ionization effects (55). In C and D, these effects were minimal.
(in the sequence 67 DQTPTTR 73 ). The sequential mechanism is summarized in Fig. 4A.
Kinetic Models Describe MEK6 and p38␣ Phosphorylation-Several questions arise from the shapes of the progress curves. Does the flux of the reactions track through the apparent observed phospho-intermediates? Further, what is the origin of the distinctive shapes of the curves at the two levels? Does this arise from the existence of different kinetic steps, or is there just a variation in the rates of the same steps? Also, is the reaction specificity due to rates of binding or catalytic events? To gain qualitative answers to these questions, kinetic models were fit to the phosphorylation progress curves. As described under "Experimental Procedures," all potential binding and catalytic steps were initially considered; then the reactions that contributed kinetically were retained in the final model. At the MAP2K level, we focused on the ASK1 phosphorylation of MEK6. At both levels of the cascade, complexes leading to different phosphorylation reactions were distinguished by underlining the amino acid in the active site (for example, the complex ASK1-MEK6/ST in Fig. 4B leads to chemistry on Thr). Further, catalysis and dissociation steps were held separate. In general, the back reactions were too slow to affect the model. At the MAPK FIGURE 4. The sequence of cascade phosphorylation events. A, signal propagation through the p38 MAP kinase cascade follows double phosphorylation reactions that apparently occur in a preferred sequence, highlighted in colored arrows. The MAP2K MEK6 is phosphorylated first on Thr-211 to give MEK6/ST* and then on Ser-207 to give MEK6/S*T*. p38␣ is phosphorylated first on Tyr-182 to give p38␣/TY* and then on Thr-180 to give p38␣/T*Y*. B and C, kinetic models were generated for ASK1 phosphorylation of MEK6 (B) and MEK6 phosphorylation of p38␣ (C). The phosphorylation site involved at each step is shown in red, and underlining indicates binding in the active site. Protein complexes leading to catalytically competent interactions were defined. This is followed by an irreversible phosphorylation step and rebinding of the monophosphorylated substrate. A second phosphorylation step is similarly modeled, leading to a final, double phosphorylated product. In both models, one set of reactions, initiated by binding MEK6/ST in (B) and p38␣/TY in (C), is unfavored and shown in gray. Only reactions contributing kinetically to the model are included, to reduce the number of variables. In B, a processive component was required (k 9 ) to account for the rapid appearance of the doubly phosphorylated MEK6/S*T*. In C, it was not necessary to separate binding and catalysis in the formation of p38␣/TY* to fit the phosphorylation data. Models were evaluated using the DynaFit software (40,56). Calculated progress curves from the models for the ASK1 phosphorylation of MEK6 (D) and MEK6 phosphorylation of p38␣ (E) are shown superimposed on the experimental data.
level, the persistence of p38␣ and p38␣/TY* was modeled as a slow loss of competence of all p38␣ proteins to be phosphorylated over the duration of the time course experiment (see Fig. 4E legend for p38␣ and see also supplemental Table S2). At the MAP2K level, no such additional terms were required to obtain a fit. The final models are given in Fig. 4, B and C, and the overlays of the calculated curves to the data are given in Fig. 4, D and E. Reactions leading to the unfavored intermediates, MEK6/S*T and p38␣/T*Y, were much slower (shown in gray in Fig. 4, B and C). The derived kinetic constants for both MEK6 and p38␣ phosphorylation are given in supplemental Tables S3 and S4, respectively. The catalytic constants from our data fit approximately with those presented by Fujioka et al. (15) (11 min Ϫ1 ) for phosphorylation of MEK and ERK in cells.
The final reaction schemes are slightly different. At the MAP2K level, it was necessary to include a term for direct conversion of ASK1-MEK6/ST* to MEK6/S*T* (k 9 ), thus accounting for the apparent processivity in Fig. 3B, where the intermediate and product appear simultaneously. In the phosphorylation of p38␣, the tyrosine phosphorylation was too fast to affect the model, allowing the binding and catalytic steps to be conflated.
Kinetic Modeling Confirms Flux through Specific Intermediates-Fitting confirms that the flux to doubly phosphorylated MEK6 occurs primarily through MEK6/ST* (with k 2 much smaller than k 1 , supplemental Table S2). Likewise, the phosphorylation of p38␣ by MEK6/S*T* was fit to the model shown in Fig. 4C, with the resulting calculated curves validating p38␣/TY* as the predominant intermediates, again with k 2 Ͻ k 1 (supplemental Table S2).
Reactions at the MAP2K but Not MAPK Level Are Partially Processive-At the MAP2K level, the reaction is partially processive, without dissociation between catalytic steps. The processivity was modeled as an additional kinetic term for the conversion from ASK1-MEK6/ST* to ASK1-MEK6/S*T* without complex dissociation (k 9 , supplemental Table S2, Fig. 4B). Models lacking a processive step resulted in a much poorer fit (data not shown). In p38␣ phosphorylation, the reaction is distributive (nonprocessive). This is apparent from the sudden appearance p38␣/TY* relative to p38␣/T*Y*. The slowest reaction in p38␣/T*Y* formation is the last step, and the first step, k 1 , is much faster (k 1 /k 5 ϳ 20), further explaining the buildup of intermediates. No improvement in the model fit was obtained by including a processive reaction term (data not shown). Thus, the difference in shapes of the two different sets of curves is significantly due to an added processive step in the ASK1 phosphorylation of MEK6.
Specific Intermediates Result from k cat as well as Binding Effects-In the ASK1 phosphorylation of MEK6, the binding constant k 1 is greater than k 2 (2 Ͼ 0.04 M Ϫ1 min Ϫ1 ). Similarly, the k cat constant k 3 is greater than k 4 (11 Ͼ 0.7 min Ϫ1 ). Apparently, both binding and catalysis contribute to the reaction specificity. In the phosphorylation of p38␣, k 1 Ͼ k 2 ; k cat for the conversion of p38␣/TY to p38a/TY* is too fast to model, giving rise to the same conclusion, that both binding and catalysis contribute to the reaction specificity.
The Rate of Reaction Increases through the Cascade-The progress curve data allow us to consider overall rates for each step in the full MAPK cascade. As a first approximation, the slowest step leading to doubly phosphorylated product is the rate-limiting step. If we only consider the dominant phosphorylation pathways, the slowest step in the formation of MEK6/ S*T* is k 1 at about 2 M Ϫ1 min Ϫ1 . The slowest step to form p38␣/T*Y* is k 5 at about 7 min Ϫ1 . Phosphorylation of ATF2 is fastest of all, at about 500 min Ϫ1 , as estimated from the initial rate of phosphorylation (Fig. 3).

DISCUSSION
Phosphorylation order in MAP kinase modules could theoretically be either random or sequential (20). Here we show the p38 MAP kinase module clearly follows the sequential model. The order of reactions is well defined at both levels of the cascade. The MAP2K MEK6 is phosphorylated primarily through the intermediate MEK6/ST*. The phosphorylation of p38␣ goes through a dominant intermediate, p38␣/TY*. This is the same order that occurs in the MAPK ERK2 (18,21). At both levels, the order is apparently determined by a difference in initial rates leading to different chemical outcomes (k 1 /k 2 Ͼ 1). We also confirmed that ATF2 is phosphorylated through the intermediate ATF2/TT* by p38␣ as reported by Waas et al. (31).
The basis for the observed order of phosphorylation events at either level of the cascade is not understood in molecular detail. Our earlier studies of MEK6-TAO2 interaction showed that mutation of either Ile-212 or Asp-206 in the activation loop of MEK6 (in the sequence 206 DSVAKTI 212 ) reduced phosphorylation by TAO2 (28). Mutational analysis of TAO2 revealed a cluster of basic residues (Lys-120, Arg-221, and Lys-222) that is important for activity (marked as ϩ in Fig. 1). This led to a model in which the isoleucine (Ile-212) in MEK6 binds to the Pϩ1 pocket (46) of TAO2, and the aspartic acid (Asp-206) binds the basic cluster, thus dictating the reaction order.
No structures or molecular models are available to indicate why tyrosine in the MAPK activation loop sequence TXY is phosphorylated first by MAP2Ks. However, the progress curves and model fitting presented here show that MEK6 binding to p38␣ at Tyr-185 leads to rapid catalysis in comparison with p38␣ Thr-183 binding. This is also the case in MAP kinase phosphatases, which dephosphorylate first on Tyr (47,48). It is also unknown how ATF2 is phosphorylated in a preferred order. Further data will be required to fully understand the underlying mechanisms driving the reaction orders established here.
At the MAPK level, the role of a phosphotyrosine intermediate has been in dispute, with the suggestion that either it is required for activity (49) or it is a stepping stone in MAPK activation (23). Here we have shown that the p38␣ phosphorylation occurs through p38␣/TY*. This result offers support to the idea that phospho-tyrosine formation is a stepping stone in MAPK activation. The stepping stone idea gains strength from the lack of Tyr residues in the ERK2 homolog ERK3 (4), suggesting that the ERK3/S* (occupying the homologous position to Thr in ERK2 and p38␣) is all that is required for activity (50) (supplemental Fig. S3). It is interesting that the overall design of the cascade represents an excursion from Ser/Thr to Tyr chemistry. The cascade starts with a Ser/Thr MAP3K and ends as a Ser/Thr MAPK but requires Tyr kinase activity in the intermediate MAP2K.
We further observe both distributive (dissociative) and processive reactions. We observe processivity at the MAP2K level, with the appearance of the intermediate and product occurring simultaneously. The shapes of the progress curves for the appearance of MEK6/ST* and MEK6/S*T* are reminiscent of processive phosphorylation reactions induced by high viscosity reported by Aoki et al. (13) for ERK2. The phosphorylation of p38␣ is clearly distributive, with significant amounts of the intermediate forming before the final product. The kinetic modeling suggests that the last step to form p38a/T*Y* is slow, thus contributing to the buildup of intermediates. The presence of intermediates should contribute to the sigmoid properties of the cascade by slowing the formation of final product (43) and by creating competing binding interactions for the upstream enzyme (51).
Because of the differences in experimental approaches, the ability to compare rate constants with published cellular assays is limited. Nevertheless, the k cat values at the MAP2K level are comparable. We derived a k cat of 11 min Ϫ1 , which is in close agreement to the results of Fujioka et al. (15). On the other hand, they report significantly faster on rates than we observed at the MAP2K level. This could be due to our use of MAP3K kinase domains instead of full-length proteins and the lack of other cellular scaffolds.
Many multiple phosphorylation events in signal transduction are being discovered, and in a very few cases, the order of reactions has been shown to occur in a precise sequence (31,52). It is interesting that ATF2 is also phosphorylated in a specific order, and further, this is an integrator of signals from multiple MAP kinases (53,54). Additional analysis is required to determine whether the reaction order is maintained across multiple MAPKs acting on ATF2.
This study represents the first in vitro analysis of a complete MAP kinase cascade using native kinase substrates and their upstream activators. This approach, in conjunction with mass spectrometry, generates high fidelity data for modeling cascade reaction kinetics. It is clear from the reaction schemes in Figs. 1 and 4A that maximally switch-like responses should be available if intermediates are inactive. Future work will focus on comparisons of the activities of phosphorylation intermediates across different MAP kinase cascades and studies of reconstituted cascades in vitro.