Activity of yeast orotidine-5'-phosphate decarboxylase in the absence of metals

Limits to Compensatory Mutations: Insights from Temperature-Sensitive Alleles

Previous experiments with temperature-sensitive mutants of the yeast enzyme orotidine 5'-phosphate decarboxylase (encoded in gene URA3) yielded the unexpected result that reversion occurs only through exact reversal of the original mutation (Jakubowska A, Korona R. 2009. Lack of evolutionary conservation at positions important for thermal stability in the yeast ODCase protein. Mol Biol Evol. 26(7):1431-1434.). We recreated a set of these mutations in which the codon had two nucleotide substitutions, making exact reversion much less likely. We screened these double mutants for reversion and obtained a number of compensatory mutations occurring at alternative sites in the molecule. None of these compensatory mutations fully restored protein performance. The mechanism of partial compensation is consistent with a model in which protein stabilization is additive, as the same secondary mutations can compensate different primary alternations. The distance between primary and compensatory residues precludes direct interaction between the sites. Instead, most of the compensatory mutants were clustered in proximity to the catalytic center. All of the second-site compensatory substitutions occurred at relatively conserved sites, and the amino acid replacements were to residues found at these sites in a multispecies alignment of the protein. Based on the estimated distribution of changes in Gibbs free energy among a large number of amino acid replacements, we estimate that, for most proteins, the probability that a second-site mutation would have a sufficiently large stabilizing effect to offset a temperature-sensitive mutation in the order of 10-4 or less. Hence compensation is likely to take place only for slightly destabilizing mutations because highly stabilizing mutations are exceeding rare.

Enzyme architecture: deconstruction of the enzyme-activating phosphodianion interactions of orotidine 5'-monophosphate decarboxylase

The mechanism for activation of orotidine 5'-monophosphate decarboxylase (OMPDC) by interactions of side chains from Gln215 and Try217 at a gripper loop and R235, adjacent to this loop, with the phosphodianion of OMP was probed by determining the kinetic parameters k(cat) and K(m) for all combinations of single, double, and triple Q215A, Y217F, and R235A mutations. The 12 kcal/mol intrinsic binding energy of the phosphodianion is shown to be equal to the sum of the binding energies of the side chains of R235 (6 kcal/mol), Q215 (2 kcal/mol), Y217 (2 kcal/mol), and hydrogen bonds to the G234 and R235 backbone amides (2 kcal/mol). Analysis of a triple mutant cube shows small (ca. 1 kcal/mol) interactions between phosphodianion gripper side chains, which are consistent with steric crowding of the side chains around the phosphodianion at wild-type OMPDC. These mutations result in the same change in the activation barrier to the OMPDC-catalyzed reactions of the whole substrate OMP and the substrate pieces (1-β-D-erythrofuranosyl)orotic acid (EO) and phosphite dianion. This shows that the transition states for these reactions are stabilized by similar interactions with the protein catalyst. The 12 kcal/mol intrinsic phosphodianion binding energy of OMP is divided between the 8 kcal/mol of binding energy, which is utilized to drive a thermodynamically unfavorable conformational change of the free enzyme, resulting in an increase in (k(cat))(obs) for OMPDC-catalyzed decarboxylation of OMP, and the 4 kcal/mol of binding energy, which is utilized to stabilize the Michaelis complex, resulting in a decrease in (K(m))(obs).

Proton transfer from C-6 of uridine 5'-monophosphate catalyzed by orotidine 5'-monophosphate decarboxylase: formation and stability of a vinyl carbanion intermediate and the effect of a 5-fluoro substituent

