Discovery of novel bacterial queuine salvage enzymes and pathways in human pathogens.

Queuosine (Q) is a complex tRNA modification widespread in eukaryotes and bacteria that contributes to the efficiency and accuracy of protein synthesis. Eukaryotes are not capable of Q synthesis and rely on salvage of the queuine base (q) as a Q precursor. While many bacteria are capable of Q de novo synthesis, salvage of the prokaryotic Q precursors preQ0 and preQ1 also occurs. With the exception of Escherichia coli YhhQ, shown to transport preQ0 and preQ1, the enzymes and transporters involved in Q salvage and recycling have not been well described. We discovered and characterized 2 Q salvage pathways present in many pathogenic and commensal bacteria. The first, found in the intracellular pathogen Chlamydia trachomatis, uses YhhQ and tRNA guanine transglycosylase (TGT) homologs that have changed substrate specificities to directly salvage q, mimicking the eukaryotic pathway. The second, found in bacteria from the gut flora such as Clostridioides difficile, salvages preQ1 from q through an unprecedented reaction catalyzed by a newly defined subgroup of the radical-SAM enzyme family. The source of q can be external through transport by members of the energy-coupling factor (ECF) family or internal through hydrolysis of Q by a dedicated nucleosidase. This work reinforces the concept that hosts and members of their associated microbiota compete for the salvage of Q precursors micronutrients.

Evolution of Enzyme Function and the Development of Catalytic Efficiency: Triosephosphate Isomerase, Jeremy R. Knowles, and W. John Albery.

Every reader knows that an enzyme accelerates a reaction by reducing the activation-energy barrier. However, understanding how this is achieved by the structure of the enzyme and its interactions with stable complexes and transition states and, then, using this to (re)design enzymes to catalyze novel reactions remain the "holy grail" of mechanistic enzymology. The necessary foundation is the free-energy profile that specifies the energies of the bound substate, product, and intervening intermediates as well as the transition states by which they are interconverted. When this free-energy profile is compared to that for the uncatalyzed reaction, strategies for establishing and enhancing catalysis can be identified. This Perspective reminds readers that the first free-energy profile determined for an enzyme-catalyzed reaction, that for triosephosphate isomerase, was published in Biochemistry in 1976 by Jeremy R. Knowles, W. John Albery, and co-workers. They used the profile to propose three steps of increasing "subtlety" that can be influenced by evolutionary pressure to increase the flux through the reaction coordinate: (1) "uniform binding" of the substrate, product, and intermediates; (2) "differential binding" of complexes so that these are isoenergetic (to minimize the energy of the intervening transition states); and (3) "catalysis of an elementary step" in which the transition state for the kinetically significant chemical step is stabilized so that flux can be determined by the rate of substrate binding or product dissociation. These papers continue to guide mechanistic studies of enzyme-catalyzed reactions and provide principles for the (re)design of novel enzymes.

Epoxyqueuosine Reductase QueH in the Biosynthetic Pathway to tRNA Queuosine Is a Unique Metalloenzyme.

Queuosine is a structurally unique and functionally important tRNA modification, widely distributed in eukaryotes and bacteria. The final step of queuosine biosynthesis is the reduction/deoxygenation of epoxyqueuosine to form the cyclopentene motif of the nucleobase. The chemistry is performed by the structurally and functionally characterized cobalamin-dependent QueG. However, the queG gene is absent from several bacteria that otherwise retain queuosine biosynthesis machinery. Members of the IPR003828 family (previously known as DUF208) have been recently identified as nonorthologous replacements of QueG, and this family was renamed QueH. Here, we present the structural characterization of QueH from Thermotoga maritima. The structure reveals an unusual active site architecture with a [4Fe-4S] metallocluster along with an adjacent coordinated iron metal. The juxtaposition of the cofactor and coordinated metal ion predicts a unique mechanism for a two-electron reduction/deoxygenation of epoxyqueuosine. To support the structural characterization, in vitro biochemical and genomic analyses are presented. Overall, this work reveals new diversity in the chemistry of iron/sulfur-dependent enzymes and novel insight into the last step of this widely conserved tRNA modification.

