Functional assignment of multiple catabolic pathways for D-apiose

Expanding the Enzymatic Toolbox for Carboligation: Increasing the Diversity of the 'Split' Transketolase Sequence Space

Transketolases (TKs) are thiamine diphosphate (ThDP)-dependent enzymes that catalyze the transfer of two-carbon units in a stereoselective manner, making them valuable biocatalysts for sustainable processes. Most known TKs are about 650 amino acids long; however, a second type found in Archaea and many Bacteria consists of two proteins, each of about 300 amino acids. Exploring the unique features and differences of split TKs may help in assessing their potential use in biocatalysis and for uncovering new reactivities. Additionally, it could provide valuable information on how their structure relates to their function, especially compared to full-length TKs. In this study, we significantly expanded the known repertoire of split TKs approximately 14-fold to the best of our knowledge, by identifying and providing accessions of nearly 500 putative split-TK subunit pairs. Moreover, we doubled the number of experimentally produced and tested split TKs by cloning, purifying, and testing ten candidates retrieved from genomes and in-house metagenomes. Interestingly, pQR2809 and pQR2812, derived from hyperthermophilic organisms, showed enhanced thermostability compared to other TK examples in the literature, maintaining partial activity after heating at 90 °C or 100 °C for 1 hour, respectively.

Effects of dietary l-carnosine supplementation on the growth, intestinal microbiota, and serum metabolome of fattening lambs

Dietary l-carnosine supplementation has been shown to enhance animal performance and improve meat quality. However, the mechanisms underlying the effects of l-carnosine on the physiological functions of animals have not been fully elucidated. We investigated the effects of dietary l-carnosine supplementation on growth performance, intestinal microbiota diversity, and the serum metabolome in fattening lambs to reveal the molecular mechanism underlying the effect of l-carnosine on the growth performance of sheep. Sixty 3-month-old male crossbred lambs (Dorper ♂ × Small Tail Han ♀) with an average body weight of 30 ± 5 kg were randomly divided into two groups: a control group (group C) fed a basal diet, and an experimental group (group L) fed a basal diet supplemented with 400 mg/kg of l-carnosine. At the end of the 60-day experiment, all sheep were weighed, and fecal and blood samples were collected from 12 random sheep. The fecal microbiota was analyzed using 16S rRNA sequencing, and serum metabolites were analyzed using liquid chromatography-tandem mass spectrometry. Spearman correlation analysis was employed to assess the associations between intestinal microbiota and serum metabolite biomarkers. The results showed that weight gain and daily weight gain were significantly increased in group L compared to group C (p < 0.01). The dominant phyla in the intestinal microbiota (Firmicutes and Bacteroidetes) did not significantly differ between the two groups (p > 0.05). At the genus level, the abundances of Syntrophococcus (p < 0.01) and Butyricimonas (p < 0.001) were higher, whereas those of Escherichia-Shigella and Candidatus Saccharimonas were significantly lower in group L than in group C (p < 0.05). Non-targeted metabolomics identified 68 differentially abundant biomarkers (VIP > 1, p < 0.05). The content of pyridine N-oxide glucuronide was significantly downregulated (p < 0.01), whereas those of l-histidinol, d-apiose, and isodomedin were significantly upregulated in group L versus group C (p < 0.001). Holdemania and Butyricimonas were positively correlated with l-histidine, d-apiose, and l-erythrulose (p < 0.001), whereas Butyricimonas was negatively correlated with pyridine N-oxide glucuronide (p < 0.001). This study provided new insights into the effects of l-carnosine on the intestinal microbiota and nutrient metabolism in fattening sheep that will be helpful for the future application of l-carnosine in ruminants.

Accurate Identification of Periplasmic Urea-binding Proteins by Structure- and Genome Context-assisted Functional Analysis