The exchange for deuterium of the C-6 protons of uridine 5'-monophosphate (UMP) and 5-fluorouridine 5'-monophosphate (F-UMP) catalyzed by yeast orotidine 5'-monophosphate decarboxylase (ScOMPDC) at pD 6.5-9.3 and 25 °C was monitored by (1)H NMR spectroscopy. Deuterium exchange proceeds by proton transfer from C-6 of the bound nucleotide to the deprotonated side chain of Lys-93 to give the enzyme-bound vinyl carbanion. The pD-rate profiles for k(cat) give turnover numbers for deuterium exchange into enzyme-bound UMP and F-UMP of 1.2 × 10(-5) and 0.041 s(-1), respectively, so that the 5-fluoro substituent results in a 3400-fold increase in the first-order rate constant for deuterium exchange. The binding of UMP and F-UMP to ScOMPDC results in 0.5 and 1.4 unit decreases, respectively, in the pK(a) of the side chain of the catalytic base Lys-93, showing that these nucleotides bind preferentially to the deprotonated enzyme. We also report the first carbon acid pK(a) values for proton transfer from C-6 of uridine (pK(CH) = 28.8) and 5-fluorouridine (pK(CH) = 25.1) in aqueous solution. The stabilizing effects of the 5-fluoro substituent on C-6 carbanion formation in solution (5 kcal/mol) and at ScOMPDC (6 kcal/mol) are similar. The binding of UMP and F-UMP to ScOMPDC results in a greater than 5 × 10(9)-fold increase in the equilibrium constant for proton transfer from C-6, so that ScOMPDC stabilizes the bound vinyl carbanions, relative to the bound nucleotides, by at least 13 kcal/mol. The pD-rate profile for k(cat)/K(m) for deuterium exchange into F-UMP gives the intrinsic second-order rate constant for exchange catalyzed by the deprotonated enzyme as 2300 M(-1) s(-1). This was used to calculate a total rate acceleration for ScOMPDC-catalyzed deuterium exchange of 3 × 10(10) M(-1), which corresponds to a transition-state stabilization for deuterium exchange of 14 kcal/mol. We conclude that a large portion of the total transition-state stabilization for the decarboxylation of orotidine 5'-monophosphate can be accounted for by stabilization of the enzyme-bound vinyl carbanion intermediate of the stepwise reaction.

Mechanism of OMP decarboxylation in orotidine 5'-monophosphate decarboxylase

Despite extensive experimental and theoretical studies, the detailed catalytic mechanism of orotidine 5'-monophosphate decarboxylase (ODCase) remains controversial. In particular simulation studies using high level quantum mechanics have failed to reproduce experimental activation free energy. One common feature of many previous simulations is that there is a water molecule in the vicinity of the leaving CO2 group whose presence was only observed in the inhibitor bound complex of ODCase/BMP. Various roles have even been proposed for this water molecule from the perspective of stabilizing the transition state and/or intermediate state. We hypothesize that this water molecule is not present in the active ODCase/OMP complex. Based on QM/MM minimum free energy path simulations with accurate density functional methods, we show here that in the absence of this water molecule the enzyme functions through a simple direct decarboxylation mechanism. Analysis of the interactions in the active site indicates multiple factors contributing to the catalysis, including the fine-tuned electrostatic environment of the active site and multiple hydrogen-bonding interactions. To understand better the interactions between the enzyme and the inhibitor BMP molecule, simulations were also carried out to determine the binding free energy of this special water molecule in the ODCase/BMP complex. The results indicate that the water molecule in the active site plays a significant role in the binding of BMP by contributing approximately -3 kcal/mol to the binding free energy of the complex. Therefore, the complex of BMP plus a water molecule, instead of the BMP molecule alone, better represents the tight binding transition state analogue of ODCase. Our simulation results support the direct decarboxylation mechanism and highlight the importance of proper recognition of protein bound water molecules in the protein-ligand binding and the enzyme catalysis.

The effective molarity of the substrate phosphoryl group in the transition state for yeast OMP decarboxylase

The second order rate constant (k(cat)/K(m)) for decarboxylation of orotidine by yeast OMP decarboxylase (ODCase), measured by trapping (14)CO(2) released during the reaction, is 2 x 10(-4)M(-1)s(-1). This very low activity may be compared with a value of 3 x 10(7)M(-1)s(-1) for the action of yeast OMP decarboxylase on the normal substrate OMP. Both activities are strongly inhibited by 6-hydroxy UMP (BMP), and abrogated by mutation of Asp-96 to alanine. These results, in conjunction with the binding affinity of inorganic phosphate as a competitive inhibitor (K(i)=7 x 10(-4)M), imply an effective concentration of 1.1 x 10(9)M for the substrate phosphoryl group in stabilizing the transition state for enzymatic decarboxylation of OMP. The observed difference in rate (1.5 x 10(11)-fold) is the largest effect of a simple substituent that appears to have been reported for an enzyme reaction.

