4-oxalocrotonate Tautomerase Research Articles

Cloning and Nucleotide Sequences of Carbazole Degradation Genes from Marine Bacterium Neptuniibacter sp. Strain CAR-SF

The marine bacterium Neptuniibacter sp. strain CAR-SF utilizes carbazole as its sole carbon and nitrogen sources. Two sets of clustered genes related to carbazole degradation, the upper and lower pathways, were obtained. The marine bacterium genes responsible for the upper carbazole degradation pathway, carAa, carBa, carBb, and carC, encode the terminal oxygenase component of carbazole 1,9a-dioxygenase, the small and large subunits of the meta-cleavage enzyme, and the meta-cleavage compound hydrolase, respectively. The genes involved in the lower degradation pathway encode the anthranilate dioxygenase large and small subunit AntA and AntB, anthranilate dioxygenase reductase AntC, 4-oxalocrotonate tautomerase, and catechol 2,3-dioxygenase. Reverse transcription-polymerase chain reaction confirmed the involvement of the isolated genes in carbazole degradation. Escherichia coli cells transformed with the CarAa of strain CAR-SF required ferredoxin and ferredoxin reductase for biotransformation of carbazole. Although carAc, which encodes the ferredoxin component of carbazole 1,9a-dioxygenase, was not found immediately downstream of carAaBaBbC, the carAc-like gene may be located elsewhere based on Southern hybridization. This is the first report of genes involved in carbazole degradation isolated from a marine bacterium.

Uncovering the Protocatechuate 2,3-Cleavage Pathway Genes

Paenibacillus sp. (formerly Bacillus macerans) strain JJ-1b is able to grow on 4-hydroxybenzoate (4HB) as a sole source of carbon and energy and is known to degrade 4HB via the protocatechuate (PCA) 2,3-cleavage pathway. However, none of the genes involved in this pathway have been identified. In this study, we identified and characterized the JJ-1b genes for the 4HB catabolic pathway via the PCA 2,3-cleavage pathway, which consisted of praR and praABEGFDCHI. Based on the enzyme activities of cell extracts of Escherichia coli carrying praI, praA, praH, praB, praC, and praD, these genes were found to code for 4HB 3-hydroxylase, PCA 2,3-dioxygenase, 5-carboxy-2-hydroxymuconate-6-semialdehyde decarboxylase, 2-hydroxymuconate-6-semialdehyde dehydrogenase, 4-oxalocrotonate (OCA) tautomerase, and OCA decarboxylase, respectively, which are involved in the conversion of 4HB into 2-hydroxypenta-2,4-dienoate (HPD). The praE, praF, and praG gene products exhibited 45 to 61% amino acid sequence identity to the corresponding enzymes responsible for the catabolism of HPD to pyruvate and acetyl coenzyme A. The deduced amino acid sequence of praR showed similarity with those of IclR-type transcriptional regulators. Reverse transcription-PCR analysis revealed that praABEGFDCHI constitute an operon, and these genes were expressed during the growth of JJ-1b on 4HB and PCA. praR-praABEGFDCHI conferred the ability to grow on 4HB to E. coli, suggesting that praEGF were functional for the conversion of HPD to pyruvate and acetyl coenzyme A. A promoter analysis suggested that praR encodes a repressor of the pra operon.

Critical role of substrate conformational change in the proton transfer process catalyzed by 4-oxalocrotonate tautomerase.

