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

Unlocking a Sustainable Future for Plastics: A Chemical-Enzymatic Pathway for Efficient Conversion of Mixed Waste to MHET and Energy-Saving PET Recycling

The heterogeneous monomers obtained from plastic waste degradation are unfavorable for PET recondensation and high-value derivative synthesis. Herein, we developed an efficient chemical-enzymatic approach to convert mixed plastic wastes into homogeneous mono-2-hydroxyethyl terephthalate (MHET) without downstream purification, benefiting from three discovered BHETases (KbEst, KbHyd, and BrevEst) in nature. Towards the mixed plastic waste, integrating the chemical K2CO3-driven glycolysis process with the BHETase depolymerization technique resulted in an MHET yield of up to 98.26 % in 40 h. Remarkably, BrevEst accomplished the highest BHET hydrolysis (~87 % efficiency in 12 h) for yielding analytical-grade MHET compared to seven state-of-the-art PET hydrolases (18 %-40 %). In an investigation combining quantum theoretical computations and experimental validations, we established a MHET-initiated PET repolymerization pathway. This shortcut approach with MHET promises to strengthen the valorization of mixed plastics, offering a substantially more efficient and energy-saving route for PET recycling.

Limitations of Current Machine-Learning Models in Predicting Enzymatic Functions for Uncharacterized Proteins

Thirty to seventy percent of proteins in any given genome have no assigned function and have been labeled as the protein "unknownme". This large knowledge gap prevents the biological community from fully leveraging the plethora of genomic data that is now available. Machine-learning approaches are showing some promise in propagating functional knowledge from experimentally characterized proteins to the correct set of isofunctional orthologs. However, they largely fail to predict enzymatic functions unseen in the training set, as shown by dissecting the predictions made for 450 enzymes of unknown function from the model bacteria Escherichia coli using the DeepECTransformer platform. Lessons from these failures can help the community develop machine-learning methods that assist domain experts in making testable functional predictions for more members of the uncharacterized proteome.

Genome mining of sulfonated lanthipeptides reveals unique cyclic peptide sulfotransferases

Although sulfonation plays crucial roles in various biological processes and is frequently utilized in medicinal chemistry to improve water solubility and chemical diversity of drug leads, it is rare and underexplored in ribosomally synthesized and post-translationally modified peptides (RiPPs). Biosynthesis of RiPPs typically entails modification of hydrophilic residues, which substantially increases their chemical stability and bioactivity, albeit at the expense of reducing water solubility. To explore sulfonated RiPPs that may have improved solubility, we conducted co-occurrence analysis of RiPP class-defining enzymes and sulfotransferase (ST), and discovered two distinctive biosynthetic gene clusters (BGCs) encoding both lanthipeptide synthetase (LanM) and ST. Upon expressing these BGCs, we characterized the structures of novel sulfonated lanthipeptides and determined the catalytic details of LanM and ST. We demonstrate that SslST-catalyzed sulfonation is leader-independent but relies on the presence of A ring formed by LanM. Both LanM and ST are promiscuous towards residues in the A ring, but ST displays strict regioselectivity toward Tyr5. The recognition of cyclic peptide by ST was further discussed. Bioactivity evaluation underscores the significance of the ST-catalyzed sulfonation. This study sets up the starting point to engineering the novel lanthipeptide STs as biocatalysts for hydrophobic lanthipeptides improvement.

Biosynthesis of Macrocyclic Peptides with C-Terminal β-Amino-α-keto Acid Groups by Three Different Metalloenzymes

Advances in genome sequencing and bioinformatics methods have identified a myriad of biosynthetic gene clusters (BGCs) encoding uncharacterized molecules. By mining genomes for BGCs containing a prevalent peptide-binding domain used for the biosynthesis of ribosomally synthesized and post-translationally modified peptides (RiPPs), we uncovered a new compound class involving modifications installed by a cytochrome P450, a multinuclear iron-dependent non-heme oxidative enzyme (MNIO, formerly DUF692), a cobalamin- and radical S-adenosyl-l-methionine-dependent enzyme (B12-rSAM), and a methyltransferase. All enzymes were functionally expressed in Burkholderia sp. FERM BP-3421. Structural characterization demonstrated that the P450 enzyme catalyzed the formation of a biaryl C-C cross-link between two Tyr residues with the B12-rSAM generating β-methyltyrosine. The MNIO transformed a C-terminal Asp residue into aminopyruvic acid, while the methyltransferase acted on the β-carbon of this α-keto acid. Exciton-coupled circular dichroism spectroscopy and microcrystal electron diffraction (MicroED) were used to elucidate the stereochemical configuration of the atropisomer formed upon biaryl cross-linking. To the best of our knowledge, the MNIO featured in this pathway is the first to modify a residue other than Cys. This study underscores the utility of genome mining to isolate new macrocyclic RiPPs biosynthesized via previously undiscovered enzyme chemistry.

