Active site hydrophobic residues impact hydrogen tunneling differently in a thermophilic alcohol dehydrogenase at optimal versus nonoptimal temperatures

A growing body of data suggests that protein motion plays an important role in enzyme catalysis. Two highly conserved hydrophobic active site residues in the cofactor-binding pocket of ht-ADH (Leu176 and V260) have been mutated to a series of hydrophobic side chains of smaller size, as well as one deletion mutant, L176Δ. Mutations decrease k(cat) and increase K(M)(NAD(+)). Most of the observed decreases in effects on k(cat) at pH 7.0 are due to an upward shift in the optimal pH for catalysis; a simple electrostatic model is invoked that relates the change in pK(a) to the distance between the positively charged nicotinamide ring and bound substrate. Structural modeling of the L176Δ and V260A variants indicates the development of a cavity behind the nicotinamide ring without any significant perturbation of the secondary structure of the enzyme relative to that of the wild type. Primary kinetic isotope effects (KIEs) are modestly increased for all mutants. Above the dynamical transition at 30 °C for ht-ADH [Kohen, A., et al. (1999) Nature 399, 496], the temperature dependence of the KIE is seen to increase with a decrease in side chain volume at positions 176 and 260. Additionally, the relative trends in the temperature dependence of the KIE above and below 30 °C appear to be reversed for the cofactor-binding pocket mutants in relation to wild-type protein. The aggregate results are interpreted in the context of a full tunneling model of enzymatic hydride transfer that incorporates both protein conformational sampling (preorganization) and active site optimization of tunneling (reorganization). The reduced temperature dependence of the KIE in the mutants below 30 °C indicates that at low temperatures, the enzyme adopts conformations refractory to donor-acceptor distance sampling.     

Kinetic investigations of the rate-limiting step in human 12- and 15-lipoxygenase

Mammalian lipoxygenases have been implicated in several inflammatory disorders; however, the details of the kinetic mechanism are still not well understood. In this paper, human platelet 12-lipoxygenase (12-hLO) and human reticulocyte 15-lipoxygenase-1 (15-hLO) were tested with arachidonic acid (AA) and linoleic acid (LA), respectively, under a variety of changing experimental conditions, such as temperature, dissolved oxygen concentration, and viscosity. The data that are presented show that 12-hLO and 15-hLO have slower rates of product release (k(cat)) than soybean lipoxygenase-1 (sLO-1), but similar or better rates of substrate capture for the fatty acid (k(cat)/K(M)) or molecular oxygen [k(cat)/K(M(O)2)]. The primary, kinetic isotope effect (KIE) for 15-hLO with LA was determined to be temperature-independent and large ((D)k(cat) = 40 +/- 8), over the range of 10-35 degrees C, indicating that C-H bond cleavage is the sole rate-limiting step and proceeds through a tunneling mechanism. The (D)k(cat)/K(M) for 15-hLO, however, was temperature-dependent, consistent with our previous results [Lewis, E. R., Johansen, E., and Holman, T. R. (1999) J. Am. Chem. Soc. 121, 1395-1396], indicating multiple rate-limiting steps. This was confirmed by a temperature-dependent, k(cat)/K(M) solvent isotope effect (SIE), which indicated a hydrogen bond rearrangement step at low temperatures, similar to that of sLO-1 [Glickman, M. H., and Klinman, J. P. (1995) Biochemistry 34, 14077-14092]. The KIE could not be determined for 12-hLO due to its inability to efficiently catalyze LA, but the k(cat)/K(M) SIE was temperature-independent, indicating distinct rate-limiting steps from both 15-hLO and sLO-1.     

Environmentally coupled hydrogen tunneling. Linking catalysis to dynamics.

Many biological C-H activation reactions exhibit nonclassical kinetic isotope effects (KIEs). These nonclassical KIEs are too large (kH/kD > 7) and/or exhibit unusual temperature dependence such that the Arrhenius prefactor KIEs (AH/AD) fall outside of the semiclassical range near unity. The focus of this minireview is to discuss such KIEs within the context of the environmentally coupled hydrogen tunneling model. Full tunneling models of hydrogen transfer assume that protein or solvent fluctuations generate a reactive configuration along the classical, heavy-atom coordinate, from which the hydrogen transfers via nuclear tunneling. Environmentally coupled tunneling also invokes an environmental vibration (gating) that modulates the tunneling barrier, leading to a temperature-dependent KIE. These properties directly link enzyme fluctuations to the reaction coordinate for hydrogen transfer, making a quantum view of hydrogen transfer necessarily a dynamic view of catalysis. The environmentally coupled hydrogen tunneling model leads to a range of magnitudes of KIEs, which reflect the tunneling barrier, and a range of AH/AD values, which reflect the extent to which gating modulates hydrogen transfer. Gating is the primary determinant of the temperature dependence of the KIE within this model, providing insight into the importance of this motion in modulating the reaction coordinate. The potential use of variable temperature KIEs as a direct probe of coupling between environmental dynamics and the reaction coordinate is described. The evolution from application of a tunneling correction to a full tunneling model in enzymatic H transfer reactions is discussed in the context of a thermophilic alcohol dehydrogenase and soybean lipoxygenase-1.     

