Contributions Of Active Site Residues Research Articles

Structural and biochemical analysis of phosphoethanolamine methyltransferase from the pine wilt nematode Bursaphelenchus xylophilus

In free‐living and parasitic nematodes, the methylation of phosphoethanolamine to phosphocholine provides a key metabolite to sustain phospholipid biosynthesis for growth and development. Because the phosphoethanolamine methyltransferases (PMT) of nematodes are essential for normal growth and development, these enzymes are potential targets of inhibitor design. The pine wilt nematode (Bursaphelenchus xylophilus) causes extensive damage to trees used for lumber and paper in Asia. As a first step toward testing BxPMT1 as a potential nematicide target, we determined the 2.05 Å resolution x‐ray crystal structure of the enzyme as a dead‐end complex with phosphoethanolamine and S‐adenosylhomocysteine. The three‐dimensional structure of BxPMT1 served as a template for site‐directed mutagenesis to probe the contribution of active site residues to catalysis and phosphoethanolamine binding using steady‐state kinetic analysis. Biochemical analysis of the mutants identifies key residues on the β1d‐α6 loop (W123F, M126I, and Y127F) and β1e‐α7 loop (S155A, S160A, H170A, T178V, and Y180F) that form the phosphobase binding site and suggest that Tyr127 facilitates the methylation reaction in BxPMT1.

In silico evaluation of NO donor heterocyclic vasodilators as SARS-CoV-2 Mpro protein inhibitor

Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) which causes COVID-19 disease has been exponentially increasing throughout the world. The mortality rate is increasing gradually as effective treatment is unavailable to date. In silico based screening for novel testable hypotheses on SARS-CoV-2 Mpro protein to discover the potential lead drug candidate is an emerging area along with the discovery of a vaccine. Administration of NO-releasing agents, NO inducers or the NO gas itself may be useful as therapeutics in the treatment of SARS-CoV-2. In the present study, a 3D structure of SARS-CoV-2 Mpro protein was used for the rational setting of inhibitors to the binding pocket of enzyme which proposed that phenyl furoxan derivative gets efficiently dock in the target pocket. Molecular docking and molecular dynamics simulations helped to investigate possible effective inhibitor candidates bound to SARS-CoV-2 Mpro substrate binding pocket. Molecular mechanics Poisson-Boltzmann surface area (MM/PBSA) calculations revealed energetic contributions of active site residues of Mpro in binding with most stable proposed NO donor heterocyclic vasodilator inhibitor molecules. Furthermore, principal component analysis (PCA) showed that the NO donor heterocyclic inhibitor molecules 14, 16, 18 and 19 was strongly bound to catalytic core of SARS-CoV-2 Mpro protein, limiting its movement to form stable complex as like control. Thus, overall in silico investigations revealed that 5-oxopiperazine-2-carboxylic acid coupled furoxan derivatives was found to be key pharmacophore in drug design for the treatment of SARS-CoV-2, a global pandemic disease with a dual mechanism of action as NO donor and a worthwhile ligand to act as SARS-CoV-2 Mpro protein inhibitor. Communicated by Ramaswamy H. Sarma

Investigation of glutathione as a natural antioxidant and multitarget inhibitor for Alzheimer’s disease: Insights from molecular simulations

Glutathione (GSH) is the most plentiful nonprotein thiol-containing substance in all kinds of cells, and plays crucial roles in the antioxidant defense system and maintaining redox homeostasis in neurons. Extensive microscopic molecular simulations with atomistic details were performed to examine GSH interaction with three proteins involved in β-Amyloid formation (BACE-1, GSK-3β, AChE). The results show that the binding of GSH to AChE is vertical and from the side of Cys residue, while for BACE-1 it was relatively horizontal. The simultaneous study of MSD and interaction energies proved that due to the stronger interaction between GSH and the active-site of BACE-1, GSH shows the least mobility in the BACE-1 containing system, also due to no binding to GSK-3β, it shows the most mobility in this system. Analyzing the radius of gyration and the molecule's inner angle showed that GSH adopts folded and extended structures in interaction with BACE-1 and AChE, respectively. Estimating binding free energy using MMPBSA calculations proved that GSH binding to BACE-1 (-16.98 kJ/mol) is thermodynamically more favorable than AChE (-8.27 kJ/mol). Contributions of active-site residues in binding free energies revealed that GSH binding into BACE-1 and AChE was controlled more dynamically and thermodynamically, respectively. Cavity size analysis showed that the accessible volume of cavities for GSH in BACE-1 and AChE is sufficient for binding. Furthermore, “adjacent pocket” and “breathing motion” types were suggested for AChE and BACE-1, respectively.

