Model chemistry for a covalent mechanism of action of orotidine 5'-phosphate decarboxylase

The Enzymatic Decarboxylation Mechanism of 5-Carboxy Uracil: A Comprehensive Quantum Chemical Study

Dynamic regulation of DNA methylation is an important process for the control of gene expression in mammals. It is believed that in the demethylation pathway of 5-methyl cytosine, the intermediate 5-carboxy cytosine (5caC) can be actively decarboxylated alongside the substitution in the base excision repair. For the active decarboxylation of 5caC, a decarboxylase has not been identified so far. Due to the similar chemistry of the decarboxylation of 5-carboxy uracil (5caU) to uracil (U) in the pyrimidine salvage pathway catalyzed by the iso-orotate decarboxylase (IDCase), the study of this reaction might give valuable insights into the active 5caC decarboxylation process. In this work, we employ quantum chemical and molecular mechanic calculations and find that the catalytic mechanism of IDCase proceeds via a direct decarboxylation mechanism. Detailed investigations on the reaction coordinate reveal that it is a one-step mechanism with concerted proton transfer and C-C bond opening.

CRISPR/Cas9 genome editing in ergot fungus Claviceps purpurea

Claviceps purpurea is a filamentous fungus well known as a widespread plant pathogen, but it is also an important ergot alkaloid producer exploited by the pharmaceutic industry. In this work, we demonstrated that CRISPR/Cas9 can be a tool for directed mutagenesis in C. purpurea targeting pyr4 and TrpE genes encoding the orotidine 5'-phosphate decarboxylase involved in pyrimidine biosynthesis and the α-subunit of the anthranilate synthase involved in tryptophan biosynthesis, respectively. After protoplast transformation and single spore isolation, homokaryotic mutants showing uridine or tryptophan auxotrophy were selected. In all cases, insertions or insertions combined with deletions were found mostly 3 bp upstream of the PAM sequence. However, transformation efficiencies of CRISPR/Cas9 and CRISPR/Cas9 mediated homology-directed repair only slightly improved in comparison to homologous recombination-mediated knocking-out of the TrpE gene. Moreover, Trp auxotrophs were non-infectious towards rye plants likely due to a decreased production of the plant hormones auxins, which are synthesized by C. purpurea from indole-3-glycerolphosphate in Trp-dependent and Trp-independent biosynthetic pathways, and help the fungus to colonize the plant host. It was demonstrated that the CRISPR/Cas9 vector containing autonomous replicative sequence AMA1 can be fully removed by further culturing of C. purpurea on non-selective media. This method enables introducing multiple mutations in Claviceps and makes feasible metabolic engineering of industrial strains.

Orotidine Monophosphate Decarboxylase--A Fascinating Workhorse Enzyme with Therapeutic Potential

Orotidine 5'-monophosphate decarboxylase (ODCase) is known as one of the most proficient enzymes. The enzyme catalyzes the last reaction step of the de novo pyrimidine biosynthesis, the conversion from orotidine 5'-monophosphate (OMP) to uridine 5'-monophosphate. The enzyme is found in all three domains of life, Bacteria, Eukarya and Archaea. Multiple sequence alignment of 750 putative ODCase sequences resulted in five distinct groups. While the universally conserved DxKxxDx motif is present in all the groups, depending on the groups, several characteristic motifs and residues can be identified. Over 200 crystal structures of ODCases have been determined so far. The structures, together with biochemical assays and computational studies, elucidated that ODCase utilized both transition state stabilization and substrate distortion to accelerate the decarboxylation of its natural substrate. Stabilization of the vinyl anion intermediate by a conserved lysine residue at the catalytic site is considered the largest contributing factor to catalysis, while bending of the carboxyl group from the plane of the aromatic pyrimidine ring of OMP accounts for substrate distortion. A number of crystal structures of ODCases complexed with potential drug candidate molecules have also been determined, including with 6-iodo-uridine, a potential antimalarial agent.

Study of orotidine 5'-monophosphate decarboxylase in complex with the top three OMP, BMP, and PMP ligands by molecular dynamics simulation

Catalytic mechanism of orotidine 5'-monophosphate decarboxylase (OMPDC), one of the nature most proficient enzymes which provides large rate enhancement, has not been fully understood yet. A series of 30 ns molecular dynamics (MD) simulations were run on X-ray structure of the OMPDC from Saccharomyces cerevisiae in its free form as well as in complex with different ligands, namely 1-(5'-phospho-D-ribofuranosyl) barbituric acid (BMP), orotidine 5'-monophosphate (OMP), and 6-phosphonouridine 5'-monophosphate (PMP). The importance of this biological system is justified both by its high rate enhancement and its potential use as a target in chemotherapy. This work focuses on comparing two physicochemical states of the enzyme (protonated and deprotonated Asp91) and three ligands (substrate OMP, inhibitor, and transition state analog BMP and substrate analog PMP). Detailed analysis of the active site geometry and its interactions is properly put in context by extensive comparison with relevant experimental works. Our overall results show that in terms of hydrogen bond occupancy, electrostatic interactions, dihedral angles, active site configuration, and movement of loops, notable differences among different complexes are observed. Comparison of the results obtained from these simulations provides some detailed structural data for the complexes, the enzyme, and the ligands, as well as useful insights into the inhibition mechanism of the OMPDC enzyme. Furthermore, these simulations are applied to clarify the ambiguous mechanism of the OMPDC enzyme, and imply that the substrate destabilization and transition state stabilization contribute to the mechanism of action of the most proficient enzyme, OMPDC.

Design of inhibitors of ODCase