4-Oxalocrotonate tautomerase enzyme (4-OT) catalyzes the isomerization of 2-oxo-4-hexenedioate to 2-oxo-3-hexenedioate. The chemical process involves two proton transfers, one from a carbon of the substrate to the nitrogen of Pro1 and another from this nitrogen atom to a different carbon of the substrate. In this paper the isomerization has been studied using the combined quantum mechanical and molecular mechanical method with a dual-level treatment of the quantum subsystem employing the MPW1BK density functional as the higher level. Exploration of the potential energy surface shows that the process is stepwise, with a stable intermediate state corresponding to the deprotonated substrate and a protonated proline. The rate constant of the overall process has been evaluated using ensemble-averaged variational transition state theory, including the quantized vibrational motion of a primary zone of active-site atoms and a transmission coefficient based on an ensemble of optimized reaction coordinates to account for recrossing trajectories and optimized multidimensional tunneling. The two proton-transfer steps have similar free energy barriers, but the transition state associated with the first proton transfer is found to be higher in energy. The calculations show that reaction progress is coupled to a conformational change of the substrate, so it is important that the simulation allows this flexibility. The coupled conformational change is promoted by changes in the electron distribution of the substrate that take place as the proton transfers occur.

The chemical versatility of the beta-alpha-beta fold: catalytic promiscuity and divergent evolution in the tautomerase superfamily.

Tautomerase superfamily members have an amino-terminal proline and a beta-alpha-beta fold, and include 4-oxalocrotonate tautomerase (4-OT), 5-(carboxymethyl)-2-hydroxymuconate isomerase (CHMI), trans- and cis-3-chloroacrylic acid dehalogenase (CaaD and cis-CaaD, respectively), malonate semialdehyde decarboxylase (MSAD), and macrophage migration inhibitory factor (MIF), which exhibits a phenylpyruvate tautomerase (PPT) activity. Pro-1 is a base (4-OT, CHMI, the PPT activity of MIF) or an acid (CaaD, cis-CaaD, MSAD). Components of the catalytic machinery have been identified and mechanistic hypotheses formulated. Characterization of new homologues shows that these mechanisms are incomplete. 4-OT, CaaD, cis-CaaD, and MSAD also have promiscuous activities with a hydratase activity in CaaD, cis-CaaD, and MSAD, PPT activity in CaaD and cis-CaaD, and CaaD and cis-CaaD activities in 4-OT. The shared promiscuous activities provide evidence for divergent evolution from a common ancestor, give hints about mechanistic relationships, and implicate catalytic promiscuity in the emergence of new enzymes.

Quantum mechanics/molecular mechanics minimum free-energy path for accurate reaction energetics in solution and enzymes: Sequential sampling and optimization on the potential of mean force surface

To accurately determine the reaction path and its energetics for enzymatic and solution-phase reactions, we present a sequential sampling and optimization approach that greatly enhances the efficiency of the ab initio quantum mechanics/molecular mechanics minimum free-energy path (QM/MM-MFEP) method. In the QM/MM-MFEP method, the thermodynamics of a complex reaction system is described by the potential of mean force (PMF) surface of the quantum mechanical (QM) subsystem with a small number of degrees of freedom, somewhat like describing a reaction process in the gas phase. The main computational cost of the QM/MM-MFEP method comes from the statistical sampling of conformations of the molecular mechanical (MM) subsystem required for the calculation of the QM PMF and its gradient. In our new sequential sampling and optimization approach, we aim to reduce the amount of MM sampling while still retaining the accuracy of the results by first carrying out MM phase-space sampling and then optimizing the QM subsystem in the fixed-size ensemble of MM conformations. The resulting QM optimized structures are then used to obtain more accurate sampling of the MM subsystem. This process of sequential MM sampling and QM optimization is iterated until convergence. The use of a fixed-size, finite MM conformational ensemble enables the precise evaluation of the QM potential of mean force and its gradient within the ensemble, thus circumventing the challenges associated with statistical averaging and significantly speeding up the convergence of the optimization process. To further improve the accuracy of the QM/MM-MFEP method, the reaction path potential method developed by Lu and Yang [Z. Lu and W. Yang, J. Chem. Phys. 121, 89 (2004)] is employed to describe the QM/MM electrostatic interactions in an approximate yet accurate way with a computational cost that is comparable to classical MM simulations. The new method was successfully applied to two example reaction processes, the classical SN2 reaction of Cl-+CH3Cl in solution and the second proton transfer step of the reaction catalyzed by the enzyme 4-oxalocrotonate tautomerase. The activation free energies calculated with this new sequential sampling and optimization approach to the QM/MM-MFEP method agree well with results from other simulation approaches such as the umbrella sampling technique with direct QM/MM dynamics sampling, demonstrating the accuracy of the iterative QM/MM-MFEP method.