Deuterium isotope effects during carbon-hydrogen bond cleavage by trimethylamine dehydrogenase. Implications for mechanism and vibrationally assisted hydrogen tunneling in wild-type and mutant enzymes

His-172 and Tyr-169 are components of a triad in the active site of trimethylamine dehydrogenase (TMADH) comprising Asp-267, His-172, and Tyr-169. Stopped-flow kinetic studies with trimethylamine as substrate have indicated that mutation of His-172 to Gln reduces the limiting rate constant for flavin reduction approximately 10-fold (Basran, J., Sutcliffe, M. J., Hille, R., and Scrutton, N. S. (1999) Biochem. J. 341, 307-314). A kinetic isotope effect (KIE = k(H)/k(D)) accompanies flavin reduction by H172Q TMADH, the magnitude of which varies significantly with solution pH. With trimethylamine, flavin reduction by H172Q TMADH is controlled by a single macroscopic ionization (pK(a) = 6.8 +/- 0.1). This ionization is perturbed (pK(a) = 7.4 +/- 0.1) in reactions with perdeuterated trimethylamine and is responsible for the apparent variation in the KIE with solution pH. At pH 9.5, where the functional group controlling flavin reduction is fully ionized, the KIE is independent of temperature in the range 277-297 K, consistent with vibrationally assisted hydrogen tunneling during breakage of the substrate C-H bond. Y169F TMADH is approximately 4-fold more compromised than H172Q TMADH for hydrogen transfer, which occurs non-classically. Studies with Y169F TMADH suggest partial thermal excitation of substrate prior to hydrogen tunneling by a vibrationally assisted mechanism. Our studies illustrate the varied effects of compromising mutations on tunneling regimes in enzyme molecules.     

Roles of the conserved aspartate and arginine in the catalytic mechanism of an archaeal beta-class carbonic anhydrase

The roles of an aspartate and an arginine, which are completely conserved in the active sites of beta-class carbonic anhydrases, were investigated by steady-state kinetic analyses of replacement variants of the beta-class enzyme (Cab) from the archaeon Methanobacterium thermoautotrophicum. Previous kinetic analyses of wild-type Cab indicated a two-step zinc-hydroxide mechanism of catalysis in which the k(cat)/K(m) value depends only on the rate constants for the CO(2) hydration step, whereas k(cat) also depends on rate constants from the proton transfer step (K. S. Smith, N. J. Cosper, C. Stalhandske, R. A. Scott, and J. G. Ferry, J. Bacteriol. 182:6605-6613, 2000). The recently solved crystal structure of Cab shows the presence of a buffer molecule within hydrogen bonding distance of Asp-34, implying a role for this residue in the proton transport step (P. Strop, K. S. Smith, T. M. Iverson, J. G. Ferry, and D. C. Rees, J. Biol. Chem. 276:10299-10305, 2001). The k(cat)/K(m) values of Asp-34 variants were decreased relative to those of the wild type, although not to an extent which supports an essential role for this residue in the CO(2) hydration step. Parallel decreases in k(cat) and k(cat)/K(m) values for the variants precluded any conclusions regarding a role for Asp-34 in the proton transfer step; however, the k(cat) of the D34A variant was chemically rescued by replacement of 2-(N-morpholino)propanesulfonic acid buffer with imidazole at pH 7.2, supporting a role for the conserved aspartate in the proton transfer step. The crystal structure of Cab also shows Arg-36 with two hydrogen bonds to Asp-34. Arg-36 variants had both k(cat) and k(cat)/K(m) values that were decreased at least 250-fold relative to those of the wild type, establishing an essential function for this residue. Imidazole was unable to rescue the k(cat) of the R36A variant; however, partial rescue of the kinetic parameter was obtained with guanidine-HCl indicating that the guanido group of this residue is important.     

Role of active site histidines in the two half-reactions of the aryl-alcohol oxidase catalytic cycle

