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

Glycerol 3-Phosphate Dehydrogenase: Role of the Protein Conformational Change in Activation of a Readily Reversible Enzyme-Catalyzed Hydride Transfer Reaction

Kinetic parameters are reported for glycerol 3-phosphate dehydrogenase (GPDH)-catalyzed hydride transfer from the whole substrate glycerol 3-phosphate (G3P) or truncated substrate ethylene glycol (EtG) to NAD, and for activation of the hydride transfer reaction of EtG by phosphite dianion. These kinetic parameters were combined with parameters for enzyme-catalyzed hydride transfer in the microscopic reverse direction to give the reaction equilibrium constants Keq. Hydride transfer from G3P is favored in comparison to EtG because the carbonyl product of the former reaction is stabilized by hyperconjugative electron donation from the -CH2R keto substituent. The kinetic data show that the phosphite dianion provides the same 7.6 ± 0.1 kcal/mol stabilization of the transition states for enzyme-catalyzed reactions in the forward [reduction of NAD by EtG] and reverse [oxidation of NADH by glycolaldehyde] directions. The experimental evidence that supports a role for phosphite dianion in stabilizing the active closed form of the GPDH (EC) relative to the ca. 6 kcal/mol more unstable open form (EO) is summarized.

Structural and functional properties of uridine 5'-monophosphate synthase from Coffea arabica

In higher eukaryotes and plants, the last two sequential steps in the de novo biosynthesis of uridine 5'-monophosphate (UMP) are catalyzed by a bifunctional natural chimeric protein called UMP synthase (UMPS). In higher plants, UMPS consists of two naturally fused enzymes: orotate phosphoribosyltransferase (OPRTase) at N-terminal and orotidine-5'-monophosphate decarboxylase (ODCase) at C-terminal. In this work, we obtained the full functional recombinant protein UMPS from Coffea arabica (CaUMPS) and studied its structure-function relationships. A biochemical and structural characterization of a plant UMPS with its two functional domains is described together with the presentation of the first crystal structure of a plant ODCase at 1.4 Å resolution. The kinetic parameters measured of CaOPRTase and CaODCase domains were comparable to those reported. The crystallographic structure revealed that CaODCase is a dimer that conserves the typical fold observed in other ODCases from prokaryote and eukaryote with a 1-deoxy-ribofuranose-5'-phosphate molecule bound in the active site of one subunit induced a closed conformation. Our results add to the knowledge of one of the key enzymes of the de novo biosynthesis of pyrimidines in plant metabolism and open the door to future applications.

Analyzing the Reaction of Orotidine 5'-Phosphate Decarboxylase as a Way to Examine Some Key Catalytic Proposals

This study analyzes the origin of enzyme catalysis by focusing on the reaction of orotidine 5'-phosphate decarboxylase (ODCase). This reaction involves an enormous catalytic effect of 23 kcal/mol that has been attributed to reactant state destabilization associated with the use of binding energy through the so-called Circe effect. However, our early studies and subsequent key experiments have shown that the presumed effect of the binding energy (namely, the strain exerted by a bond to a phosphate group) does not contribute to the catalysis. In this study, we perform quantitative empirical valence bond calculations that reproduce the catalytic effect of ODCase and the effect of removing the phosphate side chain. The calculations demonstrate that the effect of the phosphate is due to a change in reorganization energy and should not be described as an induced fit effect. Similarly, we show that the overall catalytic effect is due to electrostatic transition state stabilization, which again reflects the smaller reorganization energy in the enzyme than in water. We also elaborate on the problems with the induced fit proposal, including the fact that it does not serve to tell us what the actual origin of the action of the catalytic effect is. In addition to the above points, we use this paper to discuss misconceptions about the meaning of the preorganization effect, as well as other misunderstandings of what is being done in consistent calculations of enzyme catalysis.

Microbial rhodoquinone biosynthesis proceeds via an atypical RquA-catalyzed amino transfer from S-adenosyl-L-methionine to ubiquinone

Rhodoquinone (RQ) is a close analogue of ubiquinone (UQ) that confers diverse bacterial and eukaryotic taxa the ability to utilize fumarate as an electron acceptor in hypoxic conditions. The RquA protein, identified in a Rhodospirillum rubrum RQ-deficient mutant, has been shown to be required for RQ biosynthesis in bacteria. In this report, we demonstrate that RquA, homologous to SAM-dependent methyltransferases, is necessary and sufficient to catalyze RQ biosynthesis from UQ in vitro. Remarkably, we show that RquA uses SAM as the amino group donor in a substitution reaction that converts UQ to RQ. In contrast to known aminotransferases, RquA does not use pyridoxal 5'-phosphate (PLP) as a coenzyme, but requires the presence of Mn²⁺ as a cofactor. As these findings reveal, RquA provides an example of a non-canonical SAM-dependent enzyme that does not catalyze methyl transfer, instead it uses SAM in an atypical amino transfer mechanism.

Role of pH in Regulating Cancer Pyrimidine Synthesis

Replication is a fundamental aspect of cancer, and replication is about reproducing all the elements and structures that form a cell. Among them are DNA, RNA, enzymes, and coenzymes. All the DNA is doubled during each S (synthesis) cell cycle phase. This means that six billion nucleic acids must be synthesized in each cycle. Tumor growth, proliferation, and mutations all depend on this synthesis. Cancer cells require a constant supply of nucleotides and other macromolecules. For this reason, they must stimulate de novo nucleotide synthesis to support nucleic acid provision. When deregulated, de novo nucleic acid synthesis is controlled by oncogenes and tumor suppressor genes that enable increased synthesis and cell proliferation. Furthermore, cell duplication must be achieved swiftly (in a few hours) and in the midst of a nutrient-depleted and hypoxic environment. This also means that the enzymes participating in nucleic acid synthesis must work efficiently. pH is a critical factor in enzymatic efficiency and speed. This review will show that the enzymatic machinery working in nucleic acid synthesis requires a pH on the alkaline side in most cases. This coincides with many other pro-tumoral factors, such as the glycolytic phenotype, benefiting from an increased intracellular pH. An increased intracellular pH is a perfect milieu for high de novo nucleic acid production through optimal enzymatic performance.