ABC transporters are ancient and ubiquitous nutrient transport systems in bacteria and play a central role in defining lifestyles. Periplasmic solute-binding proteins (SBPs) are components that deliver ligands to their translocation machinery. SBPs have diversified to bind a wide range of ligands with high specificity and affinity. However, accurate assignment of cognate ligands remains a challenging problem in SBPs. Urea metabolism plays an important role in the nitrogen cycle; anthropogenic sources account for more than half of global nitrogen fertilizer. We report identification of urea-binding proteins within a large SBP sequence family that encodes diverse functions. By combining genetic linkage between SBPs, ABC transporter components, enzymes or transcription factors, we accurately identified cognate ligands, as we verified experimentally by biophysical characterization of ligand binding and crystallographic determination of the urea complex of a thermostable urea-binding homolog. Using three-dimensional structure information, these functional assignments were extrapolated to other members in the sequence family lacking genetic linkage information, which revealed that only a fraction bind urea. Using the same combined approaches, we also inferred that other family members bind various short-chain amides, aliphatic amino acids (leucine, isoleucine, valine), γ-aminobutyrate, and as yet unknown ligands. Comparative structural analysis revealed structural adaptations that encode diversification in these SBPs. Systematic assignment of ligands to SBP sequence families is key to understanding bacterial lifestyles, and also provides a rich source of biosensors for clinical and environmental analysis, such as the thermostable urea-binding protein identified here.

The ultra-high affinity transport proteins of ubiquitous marine bacteria

SAR11 bacteria are the most abundant microorganisms in the surface ocean1 and have global biogeochemical importance2-4. To thrive in their competitive oligotrophic environment, these bacteria rely heavily on solute-binding proteins that facilitate uptake of specific substrates via membrane transporters5,6. The functions and properties of these transport proteins are key factors in the assimilation of dissolved organic matter and biogeochemical cycling of nutrients in the ocean, but they have remained largely inaccessible to experimental investigation. Here we performed genome-wide experimental characterization of all solute-binding proteins in a prototypical SAR11 bacterium, revealing specific functions and general trends in their properties that contribute to the success of SAR11 bacteria in oligotrophic environments. We found that the solute-binding proteins of SAR11 bacteria have extremely high binding affinity (dissociation constant >20 pM) and high binding specificity, revealing molecular mechanisms of oligotrophic adaptation. Our functional data have uncovered new carbon sources for SAR11 bacteria and enable accurate biogeographical analysis of SAR11 substrate uptake capabilities throughout the ocean. This study provides a comprehensive view of the substrate uptake capabilities of ubiquitous marine bacteria, providing a necessary foundation for understanding their contribution to assimilation of dissolved organic matter in marine ecosystems.

Double-duty isomerases: a case study of isomerization-coupled enzymatic catalysis

Enzymes can usually be unambiguously assigned to one of seven classes specifying the basic chemistry of their catalyzed reactions. Less frequently, two or more reaction classes are catalyzed by a single enzyme within one active site. Two examples are an isomerohydrolase and an isomero-oxygenase that catalyze isomerization-coupled reactions crucial for production of vision-supporting 11-cis-retinoids. In these enzymes, isomerization is obligately paired and mechanistically intertwined with a second reaction class. A handful of other enzymes carrying out similarly coupled isomerization reactions have been described, some of which have been subjected to detailed structure-function analyses. Herein we review these rarefied enzymes, focusing on the mechanistic and structural basis of their reaction coupling with the goal of revealing catalytic commonalities.

In vitro characterization of alternative l-threonate and d-erythronate catabolic pathways

l-threonate is the metabolite of vitamin C, while d-erythronate is the metabolite of N-acetyl-d-glucosamine, the nutritional supplement for joint health. They are widely distributed in the environment and human biofluids. Nevertheless, the catabolisms of l-threonate and d-erythronate are sparsely reported. Here we explored the functional diversity of an acid sugar kinase family (Pfam families PF07005-PF17042), and discovered a novel 2-oxo-tetronate kinase. The conserved genome neighborhood of the 2-oxo-tetronate kinase encodes members of class-II fructose-bisphosphate aldolase family (F_bP_aldolase, PF01116) and a dehydrogenase family (PF03446-PF14833). Instructed by this analysis, we experimentally verified that these enzymes are capable of degrading l-threonate into dihydroxyacetone phosphate (DHAP) in Arthrobacter sp. ZBG10, Clostridium scindens ATCC 35704, and Pseudonocardia dioxanivorans ATCC 55486. Meanwhile, a convergent catabolic pathway for d-erythronate was characterized in P. dioxanivorans ATCC 55486. Moreover, the phylogenetic distribution analysis indicates that the biological range of the identified l-threonate and d-erythronate catabolic pathways appears to extend mostly to members of the Actinomycetota, Cyanobacteriota, Bacillota, Pseudomonadota, and Bacteroidota phyla.

Discovery of a non-canonical prototype long-chain monoacylglycerol lipase through a structure-based endogenous reaction intermediate complex