The crystal structure of aryl-alcohol oxidase (AAO), a flavoenzyme involved in lignin degradation, reveals two active-site histidines, whose role in the two enzyme half-reactions was investigated. The redox state of flavin during turnover of the variants obtained show a stronger histidine involvement in the reductive than in the oxidative half-reaction. This was confirmed by the k(cat)/K(m(Al)) and reduction constants that are 2-3 orders of magnitude decreased for the His546 variants and up to 5 orders for the His502 variants, while the corresponding O(2) constants only decreased up to 1 order of magnitude. These results confirm His502 as the catalytic base in the AAO reductive half-reaction. The solvent kinetic isotope effect (KIE) revealed that hydroxyl proton abstraction is partially limiting the reaction, while the α-deuterated alcohol KIE showed a stereoselective hydride transfer. Concerning the oxidative half-reaction, directed mutagenesis and computational simulations indicate that only His502 is involved. Quantum mechanical/molecular mechanical (QM/MM) reveals an initial partial electron transfer from the reduced FADH(-) to O(2), without formation of a flavin-hydroperoxide intermediate. Reaction follows with a nearly barrierless His502H(+) proton transfer that decreases the triplet/singlet gap. Spin inversion and second electron transfer, concomitant with a slower proton transfer from flavin N5, yields H(2)O(2). No solvent KIE was found for O(2) reduction confirming that the His502 proton transfer does not limit the oxidative half-reaction. However, the small KIE on k(cat)/K(m(Ox)), during steady-state oxidation of α-deuterated alcohol, suggests that the second proton transfer from N5H is partially limiting, as predicted by the QM/MM simulations.     

Reaction of vascular adhesion protein-1 (VAP-1) with primary amines: mechanistic insights from isotope effects and quantitative structure-activity relationships

Human vascular adhesion protein-1 (VAP-1) is an endothelial copper-dependent amine oxidase involved in the recruitment and extravasation of leukocytes at sites of inflammation. VAP-1 is an important therapeutic target for several pathological conditions. We expressed soluble VAP-1 in HEK293 EBNA1 cells at levels suitable for detailed mechanistic studies with model substrates. Using the model substrate benzylamine, we analyzed the steady-state kinetic parameters of VAP-1 as a function of solution pH. We found two macroscopic pK(a) values that defined a bell-shaped plot of turnover number k(cat,app) as a function of pH, representing ionizable groups in the enzyme-substrate complex. The dependence of (k(cat)/K(m))(app) on pH revealed a single pK(a) value (～9) that we assigned to ionization of the amine group in free benzylamine substrate. A kinetic isotope effect (KIE) of 6 to 7.6 on (k(cat)/K(m))(app) over the pH range of 6 to 10 was observed with d(2)-benzylamine. Over the same pH range, the KIE on k(cat) was found to be close to unity. The unusual KIE values on (k(cat)/K(m))(app) were rationalized using a mechanistic scheme that includes the possibility of multiple isotopically sensitive steps. We also report the analysis of quantitative structure-activity relationships (QSAR) using para-substituted protiated and deuterated phenylethylamines. With phenylethylamines we observed a large KIE on k(cat,app) (8.01 ± 0.28 with phenylethylamine), indicating that C-H bond breakage is limiting for 2,4,5-trihydroxyphenylalanine quinone reduction. Poor correlations were observed between steady-state rate constants and QSAR parameters. We show the importance of combining KIE, QSAR, and structural studies to gain insight into the complexity of the VAP-1 steady-state mechanism.     

Vibrationally enhanced hydrogen tunneling in the Escherichia coli thymidylate synthase catalyzed reaction

The enzyme thymidylate synthase (TS) catalyzes a complex reaction that involves forming and breaking at least six covalent bonds. The physical nature of the hydride transfer step in this complex reaction cascade has been studied by means of isotope effects and their temperature dependence. Competitive kinetic isotope effects (KIEs) on the second-order rate constant (V/K) were measured over a temperature range of 5-45 degrees C. The observed H/T ((T)V/K(H)) and D/T ((T)V/K(D)) KIEs were used to calculate the intrinsic KIEs throughout the temperature range. The Swain-Schaad relationships between the H/T and D/T V/K KIEs revealed that the hydride transfer step is the rate-determining step at the physiological temperature of Escherichia coli (20-30 degrees C) but is only partly rate-determining at elevated and reduced temperatures. H/D KIE on the first-order rate constant k(cat) ((D)k = 3.72) has been previously reported [Spencer et al. (1997) Biochemistry 36, 4212-4222]. Additionally, the Swain-Schaad relationships between that (D)k and the V/K KIEs reported here suggested that at 20 degrees C the hydride transfer step is the rate-determining step for both rate constants. Intrinsic KIEs were calculated here and were found to be virtually temperature independent (DeltaE(a) = 0 within experimental error). The isotope effects on the preexponential Arrhenius factors for the intrinsic KIEs were A(H)/A(T) = 6.8 +/- 2.8 and A(D)/A(T) = 1.9 +/- 0.25. Both effects are significantly above the semiclassical (no-tunneling) predicted values and indicate a contribution of quantum mechanical tunneling to this hydride transfer reaction. Tunneling correction to transition state theory would predict that these isotope effects on activation parameters result from no energy of activation for all isotopes. Yet, initial velocity measurements over the same temperature range indicate cofactor inhibition and result in significant activation energy on k(cat) (4.0 +/- 0.1 kcal/mol). Taken together, the temperature-independent KIEs, the large isotope effects on the preexponential Arrhenius factors, and a significant energy of activation all suggest vibrationally enhanced hydride tunneling in the TS-catalyzed reaction.     