Implications for an imidazol-2-yl carbene intermediate in the rhodanase-catalyzed C-S bond formation reaction of anaerobic ergothioneine biosynthesis

In the anaerobic ergothioneine biosynthetic pathway, a rhodanese domain containing enzyme (EanB) activates tne hercynine's sp² ε-C-H Dona ana replaces it with a C-S bond to produce ergothioneine. The key intermediate for this trans-sulfuration reaction is the Cys412 persulfide. Substitution of the EanB-Cys412 persulfide with a Cys412 perselenide does not yield the selenium analog of ergothioneine, selenoneine. However, in deuterated buffer, the perselenide-modified EanB catalyzes the deuterium exchange between hercynine's sp² ε-C-H bond and D₂O. Results from QM/MM calculations suggest that the reaction involves a carbene intermediate and that Tyr353 plays a key role. We hypothesize that modulating the pK_a of Tyr353 will affect the deuterium-exchange rate. Indeed, the 3,5-difluoro tyrosine containing EanB catalyzes the deuterium exchange reaction with k _ex of ~10-fold greater than the wild-type EanB (EanB_WT). With regards to potential mechanisms, these results support the involvement of a carbene intermediate in EanB-catalysis, rendering EanB as one of the few carbene-intermediate involving enzymatic systems.

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.

A Directed Evolution System for Lysine Deacetylases

Lysine acetylation is a ubiquitous modification permeating the proteomes of organisms from all domains of life. Lysine deacetylases (KDACs) reverse this modification by following two fundamentally different enzymatic mechanisms, which differ mainly by the need for NAD⁺ as stoichiometric co-substrate. KDACs are often found as catalytic subunit in protein complexes involved in cell cycle regulation, chromatin organization and transcription. Their promiscuity with respect to sequence context and type of lysine acylation convolutes the network of functional and physical connections.Here we present an efficient selection method for KDACs in E. coli, which allows for the creation of acyl-type specific KDAC variants, which greatly facilitate the investigation of their physiological function . The selection system builds on the incorporation of acylated lysines by genetic code expansion in reporter enzymes with essential lysine residues. We describe the creation of KDAC mutant libraries by saturation mutagenesis of active site residues, the isolation of individual mutants from this library using the selection system, and their biochemical characterization with acylated firefly luciferase.

Design and Applications of Water-Soluble Coordination Cages

Compartmentalization of the aqueous space within a cell is necessary for life. In similar fashion to the nanometer-scale compartments in living systems, synthetic water-soluble coordination cages (WSCCs) can isolate guest molecules and host chemical transformations. Such cages thus show promise in biological, medical, environmental, and industrial domains. This review highlights examples of three-dimensional synthetic WSCCs, offering perspectives so as to enhance their design and applications. Strategies are presented that address key challenges for the preparation of coordination cages that are soluble and stable in water. The peculiarities of guest binding in aqueous media are examined, highlighting amplified binding in water, changing guest properties, and the recognition of specific molecular targets. The properties of WSCC hosts associated with biomedical applications, and their use as vessels to carry out chemical reactions in water, are also presented. These examples sketch a blueprint for the preparation of new metal-organic containers for use in aqueous solution, as well as guidelines for the engineering of new applications in water.

Ground-State Destabilization by Active-Site Hydrophobicity Controls the Selectivity of a Cofactor-Free Decarboxylase

Bacterial arylmalonate decarboxylase (AMDase) and evolved variants have become a valuable tool with which to access both enantiomers of a broad range of chiral arylaliphatic acids with high optical purity. Yet, the molecular principles responsible for the substrate scope, activity, and selectivity of this enzyme are only poorly understood to date, greatly hampering the predictability and design of improved enzyme variants for specific applications. In this work, empirical valence bond and metadynamics simulations were performed on wild-type AMDase and variants thereof to obtain a better understanding of the underlying molecular processes determining reaction outcome. Our results clearly reproduce the experimentally observed substrate scope and support a mechanism driven by ground-state destabilization of the carboxylate group being cleaved by the enzyme. In addition, our results indicate that, in the case of the nonconverted or poorly converted substrates studied in this work, increased solvent exposure of the active site upon binding of these substrates can disturb the vulnerable network of interactions responsible for facilitating the AMDase-catalyzed cleavage of CO2. Finally, our results indicate a switch from preferential cleavage of the pro-(R) to the pro-(S) carboxylate group in the CLG-IPL variant of AMDase for all substrates studied. This appears to be due to the emergence of a new hydrophobic pocket generated by the insertion of the six amino acid substitutions, into which the pro-(S) carboxylate binds. Our results allow insight into the tight interaction network determining AMDase selectivity, which in turn provides guidance for the identification of target residues for future enzyme engineering.

Crystal structure of cis-aconitate decarboxylase reveals the impact of naturally occurring human mutations on itaconate synthesis