ODCase is a highly proficient enzyme responsible for the decarboxylation of orotidine monophosphate to generate uridine monophosphate. ODCase has attracted early attention due to its interesting mechanism of catalysis. In order to exploit therapeutic advantages due to the inhibition of ODCase, one must have selective inhibitors of this enzyme from the pathogen, or a dysregulated molecular mechanism involving ODCase. ODCase inhibitors have potential applications as anticancer agents, antiviral agents, antimalarial agents and potentially act against other parasitic diseases. A variety of C6-substituted uridine monophosphate derivatives have shown excellent inhibition of ODCase. 6-iodouridine is a potent inhibitor of the malaria parasite, and its monophosphate form covalently inhibits ODCase. A variety of inhibitors of ODCase with potential applications as therapeutic agents are discussed in this review.

Role of a guanidinium cation-phosphodianion pair in stabilizing the vinyl carbanion intermediate of orotidine 5'-phosphate decarboxylase-catalyzed reactions

The side chain cation of Arg235 provides a 5.6 and 2.6 kcal/mol stabilization of the transition states for orotidine 5'-monophosphate (OMP) decarboxylase (OMPDC) from Saccharomyces cerevisiae catalyzed reactions of OMP and 5-fluoroorotidine 5'-monophosphate (FOMP), respectively, a 7.2 kcal/mol stabilization of the vinyl carbanion-like transition state for enzyme-catalyzed exchange of the C-6 proton of 5-fluorouridine 5'-monophosphate (FUMP), but no stabilization of the transition states for enzyme-catalyzed decarboxylation of truncated substrates 1-(β-d-erythrofuranosyl)orotic acid and 1-(β-d-erythrofuranosyl) 5-fluorouracil. These observations show that the transition state stabilization results from formation of a protein cation-phosphodianion pair, and that there is no detectable stabilization from an interaction between the side chain and the pyrimidine ring of substrate. The 5.6 kcal/mol side chain interaction with the transition state for the decarboxylation reaction is 50% of the total 11.2 kcal/mol transition state stabilization by interactions with the phosphodianion of OMP, whereas the 7.2 kcal/mol side chain interaction with the transition state for the deuterium exchange reaction is a larger 78% of the total 9.2 kcal/mol transition state stabilization by interactions with the phosphodianion of FUMP. The effect of the R235A mutation on the enzyme-catalyzed deuterium exchange is expressed predominantly as a change in the turnover number kex, whereas the effect on the enzyme-catalyzed decarboxylation of OMP is expressed predominantly as a change in the Michaelis constant Km. These results are rationalized by a mechanism in which the binding of OMP, compared with that for FUMP, provides a larger driving force for conversion of OMPDC from an inactive open conformation to a productive, active, closed conformation.

Stabilities of Uracil and Pyridone-Based Carbanions: A Systematic Study in the Gas Phase and Solution and Implications for the Mechanism of Orotidine-5'-Monophosphate Decarboxylase

The stabilities of the C6-centered carbanions derived from 1,3-dimethyluracil, N-methyl-2-pyridone, and N-methyl-4-pyridone were systematically investigated in the gas phase and in DMSO and water solutions. The stabilities of the carbanions in the gas phase and DMSO were directly measured through their reactions with carbon acids with known proton affinity or pK_a values. The stabilities of the carbanions in DMSO were also probed through their kinetic isotope effects of protonation over deuteriation using acids with different acidity. The stabilities of the carbanions in water were determined through the rates of hydrogen-deuterium exchange reactions of the corresponding conjugate acids. The carbanions derived from the two pyridones were found to have the same stability, whereas the carbanion derived from 1,3-dimethyluracil was more stable. The order of the stability of the carbanions showed no correlation with the decarboxylation rates of their corresponding carboxylic acids. The implications of the results for the mechanism of orotidine-5'-monophosphate decarboxylase (ODCase) are discussed.

Quantum and classical simulations of orotidine monophosphate decarboxylase: support for a direct decarboxylation mechanism

Orotidine 5'-monophosphate (OMP) decarboxylase (ODCase) catalyzes the decarboxylation of OMP to uridine 5'-monophosphate (UMP). Numerous studies of this reaction have suggested a plethora of mechanisms including covalent addition, ylide or carbene formation, and concerted or stepwise protonation. Recent experiments and simulations present strong evidence for a direct decarboxylation mechanism, although direct comparison between experiment and theory is still lacking. In the current work we present hybrid quantum mechanics-molecular mechanics simulations that address the detailed decarboxylation mechanisms for OMP and 5-fluoro-OMP by ODCase. Multidimensional potentials of mean force are computed as functions of structural progress coordinates for the Methanobacterium thermoautotrophicum ODCase reaction: the decarboxylation reaction coordinate, an orbital rehybridization coordinate, and the proton transfer coordinate between Lys72 and the substrate. The computed free energy profiles are in accord with the available experimental data. To facilitate further direct comparison with experiment, we compute the kinetic isotope effects (KIEs) for the enzyme-catalyzed reactions using a mass-perturbation-based path-integral method. The computed KIE provide further support for a direct decarboxylation mechanism. In agreement with experiment, the data suggest a role for Lys72 in stabilizing the transition state in the catalysis of OMP and, to a somewhat lesser extent, in 5-fluoro-OMP.

Using catalytic atom maps to predict the catalytic functions present in enzyme active sites