Probes of hydrogen tunneling with horse liver alcohol dehydrogenase at subzero temperatures

The temperature dependence of steady-state kinetics has been studied with horse liver alcohol dehydrogenase (HLADH) using protonated and deuterated benzyl alcohol as substrates in methanol/water mixtures between +3 and -50 degrees C. Additionally, the competitive isotope effects, k(H)/k(T) and k(D)/k(T), were measured. The studies indicate increasing kinetic complexity for wild-type HLADH at subzero temperatures. Consistent with earlier findings at 25 degrees C [Bahnson et al. (1993) Biochemistry 31, 5503], the F93W mutant shows much less kinetic complexity than the wild-type enzyme between 3 and -35 degrees C. An analysis of noncompetitive deuterium isotope effects and competitive tritium isotope effects leads to the conclusion that the reaction of F93W involves substantial hydrogen tunneling down to -35 degrees C. The effect of methanol on kinetic properties for the F93W mutant was analyzed, showing a dependence of competitive KIEs on the NAD(+) concentration. This indicates a more random bi--bi kinetic mechanism, in comparison to an ordered bi-bi kinetic mechanism in water. Although MeOH also affects the magnitude of the reaction rates and, to some extent, the observed KIEs, the ratio of ln k(H)/k(T) to ln k(D)/k(T) for primary isotope effects has not changed in methanol, and we conclude little or no change in kinetic complexity. Importantly, the degree of tunneling, as shown from the relationship between the secondary k(H)/k(T) and k(D)/k(T) values, is the same in water and MeOH/water mixtures, implicating similar trajectories for H transfer in both solvents. In a recent study of a thermophilic alcohol dehydrogenase [Kohen et al. (1999) Nature 399, 496], it was shown that decreases in temperatures below a transition temperature lead to decreased tunneling. This arises because of a change in protein dynamics below a break point in enzyme activity [Kohen et al. (2000) J. Am. Chem. Soc. 122, 10738-10739]. For the mesophilic HLADH described herein, an opposite trend is observed in which tunneling increases at subzero temperatures. These differences are attributed to inherent differences in tunneling probabilities between 0 and 100 degrees C vs subzero temperatures, as opposed to fundamental differences in protein structure for enzymes from mesophilic vs thermophilic sources. We propose that future investigations of the relationship between protein flexibility and hydrogen tunneling are best approached using enzymes from thermophilic sources.     

Studies of the enzymic mechanism of Candida tenuis xylose reductase (AKR 2B5): X-ray structure and catalytic reaction profile for the H113A mutant

Xylose reductase from the yeast Candida tenuis (CtXR) is a family 2 member of the aldo-keto reductase (AKR) superfamily of proteins and enzymes. Active site His-113 is conserved among AKRs, but a unified mechanism of how it affects catalytic activity is outstanding. We have replaced His-113 by alanine using site-directed mutagenesis, determined a 2.2 A structure of H113A mutant bound to NADP(+), and compared catalytic reaction profiles of NADH-dependent reduction of different aldehydes catalyzed by the wild type and the mutant. Deuterium kinetic isotope effects (KIEs) on k(cat) and k(cat)/K(m xylose) show that, relative to the wild type, the hydride transfer rate constant (k(7) approximately 0.16 s(-1)) has decreased about 1000-fold in H113A whereas xylose binding was not strongly affected. No solvent isotope effect was seen on k(cat) and k(cat)/K(m xylose) for H113A, suggesting that proton transfer has not become rate-limiting as a result of the mutation. The pH profiles of log(k(cat)/K(m xylose)) for the wild type and H113A decreased above apparent pK(a) values of 8.85 and 7.63, respectively. The DeltapK(a) of -1.2 pH units likely reflects a proximally disruptive character of the mutation, affecting the position of Asp-50. A steady-state kinetic analysis for H113A-catalyzed reduction of a homologous series of meta-substituted benzaldehyde derivatives was carried out, and quantitative structure-reactivity correlations were used to factor the observed kinetic substituent effect on k(cat) and k(cat)/K(m aldehyde) into an electronic effect and bonding effects (which are lacking in the wild type). Using the Hammett sigma scale, electronic parameter coefficients (rho) of +0.64 (k(cat)) and +0.78 (k(cat)/K(m aldehyde)) were calculated and clearly differ from rho(k(cat)/K(aldehyde)) and rho(k(cat)) values of +1.67 and approximately 0.0, respectively, for the wild-type enzyme. Hydride transfer rate constants of H113A, calculated from kinetic parameters and KIE data, display a substituent dependence not seen in the corresponding wild-type enzyme rate constants. An enzymic mechanism is proposed in which His-113, through a hydrogen bond from Nepsilon2 to aldehyde O1, assists in catalysis by optimizing the C=O bond charge separation and orbital alignment in the ternary complex.     