A bifunctional salvage pathway for two distinct S-adenosylmethionine by-products that is widespread in bacteria, including pathogenic Escherichia coli.

S-adenosyl-l-methionine (SAM) is a necessary cosubstrate for numerous essential enzymatic reactions including protein and nucleotide methylations, secondary metabolite synthesis and radical-mediated processes. Radical SAM enzymes produce 5'-deoxyadenosine, and SAM-dependent enzymes for polyamine, neurotransmitter and quorum sensing compound synthesis produce 5'-methylthioadenosine as by-products. Both are inhibitory and must be addressed by all cells. This work establishes a bifunctional oxygen-independent salvage pathway for 5'-deoxyadenosine and 5'-methylthioadenosine in both Rhodospirillum rubrum and Extraintestinal Pathogenic Escherichia coli. Homologous genes for this pathway are widespread in bacteria, notably pathogenic strains within several families. A phosphorylase (Rhodospirillum rubrum) or separate nucleoside and kinase (Escherichia coli) followed by an isomerase and aldolase sequentially function to salvage these two wasteful and inhibitory compounds into adenine, dihydroxyacetone phosphate and acetaldehyde or (2-methylthio)acetaldehyde during both aerobic and anaerobic growth. Both SAM by-products are metabolized with equal affinity during aerobic and anaerobic growth conditions, suggesting that the dual-purpose salvage pathway plays a central role in numerous environments, notably the human body during infection. Our newly discovered bifunctional oxygen-independent pathway, widespread in bacteria, salvages at least two by-products of SAM-dependent enzymes for carbon and sulfur salvage, contributing to cell growth.

Characterization of an l-Ascorbate Catabolic Pathway with Unprecedented Enzymatic Transformations.

l-Ascorbate (vitamin C) is ubiquitous in both our diet and the environment. Here we report that Ralstonia eutropha H16 (Cupriavidus necator ATCC 17699) uses l-ascorbate as sole carbon source via a novel catabolic pathway. RNaseq identified eight candidate catabolic genes, sequence similarity networks, and genome neighborhood networks guided predictions for function of the encoded proteins, and the predictions were confirmed by in vitro assays and in vivo growth phenotypes of gene deletion mutants. l-Ascorbate, a lactone, is oxidized and ring-opened by enzymes in the cytochrome b561 and gluconolactonase families, respectively, to form 2,3-diketo-l-gulonate. A protein predicted to have a WD40-like fold catalyzes an unprecedented benzilic acid rearrangement involving migration of a carboxylate group to form 2-carboxy-l-lyxonolactone; the lactone is hydrolyzed by a member of the amidohydrolase superfamily to yield 2-carboxy-l-lyxonate. A member of the PdxA family of oxidative decarboxylases catalyzes a novel decarboxylation that uses NAD+ catalytically. The product, l-lyxonate, is catabolized to α-ketoglutarate by a previously characterized pathway. The pathway is found in hundreds of bacteria, including the pathogens Pseudomonas aeruginosa and Acinetobacter baumannii.

Functional assignment of multiple catabolic pathways for D-apiose.