Structure-Based Design of Tetrahydroisoquinoline-Based Hydroxamic Acid Derivatives Inhibiting Human Histone Deacetylase 8

We have designed new human histone deacetylase 8 (HDAC8) inhibitors using structure-based molecular design. 3D models of HDAC8–inhibitor complexes were prepared by in situ modification of the crystal structure of HDAC8 co-crystallized with the hydroxamic acid suberoylanilide (SAHA) and a training set (TS) of tetrahydroisoquinoline-based hydroxamic acid derivatives (DAHTs) with known inhibitory potencies. A QSAR model was elaborated for the TS yielding a linear correlation between the computed Gibbs free energies (GFE) of HDAC8–DAHTs complexation (∆∆Gcom) and observed half-maximal enzyme inhibitory concentrations (IC50exp). From this QSAR model a 3D-QSAR pharmacophore (PH4) was generated. Structural information derived from the 3D model and breakdown of computed HDAC8–DAHTs interaction energies up to individual active site residue contributions helped us to design new more potent HDAC8 inhibitors. We obtained a reasonable agreement ∆∆Gcom and values: (pIC50exp = – 0.0376 × ∆∆Gcom + 7.4605, R2 = 0.89). Similar agreement was established for the PH4 model (pIC50exp = 0.8769 × + 0.7854, R2 = 0.87). A comparative analysis of the contributions of active site residues guided the choice of fragments used in designing a virtual combinatorial library (VCL) of DAHT analogs. The VCL of more than 17 thousand DAHTs was screened by the PH4 and furnished 229 new DAHTs. The best-designed analog displayed predicted inhibitory potency up to 110 times higher than that of DAHT1 (IC50exp = 0.047 µM). Predicted pharmacokinetic profiles of the new analogs were compared to current per oral anticancer compounds. This computational approach, which combines molecular modelling, pharmacophore generation, analysis of HDAC8–DAHTs interaction energies and virtual screening of a combinatorial library of DAHTs resulted in a set of proposed new HDAC8 inhibitors. It can thus direct medicinal chemists in their search for new anticancer agents.

Structural and biochemical analysis of phosphoethanolamine methyltransferase from the pine wilt nematode Bursaphelenchus xylophilus

In free-living and parasitic nematodes, the methylation of phosphoethanolamine to phosphocholine provides a key metabolite to sustain phospholipid biosynthesis for growth and development. Because the phosphoethanolamine methyltransferases (PMT) of nematodes are essential for normal growth and development, these enzymes are potential targets of inhibitor design. The pine wilt nematode (Bursaphelenchus xylophilus) causes extensive damage to trees used for lumber and paper in Asia. As a first step toward testing BxPMT1 as a potential nematicide target, we determined the 2.05 Å resolution x-ray crystal structure of the enzyme as a dead-end complex with phosphoethanolamine and S-adenosylhomocysteine. The three-dimensional structure of BxPMT1 served as a template for site-directed mutagenesis to probe the contribution of active site residues to catalysis and phosphoethanolamine binding using steady-state kinetic analysis. Biochemical analysis of the mutants identifies key residues on the β1d-α6 loop (W123F, M126I, and Y127F) and β1e-α7 loop (S155A, S160A, H170A, T178V, and Y180F) that form the phosphobase binding site and suggest that Tyr127 facilitates the methylation reaction in BxPMT1.

Regioselectivity of hyoscyamine 6β-hydroxylase-catalysed hydroxylation as revealed by high-resolution structural information and QM/MM calculations.

Hyoscyamine 6β-hydroxylase (H6H) is a bifunctional non-heme 2-oxoglutarate/Fe2+-dependent dioxygenase that catalyzes the two final steps in the biosynthesis of scopolamine. Based on high resolution crystal structures of H6H from Datura metel, detailed information on substrate binding was obtained that provided insights into the onset of the enzymatic process. In particular, the role of two prominent residues was revealed - Glu-116 that interacts with the tertiary amine located on the hyoscyamine tropane moiety and Tyr-326 that forms CH-π hydrogen bonds with the hyoscyamine phenyl ring. The structures were used as the basis for QM/MM calculations that provided an explanation for the regioselectivity of the hydroxylation reaction on the hyoscyamine tropane moiety (C6 vs. C7) and quantified contributions of active site residues to respective barrier heights.