Proposal for a hydrogen bond network in the active site of the prototypic gamma-class carbonic anhydrase

The crystal structure of Cam, the prototypic gamma-class carbonic anhydrase, reveals active site residues Gln75, Asn73, and Asn 202 previously hypothesized to participate in catalysis. These potential roles were investigated for the first time by kinetic analyses of site-specific replacement variants of the zinc and cobalt forms of Cam. Gln75 replacement variants showed large decreases in k(cat)/K(m) relative to wild-type. Further, the Gln75 variants showed a loss of the pK(a) in pH versus k(cat)/K(m) profiles previously attributed to ionization of the metal-bound water yielding the hydroxyl group attacking CO(2). These results support the previously proposed role for Gln75 in hydrogen bonding with the catalytic hydroxyl orienting it for attack on CO(2). Kinetic analyses of Asn73 variants were consistent with a role in hydrogen bonding with Gln75 to position it for optimal interaction with the catalytic hydroxyl. Kinetic analyses of Asn202 variants showed substantial decreases in k(cat)/K(m) relative to the wild-type enzyme supporting the previously hypothesized role in polarizing CO(2) and facilitating attack from the metal-bound hydroxyl. On the basis of results presented here, and previously reported structural analyses, we present a catalytic mechanism involving Gln75, Asn73, and Asn202 that also suggests a role for Glu62 not previously recognized. Finally, the results suggest that the gamma-, beta-, and alpha-class carbonic anhydrases each independently evolved variations of a fundamental hydrogen bond network essential for catalysis.     

Nature of hydrogen transfer in soybean lipoxygenase 1: separation of primary and secondary isotope effects

Previous measurements of the kinetics of oxidation of linoleic acid by soybean lipoxygenase 1 have indicated very large deuterium isotope effects, but have not been able to distinguish the primary isotope effect from the alpha-secondary effect. To address this question, singly deuterated linoleic acid was prepared, and enantiomerically resolved using the enzyme itself. Noncompetitive measurements of the primary deuterium isotope effect give a value of ca. 40 which is temperature-independent. The enthalpy of activation is low and isotope-independent, and there is a large isotope effect on the Arrhenius prefactor. A very large apparent secondary isotope effect (ca. 2.1) is measured with deuterium in the primary position, but a greatly reduced value (1.1) is observed with protium in the primary position. Mutagenesis of the active site leads to a significant reduction in k(cat) and perturbed isotope effects, in particular, a secondary effect of 5.6 when deuterium is in the primary position. The anomalous secondary isotope effects are shown to arise from imperfect stereoselectivity of hydrogen abstraction which, for the mutant, is attributed to a combination of inverse substrate binding and increased flexibility at the reactive carbon. After correction, a very large primary (76-84) and small secondary (1.1-1.2) kinetic isotope effects are calculated for both mutant and wild-type enzymes. The weight of the evidence is taken to favor hydrogen tunneling as the primary mechanism of hydrogen transfer.     

Probing functional perfection in substructures of ribonuclease T1: double combinatorial random mutagenesis involving Asn43, Asn44, and Glu46 in the guanine binding loop

Combinatorial random mutageneses involving either Asn43 with Asn44 (set 1) or Glu46 with an adjacent insertion (set 2) were undertaken to explore the functional perfection of the guanine recognition loop of ribonuclease T(1) (RNase T(1)). Four hundred unique recombinants were screened in each set for their ability to enhance enzyme catalysis of RNA cleavage. After a thorough selection procedure, only six variants were found that were either as active or more active than wild type which included substitutions of Asn43 by Gly, His, Leu, or Thr, an unplanned Tyr45Ser substitution and Glu46Pro with an adjacent Glu47 insertion. Asn43His-RNase T(1) has the same loop sequence as that for RNases Pb(1) and Fl(2). None of the most active mutants were single substitutions at Asn44 or double substitutions at Asn43 and Asn44. A total of 13 variants were purified, and these were subjected to kinetic analysis using RNA, GpC, and ApC as substrates. Modestly enhanced activities with GpC and RNA involved both k(cat) and K(M) effects. Mutants having low activity with GpC had proportionately even lower relative activity with RNA. Asn43Gly-RNase T(1) and all five of the purified mutants in set 2 exhibited similar values of k(cat)/K(M) for ApC which were the highest observed and about 10-fold that for wild type. The specificity ratio [(k(cat)/K(M))(GpC)/(k(cat)/K(M))(ApC)] varied over 30 000-fold including a 10-fold increase [Asn43His variant; mainly due to a low (k(cat)/K(M))(ApC)] and a 3000-fold decrease (Glu46Ser/(insert)Gly47 variant; mainly due to a low (k(cat)/K(M))(GpC)) as compared with wild type. It is interesting that k(cat) (GpC) for the Tyr45Ser variant was almost 4-fold greater than for wild type and that Pro46/(insert)Glu47 RNase T(1) is 70-fold more active than the permuted variant (insert)Pro47-RNase T(1) which has a conserved Glu46. In any event, the observation that only 6 out of 800 variants surveyed had wild-type activity supports the view that functional perfection of the guanine recognition loop of RNase T(1) has been achieved.     