The identification and characterization of enzyme function is largely lacking behind the rapidly increasing availability of large numbers of sequences and associated high-resolution structures. This is often hampered by lack of knowledge on in vivo relevant substrates. Here, we present a case study of a high-resolution structure of an unusual orphan lipase in complex with an endogenous C18 monoacylglycerol ester reaction intermediate from the expression host, which is insoluble under aqueous conditions and thus not accessible for studies in solution. The data allowed its functional characterization as a prototypic long-chain monoacylglycerol lipase, which uses a minimal lid domain to position the substrate through a hydrophobic tunnel directly to the enzyme's active site. Knowledge about the molecular details of the substrate binding site allowed us to modulate the enzymatic activity by adjusting protein/substrate interactions, demonstrating the potential of our findings for future biotechnology applications.

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.

A Novel Pathway for Biosynthesis of the Herbicidal Phosphonate Natural Product Phosphonothrixin Is Widespread in Actinobacteria

Phosphonothrixin is an herbicidal phosphonate natural product with an unusual, branched carbon skeleton. Bioinformatic analyses of the ftx gene cluster, which is responsible for synthesis of the compound, suggest that early steps of the biosynthetic pathway, up to production of the intermediate 2,3-dihydroxypropylphosphonic acid (DHPPA) are identical to those of the unrelated phosphonate natural product valinophos. This conclusion was strongly supported by the observation of biosynthetic intermediates from the shared pathway in spent media from two phosphonothrixin producing strains. Biochemical characterization of ftx-encoded proteins confirmed these early steps, as well as subsequent steps involving the oxidation of DHPPA to 3-hydroxy-2-oxopropylphosphonate and its conversion to phosphonothrixin by the combined action of an unusual heterodimeric, thiamine-pyrophosphate (TPP)-dependent ketotransferase and a TPP-dependent acetolactate synthase. The frequent observation of ftx-like gene clusters within actinobacteria suggests that production of compounds related to phosphonothrixin is common within these bacteria. IMPORTANCE Phosphonic acid natural products, such as phosphonothrixin, have great potential for biomedical and agricultural applications; however, discovery and development of these compounds requires detailed knowledge of the metabolism involved in their biosynthesis. The studies reported here reveal the biochemical pathway phosphonothrixin production, which enhances our ability to design strains that overproduce this potentially useful herbicide. This knowledge also improves our ability to predict the products of related biosynthetic gene clusters and the functions of homologous enzymes.

Enzymatic Conversion of CO2: From Natural to Artificial Utilization

Enzymatic carbon dioxide fixation is one of the most important metabolic reactions as it allows the capture of inorganic carbon from the atmosphere and its conversion into organic biomass. However, due to the often unfavorable thermodynamics and the difficulties associated with the utilization of CO2, a gaseous substrate that is found in comparatively low concentrations in the atmosphere, such reactions remain challenging for biotechnological applications. Nature has tackled these problems by evolution of dedicated CO2-fixing enzymes, i.e., carboxylases, and embedding them in complex metabolic pathways. Biotechnology employs such carboxylating and decarboxylating enzymes for the carboxylation of aromatic and aliphatic substrates either by embedding them into more complex reaction cascades or by shifting the reaction equilibrium via reaction engineering. This review aims to provide an overview of natural CO2-fixing enzymes and their mechanistic similarities. We also discuss biocatalytic applications of carboxylases and decarboxylases for the synthesis of valuable products and provide a separate summary of strategies to improve the efficiency of such processes. We briefly summarize natural CO2 fixation pathways, provide a roadmap for the design and implementation of artificial carbon fixation pathways, and highlight examples of biocatalytic cascades involving carboxylases. Additionally, we suggest that biochemical utilization of reduced CO2 derivates, such as formate or methanol, represents a suitable alternative to direct use of CO2 and provide several examples. Our discussion closes with a techno-economic perspective on enzymatic CO2 fixation and its potential to reduce CO2 emissions.

Microbial life in 25-m-deep boreholes in ancient permafrost illuminated by metagenomics

