The Complete Pathway for ERK2-catalyzed Reaction
Zhi-Xin Wang,Jia-Wei Wu
DOI: https://doi.org/10.1074/jbc.m703161200
IF: 5.485
2007-01-01
Journal of Biological Chemistry
Abstract:In the present study, the enzymatic mechanism of ERK2 is re-examined by a combination of steady-state kinetic studies in the absence and presence of viscosogenic agents. Kinetic studies carried out in various concentrations of sucrose revealed that both kcat and kcat/Km for either ATP or EtsΔ138 were highlysensitive to solvent viscosity, suggesting that the rapid equilibrium assumption is not valid for the phosphorylation of protein substrate by ERK2. Furthermore, the kinetic analysis with the minimal random Bi Bi reaction mechanism is shown to be inconsistent with the principle of the detailed balance. This inconsistent calculation strongly suggests that there is isomerization of the enzyme-substrate ternary complex. The viscosity-dependent steady-state kinetic data are combined to establish a kinetic mechanism for the ERK2-catalyzed reaction that predicts initial reaction velocities under varying concentrations of ATP and substrate. These results complement previous structure-function studies of mitogen-activated protein kinases and provide important insight for mechanistic interpretation of the kinase functions. In the present study, the enzymatic mechanism of ERK2 is re-examined by a combination of steady-state kinetic studies in the absence and presence of viscosogenic agents. Kinetic studies carried out in various concentrations of sucrose revealed that both kcat and kcat/Km for either ATP or EtsΔ138 were highlysensitive to solvent viscosity, suggesting that the rapid equilibrium assumption is not valid for the phosphorylation of protein substrate by ERK2. Furthermore, the kinetic analysis with the minimal random Bi Bi reaction mechanism is shown to be inconsistent with the principle of the detailed balance. This inconsistent calculation strongly suggests that there is isomerization of the enzyme-substrate ternary complex. The viscosity-dependent steady-state kinetic data are combined to establish a kinetic mechanism for the ERK2-catalyzed reaction that predicts initial reaction velocities under varying concentrations of ATP and substrate. These results complement previous structure-function studies of mitogen-activated protein kinases and provide important insight for mechanistic interpretation of the kinase functions. The mitogen-activated protein (MAP) 3The abbreviations used are: MAPmitogen-activated proteinERKextracellular signal-regulated kinaseJNKc-Jun N-terminal kinaseMEKMAP kinase/ERK kinaseMOPS4-morpholinepropanesulfonic acid. 3The abbreviations used are: MAPmitogen-activated proteinERKextracellular signal-regulated kinaseJNKc-Jun N-terminal kinaseMEKMAP kinase/ERK kinaseMOPS4-morpholinepropanesulfonic acid. kinases play a central role in signaling pathways in response to extracellular stimuli such as growth factors, cytokines, and various forms of environmental stress. The MAP kinase family comprises the ERK, JNK, and p38 subfamilies of kinases. All participate in three-tier protein phosphorylation cascades that serve to mediate diverse cellular processes. MAP kinase cascades are conserved in organisms ranging from yeast to human. A typical MAP kinase cascade consists of at least three kinases, a MAP kinase, a MAP kinase/ERK kinase (MEK) that activates the MAP kinase, and a MEK kinase that activates the MEK. Although ERKs respond robustly to growth factors and certain hormones, JNK kinases and p38 kinases react primarily to cytokines and stress stimuli. After activation, each MAP kinase phosphorylates a distinct spectrum of substrates, which include key regulatory enzymes, cytoskeletal proteins, nuclear receptors, regulators of apoptosis, and many transcription factors. Such an array of substrates is consistent with the observation that MAP kinases are involved in many critical cell functions (1Lewis T.S. Shapiro P.S. Ahn N.G. Adv. Cancer Res. 1998; 74: 49-139Crossref PubMed Google Scholar, 2Widmann C. Gibson S. Jarpe M.B. Johnson G.L. Physiol. Rev. 1999; 79: 143-180Crossref PubMed Scopus (2255) Google Scholar, 3Cobb M.H. Prog. Biophys. Mol. Biol. 1999; 71: 479-500Crossref PubMed Scopus (758) Google Scholar, 4Chen Z. Gibson T.B. Robinson F. Silvestro L. Pearson G. Xu B. Wright A. Vanderbilt C. Cobb M.H. Chem. Rev. 2001; 101: 2449-2476Crossref PubMed Scopus (779) Google Scholar). mitogen-activated protein extracellular signal-regulated kinase c-Jun N-terminal kinase MAP kinase/ERK kinase 4-morpholinepropanesulfonic acid. mitogen-activated protein extracellular signal-regulated kinase c-Jun N-terminal kinase MAP kinase/ERK kinase 4-morpholinepropanesulfonic acid. Because of their central role in signal transduction, the regulation and function of MAP kinases has been the subject of intense interest since their discovery. The prototypic member of the MAP kinase family is extracellular-regulated kinase 2 (ERK2). As a result, abundant molecular biological, biochemical, and structural information has occurred on ERK2. Like all MAP kinases, ERK2 catalyzes the transfer of the γ-phosphate of ATP to serine or threonine residues that precede proline (5Clark-Lewis I. Sanghera J.S. Pelech S.L. J. Biol. Chem. 1991; 266: 15180-15184Abstract Full Text PDF PubMed Google Scholar). Although this specificity is critical, further specificity is thought to be conferred through substrate docking domains where interactions are made between ERK2 and the substrate at sites on the enzyme that can lie a considerable distance from the active site (6Tanoue T. Adachi M. Moriguchi T. Nishida E. Nat. Cell Biol. 2000; 2: 110-116Crossref PubMed Scopus (671) Google Scholar, 7Seidel J.J. Graves B.J. Genes Dev. 2002; 16: 127-137Crossref PubMed Scopus (145) Google Scholar). Although progress in identifying these sequences and demonstrating their involvement in mediating specificity has advanced our knowledge considerably, precisely how they contribute to regulation and catalysis is still unclear. The mechanism of ERK2-catalyzed reaction has been elucidated by initial velocity and product inhibition studies (8Waas W.F. Dalby K.N. J. Biol. Chem. 2002; 277: 12532-12540Abstract Full Text Full Text PDF PubMed Scopus (39) Google Scholar). The substrate used in this study was EtsΔ138, a model protein substrate derived from the N terminus of Ets-1, which is currently the most thoroughly characterized protein substrate of ERK2. As shown by Seidel and Graves (7Seidel J.J. Graves B.J. Genes Dev. 2002; 16: 127-137Crossref PubMed Scopus (145) Google Scholar), it is an excellent surrogate for the full-length protein. In the presence of MgATP2– and excess Mg2+, ERK2 phosphorylates EtsΔ138 exclusively on Thr-38 to a stoichiometry of 1 mol/mol. The kinetic dependence on substrate concentration and product inhibition patterns was consistent with a rapid equilibrium random order ternary complex mechanism. Although product inhibition data are not definitive evidence for a rapid equilibrium mechanism, the approximation of steady-state data to a rapid equilibrium model suggests that the dissociation rate constants from the ternary complex for both substrates must be of greater magnitude than the catalytic constant. On the other hand, Prowse et al. (9Prowse C.N. Hagopian J.C. Cobb M.H. Ahn N.G. Lew J. Biochemistry. 2000; 39: 6258-6266Crossref PubMed Scopus (51) Google Scholar) employed solvent viscosometric techniques to characterize the catalytic reaction pathway of ERK2 with respect to the phosphorylation of myelin basic protein. In this study the investigators showed that the catalytic mechanism for this reaction consists of the diffusion-limited binding of myelin basic protein (k–2 << k3) and the rapid equilibrium binding of ATP k−2'≫k3, followed by phosphoryl group transfer, which is slow in comparison with the release of products. This is in contrast to the kinetic mechanism established for the phosphorylation of EtsΔ138 by ERK2, in which the dissociation rate constants of EtsΔ138, k–2, is of similar magnitude to the catalytic constant, k3. Furthermore, as ATP binds deep within the active site cleft of protein kinases (10Akamine Madhusudan Nat. Struct. Biol. 2002; 9: 273-277Crossref PubMed Scopus (188) Google Scholar), it is conceivable that such binding is prevented when protein substrate is bound. Because the rapid equilibrium ordered mechanism with EB and EP abortive complexes is distinguishable by initial velocity and product inhibition studies from the random addition case, further studies are required to address this point (11Frieden C. Biochem. Biophys. Res. Commun. 1976; 68: 914-917Crossref PubMed Scopus (29) Google Scholar). In the present study, the enzymatic mechanism of ERK2 is re-examined by a combination of steady-state kinetic studies in the absence and presence of viscosogenic agents. Our results showed that 1) the reaction pathway for the turnover of transcription factor EtsΔ138 occurs through two partially rate-limiting steps: one step correlates with phosphorylation of EtsΔ138 on the enzyme, and the second step is related to product release; 2) the dissociation rate constant of EtsΔ138 is much less than the catalytic constant, and therefore the rapid equilibrium assumption is not valid for the phosphorylation of EtsΔ138 by ERK2; and 3) the kinetic analysis with the minimal steady-state random mechanism was shown to be inconsistent with the principle of the detailed balance, suggesting that there is isomerization of the enzyme-substrate ternary complex. These results complement previous structure-function studies of MAP kinases and provide important insight for mechanistic interpretation of the kinase functions. Materials-ATP, phospho(enol) pyruvate, nicotinamide adenine dinucleotide, reduced NADH, lactate dehydrogenase, and pyruvate kinase, were purchased from Sigma. Sucrose was purchased from Invitrogen. MOPS was purchased from Roche Applied Science. Other reagents were local products of analytical grade used without further purification. Double-deionized water was used throughout. Expression and Purification of Proteins-The plasmid pETHis6-ERK2-MEK1(R4F) (a generous gift of Dr. Melanie Cobb) was used to coexpress a constitutively active MEK1 and an N-terminal His6-tagged ERK2 in Escherichia coli BL21(DE3). The expression and purification of bisphosphorylated ERK2 were carried out following the procedure described by Wilsbacher and Cobb (12Wilsbacher J.L. Cobb M.H. Methods Enzymol. 2001; 332: 387-400Crossref PubMed Scopus (23) Google Scholar). After the final fast protein liquid chromatography step, ∼3 mg of ERK2/pTpY were obtained from 6 liters of culture. Electrospray ionization mass spectrometry analysis confirmed that the purified ERK2/pTpY was effectively homogenous and was phosphorylated stoichiometrically to a ratio of 2 mol of phosphate/mole of ERK2. Plasmid used to express Ets-1 (1–138) has been reported previously. The N-terminal His6-tagged Ets-1 (EtsΔ138) was expressed in E. coli BL21(DE3) and purified using standard procedures of nickel chelate affinity chromatography (13Waas W.F. Dalby K.N. Protein Expression Purif. 2001; 23: 191-197Crossref PubMed Scopus (21) Google Scholar). The protein purity was over 90% as judged by SDS-PAGE. The purified proteins were made to 20% glycerol and stored at –80 °C. Enzyme Assays for ERK2-Enzyme activity of ERK2 was determined with the EtsΔ138 as a substrate using a coupled spectrophotometric assay (14Cook P.F. Neville M.E. Biochemistry. 1982; 21: 5794-5799Crossref PubMed Scopus (346) Google Scholar). The standard assay for ERK2 was carried out at 30 °C in 1.8 ml of reaction mixture containing 100 mm MOPS buffer, pH 7.0, 200 mm NaCl, 20 mm MgCl2, 200 μm NADH, 1 mm phospho(enol) pyruvate, 18 units/ml lactate dehydrogenase, 7.5 units/ml pyruvate kinase, and different concentrations of substrates. The reactions were initiated by the addition of ERK2 to the reaction mixture. Progress of the reaction was monitored continuously by following the formation of NAD+ at 340 nm, on a PerkinElmer Life Sciences Lambda 45 spectrophotometer equipped with a magnetic stirrer in the cuvette holder. The concentrations of ADP formed in ERK2-catalyzed reaction were determined using an extinction coefficient for NADH of 6220 cm–1 m–1 at 340 nm. The concentration of ERK2 was determined spectrophotometrically based on an extinction coefficient (ϵ = 44,825 cm–1m–1) calculated from its primary amino acid sequence and verified by the method of Bradford with bovine serum albumin as a standard. The concentration of EtsΔ138 was determined by turnover with the enzyme under conditions of limiting EtsΔ138 at a fixed low concentration of ATP (17.5 μm). Solvent Viscosometric Studies-Steady-state assays were carried out as described above in buffer containing varied sucrose, corresponding to relative solvent viscosities varying between 1 and 2.92. The relative solvent viscosities (ηrel) were determined from the equation ηrel=(t/t0)(ρ/ρ0) (Eq. 1) where ηrel is the solvent viscosity relative to the buffer containing no viscosogen, t is the solvent transit time measured by an Ostwald capillary viscometer, and ρ is the solvent specific gravity (15Freifelder D. Physical Biochemistry: Applications to Biochemistry and Molecular Biology. 1982; Google Scholar). Superscript 0 indicates conditions of the buffer with no viscosogen. Effects of Viscosogens on the Steady-state Kinetic Parameters - Scheme 1 describes the minimal kinetic mechanism for an enzyme-catalyzed two-substrate, two-product reaction, where k3 is the rate constant for the functional group transfer step, k4 the rate constant for products release, E represents enzyme, and A and B represent ATP and the protein substrate, respectively. When the concentration of one substrate is infinite, the enzyme can be treated as “one-substrate” enzyme for kinetic purposes, and the steady-state rate equation at saturating A can be written as follows. v=c2[B][E]0c8+c10[B]=kcatB[B][E]0KmB+[B] (Eq. 2) where kcatB=k3k4k3+k4 (Eq. 3) KmB=k4k3+k4k−2+k3k2 (Eq. 4) Combining Equations 3 and 4, we obtain kcatBKmB=k2k3k−2+k3 (Eq. 5) Because k2, k–2, and k4 are presumed to represent diffusion-controlled processes, these rate constants will be dependent on the relative viscosity of the solvent. A comprehensive explanation of the use of solvent viscosogens in enzymatic analysis is presented by Gutfreund (16Gutfreund H. Kinetic for the Life Sciences. 1995; : 248-261Google Scholar), and its application to protein kinases is reviewed by Adams (17Adams J.A. Chem. Rev. 2001; 101: 2271-2290Crossref PubMed Scopus (481) Google Scholar). The relative solution viscosity (ηrel) can be equated to the ratio of the microscopic rate constants, ηrel=k20/k2=k−20/k−2=k40/k4(Eq. 6) where k20 is the association rate constant of substrate, and k−20 and k40 are the dissociation rate constants of substrate and product, respectively, when ηrel = 1. By taking separately the ratio of Equations 3 and 5 in the absence and presence of viscosogen according to Equation 6, two linear equations are derived upon plotting kcatB0/kcatB and kcat/KmB0/kcat/KmB as a function of ηrel. The slopes of these lines are given by KcatBη=K3K3+K40(Eq. 7) and KcatBKmBη=K3K−20+K3(Eq. 8) Thus, (k(kcatB)η and (kcatB/KmB)η can vary between 0 and 1, the limit of diffusion control. Similarly at high [B] v=C3AE0C9+C11A=KcatAAE0KmA+A(Eq. 9) where KcatA=K3K4K3+K4(Eq. 10) and KmA=K4K3+K4K−2'+K3K2'(Eq. 11) Saturating ATP and Varying EtsΔ138 Concentrations-To characterize the steady-state kinetic properties of ERK2, the kinase activity of ERK2 were measured by following ADP formation using an enzyme-coupled spectrophotometric assay. All of the initial velocity studies were carried out using a continuous assay in which the linearity of the initial portion of the reaction profile was carefully monitored. Because the spontaneous ATP hydrolysis by ERK2 occurs in parallel to the phosphorylation reaction, the assays were corrected for the ATPase effect, which is taken as the observed absorption changes in the absence of EtsΔ138. We first attempted to define the diffusion-controlled steps in the catalytic reaction of ERK2 by examining the steady-state phosphorylation of EtsΔ138 in the presence of sucrose. The initial velocity data were collected under conditions of saturating ATP (2.5 mm), varied EtsΔ138, and varied relative solvent viscosities ranging from 1 to 2.92. When the initial velocity data were plotted versus EtsΔ138 concentrations, typical Michaelis-Menten behavior was observed in the absence and presence of the viscosogen. Fig. 1 shows the effect of increasing solvent viscosity on the steady-state kinetic parameters for the phosphorylation of EtsΔ138 by ERK2. It can be seen from this figure that both kcatB and kcatB/KmB are linear functions of the relative solvent viscosity. The slope of a plot of KcatB0 versus ηrel,(kcatB)η, which indicates the sucrose effect on EtsΔ138 phosphorylation, was determined to be 0.49 ± 0.02. From these values, the rates of product release and phosphoryl group transfer can then be calculated as K40=KcatB0/kcatBη=37.76s−1 and K3=KcatB0/1−kcatBη=36.27s−1, respectively. The dissociation of products, ADP and phospho-EtsΔ138, is combined in K40, and the viscosity method cannot discern which product release step controls K40 without other kinetic experiments. By comparison, a nearly full viscosity effect was observed on kcatB/KmB. The slope of a plot of (kcat/KmB)0/(kcat/KmB) versus ηrel,(kcatB/KmB)η, was determined to be 0.92 ± 0.05. Because the value of k3 is known, the dissociation rate constant of EtsΔ138 to the E-ATP complex can then be calculated as K−20=K31−KcatB/KmB)η/KcatB/KmBη=3.15s−1. The corresponding association rate constant can then be calculated as K20=KcatB/KmB0/KcatB/KmBη=1.90s−1. Saturating EtsΔ138 and Varying ATP Concentrations-We next tested the effect of increased solvent viscosity on the catalytic reaction mechanism when the concentration of ATP was varied and the amount of EtsΔ138 was fixed. The effect of increasing solvent viscosity on the steady-state kinetic parameters for the phosphorylation of EtsΔ138 by ERK2 is shown in Fig. 2. At saturating EtsΔ138 concentration (440 μm), the viscosity effect on kcat was close ((kcat)η = 0.42 ± 0.01) to that observed when the EtsΔ138 concentration was varied under saturating ATP. Similarly, a partial viscosity effect on kcatA/KmA at high EtsΔ138 concentration was observed ((kcatA/KmA)η = 0.69 ± 0.01). If an ordered mechanism (with ATP binding first) is assumed, the value for kcatA/KmA in such a mechanism would be equal to k1, and the viscosity effect on the parameter would be predicted as (kcatA/KmA)η = 1 at saturating concentrations of EtsΔ138. Our viscosity data, apparently, conflict with such an ordered binding mechanism. Taken together, the results from viscosity experiments are consistent with random addition of ATP and EtsΔ138 to the ERK2 active site. In this model, the substrate ATP combines with the enzyme-EtsΔ138 complex to form the active ternary complex by the second order association rate constant, k2', and the dissociation rate constant, k−2'. Therefore, the association and dissociation rate constants of ATP from the enzyme-EtsΔ138 complex can be estimated as k2'o=KcatA/KmA0/KcatA/KmAη=0.10μM−1S−1 and K−2'o=k3[1−(kcatA/KmA)η]/(kcatA/KmA)η=15.54 s−1, respectively. Note that the full or partial viscosity effects on kcat/Km indicate that the dissociation rates of ATP and EtsΔ138 from the E-ATP-EtsΔ138 complex are not much faster than the phosphoryl group transfer of ATP to the EtsΔ138. Thus, the ERK2-catalyzed reaction has a random, but not truly rapid equilibrium random mechanism. Steady-state Kinetics of the ERK2-catalyzed EtsΔ138 Phosphorylation Reaction-To analyze the kinetic mechanism of ERK2 catalyzed reaction further, 12 data sets of initial velocities were obtained over a wide range of ATP and EtsΔ138 concentrations. Fig. 3 shows double reciprocal plots of the initial velocity versus either EtsΔ138 or ATP concentrations (for clarity, only six data sets are shown). It can be seen from Fig. 3 that the Lineweaver-Burk plots were linear and displayed a nonparallel pattern of intersecting lines. Thus, a ping-pong mechanism (where a parallel line pattern is characteristic) can be excluded from such analysis. As mentioned above, the simplest explanation for the kinetic observation is the random Bi Bi mechanism. It is therefore reasonable to start the analysis of our experimental data with this mechanism first. The exact steady-state rate equation for the general sequential mechanism is too complicated for practical use. As pointed by Dalziel (18Dalziel K. Trans. Faraday Soc. 1958; 545: 1247-1253Crossref Scopus (30) Google Scholar), however, if substrate concentrations are sufficiently small to satisfy conditions of k–1 >> k [B] and k−1'≫k2'A, the intial rate equation for the steady-state random mechanism reduces to the simple form v=K3K4K3+K4E0ABK4K−2+K−2'+K3K−1K−1'+K1K−1'A+K−1K1'BK3+K4K2'K−1K1'+K2K1K−1'+AB(Eq. 12) which is identical to the rate equation obtained by Seubert et al. with the partial equilibrium assumption (19Seubert P.A. Renosto F. Knudson P. Segel I.H. Arch. Biochem. Biophys. 1985; 240: 509-523Crossref PubMed Scopus (36) Google Scholar). Therefore, the condition of k–1 >> k2[B] and K−1'≫K2'A is equivalent to the assumption that the binary enzyme-substrate complexes remain in equilibrium with the free enzyme and substrate. Because the values of k3, K40,K20,K−20,K2'0, and K−2'0 have been determined separately by solvent viscosometric analysis, we can fix them to constant values in curve fitting. The data sets of initial velocities were analyzed with Equation 12 by global fitting procedures, and the kinetic parameters were determined to be KA=K−10/K10=32.26±4.22 μm, and KB=K−10/K10=3.57±0.54μM, respectively. However, according to the principle of detailed balance, the rate constants must satisfy the following condition (20Hammes G.G. Thermodynamics and Kinetic for the Biological Sciences. 1995; Google Scholar). K1K2K−2'K−1'=K1'K2'K−2K−1(Eq. 13) or K2K−2'KB=K2'K−2KA(Eq. 14) It is important to confirm that this principle is obeyed whenever reaction cycles are present. Obviously, our kinetic analysis is inconsistent with the principles of detailed balance, and thus, the calculation of rate constants is invalid. Thus, the minimal random mechanism is not applicable to our experimental data. These inconsistent calculations strongly suggest that there is isomerization of the ERK2-EtsΔ138 complex. Kinetic Analysis with an Iso Random Bi Bi Mechanism-The steady-state kinetic data and the effects of solvent viscosity on the phosphorylation of EtsΔ138 can be combined to outline a complete kinetic mechanism for the ERK2-catalyzed reaction. Scheme 2 describes the sequential step processing of ATP and EtsΔ138 under all substrate concentrations. The steady-state rate equation at saturating A is given by the following. KcatB=K3K4KcK3Kc+K4Kc+K−c(Eq. 15) KmB=K4KcK−2+K3K2K3Kc+K4Kc+K−c(Eq. 16) and KcatBKmB=K2K3K−2+K3(Eq. 17) Applying the same viscosity constraints as for Equations 15 and 17, the slopes of the viscosity-dependent steady-state kinetic parameters are KcatBη=K3KcK3Kc+K40Kc+K−c(Eq. 18) and kcatBKmBη=k3k−20+k3 (Eq. 19) Because (kcatB)η = 0.49 ± 0.02 and (kcat/KmB)η = 0.92 ± 0.05, the rate constants of product release and substrate association can be calculated as k40=kcatB0/(kcatB)η=37.76±3.83 s−1 and k20=(kcatB0/KmB0)/(kcat/KmB)η=1.90±0.56 s−1, respectively. When substrate B is saturating, we have kcatA=k3k4kck3kc+k4(kc+k−c+k3) (Eq. 20) KmA=K4K−2'K−c+K3+K3KcK2'K3Kc+K4Kc+K−c+K3(Eq. 21) kcatAKmA=K2'K3KcK−2'K−c+K3+K3Kc(Eq. 22) The slopes of the viscosity-dependent steady-state kinetic parameters are given by the following. (kcatA)η=k3kck3kc+k40(kc+k−c+k3)(Eq. 23) kcatAKmAη=k3kcK−2'0K−c+K3+K3Kc(Eq. 24) Similarly, from the value of (kcatA/KmA)η = 0.69 ± 0.01, the association rate constant of ATP can be estimated to be (kcatA/KmA)0/(kcatA/KmA)η = 0.10 ± 0.01 μ–1 s–1. In addition, from Equations 18 and 23, we have the following. 1(kcatA)η−1(kcatB)η=k40kc (Eq. 25) Hence kc=k40(kcatB)η(kcatA)η(kcatB)η−(kcatA)η=111.01±50.43 (Eq. 26) In this case, however, the values of k3, k−2'0, and k−20 cannot be determined directly from the viscosity experiment. Instead, the following relationships between the rate constants can be obtained. k3=kcatB01−(kcatB)ηkc+k−ckc=36.27kc+k−ckc (Eq. 27) k−20=k3{1−(kcatB/KmB)η}(kcatB/KmB)η=0.087k3 (Eq. 28) k−2′0={1−(kcatA/KmA)η}(kcatA/KmA)ηk3kck−c+k3=0.43k3kck−c+k3 (Eq. 29) To determine the values of other kinetic parameters in Scheme 2, the data sets of initial velocities shown in Fig. 3 must be reanalyzed according to the rate equation derived from this model. The minimum assumption that needs to be made to obtain a velocity equation without squared terms is that the reactions forming the EA and EB complexes are at equilibrium. The method of Cha (21Cha S. J. Biol. Chem. 1968; 243: 820-825Abstract Full Text PDF PubMed Google Scholar) can then be used to solve for the concentrations of various enzyme species. The resulting expression for the initial rate of reaction as a function of the concentrations of substrates is given by the following. ν=kc[E]0[A][B]α{KAKB+KB[A]+KA[B]}+β[A][B] (Eq. 30) where kc=k3k4k2KBkc+k−2'+k−2'KAKc(Eq. 31) α=k4[k−2'(k−2+k−c+k3)+kc(k−2+k3)] (Eq. 32) β=k2'KA[k4(k−2+k−c+k3)+kc(k3+k4)]+k2KB[k4k−c+(kc+k−2')(k3+k4)] (Eq. 33) The required microscopic reversibility condition for Scheme 2 is then given by the following. KB=k2'k−2kcKAk2k−2'k−c(Eq. 34) Therefore, by substituting Equations 27, 28, 29 and 34 into Equation 30 and treating the parameters, k40,k20,k2'0, and kc as fixed constants, all of the data (12 data sets) were simultaneously analyzed with Equation 30. The resulting family of curves shown in Fig. 3 represents a best fit of the model to all data points simultaneously constrained, from which the kinetic parameters were determined to be KA = 107.21 ± 0.87 μm, and k–c = 3.93 ± 0.07 s–1, respectively. It can be seen from this figure that the experimental data fit very well to this model with standard errors less than 2%, indicating that the proposed model (Scheme 2) satisfactorily described the reaction. Thus, the values of k3, k−20, k−2'0, and KB can be calculated from Equations 27, 28, 29 and 34. All of the kinetic parameters for the mechanism of Scheme 2 are listed in Table 1.TABLE 1Kinetic constants for phosphorylation of EtsΔ 138 by ERK2Kinetic parameterValueKATP107.21 μmKEtsΔ13812.05 μmK201.90 μm-1 s-1K−203.27 s-1K2'00.10 μm-1 s-1K−2'043.21 s-1kc111.01 s-1k-c3.93 s-1k337.55 s-1K4037.76 s-1 Open table in a new tab It should be indicated that the isomerization step invoked in the reaction pathway of binding the peptide substrate first could be either before (in the binary complex) or after (in the ternary complex) ATP binding. If an intermediate is formed before the central ternary complex, the initial rate equation is first degree with respect to [B] and secondary degree to [A]. It predict that the plots of 1/v versus 1/[B] are linear and plots of 1/v versus 1/[A] are theoretically nonlinear. These random Bi Bi mechanisms can be eliminated readily, because no evidence for this type behavior has been observed. The sequential random two-substrate mechanism is of great importance in enzyme kinetics and has been discussed by several authors (18Dalziel K. Trans. Faraday Soc. 1958; 545: 1247-1253Crossref Scopus (30) Google Scholar, 22Ferdinand W. Biochem. J. 1966; 98: 278-283Crossref PubMed Scopus (130) Google Scholar, 23Pettersson G. Acta Chem. Scand. 1969; 23: 2717-2726Crossref Google Scholar). A complete steady-state treatment for the random ternary complex mechanism yields a complicated velocity equation containing [A]2, [B]2, [A]2[B], and [A][B]2 terms. It has, in fact, often been stated in the literature that random order mechanisms (with exception for the rapid equilibrium case) always yield nonlinear Lineweaver-Burk plots, and a large number of two-substrate enzyme systems have been claimed to operate by a compulsory order or a rapid equilibrium mechanism merely because they obey a hyperbolic type of rate equation. Some years ago, Pettersson (23Pettersson G. Acta Chem. Scand. 1969; 23: 2717-2726Crossref Google Scholar) carried out detailed theoretical studies on random order steady-state systems and showed that the degree to which reciprocal plots deviate from linearity in systems that obey a steady-state random order mechanism will depend on the values of individual rate constants relative to the concentrations of the two substrates. At relatively low substrate concentrations, linear Lineweaver-Burk relationships will always be obtained. Gulbinsky and Cleland (24Gulbinsky J. Cleland W.W. Biochemistry. 1968; 2: 566-575Crossref Scopus (63) Google Scholar) have also shown by computer simulation studies of the random order two-substrate mechanism that the initial velocity patterns did not deviate significantly from linearity under some circumstances in which equilibrium conditions did not apply. The present investigation confirmed experimentally that although the dissociation rate for EtsΔ138 (k–2) is 10-fold slower than that for the phosphoryl group transfer step (k3) and therefore the rapid equilibrium assumption is not valid, linear Lineweaver-Burk relationships were indeed observed when the kinetic data of ERK2-catalyzed reaction were plotted in the traditional double-reciprocal forms. A thorough understanding of the mechanism of protein phosphorylation catalyzed by protein kinases requires the unambiguous identification of individual kinetic steps comprising the catalytic pathway. In this report we have deduced from the studies of initial velocity and viscosity that ERK2 conforms to an iso random mechanism for the phosphorylation of EtsΔ138. All of the rate constant values to each step in the overall kinetic mechanism for the phosphorylation of EtsΔ138 have been determined from steady-state kinetic data. The concentration of ATP in mammalian cells is ∼1 mm, which is 10-fold higher than the dissociation constant for ATP and the free enzyme, KA (25Gribble F.M. Loussouarn G. Tucker S.J. Zhao C. Nichols C.G. J. Biol. Chem. 2000; 275: 30046-30049Abstract Full Text Full Text PDF PubMed Scopus (237) Google Scholar). Thus, although the mechanism is formally random sequential, there is a preferred binding order under limiting amounts of protein substrate with ATP binding first, because 1 mm ATP is capable of trapping 90% of the free enzyme. However, at higher protein substrate concentrations (say 10-fold higher than KB), binding of the protein substrate as the first substrate is physiological relevant, and ATP binds to form an unproductive enzyme-substrate conformer that can reversibly covert to productive ternary complex. Currently, the molecular details of ERK2-substrate interactions remain unclear, although there is a large body of evidence suggesting that ERK2 facilitates substrate recognition, at least in part, through interactions that are discrete from those at the active site (6Tanoue T. Adachi M. Moriguchi T. Nishida E. Nat. Cell Biol. 2000; 2: 110-116Crossref PubMed Scopus (671) Google Scholar). Because the three-dimensional structure of ERK2 cocrystallized with a protein substrate or substrate analogue has not been solved, there is little mechanistic information available to ultimately correlate the structural data with the substrate recognition properties of this enzyme. Some years ago, Adams and Taylor (26Adams J.A. Taylor S.S. Biochemistry. 1992; 31: 8516-8522Crossref PubMed Scopus (159) Google Scholar) proposed a similar kinetic mechanism for the catalytic subunit of cAMP-dependent protein kinase. The x-ray crystal structure of the binary enzyme-peptide complex of the catalytic subunit of cAPK indicates that ATP binds deep within the active site groove with peptide directly lining the periphery. Therefore, they concluded that ATP prebinding directs the productive association of substrate, and in the presence of a large substrate pool, bound peptide substrate adjusts its conformation upon ATP binding to position the γ phosphate for delivery. It is likely that ERK2 has properties similar to those of the catalytic subunit of cAPK.