Catalytic atom maps (CAMs) are minimal models of enzyme active sites. The structures in the Protein Data Bank (PDB) were examined to determine if proteins with CAM-like geometries in their active sites all share the same catalytic function. We combined the CAM-based search protocol with a filter based on the weighted contact number (WCN) of the catalytic residues, a measure of the "crowdedness" of the microenvironment around a protein residue. Using this technique, a CAM based on the Ser-His-Asp catalytic triad of trypsin was able to correctly identify catalytic triads in other enzymes within 0.5 Å rmsd of the CAM with 96% accuracy. A CAM based on the Cys-Arg-(Asp/Glu) active site residues from the tyrosine phosphatase active site achieved 89% accuracy in identifying this type of catalytic functionality. Both of these CAMs were able to identify active sites across different fold types. Finally, the PDB was searched to locate proteins with catalytic functionality similar to that present in the active site of orotidine 5'-monophosphate decarboxylase (ODCase), whose mechanism is not known with certainty. A CAM, based on the conserved Lys-Asp-Lys-Asp tetrad in the ODCase active site, was used to search the PDB for enzymes with similar active sites. The ODCase active site has a geometry similar to that of Schiff base-forming Class I aldolases, with lowest aldolase rmsd to the ODCase CAM at 0.48 Å. The similarity between this CAM and the aldolase active site suggests that ODCase has the correct catalytic functionality present in its active site for the generation of a nucleophilic lysine.

Chemical models and their mechanistic implications for the transformation of 6-cyanouridine 5'-monophosphate catalyzed by orotidine 5'-monophosphate decarboxylase

The reactions of 6-cyano-1,3-dimethyluracil have been studied as chemical models to illustrate the mechanism for the transformation of 6-cyanouridine 5'-monophosphate (6-CN-UMP) to barbiturate ribonucleoside 5'-monophosphate (BMP) catalyzed by orotidine 5'-monophosphate decarboxylase (ODCase). The results suggest that the Asp residue in the ODCase active site plays the role of a general base in the transformation.

An examination of the relationship between active site loop size and thermodynamic activation parameters for orotidine 5'-monophosphate decarboxylase from mesophilic and thermophilic organisms

Closure of the active site phosphate gripper loop of orotidine 5'-monophosphate decarboxylase from Saccharomyces cerevisiae (ScOMPDC) over the bound substrate orotidine 5'-monophosphate (OMP) activates the bound substrate for decarboxylation by at least 10(4)-fold [Amyes, T. L., Richard, J. P., and Tait, J. J. (2005) J. Am. Chem. Soc. 127, 15708-15709]. The 19-residue phosphate gripper loop of the mesophilic ScOMPDC is much larger than the nine-residue loop at the ortholog from the thermophile Methanothermobacter thermautotrophicus (MtOMPDC). This difference in loop size results in a small decrease in the total intrinsic phosphate binding energy of the phosphodianion group of OMP from 11.9 to 11.6 kcal/mol, along with a modest decrease in the extent of activation by phosphite dianion of decarboxylation of the truncated substrate 1-(beta-D-erythrofuranosyl)orotic acid. The activation parameters DeltaH(double dagger) and DeltaS(double dagger) for k(cat) for decarboxylation of OMP are 3.6 kcal/mol and 10 cal K(-1) mol(-1) more positive, respectively, for MtOMPDC than for ScOMPDC. We suggest that these differences are related to the difference in the size of the active site loops at the mesophilic ScOMPDC and the thermophilic MtOMPDC. The greater enthalpic transition state stabilization available from the more extensive loop-substrate interactions for the ScOMPDC-catalyzed reaction is largely balanced by a larger entropic requirement for immobilization of the larger loop at this enzyme.

Mechanism of the orotidine 5'-monophosphate decarboxylase-catalyzed reaction: effect of solvent viscosity on kinetic constants

Orotidine 5'-monophosphate decarboxylase (OMPDC) is an exceptionally proficient catalyst: the rate acceleration (k(cat)/k(non)) is 7.1 x 10(16), and the proficiency [(k(cat)/K(M))/k(non)] is 4.8 x 10(22) M(-1). The structural basis for the large rate acceleration and proficiency is unknown, although the mechanism has been established to involve a stabilized carbanion intermediate. To provide reaction coordinate context for interpretation of the values of k(cat), k(cat)/K(M), and kinetic isotope effects, we investigated the effect of solvent viscosity on k(cat) and k(cat)/K(M) for the OMPDCs from Methanothermobacter thermautotrophicus (MtOMPDC) and Saccharomyces cerevisiae (ScOMPDC). For MtOMPDC, we used not only the natural OMP substrate but also a catalytically impaired mutant (D70N) and a more reactive substrate (FOMP); for ScOMPDC, we used OMP and FOMP. With MtOMPDC and OMP, k(cat) is independent of solvent viscosity, indicating that decarboxylation is fully rate-determining; k(cat)/K(M) displays a fractional dependence of solvent viscosity, suggesting that both substrate binding and decarboxylation determine this kinetic constant. For ScOMPDC with OMP, we observed that both k(cat) and k(cat)/K(M) are fractionally dependent on solvent viscosity, suggesting that the rates of substrate binding, decarboxylation, and product dissociation are similar. Consistent with these interpretations, for both enzymes with FOMP, the increases in the values of k(cat) and k(cat)/K(M) are much less than expected based on the ability of the 5-fluoro substituent to stabilize the anionic intermediate; i.e., substrate binding and product dissociation mask the kinetic effects of stabilization of the intermediate by the substituent.

Lys314 is a nucleophile in non-classical reactions of orotidine-5'-monophosphate decarboxylase

