Oxidative opening of the aromatic ring: Tracing the natural history of a large superfamily of dioxygenase domains and their relatives

Analysis of Rhizonin Biosynthesis Reveals Origin of Pharmacophoric Furylalanine Moieties in Diverse Cyclopeptides

Rhizonin A and B are hepatotoxic cyclopeptides produced by bacterial endosymbionts (Mycetohabitans endofungorum) of the fungus Rhizopus microsporus. Their toxicity critically depends on the presence of 3-furylalanine (Fua) residues, which also occur in pharmaceutically relevant cyclopeptides of the endolide and bingchamide families. The biosynthesis and incorporation of Fua by non-ribosomal peptide synthetases (NRPS), however, has remained elusive. By genome sequencing and gene inactivation we elucidated the gene cluster responsible for rhizonin biosynthesis. A suite of isotope labeling experiments identified tyrosine and l-DOPA as Fua precursors and provided the first mechanistic insight. Bioinformatics, mutational analysis and heterologous reconstitution identified dioxygenase RhzB as necessary and sufficient for Fua formation. RhzB is a novel type of heme-dependent aromatic oxygenases (HDAO) that enabled the discovery of the bingchamide biosynthesis gene cluster through genome mining.

EFI-EST, EFI-GNT, and EFI-CGFP: Enzyme Function Initiative (EFI) Web Resource for Genomic Enzymology Tools

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.

Large-scale phylogenomics of aquatic bacteria reveal molecular mechanisms for adaptation to salinity

The crossing of environmental barriers poses major adaptive challenges. Rareness of freshwater-marine transitions separates the bacterial communities, but how these are related to brackish counterparts remains elusive, as do the molecular adaptations facilitating cross-biome transitions. We conducted large-scale phylogenomic analysis of freshwater, brackish, and marine quality-filtered metagenome-assembled genomes (11,248). Average nucleotide identity analyses showed that bacterial species rarely existed in multiple biomes. In contrast, distinct brackish basins cohosted numerous species, but their intraspecific population structures displayed clear signs of geographic separation. We further identified the most recent cross-biome transitions, which were rare, ancient, and most commonly directed toward the brackish biome. Transitions were accompanied by systematic changes in amino acid composition and isoelectric point distributions of inferred proteomes, which evolved over millions of years, as well as convergent gains or losses of specific gene functions. Therefore, adaptive challenges entailing proteome reorganization and specific changes in gene content constrains the cross-biome transitions, resulting in species-level separation between aquatic biomes.

How Single Amino Acid Substitutions Can Disrupt a Protein Hetero-Dimer Interface: Computational and Experimental Studies of the LigAB Dioxygenase from Sphingobium sp. Strain SYK-6

Protocatechuate 4,5-dioxygenase (LigAB) is a heterodimeric enzyme that catalyzes the dioxygenation of multiple lignin derived aromatic compounds. The active site of LigAB is at the heterodimeric interface, with specificity conferred by the alpha subunit and catalytic residues contributed by the beta subunit. Previous research has indicated that the phenylalanine at the 103 position of the alpha subunit (F103α) controls selectivity for the C5 position of the aromatic substrates, and mutations of this residue can enhance the rate of catalysis for substrates with larger functional groups at this position. While several of the mutations to this position (Valine, V; Threonine, T; Leucine, L; and Histidine, H) were catalytically active, other mutations (Alanine, A; and Serine, S) were found to have reduced dimer interface affinity, leading to challenges in copurifing the catalytically active enzyme complex under high salt conditions. In this study, we aimed to experimentally and computationally interrogate residues at the dimer interface to discern the importance of position 103α for maintaining the integrity of the heterodimer. Molecular dynamic simulations and electrophoretic mobility assays revealed a preference for nonpolar/aromatic amino acids in this position, suggesting that while substitutions to polar amino acids may produce a dioxygenase with a useful substrate utilization profile, those considerations may be off-set by potential destabilization of the catalytically active oligomer. Understanding the dimerization of LigAB provides insight into the multimeric proteins within the largely uncharacterized superfamily and characteristics to consider when engineering proteins that can degrade lignin efficiently. These results shed light on the challenges associated with engineering proteins for broader substrate specificity.