Simple Formulas for Improved Point-Charge Electrostatics in Classical Force Fields and Hybrid Quantum Mechanical/Molecular Mechanical Embedding.

We present a simple damping scheme for point-charge electrostatics that could be used directly in classical force fields. The approach acts at the charge (or monopole) level only and allows the inclusion of short-range electrostatic penetration effects at a very low cost. Results are compared with density functional theory Coulomb intermolecular interaction energies and with several other methods such as distributed multipoles, damped distributed multipoles, and transferable Hermite-Gaussian densities. Realistic trends in the interactions are observed for atom-centered Mertz-Kollman corrected point-charge distributions. The approach allows increasing the selectivity of parameters in the case of metal complexes. In addition, two QM/MM calculations are presented where the damping function is employed to include the MM atoms located at the QM/MM boundary. The first calculation corresponds to the gas-phase proton transfer of aspartic acid through water and the second is the first step of the reaction catalyzed by the 4-oxalocrotonate tautomerase (4OT) enzyme. First, improved agreement is observed when using the damping approach compared with the conventional excluded charge method or when including all charges in the calculation. Second, in the case of 4OT, the damped charge approach is in agreement with previous calculations, whereas including all charges gives a significantly higher energy barrier. In both cases, no reparameterization of the van der Waals part of the force field was performed.

Isolation of the phe-operon from G. stearothermophilus comprising the phenol degradative meta-pathway genes and a novel transcriptional regulator

BackgroundGeobacillus stearothermophilus is able to utilize phenol as a sole carbon source. A DNA fragment encoding a phenol hydroxylase catalyzing the first step in the meta-pathway has been isolated previously. Based on these findings a PCR-based DNA walk was performed initially to isolate a catechol 2,3-dioxygenase for biosensoric applications but was continued to elucidate the organisation of the genes encoding the proteins for the metabolization of phenol.ResultsA 20.2 kb DNA fragment was isolated as a result of the DNA walk. Fifteen open reading frames residing on a low-copy megaplasmid were identified. Eleven genes are co-transcribed in one polycistronic mRNA as shown by reverse transcription-PCR. Ten genes encode proteins, that are directly linked with the meta-cleavage pathway. The deduced amino acid sequences display similarities to a two-component phenol hydroxylase, a catechol 2,3-dioxygenase, a 4-oxalocrotonate tautomerase, a 2-oxopent-4-dienoate hydratase, a 4-oxalocrotonate decarboxylase, a 4-hydroxy-2-oxovalerate aldolase, an acetaldehyde dehydrogenase, a plant-type ferredoxin involved in the reactivation of extradiol dioxygenases and a novel regulatory protein. The only enzymes missing for the complete mineralization of phenol are a 2-hydroxymuconic acid-6-semialdehyde hydrolase and/or 2-hydroxymuconic acid-6-semialdehyde dehydrogenase.ConclusionResearch on the bacterial degradation of aromatic compounds on a sub-cellular level has been more intensively studied in gram-negative organisms than in gram-positive bacteria. Especially regulatory mechanisms in gram-positive (thermophilic) prokaryotes remain mostly unknown. We isolated the first complete sequence of an operon from a thermophilic bacterium encoding the meta-pathway genes and analyzed the genetic organization. Moreover, the first transcriptional regulator of the phenol metabolism in gram-positive bacteria was identified. This is a first step to elucidate regulatory mechanisms that are likely to be distinct from modes described for gram-negative bacteria.

Sequential quadratic programming method for determining the minimum energy path