Discrimination of cognate and noncognate substrates at the active site of class I lysyl-tRNA synthetase

The aminoacyl-tRNA synthetases are divided into two unrelated structural classes, with lysyl-tRNA synthetase (LysRS) being the only enzyme represented in both classes. On the basis of the structure of l-lysine complexed with Pyrococcus horikoshii class I LysRS (LysRS1) and homology to glutamyl-tRNA synthetase (GluRS), residues implicated in amino acid recognition and noncognate substrate discrimination were systematically replaced in Borrelia burgdorferi LysRS1. The catalytic efficiency of steady-state aminoacylation (k(cat)/K(M)) with lysine by LysRS1 variants fell by 1-4 orders of magnitude compared to that of the wild type. Disruption of putative hydrogen bonding interactions through replacement of G29, T31, and Y269 caused up to 1500-fold reductions in k(cat)/K(M), similar to changes previously observed for comparable variants of class II LysRS (LysRS2). Replacements of W220 and H242, both of which are implicated in hydrophobic interactions with the side chain of lysine, resulted in more dramatic changes with up to 40000-fold reductions in k(cat)/K(M) observed. This indicates that the more compact LysRS1 active site employs both electrostatic and hydrophobic interactions during lysine discrimination, explaining the ability of LysRS1 to discriminate against noncognate substrates accepted by LysRS2. Several of the LysRS1 variants were found to be more specific than the wild type with respect to noncognate amino acid recognition but less efficient in cognate aminoacylation. This indicates that LysRS1 compromises between efficient catalysis and substrate discrimination, in contrast to LysRS2 which is considerably more effective in catalysis but is less specific than its class I counterpart.     

Loop Residues and Catalysis in OMP Synthase

Residue-to-alanine mutations and a two-amino acid deletion have been made in the highly conserved catalytic loop (residues 100-109) of Salmonella typhimurium OMP synthase (orotate phosphoribosyltransferase, EC 2.4.2.10). As described previously, the K103A mutant enzyme exhibited a 10(4)-fold decrease in k(cat)/K(M) for PRPP; the K100A enzyme suffered a 50-fold decrease. Alanine mutations at His105 and Glu107 produced 40- and 7-fold decreases in k(cat)/K(M), respectively, and E101A, D104A, and G106A were slightly faster than the wild-type (WT) in terms of k(cat), with minor effects on k(cat)/K(M). Equilibrium binding of OMP or PRPP in binary complexes was affected little by loop mutation, suggesting that the energetics of ground-state binding have little contribution from the catalytic loop, or that a favorable binding energy is offset by costs of loop reorganization. Pre-steady-state kinetics for mutants showed that K103A and E107A had lost the burst of product formation in each direction that indicated rapid on-enzyme chemistry for WT, but that the burst was retained by H105A. Δ102Δ106, a loop-shortened enzyme with Ala102 and Gly106 deleted, showed a 10(4)-fold reduction of k(cat) but almost unaltered K(D) values for all four substrate molecules. The 20% (i.e., 1.20) intrinsic [1'-(3)H]OMP kinetic isotope effect (KIE) for WT is masked because of high forward and reverse commitment factors. K103A failed to express intrinsic KIEs fully (1.095 ± 0.013). In contrast, H105A, which has a smaller catalytic lesion, gave a [1'-(3)H]OMP KIE of 1.21 ± 0.0005, and E107A (1.179 ± 0.0049) also gave high values. These results are interpreted in the context of the X-ray structure of the complete substrate complex for the enzyme [Grubmeyer, C., Hansen, M. R., Fedorov, A. A., and Almo, S. C. (2012) Biochemistry 51 (preceding paper in this issue, DOI 10.1021/bi300083p )]. The full expression of KIEs by H105A and E107A may result from a less secure closure of the catalytic loop. The lower level of expression of the KIE by K103A suggests that in these mutant proteins the major barrier to catalysis is successful closure of the catalytic loop, which when closed, produces rapid and reversible catalysis.     