Subtle sequence differences between two interacting σ54 -dependent regulators lead to different activation mechanisms

In the strictly anaerobic nitrate reducing bacterium Aromatoleum anaerobium, degradation of 1,3-dihydroxybenzene (1,3-DHB, resorcinol) is controlled by two bacterial enhancer-binding proteins (bEBPs), RedR1 and RedR2, which regulate the transcription of three σ⁵⁴ -dependent promoters controlling expression of the pathway. RedR1 and RedR2 are identical over their length except for their N-terminal tail which differ in sequence and length (six and eight residues, respectively), a single change in their N-terminal domain (NTD), and nine non-identical residues in their C-terminal domain (CTD). Their NTD is composed of a GAF and a PAS domain connected by a linker helix. We show that each regulator is controlled by a different mechanism: whilst RedR1 responds to the classical NTD-mediated negative regulation that is released by the presence of its effector, RedR2 activity is constitutive and controlled through interaction with BtdS, an integral membrane subunit of hydroxyhydroquinone dehydrogenase carrying out the second step in 1,3-DHB degradation. BtdS sequesters the RedR2 regulator to the membrane through its NTD, where a four-Ile track in the PAS domain, interrupted by a Thr in RedR1, and the N-terminal tail are involved. The presence of 1,3-DHB, which is metabolized to hydroxybenzoquinone, releases RedR2 from the membrane. Most bEBPs assemble into homohexamers to activate transcription; we show that hetero-oligomer formation between RedR1 and RedR2 is favoured over homo-oligomers. However, either an NTD-truncated version of RedR1 or a full-length RedR2 are capable of promoter activation on their own, suggesting they should assemble into homohexamers in vivo. We show that promoter DNA behaves as an allosteric effector through binding the CTD to control ΔNTD-RedR1 multimerization and activity. Overall, the regulation of the 1,3-DHB anaerobic degradation pathway can be described as a novel mode of bEBP activation and assembly.

Reannotation of Fly Amanita l-DOPA Dioxygenase Gene Enables Its Cloning and Heterologous Expression

The l-DOPA dioxygenase of Amanita muscaria (AmDODA) participates in the biosynthesis of betalain- and hygroaurin-type natural pigments. AmDODA is encoded by the dodA gene, whose DNA sequence was inferred from cDNA and gDNA libraries almost 30 years ago. However, reports on its heterologous expression rely on either the original 5'-truncated cDNA plasmid or artificial gene synthesis. We provide unequivocal evidence that the heterologous expression of AmDODA from A. muscaria specimens is not possible by using the coding sequence previously inferred for dodA. Here, we rectify and reannotate the full-length coding sequence for AmDODA and express a 205-aa His-tagged active enzyme, which was used to produce the l-DOPA hygroaurin, a rare fungal pigment. Moreover, AmDODA and other isozymes from bacteria were submitted to de novo folding using deep learning algorithms, and their putative active sites were inferred and compared. The wide catalytic pocket of AmDODA and the presence of the His-His-His and His-His-Asp motifs can provide insight into the dual cleavage of l-DOPA at positions 2,3 and 4,5 as per the mechanism proposed for nonheme dioxygenases.

Characterization of Protocatechuate 4,5-Dioxygenase from Pseudarthrobacter phenanthrenivorans Sphe3 and In Situ Reaction Monitoring in the NMR Tube