Contribution of Active Site Residues to Substrate Hydrolysis by USP2: Insights into Catalysis by Ubiquitin Specific Proteases

The ubiquitin-specific protease (USP) structural class represents the largest and most diverse family of deubiquitinating enzymes (DUBs). Many USPs assume important biological roles and emerge as potential targets for therapeutic intervention. A clear understanding of USP catalytic mechanism requires a functional evaluation of the proposed key active site residues. Crystallographic data of ubiquitin aldehyde adducts of USP catalytic cores provided structural details on the catalytic triad residues, namely the conserved Cys and His, and a variable putative third residue, and inferred indirect structural roles for two other conserved residues (Asn and Asp), in stabilizing via a bridging water molecule the oxyanion of the tetrahedral intermediate (TI). We have expressed the catalytic domain of USP2 and probed by site-directed mutagenesis the role of these active site residues in the hydrolysis of peptide and isopeptide substrates, including a synthetic K48-linked diubiquitin substrate for which a label-free, mass spectrometry based assay has been developed to monitor cleavage. Hydrolysis of ubiquitin-AMC, a model substrate, was not affected by the mutations. Molecular dynamics simulations of USP2, free and complexed with the TI of a bona fide isopeptide substrate, were carried out. We found that Asn271 is structurally poised to directly stabilize the oxyanion developed in the acylation step, while being structurally supported by the adjacent absolutely conserved Asp575. Mutagenesis data functionally confirmed this structural role independent of the nature (isopeptide vs peptide) of the bond being cleaved. We also found that Asn574, structurally located as the third member of the catalytic triad, does not fulfill this role functionally. A dual supporting role is inferred from double-point mutation and structural data for the absolutely conserved residue Asp575, in oxyanion hole formation, and in maintaining the correct alignment and protonation of His557 for catalytic competency.

Divergence of Chemical Function in the Alkaline Phosphatase Superfamily: Structure and Mechanism of the P−C Bond Cleaving Enzyme Phosphonoacetate Hydrolase

Phosphonates constitute a class of natural products that mimic the properties of the more common organophosphate ester metabolite yet are not readily degraded owing to the direct linkage of the phosphorus atom to the carbon atom. Phosphonate hydrolases have evolved to allow bacteria to utilize environmental phosphonates as a source of carbon and phosphorus. The work reported in this paper examines one such enzyme, phosphonoacetate hydrolase. By using a bioinformatic approach, we circumscribed the biological range of phosphonoacetate hydrolase to a select group of bacterial species from different classes of Proteobacteria. In addition, using gene context, we identified a novel 2-aminoethylphosphonate degradation pathway in which phosphonoacetate hydrolase is a participant. The X-ray structure of phosphonoformate-bound phosphonoacetate hydrolase was determined to reveal that this enzyme is most closely related to nucleotide pyrophosphatase/diesterase, a promiscuous two-zinc ion metalloenzyme of the alkaline phosphatase enzyme superfamily. The X-ray structure and metal ion specificity tests showed that phosphonoacetate hydrolase is also a two-zinc ion metalloenzyme. By using site-directed mutagenesis and (32)P-labeling strategies, the catalytic nucleophile was shown to be Thr64. A structure-guided, site-directed mutation-based inquiry of the catalytic contributions of active site residues identified Lys126 and Lys128 as the most likely candidates for stabilization of the aci-carboxylate dianion leaving group. A catalytic mechanism is proposed which combines Lys12/Lys128 leaving group stabilization with zinc ion activation of the Thr64 nucleophile and the substrate phosphoryl group.

Probing the Molecular Basis of Substrate Specificity, Stereospecificity, and Catalysis in the Class II Pyruvate Aldolase, BphI

