An alternative explanation for the catalytic proficiency 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.

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.

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.

Ground state destabilization from a positioned general base in the ketosteroid isomerase active site

We compared the binding affinities of ground state analogues for bacterial ketosteroid isomerase (KSI) with a wild-type anionic Asp general base and with uncharged Asn and Ala in the general base position to provide a measure of potential ground state destabilization that could arise from the close juxtaposition of the anionic Asp and hydrophobic steroid in the reaction's Michaelis complex. The analogue binding affinity increased ~1 order of magnitude for the Asp38Asn mutation and ~2 orders of magnitude for the Asp38Ala mutation, relative to the affinity with Asp38, for KSI from two sources. The increased level of binding suggests that the abutment of a charged general base and a hydrophobic steroid is modestly destabilizing, relative to a standard state in water, and that this destabilization is relieved in the transition state and intermediate in which the charge on the general base has been neutralized because of proton abstraction. Stronger binding also arose from mutation of Pro39, the residue adjacent to the Asp general base, consistent with an ability of the Asp general base to now reorient to avoid the destabilizing interaction. Consistent with this model, the Pro mutants reduced or eliminated the increased level of binding upon replacement of Asp38 with Asn or Ala. These results, supported by additional structural observations, suggest that ground state destabilization from the negatively charged Asp38 general base provides a modest contribution to KSI catalysis. They also provide a clear illustration of the well-recognized concept that enzymes evolve for catalytic function and not, in general, to maximize ground state binding. This ground state destabilization mechanism may be common to the many enzymes with anionic side chains that deprotonate carbon acids.

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.

Quantum Chemical Modeling of Enzymatic Reactions: The Case of Decarboxylation

We present a systematic study of the decarboxylation step of the enzyme aspartate decarboxylase with the purpose of assessing the quantum chemical cluster approach for modeling this important class of decarboxylase enzymes. Active site models ranging in size from 27 to 220 atoms are designed, and the barrier and reaction energy of this step are evaluated. To model the enzyme surrounding, homogeneous polarizable medium techniques are used with several dielectric constants. The main conclusion is that when the active site model reaches a certain size, the solvation effects from the surroundings saturate. Similar results have previously been obtained from systematic studies of other classes of enzymes, suggesting that they are of a quite general nature.

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.

Structure-activity relationships of orotidine-5'-monophosphate decarboxylase inhibitors as anticancer agents

A series of 6-substituted and 5-fluoro-6-substituted uridine derivatives were synthesized and evaluated for their potential as anticancer agents. The designed molecules were synthesized from either fully protected uridine or the corresponding 5-fluorouridine derivatives. The mononucleotide derivatives were used for enzyme inhibition investigations against ODCase. Anticancer activities of all the synthesized derivatives were evaluated using the nucleoside forms of the inhibitors. 5-Fluoro-UMP was a very weak inhibitor of ODCase. 6-Azido-5-fluoro and 5-fluoro-6-iodo derivatives are covalent inhibitors of ODCase, and the active site Lys145 residue covalently binds to the ligand after the elimination of the 6-substitution. Among the synthesized nucleoside derivatives, 6-azido-5-fluoro, 6-amino-5-fluoro, and 6-carbaldehyde-5-fluoro derivatives showed potent anticancer activities in cell-based assays against various leukemia cell lines. On the basis of the overall profile, 6-azido-5-fluoro and 6-amino-5-fluoro uridine derivatives exhibited potential for further investigations.

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.

A substantial oxygen isotope effect at O2 in the OMP decarboxylase reaction: mechanistic implications