The current study aims at the functional and kinetic characterization of protocatechuate (PCA) 4,5-dioxygenase (PcaA) from Pseudarthrobacter phenanthrenivorans Sphe3. This is the first single subunit Type II dioxygenase characterized in Actinobacteria. RT-PCR analysis demonstrated that pcaA and the adjacent putative genes implicated in the PCA meta-cleavage pathway comprise a single transcriptional unit. The recombinant PcaA is highly specific for PCA and exhibits Michaelis-Menten kinetics with Km and Vmax values of 21 ± 1.6 μM and 44.8 ± 4.0 U × mg-1, respectively, in pH 9.5 and at 20 °C. PcaA also converted gallate from a broad range of substrates tested. The enzymatic reaction products were identified and characterized, for the first time, through in situ biotransformation monitoring inside an NMR tube. The PCA reaction product demonstrated a keto-enol tautomerization, whereas the gallate reaction product was present only in the keto form. Moreover, the transcriptional levels of pcaA and pcaR (gene encoding a LysR-type regulator of the pathway) were also determined, showing an induction when cells were grown on PCA and phenanthrene. Studying key enzymes in biodegradation pathways is significant for bioremediation and for efficient biocatalysts development.

Shape-function of a novel metapyrocatechase, RW4-MPC: Metagenomics to SAXS data based insight into deciphering regulators of function

The oxygenases have attracted considerable attention in enzyme-mediated bioremediation of xenobiotic compounds due to their high specificity, cost-effectiveness, and targeted field applications. Here, we performed a functional metagenomics approach to cope with cultivability limitations to isolate a novel extradiol dioxygenase. Fosmid clone harboring dioxygenase gene was sequenced and analyzed by bioinformatics tools. One ring-cleaving dioxygenase RW4-MPC (metapyrocatechase) was purified and characterized to examine its degradation efficiency. The RW4-MPC was significantly active in the temperature and pH range of 5 to 40 °C, and 7-10, respectively, with an optimum temperature of 25 °C and pH 8. To gain insight into observed differential activity, Small-Angle X-ray Scattering (SAXS) data of the protein samples were analyzed, which brought forth that the RW4-MPC molecules form tight globular tetramers in solution. This native association was stable till 35 °C, and protein started to associate at higher temperatures, explaining heat-induced loss of function. Similarly, RW4-MPC aggregated or lost globular profile below pH 7 or at pH 10, respectively. The kinetic parameters showed the six folds high catalytic efficiency of RW4-MPC towards 2,3-dihydroxy biphenyl than catechol and its derivatives. RW4-MPC molecules showed remarkable retention of functionality in hypersaline conditions with more than 70% activity in a buffer having 3 M NaCl concentration. In concordance, SAXS data analysis showed retention of functional shape profile in hypersaline conditions. The halotolerant and oxygen insensitive nature of this enzyme makes it a potential candidate for bioremediation.

L-DOPA Dioxygenase Activity on 6-Substituted Dopamine Analogues

Dioxygenase enzymes are essential protein catalysts for the breakdown of catecholic rings, structural components of plant woody tissue. This powerful chemistry is used in nature to make antibiotics and other bioactive materials or degrade plant material, but we have a limited understanding of the breadth and depth of substrate space for these potent catalysts. Here we report steady-state and pre-steady-state kinetic analysis of dopamine derivatives substituted at the 6-position as substrates of L-DOPA dioxygenase, and an analysis of that activity as a function of the electron-withdrawing nature of the substituent. Steady-state and pre-steady-state kinetic data demonstrate the dopamines are impaired in binding and catalysis with respect to the cosubstrate molecular oxygen, which likely afforded spectroscopic observation of an early reaction intermediate, the semiquinone of dopamine. The reaction pathway of dopamine in the pre-steady state is consistent with a nonproductive mode of binding of oxygen at the active site. Despite these limitations, L-DOPA dioxygenase is capable of binding all of the dopamine derivatives and catalyzing multiple turnovers of ring cleavage for dopamine, 6-bromodopamine, 6-carboxydopamine, and 6-cyanodopamine. 6-Nitrodopamine was a single-turnover substrate. The variety of substrates accepted by the enzyme is consistent with an interplay of factors, including the capacity of the active site to bind large, negatively charged groups at the 6-position and the overall oxidizability of each catecholamine, and is indicative of the utility of extradiol cleavage in semisynthetic and bioremediation applications.

How Do Shipworms Eat Wood? Screening Shipworm Gill Symbiont Genomes for Lignin-Modifying Enzymes