Colocation of the genes encoding ABC, TRAP, and TCT transport systems and catabolic pathways for the transported ligand provides a strategy for discovering novel microbial enzymes and pathways. We screened solute-binding proteins (SBPs) for ABC transport systems and identified three that bind D-apiose, a branched pentose in the cell walls of higher plants. Guided by sequence similarity networks (SSNs) and genome neighborhood networks (GNNs), the identities of the SBPs enabled the discovery of four catabolic pathways for D-apiose with eleven previously unknown reactions. The new enzymes include D-apionate oxidoisomerase, which catalyzes hydroxymethyl group migration, as well as 3-oxo-isoapionate-4-phosphate decarboxylase and 3-oxo-isoapionate-4-phosphate transcarboxylase/hydrolase, which are RuBisCO-like proteins (RLPs). The web tools for generating SSNs and GNNs are publicly accessible ( http://efi.igb.illinois.edu/efi-est/ ), so similar 'genomic enzymology' strategies for discovering novel pathways can be used by the community.

Novel Metabolic Pathways and Regulons for Hexuronate Utilization in Proteobacteria.

We used comparative genomics to reconstruct d-galacturonic and d-glucuronic acid catabolic pathways and associated transcriptional regulons involving the tripartite ATP-independent periplasmic (TRAP) family transporters that bind hexuronates in proteobacteria. The reconstructed catabolic network involves novel transcription factors, catabolic enzymes, and transporters for utilization of both hexuronates and aldarates (d-glucarate and meso-galactarate). The reconstructed regulons for a novel GntR family transcription factor, GguR, include the majority of hexuronate/aldarate utilization genes in 47 species from the Burkholderiaceae, Comamonadaceae, Halomonadaceae, and Pseudomonadaceae families. GudR, GulR, and UdhR are additional local regulators of some hexuronate/aldarate utilization genes in some of the above-mentioned organisms. The predicted DNA binding motifs of GguR and GudR regulators from Ralstonia pickettii and Polaromonas were validated by in vitro binding assays. Genes from the GulR- and GguR-controlled loci were differentially expressed in R. pickettii grown on hexuronates and aldarates. By a combination of bioinformatics and experimental techniques we identified a novel variant of the oxidative pathway for hexuronate utilization, including two previously uncharacterized subfamilies of lactone hydrolases (UxuL and UxuF). The genomic context of respective genes and reconstruction of associated pathways suggest that both enzymes catalyze the conversion of d-galactaro- and d-glucaro-1,5-lactones to the ring-opened aldarates. The activities of the purified recombinant enzymes, UxuL and UxuF, from four proteobacterial species were directly confirmed and kinetically characterized. The inferred novel aldarate-specific transporter from the tripartite tricarboxylate transporter (TTT) family transporter TctC was confirmed to bind d-glucarate in vitro This study expands our knowledge of bacterial carbohydrate catabolic pathways by identifying novel families of catabolic enzymes, transcriptional regulators, and transporters.IMPORTANCE Hexuronate catabolic pathways and their transcriptional networks are highly variable among different bacteria. We identified novel transcriptional regulators that control the hexuronate and aldarate utilization genes in four families of proteobacteria. By regulon reconstruction and genome context analysis we identified several novel components of the common hexuronate/aldarate utilization pathways, including novel uptake transporters and catabolic enzymes. Two novel families of lactonases involved in the oxidative pathway of hexuronate catabolism were characterized. Novel transcriptional regulons were validated via in vitro binding assays and gene expression studies with Polaromonas and Ralstonia species. The reconstructed catabolic pathways are interconnected with each other metabolically and coregulated via the GguR regulons in proteobacteria.

The EFI Web Resource for Genomic Enzymology Tools: Leveraging Protein, Genome, and Metagenome Databases to Discover Novel Enzymes and Metabolic Pathways.

The assignment of functions to uncharacterized proteins discovered in genome projects requires easily accessible tools and computational resources for large-scale, user-friendly leveraging of the protein, genome, and metagenome databases by experimentalists. This article describes the web resource developed by the Enzyme Function Initiative (EFI; accessed at https://efi.igb.illinois.edu/ ) that provides "genomic enzymology" tools ("web tools") for (1) generating sequence similarity networks (SSNs) for protein families (EFI-EST); (2) analyzing and visualizing genome context of the proteins in clusters in SSNs (in genome neighborhood networks, GNNs, and genome neighborhood diagrams, GNDs) (EFI-GNT); and (3) prioritizing uncharacterized SSN clusters for functional assignment based on metagenome abundance (chemically guided functional profiling, CGFP) (EFI-CGFP). The SSNs generated by EFI-EST are used as the input for EFI-GNT and EFI-CGFP, enabling easy transfer of information among the tools. The networks are visualized and analyzed using Cytoscape, a widely used desktop application; GNDs and CGFP heatmaps summarizing metagenome abundance are viewed within the tools. We provide a detailed example of the integrated use of the tools with an analysis of glycyl radical enzyme superfamily (IPR004184) found in the human gut microbiome. This analysis demonstrates that (1) SwissProt annotations are not always correct, (2) large-scale genome context analyses allow the prediction of novel metabolic pathways, and (3) metagenome abundance can be used to identify/prioritize uncharacterized proteins for functional investigation.

Discovery of new enzymes and metabolic pathways by using structure and genome context.

Assigning valid functions to proteins identified in genome projects is challenging: overprediction and database annotation errors are the principal concerns. We and others are developing computation-guided strategies for functional discovery with 'metabolite docking' to experimentally derived or homology-based three-dimensional structures. Bacterial metabolic pathways often are encoded by 'genome neighbourhoods' (gene clusters and/or operons), which can provide important clues for functional assignment. We recently demonstrated the synergy of docking and pathway context by 'predicting' the intermediates in the glycolytic pathway in Escherichia coli. Metabolite docking to multiple binding proteins and enzymes in the same pathway increases the reliability of in silico predictions of substrate specificities because the pathway intermediates are structurally similar. Here we report that structure-guided approaches for predicting the substrate specificities of several enzymes encoded by a bacterial gene cluster allowed the correct prediction of the in vitro activity of a structurally characterized enzyme of unknown function (PDB 2PMQ), 2-epimerization of trans-4-hydroxy-L-proline betaine (tHyp-B) and cis-4-hydroxy-D-proline betaine (cHyp-B), and also the correct identification of the catabolic pathway in which Hyp-B 2-epimerase participates. The substrate-liganded pose predicted by virtual library screening (docking) was confirmed experimentally. The enzymatic activities in the predicted pathway were confirmed by in vitro assays and genetic analyses; the intermediates were identified by metabolomics; and repression of the genes encoding the pathway by high salt concentrations was established by transcriptomics, confirming the osmolyte role of tHyp-B. This study establishes the utility of structure-guided functional predictions to enable the discovery of new metabolic pathways.

Assignment of function to a domain of unknown function: DUF1537 is a new kinase family in catabolic pathways for acid sugars.

Using a large-scale "genomic enzymology" approach, we (i) assigned novel ATP-dependent four-carbon acid sugar kinase functions to members of the DUF1537 protein family (domain of unknown function; Pfam families PF07005 and PF17042) and (ii) discovered novel catabolic pathways for d-threonate, l-threonate, and d-erythronate. The experimentally determined ligand specificities of several solute binding proteins (SBPs) for TRAP (tripartite ATP-independent permease) transporters for four-carbon acids, including d-erythronate and l-erythronate, were used to constrain the substrates for the catabolic pathways that degrade the SBP ligands to intermediates in central carbon metabolism. Sequence similarity networks and genome neighborhood networks were used to identify the enzyme components of the pathways. Conserved genome neighborhoods encoded SBPs as well as permease components of the TRAP transporters, members of the DUF1537 family, and a member of the 4-hydroxy-l-threonine 4-phosphate dehydrogenase (PdxA) oxidative decarboxylase, class II aldolase, or ribulose 1,5-bisphosphate carboxylase/oxygenase, large subunit (RuBisCO) superfamily. Because the characterized substrates of members of the PdxA, class II aldolase, and RuBisCO superfamilies are phosphorylated, we postulated that the members of the DUF1537 family are novel ATP-dependent kinases that participate in catabolic pathways for four-carbon acid sugars. We determined that (i) the DUF1537/PdxA pair participates in a pathway for the conversion of d-threonate to dihydroxyacetone phosphate and CO2 and (ii) the DUF1537/class II aldolase pair participates in pathways for the conversion of d-erythronate and l-threonate (epimers at carbon-3) to dihydroxyacetone phosphate and CO2 The physiological importance of these pathways was demonstrated in vivo by phenotypic and genetic analyses.

Enzyme Architecture: Breaking Down the Catalytic Cage that Activates Orotidine 5'-Monophosphate Decarboxylase for Catalysis.

We report the results of a study of the catalytic role of a network of four interacting amino acid side chains at yeast orotidine 5'-monophosphate decarboxylase ( ScOMPDC), by the stepwise replacement of all four side chains. The H-bond, which links the -CH2OH side chain of S154 from the pyrimidine umbrella loop of ScOMPDC to the amide side chain of Q215 in the phosphodianion gripper loop, creates a protein cage for the substrate OMP. The role of this interaction in optimizing transition state stabilization from the dianion gripper side chains Q215, Y217, and R235 was probed by determining the kinetic parameter kcat/ Km for 16 enzyme variants, which include all combinations of single, double, triple, and quadruple S154A, Q215A, Y217F, and R235A mutations. The effects of consecutive Q215A, Y217F, and R235A mutations on Δ G⧧ for wild-type enzyme-catalyzed decarboxylation sum to 11.6 kcal/mol, but to only 7.6 kcal/mol when starting from S154A mutant. This shows that the S154A mutation results in a (11.6-7.6) = 4.0 kcal/mol decrease in transition state stabilization from interactions with Q215, Y217, and R235. Mutant cycles show that ca. 2 kcal/mol of this 4 kcal/mol effect is from the direct interaction between the S154 and Q215 side chains and that ca. 2 kcal/mol is from a tightening in the stabilizing interactions of the Y217 and R235 side chains. The sum of the effects of individual A154S, A215Q, F217Y and A235R substitutions at the quadruple mutant of ScOMPDC to give the corresponding triple mutants, 5.6 kcal/mol, is much smaller than 16.0 kcal/mol, the sum of the effects of the related four substitutions in wild-type ScOMPDC to give the respective single mutants. The small effect of substitutions at the quadruple mutant is consistent with a large entropic cost to holding the flexible loops of ScOMPDC in the active closed conformation.

Genomic Enzymology: Web Tools for Leveraging Protein Family Sequence-Function Space and Genome Context to Discover Novel Functions.

The exponentially increasing number of protein and nucleic acid sequences provides opportunities to discover novel enzymes, metabolic pathways, and metabolites/natural products, thereby adding to our knowledge of biochemistry and biology. The challenge has evolved from generating sequence information to mining the databases to integrating and leveraging the available information, i.e., the availability of "genomic enzymology" web tools. Web tools that allow identification of biosynthetic gene clusters are widely used by the natural products/synthetic biology community, thereby facilitating the discovery of novel natural products and the enzymes responsible for their biosynthesis. However, many novel enzymes with interesting mechanisms participate in uncharacterized small-molecule metabolic pathways; their discovery and functional characterization also can be accomplished by leveraging information in protein and nucleic acid databases. This Perspective focuses on two genomic enzymology web tools that assist the discovery novel metabolic pathways: (1) Enzyme Function Initiative-Enzyme Similarity Tool (EFI-EST) for generating sequence similarity networks to visualize and analyze sequence-function space in protein families and (2) Enzyme Function Initiative-Genome Neighborhood Tool (EFI-GNT) for generating genome neighborhood networks to visualize and analyze the genome context in microbial and fungal genomes. Both tools have been adapted to other applications to facilitate target selection for enzyme discovery and functional characterization. As the natural products community has demonstrated, the enzymology community needs to embrace the essential role of web tools that allow the protein and genome sequence databases to be leveraged for novel insights into enzymological problems.

Investigating the role of a backbone to substrate hydrogen bond in OMP decarboxylase using a site-specific amide to ester substitution.

Hydrogen bonds between backbone amide groups of enzymes and their substrates are often observed, but their importance in substrate binding and/or catalysis is not easy to investigate experimentally. We describe the generation and kinetic characterization of a backbone amide to ester substitution in the orotidine 5'-monophosphate (OMP) decarboxylase from Methanobacter thermoautotrophicum (MtOMPDC) to determine the importance of a backbone amide-substrate hydrogen bond. The MtOMPDC-catalyzed reaction is characterized by a rate enhancement (∼10(17)) that is among the largest for enzyme-catalyzed reactions. The reaction proceeds through a vinyl anion intermediate that may be stabilized by hydrogen bonding interaction between the backbone amide of a conserved active site serine residue (Ser-127) and oxygen (O4) of the pyrimidine moiety and/or electrostatic interactions with the conserved general acidic lysine (Lys-72). In vitro translation in conjunction with amber suppression using an orthogonal amber tRNA charged with L-glycerate ((HO)S) was used to generate the ester backbone substitution (S127(HO)S). With 5-fluoro OMP (FOMP) as substrate, the amide to ester substitution increased the value of Km by ∼1.5-fold and decreased the value of kcat by ∼50-fold. We conclude that (i) the hydrogen bond between the backbone amide of Ser-127 and O4 of the pyrimidine moiety contributes a modest factor (∼10(2)) to the 10(17) rate enhancement and (ii) the stabilization of the anionic intermediate is accomplished by electrostatic interactions, including its proximity of Lys-72. These conclusions are in good agreement with predictions obtained from hybrid quantum mechanical/molecular mechanical calculations.

ATP-binding Cassette (ABC) Transport System Solute-binding Protein-guided Identification of Novel d-Altritol and Galactitol Catabolic Pathways in Agrobacterium tumefaciens C58.

Innovations in the discovery of the functions of uncharacterized proteins/enzymes have become increasingly important as advances in sequencing technology flood protein databases with an exponentially growing number of open reading frames. This study documents one such innovation developed by the Enzyme Function Initiative (EFI; U54GM093342), the use of solute-binding proteins for transport systems to identify novel metabolic pathways. In a previous study, this strategy was applied to the tripartite ATP-independent periplasmic transporters. Here, we apply this strategy to the ATP-binding cassette transporters and report the discovery of novel catabolic pathways for d-altritol and galactitol in Agrobacterium tumefaciens C58. These efforts resulted in the description of three novel enzymatic reactions as follows: 1) oxidation of d-altritol to d-tagatose via a dehydrogenase in Pfam family PF00107, a previously unknown reaction; 2) phosphorylation of d-tagatose to d-tagatose 6-phosphate via a kinase in Pfam family PF00294, a previously orphan EC number; and 3) epimerization of d-tagatose 6-phosphate C-4 to d-fructose 6-phosphate via a member of Pfam family PF08013, another previously unknown reaction. The epimerization reaction catalyzed by a member of PF08013 is especially noteworthy, because the functions of members of PF08013 have been unknown. These discoveries were assisted by the following two synergistic bioinformatics web tools made available by the Enzyme Function Initiative: the EFI-Enzyme Similarity Tool and the EFI-Genome Neighborhood Tool.