A new method, referred to as the sequential quadratic programming method, is presented for determining minimum energy paths. The method is based on minimizing the points representing the path in the subspace perpendicular to the tangent of the path while using a penalty term to prevent kinks from forming. Rather than taking one full step, the minimization is divided into a number of sequential steps on an approximate quadratic surface. The resulting method can efficiently determine the reaction mechanism, from which transition state can be easily identified and refined with other methods. To improve the resolution of the path close to the transition state, points are clustered close to this region with a reparametrization scheme. The usefulness of the algorithm is demonstrated for the Muller-Brown potential, amide hydrolysis, and an 89 atom cluster taken from the active site of 4-oxalocrotonate tautomerase for the reaction which catalyzes 2-oxo-4-hexenedioate to the intermediate 2-hydroxy-2,4-hexadienedioate.

Kinetic and Stereochemical Analysis of YwhB, a 4-Oxalocrotonate Tautomerase Homologue in Bacillus subtilis: Mechanistic Implications for the YwhB- and 4-Oxalocrotonate Tautomerase-Catalyzed Reactions

YwhB, a 4-oxalocrotonate tautomerase (4-OT) homologue in Bacillus subtilis, has no known biological role, and the gene has no apparent genomic context. The kinetic and stereochemical properties of YwhB have been examined using available enol and dienol compounds. The kinetic analysis shows that YwhB has a relatively nonspecific 1,3- and 1,5-keto-enol tautomerase activity, with the former activity prevailing. Replacement of Pro-1 or Arg-11 with an alanine significantly reduces or abolishes these activities, implicating both residues as critical ones for the activities. In D2O, ketonization of two monoacid substrates (2-hydroxy-2,4-pentadienoate and phenylenolpyruvate) produces a mixture of stereoisomers {2-keto-3-[2H]-4-pentenoate and 3-[2H]-phenylpyruvate}, where the (3R)-isomers predominate. Ketonization of 2-hydroxy-2,4-hexadienedioate, a diacid, in D2O affords mostly the opposite enantiomer, (3S)-2-oxo-[3-2H]-4-hexenedioate. The mono- and diacids apparently bind in different orientations in the active site of YwhB, but the highly stereoselective nature of the YwhB reaction using a diacid suggests that the biological substrate for YwhB may be a diacid. Moreover, of the three dienols examined, 1,3- and 1,5-keto-enol tautomerization reactions are only observed for 2-hydroxy-2,4-hexadienedioate, indicating that the C-3 and C-5 positions are accessible for protonation in this compound. Incubation of 4-OT with 2-hydroxy-2,4-hexadienedioate in D2O results in a racemic mixture of 2-oxo-[3-2H]-4-hexenedioate, suggesting that 4-OT may not catalyze a 1,3-keto-enol tautomerization reaction using this dienol. It has previously been shown that 4-OT catalyzes the near stereospecific conversion of 2-hydroxy-2,4-hexadienedioate to (5S)-[5-2H]-2-oxo-3-hexenedioate in D2O. Taken together, these observations suggest that 4-OT might function as a 1,5-keto-enol tautomerase using 2-hydroxy-2,4-hexadienedioate.

Quantum chemical modeling of enzymatic reactions: The case of 4-oxalocrotonate tautomerase

The reaction mechanism of 4-oxalocrotonate tautomerase (4-OT) is studied using the density functional theory method B3LYP. This enzyme catalyzes the isomerisation of unconjugated α-keto acids to their conjugated isomers. Two different quantum chemical models of the active site are devised and the potential energy curves for the reaction are computed. The calculations support the proposed reaction mechanism in which Pro-1 acts as a base to shuttle a proton from the C3 to the C5 position of the substrate. The first step (proton transfer from C3 to proline) is shown to be the rate-limiting step. The energy of the charge-separated intermediate (protonated proline-deprotonated substrate) is calculated to be quite low, in accordance with measured p K a values. The results of the two models are used to evaluate the methodology employed in modeling enzyme active sites using quantum chemical cluster models.

Substrate Orientation in 4-Oxalocrotonate Tautomerase and Its Effect on QM/MM Energy Profiles