Shipworms are ecologically and economically important mollusks that feed on woody plant material (lignocellulosic biomass) in marine environments. Digestion occurs in a specialized cecum, reported to be virtually sterile and lacking resident gut microbiota. Wood-degrading CAZymes are produced both endogenously and by gill endosymbiotic bacteria, with extracellular enzymes from the latter being transported to the gut. Previous research has predominantly focused on how these animals process the cellulose component of woody plant material, neglecting the breakdown of lignin – a tough, aromatic polymer which blocks access to the holocellulose components of wood. Enzymatic or non-enzymatic modification and depolymerization of lignin has been shown to be required in other wood-degrading biological systems as a precursor to cellulose deconstruction. We investigated the genomes of five shipworm gill bacterial symbionts obtained from the Joint Genome Institute Integrated Microbial Genomes and Microbiomes Expert Review for the production of lignin-modifying enzymes, or ligninases. The genomes were searched for putative ligninases using the Joint Genome Institute’s Function Profile tool and blastp analyses. The resulting proteins were then modeled using SWISS-MODEL. Although each bacterial genome possessed at least four predicted ligninases, the percent identities and protein models were of low quality and were unreliable. Prior research demonstrates limited endogenous ability of shipworms to modify lignin at the chemical/molecular level. Similarly, our results reveal that shipworm bacterial gill-symbiont enzymes are unlikely to play a role in lignin modification during lignocellulose digestion in the shipworm gut. This suggests that our understanding of how these keystone organisms digest and process lignocellulose is incomplete, and further research into non-enzymatic and/or other unknown mechanisms for lignin modification is required.

Potential applications of halophilic microorganisms for biological treatment of industrial process brines contaminated with aromatics

Saline wastewater contaminated with aromatic compounds can be frequently found in various industrial sectors. Those compounds need to be degraded before reuse of wastewater in other process steps or release to the environment. Halophiles have been reported to efficiently degrade aromatics, but their application to treat industrial wastewater is rare. Halophilic processes for industrial wastewater treatment need to satisfy certain requirements: a continuous process mode, low operational expenditures, suitable reactor systems and a monitoring and control strategy. The aim of this review is to provide an overview of halophilic microorganisms, principles of aromatic biodegradation, and sources of saline wastewater containing aromatics and other contaminants. Finally, process examples for halophilic wastewater treatment and potential process monitoring strategies are discussed. To further illustrate the significant potential of halophiles for saline wastewater treatment and to facilitate development of ready-to-implement processes, future research should focus on scale-up and innovative process monitoring and control strategies.

The evolution of oxygen-utilizing enzymes suggests early biosphere oxygenation

Production of molecular oxygen was a turning point in the Earth's history. The geological record indicates the Great Oxidation Event, which marked a permanent transition to an oxidizing atmosphere around 2.4 Ga. However, the degree to which oxygen was available to life before oxygenation of the atmosphere remains unknown. Here, phylogenetic analysis of all known oxygen-utilizing and -producing enzymes (O₂-enzymes) indicates that oxygen became widely available to living organisms well before the Great Oxidation Event. About 60% of the O₂-enzyme families whose birth can be dated appear to have emerged at the separation of terrestrial and marine bacteria (22 families, compared to two families assigned to the last universal common ancestor). This node, dubbed the last universal oxygen ancestor, coincides with a burst of emergence of both oxygenases and other oxidoreductases, thus suggesting a wider availability of oxygen around 3.1 Ga.

Jumbo Phages: A Comparative Genomic Overview of Core Functions and Adaptions for Biological Conflicts

