Structures and mechanisms of Nudix hydrolases.

Nudix hydrolases catalyze the hydrolysis of nucleoside diphosphates linked to other moieties, X, and contain the sequence motif or Nudix box, GX(5)EX(7)REUXEEXGU. The mechanisms of Nudix hydrolases are highly diverse in the position on the substrate at which nucleophilic substitution occurs, and in the number of required divalent cations. While most proceed by associative nucleophilic substitutions by water at specific internal phosphorus atoms of a diphosphate or polyphosphate chain, members of the GDP-mannose hydrolase sub-family catalyze dissociative nucleophilic substitutions, by water, at carbon. The site of substitution is likely determined by the positions of the general base and the entering water. The rate accelerations or catalytic powers of Nudix hydrolases range from 10(9)- to 10(12)-fold. The reactions are accelerated 10(3)-10(5)-fold by general base catalysis by a glutamate residue within, or beyond the Nudix box, or by a histidine beyond the Nudix box. Lewis acid catalysis, which contributes 10(3)-10(5)-fold to the rate acceleration, is provided by one, two, or three divalent cations. One divalent cation is coordinated by two or three conserved residues of the Nudix box, the initial glycine and one or two glutamate residues, together with a remote glutamate or glutamine ligand from beyond the Nudix box. Some Nudix enzymes require one (MutT) or two additional divalent cations (Ap(4)AP), to neutralize the charge of the polyphosphate chain, to help orient the attacking hydroxide or oxide nucleophile, and/or to facilitate the departure of the anionic leaving group. Additional catalysis (10-10(3)-fold) is provided by the cationic side chains of lysine and arginine residues and by H-bond donation by tyrosine residues, to orient the general base, or to promote the departure of the leaving group. The overall rate accelerations can be explained by both independent and cooperative effects of these catalytic components.

Mechanistic implications of methylglyoxal synthase complexed with phosphoglycolohydroxamic acid as observed by X-ray crystallography and NMR spectroscopy.

Methylglyoxal synthase (MGS) and triosephosphate isomerase (TIM) share neither sequence nor structural similarities, yet the reactions catalyzed by both enzymes are similar, in that both initially convert dihydroxyacetone phosphate to a cis-enediolic intermediate. This enediolic intermediate is formed from the abstraction of the pro-S C3 proton of DHAP by Asp-71 of MGS or the pro-R C3 proton of DHAP by Glu-165 of TIM. MGS then catalyzes the elimination of phosphate from this enediolic intermediate to form the enol of methylglyoxal, while TIM catalyzes proton donation to C2 to form D-glyceraldehyde phosphate. A competitive inhibitor of TIM, phosphoglycolohydroxamic acid (PGH) is found to be a tight binding competitive inhibitor of MGS with a K(i) of 39 nM. PGH's high affinity for MGS may be due in part to a short, strong hydrogen bond (SSHB) from the NOH of PGH to the carboxylate of Asp-71. Evidence for this SSHB is found in X-ray, 1H NMR, and fractionation factor data. The X-ray structure of the MGS homohexamer complexed with PGH at 2.0 A resolution shows this distance to be 2.30-2.37 +/- 0.24 A. 1H NMR shows a PGH-dependent 18.1 ppm signal that is consistent with a hydrogen bond length of 2.49 +/- 0.02 A. The D/H fractionation factor (phi = 0.43 +/- 0.02) is consistent with a hydrogen bond length of 2.53 +/- 0.01 A. Further, 15N NMR suggests a significant partial positive charge on the nitrogen atom of bound PGH, which could strengthen hydrogen bond donation to Asp-71. Both His-98 and His-19 are uncharged in the MGS-PGH complex on the basis of the chemical shifts of their Cdelta and C(epsilon) protons. The crystal structure reveals that Asp-71, on the re face of PGH, and His-19, on the si face of PGH, both approach the NO group of the analogue, while His-98, in the plane of PGH, approaches the carbonyl oxygen of the analogue. The phosphate group of PGH accepts nine hydrogen bonds from seven residues and is tilted out of the imidate plane of PGH toward the re face. Asp-71 and phosphate are thus positioned to function as the base and leaving group, respectively, in a concerted suprafacial 1,4-elimination of phosphate from the enediolic intermediate in the second step of the MGS reaction. Combined, these data suggest that Asp-71 is the one base that initially abstracts the C3 pro-S proton from DHAP and subsequently the 3-OH proton from the enediolic intermediate. This mechanism is compared to an alternative TIM-like mechanism for MGS, and the relative merits of both mechanisms are discussed.

Short, strong hydrogen bonds at the active site of human acetylcholinesterase: proton NMR studies.