cis-Aconitate decarboxylase (CAD, also known as ACOD1 or Irg1) converts cis-aconitate to itaconate and plays central roles in linking innate immunity with metabolism and in the biotechnological production of itaconic acid by Aspergillus terreus We have elucidated the crystal structures of human and murine CADs and compared their enzymological properties to CAD from A. terreus Recombinant CAD is fully active in vitro without a cofactor. Murine CAD has the highest catalytic activity, whereas Aspergillus CAD is best adapted to a more acidic pH. CAD is not homologous to any known decarboxylase and appears to have evolved from prokaryotic enzymes that bind negatively charged substrates. CADs are homodimers, the active center is located in the interface between 2 distinct subdomains, and structural modeling revealed conservation in zebrafish and Aspergillus We identified 8 active-site residues critical for CAD function and rare naturally occurring human mutations in the active site that abolished CAD activity, as well as a variant (Asn152Ser) that increased CAD activity and is common (allele frequency 20%) in African ethnicity. These results open the way for 1) assessing the potential impact of human CAD variants on disease risk at the population level, 2) developing therapeutic interventions to modify CAD activity, and 3) improving CAD efficiency for biotechnological production of itaconic acid.

Biologically generated carbon dioxide: nature's versatile chemical strategies for carboxy lyases

Covering: up to 2019Metabolic production of CO₂ is natural product chemistry on a mammoth scale. Just counting humans, among all other respiring organisms, the seven billion people on the planet exhale about 3 billion tons of CO₂ per year. Essentially all of the biogenic CO₂ arises by action of discrete families of decarboxylases. The mechanistic routes to CO₂ release from carboxylic acid metabolites vary with the electronic demands and structures of specific substrates and illustrate the breadth of chemistry employed for C-COO (C-C bond) disconnections. Most commonly decarboxylated are α-keto acid and β-keto acid substrates, the former requiring thiamin-PP as cofactor, the latter typically cofactor-free. The extensive decarboxylation of amino acids, e.g. to neurotransmitter amines, is synonymous with the coenzyme form of vitamin B₆, pyridoxal-phosphate, although covalent N-terminal pyruvamide residues serve in some amino acid decarboxylases. All told, five B vitamins (B₁, B₂, B₃, B₆, B₇), ATP, S-adenosylmethionine, manganese and zinc ions are pressed into service for specific decarboxylase catalyses. There are additional cofactor-independent decarboxylases that operate by distinct chemical routes. Finally, while most decarboxylases use heterolytic ionic mechanisms, a small number of decarboxylases carry out radical pathways.

Functional classification of protein structures by local structure matching in graph representation

As a result of high‐throughput protein structure initiatives, over 14,400 protein structures have been solved by Structural Genomics (SG) centers and participating research groups. While the totality of SG data represents a tremendous contribution to genomics and structural biology, reliable functional information for these proteins is generally lacking. Better functional predictions for SG proteins will add substantial value to the structural information already obtained. Our method described herein, Graph Representation of Active Sites for Prediction of Function (GRASP‐Func), predicts quickly and accurately the biochemical function of proteins by representing residues at the predicted local active site as graphs rather than in Cartesian coordinates. We compare the GRASP‐Func method to our previously reported method, Structurally Aligned Local Sites of Activity (SALSA), using the Ribulose Phosphate Binding Barrel (RPBB), 6‐Hairpin Glycosidase (6‐HG), and Concanavalin A‐like Lectins/Glucanase (CAL/G) superfamilies as test cases. In each of the superfamilies, SALSA and the much faster method GRASP‐Func yield similar correct classification of previously characterized proteins, providing a validated benchmark for the new method. In addition, we analyzed SG proteins using our SALSA and GRASP‐Func methods to predict function. Forty‐one SG proteins in the RPBB superfamily, nine SG proteins in the 6‐HG superfamily, and one SG protein in the CAL/G superfamily were successfully classified into one of the functional families in their respective superfamily by both methods. This improved, faster, validated computational method can yield more reliable predictions of function that can be used for a wide variety of applications by the community.

A reevaluation of the origin of the rate acceleration for enzyme-catalyzed hydride transfer

There is no consensus of opinion on the origin of the large rate accelerations observed for enzyme-catalyzed hydride transfer. The interpretation of recent results from studies on hydride transfer reactions catalyzed by alcohol dehydrogenase (ADH) focus on the proposal that the effective barrier height is reduced by quantum-mechanical tunneling through the energy barrier. This interpretation contrasts sharply with the notion that enzymatic rate accelerations are obtained through direct stabilization of the transition state for the nonenzymatic reaction in water. The binding energy of the dianion of substrate DHAP provides 11 kcal mol^-1 stabilization of the transition state for the hydride transfer reaction catalyzed by glycerol-3-phosphate dehydrogenase (GPDH). We summarize evidence that the binding interactions between (GPDH) and dianion activators are utilized directly for stabilization of the transition state for enzyme-catalyzed hydride transfer. The possibility is considered, and then discounted, that these dianion binding interactions are utilized for the stabilization of a tunnel ready state (TRS) that enables efficient tunneling of the transferred hydride through the energy barrier, and underneath the energy maximum for the transition state. It is noted that the evidence to support the existence of a tunnel-ready state for the hydride transfer reactions catalyzed by ADH is ambiguous. We propose that the rate acceleration for ADH is due to the utilization of the binding energy of the cofactor NAD+/NADH in the stabilization of the transition state for enzyme-catalyzed hydride transfer.

Isolation of uracil auxotroph mutants of coral symbiont alga for symbiosis studies