Orotidine-5'-monophosphate decarboxylase (OMPD) catalyzes the decarboxylation of orotidine-5'-monophosphate (OMP) to uridine-5'-monophosphate (UMP) in an extremely proficient manner. The reaction does not require any cofactors and proceeds by an unknown mechanism. In addition to decarboxylation, OMPD is able to catalyze other reactions. We show that several C6-substituted UMP derivatives undergo hydrolysis or substitution reactions that depend on a lysine residue (Lys314) in the OMPD active site. 6-Cyano-UMP is converted to UMP, and UMP derivatives with good leaving groups inhibit OMPD by a suicide mechanism in which Lys314 covalently binds to the substrate. These non-classical reactivities of human OMPD were characterized by cocrystallization and freeze-trapping experiments with wild-type OMPD and two active-site mutants by using substrate and inhibitor nucleotides. The structures show that the C6-substituents are not coplanar with the pyrimidine ring. The extent of this substrate distortion is a function of the substituent geometry. Structure-based mechanisms for the reaction of 6-substituted UMP derivatives are extracted in accordance with results from mutagenesis, mass spectrometry, and OMPD enzyme activity. The Lys314-based mechanisms explain the chemodiversity of OMPD, and offer a strategy to design mechanism-based inhibitors that could be used for antineoplastic purposes for example.

Structural characterization of the molecular events during a slow substrate-product transition in orotidine 5'-monophosphate decarboxylase

Crystal structures of substrate-product complexes of Methanobacterium thermoautotrophicum orotidine 5'-monophosphate decarboxylase, obtained at various steps in its catalysis of the unusual transformation of 6-cyano-uridine 5'-monophosphate (UMP) into barbituric acid ribosyl monophosphate, show that the cyano substituent of the substrate, when bound to the active site, is first bent significantly from the plane of the pyrimidine ring and then replaced by an oxygen atom. Although the K72A and D70A/K72A mutants are either catalytically impaired or even completely inactive, they still display bending of the C6 substituent. Interestingly, high-resolution structures of the D70A and D75N mutants revealed a covalent bond between C6 of UMP and the Lys72 side chain after the -CN moiety's release. The same covalent bond was observed when the native enzyme was incubated with 6-azido-UMP and 6-iodo-UMP; in contrast, the K72A mutant transformed 6-iodo-UMP to barbituric acid ribosyl 5'-monophosphate. These results demonstrate that, given a suitable environment, native orotidine 5'-monophosphate decarboxylase and several of its mutants are not restricted to the physiologically relevant decarboxylation; they are able to catalyze even nucleophilic substitution reactions but consistently maintain distortion on the C6 substituent as an important feature of catalysis.

Development and application of ab initio QM/MM methods for mechanistic simulation of reactions in solution and in enzymes

Determining the free energies and mechanisms of chemical reactions in solution and enzymes is a major challenge. For such complex reaction processes, combined quantum mechanics/molecular mechanics (QM/MM) method is the most effective simulation method to provide an accurate and efficient theoretical description of the molecular system. The computational costs of ab initio QM methods, however, have limited the application of ab initio QM/MM methods. Recent advances in ab initio QM/MM methods allowed the accurate simulation of the free energies for reactions in solution and in enzymes and thus paved the way for broader application of the ab initio QM/MM methods. We review here the theoretical developments and applications of the ab initio QM/MM methods, focusing on the determination of reaction path and the free energies of the reaction processes in solution and enzymes.

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.

Structures of the human orotidine-5'-monophosphate decarboxylase support a covalent mechanism and provide a framework for drug design

UMP synthase (UMPS) catalyzes the last two steps of de novo pyrimidine nucleotide synthesis and is a potential cancer drug target. The C-terminal domain of UMPS is orotidine-5'-monophosphate decarboxylase (OMPD), a cofactor-less yet extremely efficient enzyme. Studies of OMPDs from micro-organisms led to the proposal of several noncovalent decarboxylation mechanisms via high-energy intermediates. We describe nine crystal structures of human OMPD in complex with substrate, product, and nucleotide inhibitors. Unexpectedly, simple compounds can replace the natural nucleotides and induce a closed conformation of OMPD, defining a tripartite catalytic site. The structures outline the requirements drugs must meet to maximize therapeutic effects and minimize cross-species activity. Chemical mimicry by iodide identified a CO(2) product binding site. Plasticity of catalytic residues and a covalent OMPD-UMP complex prompt a reevaluation of the prevailing decarboxylation mechanism in favor of covalent intermediates. This mechanism can also explain the observed catalytic promiscuity of OMPD.

QM/MM metadynamics study of the direct decarboxylation mechanism for orotidine-5'-monophosphate decarboxylase using two different QM regions: acceleration too small to explain rate of enzyme catalysis

Despite decades of study, the mechanism by which orotidine-5'-monophosphate decarboxylase (ODCase) catalyzes the decarboxylation of orotidine monophosphate remains unresolved. A computational investigation of the direct decarboxylation mechanism has been performed using mixed quantum mechanical/molecular mechanical (QM/MM) dynamics simulations. The study was performed with the program CP2K that integrates classical dynamics and ab initio dynamics based on the Born-Oppenheimer approach. Two different QM regions were explored. The free energy barriers for direct decarboxylation of orotidine-5'-monophosphate (OMP) in solution and in the enzyme (using the larger QM region) were determined with the metadynamics method to be 40 and 33 kcal/mol, respectively. The calculated change in activation free energy (DeltaDeltaG++) on going from solution to the enzyme is therefore -7 kcal/mol, far less than the experimental change of -23 kcal/ mol (for k(cat.)/k(uncat.): Radzicka, A.; Wolfenden, R., Science 1995, 267, 90-92). These results do not support the direct decarboxylation mechanism that has been proposed for the enzyme. However, in the context of QM/MM calculations, it was found that the size of the QM region has a dramatic effect on the calculated reaction barrier.