Orotidine-5'-monophosphate decarboxylase (OMP decarboxylase, ODCase) catalyzes the decarboxylation of orotidine-5'-monophosphate (OMP) to uridine-5'-monophosphate (UMP). Despite extensive enzymological, structural, and computational studies, the mechanism of ODCase remains incompletely characterized. Herein, carbon kinetic isotope effects were measured for both the natural abundance substrate and a substrate mixture synthesized for the purpose of carrying out the remote double label isotope effect procedure, with O2 of the substrate as the remote position. The carbon kinetic isotope effect on enzymatic decarboxylation of this substrate mix was measured to be 1.0199 +/- 0.0007, compared to the value of 1.0289 +/- 0.0009 for natural abundance OMP, revealing an (18)O2 isotope effect of 0.991 +/- 0.001. This value equates to an intrinsic isotope effect of approximately 0.983, using a calculated commitment factor derived from previous isotope effect data. The measured (18)O2 isotope effect requires a mechanism with one or more enzymatic processes, including binding and/or chemistry, that contribute to this substantial inverse isotope effect. (18)O2 kinetic isotope effects were calculated for four proposed mechanisms: decarboxylation preceded by proton transfer to 1) O2; 2) O4; and 3) C5; and 4) decarboxylation without a preceding protonation step. A mechanism involving no pre-decarboxylation step does not appear to have any steps with the necessary substantial inverse (18)O2 effect, thus calling into question any mechanism involving simple direct decarboxylation. Protonation at O2, O4, or C5 are all calculated to proceed with inverse (18)O2 effects, and could contribute to the experimentally measured value. Recent crystal structures indicate that O2 of the substrate appears to be involved in an intricate bonding arrangement involving the substrate phosphoryl group, an enzyme Gln side chain, and a bound water molecule; this interaction likely contributes to the observed isotope effect.

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.

Free energies of chemical reactions in solution and in enzymes with ab initio quantum mechanics/molecular mechanics methods

Combined quantum mechanics/molecular mechanics (QM/MM) methods provide an accurate and efficient energetic description of complex chemical and biological systems, leading to significant advances in the understanding of chemical reactions in solution and in enzymes. Here we review progress in QM/MM methodology and applications, focusing on ab initio QM-based approaches. Ab initio QM/MM methods capitalize on the accuracy and reliability of the associated quantum-mechanical approaches, however, at a much higher computational cost compared with semiempirical quantum-mechanical approaches. Thus reaction-path and activation free-energy calculations based on ab initio QM/MM methods encounter unique challenges in simulation timescales and phase-space sampling. This review features recent developments overcoming these challenges and enabling accurate free-energy determination for reaction processes in solution and in enzymes, along with applications.

Indiscriminate binding by orotidine 5'-phosphate decarboxylase of uridine 5'-phosphate derivatives with bulky anionic c6 substituents

Orotidine 5'-phosphate (OMP) decarboxylase appears to act upon its substrate without the intervention of metals or other cofactors and without the formation of covalent bonds between the enzyme and the substrate. Crystallographic information indicates that substrate binding forces the substrate's scissile carboxylate group into the neighborhood of several charged groups at the active site. It has been proposed that binding might result in electrostatic stress at the substrate's C6 carboxylate group in such a way as to promote decarboxylation by destabilizing the enzyme-substrate complex in its ground state. If that were the case, one would expect uridine 5'-phosphate (UMP) derivatives with bulky anionic substituents at C6 to be bound weakly compared with UMP, which is unsubstituted at C6. Here, we describe the formation of anionic 5,6-dihydro-6-sulfonyl derivatives by spontaneous addition of sulfite to UMP and to OMP. These sulfite addition reactions, which are slowly reversible and are not catalyzed by the enzyme, result in the appearance of one (or, in the case of OMP, two) bulky anionic substituents at the 6-carbon atom of UMP. These inhibitors are bound with affinities that surpass the binding affinity of UMP. We are led to infer that the active site of OMP decarboxylase is remarkably accommodating in the neighborhood of C6. These are not the properties that one would expect of an active site with a rigid structure that imposes sufficient electrostatic stress on the substrate to produce a major advancement along the reaction coordinate.

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.

A proficient enzyme: insights on the mechanism of orotidine monophosphate decarboxylase from computer simulations