Genome Mining for New Enzyme Chemistry

A revolution in the field of biocatalysis has enabled scalable access to compounds of high societal values using enzymes. The construction of biocatalytic routes relies on the reservoir of available enzymatic transformations. A review of uncharacterized proteins predicted from genomic sequencing projects shows that a treasure trove of enzyme chemistry awaits to be uncovered. This Review highlights enzymatic transformations discovered through various genome mining methods and showcases their potential future applications in biocatalysis.

Biosynthesis and function of 7-deazaguanine derivatives in bacteria and phages

SUMMARYDeazaguanine 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.

Divergent Biosynthesis of Bridged Polycyclic Sesquiterpenoids by a Minimal Fungal Biosynthetic Gene Cluster

The bridged polycyclic sesquiterpenoids derived from sativene, isosativene, and longifolene have unique structures, and many chemical synthesis approaches with at least 10 steps have been reported. However, their biosynthetic pathway remains undescribed. A minimal biosynthetic gene cluster (BGC), named bip, encoding a sesquiterpene cyclase (BipA) and a cytochrome P450 (BipB) is characterized to produce such complex sesquiterpenoids with multiple carbon skeletons based on enzymatic assays, heterologous expression, and precursor experiments. BipA is demonstrated as a versatile cyclase with (-)-sativene as the dominant product and (-)-isosativene and (-)-longifolene as minor ones. BipB is capable of hydroxylating different enantiomeric sesquiterpenes, such as (-)-longifolene and (+)-longifolene, at C-15 and C-14 in turn. The C-15- or both C-15- and C-14-hydroxylated products are then further oxidized by unclustered oxidases, resulting in a structurally diverse array of sesquiterpenoids. Bioinformatic analysis reveals the BipB homologues as a discrete clade of fungal sesquiterpene P450s. These findings elucidate the concise and divergent biosynthesis of such intricate bridged polycyclic sesquiterpenoids, offer valuable biocatalysts for biotransformation, and highlight the distinct biosynthetic strategy employed by nature compared to chemical synthesis.

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.

Biosynthesis of macrocyclic peptides with C-terminal β-amino-α-keto acid groups by three different metalloenzymes

Advances in genome sequencing and bioinformatics methods have identified a myriad of biosynthetic gene clusters (BGCs) encoding uncharacterized molecules. By mining genomes for BGCs containing a prevalent peptide-binding domain used for the biosynthesis of ribosomally synthesized and post-translationally modified peptides (RiPPs), we uncovered a new class involving modifications installed by a cytochrome P450, a multi-nuclear iron-dependent non-heme oxidative enzyme (MNIO, formerly DUF692), a cobalamin- and radical S-adenosyl-L-methionine-dependent enzyme (B12-rSAM), and a methyltransferase. All enzymes encoded by the BGC were functionally expressed in Burkholderia sp. FERM BP-3421. Structural characterization with 2D-NMR and Marfey's method on the resulting RiPP demonstrated that the P450 enzyme catalyzed the formation of a biaryl C-C crosslink between two Tyr residues with the B12-rSAM generating β-methyltyrosine. The MNIO transformed a C-terminal Asp residue into aminopyruvic acid while the methyltransferase acted on the β-carbon of the α-keto acid. Exciton-coupled circular dichroism spectroscopy and microcrystal electron diffraction (MicroED) were used to elucidate the stereochemical configurations of the atropisomer that formed upon biaryl crosslinking. The conserved Cys residue in the precursor peptide was not modified as in all other characterized MNIO-containing BGCs; However, mutational analyses demonstrated that it was essential for the MNIO activity on the C-terminal Asp. To the best of our knowledge, the MNIO featured in this pathway is the first to modify a residue other than Cys. This study underscores the utility of genome mining to discover new macrocyclic RiPPs and that RiPPs remain a significant source of previously undiscovered enzyme chemistry.

Core-dependent post-translational modifications guide the biosynthesis of a new class of hypermodified peptides

The ribosomally synthesized and post-translationally modified peptide (RiPPs) class of natural products has undergone significant expansion due to the rapid growth in genome sequencing data. Using a bioinformatics approach, we identify the dehydrazoles, a novel class of hypermodified RiPPs that contain both side chain dehydration of Ser residues, and backbone heterocyclization at Ser, Thr, and Cys residues to the corresponding azol(in)es. Structure elucidation of the hypermodified peptide carnazolamide, a representative class member, shows that 18 post-translational modifications are installed by just five enzymes. Complete biosynthetic reconstitution demonstrates that dehydration is carried out by an unusual DUF4135 dehydration domain fused to a zinc-independent cyclase domain (CcaM). We demonstrate that CcaM only modifies Ser residues that precede an azole in the core peptide. As heterocyclization removes the carbonyl following the Ser residue, CcaM likely catalyzes dehydration without generating an enolate intermediate. Additionally, CcaM does not require the leader peptide, and this core-dependence effectively sets the order for the biosynthetic reactions. Biophysical studies demonstrate direct binding of azoles to CcaM consistent with this azole moiety-dependent dehydration. Bioinformatic analysis reveals more than 50 related biosynthetic gene clusters that contain additional catalysts that may produce structurally diverse scaffolds.