Stability of the 6-carbanion of uracil analogues: mechanistic implications for model reactions of orotidine-5'-monophosphate decarboxylase

[Structure: see text] The pKa's of the 6-CH groups of N-methyl-2-pyridone and N-methyl-4-pyridone in aqueous solution were determined. No correlation between the stability of the carbanions and the rate of decarboxylation of the corresponding carboxylic acids was found.

Carbanions from decarboxylation of orotate analogs: stability and mechanistic implications

The pKa's of the 6-CH groups of 1,3-dimethyluracil, N-methyl-2-pyridone, and N-methyl-4-pyridone were determined through their reactions with bases derived from carbon acids with known pKa and the reactions of their corresponding carbanions with the carbon acids. No correlation between the stability of the carbanions and the rate of decarboxylation of corresponding carboxylic acids was found.

Theoretical Studies of the Effect of Thio Substitution on Orotidine Monophosphate Decarboxylase Substrates

[reaction: see text] The effect of replacing carbonyl oxygens with sulfur in a series of orotidine 5'-monophosphate decarboxylase (ODCase) substrates was studied computationally. Previous experimental results indicate that while 2-thio-orotidine 5'-monophosphate (2-thio-OMP) is a poor substrate for ODCase, 4-thio-orotidine 5'-monophosphate (4-thio-OMP) binds to ODCase, and the resultant k(cat) is measurable. Energetics calculations on 2-thio-1-methyl-orotate and 4-thio-1-methyl-orotate (as models for the 2- and 4-thio-OMPs) indicate that mechanisms involving proton transfer to the 2- or 4-site, regardless of substrate and regardless of whether the 2- or 4-position is a carbonyl or thiocarbonyl, are energetically favorable, as compared to direct decarboxylation without proton transfer. Proton transfer to the 4-site during decarboxylation is found to be energetically more favorable than 2-protonation. Each thiocarbonyl is also found to be more basic than its carbonyl counterpart. Therefore, if 2- or 4-proton transfer is the operative catalytic pathway, energetics alone would not explain why 2-thio-orotidine 5'-monophosphate is a poor ODCase substrate. Conformational preferences for a series of ODCase substrates were also examined computationally. Specifically, the energies and Boltzmann probabilities of the conformers resulting from rotation about the C1'-N1 bond (O4'-C1'-N1-C2 rotation from 0 degrees to 360 degrees ) were calculated. It was found that a calculated preference for the syn versus the anti nucleoside conformation correlates to an experimentally better substrate: the OMP and 4-thio-OMP models show a preference for syn conformations, whereas the 2-thio-OMP (the only substrate of the three OMPs that is experimentally found to bind poorly) model shows a preference for an anti conformation. The same rough correlation was found for a series of ODCase inhibitors; that is, a preference for the syn conformation correlates to a better inhibitor. This result is of interest and points to the possibility that the ability for a substrate to bind well to ODCase may be related to its tendency to favor the syn conformation.

Modest catalysis of the decarboxylation of orotate by hydrogen bonding: a theoretical model for orotidine- 5' -monophosphate decarboxylase

As a model for interactions present in the active site of orotidine-5'-monophosphate decarboxylase (ODCase), the effect of hydrogen bonds to the carbonyl groups (O-2 and O-4) of orotic acid and its decarboxylation product was probed with ab initio calculations. We have found that the transition state/carbanion intermediate is a better proton receptor and therefore, the hydrogen bonds can be a modest source of catalysis. Comparison of the calculated data with results from site-directed mutagenesis provides some insights into the polarity of the active site.

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.

The acidity of uracil and uracil analogs in the gas phase: four surprisingly acidic sites and biological implications

The gas phase acidities of a series of uracil derivatives (1-methyluracil, 3-methyluracil, 6-methyluracil, 5,6-dimethyluracil, and 1,3-dimethyluracil) have been bracketed to provide an understanding of the intrinsic reactivity of uracil. The experiments indicate that in the gas phase, uracil has four sites more acidic than water. Among the uracil analogs, the N1-H sites have deltaH(acid) values of 331-333 kcal mol(-1); the acidity of the N3 sites fall between 347-352 kcal mol(-1). The vinylic C6 in 1-methyluracil and 3-methyluracil brackets to 363 kcal mol(-1), and 369 kcal mol(-1) in 1,3-dimethyluracil; the C5 of 1,3-dimethyluracil brackets to 384 kcal mol(-1). Calculations conducted at B3LYP/6-31+G* are in agreement with the experimental values. The bracketing of several of these sites involved utilization of an FTMS protocol to measure the less acidic site in a molecule that has more than one acidic site, establishing the generality of this method. In molecules with multiple acidic sites, only the two most acidic sites were bracketable, which is attributable to a kinetic effect. The measured acidities are in direct contrast to in solution, where the two most acidic sites of uracil (N1 and N3) are indifferentiable. The vinylic C6 site is also particularly acidic, compared to acrolein and pyridine. The biological implications of these results, particularly with respect to enzymes for which uracil is a substrate, are discussed.

Substrate binding induces domain movements in orotidine 5'-monophosphate decarboxylase