On the importance of being zwitterionic: enzymatic catalysis of decarboxylation and deprotonation of cationic carbon

Carbanion ylides are strongly stabilized by electrostatic interactions between opposing charges at neighboring atoms and this stabilizing electrostatic interaction increases with decreasing dielectric constant of the medium through which the charges interact. Consequently, there is a large increase in the thermodynamic driving force, with decreasing dielectric constant of the reaction medium, for deprotonation of cationic carbon acids and decarboxylation to form related ylides. This favors catalysis of the formation of unstable ylides at enzyme active sites of low dielectric constant. A brief survey of enzymes that catalyze deprotonation of cationic carbon acids and related decarboxylation reactions shows catalysis generally occurs for substrates that are bound in a deep pocket on the protein, with an apparent dielectric constant that is much lower than for the solvent water. In several cases, proton transfer is to a catalytic residue that is relatively weakly solvated in water. We suggest that there is a strong advantage for evolution of protein catalysts that utilize weakly solvated basic side chains which are relatively easily buried in nonpolar active sites that are favorable for zwitterion formation.

Catalytic proficiency: the unusual case of OMP decarboxylase

Enzymes are called upon to differ greatly in the difficulty of the tasks that they perform. The catalytic proficiency of an enzyme can be evaluated by comparing the second-order rate constant (kcat/Km) with the rate of the spontaneous reaction in neutral solution in the absence of a catalyst. The proficiencies of enzymes, measured in this way, are matched by their affinity constants for the altered substrate in the transition state. These values vary from approximately approximately 10(9) M(-1) for carbonic anhydrase to approximately 10(23) M(-1) for yeast orotidine 5'-phosphate decarboxylase (ODCase). ODCase turns its substrate over with a half-time of 18 ms, in a reaction that proceeds in its absence with a half-time of 78 million years in neutral solution. ODCase differs from other decarboxylases in that its catalytic activity does not depend on the presence of metals or other cofactors, or on the formation of a covalent bond to the substrate. Several mechanisms of transition state stabilization are considered in terms of ODCase crystal structures observed in the presence and absence of bound analogs of the substrate, transition state, and product. Very large connectivity effects are indicated by the results of experiments testing how transition state stabilization is affected by the truncation of binding determinants of the substrate and the active site.

Crystal structures of inhibitor complexes reveal an alternate binding mode in orotidine-5'-monophosphate decarboxylase

The crystal structures of the enzyme orotidine-5'-monophosphate decarboxylase from Methanobacterium thermoautotrophicum complexed with its product UMP and the inhibitors 6-hydroxyuridine 5'-phosphate (BMP), XMP, and CMP are reported. A mutant version of the protein, in which four residues of the flexible phosphate-binding loop (180)Gly-Gly(190) were removed and Arg(203) was replaced by alanine, was also analyzed. The XMP and CMP complexes reveal a ligand-binding mode that is distinct from the one identified previously with the aromatic rings located outside the binding pocket. A potential pathway for ligand binding is discussed.

Divergent evolution of enzymatic function: mechanistically diverse superfamilies and functionally distinct suprafamilies

The protein sequence and structure databases are now sufficiently representative that strategies nature uses to evolve new catalytic functions can be identified. Groups of divergently related enzymes whose members catalyze different reactions but share a common partial reaction, intermediate, or transition state (mechanistically diverse superfamilies) have been discovered, including the enolase, amidohydrolase, thiyl radical, crotonase, vicinal-oxygen-chelate, and Fe-dependent oxidase superfamilies. Other groups of divergently related enzymes whose members catalyze different overall reactions that do not share a common mechanistic strategy (functionally distinct suprafamilies) have also been identified: (a) functionally distinct suprafamilies whose members catalyze successive transformations in the tryptophan and histidine biosynthetic pathways and (b) functionally distinct suprafamilies whose members catalyze different reactions in different metabolic pathways. An understanding of the structural bases for the catalytic diversity observed in super- and suprafamilies may provide the basis for discovering the functions of proteins and enzymes in new genomes as well as provide guidance for in vitro evolution/engineering of new enzymes.