Genome Mining of Linaridins Provides Insights into the Widely Distributed LinC Oxidoreductases

Linaridins are a family of underexplored ribosomally synthesized and post-translationally modified peptides despite the prevalence of their biosynthetic gene clusters (BGCs) in microbial genomes, as shown by bioinformatic studies. Our genome mining efforts reveal that 96 putative oxidoreductase genes, namely, LinC, are encoded in linaridin BGCs. We heterologously expressed two such LinC-containing linaridin BGCs, yan and ydn, from Streptomyces yunnanensis and obtained three new linaridins, named yunnanaridins A-C (1-3). Their structures are characterized by Z-configurations of the dehydrobutyrines and the presence of a variety of epimerized amino acid residues. Yunnanaridin A (1) is the sixth member of the family of type-B linaridins, whereas yunnanaridins B (2) and C (3) represent the first examples of expressed type-C linaridins. Interestingly, heterologous expression of the same BGCs with LinC in-frame knockouts produced the same compounds. This work expands the structural diversity of linaridins and provides evidence for the notion that the widespread LinCs may not be involved in linaridin biosynthesis.

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.

Discovery and mechanism-guided engineering of BHET hydrolases for improved PET recycling and upcycling

Although considerable research achievements have been made to address the plastic crisis using enzymes, their applications are limited due to incomplete degradation and low efficiency. Herein, we report the identification and subsequent engineering of BHETases, which have the potential to improve the efficiency of PET recycling and upcycling. Two BHETases (ChryBHETase and BsEst) are identified from the environment via enzyme mining. Subsequently, mechanism-guided barrier engineering is employed to yield two robust and thermostable ΔBHETases with up to 3.5-fold enhanced kcat/KM than wild-type, followed by atomic resolution understanding. Coupling ΔBHETase into a two-enzyme system overcomes the challenge of heterogeneous product formation and results in up to 7.0-fold improved TPA production than seven state-of-the-art PET hydrolases, under the conditions used here. Finally, we employ a ΔBHETase-joined tandem chemical-enzymatic approach to valorize 21 commercial post-consumed plastics into virgin PET and an example chemical (p-phthaloyl chloride) for achieving the closed-loop PET recycling and open-loop PET upcycling.

Division of developmental phases of freshwater leech Whitmania pigra and key genes related to neurogenesis revealed by whole genome and transcriptome analysis

The freshwater leech Whitmania pigra (W. pigra) Whitman (Annelida phylum) is a model organism for neurodevelopmental studies. However, molecular biology research on its embryonic development is still scarce. Here, we described a series of developmental stages of the W. pigra embryos and defined five broad stages of embryogenesis: cleavage stages, blastocyst stage, gastrula stage, organogenesis and refinement, juvenile. We obtained a total of 239.64 Gb transcriptome data of eight representative developmental phases of embryos (from blastocyst stage to maturity), which was then assembled into 21,482 unigenes according to our reference genome sequenced by single-molecule real-time (SMRT) long-read sequencing. We found 3114 genes differentially expressed during the eight phases with phase-specific expression pattern. Using a comprehensive transcriptome dataset, we demonstrated that 57, 49 and 77 DEGs were respectively related to morphogenesis, signal pathways and neurogenesis. 49 DEGs related to signal pathways included 30 wnt genes, 14 notch genes, and 5 hedgehog genes. In particular, we found a cluster consisting of 7 genes related to signal pathways as well as synapses, which were essential for regulating embryonic development. Eight genes cooperatively participated in regulating neurogenesis. Our results reveal the whole picture of W. pigra development mechanism from the perspective of transcriptome and provide new clues for organogenesis and neurodevelopmental studies of Annelida species.

An informatic workflow for the enhanced annotation of excretory/secretory proteins of Haemonchus contortus