Coral reef ecosystems rely on stable symbiotic relationship between the dinoflagellate Symbiodinium spp. and host cnidarian animals. The collapse of such symbiosis could cause coral ‘bleaching’ and subsequent host death. Despite huge interest on Symbiodinium, lack of mutant strains and readily available genetic tools have hampered molecular research. A major issue was the tolerance to marker antibiotics. Here, we isolated Symbiodinium mutants requiring uracil for growth, and hence, useful in transformation screening. We cultured Symbiodinium spp. cells in the presence of 5-fluoroorotic acid (5FOA), which inhibits the growth of cells expressing URA3 encoding orotidine-5′-monophosphate decarboxylase, and isolated cells that require uracil for growth. Sequence analyses and genetic complementation tests using yeast demonstrated that one of the mutant cell lines had a point mutation in URA3, resulting in a splicing error at an unusual exon–intron junction, and consequently, loss of enzyme activity. This mutant could maintain a symbiotic relationship with the model sea anemone Exaiptasia pallida only in sea water containing uracil. Results show that the URA3 mutant will be a useful tool for screening Symbiodinium transformants, both ex and in hospite, as survival in the absence of uracil is possible only upon successful introduction of URA3.

Supramolecular Assembly of Uridine Monophosphate (UMP) and Thymidine Monophosphate (TMP) with a Dinuclear Copper(II) Receptor

Understanding the intermolecular interactions between nucleotides and artificial receptors is crucial to understanding the role of nucleic acids in living systems. However, direct structural evidence showing precise interactions and bonding features of a nucleoside monophosphate (NMP) with a macrocycle-based synthetic molecule has not been provided so far. Herein, we present two novel crystal structures of uridine monophosphate (UMP) and thymidine monophosphate (TMP) complexes with a macrocycle-based dinuclear receptor. Structural characterization of these complexes reveals that the receptor recognizes UMP through coordinate-covalent interactions with phosphates and π-π stackings with nucleobases and TMP through coordinate-covalent interactions with phosphate groups. Furthermore, the receptor has been shown to effectively bind nucleoside monophosphates in the order of GMP > AMP > UMP > TMP > CMP in water at physiological pH, as investigated by an indicator displacement assay.

Primary Deuterium Kinetic Isotope Effects From Product Yields: Rationale, Implementation, and Interpretation

A simple and convenient method is described to determine primary deuterium kinetic isotope effects (1°DKIEs) on reactions where the hydron incorporated into the reaction product is derived from solvent water. The 1°DKIE may be obtained by ¹H NMR analyses as the ratio of the yields of H- and D-labeled products from a reaction in 50:50 (v/v) HOH/DOD. The procedures for these ¹H NMR analyses are reviewed. This product deuterium isotope effect (PDIE) is defined as 1/ϕ_EL for fractionation of hydrons between solvent and the transition state for the reaction examined. When the solvent is not the direct hydron donor, it is necessary to correct the PDIE for the fractionation factor ϕ_EL for partitioning of the hydron between the solvent and the direct donor EL. This method was used to determine the 1°DKIE on decarboxylation reactions catalyzed by wild-type orotidine 5'-monophosphate decarboxylase (OMPDC) and by mutants of OMPDC, and then in the determination of the 1°DKIE on the decarboxylation reaction catalyzed by 5-carboxyvanillate decarboxylase. The experimental procedures used in studies on OMPDC and the rationale for these procedures are described.

Orotidine-Containing RNA: Implications for the Hierarchical Selection (Systems Chemistry Emergence) of RNA

The prebiotic synthesis of canonical nucleobases from HCN is a cornerstone for the RNA world hypothesis. However, their role in the primordial pathways to RNA is still debated. The very same process starting from HCN also gives rise to orotic acid, which (via orotidine) plays a crucial role in extant biology in the de novo synthesis of uridine and cytidine, the informational base-pairs in RNA. However, orotidine itself is absent in RNA. Given the prebiotic and biological relevance of orotic acid vis-à-vis uracil, we investigated orotidine-containing RNA oligonucleotides and show that they have severely compromised base-pairing properties. While not unexpected, these results suggest that the emergence of extant RNA cannot just be a consequence of the plausible prebiotic formation of its chemical constituents/building blocks. In combination with other investigations on alternative prebiotic nucleobases, sugars, and linkers, these findings imply that the selection of the components of extant RNA occurred at a higher hierarchical level of an oligomer/polymer based on its functional properties-pointing to a systems chemistry emergence of RNA from a library of precursors.

Organization of Nucleotides in Different Environments and the Formation of Pre-Polymers

RNA is a linear polymer of nucleotides linked by a ribose-phosphate backbone. Polymerization of nucleotides occurs in a condensation reaction in which phosphodiester bonds are formed. However, in the absence of enzymes and metabolism there has been no obvious way for RNA-like molecules to be produced and then encapsulated in cellular compartments. We investigated 5′-adenosine monophosphate (AMP) and 5′-uridine monophosphate (UMP) molecules confined in multi-lamellar phospholipid bilayers, nanoscopic films, ammonium chloride salt crystals and Montmorillonite clay, previously proposed to promote polymerization. X-ray diffraction was used to determine whether such conditions imposed a degree of order on the nucleotides. Two nucleotide signals were observed in all matrices, one corresponding to a nearest neighbour distance of 4.6 Å attributed to nucleotides that form a disordered, glassy structure. A second, smaller distance of 3.4 Å agrees well with the distance between stacked base pairs in the RNA backbone, and was assigned to the formation of pre-polymers, i.e., the organization of nucleotides into stacks of about 10 monomers. Such ordering can provide conditions that promote the nonenzymatic polymerization of RNA strands under prebiotic conditions. Experiments were modeled by Monte-Carlo simulations, which provide details of the molecular structure of these pre-polymers.

One-carbon chemistry of oxalate oxidoreductase captured by X-ray crystallography