Cholinesterases use a Glu-His-Ser catalytic triad to enhance the nucleophilicity of the catalytic serine. We have previously shown by proton NMR that horse serum butyryl cholinesterase, like serine proteases, forms a short, strong hydrogen bond (SSHB) between the Glu-His pair upon binding mechanism-based inhibitors, which form tetrahedral adducts, analogous to the tetrahedral intermediates in catalysis [Viragh, C., et al. (2000) Biochemistry 39, 16200-16205]. We now extend these studies to human acetylcholinesterase, a 136 kDa homodimer. The free enzyme at pH 7.5 shows a proton resonance at 14.4 ppm assigned to an imidazole NH of the active-site histidine, but no deshielded proton resonances between 15 and 21 ppm. Addition of a 3-fold excess of the mechanism-based inhibitor m-(N,N,N-trimethylammonio)trifluoroacetophenone (TMTFA) induced the complete loss of the 14.4 ppm signal and the appearance of a broad, deshielded resonance of equal intensity with a chemical shift delta of 17.8 ppm and a D/H fractionation factor phi of 0.76 +/- 0.10, consistent with a SSHB between Glu and His of the catalytic triad. From an empirical correlation of delta with hydrogen bond lengths in small crystalline compounds, the length of this SSHB is 2.62 +/- 0.02 A, in agreement with the length of 2.63 +/- 0.03 A, independently obtained from phi. Upon addition of a 3-fold excess of the mechanism-based inhibitor 4-nitrophenyl diethyl phosphate (paraoxon) to the free enzyme at pH 7.5, and subsequent deethylation, two deshielded resonances of unequal intensity appeared at 16.6 and 15.5 ppm, consistent with SSHBs with lengths of 2.63 +/- 0.02 and 2.65 +/- 0.02 A, respectively, suggesting conformational heterogeneity of the active-site histidine as a hydrogen bond donor to either Glu-327 of the catalytic triad or to Glu-199, also in the active site. Conformational heterogeneity was confirmed with the methylphosphonate ester anion adduct of the active-site serine, which showed two deshielded resonances of equal intensity at 16.5 and 15.8 ppm with phi values of 0.47 +/- 0.10 and 0.49 +/- 0.10 corresponding to average hydrogen bond lengths of 2.59 +/- 0.04 and 2.61 +/- 0.04 A, respectively. Similarly, lowering the pH of the free enzyme to 5.1 to protonate the active-site histidine (pK(a) = 6.0 +/- 0.4) resulted in the appearance of two deshielded resonances, at 17.7 and 16.4 ppm, consistent with SSHBs with lengths of 2.62 +/- 0.02 and 2.63 +/- 0.02 A, respectively. The NMR-derived distances agree with those found in the X-ray structures of the homologous acetylcholinesterase from Torpedo californica complexed with TMTFA (2.66 +/- 0.28 A) and sarin (2.53 +/- 0.26 A) and at low pH (2.52 +/- 0.25 A). However, the order of magnitude greater precision of the NMR-derived distances establishes the presence of SSHBs at the active site of acetylcholinesterase, and detect conformational heterogeneity of the active-site histidine. We suggest that the high catalytic power of cholinesterases results in part from the formation of a SSHB between Glu and His of the catalytic triad.

The structural basis for the perturbed pKa of the catalytic base in 4-oxalocrotonate tautomerase: kinetic and structural effects of mutations of Phe-50.

The amino-terminal proline of 4-oxalocrotonate tautomerase (4-OT) functions as the general base catalyst in the enzyme-catalyzed isomerization of beta,gamma-unsaturated enones to their alpha,beta-isomers because of its unusually low pK(a) of 6.4 +/- 0.2, which is 3 units lower than that of the model compound, proline amide. Recent studies show that this abnormally low pK(a) is not due to the electrostatic effects of nearby cationic residues (Arg-11, Arg-39, and Arg-61) [Czerwinski, R. M., Harris, T. K., Johnson, Jr., W. H., Legler, P. M., Stivers, J. T., Mildvan, A. S., and Whitman, C. P. (1999) Biochemistry 38, 12358-12366]. Hence, it may result solely from a low local dielectric constant of 14.7 +/- 0.8 at the otherwise hydrophobic active site. Support for this mechanism comes from the study of mutants of the active site Phe-50, which is 5.8 A from Pro-1 and is one of 12 apolar residues within 9 A of Pro-1. Replacing Phe-50 with Tyr does not significantly alter k(cat) or K(m) and results in a pK(a) of 6.0 +/- 0.1 for Pro-1 as determined by (15)N NMR spectroscopy, comparable to that observed for wild type. (1)H-(15)N HSQC and 3D (1)H-(15)N NOESY HSQC spectra of the F50Y mutant demonstrate its conformation to be very similar to that of the wild-type enzyme. In the F50Y mutant, the pK(a) of Tyr-50 is increased by two units from that of a model compound N-acetyl-tyrosine amide to 12.2 +/- 0.3, as determined by UV and (1)H NMR titrations, yielding a local dielectric constant of 13.4 +/- 1.7, in agreement with the value of 13.7 +/- 0.3 determined from the decreased pK(a) of Pro-1 in this mutant. In the F50A mutant, the pK(a) of Pro-1 is 7.3 +/- 0.1 by (15)N NMR titration, comparable to the pK(a) of 7.6 +/- 0.2 found in the pH vs k(cat)/K(m) rate profile, and is one unit greater than that of the wild-type enzyme, indicating an increase in the local dielectric constant to a value of 21.2 +/- 2.6. A loss of structure of the beta-hairpin from residues 50 to 57, which covers the active site, and is the site of the mutation, is indicated by the disappearance in the F50A mutant of four interstrand NOEs and one turn NOE found in wild-type 4-OT. (1)H-(15)N HSQC spectra of the F50A mutant reveal widespread and large changes in the backbone (15)N and NH chemical shifts including those of Gly residues 48, 51, 53, and 54 causing their loss of dispersion at 23 degrees C and their disappearance at 43 degrees C due to rapid exchange with solvent. These observations confirm that the active site of the F50A mutant is more accessible to the external aqueous environment, causing an increase in the local dielectric constant and in the pK(a) of Pro-1. In addition, the F50A mutation decreased k(cat) 167-fold and increased K(m) 11-fold from those of the wild-type enzyme, suggesting an important role for the hydrophobic environment in catalysis, beyond that of decreasing the pK(a) of Pro-1. The F50I and F50V mutations destabilize the protein and decrease k(cat) by factors of 58 and 1.6, and increase K(m) by 3.3- and 3.8-fold, respectively.

NMR evidence for a short, strong hydrogen bond at the active site of a cholinesterase.