Orotidine 5'-monophosphate decarboxylase (ODCase) catalyses the decarboxylation of orotidine 5'-monophosphate to uridine 5'-monophosphate (UMP). We have earlier determined the structure of ODCase from Escherichia coli complexed with the inhibitor 1-(5'-phospho-beta-d-ribofuranosyl)barbituric acid (BMP); here we present the 2.5 A structure of the uncomplexed apo enzyme, determined from twinned crystals. A structural analysis and comparison of the two structures of the E. coli enzyme show that binding of the inhibitor is accompanied by significant domain movements of approximately 12 degrees around a hinge that crosses the active site. Hence, the ODCase dimer, which contains two active sites, may be divided in three domains: a central domain that is fixed, and two lids which independently move 12 degrees upon binding. Corresponding analyses, presented herein, of the two Saccharomyces cerevisiae ODCase structures (with and without BMP) and the Methanobacterium thermoautotrophicum ODCase structures (with and without 6-aza UMP) show very similar, but somewhat smaller domain movements. The domain movements seem to be initiated by the phosphoryl binding to the enzyme and can explain why the binding of the phosphoryl group is essential for the catalytic function.

An alternative explanation for the catalytic proficiency of orotidine 5'-phosphate decarboxylase

Orotidine 5'-phosphate decarboxylase (ODCase) is the most proficient enzyme known, enhancing the rate of decarboxylation of orotidine 5'-phosphate (OMP) by a factor of 10(17), which corresponds to a DeltaDeltaG++ of approximately 24 kcal/mol. Ground-state destabilization through local electrostatic stress has been recently proposed as the basis of catalytic rate enhancement for a mechanism that is the same as in solution. We have carried out gas-phase ab initio quantum mechanical calculations combined with a free energy method, a continuum solvent model, and molecular dynamics simulations to assess an alternative mechanism. Although we are not able to reproduce the experimentally observed DeltaDeltaG++ quantitatively, we present evidence that this DeltaDeltaG++ is very large, in the range found experimentally. We thus conclude that the preferred mechanism may well be different from that in solution, involving an equilibrium pre-protonation of OMP C5 by a catalytic lysine residue that greatly reduces the barrier to subsequent decarboxylation.

Theoretical studies of mechanisms and kinetic isotope effects on the decarboxylation of orotic acid and derivatives

The mechanism of orotidine 5'-monophosphate decarboxylase was studied computationally by using the decarboxylation of orotic acid analogues as model systems. These calculations indicate that mechanisms involving proton transfer to the 2-oxygen or the 4-oxygen are energetically favorable, as compared to direct decarboxylation without proton transfer, for a series of model compounds where N1 is substituted with respectively H, CH(3), and a tetrahydrofuran moiety. Proton transfer to the 4-oxygen during decarboxylation is found to be energetically more favorable than 2-protonation, which is attributable to both the 4-oxygen site being more basic and an apparent intrinsic preference for the 4-protonation pathway. (15)N isotope effect calculations were also conducted, and compared to experimental (15)N isotope effects previously measured at N1 by Rishavy and Cleland (Biochemistry 2000, 39, 4569-4574). The theoretical isotope effects establish, for the first time, that the experimental (15)N isotope effect is consistent with decarboxylation without protonation, as well as with decarboxylation with protonation, at either O2 or at O4. Furthermore, we propose herein an isotope measurement that could potentially distinguish among mechanisms involving protonation from those that do not involve proton transfer.

The structural basis for the remarkable catalytic proficiency of orotidine 5'-monophosphate decarboxylase

The three-dimensional structures of orotidine 5'-monophosphate decarboxylases from four different organisms have been determined by X-ray crystallography. The structures reveal an active site in which the pyrimidine base and phosphate groups are rigidly held in place. Surprisingly, both pyrimidine carbonyl groups are hydrogen bonded to amide groups, rather than to strong active site acids, as was previously predicted. The positioning of a conserved aspartate sidechain close to the substrate carboxylate and a conserved lysine ammonium group close to the C6 of the pyrimidine suggests a novel mechanism to explain the extreme catalytic proficiency of this enzyme.

Structural basis for the catalytic mechanism of a proficient enzyme: orotidine 5'-monophosphate decarboxylase

Orotidine 5'-monophosphate decarboxylase (ODCase) catalyzes the decarboxylation of orotidine 5'-monophosphate, the last step in the de novo synthesis of uridine 5'-monophosphate. ODCase is a very proficient enzyme [Radzicka, A., and Wolfenden, R. (1995) Science 267, 90-93], enhancing the reaction rate by a factor of 10(17). This proficiency has been enigmatic, since it is achieved without metal ions or cofactors. Here we present a 2.5 A resolution structure of ODCase complexed with the inhibitor 1-(5'-phospho-beta-D-ribofuranosyl)barbituric acid. It shows a closely packed dimer composed of two alpha/beta-barrels with two shared active sites. The orientation of the orotate moiety of the substrate is unambiguously deduced from the structure, and previously proposed catalytic mechanisms involving protonation of O2 or O4 can be ruled out. The proximity of the OMP carboxylate group with Asp71 appears to be instrumental for the decarboxylation of OMP, either through charge repulsion or through the formation of a very short O.H.O hydrogen bond between the two carboxylate groups.

The crystal structure and mechanism of orotidine 5'-monophosphate decarboxylase