Enzyme Function Initiative-Enzyme Similarity Tool (EFI-EST): A web tool for generating protein sequence similarity networks.

The Enzyme Function Initiative, an NIH/NIGMS-supported Large-Scale Collaborative Project (EFI; U54GM093342; http://enzymefunction.org/), is focused on devising and disseminating bioinformatics and computational tools as well as experimental strategies for the prediction and assignment of functions (in vitro activities and in vivo physiological/metabolic roles) to uncharacterized enzymes discovered in genome projects. Protein sequence similarity networks (SSNs) are visually powerful tools for analyzing sequence relationships in protein families (H.J. Atkinson, J.H. Morris, T.E. Ferrin, and P.C. Babbitt, PLoS One 2009, 4, e4345). However, the members of the biological/biomedical community have not had access to the capability to generate SSNs for their "favorite" protein families. In this article we announce the EFI-EST (Enzyme Function Initiative-Enzyme Similarity Tool) web tool (http://efi.igb.illinois.edu/efi-est/) that is available without cost for the automated generation of SSNs by the community. The tool can create SSNs for the "closest neighbors" of a user-supplied protein sequence from the UniProt database (Option A) or of members of any user-supplied Pfam and/or InterPro family (Option B). We provide an introduction to SSNs, a description of EFI-EST, and a demonstration of the use of EFI-EST to explore sequence-function space in the OMP decarboxylase superfamily (PF00215). This article is designed as a tutorial that will allow members of the community to use the EFI-EST web tool for exploring sequence/function space in protein families.

Rate and Equilibrium Constants for an Enzyme Conformational Change during Catalysis by Orotidine 5'-Monophosphate Decarboxylase.

The caged complex between orotidine 5'-monophosphate decarboxylase (ScOMPDC) and 5-fluoroorotidine 5'-monophosphate (FOMP) undergoes decarboxylation ∼300 times faster than the caged complex between ScOMPDC and the physiological substrate, orotidine 5'-monophosphate (OMP). Consequently, the enzyme conformational changes required to lock FOMP at a protein cage and release product 5-fluorouridine 5'-monophosphate (FUMP) are kinetically significant steps. The caged form of ScOMPDC is stabilized by interactions between the side chains from Gln215, Tyr217, and Arg235 and the substrate phosphodianion. The control of these interactions over the barrier to the binding of FOMP and the release of FUMP was probed by determining the effect of all combinations of single, double, and triple Q215A, Y217F, and R235A mutations on kcat/Km and kcat for turnover of FOMP by wild-type ScOMPDC; its values are limited by the rates of substrate binding and product release, respectively. The Q215A and Y217F mutations each result in an increase in kcat and a decrease in kcat/Km, due to a weakening of the protein-phosphodianion interactions that favor fast product release and slow substrate binding. The Q215A/R235A mutation causes a large decrease in the kinetic parameters for ScOMPDC-catalyzed decarboxylation of OMP, which are limited by the rate of the decarboxylation step, but much smaller decreases in the kinetic parameters for ScOMPDC-catalyzed decarboxylation of FOMP, which are limited by the rate of enzyme conformational changes. By contrast, the Y217A mutation results in large decreases in kcat/Km for ScOMPDC-catalyzed decarboxylation of both OMP and FOMP, because of the comparable effects of this mutation on rate-determining decarboxylation of enzyme-bound OMP and on the rate-determining enzyme conformational change for decarboxylation of FOMP. We propose that kcat = 8.2 s(-1) for decarboxylation of FOMP by the Y217A mutant is equal to the rate constant for cage formation from the complex between FOMP and the open enzyme, that the tyrosyl phenol group stabilizes the closed form of ScOMPDC by hydrogen bonding to the substrate phosphodianion, and that the phenyl group of Y217 and F217 facilitates formation of the transition state for the rate-limiting conformational change. An analysis of kinetic data for mutant enzyme-catalyzed decarboxylation of OMP and FOMP provides estimates for the rate and equilibrium constants for the conformational change that traps FOMP at the enzyme active site.

Experimental strategies for functional annotation and metabolism discovery: targeted screening of solute binding proteins and unbiased panning of metabolomes.

The rate at which genome sequencing data is accruing demands enhanced methods for functional annotation and metabolism discovery. Solute binding proteins (SBPs) facilitate the transport of the first reactant in a metabolic pathway, thereby constraining the regions of chemical space and the chemistries that must be considered for pathway reconstruction. We describe high-throughput protein production and differential scanning fluorimetry platforms, which enabled the screening of 158 SBPs against a 189 component library specifically tailored for this class of proteins. Like all screening efforts, this approach is limited by the practical constraints imposed by construction of the library, i.e., we can study only those metabolites that are known to exist and which can be made in sufficient quantities for experimentation. To move beyond these inherent limitations, we illustrate the promise of crystallographic- and mass spectrometric-based approaches for the unbiased use of entire metabolomes as screening libraries. Together, our approaches identified 40 new SBP ligands, generated experiment-based annotations for 2084 SBPs in 71 isofunctional clusters, and defined numerous metabolic pathways, including novel catabolic pathways for the utilization of ethanolamine as sole nitrogen source and the use of d-Ala-d-Ala as sole carbon source. These efforts begin to define an integrated strategy for realizing the full value of amassing genome sequence data.

A unique cis-3-hydroxy-l-proline dehydratase in the enolase superfamily.

The genome of Labrenzia aggregata IAM 12614 encodes an uncharacterized member of the muconate lactonizing enzyme (MLE) subgroup of the enolase superfamily (UniProt ID A0NXQ8 ). The gene encoding A0NXQ8 is located between genes that encode members of the proline racemase superfamily, 4R-hydroxyproline 2-epimerase (UniProt ID A0NXQ7 ; 4HypE) and trans-3-hydroxy-l-proline dehydratase (UniProt ID A0NXQ9 ; t3LHypD). A0NXQ8 was screened with a library of proline analogues; two reactions were observed with cis-3-hydroxy-l-proline (c3LHyp), competing 2-epimerization to trans-3-hydroxy-d-proline (1,1-proton transfer) and dehydration to Δ(1)-pyrroline-2-carboxylate (ß-elimination; c3LHyp dehydratase), with eventual total dehydration. The genome context encoding A0NXQ8 both (1) confirms its novel c3LHyp dehydratase function and (2) provides evidence for metabolic pathways that allow L. aggregata to utilize several isomeric 3- and 4-hydroxyprolines as sole carbon sources.

A General Strategy for the Discovery of Metabolic Pathways: d-Threitol, l-Threitol, and Erythritol Utilization in Mycobacterium smegmatis.

We describe a general integrated bioinformatic and experimental strategy to discover the in vitro enzymatic activities and in vivo functions (metabolic pathways) of uncharacterized enzymes discovered in microbial genome projects using the ligand specificities of the solute binding proteins (SBPs) for ABC transporters. Using differential scanning fluorimetry, we determined that the SBP for an ABC transporter encoded by the genome of Mycobacterium smegmatis is stabilized by d-threitol. Using sequence similarity networks and genome neighborhood networks to guide selection of target proteins for pathway enzymes, we applied both in vitro and in vivo experimental approaches to discover novel pathways for catabolism of d-threitol, l-threitol, and erythritol.

Homology models guide discovery of diverse enzyme specificities among dipeptide epimerases in the enolase superfamily.

The rapid advance in genome sequencing presents substantial challenges for protein functional assignment, with half or more of new protein sequences inferred from these genomes having uncertain assignments. The assignment of enzyme function in functionally diverse superfamilies represents a particular challenge, which we address through a combination of computational predictions, enzymology, and structural biology. Here we describe the results of a focused investigation of a group of enzymes in the enolase superfamily that are involved in epimerizing dipeptides. The first members of this group to be functionally characterized were Ala-Glu epimerases in Eschericiha coli and Bacillus subtilis, based on the operon context and enzymological studies; these enzymes are presumed to be involved in peptidoglycan recycling. We have subsequently studied more than 65 related enzymes by computational methods, including homology modeling and metabolite docking, which suggested that many would have divergent specificities;, i.e., they are likely to have different (unknown) biological roles. In addition to the Ala-Phe epimerase specificity reported previously, we describe the prediction and experimental verification of: (i) a new group of presumed Ala-Glu epimerases; (ii) several enzymes with specificity for hydrophobic dipeptides, including one from Cytophaga hutchinsonii that epimerizes D-Ala-D-Ala; and (iii) a small group of enzymes that epimerize cationic dipeptides. Crystal structures for certain of these enzymes further elucidate the structural basis of the specificities. The results highlight the potential of computational methods to guide experimental characterization of enzymes in an automated, large-scale fashion.