Cholinesterases (ChE), use a Glu-His-Ser catalytic triad to enhance the nucleophilicity of the catalytic serine. It has been shown that serine proteases, which employ an Asp-His-Ser catalytic triad for optimal catalytic efficiency, decrease the hydrogen bonding distance between the Asp-His pair to form a short, strong hydrogen bond (SSHB) upon binding mechanism-based inhibitors, which form tetrahedral Ser-adducts, analogous to the tetrahedral intermediates in catalysis, or at low pH when the histidine is protonated [Cassidy, C. S., Lin, J., Frey, P. A. (1997) Biochemistry 36, 4576-4584]. Two types of mechanism-based inhibitors were bound to pure equine butyrylcholinesterase (BChE), a 364 kDa homotetramer, and the complexes were studied by (1)H NMR at 600 MHz and 25-37 degrees C. The downfield region of the (1)H NMR spectrum of free BChE at pH 7.5 showed a broad, weak, deshielded resonance with a chemical shift, delta = 16.1 ppm, ascribed to a small amount of the histidine-protonated form. Upon addition of a 3-fold excess of diethyl 4-nitrophenyl phosphate (paraoxon) and subsequent dealkylation, the broad 16.1 ppm resonance increased in intensity 4.7-fold, and yielded a D/H fractionation factor phi = 0.72+/-0.10 consistent with a SSHB between Glu and His of the catalytic triad. From an empirical correlation of delta with hydrogen-bond length in small crystalline compounds, the length of this SSBH is 2.64+/-0.04 A, in agreement with the length of 2.62+/-0.02 A independently obtained from phi. The addition of a 3-fold excess of m-(N,N, N-trimethylammonio)trifluoroacetophenone to BChE yielded no signal at 16.1 ppm, and a 640 Hz broad, highly deshielded proton resonance with a chemical shift delta = 18.1 ppm and a D/H fractionation factor phi = 0.63+/-0.10, also consistent with a SSHB. The length of this SSHB is calculated to be 2.62+/-0.04 A from delta and 2.59+/-0.03 A from phi. These NMR-derived distances agree with those found in the X-ray structures of the homologous acetylcholinesterase complexed with the same mechanism-based inhibitors, 2.60+/-0.22 and 2.66+/-0.28 A. However, the order of magnitude greater precision of the NMR-derived distances establish the presence of SSHBs. We suggest that ChEs achieve their remarkable catalytic power in ester hydrolysis, in part, due to the formation of a SSHB between Glu and His of the catalytic triad.

GDP-mannose mannosyl hydrolase catalyzes nucleophilic substitution at carbon, unlike all other Nudix hydrolases.

GDP-mannose mannosyl hydrolase (GDPMH) from Escherichia coli is a 36. 8 kDa homodimer which, in the presence of Mg(2+), catalyzes the hydrolysis of GDP-alpha-D-mannose or GDP-alpha-D-glucose to yield sugar and GDP. On the basis of its amino acid sequence, GDPMH is a member of the Nudix family of enzymes which catalyze the hydrolysis of nucleoside diphosphate derivatives by nucleophilic substitution at phosphorus. However, GDPMH has a sequence rearrangement (RE to ER) in the conserved Nudix motif and is missing a Glu residue characteristic of the Nudix signature sequence. By (1)H NMR, the initial hydrolysis product of GDP-alpha-D-glucose is beta-D-glucose, indicating nucleophilic substitution with inversion at C1' of glucose. Substitution at carbon was confirmed by two-dimensional (1)H-(13)C HSQC spectra of the products of hydrolysis in 48.4% (18)O-labeled water which showed an additional C1' resonance of beta-D-glucose with a typical upfield (18)O isotope shift of 18 ppb and an intensity of 47.6% of the total signal. No (18)O isotope-shifted resonances (<4%) were found in the (31)P NMR spectrum of the GDP product. Thus, unlike all other Nudix enzymes studied so far, GDPMH catalyzes nucleophilic substitution at carbon rather than at phosphorus. A small solvent kinetic deuterium isotope effect on k(cat) of 1.76 +/- 0.25, independent of pH over the range of 6.0-9.3, suggests that the deprotonation of water may be part of the rate-limiting step.

Metal requirements of a diadenosine pyrophosphatase from Bartonella bacilliformis: magnetic resonance and kinetic studies of the role of Mn2+.

Recombinant IalA protein from Bartonella bacilliformis is a monomeric adenosine 5'-tetraphospho-5'-adenosine (Ap4A) pyrophosphatase of 170 amino acids that catalyzes the hydrolysis of Ap4A, Ap5A, and Ap6A by attack at the delta-phosphorus, with the departure of ATP as the leaving group [Cartwright et al. (1999) Biochem. Biophys. Res. Commun. 256, 474-479]. When various divalent cations were tested over a 300-fold concentration range, Mg2+, Mn2+, and Zn2+ ions were found to activate the enzyme, while Ca2+ did not. Sigmoidal activation curves were observed with Mn2+ and Mg2+ with Hill coefficients of 3.0 and 1.6 and K0.5 values of 0.9 and 5.3 mM, respectively. The substrate M2+ x Ap4A showed hyperbolic kinetics with Km values of 0.34 mM for both Mn2+ x Ap4A and Mg2+ x Ap4A. Direct Mn2+ binding studies by electron paramagnetic resonance (EPR) and by the enhancement of the longitudinal relaxation rate of water protons revealed two Mn2+ binding sites per molecule of Ap4A pyrophosphatase with dissociation constants of 1.1 mM, comparable to the kinetically determined K0.5 value of Mn2+. The enhancement factor of the longitudinal relaxation rate of water protons due to bound Mn2+ (epsilon b) decreased with increasing site occupancy from a value of 12.9 with one site occupied to 3.3 when both are occupied, indicating site-site interaction between the two enzyme-bound Mn2+ ions. Assuming the decrease in epsilon(b) to result from cross-relaxation between the two bound Mn2+ ions yields an estimated distance of 5.9 +/- 0.4 A between them. The substrate Ap4A binds one Mn2+ (Kd = 0.43 mM) with an epsilon b value of 2.6, consistent with the molecular weight of the Mn2+ x Ap4A complex. Mg2+ binding studies, in competition with Mn2+, reveal two Mg2+ binding sites on the enzyme with Kd values of 8.6 mM and one Mg2+ binding site on Ap4A with a Kd of 3.9 mM, values that are comparable to the K0.5 for Mg2+. Hence, with both Mn2+ and Mg2+, a total of three metal binding sites were found-two on the enzyme and one on the substrate-with dissociation constants comparable to the kinetically determined K0.5 values, suggesting a role in catalysis for three bound divalent cations. Ca2+ does not activate Ap4A pyrophosphatase but inhibits the Mn2+-activated enzyme competitively with a Ki = 1.9 +/- 1.3 mM. Ca2+ binding studies, in competition with Mn2+, revealed two sites on the enzyme with dissociation constants (4.3 +/- 1.3 mM) and one on Ap4A with a dissociation constant of 2.1 mM. These values are similar to its Ki suggesting that inhibition by Ca2+ results from the complete displacement of Mn2+ from the active site. Unlike the homologous MutT pyrophosphohydrolase, which requires only one enzyme-bound divalent cation in an E x M2+ x NTP x M2+ complex for catalytic activity, Ap4A pyrophosphatase requires two enzyme-bound divalent cations that function in an active E x (M2+)2 x Ap4A x M2+ complex.