This study describes the composition and potential metabolic adaptation of microbial communities in northeastern Siberia, a repository of the oldest permafrost in the Northern Hemisphere. Samples of contrasting depth (1.75 to 25.1 m below surface), age (from ~ 10 kyr to 1.1 Myr) and salinity (from low 0.1-0.2 ppt and brackish 0.3-1.3 ppt to saline 6.1 ppt) were collected from freshwater permafrost (FP) of borehole AL1_15 on the Alazeya River, and coastal brackish permafrost (BP) overlying marine permafrost (MP) of borehole CH1_17 on the East Siberian Sea coast. To avoid the limited view provided with culturing work, we used 16S rRNA gene sequencing to show that the biodiversity decreased dramatically with permafrost age. Nonmetric multidimensional scaling (NMDS) analysis placed the samples into three groups: FP and BP together (10-100 kyr old), MP (105-120 kyr old), and FP (> 900 kyr old). Younger FP/BP deposits were distinguished by the presence of Acidobacteriota, Bacteroidota, Chloroflexota_A, and Gemmatimonadota, older FP deposits had a higher proportion of Gammaproteobacteria, and older MP deposits had much more uncultured groups within Asgardarchaeota, Crenarchaeota, Chloroflexota, Patescibacteria, and unassigned archaea. The 60 recovered metagenome-assembled genomes and un-binned metagenomic assemblies suggested that despite the large taxonomic differences between samples, they all had a wide range of taxa capable of fermentation coupled to nitrate utilization, with the exception of sulfur reduction present only in old MP deposits.

Mechanism of Sugar Ring Contraction and Closure Catalyzed by UDP-d-apiose/UDP-d-xylose Synthase (UAXS)

Uridine diphosphate (UDP)-apiose/UDP-xylose synthase (UAXS) is a member of the short-chain dehydrogenase/reductase superfamily (SDR), which catalyzes the ring contraction and closure of UDP-d-glucuronic acid (UDP-GlcA), affording UDP-apiose and UDP-xylose. UAXS is a special enzyme that integrates ring-opening, decarboxylation, rearrangement, and ring closure/contraction in a single active site. Recently, the ternary complex structure of UAXS was crystallized from Arabidopsis thaliana. In this work, to gain insights into the detailed formation mechanism of UDP-apiose and UDP-xylose, an enzyme-substrate reactant model has been constructed and quantum mechanical/molecular mechanical (QM/MM) calculations have been performed. Our calculation results reveal that the reaction starts from the C4-OH oxidation, which is accompanied by the conformational transformation of the sugar ring from chair type to boat type. The sugar ring-opening is prior to decarboxylation, and the deprotonation of the C2-OH group is the prerequisite for sugar ring-opening. Moreover, the keto-enol tautomerization of the decarboxylated intermediate is a necessary step for ring closure/contraction. Based on our calculation results, more UDP-apiose product was expected, which is in line with the experimental observation. Three titratable residues, Tyr185, Cys100, and Cys140, steer the reaction by proton transfer from or to UDP-GlcA, and Arg182, Glu141, and D337 constitute a proton conduit for sugar C2-OH deprotonation. Although Thr139 and Tyr105 are not directly involved in the enzymatic reaction, they are responsible for promoting the catalysis by forming hydrogen-bonding interactions with GlcA. Our calculations may provide useful information for understanding the catalysis of the SDR family.

Comparative Metagenomics of the Active Layer and Permafrost from Low-Carbon Soil in the Canadian High Arctic

Approximately 87% of the Arctic consists of low-organic carbon mineral soil, but knowledge of microbial activity in low-carbon permafrost (PF) and active layer soils remains limited. This study investigated the taxonomic composition and genetic potential of microbial communities at contrasting depths of the active layer (5, 35, and 65 cm below surface, bls) and PF (80 cm bls). We showed microbial communities in PF to be taxonomically and functionally different from those in the active layer. 16S rRNA gene sequence analysis revealed higher biodiversity in the active layer than in PF, and biodiversity decreased significantly with depth. The reconstructed 91 metagenome-assembled genomes showed that PF was dominated by heterotrophic, fermenting Bacteroidota using nitrite as their main electron acceptor. Prevalent microbes identified in the active layer belonged to bacterial taxa, gaining energy via aerobic respiration. Gene abundance in metagenomes revealed enrichment of genes encoding the plant-derived polysaccharide degradation and metabolism of nitrate and sulfate in PF, whereas genes encoding methane/ammonia oxidation, cold-shock protein, and two-component systems were generally more abundant in the active layer, particularly at 5 cm bls. The results of this study deepen our understanding of the low-carbon Arctic soil microbiome and improve prediction of the impacts of thawing PF.

On the diversity of F420 -dependent oxidoreductases: A sequence- and structure-based classification