Thiamine pyrophosphate (TPP)-dependent oxalate oxidoreductase (OOR) metabolizes oxalate, generating two molecules of CO2 and two low-potential electrons, thus providing both the carbon and reducing equivalents for operation of the Wood-Ljungdahl pathway of acetogenesis. Here we present structures of OOR in which two different reaction intermediate bound states have been trapped: the covalent adducts between TPP and oxalate and between TPP and CO2. These structures, along with the previously determined structure of substrate-free OOR, allow us to visualize how active site rearrangements can drive catalysis. Our results suggest that OOR operates via a bait-and-switch mechanism, attracting substrate into the active site through the presence of positively charged and polar residues, and then altering the electrostatic environment through loop and side chain movements to drive catalysis. This simple but elegant mechanism explains how oxalate, a molecule that humans and most animals cannot break down, can be used for growth by acetogenic bacteria.

Nucleotides, Nucleosides, and Nucleobases

We review literature on the metabolism of ribo- and deoxyribonucleotides, nucleosides, and nucleobases in Escherichia coli and Salmonella,including biosynthesis, degradation, interconversion, and transport. Emphasis is placed on enzymology and regulation of the pathways, at both the level of gene expression and the control of enzyme activity. The paper begins with an overview of the reactions that form and break the N-glycosyl bond, which binds the nucleobase to the ribosyl moiety in nucleotides and nucleosides, and the enzymes involved in the interconversion of the different phosphorylated states of the nucleotides. Next, the de novo pathways for purine and pyrimidine nucleotide biosynthesis are discussed in detail.Finally, the conversion of nucleosides and nucleobases to nucleotides, i.e.,the salvage reactions, are described. The formation of deoxyribonucleotides is discussed, with emphasis on ribonucleotidereductase and pathways involved in fomation of dUMP. At the end, we discuss transport systems for nucleosides and nucleobases and also pathways for breakdown of the nucleobases.

Cooperative Electrostatic Interactions Drive Functional Evolution in the Alkaline Phosphatase Superfamily

It is becoming widely accepted that catalytic promiscuity, i.e., the ability of a single enzyme to catalyze the turnover of multiple, chemically distinct substrates, plays a key role in the evolution of new enzyme functions. In this context, the members of the alkaline phosphatase superfamily have been extensively studied as model systems in order to understand the phenomenon of enzyme multifunctionality. In the present work, we model the selectivity of two multiply promiscuous members of this superfamily, namely the phosphonate monoester hydrolases from Burkholderia caryophylli and Rhizobium leguminosarum. We have performed extensive simulations of the enzymatic reaction of both wild-type enzymes and several experimentally characterized mutants. Our computational models are in agreement with key experimental observables, such as the observed activities of the wild-type enzymes, qualitative interpretations of experimental pH-rate profiles, and activity trends among several active site mutants. In all cases the substrates of interest bind to the enzyme in similar conformations, with largely unperturbed transition states from their corresponding analogues in aqueous solution. Examination of transition-state geometries and the contribution of individual residues to the calculated activation barriers suggest that the broad promiscuity of these enzymes arises from cooperative electrostatic interactions in the active site, allowing each enzyme to adapt to the electrostatic needs of different substrates. By comparing the structural and electrostatic features of several alkaline phosphatases, we suggest that this phenomenon is a generalized feature driving selectivity and promiscuity within this superfamily and can be in turn used for artificial enzyme design.

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.

Synthesis of orotidine by intramolecular nucleosidation

An intramolecular nucleosidation approach, using orotate as a leaving group, provides easy access to orotidine in high yields.

Thermoadaptation-directed enzyme evolution in an error-prone thermophile derived from Geobacillus kaustophilus HTA426

Thermostability is an important property of enzymes utilized for practical applications because it allows long-term storage and use as catalysts. In this study, we constructed an error-prone strain of the thermophile Geobacillus kaustophilus HTA426 and investigated thermoadaptation-directed enzyme evolution using the strain. A mutation frequency assay using the antibiotics rifampin and streptomycin revealed that G. kaustophilus had substantially higher mutability than Escherichia coli and Bacillus subtilis. The predominant mutations in G. kaustophilus were A · T→G · C and C · G→T · A transitions, implying that the high mutability of G. kaustophilus was attributable in part to high-temperature-associated DNA damage during growth. Among the genes that may be involved in DNA repair in G. kaustophilus, deletions of the mutSL, mutY, ung, and mfd genes markedly enhanced mutability. These genes were subsequently deleted to construct an error-prone thermophile that showed much higher (700- to 9,000-fold) mutability than the parent strain. The error-prone strain was auxotrophic for uracil owing to the fact that the strain was deficient in the intrinsic pyrF gene. Although the strain harboring Bacillus subtilis pyrF was also essentially auxotrophic, cells became prototrophic after 2 days of culture under uracil starvation, generating B. subtilis PyrF variants with an enhanced half-denaturation temperature of >10°C. These data suggest that this error-prone strain is a promising host for thermoadaptation-directed evolution to generate thermostable variants from thermolabile enzymes.

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.

Substrate distortion contributes to the catalysis of orotidine 5'-monophosphate decarboxylase