Mutational, kinetic, and NMR studies of the roles of conserved glutamate residues and of lysine-39 in the mechanism of the MutT pyrophosphohydrolase.

The MutT enzyme catalyzes the hydrolysis of nucleoside triphosphates (NTP) to NMP and PP(i) by nucleophilic substitution at the rarely attacked beta-phosphorus. The solution structure of the quaternary E-M(2+)-AMPCPP-M(2+) complex indicated that conserved residues Glu-53, -56, -57, and -98 are at the active site near the bound divalent cation possibly serving as metal ligands, Lys-39 is positioned to promote departure of the NMP leaving group, and Glu-44 precedes helix I (residues 47-59) possibly stabilizing this helix which contributes four catalytic residues to the active site [Lin, J. , Abeygunawardana, C., Frick, D. N., Bessman, M. J., and Mildvan, A. S. (1997) Biochemistry 36, 1199-1211]. To test these proposed roles, the effects of mutations of each of these residues on the kinetic parameters and on the Mn(2+), Mg(2+), and substrate binding properties were examined. The largest decreases in k(cat) for the Mg(2+)-activated enzyme of 10(4.7)- and 10(2.6)-fold were observed for the E53Q and E53D mutants, respectively, while 97-, 48-, 25-, and 14-fold decreases were observed for the E44D, E56D, E56Q, and E44Q mutations, respectively. Smaller effects on k(cat) were observed for mutations of Glu-98 and Lys-39. For wild type MutT and its E53D and E44D mutants, plots of log(k(cat)) versus pH exhibited a limiting slope of 1 on the ascending limb and then a hump, i.e., a sharply defined maximum near pH 8 followed by a plateau, yielding apparent pK(a) values of 7.6 +/- 0.3 and 8.4 +/- 0.4 for an essential base and a nonessential acid catalyst, respectively, in the active quaternary MutT-Mg(2+)-dGTP-Mg(2+) complex. The pK(a) of 7.6 is assigned to Glu-53, functioning as a base catalyst in the active quaternary complex, on the basis of the disappearance of the ascending limb of the pH-rate profile of the E53Q mutant, and its restoration in the E53D mutant with a 10(1.9)-fold increase in (k(cat))(max). The pK(a) of 8.4 is assigned to Lys-39 on the basis of the disappearance of the descending limb of the pH-rate profile of the K39Q mutant, and the observation that removal of the positive charge of Lys-39, by either deprotonation or mutation, results in the same 8.7-fold decrease in k(cat). Values of k(cat) of both wild type MutT and the E53Q mutant were independent of solvent viscosity, indicating that a chemical step is likely to be rate-limiting with both. A liganding role for Glu-53 and Glu-56, but not Glu-98, in the binary E-M(2+) complex is indicated by the observation that the E53Q, E53D, E56Q, and E56D mutants bound Mn(2+) at the active site 36-, 27-, 4.7-, and 1.9-fold weaker, and exhibited 2.10-, 1.50-, 1.12-, and 1.24-fold lower enhanced paramagnetic effects of Mn(2+), respectively, than the wild type enzyme as detected by 1/T(1) values of water protons, consistent with the loss of a metal ligand. However, the K(m) values of Mg(2+) and Mn(2+) indicate that Glu-56, and to a lesser degree Glu-98, contribute to metal binding in the active quaternary complex. Mutations of the more distant but conserved residue Glu-44 had little effect on metal binding or enhancement factors in the binary E-M(2+) complexes. Two-dimensional (1)H-(15)N HSQC and three-dimensional (1)H-(15)N NOESY-HSQC spectra of the kinetically damaged E53Q and E56Q mutants showed largely intact proteins with structural changes near the mutated residues. Structural changes in the kinetically more damaged E44D mutant detected in (1)H-(15)N HSQC spectra were largely limited to the loop I-helix I motif, suggesting that Glu-44 stabilizes the active site region. (1)H-(15)N HSQC titrations of the E53Q, E56Q, and E44D mutants with dGTP showed changes in chemical shifts of residues lining the active site cleft, and revealed tighter nucleotide binding by these mutants, indicating an intact substrate binding site. (ABSTRACT TRUNCATED)

Kinetic, stereochemical, and structural effects of mutations of the active site arginine residues in 4-oxalocrotonate tautomerase.

Three arginine residues (Arg-11, Arg-39, Arg-61) are found at the active site of 4-oxalocrotonate tautomerase in the X-ray structure of the affinity-labeled enzyme [Taylor, A. B., Czerwinski, R. M., Johnson, R. M., Jr., Whitman, C. P., and Hackert, M. L. (1998) Biochemistry 37, 14692-14700]. The catalytic roles of these arginines were examined by mutagenesis, kinetic, and heteronuclear NMR studies. With a 1,6-dicarboxylate substrate (2-hydroxymuconate), the R61A mutation showed no kinetic effects, while the R11A mutation decreased k(cat) 88-fold and increased K(m) 8.6-fold, suggesting both binding and catalytic roles for Arg-11. With a 1-monocarboxylate substrate (2-hydroxy-2,4-pentadienoate), no kinetic effects of the R11A mutation were found, indicating that Arg-11 interacts with the 6-carboxylate of the substrate. The stereoselectivity of the R11A-catalyzed protonation at C-5 of the dicarboxylate substrate decreased, while the stereoselectivity of protonation at C-3 of the monocarboxylate substrate increased in comparison with wild-type 4-OT, indicating the importance of Arg-11 in properly orienting the dicarboxylate substrate by interacting with the charged 6-carboxylate group. With 2-hydroxymuconate, the R39A and R39Q mutations decreased k(cat) by 125- and 389-fold and increased K(m) by 1.5- and 2.6-fold, respectively, suggesting a largely catalytic role for Arg-39. The activity of the R11A/R39A double mutant was at least 10(4)-fold lower than that of the wild-type enzyme, indicating approximate additivity of the effects of the two arginine mutants on k(cat). For both R11A and R39Q, 2D (1)H-(15)N HSQC and 3D (1)H-(15)N NOESY-HSQC spectra showed chemical shift changes mainly near the mutated residues, indicating otherwise intact protein structures. The changes in the R39Q mutant were mainly in the beta-hairpin from residues 50 to 57 which covers the active site. HSQC titration of R11A with the substrate analogue cis, cis-muconate yielded a K(d) of 22 mM, 37-fold greater than the K(d) found with wild-type 4-OT (0.6 mM). With the R39Q mutant, cis, cis-muconate showed negative cooperativity in active site binding with two K(d) values, 3.5 and 29 mM. This observation together with the low K(m) of 2-hydroxymuconate (0.47 mM) suggests that only the tight binding sites function catalytically in the R39Q mutant. The (15)Nepsilon resonances of all six Arg residues of 4-OT were assigned, and the assignments of Arg-11, -39, and -61 were confirmed by mutagenesis. The binding of cis,cis-muconate to wild-type 4-OT upshifts Arg-11 Nepsilon (by 0.05 ppm) and downshifts Arg-39 Nepsilon (by 1.19 ppm), indicating differing electronic delocalizations in the guanidinium groups. A mechanism is proposed in which Arg-11 interacts with the 6-carboxylate of the substrate to facilitate both substrate binding and catalysis and Arg-39 interacts with the 1-carboxylate and the 2-keto group of the substrate to promote carbonyl polarization and catalysis, while Pro-1 transfers protons from C-3 to C-5. This mechanism, together with the effects of mutations of catalytic residues on k(cat), provides a quantitative explanation of the 10(7)-fold catalytic power of 4-OT. Despite its presence in the active site in the crystal structure of the affinity-labeled enzyme, Arg-61 does not play a significant role in either substrate binding or catalysis.