Major advances in genomic and associated technologies have demanded reliable bioinformatic tools and workflows for the annotation of genes and their products via comparative analyses using well-curated reference data sets, accessible in public repositories. However, the accurate in silico annotation of molecules (proteins) encoded in organisms (e.g., multicellular parasites) which are evolutionarily distant from those for which these extensive reference data sets are available, including invertebrate model organisms (e.g., Caenorhabditis elegans - free-living nematode, and Drosophila melanogaster - the vinegar fly) and vertebrate species (e.g., Homo sapiens and Mus musculus), remains a major challenge. Here, we constructed an informatic workflow for the enhanced annotation of biologically-important, excretory/secretory (ES) proteins ("secretome") encoded in the genome of a parasitic roundworm, called Haemonchus contortus (commonly known as the barber's pole worm). We critically evaluated the performance of five distinct methods, refined some of them, and then combined the use of all five methods to comprehensively annotate ES proteins, according to gene ontology, biological pathways and/or metabolic (enzymatic) processes. Then, using optimised parameter settings, we applied this workflow to comprehensively annotate 2591 of all 3353 proteins (77.3%) in the secretome of H. contortus. This result is a substantial improvement (10-25%) over previous annotations using individual, "off-the-shelf" algorithms and default settings, indicating the ready applicability of the present, refined workflow to gene/protein sequence data sets from a wide range of organisms in the Tree-of-Life.

Expanding the viewpoint: Leveraging sequence information in enzymology

The use of protein sequence to inform enzymology in terms of structure, mechanism, and function has burgeoned over the past two decades. Referred to as genomic enzymology, the utilization of bioinformatic tools such as sequence similarity networks and phylogenetic analyses has allowed the identification of new substrates and metabolites, novel pathways, and unexpected reaction mechanisms. The holistic examination of superfamilies can yield insight into the origins and paths of evolution of enzymes and the range of their substrates and mechanisms. Herein, we highlight advances in the use of genomic enzymology to address problems which the in-depth analyses of a single enzyme alone could not enable.

Amino acid (acyl carrier protein) ligase-associated biosynthetic gene clusters reveal unexplored biosynthetic potential

This study aimed to evaluate the postulated cellular function of a novel family of amino acid (acyl carrier protein) ligases (AALs) in natural product biosynthesis. Here, we analyzed the manually curated, putative, aal-associated natural product biosynthetic gene clusters (NP BGCs) using two computational platforms for NP prediction, antiSMASH-BiG-SCAPE-CORASON and DeepBGC. The detected BGCs included a diversity of type I polyketide/nonribosomal peptide (PKS/NRPS) hybrid BGCs, exemplified by the guadinomine BGC, which suggested a dedicated function of AALs in the biosynthesis of rare (2S)-aminomalonyl-ACP extension units. Besides modular PKS/NRPSs and NRPSs, AAL-associated BGCs were predicted to assemble arylpolyenes, ladderane lipids, phosphonates, aminoglycosides, β-lactones, and thioamides of both nonribosomal and ribosomal origins. Additionally, we revealed a frequent association of AALs with putative, seldom observed transglutaminase-like and BtrH-like transferases of the cysteine protease superfamily, which may form larger families of ACP-dependent amide bond catalysts used in NP synthesis. Our results disclosed an exceptional chemical novelty and biosynthetic potential of the AAL-associated BGCs in NP biosynthesis. The presented in silico evidence supports the initial hypothesis and provides an important foundation for future experimental studies aimed at discovering novel pharmaceutically relevant active compounds.

YcaO-mediated ATP-dependent peptidase activity in ribosomal peptide biosynthesis

YcaO enzymes catalyze ATP-dependent post-translation modifications on peptides, including the installation of (ox/thi)azoline, thioamide and/or amidine moieties. Here we demonstrate that, in the biosynthesis of the bis-methyloxazolic alkaloid muscoride A, the YcaO enzyme MusD carries out both ATP-dependent cyclodehydration and peptide bond cleavage, which is a mechanism unprecedented for such a reaction. YcaO-catalyzed modifications are proposed to occur through a backbone O-phosphorylated intermediate, but this mechanism remains speculative. We report, to our knowedge, the first characterization of an acyl-phosphate species consistent with the proposed mechanism for backbone amide activation. The 3.1-Å-resolution cryogenic electron microscopy structure of MusD along with biochemical analysis allow identification of residues that enable peptide cleavage reaction. Bioinformatics analysis identifies other cyanobactin pathways that may deploy bifunctional YcaO enzymes. Our structural, mutational and mechanistic studies expand the scope of modifications catalyzed by YcaO proteins to include peptide hydrolysis and provide evidence for a unifying mechanism for the catalytically diverse outcomes.

General strategies for using amino acid sequence data to guide biochemical investigation of protein function

The rapid increase of '-omics' data warrants the reconsideration of experimental strategies to investigate general protein function. Studying individual members of a protein family is likely insufficient to provide a complete mechanistic understanding of family functions, especially for diverse families with thousands of known members. Strategies that exploit large amounts of available amino acid sequence data can inspire and guide biochemical experiments, generating broadly applicable insights into a given family. Here we review several methods that utilize abundant sequence data to focus experimental efforts and identify features truly representative of a protein family or domain. First, coevolutionary relationships between residues within primary sequences can be successfully exploited to identify structurally and/or functionally important positions for experimental investigation. Second, functionally important variable residue positions typically occupy a limited sequence space, a property useful for guiding biochemical characterization of the effects of the most physiologically and evolutionarily relevant amino acids. Third, amino acid sequence variation within domains shared between different protein families can be used to sort a particular domain into multiple subtypes, inspiring further experimental designs. Although generally applicable to any kind of protein domain because they depend solely on amino acid sequences, the second and third approaches are reviewed in detail because they appear to have been used infrequently and offer immediate opportunities for new advances. Finally, we speculate that future technologies capable of analyzing and manipulating conserved and variable aspects of the three-dimensional structures of a protein family could lead to broad insights not attainable by current methods.

Flavonoids' Dual Benefits in Gastrointestinal Cancer and Diabetes: A Potential Treatment on the Horizon?

Diabetes and gastrointestinal cancers (GI) are global health conditions with a massive burden on patients' lives worldwide. The development of both conditions is influenced by several factors, such as diet, genetics, environment, and infection, which shows a potential link between them. Flavonoids are naturally occurring phenolic compounds present in fruits and vegetables. Once ingested, unabsorbed flavonoids reaching the colon undergo enzymatic modification by the gut microbiome to facilitate absorption and produce ring fission products. The metabolized flavonoids exert antidiabetic and anti-GI cancer properties, targeting major impaired pathways such as apoptosis and cellular proliferation in both conditions, suggesting the potentially dual effects of flavonoids on diabetes and GI cancers. This review summarizes the current knowledge on the impact of flavonoids on diabetes and GI cancers in four significant pathways. It also addresses the synergistic effects of selected flavonoids on both conditions. While this is an intriguing approach, more studies are required to better understand the mechanism of how flavonoids can influence the same impaired pathways with different outcomes depending on the disease.

The Power of Biocatalysts for Highly Selective and Efficient Phosphorylation Reactions

Reactions involving the transfer of phosphorus-containing groups are of key importance for maintaining life, from biological cells, tissues and organs to plants, animals, humans, ecosystems and the whole planet earth. The sustainable utilization of the nonrenewable element phosphorus is of key importance for a balanced phosphorus cycle. Significant advances have been achieved in highly selective and efficient biocatalytic phosphorylation reactions, fundamental and applied aspects of phosphorylation biocatalysts, novel phosphorylation biocatalysts, discovery methodologies and tools, analytical and synthetic applications, useful phosphoryl donors and systems for their regeneration, reaction engineering, product recovery and purification. Biocatalytic phosphorylation reactions with complete conversion therefore provide an excellent reaction platform for valuable analytical and synthetic applications.

Class V Lanthipeptide Cyclase Directs the Biosynthesis of a Stapled Peptide Natural Product

Lanthipeptides are a class of cyclic peptides characterized by the presence of one or more lanthionine (Lan) or methyllanthionine (MeLan) thioether rings. These cross-links are produced by α,β-unsaturation of Ser or Thr residues in peptide substrates by dehydration, followed by a Michael-type conjugate addition of Cys residues onto the dehydroamino acids. Lanthipeptides may be broadly classified into at least five different classes, and the biosynthesis of classes I-IV lanthipeptides requires catalysis by LanC cyclases that control both the site-specificity and the stereochemistry of the conjugate addition. In contrast, there are no current examples of LanCs that occur in class V biosynthetic clusters, despite the presence of lanthionine rings in these compounds. In this work, bioinformatics-guided co-occurrence analysis identifies more than 240 putative class V lanthipeptide clusters that contain a LanC cyclase. Reconstitution studies demonstrate that the cyclase-catalyzed product is notably distinct from the product formed spontaneously. Stereochemical analysis shows that the cyclase diverts the final product to a configuration that is distinct from one that is energetically favored. Structural characterization of the final product by multi-dimensional NMR spectroscopy reveals that it forms a helical stapled peptide. Mutational analysis identified a plausible order for cyclization and suggests that enzymatic rerouting to the final structure is largely directed by the construction of the first lanthionine ring. These studies show that lanthipeptide cyclases are needed for the biosynthesis of some constrained peptides, the formations of which would otherwise be energetically unfavored.

RaS-RiPPs in Streptococci and the Human Microbiome

Radical S-adenosylmethionine (RaS) enzymes have quickly advanced to one of the most abundant and versatile enzyme superfamilies known. Their chemistry is predicated upon reductive homolytic cleavage of a carbon–sulfur bond in cofactor S-adenosylmethionine forming an oxidizing carbon-based radical, which can initiate myriad radical transformations. An emerging role for RaS enzymes is their involvement in the biosynthesis of ribosomally synthesized and post-translationally modified peptides (RiPPs), a natural product family that has become known as RaS-RiPPs. These metabolites are especially prevalent in human and mammalian microbiomes because the complex chemistry of RaS enzymes gives rise to correspondingly complex natural products with minimal cellular energy and genomic fingerprint, a feature that is advantageous in microbes with small, host-adapted genomes in competitive environments. Herein, we review the discovery and characterization of RaS-RiPPs from the human microbiome with a focus on streptococcal bacteria. We discuss the varied chemical modifications that RaS enzymes introduce onto their peptide substrates and the diverse natural products that they give rise to. The majority of RaS-RiPPs remain to be discovered, providing an intriguing avenue for future investigations at the intersection of metalloenzymology, chemical ecology, and the human microbiome.

A roadmap for the functional annotation of protein families: a community perspective

Over the last 25 years, biology has entered the genomic era and is becoming a science of ‘big data’. Most interpretations of genomic analyses rely on accurate functional annotations of the proteins encoded by more than 500 000 genomes sequenced to date. By different estimates, only half the predicted sequenced proteins carry an accurate functional annotation, and this percentage varies drastically between different organismal lineages. Such a large gap in knowledge hampers all aspects of biological enterprise and, thereby, is standing in the way of genomic biology reaching its full potential. A brainstorming meeting to address this issue funded by the National Science Foundation was held during 3–4 February 2022. Bringing together data scientists, biocurators, computational biologists and experimentalists within the same venue allowed for a comprehensive assessment of the current state of functional annotations of protein families. Further, major issues that were obstructing the field were identified and discussed, which ultimately allowed for the proposal of solutions on how to move forward.

A Silent Biosynthetic Gene Cluster from a Methanotrophic Bacterium Potentiates Discovery of a Substrate Promiscuous Proteusin Cyclodehydratase

Natural product-encoding biosynthetic gene clusters (BGCs) within microbial genomes far outnumber the known natural products; chemical products from such BGCs remain cryptic. These silent BGCs hold promise not only for the elaboration of new natural products but also for the discovery of useful biosynthetic enzymes. Here, we describe a genome mining strategy targeted toward the discovery of substrate promiscuous natural product biosynthetic enzymes. In the genome of the methanotrophic bacterium Methylovulum psychrotolerans Sph1^T, we discover a transcriptionally silent natural product BGC that encoded numerous ribosomally synthesized and post-translationally modified peptide (RiPP) natural products. These cryptic RiPP natural products were accessed using heterologous expression of the substrate peptide and biosynthetic enzyme-encoded genes. In line with our genome mining strategy, the RiPP biosynthetic enzymes in this BGC were found to be substrate promiscuous, which allowed us to use them in a combinatorial fashion with a similarly substrate-tolerant cyanobactin biosynthetic enzyme to introduce head-to-tail macrocyclization in the proteusin family of RiPP natural products.

Catalysts for the Enzymatic Lipidation of Peptides

Biologically active peptides are a major growing class of drugs, but their therapeutic potential is constrained by several limitations including bioavailability and poor pharmacokinetics. The attachment of functional groups like lipids has proven to be a robust and effective strategy for improving their therapeutic potential. Biochemical and bioactivity-guided screening efforts have identified the cyanobactins as a large class of ribosomally synthesized and post-translationally modified peptides (RiPPs) that are modified with lipids. These lipids are attached by the F superfamily of peptide prenyltransferase enzymes that utilize 5-carbon (prenylation) or 10-carbon (geranylation) donors. The chemical structures of various cyanobactins initially showed isoprenoid attachments on Ser, Thr, or Tyr. Biochemical characterization of the F prenyltransferases from the corresponding clusters shows that the different enzymes have different acceptor residue specificities but are otherwise remarkably sequence tolerant. Hence, these enzymes are well suited for biotechnological applications. The crystal structure of the Tyr O-prenyltransferase PagF reveals that the F enzyme shares a domain architecture reminiscent of a canonical ABBA prenyltransferase fold but lacks secondary structural elements necessary to form an enclosed active site. Binding of either cyclic or linear peptides is sufficient to close the active site to allow for productive catalysis, explaining why these enzymes cannot use isolated amino acids as substrates.Almost all characterized isoprenylated cyanobactins are modified with 5-carbon isoprenoids. However, chemical characterization demonstrates that the piricyclamides are modified with a 10-carbon geranyl moiety, and in vitro reconstitution of the corresponding PirF shows that the enzyme is a geranyltransferase. Structural analysis of PirF shows an active site nearly identical with that of the PagF prenyltransferase but with a single amino acid substitution. Of note, mutation at this residue in PagF or PirF can completely switch the isoprenoid donor specificity of these enzymes. Recent efforts have resulted in significant expansion of the F family with enzymes identified that can carry out C-prenylations of Trp, N-prenylations of Trp, and bis-N-prenylations of Arg. Additional genome-guided efforts based on the sequence of F enzymes identify linear cyanobactins that are α-N-prenylated and α-C-methylated by a bifunctional prenyltransferase/methyltransferase fusion and a bis-α-N- and α-C-prenylated linear peptide. The discovery of these different classes of prenyltransferases with diverse acceptor residue specificities expands the biosynthetic toolkit for enzymatic prenylation of peptide substrates.In this Account, we review the current knowledge scope of the F family of peptide prenyltransferases, focusing on the biochemical, structure-function, and chemical characterization studies that have been carried out in our laboratories. These enzymes are easily amenable for diversity-oriented synthetic efforts as they can accommodate substrate peptides of diverse sequences and are thus attractive catalysts for use in synthetic biology approaches to generate high-value peptidic therapeutics.

Fungal dye-decolorizing peroxidase diversity: roles in either intra- or extracellular processes

Abstract

Fungal dye-decolorizing peroxidases (DyPs) have found applications in the treatment of dye-contaminated industrial wastes or to improve biomass digestibility. Their roles in fungal biology are uncertain, although it has been repeatedly suggested that they could participate in lignin degradation and/or modification. Using a comprehensive set of 162 fully sequenced fungal species, we defined seven distinct fungal DyP clades on basis of a sequence similarity network. Sequences from one of these clades clearly diverged from all others, having on average the lower isoelectric points and hydropathy indices, the highest number of N-glycosylation sites, and N-terminal sequence peptides for secretion. Putative proteins from this clade are absent from brown-rot and ectomycorrhizal species that have lost the capability of degrading lignin enzymatically. They are almost exclusively present in white-rot and other saprotrophic Basidiomycota that digest lignin enzymatically, thus lending support for a specific role of DyPs from this clade in biochemical lignin modification. Additional nearly full-length fungal DyP genes were isolated from the environment by sequence capture by hybridization; they all belonged to the clade of the presumably secreted DyPs and to another related clade. We suggest focusing our attention on the presumably intracellular DyPs from the other clades, which have not been characterized thus far and could represent enzyme proteins with novel catalytic properties.

Key points

• A fungal DyP phylogeny delineates seven main sequence clades.

• Putative extracellular DyPs form a single clade of Basidiomycota sequences.

• Extracellular DyPs are associated to white-rot fungi.

Supplementary Information

The online version contains supplementary material available at 10.1007/s00253-022-11923-0.

RadicalSAM.org: A Resource to Interpret Sequence-Function Space and Discover New Radical SAM Enzyme Chemistry

The radical SAM superfamily (RSS), arguably the most functionally diverse enzyme superfamily, is also one of the largest with ~700K members currently in the UniProt database. The vast majority of the members have uncharacterized enzymatic activities and metabolic functions. In this Perspective, we describe RadicalSAM.org, a new web-based resource that enables a user-friendly genomic enzymology strategy to explore sequence-function space in the RSS. The resource attempts to enable identification of isofunctional groups of radical SAM enzymes using sequence similarity networks (SSNs) and the genome context of the bacterial, archaeal, and fungal members provided by genome neighborhood diagrams (GNDs). Enzymatic activities and in vivo functions frequently can be inferred from genome context given the tendency for genes of related function to be clustered. We invite the scientific community to use RadicalSAM.org to (i) guide their experimental studies to discover new enzymatic activities and metabolic functions, (ii) contribute experimentally verified annotations to RadicalSAM.org to enhance the ability to predict novel activities and functions, and (iii) provide suggestions for improving this resource.

Biodegradation of herbicides by a plant nonheme iron dioxygenase: mechanism and selectivity of substrate analogues

Aryloxyalkanoate dioxygenases are unique herbicide biodegrading nonheme iron enzymes found in plants and hence, from environmental and agricultural point of view they are important and valuable. However, they often are substrate specific and little is known on the details of the mechanism and the substrate scope. To this end, we created enzyme models and calculate the mechanism for 2,4-dichlorophenoxyacetic acid biodegradation and 2-methyl substituted analogs by density functional theory. The work shows that the substrate binding is tight and positions the aliphatic group close to the metal center to enable a chemoselective reaction mechanism to form the C 2 -hydroxy products, whereas the aromatic hydroxylation barriers are well higher in energy. Subsequently, we investigated the metabolism of R - and S -methyl substituted inhibitors and show that these do not react as efficiently as 2,4-dichlorophenoxyacetic acid substrate due to stereochemical clashes in the active site and particularly for the R -isomer give high rebound barriers.

The Birth of Genomic Enzymology: Discovery of the Mechanistically Diverse Enolase Superfamily

Enzymes are categorized into superfamilies by sequence, structural, and mechanistic similarities. The evolutionary implications can be profound. Until the mid-1990s, the approach was fragmented largely due to limited sequence and structural data. However, in 1996, Babbitt et al. published a paper in Biochemistry that demonstrated the potential power of mechanistically diverse superfamilies to identify common ancestry, predict function, and, in some cases, predict specificity. This Perspective describes the findings of the original work and reviews the current understanding of structure and mechanism in the founding family members. The outcomes of the genomic enzymology approach have reached far beyond the functional assignment of members of the enolase superfamily, inspiring the study of superfamilies and the adoption of sequence similarity networks and genome context and yielding fundamental insights into enzyme evolution.

Rid7C, a member of the YjgF/YER057c/UK114 (Rid) protein family, is a novel endoribonuclease that regulates the expression of a specialist RNA polymerase involved in differentiation in Nonomuraea gerenzanensis

The YjgF/YER057c/UK114 (Rid) is a protein family breadth conserved in all domains of life and includes the widely distributed archetypal RidA (YjgF) subfamily and seven other subfamilies (Rid1 to Rid7). Among these subfamilies, RidA is the only family to have been biochemically well characterized and is involved in the deamination of the reactive enamine/imine intermediates. In this study, we have characterized a protein of the Rid7 subfamily, named Rid7C, in Nonomuraea gerenzanensis, an actinomycete that is characterized by the presence of two types of RNA polymerases. This is due to the co-existence in its genome of two RNAP β chain-encoding genes: rpoB(S) (the wild-type rpoB gene) and rpoB(R) (a specialist, mutant-type rpoB gene) that controls A40926 antibiotic production and a wide range of metabolic adaptive behaviors. Here, we found that expression of rpoB(R) is regulated post-transcriptionally by RNA processing in the 5'-UTR of rpoB(R) mRNA, and that the endoribonuclease activity of Rid7C is responsible for mRNA processing thereby overseeing several tracts of morphological and biochemical differentiation. We also provide evidence that Rid7C may be associated with ribonuclease P M1 RNA, although M1 RNA is not required for rpoB(R) mRNA processing in vitro, and that Rid7C endoribonuclease activity is inhibited by A40926 suggesting the existence of a negative feedback loop on A40926 production, and a role of the endogenous synthesis of A40926 in the modulation of biochemical differentiation in this microorganism. Importance The YjgF/YER057c/UK114 family includes many proteins with diverse functions involved in detoxification, RNA maturation, and control of mRNA translation. We found that Rid7C is an endoribonuclease that is involved in processing of rpoB(R) mRNA, coding for a specialized RNA polymerase beta subunit that oversees morphological differentiation and A40926 antibiotic production in Nonomuraea gerenzanensis. Rid7C-mediated processing promotes rpoB(R) mRNA translation and antibiotic production, while Rid7C endoribonuclease activity is inhibited by A40926 suggesting a role of the endogenous synthesis of A40926 in modulation of biochemical differentiation in this microorganism. Finally, we show that recombinant Rid7C co-purified with M1 RNA (the RNA subunit of ribonuclease P) from Escherichia coli extract, suggesting a functional interaction between Rid7C and M1 RNA activities.

How do I get the most out of my protein sequence using bioinformatics tools?

Biochemical and biophysical experiments are essential for uncovering the three-dimensional structure and biological role of a protein of interest. However, meaningful predictions can frequently also be made using bioinformatics resources that transfer knowledge from a well studied protein to an uncharacterized protein based on their evolutionary relatedness. These predictions are helpful in developing specific hypotheses to guide wet-laboratory experiments. Commonly used bioinformatics resources include methods to identify and predict conserved sequence motifs, protein domains, transmembrane segments, signal sequences, and secondary as well as tertiary structure. Here, several such methods available through the MPI Bioinformatics Toolkit (https://toolkit.tuebingen.mpg.de) are described and how their combined use can provide meaningful information on a protein of unknown function is demonstrated. In particular, the identification of homologs of known structure using HHpred, internal repeats using HHrepID, coiled coils using PCOILS and DeepCoil, and transmembrane segments using Quick2D are focused on.

Discovery and mining of enzymes from the human gut microbiome

Advances in technological and bioinformatics approaches have led to the generation of a plethora of human gut metagenomic datasets. Metabolomics has also provided substantial data regarding the small metabolites produced and modified by the microbiota. Comparatively, the microbial enzymes mediating the transformation of metabolites have not been intensively investigated. Here, we discuss the recent efforts and technologies used for discovering and mining enzymes from the human gut microbiota. The wealth of knowledge on metabolites, reactions, genome sequences, and structures of proteins, may drive the development of strategies for enzyme mining. Ongoing efforts to annotate gut microbiota enzymes will explain catalytic mechanisms that may guide the clinical applications of the gut microbiome for diagnostic and therapeutic purposes.

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.