The tautomerization of 2-oxo-4E-hexendioate by 4-oxalocrotonate tautomerase has been studied by quantum mechanical/molecular mechanical (QM/MM) methods using three models, A-C, with different substrate orientations. The computed QM/MM energy profiles are rather different. Various energy partitioning analyses indicate the origin of these differences and the role of the active site residues for different substrate orientations. The proposed new model C is preferred over the previously used models A and B because it combines favorable substrate binding geometries with reasonable barriers and is consistent with the experimental evidence from mutation studies concerning the catalytic ability of specific residues in the binding site, especially R11'.

Composite sequence proteomic analysis of protein biomarkers of Campylobacter coli, C. lari and C. concisus for bacterial identification

Protein biomarkers observed in the matrix-assisted laser desorption/ionization time-of-flight mass spectra (MALDI-TOF-MS) of cell lysates of three strains of Campylobacter coli, two strains of C. lari and one strain of C. concisus have been identified by 'bottom-up' proteomic techniques. The significant findings are as follows. First, the protein biomarkers identified were: PhnA-related protein, 4-oxalocrotonate tautomerase (DmpI)-related protein, NifU-like protein, cytochrome c, DNA-binding protein HU, 10 kDa chaperonin, thioredoxin, as well as several conserved hypothetical and ribosomal proteins. Second, variations in the biomarker ion m/z in MALDI-TOF-MS spectra across species and strains are the result of variations in the amino acid sequence of the protein due to non-synonymous mutations of the biomarker gene. Third, the most common post-translational modifications (PTMs) were the removal of the N-terminal methionine and N-terminal signal peptides. However, in the case of the NifU protein (an iron-sulfur cluster transport protein), post-translational cleavage occurred from the C-terminus. Fourth, only the genomes of the C. coli strain RM2228 and C. lari strain RM2100 have been sequenced; thus, proteomic identification of the proteins of the other strains in this study relied upon sequence homology to the genomic sequence of these strains as well as the genomes of sequences of other Campylobacter strains. In some cases, the determination of the full amino acid sequence of a protein biomarker from a genomically non-sequenced strain was accomplished by combining non-overlapping partial sequences from proteomic identifications of genomically-sequenced strains that were of the same species (or of a different species) to that of the non-sequenced strain. The accuracy of this composite sequence was confirmed by both MS and MS/MS. It was necessary, in some cases, to perform de novo sequencing on 'gaps' in the composite sequence that were not homologous to any genomically-sequenced strain. In order to validate the composite sequence approach, composite sequences were further confirmed by subsequent DNA sequencing of the biomarker gene. Thus, using the composite sequence approach, it was possible to determine the full amino acid sequence of an unknown protein from a genomically non-sequenced bacterial strain without the necessity of either sequencing the biomarker gene or performing full de novo MS/MS sequencing. The sequence obtained could then be used as a strain-specific biomarker for analysis by 'top-down' proteomics techniques.

Metabolism of 4-Amino-3-hydroxybenzoic Acid byBordetellasp. Strain 10d: A Different ModifiedMeta-Cleavage Pathway for 2-Aminophenols

Bordetella sp. strain 10d metabolizes 4-amino-3-hydroxybenzoic acid via 2-hydroxymuconic 6-semialdehyde. Cell extracts from 4-amino-3-hydroxybenzoate-grown cells showed high NAD(+)-dependent 2-hydroxymuconic 6-semialdehyde dehydrogenase, 4-oxalocrotonate tautomerase, 4-oxalocrotonate decarboxylase, and 2-oxopent-4-enoate hydratase activities, but no 2-hydroxymuconic 6-semialdehyde hydrolase activity. These enzymes involved in 4-amino-3-hydroxybenzoate metabolism were purified and characterized. When 2-hydroxymuconic 6-semialdehyde was used as substrate in a reaction mixture containing NAD(+) and cell extracts from 4-amino-3-hydroxybenzoate-grown cells, 4-oxalocrotonic acid, 2-oxopent-4-enoic acid, and 4-hydroxy-2-oxovaleric acid were identified as intermediates, and pyruvic acid was identified as the final product. A complete pathway for the metabolism of 4-amino-3-hydroxybenzoic acid in strain 10d is proposed. Strain 10d metabolized 2-hydroxymuconic 6-semialdehyde derived from 4-amino-3-hydroxybenzoic acid via a dehydrogenative route, not via a hydrolytic route. This proposed metabolic pathway differs considerably from the modified meta-cleavage pathway of 2-aminophenol and those previously reported for methyl- and chloro-derivatives.