The F₄₂₀ deazaflavin cofactor is an intriguing molecule as it structurally resembles the canonical flavin cofactor, although behaves as a nicotinamide cofactor due to its obligate hydride-transfer reactivity and similar low redox potential. Since its discovery, numerous enzymes relying on it have been described. The known deazaflavoproteins are taxonomically restricted to Archaea and Bacteria. The biochemistry of the deazaflavoenzymes is diverse and they exhibit great structural variability. In this study a thorough sequence and structural homology evolutionary analysis was performed in order to generate an overarching classification of the F₄₂₀ -dependent oxidoreductases. Five different deazaflavoenzyme Classes (I-V) are described according to their structural folds as follows: Class I encompassing the TIM-barrel F₄₂₀ -dependent enzymes; Class II including the Rossmann fold F₄₂₀ -dependent enzymes; Class III comprising the β-roll F₄₂₀ -dependent enzymes; Class IV which exclusively gathers the SH3 barrel F₄₂₀ -dependent enzymes and Class V including the three layer ββα sandwich F₄₂₀ -dependent enzymes. This classification provides a framework for the identification and biochemical characterization of novel deazaflavoenzymes.

Deep ocean metagenomes provide insight into the metabolic architecture of bathypelagic microbial communities

The deep sea, the largest ocean's compartment, drives planetary-scale biogeochemical cycling. Yet, the functional exploration of its microbial communities lags far behind other environments. Here we analyze 58 metagenomes from tropical and subtropical deep oceans to generate the Malaspina Gene Database. Free-living or particle-attached lifestyles drive functional differences in bathypelagic prokaryotic communities, regardless of their biogeography. Ammonia and CO oxidation pathways are enriched in the free-living microbial communities and dissimilatory nitrate reduction to ammonium and H2 oxidation pathways in the particle-attached, while the Calvin Benson-Bassham cycle is the most prevalent inorganic carbon fixation pathway in both size fractions. Reconstruction of the Malaspina Deep Metagenome-Assembled Genomes reveals unique non-cyanobacterial diazotrophic bacteria and chemolithoautotrophic prokaryotes. The widespread potential to grow both autotrophically and heterotrophically suggests that mixotrophy is an ecologically relevant trait in the deep ocean. These results expand our understanding of the functional microbial structure and metabolic capabilities of the largest Earth aquatic ecosystem.

Discovery of new enzymatic functions and metabolic pathways using genomic enzymology web tools

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.

Discovery of novel pathways for carbohydrate metabolism

Closing the gap between the increasing availability of complete genome sequences and the discovery of novel enzymes in novel metabolic pathways is a significant challenge. Here, we review recent examples of assignment of in vitro enzymatic activities and in vivo metabolic functions to uncharacterized proteins, with a focus on enzymes and metabolic pathways involved in the catabolism and biosynthesis of monosaccharides and polysaccharides. The most effective approaches are based on analyses of sequence-function space in protein families that provide clues for the predictions of the functions of the uncharacterized enzymes. As summarized in this Opinion, this approach allows the discovery of the catabolism of new molecules, new pathways for common molecules, and new enzymatic chemistries.

Evolutionary trends in RuBisCO kinetics and their co-evolution with CO2 concentrating mechanisms

RuBisCO-catalyzed CO₂ fixation is the main source of organic carbon in the biosphere. This enzyme is present in all domains of life in different forms (III, II, and I) and its origin goes back to 3500 Mya, when the atmosphere was anoxygenic. However, the RuBisCO active site also catalyzes oxygenation of ribulose 1,5-bisphosphate, therefore, the development of oxygenic photosynthesis and the subsequent oxygen-rich atmosphere promoted the appearance of CO₂ concentrating mechanisms (CCMs) and/or the evolution of a more CO₂ -specific RuBisCO enzyme. The wide variability in RuBisCO kinetic traits of extant organisms reveals a history of adaptation to the prevailing CO₂ /O₂ concentrations and the thermal environment throughout evolution. Notable differences in the kinetic parameters are found among the different forms of RuBisCO, but the differences are also associated with the presence and type of CCMs within each form, indicative of co-evolution of RuBisCO and CCMs. Trade-offs between RuBisCO kinetic traits vary among the RuBisCO forms and also among phylogenetic groups within the same form. These results suggest that different biochemical and structural constraints have operated on each type of RuBisCO during evolution, probably reflecting different environmental selective pressures. In a similar way, variations in carbon isotopic fractionation of the enzyme point to significant differences in its relationship to the CO₂ specificity among different RuBisCO forms. A deeper knowledge of the natural variability of RuBisCO catalytic traits and the chemical mechanism of RuBisCO carboxylation and oxygenation reactions raises the possibility of finding unrevealed landscapes in RuBisCO evolution.