The crystal structure of Bacillus subtilis orotidine 5'-monophosphate (OMP) decarboxylase with bound uridine 5'-monophosphate has been determined by multiple wavelength anomalous diffraction phasing techniques and refined to an R-factor of 19.3% at 2.4 A resolution. OMP decarboxylase is a dimer of two identical subunits. Each monomer consists of a triosephosphate isomerase barrel and contains an active site that is located across one end of the barrel and near the dimer interface. For each active site, most of the residues are contributed by one monomer with a few residues contributed from the adjacent monomer. The most highly conserved residues are located in the active site and suggest a novel catalytic mechanism for decarboxylation that is different from any previously proposed OMP decarboxylase mechanism. The uridine 5'-monophosphate molecule is bound to the active site such that the phosphate group is most exposed and the C5-C6 edge of the pyrimidine base is most buried. In the proposed catalytic mechanism, the ground state of the substrate is destabilized by electrostatic repulsion between the carboxylate of the substrate and the carboxylate of Asp60. This repulsion is reduced in the transition state by shifting negative charge from the carboxylate to C6 of the pyrimidine, which is close to the protonated amine of Lys62. We propose that the decarboxylation of OMP proceeds by an electrophilic substitution mechanism in which decarboxylation and carbon-carbon bond protonation by Lys62 occur in a concerted reaction.

The mechanism of orotidine 5'-monophosphate decarboxylase: catalysis by destabilization of the substrate

The mechanism of orotidine 5'-monophosphate decarboxylase (OMP decarboxylase, ODCase) was studied using the decarboxylation of orotic acid analogues as a model system. The rate of decarboxylation of 1,3-dimethylorotic acid and its analogues as well as the stability of their corresponding carbanion intermediates was determined. The results have shown that the stability of the carbanion intermediate is not a critical factor in the rate of decarboxylation. On the other hand, the reaction rate is largely dependent on the equilibrium constant for the formation of a zwitterion. Based on these results, we have proposed a new mechanism in which ODCase catalyzes the decarboxylation of OMP by binding the substrate in a zwitterionic form and providing a destabilizing environment for the carboxylate group of OMP.

No metal cofactor in orotidine 5'-monophosphate decarboxylase

Orotidine 5'-monophosphate decarboxylase (OMP decarboxylase, ODCase) is an important enzyme that catalyzes the final step of de novo pyrimidine nucleotide biosynthesis. The mechanism of this unique enzyme and whether metal ions play any role in catalysis have been topics of intense research interest. In this report, the role of Zn in ODCase was reexamined. Atomic absorption (AA) and X-ray absorption (XAS) spectroscopic studies did not detect zinc in active enzyme samples at high concentration. The XAS results also indicated the absence of other transition metal ions in ODCase.

A unique catalytic and inhibitor-binding role for Lys93 of yeast orotidylate decarboxylase

The presence of a proton-donating catalytic amino acid side chain in orotidylate decarboxylase (ODCase) was sought by site-directed mutagenesis. Replacement of yeast ODCase Lys93 with a cysteine resulted in a mutant protein (K93C) with no measurable activity, representing a decrease in activity by a factor of, at most, 2 x 10(-8) times the activity of the wild-type enzyme. Treatment of this mutant protein with 2-bromoethylamine, designed to append Cys93 to yield S-(2-aminoethyl)cysteine, restored activity by a factor of at least 5 x 10(5) over the untreated mutant protein. Activity could not be restored by treatment with other brominated reagents designed to replace the epsilon-amino of S-(2-aminoethyl)Cys93 with a different functional group. The overall architecture of the K93C protein was not significantly changed, as judged by the similar dimerization properties (in the absence of ligands) of the mutant enzyme compared to the wild-type enzyme. The binding affinity of the substrate orotidylate was not measurably changed by the mutation, indicating that Lys93 has an essential role in catalysis which is mechanistically distinguishable from substrate binding. Apparently the mutation removes an integral portion of the active site and does not drastically affect the structural or substrate binding properties. However, the affinities of the mutant protein for the competitive inhibitors 6-azauridylate (6-azaUMP) and UMP are significantly altered from the pattern seen with the wild-type enzyme. The K93C protein has an affinity for the neutral ligand UMP which is greater than that for the anionic 6-azaUMP, in clear contrast to the preference for 6-azaUMP displayed by the wild-type enzyme. Lys93 is apparently critical for catalysis of the substrate to product and for the binding of anionic inhibitors; the data are discussed in terms of previously existing models for transition-state analogue inhibitor binding and catalysis.

Orotidylate decarboxylase: insights into the catalytic mechanism from substrate specificity studies

Pyrimidine nucleotides were tested as substrates for pure yeast orotidylate decarboxylase in an attempt to gain insight into the nature of the catalytic mechanism of the enzyme. Substitutions of the 5-position in the pyrimidine ring of the orotidylate substrate resulted in compounds that are either excellent inhibitors or substrates of the enzyme. The 5-bromo- and 5-chloroorotidylates are potent inhibitors while the 5-fluoro derivative is a good substrate with a turnover number 30 times that observed with orotidylate. When carbon 5 of the pyrimidine ring is replaced by nitrogen in 5-azaorotidylate, the resulting compound is unstable in solution with a half-life of 25 min at pH 6. However, studies with freshly generated 5-azaorotidylate show that an enzyme-dependent reaction occurs, presumably decarboxylation. This enzyme reaction follows simple Michaelis-Menten kinetics. Because the 5-aza group is not electrophilic, an enzyme mechanism utilizing a nucleophilic addition of the enzyme at the 5-position is ruled out. We also present studies that are not compatible with a mechanism requiring the formation of a Schiff's base prior to decarboxylation. The enzyme is tolerant of modest substitution at the 4-position, for the 4-keto group can be replaced with a thioketone. However, no catalysis is observed when the same substitution is made at the 2-position. Similarities in the substrate specificity of orotate phosphoribosyltransferase and orotidylate decarboxylase led us to compare the amino acid sequences of the two enzymes; significant (20%) sequence homology was observed.(ABSTRACT TRUNCATED AT 250 WORDS)

Orotidylate decarboxylase of yeast and man