Understanding the Enzymatic Activity of 4-Oxalocrotonate Tautomerase and Its Mutant Analogues: A Computational Study

The effect of replacing arginine residues (Arg) with citrulline residues (Cit) in the binding site of 4-oxalocrotonate tautomerase (4-OT) was investigated with force field molecular dynamics and hybrid quantum mechanics/molecular mechanics studies. It is found that the Arg61Cit mutation has only minor effects on the k(cat) and K(M) values determined experimentally because of the flexibility of this residue. The decrease in k(cat) and increase in K(M) for the Arg11Cit and Arg39Cit mutations are due to the disruption of the binding site, which arises from repulsive interactions with neighboring residues. The results of this investigation shed new light on the effects of mutations in the binding site of 4-OT and consequently on how the enzyme binds the active substrate.

Evolution of Enzymatic Activity in the Tautomerase Superfamily: Mechanistic and Structural Consequences of the L8R Mutation in 4-Oxalocrotonate Tautomerase,

4-Oxalocrotonate tautomerase (4-OT) and trans-3-chloroacrylic acid dehalogenase (CaaD) are members of the tautomerase superfamily, a group of structurally homologous proteins that share a beta-alpha-beta fold and a catalytic amino-terminal proline. 4-OT, from Pseudomonas putida mt-2, catalyzes the conversion of 2-oxo-4-hexenedioate to 2-oxo-3-hexenedioate through the dienol intermediate 2-hydroxymuconate in a catabolic pathway for aromatic hydrocarbons. CaaD, from Pseudomonas pavonaceae 170, catalyzes the hydrolytic dehalogenation of trans-3-chloroacrylate in the trans-1,3-dichloropropene degradation pathway. Both reactions may involve an arginine-stabilized enediolate intermediate, a capability that may partially account for the low-level CaaD activity of 4-OT. Two active-site residues in 4-OT, Leu-8 and Ile-52, have now been mutated to the positionally conserved and catalytic ones in CaaD, alphaArg-8, and alphaGlu-52. The L8R and L8R/I52E mutants show improved CaaD activity (50- and 32-fold increases in k(cat)/K(m), respectively) and diminished 4-OT activity (5- and 1700-fold decreases in k(cat)/K(m), respectively). The increased efficiency of L8R-4-OT for the CaaD reaction stems primarily from an 8.8-fold increase in k(cat), whereas that of the L8R/I52E mutant is due largely to a 23-fold decrease in K(m). The presence of the additional arginine residue in the active site of L8R-4-OT does not alter the pK(a) of the Pro-1 amino group from that measured for the wild type (6.5 +/- 0.1 versus 6.4 +/- 0.2). Moreover, the crystal structure of L8R-4-OT is comparable to that of the wild type. Hence, the enhanced CaaD activity of L8R-4-OT is likely due to the additional arginine residue that can participate in substrate binding and/or stabilization of the putative enediolate intermediate. The results also suggest that the evolution of new functions within the tautomerase superfamily could be quite facile, requiring only a few strategically placed active-site mutations.

A persistent pesticide residue and the unusual catalytic proficiency of a dehalogenating enzyme