BphI, a pyruvate-specific class II aldolase found in the polychlorinated biphenyls (PCBs) degradation pathway, catalyzes the reversible C-C bond cleavage of (4S)-hydroxy-2-oxoacids to form pyruvate and an aldehyde. Mutations were introduced into bphI to probe the contribution of active site residues to substrate recognition and catalysis. In contrast to the wild-type enzyme that has similar specificities for acetaldehyde and propionaldehyde, the L87A variant exhibited a 40-fold preference for propionaldehyde over acetaldehyde. The specificity constant of the L89A variant in the aldol addition reaction using pentaldehyde is increased ∼50-fold, making it more catalytically efficient for pentaldehyde utilization compared to the wild-type utilization of the natural substrate, acetaldehyde. Replacement of Tyr-290 with phenylalanine or serine resulted in a loss of stereochemical control as the variants were able to utilize substrates with both R and S configurations at C4 with similar kinetic parameters. Aldol cleavage and pyruvate α-proton exchange activity were undetectable in the R16A variant, supporting the role of Arg-16 in stabilizing a pyruvate enolate intermediate. The pH dependence of the enzyme is consistent with a single deprotonation by a catalytic base with pK(a) values of approximately 7. In H20A and H20S variants, pH profiles show the dependence of enzyme activity on hydroxide concentration. On the basis of these results, a catalytic mechanism is proposed.

Contributions of active site residues to cofactor binding and catalysis of 3α-hydroxysteroid dehydrogenase/carbonyl reductase

3α-Hydroxysteroid dehydrogenase/carbonyl reductase reversely catalyzes the oxidation of androsterone with NAD+ to form androstanedione and NADH. In this study, we investigated the function of active site residues N86, Y155, and K159 in NADH binding and catalysis in the reduction of androstanedione, using site-directed mutagenesis, steady-state kinetics, fluorescence quenching, and anisotropy measurements. The N86A, Y155F, and K159A mutant enzymes decreased the catalytic constant by 37- to 220-fold and increased the dissociation constant by 3- to 75-fold, respectively. Binding of NADH with wild-type and mutant enzymes caused different levels of fluorescence resonance energy transfer, implying a different orientation of nicotinamide ring versus W173. In addition, the enzyme-bound NADH decreased the fluorescence anisotropy value in the order WT>N86A>Y155F>K159A, indicating an increase in the mobility of the bound NADH for the mutants. Data suggest that hydrogen bonding with the hydroxyl group of nicotinamide ribose by K159 and Y155 is important to maintain the orientation of NADH and contributes greatly to the transition-state binding energy to facilitate the catalysis. N86 is important for stabilizing the position of K159. Substitution of alanine for N86 has a minor effect on NADH binding through K159, resulting in a slight increase in the mobility of the bound NADH and decreases in affinity and catalytic constant.

Theoretical Determination of the Redox Potentials of NRH:Quinone Oxidoreductase 2 Using Quantum Mechanical/Molecular Mechanical Simulations

NRH:quinone oxidoreductase 2 (NQO2) is a flavoenzyme that catalyzes a one-step two-electron reduction of quinones. During this enzyme catalysis, the 7,8-dimethyl isoalloxazine (flavin) ring of the enzyme-bound cofactor, flavin adenine dinucleotide (FAD), shuttles between reduced and oxidized states as the enzyme passes through multiple cycles of binding/release of alternate substrates. These redox changes in NQO2, however, lead to unequal charge separation between the flavin ring and the active site, which must be stabilized by reorganization of the surrounding protein matrix. In this study, we have used a combined quantum mechanical/molecular mechanical method to simulate the electron and proton addition reactions of the flavin-bound NQO2. We have computed the redox potentials and pK(a)'s of the enzyme-bound flavin. The present work demonstrates that upon reduction, the NQO2 active site stabilizes the flavin anionic hydroquinone state. Simulation data has also allowed quantitative estimation of the electrostatic contributions of active site residues. Their significance in oscillatory redox transition of this flavoenzyme is discussed.

Comparative residue interaction analysis (CoRIA): a 3D-QSAR approach to explore the binding contributions of active site residues with ligands

A novel approach termed comparative residue-interaction analysis (CoRIA), emphasizing the trends and principles of QSAR in a ligand-receptor environment has been developed to analyze and predict the binding affinity of enzyme inhibitors. To test this new approach, a training set of 36 COX-2 inhibitors belonging to nine families was selected. The putative binding (bioactive) conformations of inhibitors in the COX-2 active site were searched using the program DOCK. The docked configurations were further refined by a combination of Monte Carlo and simulated annealing methods with the Affinity program. The non-bonded interaction energies of the inhibitors with the individual amino acid residues in the active site were then computed. These interaction energies, plus specific terms describing the thermodynamics of ligand-enzyme binding, were correlated to the biological activity with G/PLS. The various QSAR models obtained were validated internally by cross validation and boot strapping, and externally using a test set of 13 molecules. The QSAR models developed on the CoRIA formalism were robust with good r (2), q (2) and r (pred) (2) values. The major highlights of the method are: adaptation of the QSAR formalism in a receptor setting to answer both the type (qualitative) and the extent (quantitative) of ligand-receptor binding, and use of descriptors that account for the complete thermodynamics of the ligand-receptor binding. The CoRIA approach can be used to identify crucial interactions of inhibitors with the enzyme at the residue level, which can be gainfully exploited in optimizing the inhibitory activity of ligands. Furthermore, it can be used with advantage to guide point mutation studies. As regards the COX-2 dataset, the CoRIA approach shows that improving Coulombic interaction with Pro528 and reducing van der Waals interaction with Tyr385 will improve the binding affinity of inhibitors.