Role of arginine 59 in the gamma-class carbonic anhydrases.

The functional role of the highly conserved active site Arg 59 in the prototype of the gamma-class carbonic anhydrase Cam (carbonic anhydrase from Methanosarcina thermophila) was investigated. Variants (R59A, -C, -E, -H, -K, -M, and -Q) were prepared by site-directed mutagenesis and characterized by size exclusion chromatography (SEC), circular dichroism (CD) spectroscopy, and stopped-flow kinetic analyses. CD spectra indicated similar secondary structures for the wild type and the R59A and -K variants, independent of nondenaturing concentrations of guanidine hydrochloride (GdnHCl). SEC indicated that all variants purified as homotrimers like the wild type. SEC also revealed that the R59A and -K variants unfolded at > or = 1.5 M GdnHCl, compared to 3.0 M GdnHCl for the wild type. These results indicate that Arg 59 contributes to the thermodynamic stability of the Cam trimer. The R59K variant had k(cat) and k(cat)/K(m) values that were 8 and 5% of the wild-type values, respectively, while all other variants had k(cat) and k(cat)/K(m) values 10-100-fold lower than those of the wild type. The R59A, -C, -E, -M, and -Q variants exhibited 4-63-fold increases in k(cat) and 9-120-fold increases in k(cat)/K(m) upon addition of 100 mM GdnHCl, with the largest increases observed for the R59A variant, which was comparable to the R59K variant. The kinetic results indicate that a positive charge at position 59 is essential for the CO(2) hydration step of the overall catalytic mechanism.     

A rearranging ligand enables allosteric control of catalytic activity in copper-containing nitrite reductase

In Cu-containing nitrite reductase from Alcaligenes faecalis S-6 the axial methionine ligand of the type-1 site was replaced (M150G) to make the copper ion accessible to external ligands that might affect the enzyme's catalytic activity. The type-1 site optical spectrum of M150G (A(460)/A(600)=0.71) differs significantly from that of the native nitrite reductase (A(460)/A(600)=1.3). The midpoint potential of the type-1 site of nitrite reductase M150G (E(M)=312(+/-5)mV versus hydrogen) is higher than that of the native enzyme (E(M)=213(+/-5)mV). M150G has a lower catalytic activity (k(cat)=133(+/-6)s(-1)) than the wild-type nitrite reductase (k(cat)=416(+/-10)s(-1)). The binding of external ligands to M150G restores spectral properties, midpoint potential (E(M)<225mV), and catalytic activity (k(cat)=374(+/-28)s(-1)). Also the M150H (A(460)/A(600)=7.7, E(M)=104(+/-5)mV, k(cat)=0.099(+/-0.006)s(-1)) and M150T (A(460)/A(600)=0.085, E(M)=340(+/-5)mV, k(cat)=126(+/-2)s(-1)) variants were characterized. Crystal structures show that the ligands act as allosteric effectors by displacing Met62, which moves to bind to the Cu in the position emptied by the M150G mutation. The reconstituted type-1 site has an otherwise unaltered geometry. The observation that removal of an endogenous ligand can introduce allosteric control in a redox enzyme suggests potential for structural and functional flexibility of copper-containing redox sites.     

Thermostability of native and pegylated Myceliophthora thermophila laccase in aqueous and mixed solvents

A commercial preparation of laccase (EC 1.10.3.2), cloned from Myceliophthora thermophila and expressed in Aspergillus oryzae (MtL), was purified and modified by conjugation with poly(ethylene glycol) (M(r) = 5000) and is labeled PEG-MtL. Native enzyme was found to have a molecular mass of 80 kDa, as determined by gel filtration, and 110 kDa, by SDS-PAGE. The oxidative dimerization of 2,6-dimethoxyphenol (DMP) to produce the corresponding dibenzoquinone was catalyzed by MtL in a manner comparable to that for a diffusion-controlled reaction (k(cat)/K(M) approximately = 10(8) M(-)(1) s(-)(1) and E(a) approximately = 18 kJ M(-)(1)). PEG-MtL was found, by TNBS titration, to have blocked 54% of lysine groups; its hydrodynamic and charge properties were different from those of MtL. Catalytic efficiency (k(cat)/K(M)) of PEG-MtL was similar to that of MtL with DMP as substrate; however, k(cat)/K(M) was 2-fold reduced for the reaction in which 2',2-azinobis(3-ethylbenzothiazoline-6-sulfonic acid) (ABTS) is oxidized to form a radical cation. E(a) values were similar in both enzyme preparations when assayed in buffered solutions. Far-UV CD spectra were similar for MtL and PEG-MtL and consistent with a protein rich in beta-sheet structure with negligible content of alpha-helices. A blue shift of near-UV CD spectrum for PEG-MtL as compared to MtL was consistent with the decreased polarity of the tyrosyl side chains upon PEG conjugation. Also the blue band of the copper active site was shifted from lambda approximately 610 nm (MtL) to lambda approximately 575 nm (PEG-MtL). Scanning microcalorimetry showed small denaturation enthalpies (6.3 and 7.5 J g(-)(1) for MtL and PEG-MtL, respectively), indicating the high stability of the beta-sheet folding pattern of laccases. However, PEG-MtL proved to be more stable, its half-denaturation temperature being 2 degrees C higher than that of MtL. In 30% alcohol, pegylated laccase showed slower enzyme-activity decay rates than the unmodified enzyme; this behavior was caused by a decrease in the activation entropy of the denaturation reaction. Results can be explained by entropic stabilization by PEG conjugation because of the restricted motion of some surface amino acid side chains, which results in a more stable active site.     