The soil of potato fields in The Netherlands harbors bacteria with the ability to metabolize 3-chloroacrylic acid, generated by the degradation of a pesticide (1,3-dichloropropene) that entered the environment in 1946. From examination of rate constants at elevated temperatures, we infer that the half-time at 25 degrees C for spontaneous hydrolytic dechlorination of trans-3-chloroacrylic acid is 10,000 years, several orders of magnitude longer than half-times for spontaneous decomposition of other environmental pollutants such as 1,2-dichloroethane (72 years), paraoxon (13 months), atrazine (5 months), and aziridine (52 h). With thermodynamic parameters for activation similar to those for the spontaneous hydration of fumarate at pH 6.8, this slow reaction proceeds at a constant rate through the pH range between 2 and 12. However, at the active site of the enzyme 3-chloroacrylate dehalogenase (CaaD), isolated from a pseudomonad growing in these soils, hydrolytic dechlorination proceeds with a half-time of 0.18 s. Neither k(cat) nor k(cat)/K(m) is reduced by increasing solvent viscosity with trehalose, implying that the rate of enzymatic dechlorination is controlled by chemical events in catalysis rather than by diffusion-limited substrate binding or product release. CaaD achieves an approximately 10(12)-fold rate enhancement, matching or surpassing the rate enhancements produced by many enzymes that act on more conventional biological substrates. One of those enzymes is 4-oxalocrotonate tautomerase, with which CaaD seems to share a common evolutionary origin.

Theoretical and Experimental Determination on Two Substrates Turned over by 4-Oxalocrotonate Tautomerase

Quantum mechanical/molecular mechanical (QM/MM) calculations and experimental kinetic studies have been performed on 4-oxalocrotonate tautomerase (4OT) for two different substrates, 2-hydroxymuconate (2HM) and 2-oxo-4-hexenedioate (2o4hex). Potential (deltaE) and free energy (deltaG) paths for both steps of the reaction using both substrates were calculated to determine the free energy barriers and compared to the experimental values obtained from the kinetic studies via the transition state theory. In the first step, a proton from the hydroxyl oxygen on the second carbon of 2HM, or from the third carbon of 2o4hex, is abstracted by Pro-1. In the second step, the proton is transferred to the fifth carbon of the substrate to form the product, 2-oxo-3-hexenedioate (2o3hex). For both substrates we obtain a calculated deltaG of approximately 13 kcal/mol, in agreement with experimental determinations. The calculated free energy barrier difference deltaG2o4hex - deltaG2HM (deltadeltaG) is 0.87 kcal/mol. We obtained an experimental deltadeltaG of 0.85 kcal/mol. These results suggest that 2HM is turned over faster than 2o4hex by 4OT. However, these energy differences are so small that both 2HM and 2o4hex need to be taken into account in considering the mechanism of catalysis of 4OT.

Half-of-the-Sites Binding of Reactive Intermediates and Their Analogues to 4-Oxalocrotonate Tautomerase and Induced Structural Asymmetry of the Enzyme