Orotidine 5'-monophosphate decarboxylase (ODCase) accelerates the decarboxylation of orotidine 5'-monophosphate (OMP) to uridine 5'-monophosphate (UMP) by 17 orders of magnitude. Eight new crystal structures with ligand analogues combined with computational analyses of the enzyme's short-lived intermediates and the intrinsic electronic energies to distort the substrate and other ligands improve our understanding of the still controversially discussed reaction mechanism. In their respective complexes, 6-methyl-UMP displays significant distortion of its methyl substituent bond, 6-amino-UMP shows the competition between the K72 and C6 substituents for a position close to D70, and the methyl and ethyl esters of OMP both induce rotation of the carboxylate group substituent out of the plane of the pyrimidine ring. Molecular dynamics and quantum mechanics/molecular mechanics computations of the enzyme-substrate complex also show the bond between the carboxylate group and the pyrimidine ring to be distorted, with the distortion contributing a 10-15% decrease of the ΔΔG(⧧) value. These results are consistent with ODCase using both substrate distortion and transition-state stabilization, primarily exerted by K72, in its catalysis of the OMP decarboxylation reaction.

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.

Catalysis in Enzymatic Decarboxylations: Comparison of Selected Cofactor-dependent and Cofactor-independent Examples

This review is focused on three types of enzymes decarboxylating very different substrates: (1) Thiamin diphosphate (ThDP)-dependent enzymes reacting with 2-oxo acids; (2) Pyridoxal phosphate (PLP)-dependent enzymes reacting with α-amino acids; and (3) An enzyme with no known co-factors, orotidine 5'-monophosphate decarboxylase (OMPDC). While the first two classes have been much studied for many years, during the past decade studies of both classes have revealed novel mechanistic insight challenging accepted understanding. The enzyme OMPDC has posed a challenge to the enzymologist attempting to explain a 10¹⁷-fold rate acceleration in the absence of cofactors or even metal ions. A comparison of the available evidence on the three types of decarboxylases underlines some common features and more differences. The field of decarboxylases remains an interesting and challenging one for the mechanistic enzymologist notwithstanding the large amount of information already available.

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.

Ter-dependent stress response systems: novel pathways related to metal sensing, production of a nucleoside-like metabolite, and DNA-processing

The mode of action of the bacterial ter cluster and TelA genes, implicated in natural resistance to tellurite and other xenobiotic toxic compounds, pore-forming colicins and several bacteriophages, has remained enigmatic for almost two decades. Using comparative genomics, sequence-profile searches and structural analysis we present evidence that the ter gene products and their functional partners constitute previously underappreciated, chemical stress response and anti-viral defense systems of bacteria. Based on contextual information from conserved gene neighborhoods and domain architectures, we show that the ter gene products and TelA lie at the center of membrane-linked metal recognition complexes with regulatory ramifications encompassing phosphorylation-dependent signal transduction, RNA-dependent regulation, biosynthesis of nucleoside-like metabolites and DNA processing. Our analysis suggests that the multiple metal-binding and non-binding TerD paralogs and TerC are likely to constitute a membrane-associated complex, which might also include TerB and TerY, and feature several, distinct metal-binding sites. Versions of the TerB domain might also bind small molecule ligands and link the TerD paralog-TerC complex to biosynthetic modules comprising phosphoribosyltransferases (PRTases), ATP grasp amidoligases, TIM-barrel carbon-carbon lyases, and HAD phosphoesterases, which are predicted to synthesize novel nucleoside-like molecules. One of the PRTases is also likely to interact with RNA by means of its Pelota/Ribosomal protein L7AE-like domain. The von Willebrand factor A domain protein, TerY, is predicted to be part of a distinct phosphorylation switch, coupling a protein kinase and a PP2C phosphatase. We show, based on the evidence from numerous conserved gene neighborhoods and domain architectures, that both the TerB and TelA domains have been linked to diverse lipid-interaction domains, such as two novel PH-like and the Coq4 domains, in different bacteria, and are likely to comprise membrane-associated sensory complexes that might additionally contain periplasmic binding-protein-II and OmpA domains. We also show that the TerD and TerB domains and the TerY-associated phosphorylation system are functionally linked to many distinct DNA-processing complexes, which feature proteins with SWI2/SNF2 and RecQ-like helicases, multiple AAA+ ATPases, McrC-N-terminal domain proteins, several restriction endonuclease fold DNases, DNA-binding domains and a type-VII/Esx-like system, which is at the center of a predicted DNA transfer apparatus. These DNA-processing modules and associated genes are predicted to be involved in restriction or suicidal action in response to phages and possibly repairing xenobiotic-induced DNA damage. In some eukaryotes, certain components of the ter system appear to be recruited to function in conjunction with the ubiquitin system and calcium-signaling pathways.

Specificity in transition state binding: the Pauling model revisited

Linus Pauling proposed that the large rate accelerations for enzymes are caused by the high specificity of the protein catalyst for binding the reaction transition state. The observation that stable analogues of the transition states for enzymatic reactions often act as tight-binding inhibitors provided early support for this simple and elegant proposal. We review experimental results that support the proposal that Pauling's model provides a satisfactory explanation for the rate accelerations for many heterolytic enzymatic reactions through high-energy reaction intermediates, such as proton transfer and decarboxylation. Specificity in transition state binding is obtained when the total intrinsic binding energy of the substrate is significantly larger than the binding energy observed at the Michaelis complex. The results of recent studies that aimed to characterize the specificity in binding of the enolate oxygen at the transition state for the 1,3-isomerization reaction catalyzed by ketosteroid isomerase are reviewed. Interactions between pig heart succinyl-coenzyme A:3-oxoacid coenzyme A transferase (SCOT) and the nonreacting portions of coenzyme A (CoA) are responsible for a rate increase of 3 × 10(12)-fold, which is close to the estimated total 5 × 10(13)-fold enzymatic rate acceleration. Studies that partition the interactions between SCOT and CoA into their contributing parts are reviewed. Interactions of the protein with the substrate phosphodianion group provide an ~12 kcal/mol stabilization of the transition state for the reactions catalyzed by triosephosphate isomerase, orotidine 5'-monophosphate decarboxylase, and α-glycerol phosphate dehydrogenase. The interactions of these enzymes with the substrate piece phosphite dianion provide a 6-8 kcal/mol stabilization of the transition state for reaction of the appropriate truncated substrate. Enzyme activation by phosphite dianion reflects the higher dianion affinity for binding to the enzyme-transition state complex compared with that of the free enzyme. Evidence is presented that supports a model in which the binding energy of the phosphite dianion piece, or the phosphodianion group of the whole substrate, is utilized to drive an enzyme conformational change from an inactive open form E(O) to an active closed form E(C), by closure of a phosphodianion gripper loop. Members of the enolase and haloalkanoic acid dehalogenase superfamilies use variable capping domains to interact with nonreacting portions of the substrate and sequester the substrate from interaction with bulk solvent. Interactions of this capping domain with the phenyl group of mandelate have been shown to activate mandelate racemase for catalysis of deprotonation of α-carbonyl carbon. We propose that an important function of these capping domains is to utilize the binding interactions with nonreacting portions of the substrate to activate the enzyme for catalysis.