Dissecting a charged network at the active site of orotidine-5'-phosphate decarboxylase

The crystal structure of yeast orotidine-5'-phosphate decarboxylase in complex with the postulated transition state analog, 6-hydroxyuridine-5'-phosphate, reveals contacts between this inhibitor and a novel quartet of charged residues (Lys-59, Asp-91, Lys-93, and Asp-96) within the active site. The structure also suggests a possible interaction between O2 of the 6-hydroxyuridine-5'-phosphate pyrimidine ring and Gln-215. Here we report the results of mutagenesis of each of the charged active site residues and Gln-215. The activities of the Q215A and wild-type enzymes were equal indicating that any interactions between this residue and the pyrimidine ring are dispensable for efficient decarboxylation. For the D91A and K93A mutant enzymes, activity was reduced by more than 5 orders of magnitude and substrate binding could not be detected by isothermal calorimetry. For the D96A mutant enzyme, k(cat) was reduced by more than 5 orders of magnitude, and isothermal calorimetry indicated an 11-fold decrease in the affinity of this enzyme for the substrate in the ground state. For the K59A enzyme, k(cat) was reduced by a factor of 130, and K(m) had increased by a factor of 900. These results indicate that the integrity of the network of charged residues is essential for transition state stabilization.

Yeast orotidine-5'-phosphate decarboxylase: steady-state and pre-steady-state analysis of the kinetic mechanism of substrate decarboxylation

The catalytically active form of monofunctional yeast orotidine-5'-phosphate decarboxylase was a dimer (E(2)). The dimer equilibrium dissociation constant was 0.25 microM in 0.01 M MOPS Na(+) at pH 7.2. The bimolecular rate constant for dimer formation was 1.56 microM(-1) s(-1). The dimeric form of the enzyme was stabilized by NaCl such that the enzyme was E(2) in 100 mM NaCl at all concentrations of enzyme tested. The kinetics of binding of OMP to E(2) was governed by two ionizations (pK(1) = 6.1 and pK(2) = 7.7). From studies with substrate analogues, the higher pK was assigned to a group on the enzyme that interacted with the pyrimidinyl moiety. The value of the lower pK was dependent on the substrate analogue, which suggested that it was not exclusively the result of ionization of the phosphoryl moiety. During the decarboxylation of OMP, the fluorescence of E(2) was quenched over 20%. The enzymatic species with reduced fluorescence was a catalytically competent intermediate that had kinetic properties consistent with it being the initial enzyme-substrate complex. The stoichiometry for binding of OMP to E(2) was one OMP per enzyme monomer. The value of the first-order rate constant for conversion of the enzyme-substrate complex to free enzyme (36 s(-1)) calculated from a single turnover experiment ([E] >> [S]) was slightly greater than the value of k(cat), 20 s(-1) (corrected for stoichiometry), calculated from steady-state data. In the single turnover experiments, the enzyme was E(2)*S, whereas in the steady-state turnover the experiment enzyme was E(2)*S(2). The similarity of these values suggested that the subunits were catalytically independent such that E(2)*S(2) could be treated as E*S and that conversion of the enzyme-substrate complex to E was k(cat). Kinetic data for the approach to the steady-state with OMP and E(2) yield a bimolecular association rate complex of 62 microM(-1) s(-1)and a dissociation rate constant for E*S of 60 s(-1). The commitment to catalysis was 0.25. By monitoring the effect of carbonic anhydrase on [H(+)] changes during a single turnover experiment, the initial product of the decarboxylation reaction was shown to be CO(2) not HCO(3-). UMP was released from the enzyme concomitantly with CO(2) during the conversion of E*S to E. Furthermore, the enzyme removed an enzyme equivalent of H(+) from solvent during this step of the reaction. The bimolecular rate constants for association of 6-AzaUMP and 8-AzaXMP, substrate analogues with markedly different nucleobases, had association rate constants of 112 and 130 microM(-1) s(-1), respectively. These results suggested that the nucleobase did not contribute significantly to the success of formation of the initial enzyme-substrate complex.