Jumbo phages have attracted much attention by virtue of their extraordinary genome size and unusual aspects of biology. By performing a comparative genomics analysis of 224 jumbo phages, we suggest an objective inclusion criterion based on genome size distributions and present a synthetic overview of their manifold adaptations across major biological systems. By means of clustering and principal component analysis of the phyletic patterns of conserved genes, all known jumbo phages can be classified into three higher-order groups, which include both myoviral and siphoviral morphologies indicating multiple independent origins from smaller predecessors. Our study uncovers several under-appreciated or unreported aspects of the DNA replication, recombination, transcription and virion maturation systems. Leveraging sensitive sequence analysis methods, we identify novel protein-modifying enzymes that might help hijack the host-machinery. Focusing on host–virus conflicts, we detect strategies used to counter different wings of the bacterial immune system, such as cyclic nucleotide- and NAD⁺-dependent effector-activation, and prevention of superinfection during pseudolysogeny. We reconstruct the RNA-repair systems of jumbo phages that counter the consequences of RNA-targeting host effectors. These findings also suggest that several jumbo phage proteins provide a snapshot of the systems found in ancient replicons preceding the last universal ancestor of cellular life.

Chloramphenicol biodegradation by enriched bacterial consortia and isolated strain Sphingomonas sp. CL5.1: The reconstruction of a novel biodegradation pathway

Figuring out the comprehensive metabolic mechanism of chloramphenicol (CAP) is critical to improving CAP removal in the bioremediation process. In this study, CAP biodegradation by six consortia and isolated Sphingomonas sp. CL5.1 were systematically investigated using the combination of high-performance liquid chromatography-quadrupole time-of-flight mass spectrometry, second-generation, and third-generation sequencing technologies. The CAP-degrading capability of six consortia was enhanced while CAP mineralization rate declined after long-term enrichment. The microbial community structures of six consortia were all simplified with 69%-82% decline in species richness after continuous passages for one year. The core genera of consortia CL and CH included Sphingomonas, Cupriavidus, Burkholderia, Chryseobacterium, and Pigmentiphaga, which accounted for over 98% of the total population. Sphingomonas was discovered as a new CAP degrader that could subsist on CAP as the sole carbon, nitrogen, and energy sources. Sphingomonas sp. CL5.1 was able to completely remove 120 mg/L CAP within 48 hours with a mineralization rate of 50.4%. The presence of acetate or nitrite could inhibit CAP metabolization by strain CL5.1. Four CAP metabolic pathways were constructed, including modification of the C₃ hydroxyl group of CAP via acetylation, oxidization, dehydration and the bond cleavage between C₁ and C₂. C₃ hydroxyl group dehydration and C₁-C₂ bond-cleavage were first reported regarding to CAP biotransformation. Strain CL5.1 played a core role in the consortia and was responsible for C3 hydroxyl oxidation, C3 dehydration, and C1-C2 bond cleavage. Genomic information of strain CL5.1 revealed the further mineralization pathways of downstream product p-nitrobenzoic acid via ortho- and meta-cleavage.

Finding MEMO-Emerging Evidence for MEMO1's Function in Development and Disease

Although conserved throughout animal kingdoms, the protein encoded by the gene Mediator of ERBB2 Driven Cell Motility 1 or MEMO1, has only recently come into focus. True to its namesake, MEMO1 first emerged from a proteomic screen of molecules bound to the ERBB2 receptor and was found to be necessary for efficient cell migration upon receptor activation. While initially placed within the context of breast cancer metastasis—a pathological state that has provided tremendous insight into MEMO1′s cellular roles—MEMO1′s function has since expanded to encompass additional cancer cell types, developmental processes during embryogenesis and homeostatic regulation of adult organ systems. Owing to MEMO1′s deep conservation, a variety of model organisms have been amenable to uncovering biological facets of this multipurpose protein; facets ranging from the cellular (e.g., receptor signaling, cytoskeletal regulation, redox flux) to the organismal (e.g., mineralization and mineral homeostasis, neuro/gliogenesis, vasculogenesis) level. Although these facets emerge at the intersection of numerous biological and human disease processes, how and if they are interconnected remains to be resolved. Here, we review our current understanding of this ‘enigmatic’ molecule, its role in development and disease and open questions emerging from these previous studies.

Evolution of l-DOPA 4,5-dioxygenase activity allows for recurrent specialisation to betalain pigmentation in Caryophyllales