Protein function annotation with Structurally Aligned Local Sites of Activity (SALSAs)

Background

The prediction of biochemical function from the 3D structure of a protein has proved to be much more difficult than was originally foreseen. A reliable method to test the likelihood of putative annotations and to predict function from structure would add tremendous value to structural genomics data. We report on a new method, Structurally Aligned Local Sites of Activity (SALSA), for the prediction of biochemical function based on a local structural match at the predicted catalytic or binding site.

Results

Implementation of the SALSA method is described. For the structural genomics protein PY01515 (PDB ID 2aqw) from Plasmodium yoelii, it is shown that the putative annotation, Orotidine 5'-monophosphate decarboxylase (OMPDC), is most likely correct. SALSA analysis of YP_001304206.1 (PDB ID 3h3l), a putative sugar hydrolase from Parabacteroides distasonis, shows that its active site does not bear close resemblance to any previously characterized member of its superfamily, the Concanavalin A-like lectins/glucanases. It is noted that three residues in the active site of the thermophilic beta-1,4-xylanase from Nonomuraea flexuosa (PDB ID 1m4w), Y78, E87, and E176, overlap with POOL-predicted residues of similar type, Y168, D153, and E232, in YP_001304206.1. The substrate recognition regions of the two proteins are rather different, suggesting that YP_001304206.1 is a new functional type within the superfamily. A structural genomics protein from Mycobacterium avium (PDB ID 3q1t) has been reported to be an enoyl-CoA hydratase (ECH), but SALSA analysis shows a poor match between the predicted residues for the SG protein and those of known ECHs. A better local structural match is obtained with Anabaena beta-diketone hydrolase (ABDH), a known β-diketone hydrolase from Cyanobacterium anabaena (PDB ID 2j5s). This suggests that the reported ECH function of the SG protein is incorrect and that it is more likely a β-diketone hydrolase.

Conclusions

A local site match provides a more compelling function prediction than that obtainable from a simple 3D structure match. The present method can confirm putative annotations, identify misannotation, and in some cases suggest a more probable annotation.

Atomic resolution structure of the orotidine 5'-monophosphate decarboxylase product complex combined with surface plasmon resonance analysis: implications for the catalytic mechanism

Orotidine 5'-monophosphate decarboxylase (ODCase) accelerates the decarboxylation of its substrate by 17 orders of magnitude. One argument brought forward against steric/electrostatic repulsion causing substrate distortion at the carboxylate substituent as part of the catalysis has been the weak binding affinity of the decarboxylated product (UMP). The crystal structure of the UMP complex of ODCase at atomic resolution (1.03 Å) shows steric competition between the product UMP and the side chain of a catalytic lysine residue. Surface plasmon resonance analysis indicates that UMP binds 5 orders of magnitude more tightly to a mutant in which the interfering side chain has been removed than to wild-type ODCase. These results explain the low affinity of UMP and counter a seemingly very strong argument against a contribution of substrate distortion to the catalytic reaction mechanism of ODCase.

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.

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.

Decarboxylation mechanisms in biological system

This review examines the mechanisms propelling cofactor-independent, organic cofactor-dependent and metal-dependent decarboxylase chemistry. Decarboxylation, the removal of carbon dioxide from organic acids, is a fundamentally important reaction in biology. Numerous decarboxylase enzymes serve as key components of aerobic and anaerobic carbohydrate metabolism and amino acid conversion. In the past decade, our knowledge of the mechanisms enabling these crucial decarboxylase reactions has continued to expand and inspire. This review focuses on the organic cofactors biotin, flavin, NAD, pyridoxal 5'-phosphate, pyruvoyl, and thiamin pyrophosphate as catalytic centers. Significant attention is also placed on the metal-dependent decarboxylase mechanisms.

Emergent decarboxylase activity and attenuation of α/β-hydrolase activity during the evolution of methylketone biosynthesis in tomato

