Previous Article | Next Article ![]()
Journal of Bacteriology, April 2005, p. 2386-2394, Vol. 187, No. 7
0021-9193/05/$08.00+0 doi:10.1128/JB.187.7.2386-2394.2005
Copyright © 2005, American Society for Microbiology. All Rights Reserved.
,
Andrea Gorrell,
,
Sarah H. Lawrence,
Prabha Iyer,¶
Kerry Smith,
and
James G. Ferry*
Department of Biochemistry and Molecular Biology, Pennsylvania State University, University Park, Pennsylvania
Received 22 October 2004/ Accepted 8 December 2004
|
|
|---|
-phosphoryl group to acetate. Inspection of the crystal structure of the Methanosarcina thermophila enzyme containing only ADP revealed a solvent-accessible hydrophobic pocket formed by residues Val93, Leu122, Phe179, and Pro232 in the active site cleft, which identified a potential acetate binding site. The hypothesis that this was a binding site was further supported by alignment of all acetate kinase sequences available from databases, which showed strict conservation of all four residues, and the recent crystal structure of the M. thermophila enzyme with acetate bound in this pocket. Replacement of each residue in the pocket produced variants with Km values for acetate that were 7- to 26-fold greater than that of the wild type, and perturbations of this binding pocket also altered the specificity for longer-chain carboxylic acids and acetyl phosphate. The kinetic analyses of variants combined with structural modeling indicated that the pocket has roles in binding the methyl group of acetate, influencing substrate specificity, and orienting the carboxyl group. The kinetic analyses also indicated that binding of acetyl phosphate is more dependent on interactions of the phosphate group with an unidentified residue than on interactions between the methyl group and the hydrophobic pocket. The analyses also indicated that Phe179 is essential for catalysis, possibly for domain closure. Alignments of acetate kinase, propionate kinase, and butyrate kinase sequences obtained from databases suggested that these enzymes have similar catalytic mechanisms and carboxylic acid substrate binding sites. |
|
|---|
![]() | (1) |
![]() | (2) |
The reaction catalyzed by acetate kinase was initially described in 1944 by Lipmann (24). Following the first purification in 1954 from Escherichia coli (31), the enzyme was the subject of several investigations, yet questions concerning the catalytic mechanism and particularly the substrate binding sites have remained unanswered. Investigations with the E. coli enzyme have led to two proposals for the catalytic mechanism: direct in-line transfer of the
-phosphoryl group of ATP to acetate (6, 34) and a covalent triple-displacement mechanism involving two phosphoenzyme intermediates (35). Although the phosphorylated E. coli acetate kinase is able to transfer the phosphoryl group to acetate (39, 40), it is not kinetically competent (6). Furthermore, the phosphoenzyme has been shown to phosphorylate enzyme I of the bacterial phosphotransferase system (13) and CheY (11), a member of the flagellar motor cascade; thus, the phosphoenzyme may function only in sugar transport and chemotaxis rather than play an essential role in the catalytic mechanism of equation 1. Finally, the acetate kinase of Methanosarcina thermophila has been shown to be inhibited by a putative transition state analogue, ADP-AlF3-acetate, in which the AlF3 is proposed to mimic the meta-phosphate in a direct phosphoryl transfer mechanism (26). The E. coli and M. thermophila acetate kinases exhibit 44% sequence identity, suggesting that they utilize similar catalytic mechanisms (32).
The first crystal structure reported for acetate kinase was that of the M. thermophila enzyme containing ADP bound in the active site cleft proximal to a solvent-accessible pocket of hydrophobic residues proposed to accept the methyl group of acetate or acetyl phosphate (8). A recent crystal structure of the M. thermophila acetate kinase containing acetate, ADP, and the transition state analog AlF3 further supported the hypothesis that there is a direct in-line mechanism and showed the methyl group of acetate located in the previously proposed binding pocket (14). Kinetic analyses of active site replacement variants also supported the hypothesis that there is a direct in-line mechanism for the M. thermophila enzyme and suggested roles for conserved arginine, histidine, and glutamate residues (16, 26, 27, 32, 33). Here we report kinetic analyses of variants designed to characterize the proposed hydrophobic binding pocket which indicated that the pocket has a role in both acetate binding and specificity. Finally, our results also suggested that the conserved Phe179 residue has a role in catalysis by contributing to domain movement.
|
|
|---|
Analysis of acetate, propionate, and butyrate kinase sequences. The nonredundant protein and nucleotide sequence databases at the National Center for Biotechnology Information were searched for acetate, propionate, and butyrate kinase sequences by using the BLAST network server and the BLASTp and tBLASTn programs (3, 4). Sequences were aligned with ClustalX (38) by using a Gonnet PAM 250 weight matrix with a gap opening penalty of 10.0 and a gap extension penalty of 0.05.
Site-directed mutagenesis. Mutagenesis was performed by the oligonucleotide-directed in vitro mutagenesis method (21) with a QuikChange mutagenesis kit (Stratagene). Plasmid pML703 (23), a derivative of the expression vector pT7-7 (37) containing the M. thermophila ack gene, was the target for mutagenesis with the primers listed in Table 1. Mutations were verified by dye termination cycle sequencing by using an ABI PRISM 377 DNA sequencer (Applied Biosystems) at the Nucleic Acid Facility at Pennsylvania State University.
|
View this table: [in a new window] |
TABLE 1. Mutagenic oligonucleotide primers
|
(DE3)] and were purified as described previously (23). Protein purity was examined by sodium dodecyl sulfate (SDS)-polyacrylamide gel electrophoresis (PAGE) (22), and protein concentrations were determined by the Bradford method (7) by using the Bio-Rad dye reagent with bovine serum albumin as the standard.
Molecular masses.
The native molecular masses of the wild-type and variant enzymes were determined by gel filtration chromatography by using a Superose 12 gel filtration column (Amersham Pharmacia Biotech) calibrated with blue dextran (2,000 kDa), urease (trimer, 272 kDa; hexamer, 545 kDa), bovine serum albumin (monomer, 66 kDa; dimer, 132 kDa), chicken egg albumin (45 kDa), bovine erythrocyte carbonic anhydrase (29 kDa), and bovine milk
-lactalbumin (14.2 kDa). The column was preequilibrated with 50 mM potassium phosphate (pH 6.8) containing 150 mM KCl and was developed at a flow rate of 0.4 ml/min.
Kinetic parameters of wild-type and variant acetate kinases. The hydroxamate assay, an adaptation of the methods of Lipmann and Rose et al. (1, 24, 31), detects acetyl phosphate formation from acetate and ATP and was used to measure the initial acetate kinase activities and kinetic parameters in the forward (ADP/acetyl phosphate-producing) direction. By utilizing the pyruvate kinase (PK)-lactate dehydrogenase (LDH) coupled assay system (1), hydroxylamine was found to be an inhibitor of acetate kinase activity (14); therefore, kinetic parameters of wild-type and variant acetate kinases were determined by utilizing the PK-LDH coupled assay system. Briefly, each assay solution contained 60 mM HEPES (pH 7.0), 5 mM MgCl2, 16.7 U of PK, 36 U of LDH, 3 mM phosphoenolpyruvate, and 0.2 mM NADH along with a fixed concentration of substrate (200 mM acetate or 1 to 2 mM equimolar ATP-MgCl2) as appropriate. The wild-type or variant enzyme concentrations ranged from 0.5 to 50 µg/ml, depending upon the specific activity. Changes in absorbance at 340 nm were monitored for 1 to 5 min at 1-s intervals with a Beckman DU640 spectrophotometer. The Km and kcat values were determined by nonlinear regression data analysis fit to the Michaelis-Menten equation by using the Kaleidagraph program (Synergy Software, Reading, Pa.). The kinetic parameters reported below are averages of at least three independent trials, with kcat values independent of the variable substrate.
To determine kinetic parameters of the ATP-producing reaction, the previously described hexokinase-glucose 6-phosphatase enzyme-linked assay was used (1). Briefly, the assay mixtures contained 100 mM Tris (pH 7.4), 0.2 mM dithiothreitol, 10 mM MgCl2, 4.4 mM glucose, 1 mM NADP, 10 U of hexokinase (yeast), 10 U of glucose-6-phosphate dehydrogenase (yeast), 5 mM ADP, and different acetyl phosphate concentrations. Kinetic constants were determined by using nonlinear regression to fit data with the Kaleidagraph program (Synergy Software).
Kinetic parameters with the alternative substrates propionate and butyrate. The abilities of the wild-type and variant acetate kinases to utilize the substrates propionate and butyrate in the propionyl or butyryl phosphate reaction direction were determined by utilizing the PK-LDH coupled assay system described above, with the ATP-MgCl2 concentration kept at 2 mM and various concentrations of propionate or butyrate as appropriate for the enzyme of interest.
Cavity volume determination and modeling of propionate and butyrate in the active site pocket. Models of Val93Ala and Val93Gly variant enzymes were generated by utilizing O (17), and GROMOS96 energy minimization was performed (41). The volumes of the cavities in the wild-type, Val93Ala, and Val93Gly enzymes were determined by utilizing VOIDOO (19) with a 1.4-Å rolling sphere probe and a starting point of 32.0, 44.0, 55.8 (x, y, z). PDB, topology, parameter, and connectivity files for acetate, propionate, and butyrate were obtained from the HICUP website (18). Acetate, propionate, and butyrate were positioned in the cavities located in the wild-type, Val93Ala, and Val93Gly active sites, respectively, and GROMOS96 energy minimization was performed for each preparation (41).
|
|
|---|
![]() View larger version (17K): [in a new window] |
FIG. 1. Stereo view of the proposed acetate binding site in the acetate kinase from M. thermophila. The figure was generated with previously published coordinates (9) by using MOLSCRIPT (20) and Raster3D (25). All residues are labeled at the corresponding carbons.
|
![]() View larger version (58K): [in a new window] |
FIG. 2. Alignment of acetate kinase, propionate kinase, and butyrate kinase sequences. The deduced amino acid sequences of the enzymes were aligned by using ClustalX (38). Residues essential for acetate kinase activity (number signs) and for the hydrophobic pocket (percent signs) are indicated above the alignment. Identical (asterisks) and similar (colons or periods) residues in the sequences are indicated below the alignment. The values in parentheses are the levels of identity followed by the levels of similarity (expressed as percentages) to the acetate kinase from M. thermophila. Abbreviations: MT AK, M. thermophila acetate kinase (accession number gi:584720); LC PK, Lactobacillus collinoides propionate kinase (gi:29335735); SE PK, S. enterica subsp. enterica serovar Typhimurium propionate kinase (gi:5069465); EC PK, E. coli propionate kinase (gi:1176151); CB BK1, Clostridium acetobutylicum butyrate kinase I (gi:20137334); CP BK, Clostridium perfringens butyrate kinase (gi:4239872); CB BK2, C. acetobutylicum butyrate kinase II (gi:20137415); TT BK, Thermoanaerobacter tengcongensis butyrate kinase (gi:20517209); EF BK, Enterococcus faecalis butyrate kinase (gi:20137247).
|
Kinetic parameters of wild-type and variant acetate kinases. The kinetic constants for the wild-type and variant enzymes assayed in the direction of acetyl phosphate synthesis are shown in Table 2. All variants showed less than twofold decreases in the Km for ATP relative to the wild-type enzyme, further supporting the hypothesis that the overall structure of the active site was not affected in the variant enzymes. However, all of the variants showed increases in the Km for acetate relative to the wild type; these increases ranged from 7-fold for the Phe179Ala variant to 26-fold for the Val93Gly variant. The increases in the Km for acetate for the variants indicate that these residues contribute to acetate affinity. The kcat values decreased 8-fold or less relative to the wild type for all variants, except for a 480-fold decrease for Phe179Ala. The kinetic parameters for acetyl phosphate (Table 3) were considerably less affected than the kinetic parameters for acetate by replacements in the hydrophobic pocket (Table 2). The values for the Km for acetyl phosphate ranged from 1.5-fold (Leu122Ala) to 3.5-fold (Pro232Ala) greater than the wild-type value, with the exception of Phe179Ala, which showed a 2.3-fold decrease in the Km. The moderate changes in the Km suggest that factors other than the putative hydrophobic pocket residues contribute to acetyl phosphate affinity. While kcat changed less than 2- to 3-fold for the majority of variants, the Phe179Ala variant showed 233- and 479-fold decreases in kcat for the acetate- and acetyl phosphate-forming directions, respectively, suggesting that this residue may have a catalytic role. While substantial increases in the Km for acetate were observed for all the variants relative to the wild type (Table 2), there was relatively little change in the Km values for propionate and butyrate (Table 4). In contrast, the kcat values for the Val93Ala and Val93Gly variants with propionate and butyrate increased relative to the wild type, while the kcat values for the same variants decreased with acetate (Table 2), indicating that the relationship between pocket size and substrate size correlates with the rate of catalysis. While a Km value is only an approximation of an enzyme's affinity for a substrate, in the absence of true dissociation constants, it can provide useful information.
|
View this table: [in a new window] |
TABLE 2. Kinetic parameters of wild-type and variant acetate kinases assayed in the direction of acetyl phosphate synthesis
|
|
View this table: [in a new window] |
TABLE 3. Kinetic parameters for wild-type and variant acetate kinases assayed in the direction of ATP synthesis
|
|
View this table: [in a new window] |
TABLE 4. Kinetic parameters for wild-type and variant acetate kinases assayed in the direction of propionyl phosphate and butyryl phosphate synthesis
|
![]() View larger version (35K): [in a new window] |
FIG. 3. Solvent-accessible cavities in wild-type acetate kinase and the modeled Val93Ala and Val93Gly variants. (A, C, and E) Wild-type acetate kinase and the Val93Ala and Val93Gly variants, respectively, with the active site cavity volume masks. (B, D, and F) Acetate, propionate, and butyrate modeled into wild-type acetate kinase and the Val93Ala and Val93Gly variants, respectively. All diagrams were generated by using MOLSCRIPT (20) and Raster3D (25).
|
|
|
|---|
The results obtained from modeling acetate in the previously published wild-type acetate kinase structure were consistent with the kinetic results, suggesting that the hydrophobic pocket has a role in binding the methyl group of acetate. The model also showed the carboxyl group of acetate pointing toward the sulfate ion which is the proposed site for the
phosphate of ATP (8). Furthermore, the carboxyl group of acetate in the model is poised to interact with the guanidinium cation of Arg91 previously proposed to facilitate the binding of acetate (14, 33). The modeling results are consistent with the recently determined crystal structure of M. thermophila acetate kinase containing ADP, acetate, and the transition state analog AlF3 (14). In this structure, acetate is bound with the methyl group in the pocket and the carboxyl group pointing toward AlF3, a proposed transition state analog of the
phosphate of ATP, supporting the hypothesis that there is a direct in-line mechanism. The kinetic and modeling results presented here are consistent with this mechanism.
In addition to an essential role for the hydrophobic pocket in binding acetate, the results suggest that the size of the pocket is a determinant of substrate specificity and also is important for catalysis. As the side chain volume decreased in the position 93 variants, the kcat increased for propionate and butyrate relative to the wild type, in accord with the modeling results. Notably, the kcat of the Val93Gly variant with butyrate was 1,600-fold greater than that for the wild type. Furthermore, the kcat for the Val93Ala variant with propionate was equivalent to that for the wild-type enzyme with acetate, although the catalytic efficiency of the former was lower due to the slightly larger Km for propionate. As kinetic studies have not been performed with any propionate kinase, the Km values and catalytic efficiencies cannot be compared. These results suggest that while hydrophobic interactions appear to be the main determinant of affinity for acetate, specificity for acetate is through limiting the size of the hydrophobic pocket, determined in large part by Val93 located in the floor of the pocket. These results also suggest that the side chain size at position 93 is important for catalysis in the direction of acetyl phosphate synthesis, most likely due to positioning of the carboxyl group of acetate in proximity to the
-phosphate of ATP for nucleophilic attack. However, the kcat values for the Leu122Ala, Phe179Ala, and Pro232Ala variants with acetate and propionate were substantially lower than the values for the position 93 variants, and no activity was detected with butyrate. These results suggest that the Leu122, Phe179, and Pro232 residues have greater impact than Val93 on substrate positioning and catalysis. The wild-type M. thermophila acetate kinase had substantial propionate kinase activity, which may have physiological relevance. It was recently reported that pduW encodes propionate kinase (29) in Salmonella enterica and that acetate kinase is able to function in propionate metabolism in the absence of pduW.
If the affinity of acetate kinase for acetyl phosphate occurs through the same hydrophobic interactions as the affinity for acetate, then similar changes would be expected in the Km values for acetate and propionate for the variants. However, the increases in the Km for acetyl phosphate were severalfold less than the increases in the Km for acetate for all of the variants, suggesting that the affinity of acetyl phosphate is mediated to a greater extent by interaction with the phosphoryl group than by hydrophobic interactions with the methyl group. This interpretation is further supported by the fivefold-lower Km for acetyl phosphate than for acetate for the wild-type enzyme. Except for the Phe179Ala variant, the kcat for the variants changed less in the direction of ATP synthesis than in the direction of acetyl phosphate synthesis, which is consistent with less influence for the hydrophobic pocket on positioning acetyl phosphate for optimum catalysis. Although Arg91 was hypothesized to interact with the phosphoryl group of acetyl phosphate based upon the crystal structure complexed with ADP (8), the kinetic analysis results for variants are inconsistent with this role for Arg91 (26, 33; unpublished data); thus, the residue interacting with acetyl phosphate remains unknown. Neither propionyl phosphate nor butyryl phosphate is commercially available; thus, the effect of the pocket size on the specificity for the phosphorylated substrates has not been tested.
Members of the ASKHA superfamily are known to undergo domain closure which is required for catalysis (15, 28, 30). The catalytically essential Phe179 residue, which is conserved in all three kinases (Fig. 2), is located at the N terminus of helix
3, which is at the bridge between domains I and II (8) (Fig. 4). The ASKHA family members glycerol kinase (10, 30), hexokinase (2), and phosphoglycerate kinase (5) have a conserved glycine residue in their equivalent helices (
3 in glycerol kinase and hexokinase,
7 in phosphoglycerate kinase) that is postulated to be part of the hinge at which the catalytically essential domain closure occurs (2, 5, 10, 30). In glycerol kinase and phosphoglycerate kinase a phenylalanine immediately follows this hinge glycine, and this motif is also present in the acetate kinase and propionate kinase sequences (Fig. 2). Therefore, this Gly-Phe pair may identify a previously unrecognized conserved motif required for domain motion in the ASKHA superfamily. Phe179 likely participates in the closure of domain II down onto domain I by creating a greasy hydrophobic patch that allows the sliding shear action required for movement (10, 15). The Leu122Ala variant also lies at the domain I-domain II interface (Fig. 4), implying that it contributes to the proposed hydrophobic patch. Indeed, the Leu122Ala variant showed a reduction in kcat relative to the wild type in the direction of ATP synthesis, which is consistent with a role in addition to the proposed positioning of the carboxyl group of acetate for optimal catalysis in the direction of acetyl phosphate synthesis.
![]() View larger version (34K): [in a new window] |
FIG. 4. Overall structure of acetate kinase, highlighting Phe179 and Leu122 in the hinge region. Monomer A is blue and green (domain I, green; domain II, blue), monomer B is orange and red (domain I, red; domain II, orange), and the 3 helices are cyan (monomer A) and yellow (monomer B). Phe179 and Leu122 are shown in ball-and-stick form. The figure was generated by using MOLSCRIPT (20) and Raster3D (25).
|
Conclusions. Kinetic analyses of acetate kinase variants have resulted in identification of residues that form a hydrophobic pocket and are important for binding acetate. The results further suggest that the size of the pocket is an important determinant of substrate specificity and support the hypothesis that the acetate identified in the recently reported crystal structure (14) is located in the catalytically competent position. Furthermore, the results reported here and the recently reported crystal structure both support a direct in-line mechanism. Finally, we identified a previously unrecognized catalytically essential residue, Phe179, and hypothesize that this residue has a role in domain closure during catalysis.
C.I-S. and A.G. contributed equally to this work. ![]()
Present address: Department of Genetics and Biochemistry, Clemson University, Clemson, SC 29634-0324. ![]()
Present address: Department of Chemistry, University of Northern British Columbia, Prince George, BC V2N 4Z9, Canada. ![]()
¶ Present address: Institute for Biological Energy Alternatives, Rockville, MD 20850-3343. ![]()
|
|
|---|
W in complex with nonhydrolyzable ATP analogues reveal a putative active conformation of the enzyme as a result of domain motion. Biochemistry 38:3508-3518.[CrossRef][Medline]
| |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
Copyright © 2009 by the American Society for Microbiology. For an alternate route to Journals.ASM.org, visit: http://intl-journals.asm.org | More Info»