The mechanism for ODCase appears to involve the formation of a zwitterion of OMP and a ylid on decarboxylation. Thiamin pyrophosphate catalyzes various decarboxylation and transfer reactions involving ketone groups because the thiazolium ring with its positively charged N atom can, on the loss of a proton from the adjacent C-2, generate a ylid which adds to carbonyl groups to produce a substrate ylid. The unusual aspect, then, of the ODCase reaction is that the substrate itself becomes the ylid, presumably by gaining a proton from ODCase, which results in a positive charge on the N-1 atom of the pyrimidine ring. It is a zwitterion in the transition state which momentarily becomes a ylid on decarboxylation of OMP which then yields the product, UMP. There is no known cofactor for the ODCase reaction. It will be of interest to discover the groups on the enzyme that aid in formation of the zwitterion and the ylid. Further work on the crystal structure and on the production of altered enzymes (where specific amino acids suspected to be important for the reaction are changed) should reveal more details about this important and novel reaction.

Investigation of the enzymatic mechanism of yeast orotidine-5'-monophosphate decarboxylase using 13C kinetic isotope effects

Orotidine-5'-monophosphate decarboxylase (ODCase) from Saccharomyces cerevisiae displays an observed 13C kinetic isotope effect of 1.0247 +/- 0.0008 at 25 degrees C, pH 6.8. The observed isotope effect is sensitive to changes in the reaction medium, such as pH, temperature, or glycerol content. The value of 1.0494 +/- 0.0006 measured at pH 4.0, 25 degrees C, is not altered significantly by temperature or glycerol, and thus the intrinsic isotope effect for the reaction is apparently being observed under these conditions and decarboxylation is almost entirely rate-determining. These data require a catalytic mechanism with freely reversible binding and one in which a very limited contribution to the overall rate is made by chemical steps preceding decarboxylation; the zwitterion mechanism of Beak and Siegel [Beak, P. & Siegel, B. (1976) J. Am. Chem. Soc. 98, 3601-3606], which involves only protonation of the pyrimidine ring, is such a mechanism. With use of an intrinsic isotope effect of 1.05, a partitioning factor of less than unity is calculated for ODCase at pH 6.0, 25 degrees C. A quantitative kinetic analysis using this result excludes the possibility of an enzymatic mechanism involving covalent attachment of an enzyme nucleophile to C-5 of the pyrimidine ring. The observed isotope effect does not rise to the intrinsic value above pH 8.5; instead, the observed isotope effects at 25 degrees C plotted against pH yield an asymmetric curve that at high pH plateaus at about 1.035. These data, in conjunction with the pH profile of Vmax/km, fit a kinetic model in which an enzyme proton necessary for catalysis is titrated at high pH, thus providing evidence for the catalytic mechanism of Beak and Siegel (1976).

Crystallization of yeast orotidine 5'-monophosphate decarboxylase complexed with 1-(5'-phospho-beta-D-ribofuranosyl) barbituric acid

Using an incomplete factorial experimental design, we have identified conditions for crystallization of yeast orotidine 5'-monophosphate decarboxylase (ODCase) in an unliganded state and complexed separately to two inhibitors: 6-azauridine 5'-monophosphate (aza-UMP) and 1-(5'-phospho-beta-D-ribofuranosyl) barbituric acid (BMP). Crystals of X-ray diffraction quality have been obtained of yeast ODCase complexed with BMP, a putative transition state analog inhibitor (Ki = 8.8 x 10(-12) M). ODCase:BMP complex crystals with a hexagonal rod habit were grown from a solution initially containing 12 mg/ml ODCase (205 microM dimer) plus 450 microM BMP by microdialysis at 4 degrees C against a mother liquor which consisted of 0.1 M Na-PIPES-acetate (pH 6.4), 37.5 microM BMP, 5 mM mercaptoethanol, 1% polyethylene glycol 400, and 2.3 M ammonium sulfate. Crystals were analyzed using precession photography and were assigned to trigonal space group R32 with unit cell dimensions a = b = 115 A, c = 385 A. The crystal density is 1.245 g/cm3 indicating the presence of two ODCase: BMP complex dimers (118 kDa each) per asymmetric unit with a packing density of 2.08 A3/Da and 41% solvent content. The morphological habit of crystals of the ODCase:BMP complex changed when the initial ammonium sulfate concentration was increased in 0.05 M steps from 2.3 to 2.45 M. All of these crystals diffracted to at least 3.0 A resolution over a period of several weeks at room temperature and are isomorphous.

Orotidine-5'-monophosphate decarboxylase catalysis: kinetic isotope effects and the state of hybridization of a bound transition-state analogue

The enzymatic decarboxylation of orotidine 5'-monophosphate may proceed by an addition-elimination mechanism involving a covalently bound intermediate or by elimination of CO2 to generate a nitrogen ylide. In an attempt to distinguish between these two alternatives, 1-(phosphoribosyl)barbituric acid was synthesized with 13C at the 5-position. Interaction of this potential transition-state analogue inhibitor with yeast orotidine-5'-monophosphate decarboxylase resulted in a small (0.6 ppm) downfield displacement of the C-5 resonance, indicating no rehybridization of the kind that might have been expected to accompany 5,6-addition of an enzyme nucleophile. When the substrate orotidine 5'-monophosphate was synthesized with deuterium at C-5, no significant change in kcat (H/D = 0.99 +/- 0.06) or kcat/KM (H/D = 1.00 +/- 0.06) was found to result, suggesting that C-5 does not undergo significant changes in geometry before or during the step that determines the rate of the catalytic process. These results are consistent with a nitrogen ylide mechanism and offer no support for the intervention of covalently bound intermediates in the catalytic process.