The evolution of l-DOPA 4,5-dioxygenase activity, encoded by the gene DODA, was a key step in the origin of betalain biosynthesis in Caryophyllales. We previously proposed that l-DOPA 4,5-dioxygenase activity evolved via a single Caryophyllales-specific neofunctionalisation event within the DODA gene lineage. However, this neofunctionalisation event has not been confirmed and the DODA gene lineage exhibits numerous gene duplication events, whose evolutionary significance is unclear. To address this, we functionally characterised 23 distinct DODA proteins for l-DOPA 4,5-dioxygenase activity, from four betalain-pigmented and five anthocyanin-pigmented species, representing key evolutionary transitions across Caryophyllales. By mapping these functional data to an updated DODA phylogeny, we then explored the evolution of l-DOPA 4,5-dioxygenase activity. We find that low l-DOPA 4,5-dioxygenase activity is distributed across the DODA gene lineage. In this context, repeated gene duplication events within the DODA gene lineage give rise to polyphyletic occurrences of elevated l-DOPA 4,5-dioxygenase activity, accompanied by convergent shifts in key functional residues and distinct genomic patterns of micro-synteny. In the context of an updated organismal phylogeny and newly inferred pigment reconstructions, we argue that repeated convergent acquisition of elevated l-DOPA 4,5-dioxygenase activity is consistent with recurrent specialisation to betalain synthesis in Caryophyllales.

Highly regulated, diversifying NTP-dependent biological conflict systems with implications for the emergence of multicellularity

Social cellular aggregation or multicellular organization pose increased risk of transmission of infections through the system upon infection of a single cell. The generality of the evolutionary responses to this outside of Metazoa remains unclear. We report the discovery of several thematically unified, remarkable biological conflict systems preponderantly present in multicellular prokaryotes. These combine thresholding mechanisms utilizing NTPase chaperones (the MoxR-vWA couple), GTPases and proteolytic cascades with hypervariable effectors, which vary either by using a reverse transcriptase-dependent diversity-generating system or through a system of acquisition of diverse protein modules, typically in inactive form, from various cellular subsystems. Conciliant lines of evidence indicate their deployment against invasive entities, like viruses, to limit their spread in multicellular/social contexts via physical containment, dominant-negative interactions or apoptosis. These findings argue for both a similar operational 'grammar' and shared protein domains in the sensing and limiting of infections during the multiple emergences of multicellularity.

Structural and Biochemical Insights Into Two BAHD Acyltransferases (AtSHT and AtSDT) Involved in Phenolamide Biosynthesis

Phenolamides represent one of the largest classes of plant-specialized secondary metabolites and function in diverse physiological processes, including defense responses and development. The biosynthesis of phenolamides requires the BAHD-family acyltransferases, which transfer acyl-groups from different acyl-donors specifically to amines, the acyl-group acceptors. However, the mechanisms of substrate specificity and multisite-acylation of the BAHD-family acyltransferases remain poorly understood. In this study, we provide a structural and biochemical analysis of AtSHT and AtSDT, two representative BAHD-family members that catalyze the multisite acylation of spermidine but show different product profiles. By determining the structures of AtSHT and AtSDT and using structure-based mutagenesis, we identified the residues important for substrate recognition in AtSHT and AtSDT and hypothesized that the acyl acceptor spermidine might adopt a free-rotating conformation in AtSHT, which can undergo mono-, di-, or tri-acylation; while the spermidine molecule in AtSDT might adopt a linear conformation, which only allows mono- or di-acylation to take place. In addition, through sequence similarity network (SSN) and structural modeling analysis, we successfully predicted and verified the functions of two uncharacterized Arabidopsis BAHD acyltransferases, OAO95042.1 and NP_190301.2, which use putrescine as the main acyl-acceptor. Our work provides not only an excellent starting point for understanding multisite acylation in BAHD-family enzymes, but also a feasible methodology for predicting possible acyl acceptor specificity of uncharacterized BAHD-family acyltransferases.

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.