A structure-function study of a proton transport pathway in the gamma-class carbonic anhydrase from Methanosarcina thermophila

Four glutamate residues in the prototypic gamma-class carbonic anhydrase from Methanosarcina thermophila (Cam) were characterized by site-directed mutagenesis and chemical rescue studies. Alanine substitution indicated that an external loop residue, Glu 84, and an internal active site residue, Glu 62, are both important for CO(2) hydration activity. Two other external loop residues, Glu 88 and Glu 89, are less important for enzyme function. The two E84D and -H variants exhibited significant activity relative to wild-type activity in pH 7.5 MOPS buffer, suggesting that the original glutamate residue could be substituted with other ionizable residues with similar pK(a) values. The E84A, -C, -K, -Q, -S, and -Y variants exhibited large decreases in k(cat) values in pH 7.5 MOPS buffer, but only exhibited small changes in k(cat)/K(m). These same six variants were all chemically rescued by pH 7.5 imidazole buffer, with 23-46-fold increases in the apparent k(cat). These results are consistent with Glu 84 functioning as a proton shuttle residue. The E62D variant exhibited a 3-fold decrease in k(cat) and a 2-fold decrease in k(cat)/K(m) relative to those of the wild type in pH 7.5 MOPS buffer, while other substitutions (E62A, -C, -H, -Q, -T, and -Y) resulted in much larger decreases in both k(cat) and k(cat)/K(m). Imidazole did not significantly increase the k(cat) values and slightly decreased the k(cat)/K(m) values of most of the Glu 62 variants. These results indicate a primary preference for a carboxylate group at position 62, and support a proposed catalytic role for residue Glu 62 in the CO(2) hydration step, but do not definitively establish its role in the proton transport step.     

Contribution of active-site residues to the function of onconase, a ribonuclease with antitumoral activity

Onconase (ONC), a homologue of ribonuclease A (RNase A), is in clinical trials for the treatment of cancer. ONC possesses a conserved active-site catalytic triad, which is composed of His10, Lys31, and His97. The three-dimensional structure of ONC suggests that two additional residues, Lys9 and an N-terminal lactam formed from a glutamine residue (Pca1), could also contribute to catalysis. To determine the role of Pca1, Lys9, and Lys31 in the function of ONC, site-directed mutagenesis was used to replace each with alanine. Values of k(cat)/K(M) for the variants were determined with a novel fluorogenic substrate, which was designed to match the nucleobase specificity of ONC and gives the highest known k(cat)/K(M) value for the enzyme. The K9A and K31A variants display 10(3)-fold lower k(cat)/K(M) values than the wild-type enzyme, and a K9A/K31A double variant suffers a >10(4)-fold decrease in catalytic activity. In addition, replacing Lys9 or Lys31 eliminates the antitumoral activity of ONC. The side chains of Pca1 and Lys9 form a hydrogen bond in crystalline ONC. Replacing Pca1 with an alanine residue lowers the catalytic activity of ONC by 20-fold. Yet, replacing Pca1 in the K9A variant enzyme does not further reduce catalytic activity, revealing that the function of the N-terminal pyroglutamate residue is to secure Lys9. The thermodynamic cycle derived from k(cat)/K(M) values indicates that the Pca1...Lys9 hydrogen bond contributes 2.0 kcal/mol to the stabilization of the rate-limiting transition state during catalysis. Finally, binding isotherms with a substrate analogue indicate that Lys9 and Lys31 contribute little to substrate binding and that the low intrinsic catalytic activity of ONC originates largely from the low affinity of the enzyme for its substrate. These findings could assist the further development of ONC as a cancer chemotherapeutic.