Decarboxylation of orotidine 5'-monophosphate (Omp) to uridine 5'-monophosphate by orotidine 5'-monophosphate decarboxylase (ODCase) is currently the object of vivid debate. Here, we clarify its enzymatic activity with long time scale classical molecular dynamics and hybrid ab initio Car-Parrinello/molecular mechanics simulations. The lack of structural (experimental) information on the ground state of ODCase/Omp complex is overcome by a careful construction of the model and the analysis of three different strains of the enzyme. We find that the ODCase/substrate complex is characterized by a very stable charged network Omp-Lys-Asp-Lys-Asp, which is incompatible with the previously proposed direct decarboxylation driven by a ground-state destabilization. A direct decarboxylation induced by a transition-state electrostatic stabilization is consistent with our findings. The calculated activation free energy for the direct decarboxylation with the formation of a C6 carboanionic intermediate yields an overall rate enhancement by the enzyme (k(cat)/k(wat) = 3.5 x 10(16)) in agreement with experiments (k(cat)/k(wat) = 1.7 x 10(17)). The decarboxylation is accompanied by the movement of a fully conserved lysine residue toward the developing negative charge at the C6 position.

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.

Direct hydroxide attack is a plausible mechanism for amidase antibody 43C9

Direct hydroxide attack on the scissile carbonyl of the substrate has been suggested as a likely mechanism for esterase antibodies elicited by phosphonate haptens, which mimic the transition states for the alkaline hydrolysis of esters.1 The unique amidase activity of esterase antibody 43C9 has been attributed to nucleophilic attack by an active-site histidine residue.2 Yet, the active site of 43C9 is strikingly similar to those of other esterase antibodies, particularly 17E8. We have carried out quantum mechanical calculations, molecular dynamics simulations, and free energy calculations to assess the mechanism involving direct hydroxide attack for 43C9. Results support this mechanism and suggest that the mechanism is plausible for other antiphosphonate antibodies that catalyze the hydrolysis of (p-nitro)phenyl esters.

Hydrogen isotope tracing in the reaction of orotidine-5'-monophosphate decarboxylase

The mechanism of the enzyme orotidine-5(')-monophosphate decarboxylase (OMP decarboxylase, ODCase) is not fully characterized; some of the proposed mechanisms suggest the possibility of hydrogen rearrangement (shift from C5 to C6 or loss of H5 to solvent) during catalysis. In this study, we sought mechanistic information for the ODCase reaction by examining the extent of hydrogen exchange in the product uridine-5(')-monophosphate, in combination with ODCase, at the H5 and H6 positions. In a subsequent experiment, partially deuterated OMP was prepared, and the extent of 2H5 rearrangement or loss to solvent was examined by integration of 1H nuclear magnetic resonance signals in the substrate and the resulting enzymatically decarboxylated product. The absence of detectable hydrogen exchange in these experiments limits somewhat the possible mechanisms for ODCase catalysis.

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.

Molecular dynamic study of orotidine-5'-monophosphate decarboxylase in ground state and in intermediate state: a role of the 203-218 loop dynamics

Molecular dynamics simulations have been used to derive the structures of ground (orotidine-5'-monophosphate decarboxylase x orotidine 5'-monophosphate; ODC x OMP) and intermediate (ODC x intermediate; ODC x I(-)) states in the ODC-catalyzed decarboxylation of OMP. For comparison, a molecular dynamics simulation of the conformers of OMP dissolved in water was also studied. This structural information is unavailable from present crystal structures. The electrostatic network in the active site around the carboxylate moiety of OMP exhibits remarkable stability. The conformation of enzyme-bound OMP is very similar to the conformation of OMP in water. Thus, the proposed Circe effect mechanism for ODC catalysis is unlikely. Comparison of ground state and intermediate state structures shows that on decarboxylation C6 takes the position of the carboxylate O8. This significant movement of the ligand is accompanied by a placement of the C6 carbanion in the vicinity of the protonated Lys-93 and is enforced by a change of the 203-218 loop from an unstructured form to an ordered beta-hairpin. Previously proposed mechanisms involving protonation at O2, O4, or C5 have in common internal stabilization of the anionic intermediate by conjugation with positive charge on the pyrimidine ring. These mechanisms are not supported because there are no proton sources near O2, O4, and C5. We propose that the stabilization of intermediate ODC x I(-) is achieved by movement of the carbanion toward the external cation Lys-93 on decarboxylation and organization of the 203-218 loop. Because the intermediate and transition state are energetically similar, stabilization of the former decreases the free energy content of the latter.