The chemistry of the reaction determines the invariant amino acids during the evolution and divergence of orotidine 5'-monophosphate decarboxylase

Orotidine 5'-phosphate (OMP) decarboxylase has the largest rate enhancement for any known enzyme. For an average protein of 270 amino acids from more than 80 species, only 8 amino acids are invariant, and 7 of these correspond to ligand-binding residues in the crystal structures of the enzyme from four species. It appears that the chemistry required for catalysis determines the invariant residues for this enzyme structure. A motif of three invariant amino acids at the catalytic site (DXKXXD) is also found in the enzyme hexulose-phosphate synthase. Although the core of OMP decarboxylase is conserved, it has undergone a variety of changes in subunit size or fusion to other protein domains, such as orotate phosphoribosyltransferase, during evolution in different kingdoms. The phylogeny of OMP decarboxylase shows a unique subgroup distinct from the three kingdoms of life. The enzyme subunit size almost doubles from Archaea (average mass of 24.5 kDa) to certain fungi (average mass of 41.7 kDa). These observed changes in subunit size are produced by insertions at 12 sites, largely in loops and on the exterior of the core protein. The consensus for all sequences has a minimal size of <20 kDa.

Electrostatic stress in catalysis: structure and mechanism of the enzyme orotidine monophosphate decarboxylase

Orotidine 5'-monophosphate decarboxylase catalyzes the conversion of orotidine 5'-monophosphate to uridine 5'-monophosphate, the last step in biosynthesis of pyrimidine nucleotides. As part of a Structural Genomics Initiative, the crystal structures of the ligand-free and the6-azauridine 5'-monophosphate-complexed forms have been determined at 1.8 and 1.5 A, respectively. The protein assumes a TIM-barrel fold with one side of the barrel closed off and the other side binding the inhibitor. A unique array of alternating charges (Lys-Asp-Lys-Asp) in the active site prompted us to apply quantum mechanical and molecular dynamics calculations to analyze the relative contributions of ground state destabilization and transition state stabilization to catalysis. The remarkable catalytic power of orotidine 5'-monophosphate decarboxylase is almost exclusively achieved via destabilization of the reactive part of the substrate, which is compensated for by strong binding of the phosphate and ribose groups. The computational results are consistent with a catalytic mechanism that is characterized by Jencks's Circe effect.

Anatomy of a proficient enzyme: the structure of orotidine 5'-monophosphate decarboxylase in the presence and absence of a potential transition state analog

Orotidine 5'-phosphate decarboxylase produces the largest rate enhancement that has been reported for any enzyme. The crystal structure of the recombinant Saccharomyces cerevisiae enzyme has been determined in the absence and presence of the proposed transition state analog 6-hydroxyuridine 5'-phosphate, at a resolution of 2.1 A and 2.4 A, respectively. Orotidine 5'-phosphate decarboxylase folds as a TIM-barrel with the ligand binding site near the open end of the barrel. The binding of 6-hydroxyuridine 5'-phosphate is accompanied by protein loop movements that envelop the ligand almost completely, forming numerous favorable interactions with the phosphoryl group, the ribofuranosyl group, and the pyrimidine ring. Lysine-93 appears to be anchored in such a way as to optimize electrostatic interactions with developing negative charge at C-6 of the pyrimidine ring, and to donate the proton that replaces the carboxylate group at C-6 of the product. In addition, H-bonds from the active site to O-2 and O-4 help to delocalize negative charge in the transition state. Interactions between the enzyme and the phosphoribosyl group anchor the pyrimidine within the active site, helping to explain the phosphoribosyl group's remarkably large contribution to catalysis despite its distance from the site of decarboxylation.

Can sequence determine function?

The functional annotation of proteins identified in genome sequencing projects is based on similarities to homologs in the databases. As a result of the possible strategies for divergent evolution, homologous enzymes frequently do not catalyze the same reaction, and we conclude that assignment of function from sequence information alone should be viewed with some skepticism.