Glycan sulfation is a widespread postglycosylation modification crucial for modulating biological functions including cellular adhesion, signaling, and bacterial colonization. 6-Sulfo-beta-GlcNAcases are a class of enzyme that alters sulfation patterns. Such changes in sulfation patterns are linked to diseases such as bowel inflammation, colitis, and cancer. Despite their significance, 6-sulfo-beta-GlcNAcases, which cleave beta-linked 6-sulfo-N-acetylglucosamine (6S-GlcNAc), have been but rarely identified. This scarcity results mainly from the short, diverse, and distinctive sulfate-binding motifs required for recognition of the 6-sulfate group in 6S-GlcNAc in addition to the conserved GH20 family features. In this study, we discovered 6-sulfo-beta-GlcNAcases and assigned two novel sulfate-binding motifs by the use of comparative genomics, structural predictions, and activity-based screening. Our findings expand the known microbiota capable of degrading sulfated glycans and add significant enzymes to the tool kit for analysis and synthesis of sulfated oligosaccharides.

SUMMARY Deazaguanine modifications play multifaceted roles in the molecular biology of DNA and tRNA, shaping diverse yet essential biological processes, including the nuanced fine-tuning of translation efficiency and the intricate modulation of codon-anticodon interactions. Beyond their roles in translation, deazaguanine modifications contribute to cellular stress resistance, self-nonself discrimination mechanisms, and host evasion defenses, directly modulating the adaptability of living organisms. Deazaguanine moieties extend beyond nucleic acid modifications, manifesting in the structural diversity of biologically active natural products. Their roles in fundamental cellular processes and their presence in biologically active natural products underscore their versatility and pivotal contributions to the intricate web of molecular interactions within living organisms. Here, we discuss the current understanding of the biosynthesis and multifaceted functions of deazaguanines, shedding light on their diverse and dynamic roles in the molecular landscape of life.

Nature has evolved biosynthetic pathways to molecules possessing reactive warheads that inspired the development of many therapeutic agents, including penicillin antibiotics. Peptides armed with electrophilic warheads have proven to be particularly effective covalent inhibitors, providing essential antimicrobial, antiviral and anticancer agents. Here we provide a full characterization of the pathways that nature deploys to assemble peptides with β-lactone warheads, which are potent proteasome inhibitors with promising anticancer activity. Warhead assembly involves a three-step cryptic methylation sequence, which is likely required to reduce unfavorable electrostatic interactions during the sterically demanding β-lactonization. Amide-bond synthetase and adenosine triphosphate (ATP)-grasp enzymes couple amino acids to the β-lactone warhead, generating the bioactive peptide products. After reconstituting the entire pathway to β-lactone peptides in vitro, we go on to deliver a diverse range of analogs through enzymatic cascade reactions. Our approach is more efficient and cleaner than the synthetic methods currently used to produce clinically important warhead-containing peptides.

Glycan sulfation is an important modification supporting the functionalities of many proteins in biology. Exo- acting 6S-GlcNAcases from human microbiota are glycosidases that participate in the removal of 6-sulfo-GlcNAc from host glycans and thereby play an important role in human health and disease. Nonetheless, mechanisms underlying their ability to recognize the sulfate group remain poorly understood. Using structural and kinetic analyses, we here reveal the catalytically important amino acids directly involved in the recognition and cleavage of 6S-GlcNAc, but not of 6-phospho-GlcNAc, in BbhII from Bifidobacterium bifidum, Bt4394 from Bacteroides thetaiotaomicron, and SGL from Prevotella spp. The defining features of their sulfate recognition motifs underpin a genomic enzymological exploration of 6S-GlcNAcases to identify a wider range of human health-associated bacterial species having 6S-GlcNAcase activity. Our data provide significant insights into distinct molecular mechanisms of sulfated sugar recognition employed by 6S-GlcNAcases from both Gram-positive and Gramnegative bacteria along with valuable information for the exploration of extensive interactions between microbiota and their host glycans.