Deciphering the enzymatic mechanism of sugar ring contraction in UDP-apiose biosynthesis

D-Apiose is a C-branched pentose sugar important for plant cell wall development. Its biosynthesis as UDP-D-apiose involves decarboxylation of the UDP-D-glucuronic acid precursor coupled to pyranosyl-to-furanosyl sugar ring contraction. This unusual multistep reaction is catalyzed within a single active site by UDP-D-apiose/UDP-D-xylose synthase (UAXS). Here, we decipher the UAXS catalytic mechanism based on crystal structures of the enzyme from Arabidopsis thaliana, molecular dynamics simulations expanded by QM/MM calculations, and mutational-mechanistic analyses. Our studies show how UAXS uniquely integrates a classical catalytic cycle of oxidation and reduction by a tightly bound nicotinamide coenzyme with retro-aldol/aldol chemistry for the sugar ring contraction. They further demonstrate that decarboxylation occurs only after the sugar ring opening and identify the thiol group of Cys100 in steering the sugar skeleton rearrangement by proton transfer to and from the C3’. The mechanistic features of UAXS highlight the evolutionary expansion of the basic catalytic apparatus of short-chain dehydrogenases/reductases for functional versatility in sugar biosynthesis.

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.

Lateral Gene Transfer Shapes the Distribution of RuBisCO among Candidate Phyla Radiation Bacteria and DPANN Archaea

Ribulose-1,5-bisphosphate carboxylase/oxygenase (RuBisCO) is considered to be the most abundant enzyme on Earth. Despite this, its full diversity and distribution across the domains of life remain to be determined. Here, we leverage a large set of bacterial, archaeal, and viral genomes recovered from the environment to expand our understanding of existing RuBisCO diversity and the evolutionary processes responsible for its distribution. Specifically, we report a new type of RuBisCO present in Candidate Phyla Radiation (CPR) bacteria that is related to the archaeal Form III enzyme and contains the amino acid residues necessary for carboxylase activity. Genome-level metabolic analyses supported the inference that these RuBisCO function in a CO₂-incorporating pathway that consumes nucleotides. Importantly, some Gottesmanbacteria (CPR) also encode a phosphoribulokinase that may augment carbon metabolism through a partial Calvin–Benson–Bassham cycle. Based on the scattered distribution of RuBisCO and its discordant evolutionary history, we conclude that this enzyme has been extensively laterally transferred across the CPR bacteria and DPANN archaea. We also report RuBisCO-like proteins in phage genomes from diverse environments. These sequences cluster with proteins in the Beckwithbacteria (CPR), implicating phage as a possible mechanism of RuBisCO transfer. Finally, we synthesize our metabolic and evolutionary analyses to suggest that lateral gene transfer of RuBisCO may have facilitated major shifts in carbon metabolism in several important bacterial and archaeal lineages.

Back to the future: Why we need enzymology to build a synthetic metabolism of the future

Biology is turning from an analytical into a synthetic discipline. This is especially apparent in the field of metabolic engineering, where the concept of synthetic metabolism has been recently developed. Compared to classical metabolic engineering efforts, synthetic metabolism aims at creating novel metabolic networks in a rational fashion from bottom-up. However, while the theoretical design of synthetic metabolic networks has made tremendous progress, the actual realization of such synthetic pathways is still lacking behind. This is mostly because of our limitations in enzyme discovery and engineering to provide the parts required to build synthetic metabolism. Here I discuss the current challenges and limitations in synthetic metabolic engineering and elucidate how modern day enzymology can help to build a synthetic metabolism of the future.

Revisiting the methionine salvage pathway and its paralogues

Methionine is essential for life. Its chemistry makes it fragile in the presence of oxygen. Aerobic living organisms have selected a salvage pathway (the MSP) that uses dioxygen to regenerate methionine, associated to a ratchet-like step that prevents methionine back degradation. Here, we describe the variation on this theme, developed across the tree of life. Oxygen appeared long after life had developed on Earth. The canonical MSP evolved from ancestors that used both predecessors of ribulose bisphosphate carboxylase oxygenase (RuBisCO) and methanethiol in intermediate steps. We document how these likely promiscuous pathways were also used to metabolize the omnipresent by-products of S-adenosylmethionine radical enzymes as well as the aromatic and isoprene skeleton of quinone electron acceptors.

'Democratized' genomic enzymology web tools for functional assignment

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.