Effects of mutations of the active site arginine residues in 4-oxalocrotonate tautomerase on the pKa values of active site residues and on the pH dependence of catalysis.

The unusually low pK(a) value of the general base catalyst Pro-1 (pK(a) = 6.4) in 4-oxalocrotonate tautomerase (4-OT) has been ascribed to both a low dielectric constant at the active site and the proximity of the cationic residues Arg-11 and Arg-39 [Stivers, J. T., Abeygunawardana, C., Mildvan, A. S., Hajipour, G., and Whitman, C. P. (1996) Biochemistry 35, 814-823]. In addition, the pH-rate profiles in that study showed an unidentified protonated group essential for catalysis with a pK(a) of 9.0. To address these issues, the pK(a) values of the active site Pro-1 and lower limit pK(a) values of arginine residues were determined by direct (15)N NMR pH titrations. The pK(a) values of Pro-1 and of the essential acid group were determined independently from pH-rate profiles of the kinetic parameters of 4-OT in arginine mutants of 4-OT and compared with those of wild type. The chemical shifts of all of the Arg Nepsilon resonances in wild-type 4-OT and in the R11A and R39Q mutants were found to be independent of pH over the range 4.9-9.7, indicating that no arginine is responsible for the kinetically determined pK(a) of 9.0 for an acidic group in free 4-OT. With the R11A mutant, where k(cat)/K(m) was reduced by a factor of 10(2.9), the pK(a) of Pro-1 was not significantly altered from that of the wild-type enzyme (pK(a) = 6.4 +/- 0.2) as revealed by both direct (15)N NMR titration (pK(a) = 6.3 +/- 0.1) and the pH dependence of k(cat)/K(m) (pK(a) = 6.4 +/- 0.2). The pH-rate profiles of both k(cat)/K(m) and k(cat) for the reaction of the R11A mutant with the dicarboxylate substrate, 2-hydroxymuconate, showed humps, i.e., sharply defined maxima followed by nonzero plateaus. The humps disappeared in the reaction with the monocarboxylate substrate, 2-hydroxy-2,4-pentadienoate, indicating that, unlike the wild-type enzyme which reacts only with the dianionic form of the dicarboxylic substrate, the R11A mutant reacts with both the 6-COOH and 6-COO(-) forms, with the 6-COOH form being 12-fold more active. This reversal in the preferred ionization state of the 6-carboxyl group of the substrate that occurs upon mutation of Arg-11 to Ala provides strong evidence that Arg-11 interacts with the 6-carboxylate of the substrate. In the R39Q mutant, where k(cat)/K(m) was reduced by a factor of 10(3), the kinetically determined pK(a) value for Pro-1 was 4.6 +/- 0.2, while the ionization of Pro-1 showed negative cooperativity with an apparent pK(a) of 7.1 +/- 0.1 determined by 1D (15)N NMR. From the Hill coefficient of 0.54, it can be shown that the apparent pK(a) value of 7.1 could result most simply from the averaging of two limiting pK(a) values of 4.6 and 8.2. Mutation of Arg-39, by altering the structure of the beta-hairpin which covers the active site, could result in an increase in the solvent exposure of Pro-1, raising its upper limit pK(a) value to 8.2. In the R39A mutant, the kinetically determined pK(a) of Pro-1 was also low, 5.0 +/- 0.2, indicating that in both the R39Q and R39A mutants, only the sites with low pK(a) values were kinetically operative. With the fully active R61A mutant, the kinetically determined pK(a) of Pro-1 (pK(a) = 6.5 +/- 0.2) agreed with that of wild-type 4-OT. It is concluded that the unusually low pK(a) of Pro-1 shows little contribution from electrostatic effects of the nearby cationic Arg-11, Arg-39, and Arg-61 residues but results primarily from a site of low local dielectric constant.

Cystic fibrosis transmembrane conductance regulator: solution structures of peptides based on the Phe508 region, the most common site of disease-causing DeltaF508 mutation.