The Enzyme Function Initiative (EFI) provides a web resource with “genomic enzymology” web tools to leverage the protein (UniProt) and genome (European Nucleotide Archive; ENA; https://www.ebi.ac.uk/ena/) databases to assist the assignment of in vitro enzymatic activities and in vivo metabolic functions to uncharacterized enzymes (https://efi.igb.illinois.edu/). The tools enable (1) exploration of sequence-function space in enzyme families using sequence similarity networks (SSNs; EFI-EST), (2) easy access to genome context for bacterial, archaeal, and fungal proteins in the SSN clusters so that isofunctional families can be identified and their functions inferred from genome context (EFI-GNT); and (3) determination of the abundance of SSN clusters in NIH Human Metagenome Project metagenomes using chemically guided functional profiling (EFI-CGFP). We describe enhancements that enable SSNs to be generated from taxonomy categories, allowing higher resolution analyses of sequence-function space; we provide examples of the generation of taxonomy category-specific SSNs.

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.

The continuing expansion of protein and genome sequence databases is an opportunity to identify novel enzymes with biotechnological applications. Whether applied to enzymology, chemical biology, systems biology, and microbiology, database mining must be ‘user-friendly’ so that experimentalists can devise focused strategies to discover the in vitro activities and in vivo functions of uncharacterized enzymes. We developed a suite of genomic enzymology tools (https://efi.igb.illinois.edu/) to (1) generate sequence similarity networks (SSNs) for exploration of sequence-function space in protein families (EFI-EST) and (2) provide genome context for members of protein families (EFI-GNT). Integrated analysis of this complementary information allows to generate testable hypotheses about new functions. After a brief overview of EFI-EST and EFI-GNT, we describe applications that illustrate their use.

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 preQ(0) and preQ(1) also occurs. With the exception of Escherichia coli YhhQ, shown to transport preQ(0) and preQ(1), 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 preQ(1) 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.

Summary The examination of gene neighborhood is an integral part of comparative genomics but no tools to produce publication quality graphics of gene clusters are available. Gene Graphics is a straightforward web application for creating such visuals. Supported inputs include National Center for Biotechnology Information gene and protein identifiers with automatic fetching of neighboring information, GenBank files and data extracted from the SEED database. Gene representations can be customized for many parameters including gene and genome names, colors and sizes. Gene attributes can be copied and pasted for rapid and user‐friendly customization of homologous genes between species. In addition to Portable Network Graphics and Scalable Vector Graphics, produced representations can be exported as Tagged Image File Format or Encapsulated PostScript, formats that are standard for publication. Hands‐on tutorials with real life examples inspired from publications are available for training. Availability and implementation Gene Graphics is freely available at https://katlabs.cc/genegraphics/ and source code is hosted at https://github.com/katlabs/genegraphics. Supplementary information Supplementary data are available at Bioinformatics online.

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.

The protein databases contain an exponentially growing number of sequences as a result of the recent increase in ease and decrease in cost of genome sequencing. The rate of data accumulation far exceeds the rate of functional studies, producing an increase in genomic 'dark matter', sequences for which no precise and validated function is defined. Publicly accessible, that is 'democratized,' genomic enzymology web tools are essential to leverage the protein and genome databases for discovery of the in vitro activities and in vivo functions of novel enzymes and proteins belonging to the dark matter. In this review, we discuss the use of web tools that have proven successful for functional assignment. We also describe a mechanism for ensuring the capture of published functional data so that the quality of both curated and automated annotations transfer can be improved.

Queuosine (Q) is a complex modification of the wobble base in tRNAs with GUN anticodons. The full Q biosynthesis pathway has been elucidated in Escherichia coli. FolE, QueD, QueE and QueC are involved in the conversion of guanosine triphosphate (GTP) to 7-cyano-7-deazaguanine (preQ(0)), an intermediate of increasing interest for its central role in tRNA and DNA modification and secondary metabolism. QueF then reduces preQ(0) to 7-aminomethyl-7-deazaguanine (preQ(1)). PreQ(1) is inserted into tRNAs by tRNA guanine((34)) transglycosylase (TGT). The inserted base preQ(1) is finally matured to Q by two additional steps involving QueA and QueG or QueH. Most Eubacteria harbor the full set of Q synthesis genes and are predicted to synthesize Q de novo. However, some bacteria only encode enzymes involved in the second half of the pathway downstream of preQ(0) synthesis, including the signature enzyme TGT. Different patterns of distribution of the queF, tgt, queA and queG or queH genes are observed, suggesting preQ(0), preQ(1) or even the queuine base being salvaged in specific organisms. Such salvage pathways require the existence of specific 7-deazapurine transporters that have yet to be identified. The COG1738 family was identified as a candidate for a missing preQ(0)/preQ(1) transporter in prokaryotes, by comparative genomics analyses. The existence of Q precursor salvage was confirmed for the first time in bacteria, in vivo, through an indirect assay. The involvement of the COG1738 in salvage of a Q precursor was experimentally validated in Escherichia coli, where it was shown that the COG1738 family member YhhQ is essential for preQ(0) transport.

The reduction of epoxyqueuosine (oQ) is the last step in the synthesis of the tRNA modification queuosine (Q). While the epoxyqueuosine reductase (EC 1.17.99.6) enzymatic activity was first described 30 years ago, the encoding gene queG was only identified in Escherichia coli in 2011. Interestingly, queG is absent from a large number of sequenced genomes that harbor Q synthesis or salvage genes, suggesting the existence of an alternative epoxyqueuosine reductase in these organisms. By analyzing phylogenetic distributions, physical gene clustering, and fusions, members of the Domain of Unknown Function 208 (DUF208) family were predicted to encode for an alternative epoxyqueuosine reductase. This prediction was validated with genetic methods. The Q modification is present in Lactobacillus salivarius, an organism missing queG but harboring the duf208 gene. Acinetobacter baylyi ADP1 is one of the few organisms that harbor both QueG and DUF208, and deletion of both corresponding genes was required to observe the absence of Q and the accumulation of oQ in tRNA. Finally, the conversion oQ to Q was restored in an E. coli queG mutant by complementation with plasmids harboring duf208 genes from different bacteria. Members of the DUF208 family are not homologous to QueG enzymes, and thus, duf208 is a non-orthologous replacement of queG. We propose to name DUF208 encoding genes as queH. While QueH contains conserved cysteines that could be involved in the coordination of a Fe/S center in a similar fashion to what has been identified in QueG, no cobalamin was identified associated with recombinant QueH protein.

To synthesize the cofactor thiamin diphosphate (ThDP), plants must first hydrolyze thiamin monophosphate (ThMP) to thiamin, but dedicated enzymes for this hydrolysis step were unknown and widely doubted to exist. The classical thiamin-requiring th2-1 mutation in Arabidopsis thaliana was shown to reduce ThDP levels by half and to increase ThMP levels 5-fold, implying that the THIAMIN REQUIRING2 (TH2) gene product could be a dedicated ThMP phosphatase. Genomic and transcriptomic data indicated that TH2 corresponds to At5g32470, encoding a HAD (haloacid dehalogenase) family phosphatase fused to a TenA (thiamin salvage) family protein. Like the th2-1 mutant, an insertional mutant of At5g32470 accumulated ThMP, and the thiamin requirement of the th2-1 mutant was complemented by wild-type At5g32470 Complementation tests in Escherichia coli and enzyme assays with recombinant proteins confirmed that At5g32470 and its maize (Zea mays) orthologs GRMZM2G148896 and GRMZM2G078283 are ThMP-selective phosphatases whose activity resides in the HAD domain and that the At5g32470 TenA domain has the expected thiamin salvage activity. In vitro and in vivo experiments showed that alternative translation start sites direct the At5g32470 protein to the cytosol and potentially also to mitochondria. Our findings establish that plants have a dedicated ThMP phosphatase and indicate that modest (50%) ThDP depletion can produce severe deficiency symptoms.

Gene duplication followed by mutation is a classic mechanism of neofunctionalization, producing gene families with functional diversity. In some cases, a single point mutation is sufficient to change the substrate specificity and/or the chemistry performed by an enzyme, making it difficult to accurately separate enzymes with identical functions from homologs with different functions. Because sequence similarity is often used as a basis for assigning functional annotations to genes, non-isofunctional gene families pose a great challenge for genome annotation pipelines. Here we describe how integrating evolutionary and functional information such as genome context, phylogeny, metabolic reconstruction and signature motifs may be required to correctly annotate multifunctional families. These integrative analyses can also lead to the discovery of novel gene functions, as hints from specific subgroups can guide the functional characterization of other members of the family. We demonstrate how careful manual curation processes using comparative genomics can disambiguate subgroups within large multifunctional families and discover their functions. We present the COG0720 protein family as a case study. We also discuss strategies to automate this process to improve the accuracy of genome functional annotation pipelines.

The penultimate step of thiamin diphosphate (ThDP) synthesis in plants and many bacteria is dephosphorylation of thiamin monophosphate (ThMP). Non-specific phosphatases have been thought to mediate this step and no genes encoding specific ThMP phosphatases (ThMPases) are known. Comparative genomic analysis uncovered bacterial haloacid dehalogenase (HAD) phosphatase family genes (from subfamilies IA and IB) that cluster on the chromosome with, or are fused to, thiamin synthesis genes and are thus candidates for the missing phosphatase (ThMPase). Three typical candidates (from Anaerotruncus colihominis, Dorea longicatena and Syntrophomonas wolfei) were shown to have efficient in vivo ThMPase activity by expressing them in an Escherichia coli strain engineered to require an active ThMPase for growth. In vitro assays confirmed that these candidates all preferred ThMP to any of 45 other phosphate ester substrates tested. An Arabidopsis thaliana ThMPase homologue (At4g29530) of unknown function whose expression pattern and compartmentation fit with a role in ThDP synthesis was shown to have in vivo ThMPase activity in E. coli and to prefer ThMP to any other substrate tested. However, insertional inactivation of the At4g29530 gene did not affect growth or the levels of thiamin or its phosphates, indicating that Arabidopsis has at least one other ThMPase gene. The Zea mays orthologue of At4g29530 (GRMZM2G035134) was also shown to have ThMPase activity. These data identify HAD genes specifying the elusive ThMPase activity, indicate that ThMPases are substrate-specific rather than general phosphatases and suggest that different evolutionary lineages have recruited ThMPases independently from different branches of the HAD family.

Background Gene fusions are the most powerful type of in silico -derived functional associations. However, many fusion compilations were made when <100 genomes were available, and algorithms for identifying fusions need updating to handle the current avalanche of sequenced genomes. The availability of a large fusion dataset would help probe functional associations and enable systematic analysis of where and why fusion events occur. Results Here we present a systematic analysis of fusions in prokaryotes. We manually generated two training sets: (i) 121 fusions in the model organism Escherichia coli ; (ii) 131 fusions found in B vitamin metabolism. These sets were used to develop a fusion prediction algorithm that captured the training set fusions with only 7 % false negatives and 50 % false positives, a substantial improvement over existing approaches. This algorithm was then applied to identify 3.8 million potential fusions across 11,473 genomes. The results of the analysis are available in a searchable database at http://modelseed.org/projects/fusions/ . A functional analysis identified 3,000 reactions associated with frequent fusion events and revealed areas of metabolism where fusions are particularly prevalent. Conclusions Customary definitions of fusions were shown to be ambiguous, and a stricter one was proposed. Exploring the genes participating in fusion events showed that they most commonly encode transporters, regulators, and metabolic enzymes. The major rationales for fusions between metabolic genes appear to be overcoming pathway bottlenecks, avoiding toxicity, controlling competing pathways, and facilitating expression and assembly of protein complexes. Finally, our fusion dataset provides powerful clues to decipher the biological activities of domains of unknown function.

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.

DUF89 family proteins occur widely in both prokaryotes and eukaryotes, but their functions are unknown. Here we define three DUF89 subfamilies (I, II, and III), with subfamily II being split into stand-alone proteins and proteins fused to pantothenate kinase (PanK). We demonstrated that DUF89 proteins have metal-dependent phosphatase activity against reactive phosphoesters or their damaged forms, notably sugar phosphates (subfamilies II and III), phosphopantetheine and its S-sulfonate or sulfonate (subfamily II-PanK fusions), and nucleotides (subfamily I). Genetic and comparative genomic data strongly associated DUF89 genes with phosphoester metabolism. The crystal structure of the yeast (Saccharomyces cerevisiae) subfamily III protein YMR027W revealed a novel phosphatase active site with fructose 6-phosphate and Mg(2+) bound near conserved signature residues Asp254 and Asn255 that are critical for activity. These findings indicate that DUF89 proteins are previously unrecognized hydrolases whose characteristic in vivo function is to limit potentially harmful buildups of normal or damaged phosphometabolites.

The haloacid dehalogenase (HAD)-like enzymes comprise a large superfamily of phosphohydrolases present in all organisms. The Saccharomyces cerevisiae genome encodes at least 19 soluble HADs, including 10 uncharacterized proteins. Here, we biochemically characterized 13 yeast phosphatases from the HAD superfamily, which includes both specific and promiscuous enzymes active against various phosphorylated metabolites and peptides with several HADs implicated in detoxification of phosphorylated compounds and pseudouridine. The crystal structures of four yeast HADs provided insight into their active sites, whereas the structure of the YKR070W dimer in complex with substrate revealed a composite substrate-binding site. Although the S. cerevisiae and Escherichia coli HADs share low sequence similarities, the comparison of their substrate profiles revealed seven phosphatases with common preferred substrates. The cluster of secondary substrates supporting significant activity of both S. cerevisiae and E. coli HADs includes 28 common metabolites that appear to represent the pool of potential activities for the evolution of novel HAD phosphatases. Evolution of novel substrate specificities of HAD phosphatases shows no strict correlation with sequence divergence. Thus, evolution of the HAD superfamily combines the conservation of the overall substrate pool and the substrate profiles of some enzymes with remarkable biochemical and structural flexibility of other superfamily members.