Differential contribution of active site residues in substrate recognition sites 1 and 5 to cytochrome P450 2C8 substrate selectivity and regioselectivity.

Selected active site residues in substrate recognition sites (SRS) 1 and 5 of cytochrome P450 2C8 (CYP2C8) were mutated to the corresponding amino acids present in CYP2C9 to investigate the contribution of these positions to the unique substrate selectivity and regioselectivity of CYP2C8. The effects of mutations, singly and in combination, were assessed from changes in the kinetics of paclitaxel 6alpha-hydroxylation, a CYP2C8-specific pathway, and the tolylmethyl and ring hydroxylations of torsemide, a mixed CYP2C9/CYP2C8 substrate. Within SRS1, the single mutation S114F abolished paclitaxel 6alpha-hydroxylation, while the I113V substitution resulted in modest parallel reductions in K(m) and V(max). Mutations in SRS5 (viz., V362L, G365S, and V366L) reduced paclitaxel intrinsic clearance (V(max)/K(m)) by 88-100%. Torsemide is preferentially metabolized by CYP2C9, and it was anticipated that the mutations in CYP2C8 might increase activity. However, methyl and ring hydroxylation intrinsic clearances were either unchanged or decreased by the mutations, although hydroxylation regioselectivity was often altered relative to wild-type CYP2C8. The mutations significantly increased (28-968%) K(m) values for both torsemide methyl and ring hydroxylation but had variable effects on V(max). The effects of the combined mutations in SRS1, SRS5, and SRS1 plus SRS5 were generally consistent with the changes produced by the separate mutations. Mutation of CYP2C8 at position 359 (S359I), a site of genetic polymorphism in CYP2C9, resulted in relatively minor changes in paclitaxel- and torsemide-hydroxylase activities. The results are consistent with multiple substrate binding orientations within the CYP2C8 active site and a differential contribution of active site residues to paclitaxel and torsemide binding and turnover.