4-Oxalocrotonate tautomerase (4-OT), a homohexameric enzyme, converts the unconjugated enone, 2-oxo-4-hexenedioate (1), to the conjugated enone, 2-oxo-3-hexenedioate (3), via a dienolic intermediate, 2-hydroxymuconate (2). Pro-1 serves as the general base, and both Arg-11 and Arg-39 function in substrate binding and catalysis in an otherwise hydrophobic active site. Although 4-OT exhibits hyperbolic kinetics and no structural asymmetry either by X-ray or by NMR, inactivation by two affinity labels showed half-site stoichiometry [Stivers, J. T., et al. (1996) Biochemistry 35, 803-813; Johnson, W. H., Jr., et al. (1997) Biochemistry 36, 15724-15732], and titration of the R39Q mutant with cis,cis-muconate showed negative cooperativity [Harris, T. K., et al. (1999) Biochemistry 38, 12343-12357]. To test for anticooperativity during catalysis, 4-OT was titrated with equilibrium mixtures (> or = 81% product) of the reactive dicarboxylate or monocarboxylate intermediates, 2 or 2-hydroxy-2,4-pentadienoate (4), respectively, in three types of NMR experiments: two-dimensional 1H-15N HSQC titrations of backbone NH and of Arg N epsilonH resonances and one-dimensional 15N NMR titrations of Arg N epsilon resonances. All titrations showed substoichiometric binding of the equilibrium mixtures to 3 +/- 1 sites per hexamer with apparent dissociation constants comparable to the Km values of the intermediates. Compound 4 also bound 1 order of magnitude less tightly at another site, suggesting negative cooperativity. Consistent with negative cooperativity, asymmetry of the resulting complexes at saturating levels of 2 and 4 is indicated by splitting of the backbone NH resonances of 11 residues and 10 residues of 4-OT, respectively. The dicarboxylate competitive inhibitor, (2E)-fluoromuconate (5), with a KI of 45 +/- 7 microM, also exhibited substoichiometric binding to 3 +/- 1 sites per hexamer, with a KD of 25 +/- 18 microM, and splitting of the backbone NH resonance of L8. The monocarboxylate inhibitors (2E)- (6) and (2Z)-2-fluoro-2,4-pentadienoate (7) showed much weaker binding (KD = 3.1 +/- 1.3 mM), as well as splitting of two and five backbone NH resonances, respectively, indicating asymmetry of the complexes. The N epsilon resonances of both Arg-11 and Arg-39 were shifted downfield, and that of Pro-1N was broadened by all ligands, consistent with the major catalytic roles of these residues. Structural pathways for the site-site interactions which result in negative cooperativity are proposed on the basis of the X-ray structures of free and affinity-labeled 4-OT. Selective resonance broadenings induced by the binding of inactive analogues and active intermediates indicate residues which may be mobilized during reversible ligand binding and during catalysis, respectively.

A Designed Synthetic Analogue of 4-OT Is Specific for a Non-Natural Substrate

The substrate specificity of 4-oxalocrotonate tautomerase (4-OT) is characterized by electrostatic interactions between positively charged arginine (Arg) side chains on the enzyme and the dianionic substrate, 4-oxalocrotonate. To generate specific hydrogen-bonding interactions with a monoanionic substrate analogue, we have introduced a urea functional group into the active site by replacing arginine side chains with isosteric citrulline (Cit) residues. This design was based on the complementarity between the urea functionality of citrulline and the uncharged amide function of the substrate, as opposed to the guanidinium-carboxylate electrostatic interaction between the wild-type enzyme and the natural substrate. Indeed, the synthetic (Arg39Cit)4-OT analogue catalyzed the tautomerization of the non-natural monoamide-monoacid substrate while it was a poor catalyst for the natural diacid substrate. The specificity of (Arg39Cit)4-OT for the monoamide-monoacid substrate relative to that of the diacid substrate was found to be 740-fold greater than that of the wild-type enzyme for tautomerization of the non-natural substrate as compared with the natural one. The role of electrostatic interactions in the tautomerization of the monoamide-monoacid substrate was probed in detail with several other Arg to Cit analogues of this enzyme. This study has demonstrated that chemical manipulation of the functional groups within the active site of an enzyme can modify its catalytic activity and substrate specificity in a predictable way, suggesting that the incorporation of noncoded amino acids into proteins has great promise for the development of new enzymatic mechanisms and new binding interactions.

Reaction path determination for quantum mechanical/molecular mechanical modeling of enzyme reactions by combining first order and second order “chain-of-replicas” methods

A two-step procedure for the determination of reaction paths in enzyme systems is presented. This procedure combines two chain-of-states methods: a quantum mechanical/molecular mechanical (QM/MM) implementation of the nudged elastic band (NEB) method and a second order parallel path optimizer method both recently developed in our laboratory. In the first step, a reaction path determination is performed with the NEB method, along with a restrained minimization procedure for the MM environment to obtain a first approximation to the reaction path. In the second step, the calculated path is refined with the parallel path optimizer method. By combining these two methods the reaction paths are determined accurately, and in addition, the number of path optimization iterations are significantly reduced. This procedure is tested by calculating both steps of the isomerization of 2-oxo-4-hexenedioate by 4-oxalocrotonate tautomerase, which have been previously determined by our group. The calculated paths agree with the previously reported results and we obtain a reduction of 45%-55% in the number of path optimization cycles.