Most cases of cystic fibrosis (CF), a common inherited disease of epithelial cell origin, are caused by the deletion of Phe508 located in the first nucleotide-binding domain (NBF1) of the protein called CFTR (cystic fibrosis transmembrane conductance regulator). To gain greater insight into the structure within the Phe508 region of the wild-type protein and the change in structure that occurs when this residue is deleted, we conducted nuclear magnetic resonance (NMR) studies on representative synthetic 26 and 25 amino acid peptide segments. 2D 1H NMR studies at 600 MHz of the 26-residue peptide consisting of Met498 to Ala523 in 10% DMSO, pH 4.0, at 25 degrees C show a continuous but labile helix from Gly500 to Lys522, based on both NH-NH(i,i+1) and alphaH-NH(i,i+1) NOEs. Phe508 within this helix shows only short-range (i,

High-precision measurement of hydrogen bond lengths in proteins by nuclear magnetic resonance methods.

We have compared hydrogen bond lengths on enzymes derived with high precision (< or = +/- 0.05 A) from both the proton chemical shifts (delta) and the fractionation factors (phi) of the proton involved with those obtained from protein X-ray crystallography. Hydrogen bond distances derived from proton chemical shifts were obtained from a correlation of 59 O--H....O hydrogen bond lengths, measured by small molecule high-resolution X-ray crystallography, with chemical shifts determined by solid-state nuclear magnetic resonance (NMR) in the same crystals (McDermott A, Ridenour CF, Encyclopedia of NMR, Sussex, U.K.: Wiley, 1996:3820-3825). Hydrogen bond distances were independently obtained from fractionation factors that yield distances between the two proton wells in quartic double minimum potential functions (Kreevoy MM, Liang TM, J Am Chem Soc, 1980;102:3315-3322). The high-precision hydrogen bond distances derived from their corresponding NMR-measured proton chemical shifts and fractionation factors agree well with each other and with those reported in protein X-ray structures within the larger errors (+/-0.2-0.8 A) in distances obtained by protein X-ray crystallography. The increased precision in measurements of hydrogen bond lengths by NMR has provided insight into the contributions of short, strong hydrogen bonds to catalysis for several enzymatic reactions.

Solution structure and mechanism of the MutT pyrophosphohydrolase.

The MutT enzyme prevents errors in DNA replication by hydrolyzing mutagenic nucleotide substrates such as 8-oxo-dGTP. It does so by catalyzing nucleophilic attack at the electron rich P beta of nucleoside triphosphates. Members of this small mechanistic class of enzymes require two divalent cations per active site for activity--one coordinated by the enzyme and the other by the enzyme-bound NTP--and show low catalytic powers of 10(7)- to 10(9)-fold. The first structure of an enzyme of this class, obtained by NMR methods in solution, shows MutT to be a compact globular protein with an alpha + beta-fold. The binding of the essential divalent cation activator Mg2+ and the substrate analog Mg(2+)-AMPCPP to the MutT enzyme to form the quaternary E-Mg(2+)-AMPCPP-Mg2+ complex does not alter the global fold of the enzyme but produces localized small conformational changes at or near the metal and substrate binding sites. The adenine-ribose moiety binds in a hydrophobic cleft near 3-strands of a mixed beta-sheet, with the 6-NH2 group of the purine ring approaching the -NH2 side chain of Asn-119. With a 6-keto group, GTP would interact more favorably with Asn-119, consistent with the substrate preference of MutT for guanine over adenine nucleotides. The enzyme-bound metal is coordinated by three conserved Glu residues (Glu-56, Glu-57, and Glu-98), the backbone carbonyl of a conserved Gly residue (Gly-38), and by two water ligands. The metal-triphosphate moiety of the metal-AMPCPP complex binds in the second coordination sphere of the enzyme-bound divalent cation. One of the water ligands of the enzyme-bound metal ion is well positioned to attack P beta with inversion and to be deprotonated or oriented by Glu-53, which in turn may be oriented by Arg-52. Lys-39 is positioned to interact electrostatically with the alpha-phosphoryl group and thereby to facilitate the departure of the NMP-leaving group. Quantitatively, the 10(9)-fold rate acceleration produced by the MutT enzyme may be ascribed to catalysis by approximation and polarization of the attacking water by the enzyme-bound metal ion (> or = 10(5)-fold), activation of the NMP leaving group by Lys-39 (10-fold), charge neutralization at P beta by the nucleotide-bound divalent cation (> or = 10-fold), and orientation and/or deprotonation of the attacking water by Glu-53 (> or = 10(2)-fold). This reaction mechanism, derived from the solution structure of the quaternary MutT complex, is both qualitatively and quantitatively consistent with the results of mutagenesis studies and may well be applicable to other enzymes that catalyze nucleophilic substitution at the electron-rich P beta of NTP substrates.

Proton transfer in the mechanism of triosephosphate isomerase.

Triosephosphate isomerase (TIM) catalyzes the reversible interconversion of dihydroxyacetone phosphate (DHAP) and glyceraldehyde 3-phosphate (GAP), with Glu-165 removing the pro-R proton from C1 of DHAP and neutral His-95 polarizing the carbonyl group of the substrate. During the TIM reaction, approximately 2% of the pro-R tritium from C1 of DHAP is conserved and appears at C2 of GAP [Nickbarg, E. B., and Knowles, J. R. (1988) Biochemistry 27, 5939]. In the "classical" mechanism, 98% of the pro-R tritium exchanges with solvent from Glu-165 at the intermediate state and the remaining 2% is transferred by Glu-165 to C2 of the same substrate molecule. This intramolecular transfer of tritium is therefore predicted to be independent of DHAP concentration. On the basis of NMR detection of a strong hydrogen bond between Glu-165 and the 1-OH of an analogue of the enediol intermediate [Harris, T. K., Abeygunawardana, C., and Mildvan, A. S. (1997) Biochemistry 36, 14661], we have suggested a "criss-cross" mechanism for TIM in which Glu-165 transfers a proton from C1 of DHAP to O2 of the enediol, and subsequently from O1 of the enediol to C2 of the product GAP. Since the pro-R proton is transferred to O2 instead of C2 in the criss-cross mechanism, no intramolecular transfer of label from substrate to product would be expected to occur. However, intermolecular transfer of label could occur if the label exchanges from O2 into a group on the protein and is transferred to GAP in subsequent turnovers. The extent of intermolecular tritium transfer in the criss-cross mechanism would be predicted to be dependent on DHAP concentration. The extent of tritium transfer was studied as a function of initial DHAP concentration using DHAP highly tritiated at the pro-R position. At 50% conversion to GAP, triphasic tritium transfer behavior was found. For phase 1, between 0.03 and 0.3 mM DHAP, a constant extent of tritium transfer of 1.19 +/- 0.03% occurred. For phase 2, between 0.3 and 1.0 mM DHAP, the extent of transfer progressively increased as a function of DHAP concentration to 2.17 +/- 0.15%. For phase 3, between 1.0 and 7.0 mM DHAP, the extent of transfer slightly decreased to 1.68 +/- 0.17%. In a direct test for intermolecular isotope transfer, doubly labeled [1(R)-D, 13C3]DHAP and 13C-depleted [1(R)-H,12C3]DHAP were synthesized, mixed in equal amounts, and incubated at 1 mM total DHAP with TIM, GAP dehydrogenase, NAD+, and arsenate until 50% conversion to 3-phosphoglycerate occurred. Electrospray ionization mass spectral analysis of the stable 3-phosphoglycerate product detected an extent of 1.4 +/- 0.4% of intramolecular D transfer from [13C3]DHAP to the 13C3 product, but no intermolecular transfer (

Solution structure of Delta 5-3-ketosteroid isomerase complexed with the steroid 19-nortestosterone hemisuccinate.

The solution structure of the ketosteroid isomerase homodimer complexed with the product analogue 19-nortestosterone hemisuccinate (19-NTHS) was solved by heteronuclear multidimensional NMR methods using 1647 distance restraints, 77 dihedral angle (phi) restraints, and 67 hydrogen bond restraints per monomer. The refined secondary structure of each subunit consists of three alpha-helices, eight beta-strands, four turns, and two beta-bulges. The beta-strands form a mixed beta-sheet. One of the five proline residues, Pro-39, is cis and begins a nonclassical turn. A self-consistent ensemble of 15 tertiary/quaternary structures of the enzyme dimer-steroid complex, with no distance violations greater than 0.35 A, was generated by simulated annealing and energy minimization with the program X-PLOR. The mean pairwise RMSD of the secondary structural elements was 0.63 A for the average subunit and 1.25 A for the dimer. Within each subunit, the three alpha-helices are packed onto the concave surface of the beta-sheet with a groove between them into which the steroid binds at a site defined by 14 intermolecular distances. In the productive complex, Tyr-14, from alpha-helix 1, approaches both Asp-99 and the 3-keto group of 19-NTHS while, from beta-strand 1, the carboxylate of Asp-38 approaches the beta-face of the steroid near C4 and C6, between which it transfers a proton during catalysis. Thus the solution structure of the isomerase-steroid complex can accommodate the catalytic diad mechanism in which Asp-99 donates a hydrogen bond to Tyr-14 which in turn is hydrogen bonded to the 3-oxygen of the steroid. While direct hydrogen bonding of Asp-99 to the steroid oxygen is less likely, it cannot be excluded. All other interactions of the steroid with the enzyme are hydrophobic. The dimer interface, which is between the convex surfaces of the beta-sheets, is defined by 28 intersubunit NOEs between hydrophobic residues in the 13C-filtered NOESY-HSQC spectrum of a 13C/12C-heterolabeled dimer. Both hydrophobic and polar interactions occur at the dimer interface which contains no space that would permit additional steroid binding. Comparison of the complexed enzyme with the solution structure of the free enzyme [Wu et al. (1997) Science 276, 415-418] reveals that the three helices change position in the steroid complex, becoming more closely packed onto the concave surface of the beta-sheet, thus bringing Tyr-14 closer to Asp-99 and the substrate. Comparison of the enzyme-steroid complex in solution with the free enzyme in the crystalline state reveals similar differences between the positions of the helices.

Hydrogen bonding at the active site of delta 5-3-ketosteroid isomerase.

The solution secondary structure of the highly active Y55F/Y88F "Tyr-14-only" mutant of delta 5-3-ketosteroid isomerase complexed with 19-nortestosterone hemisuccinate has been shown to consist of three helices, a six-stranded mixed beta-sheet, and five turns. The steroid binds near the general acid, Tyr-14, on helix 1, near the general base, Asp-38, on the first strand of the beta-sheet, and on the hydrophobic face of the beta-sheet [Zhao, Q., Abeygunawardana, C., & Mildvan, A. S. (1997) Biochemistry 36, 3458-3472]. On this hydrophobic face, Asp-99 is the only polar residue. Free isomerase shows a deshielded exchangeable proton resonance at 13.1 ppm assigned to the N epsilon H of neutral His-100. Its fractionation factor (phi = 0.79) and slow exchange with solvent suggest it to be buried or involved in an H-bond. The binding of dihydroequilenin or estradiol to isomerase induces the appearance of two additional deshielded proton resonances, one at 18.2 ppm assigned to the gamma-carboxyl proton of Asp-99, and the other, at 11.6 ppm, assigned to the zeta-OH proton of Tyr-14. While mutation of Asp-99 to Ala results in the disappearance of only the resonance near 18 ppm [Wu, R. W., Ebrahemian, S., Zwrotny, M. E., Thornberg, L. D., Perez-Alverado, G. C., Brothers, P., Pollack, R. M., & Summers, M. F. (1997) Science 276, 415-418], both of these resonances disappear in mutants lacking Tyr-14, suggesting an H-bonded catalytic diad, Asp-99-COOH--Tyr14-OH--O-steroid enolate. The catalytic diad is further supported by NOEs from the beta 1 and beta 2 protons of Asp-99 to the epsilon protons of Tyr-14, and from the zeta-OH proton of Tyr-14 to the gamma-carboxyl proton of Asp-99, indicating close proximity of these two residues, and by other data from the literature. A strong, low-barrier H-bond between Asp-99 and Tyr-14 is indicated by the 6.2 ppm deshielding, low fractionation factor (phi = 0.34) and slow exchange of the resonance at 18.2 ppm. A normal H-bond between Tyr-14 and the steroid is indicated by the 1.8 ppm deshielding, fractionation factor of 0.97 and the slow exchange of the resonance at 11.6 ppm. It is suggested that the 10(4.7)-fold contribution of Tyr-14 to catalysis is made possible by strong H-bonding from Asp-99 in the catalytic diad which strengthens general acid catalysis by Tyr-14. It is also noted that highly deshielded proton resonance on enzymes between 15 and 20 ppm, assigned to low-barrier H-bonds, generally involve carboxyl groups.

NMR studies of the role of hydrogen bonding in the mechanism of triosephosphate isomerase.

Triosephosphate isomerase (TIM) catalyzes the reversible interconversion of dihydroxyacetone phosphate (DHAP) and glyceraldehyde-3-phosphate (GAP), with Glu-165 removing the pro-R proton from C1 of DHAP and neutral His-95 polarizing the carbonyl group of the substrate. TIM and its complexes with the reactive intermediate analogs, phosphoglycolic acid (PGA) and phosphoglycolohydroxamic acid (PGH), were studied by 1H NMR at 600 MHz and at low temperature (-4.8 degrees C). His-95 shows an N epsilon H resonance at 13.1 ppm which shifts to 13.3 ppm in the TIM-PGA complex and to 13.5 ppm in the TIM-PGH complex. In the TIM-PGH complex, His-95 N epsilon H shows a slow, pH-independent exchange rate with water (kex = 80 s-1 at 30 degrees C, Eact = 19 kcal/mol), which is 44-fold slower than that of an exposed histidine suggesting partial shielding from bulk solvent, and a fractionation factor phi = 0.71 +/- 0.02 consistent with its donation of a normal hydrogen bond. The formation of the TIM-PGH complex results in the appearance of several deshielded proton resonances, including one at 14.9 ppm and one at 10.9 ppm which overlaps with another resonance. The resonance at 14.9 ppm is absent and the resonance at 10.9 ppm is much weaker in the TIM complex of PGA, which lacks the hydroxamic acid (-NHOH) moiety. 15N-labeled PGH was synthesized and the NH proton of free [15N]PGH shows a single 1H-15N HMQC cross peak with delta (1H) = 10.3 ppm and delta (15N) = 168 ppm which shifts to delta (1H) = 10.9 ppm and delta (15N) = 174 ppm in the TIM complex of [15N]PGH. The 15N-1H coupling in the complex indicates covalent N-H bonding, and the deshielded delta (15N) indicates a significant contribution of the imidate resonance form of PGH. The 14.9 ppm resonance is assigned to the NOH proton of bound PGH. This resonance shows a pH-independent exchange rate with water (kex = 3900 s-1 at 30 degrees C, Eact = 8.9 kcal/mol) which may reflect the dissociation of the TIM-PGH complex, and meets the criteria for a low-barrier hydrogen bond on the basis of the significant downfield shift of 6.2 ppm from the NOH proton of the model compound acetohydroxamic acid, and a very low fractionation factor phi = 0.38 +/- 0.06. In the X-ray structure of the TIM-PGH complex [Davenport, R. C., Bash, P. A., Seaton, B. A., Karplus, M., Petsko, G. A., and Ringe, D. (1991) Biochemistry 30, 5821], the NOH proton of bound PGH is hydrogen bonded to Glu-165. A low-barrier hydrogen bond from PGH NOH to Glu-165 suggests a dual role for Glu-165 in catalysis of proton transfer not only between the C1 and C2 carbons but also between the O1 and O2 oxygens in the interconversion of DHAP and GAP in wild type TIM. Such a mechanism, together with the measured exchange rate of the His-95 N epsilon H proton with solvent protons can accommodate the classical measurements of tritium incorporation from DHAP into GAP.

Mechanisms of signaling and related enzymes.

Most enzymes involved in cell signaling, such as protein kinases, protein phosphatases, GTPases, and nucleotide cyclases catalyze nucleophilic substitutions at phosphorus. When possible, the mechanisms of such enzymes are most clearly described quantitatively in terms of how associative or dissociative they are. The mechanisms of cell signaling enzymes range from < or = 8% associative (cAMP-dependent protein kinase) to approximately 50% associative (G protein Gi alpha 1). Their catalytic powers range from 10(5.7) (p21ras) to 10(11.7) (lambda Ser-Thr protein phosphatase), usually comparable in magnitude with those of nonsignaling enzymes of the same mechanistic class. Exceptions are G proteins, which are 10(3)- to 10(5)-fold poorer catalysts than F1 and myosin ATPases. The lower catalytic powers of G proteins may be ascribed to the absence of general base catalysis, and additionally in the case of p21ras, to the absence of a catalytic Arg residue, which interacts with the transition state. From kinetic studies of mutant and metal ion substituted enzymes, the catalytic powers of cell signaling and related enzymes can be rationalized quantitatively by factors contributed by metal ion catalysis (> or = 10(5), general acid catalysis (approximately 10(3 +/- 1)), general base catalysis (approximately 10(3 +/- 1)), and transition-state stabilization by cationic and hydrogen bond donating residues (approximately 10(3 +/- 1)).

Vaccinia DNA topoisomerase I: evidence supporting a free rotation mechanism for DNA supercoil relaxation.

The Vaccinia type I topoisomerase catalyzes site-specific DNA strand cleavage and religation by forming a transient phosphotyrosyl linkage between the DNA and Tyr-274, resulting in the release of DNA supercoils. For type I topoisomerases, two mechanisms have been proposed for supercoil release: (I) a coupled mechanism termed strand passage, in which a single supercoil is removed per cleavage/religation cycle, resulting in multiple topoisomer intermediates and late product formation, or (2) an uncoupled mechanism termed free rotation, where multiple supercoils are removed per cleavage/religation cycle, resulting in few intermediates and early product formation. To determine the mechanism, single-turnover experiments were done with supercoiled plasmid DNA under conditions in which the topoisomerase cleaves predominantly at a single site per DNA molecule. The concentrations of supercoiled substrate, intermediate topoisomers, and relaxed product vs time were measured by fluorescence imaging, and the rate constants for their interconversion were determined by kinetic simulation. Few intermediates and early product formation were observed. From these data, the rate constants for cleavage (0.3 s(-1)), religation (4 s(-1)), and the cleavage equilibrium constant on the enzyme (0.075) at 22 degrees C are in reasonable agreement with those obtained with small oligonucleotide substrates, while the rotation rate of the cleaved DNA strand is fast (approximately 20 rotations/s). Thus, the average number of supercoils removed for each cleavage event greatly exceeds unity (delta n = 5) and depends on kinetic competition between religation and supercoil release, establishing a free rotation mechanism. This free rotation mechanism for a type I topoisomerase differs from the strand passage mechanism proposed for the type II enzymes.