Modeling and interactions of Aspergillus fumigatus lanosterol 14-α demethylase `A' with azole antifungals

Recent identification of the sterol 14-α demethylase genes (CYP51 A and B) from Aspergillus fumigatus and other species by Mellado et al. ( J. Clin. Microbiol. 2001, 39(7), 2431–2438), has opened up possibilities of investigating the interactions of azole antifungals with the enzyme(s) from fungi. This study describes for the first time, a model of the three-dimensional structure of A. fumigatus 14-α demethylase (AF-CYP51A), using the crystal structure of Mycobacterium tuberculosis 14-α demethylase (PDB code:1EA1) as a template. The paper also describes the various interactions between azole antifungals and the target from A. fumigatus (AF-CYP51A). Quantitative evaluation of these interactions is done using COMBINE analysis to understand contributions of active site residues to ligand activity. It also provides explanation for the activity/inactivity of different ligands for AF-CYP51A.

Contribution of active site residues to the activity and thermal stability of ribonuclease Sa.

We have used site-specific mutagenesis to study the contribution of Glu 74 and the active site residues Gln 38, Glu 41, Glu 54, Arg 65, and His 85 to the catalytic activity and thermal stability of ribonuclease Sa. The activity of Gln38Ala is lowered by one order of magnitude, which confirms the involvement of this residue in substrate binding. In contrast, Glu41Lys had no effect on the ribonuclease Sa activity. This is surprising, because the hydrogen bond between the guanosine N1 atom and the side chain of Glu 41 is thought to be important for the guanine specificity in related ribonucleases. The activities of Glu54Gln and Arg65Ala are both lowered about 1000-fold, and His85Gln is totally inactive, confirming the importance of these residues to the catalytic function of ribonuclease Sa. In Glu74Lys, k(cat) is reduced sixfold despite the fact that Glu 74 is over 15 A from the active site. The pH dependence of k(cat)/K(M) is very similar for Glu74Lys and wild-type RNase Sa, suggesting that this is not due to a change in the pK values of the groups involved in catalysis. Compared to wild-type RNase Sa, the stabilities of Gln38Ala and Glu74Lys are increased, the stabilities of Glu41Lys, Glu54Gln, and Arg65Ala are decreased and the stability of His85Gln is unchanged. Thus, the active site residues in the ribonuclease Sa make different contributions to the stability.

Contributions of active site residues to the partial and overall catalytic activities of human S-adenosylhomocysteine hydrolase.

Residues glutamate 156 (E156), aspartate 190 (D190), asparagine 181 (N181), lysine 186 (K186), and asparagine 191 (N191) in the active site of S-adenosylhomocysteine (AdoHcy) hydrolase have been mutated to alanine (A). AdoHcy hydrolase achieves catalysis of AdoHcy hydrolysis to adenosine (Ado) and homocysteine (Hcy) by means of a redox partial reaction (3'-oxidation of AdoHcy at the beginning and 3'-reduction of Ado at the end of the catalytic cycle) spanning an elimination/addition partial reaction (elimination of Hcy from the oxidized substrate and addition of water to generate the oxidized product), with the enzyme in an open NAD(+) form in the ligand-free state and in a closed NADH form during the elimination/addition partial reaction. Mutation K186A reduces the rate of a model enzymatic reaction for the redox partial reaction by a factor of 280000 and the rate of a model reaction for the elimination/addition partial reaction by a factor of 24000, consistent with a primary catalytic role in both partial reactions as a proton donor/acceptor at the 3'-OH/3'-keto center. Secondary roles for N181 and N191 in localizing the flexible side chain of K186 in a catalytically effective position are supported by rate reduction factors for N181A of 2500 (redox) and 240 (elimination/addition) and for N191A of 730 (redox) and 340 (elimination/addition). A role of D190 in orienting the substrate for effective transition-state stabilization is consistent with rate reduction factors of 1300 (redox) and 30 (elimination/addition) for D190A. Residue E156 may act to maintain K186 in the desired protonation state: rate deduction factors are 1100 (redox) and 70 (elimination/addition). The mutational increases in free energy barriers for k(cat)/K(M) are described by a linear combination of the effects for the partial reactions with the coefficients equal to the fractional degree that each partial reaction determines the rate for k(cat)/K(M). A similar linear equation for k(cat) overestimates the barrier increase by a uniform 5 kJ/mol, probably reflecting reactant-state stabilization by the wild-type enzyme that is abolished by the mutations.

Lactate monooxygenase. III. Additive contributions of active site residues to catalytic efficiency and stabilization of an anionic transition state.

Lactate monooxygenase catalyzes the conversion of L-lactate to acetate, CO2, and water with incorporation of molecular oxygen. Several amino acid residues of lactate monooxygenase had been postulated to interact in specific ways with the bound substrate (Giegel, D. A., Williams, C. H., Jr., and Massey, V. (1990) J. Biol. Chem. 265, 6626-6632). Tyrosine 44 and arginine 293 were proposed to form a hydrogen bond and a salt bridge to the carboxyl-moiety of lactate. Tyrosine 152 was suggested to form a hydrogen bond to the alpha-hydroxyl group and could be involved in stabilizing a transient carbanionic intermediate of the substrate. The tyrosine residues were replaced with phenylalanines (Y44F, Y152F), and arginine 293 was mutated to a lysine (R293K). In all cases catalysis was significantly decreased; however, the binding affinity for L-lactate did not decrease. Instead, the Kd measured for Y152F was 10-fold lower than that for the wild type enzyme. The products of turnover with Y152F were similar to those with wild type enzyme, with 70-80% of the reaction proceeding to form acetate, CO2, and H2O. The catalytic reactions with both Y44F and R293K were substantially uncoupled, with between 60 and 80% of the catalytic turnover forming pyruvate and H2O2. For all mutant forms the reoxidation of enzyme with oxygen in the absence of pyruvate occurred at a rate similar to that measured for the wild type enzyme. The most important effect of the mutations was in the ability to stabilize the transition state analog oxalate. A linear relationship was found between the rate of reduction of the enzyme flavin and the dissociation constant for the binding of oxalate, demonstrating that many individual residues contribute to the lowering of the energy of the transition state, in addition to specific functions being assignable to some specific residues.