Specialized methylketone-containing metabolites accumulate in certain plants, in particular wild tomatoes in which they serve as toxic compounds against chewing insects. In Solanum habrochaites f. glabratum, methylketone biosynthesis occurs in the plastids of glandular trichomes and begins with intermediates of de novo fatty acid synthesis. These fatty-acyl intermediates are converted via sequential reactions catalyzed by Methylketone Synthase2 (MKS2) and MKS1 to produce the n-1 methylketone. We report crystal structures of S. habrochaites MKS1, an atypical member of the α/β-hydrolase superfamily. Sequence comparisons revealed the MKS1 catalytic triad, Ala-His-Asn, as divergent to the traditional α/β-hydrolase triad, Ser-His-Asp. Determination of the MKS1 structure points to a novel enzymatic mechanism dependent upon residues Thr-18 and His-243, confirmed by biochemical assays. Structural analysis further reveals a tunnel leading from the active site consisting mostly of hydrophobic residues, an environment well suited for fatty-acyl chain binding. We confirmed the importance of this substrate binding mode by substituting several amino acids leading to an alteration in the acyl-chain length preference of MKS1. Furthermore, we employ structure-guided mutagenesis and functional assays to demonstrate that MKS1, unlike enzymes from this hydrolase superfamily, is not an efficient hydrolase but instead catalyzes the decarboxylation of 3-keto acids.

Mechanism of the orotidine 5'-monophosphate decarboxylase-catalyzed reaction: importance of residues in the orotate binding site

The reaction catalyzed by orotidine 5'-monophosphate decarboxylase (OMPDC) is accompanied by exceptional values for rate enhancement (k(cat)/k(non) = 7.1 × 10(16)) and catalytic proficiency [(k(cat)/K(M))/k(non) = 4.8 × 10(22) M(-1)]. Although a stabilized vinyl carbanion/carbene intermediate is located on the reaction coordinate, the structural strategies by which the reduction in the activation energy barrier is realized remain incompletely understood. This laboratory recently reported that "substrate destabilization" by Asp 70 in the OMPDC from Methanothermobacter thermoautotrophicus (MtOMPDC) lowers the activation energy barrier by ∼5 kcal/mol (contributing ~2.7 × 10(3) to the rate enhancement) [Chan, K. K., Wood, B. M., Fedorov, A. A., Fedorov, E. V., Imker, H. J., Amyes, T. L., Richard, J. P., Almo, S. C., and Gerlt, J. A. (2009) Biochemistry 48, 5518-5531]. We now report that substitutions of hydrophobic residues in a pocket proximal to the carboxylate group of the substrate (Ile 96, Leu 123, and Val 155) with neutral hydrophilic residues decrease the value of k(cat) by as much as 400-fold but have a minimal effect on the value of k(ex) for exchange of H6 of the FUMP product analogue with solvent deuterium; we hypothesize that this pocket destabilizes the substrate by preventing hydration of the substrate carboxylate group. We also report that substitutions of Ser 127 that is proximal to O4 of the orotate ring decrease the value of k(cat)/K(M), with the S127P substitution that eliminates hydrogen bonding interactions with O4 producing a 2.5 × 10(6)-fold reduction; this effect is consistent with delocalization of the negative charge of the carbanionic intermediate on O4 that produces an anionic carbene intermediate and thereby provides a structural strategy for stabilization of the intermediate. These observations provide additional information about the identities of the active site residues that contribute to the rate enhancement and, therefore, insights into the structural strategies for catalysis.

Determination of the amino acid sequence requirements for catalysis by the highly proficient orotidine monophosphate decarboxylase

Orotidine 5'-monophosphate decarboxylase (ODCase) catalyzes the decarboxylation of orotidine 5'-monophosphate to uridine 5'-monophosphate during pyrimidine nucleotide biosynthesis. This enzyme is one of the most proficient known, exhibiting a rate enhancement of over 17 orders of magnitude over the uncatalyzed rate. An interesting question is whether the high proficiency of ODCase is associated with a highly optimized sequence of active site residues. This question was addressed by randomizing 24 residue positions in and around the active site of the E. coli ODCase (pyrF) by site-directed mutagenesis. The libraries of mutants were selected for function from a multicopy plasmid or by single-copy replacement at the pyrF locus on the E. coli chromosome. Stringent sequence requirements for function were found for the mutants expressed from the chromosomal pyrF locus. Six positions were not tolerant of substitutions and several others accepted very limited substitutions. In contrast, all positions could be substituted to some extent when the library mutants were expressed from a multicopy plasmid. For the conserved quartet of charged residues Lys44-Asp71-Lys73-Asp76, a cysteine substitution was found to provide function at positions 71 and 76. A lower pK(a) for both cysteine mutants supports a mechanism whereby the thiolate group of cysteine substitutes for the negatively charged aspartate side chain. The partial function mutants such as D71C and D76C exhibit reduced catalytic efficiency relative to wild type but nevertheless provide a rate enhancement of 15 orders of magnitude over the uncatalyzed rate indicating the catalytic proficiency of the enzyme is robust and tolerant of mutation.

Polarization in the structures of uracil and thiouracils: Implication for binding with orotidine 5'-monophosphate decarboxylase

The structures of the uracil and thiouracils were examined using NMR spectroscopy and crystal structure data when available. The relationships between the extent of polarization and the C5-C6 bond length as well as the H5-H6 coupling constants were probed. It was found that the bond length and coupling constants correlate well with the proton affinities at the carbonyl or thiocarbonyl groups at C4 but not C2. The possible implication in the tighter binding of thiouracil based nucleotides to orotidine-5'-monophosphate decarboxylase was discussed.

Accelerated hydrolysis of α-halo and α-cyano pyridinium relative to uracil derivatives: a model for ODCase-catalyzed hydrolysis of 6-cyanoUMP

α-Halo and α-cyano pyridiniums were found to undergo facile hydrolysis, in contrast to the sluggish reactions of corresponding uracils. The greatly enhanced rates found with pyridinium compounds have indicated a possible source of the rate acceleration seen in the hydrolysis of 6-cyanouridine 5'-monophosphate catalyzed by orotidine 5'-monophosphate decarboxylase.