Evolutionary, computational, and biochemical studies of the salicylaldehyde dehydrogenases in the naphthalene degradation pathway

An animal charcoal contaminated cottage industry soil highlighted by halophilic archaea dominance and decimation of bacteria

An animal charcoal contaminated cottage industry soil in Lagos, Nigeria (ACGT) was compared in an ex post facto study with a nearby unimpacted soil (ACGC). Hydrocarbon content was higher than regulatory limits in ACGT (180.2 mg/kg) but lower in ACGC (19.28 mg/kg). Heavy metals like nickel, cadmium, chromium and lead were below detection limit in ACGC. However, all these metals, except cadmium, were detected in ACGT, but at concentrations below regulatory limits. Furthermore, copper (253.205 mg/kg) and zinc (422.630 mg/kg) were above regulatory limits in ACGT. Next generation sequencing revealed that the procaryotic community was dominated by bacteria in ACGC (62%) while in ACGT archaea dominated (76%). Dominant phyla in ACGC were Euryarchaeota (37%), Pseudomonadota (16%) and Actinomycetota (12%). In ACGT it was Euryarchaeota (76%), Bacillota (9%), Pseudomonadota (7%) and Candidatus Nanohaloarchaeota (5%). Dominant Halobacteria genera in ACGT were Halobacterium (16%), Halorientalis (16%), unranked halophilic archaeon (13%) Salarchaeum (6%) and Candidatus Nanohalobium (5%), whereas ACGC showed greater diversity dominated by bacterial genera Salimicrobium (7%) and Halomonas (3%). Heavy metals homeostasis genes, especially for copper, were fairly represented in both soils but with bacterial taxonomic affiliations. Sites like ACGT, hitherto poorly studied and understood, could be sources of novel bioresources.

Gut microbiota metabolically mediate intestinal helminth infection in zebrafish

Intestinal helminth parasite (IHP) infection induces alterations in the composition of microbial communities across vertebrates, although how gut microbiota may facilitate or hinder parasite infection remains poorly defined. In this work, we utilized a zebrafish model to investigate the relationship between gut microbiota, gut metabolites, and IHP infection. We found that extreme disparity in zebrafish parasite infection burden is linked to the composition of the gut microbiome and that changes in the gut microbiome are associated with variation in a class of endogenously produced signaling compounds, N-acylethanolamines, that are known to be involved in parasite infection. Using a statistical mediation analysis, we uncovered a set of gut microbes whose relative abundance explains the association between gut metabolites and infection outcomes. Experimental investigation of one of the compounds in this analysis reveals salicylaldehyde, which is putatively produced by the gut microbe Pelomonas, as a potent anthelmintic with activity against Pseudocapillaria tomentosa egg hatching, both in vitro and in vivo. Collectively, our findings underscore the importance of the gut microbiome as a mediating agent in parasitic infection and highlight specific gut metabolites as tools for the advancement of novel therapeutic interventions against IHP infection.

IMPORTANCE

Intestinal helminth parasites (IHPs) impact human health globally and interfere with animal health and agricultural productivity. While anthelmintics are critical to controlling parasite infections, their efficacy is increasingly compromised by drug resistance. Recent investigations suggest the gut microbiome might mediate helminth infection dynamics. So, identifying how gut microbes interact with parasites could yield new therapeutic targets for infection prevention and management. We conducted a study using a zebrafish model of parasitic infection to identify routes by which gut microbes might impact helminth infection outcomes. Our research linked the gut microbiome to both parasite infection and to metabolites in the gut to understand how microbes could alter parasite infection. We identified a metabolite in the gut, salicylaldehyde, that is putatively produced by a gut microbe and that inhibits parasitic egg growth. Our results also point to a class of compounds, N-acyl-ethanolamines, which are affected by changes in the gut microbiome and are linked to parasite infection. Collectively, our results indicate the gut microbiome may be a source of novel anthelmintics that can be harnessed to control IHPs.

Gut microbiota metabolically mediate intestinal helminth infection in Zebrafish

Intestinal helminth parasite (IHP) infection induces alterations in the composition of microbial communities across vertebrates, although how gut microbiota may facilitate or hinder parasite infection remains poorly defined. In this work we utilized a zebrafish model to investigate the relationship between gut microbiota, gut metabolites, and IHP infection. We found that extreme disparity in zebrafish parasite infection burden is linked to the composition of the gut microbiome, and that changes in the gut microbiome are associated with variation in a class of endogenously-produced signaling compounds, N-acylethanolamines, that are known to be involved in parasite infection. Using a statistical mediation analysis, we uncovered a set of gut microbes whose relative abundance explains the association between gut metabolites and infection outcomes. Experimental investigation of one of the compounds in this analysis reveals salicylaldehyde, which is putatively produced by the gut microbe Pelomonas, as a potent anthelmintic with activity against Pseudocapillaria tomentosa egg hatching, both in vitro and in vivo. Collectively, our findings underscore the importance of the gut microbiome as a mediating agent in parasitic infection and highlights specific gut metabolites as tools for the advancement of novel therapeutic interventions against IHP infection.

Cultivation and characterization of functional-yet-uncultivable phenanthrene degraders by stable-isotope-probing and metagenomic-binning directed cultivation (SIP-MDC)

High-throughput identification and cultivation of functional-yet-uncultivable microorganisms is a fundamental goal in environmental microbiology. It remains as a critical challenge due to the lack of routine and effective approaches. Here, we firstly proposed an approach of stable-isotope-probing and metagenomic-binning directed cultivation (SIP-MDC) to isolate and characterize the active phenanthrene degraders from petroleum-contaminated soils. From SIP and metagenome, we assembled 13 high-quality metagenomic bins from 13C-DNA, and successfully obtained the genome of an active PHE degrader Achromobacter (genome-MB) from 13C-DNA metagenomes, which was confirmed by gyrB gene comparison and average nucleotide/amino identity (ANI/AAI), as well as the quantification of PAH dioxygenase and antibiotic resistance genes. Thereinto, we modified the traditional cultivation medium with antibiotics and specific growth factors (e.g., vitamins and metals), and separated an active phenanthrene degrader Achromobacter sp. LJB-25 via directed isolation. Strain LJB-25 could degrade phenanthrene and its identity was confirmed by ANI/AAI values between its genome and genome-MB (>99 %). Our results hinted at the feasibility of SIP-MDC to identify, isolate and cultivate functional-yet-uncultivable microorganisms (active phenanthrene degraders) from their natural habitats. Our findings developed a state-of-the-art SIP-MDC approach, expanded our knowledge on phenanthrene biodegradation mechanisms, and proposed a strategy to mine functional-yet-uncultivable microorganisms.

In-depth comparative transcriptome analysis of Purpureocillium sp. CB1 under cadmium stress

Fungal bioremediation is a very attractive tool to cope with environmental pollution. We aimed to decipher the cadmium (Cd) response of Purpureocillium sp. CB1, isolated from polluted soil, at transcriptome level by RNA-sequencing (RNA-seq). We used 500 and 2500 mg/L of Cd²⁺ concentrations at two time points (t_6;36). RNA-seq determined 620 genes that were co-expressed in all samples. The highest number of differentially expressed genes (DEGs) was obtained within the first six h of exposure to 2500 mg/L of Cd²⁺. Several genes encoding transcriptional regulators, transporters, heat shock proteins, and oxidative stress-related genes were differentially expressed under Cd²⁺ stress. Remarkably, the genes that encode salicylate hydroxylase, which is involved in naphthalene biodegradation pathway, were significantly overexpressed. Utilization of diesel as the sole carbon source by CB1 even in the presence of Cd²⁺ supported concomitant upregulation of hydrocarbon degradation pathway genes. Furthermore, leucinostatin-related gene expression levels increased under Cd²⁺ stress. In addition, leucinostatin extracts from Cd²⁺-treated CB1 cultures showed higher antifungal activity than the control. Notably, Cd²⁺ in CB1 was mainly found as bound to the cell wall, thus confirming its adsorption potential. Cd²⁺ stress slightly reduced growth and led to mycelial malformation due to Cd²⁺ adsorption, especially at a concentration of 2500 mg/L at t₃₆. A strong correlation was recorded between RNA-seq and reverse-transcriptase-quantitative polymerase chain reaction (RT-qPCR) data. In conclusion, the study represents the first transcriptome analysis of Purpureocillium sp. under Cd²⁺ stress, providing insights into the primary targets for rational engineering to construct strains with remarkable bioremediation potency. KEY POINTS: • Upregulation of genes encoding salicylate hydroxylases under Cd²⁺ stress • Maximum Cd²⁺ adsorption at 500 mg/L at t₃₆ as tightly bound to the cell wall • Concordant bioremediation potential of CB1 on Cd²⁺ and diesel.

Molecular Docking of Lac_CB10: Highlighting the Great Potential for Bioremediation of Recalcitrant Chemical Compounds by One Predicted Bacteroidetes CopA-Laccase

Laccases are multicopper oxidases (MCOs) with a broad application spectrum, particularly in second-generation ethanol biotechnology and the bioremediation of xenobiotics and other highly recalcitrant compounds. Synthetic pesticides are xenobiotics with long environmental persistence, and the search for their effective bioremediation has mobilized the scientific community. Antibiotics, in turn, can pose severe risks for the emergence of multidrug-resistant microorganisms, as their frequent use for medical and veterinary purposes can generate constant selective pressure on the microbiota of urban and agricultural effluents. In the search for more efficient industrial processes, some bacterial laccases stand out for their tolerance to extreme physicochemical conditions and their fast generation cycles. Accordingly, to expand the range of effective approaches for the bioremediation of environmentally important compounds, the prospection of bacterial laccases was carried out from a custom genomic database. The best hit found in the genome of Chitinophaga sp. CB10, a Bacteroidetes isolate obtained from a biomass-degrading bacterial consortium, was subjected to in silico prediction, molecular docking, and molecular dynamics simulation analyses. The putative laccase CB10_180.4889 (Lac_CB10), composed of 728 amino acids, with theoretical molecular mass values of approximately 84 kDa and a pI of 6.51, was predicted to be a new CopA with three cupredoxin domains and four conserved motifs linking MCOs to copper sites that assist in catalytic reactions. Molecular docking studies revealed that Lac_CB10 had a high affinity for the molecules evaluated, and the affinity profiles with multiple catalytic pockets predicted the following order of decreasing thermodynamically favorable values: tetracycline (-8 kcal/mol) > ABTS (-6.9 kcal/mol) > sulfisoxazole (-6.7 kcal/mol) > benzidine (-6.4 kcal/mol) > trimethoprim (-6.1 kcal/mol) > 2,4-dichlorophenol (-5.9 kcal/mol) mol. Finally, the molecular dynamics analysis suggests that Lac_CB10 is more likely to be effective against sulfisoxazole-like compounds, as the sulfisoxazole-Lac_CB10 complex exhibited RMSD values lower than 0.2 nm, and sulfisoxazole remained bound to the binding site for the entire 100 ns evaluation period. These findings corroborate that LacCB10 has a high potential for the bioremediation of this molecule.

Bioremediation potential of laccase for catalysis of glyphosate, isoproturon, lignin, and parathion: Molecular docking, dynamics, and simulation

The present study aimed to investigate the catalytic degradation produced by laccase in the detoxification of glyphosate, isoproturon, lignin polymer, and parathion. We explored laccase-glyphosate, laccase-lignin polymer, laccase-isoproturon, and laccase-parathion using molecular docking (MD) and molecular dynamics simulation (MDS) approaches. The results suggest that laccase interacts well with glyphosate, lignin polymer, isoproturon, and parathion during biodegradation. We calculated the root mean square deviations (RMSD) of laccase-glyphosate, laccase-lignin polymer, laccase-isoproturon, and laccase-parathion as 0.24 ± 0.02, 0.59 ± 0.32, 0.43 ± 0.07, and 0.43 ± 0.06 nm, respectively. In an aqueous solution, the stability of laccase with glyphosate, lignin polymer, isoproturon, and parathion is mediated through the formation of hydrophobic interactions, hydrogen bonds, and van der Waals interactions. The presence of xenobiotic toxic compounds in the active site changed the conformation of laccase. MDS of the laccase-substrate complexes confirmed their stability during catalytic degradation. Laccase assay results confirmed that the degradation of syringol, dihydroconiferyl alcohol, guaiacol, parathion, isoproturon, and glyphosate were 100%, 99.31%, 95.69%, 60.96%, 54.51%, and 48.34% within 2 h, respectively. Taken together, we describe a novel method to understand the molecular-level biodegradation of xenobiotic compounds through laccase and its potential application in contaminant removal.

Hexaconazole exposure ravages biosynthesis pathway of steroid hormones: revealed by molecular dynamics and interaction

Widespread application of hexaconazole for agriculture purpose poses a threat to human health by disrupting normal endocrine homeostasis. To avoid adverse health effects on human, it is crucial to identify the effects of hexaconazole on key enzymes responsible for steroidal hormone synthesis. In view of this, present study was conducted to investigate the interaction mechanisms of hexaconazole with key enzymes in comparison with their food drug administration (FDA) approved inhibitor by molecular docking and molecular dynamics simulations. Results indicate that hexaconazole contacts with the active site of the key enzymes required for steroidal hormonal synthesis. Results pertaining to root-mean-square deviation, root-mean-square calculation, radius of gyration, hydrogen bonding and solvent accessible surface area exhibited that the interaction pattern and stability of interaction of hexaconazole was similar to enzyme specific inhibitor. In addition, ligand and enzyme complex interaction energy of hexaconazole was almost similar to key enzyme and FDA-approved enzyme specific inhibitor complex. This study offers a molecular level of understanding of hexaconazole with different enzymes required for steroidal hormonal synthesis. Findings of the study clearly suggest that hexaconazole has efficacy to stably interact with various enzyme required to progress the pathway of hormonal synthesis. If incessant exposure of hexaconazole occurs during agricultural work it may lead to ravage hormonal synthesis or potent endocrine disruption. The result of binding energy and complex interaction energy is depicted in the graphical abstract.

Unveiling metabolic characteristics of an uncultured Gammaproteobacterium responsible for in situ PAH biodegradation in petroleum polluted soil

Exploring the metabolic characteristics of indigenous PAH degraders is critical to understanding the PAH bioremediation mechanism in the natural environment. While stable-isotopic probing (SIP) is a viable method to identify functional microorganisms in complex environments, the metabolic characteristics of uncultured degraders are still elusive. Here, we investigated the naphthalene (NAP) biodegradation of petroleum polluted soils by combining SIP, amplicon sequencing and metagenome binning. Based on the SIP and amplicon sequencing results, an uncultured Gammaproteobacterium sp. was identified as the key NAP degrader. Additionally, the assembled genome of this uncultured degrader was successfully obtained from the ¹³ C-DNA metagenomes by matching its 16S rRNA gene with the SIP identified OTU sequence. Meanwhile, a number of NAP degrading genes encoding naphthalene/PAH dioxygenases were identified in this genome, further confirming the direct involvement of this indigenous degrader in the NAP degradation. The degrader contained genes related to the metabolisms of several carbon sources, energy substances and vitamins, illuminating potential reasons for why microorganisms cannot be cultivated and finally realize their cultivation. Our findings provide novel information on the mechanisms of in situ PAH biodegradation and add to our current knowledge on the cultivation of non-culturable microorganisms by combining both SIP and metagenome binning.

Binding interaction of glyphosate with glyphosate oxidoreductase and C-P lyase: Molecular docking and molecular dynamics simulation studies

Widespread application of glyphosate poses a threat to living organisms. Microbial strains are able to degrade glyphosate via contrasting metabolic pathways with the help of enzymes. Glyphosate oxidoreductase (GOX) and C-P lyase are the key enzymes for the biodegradation of glyphosate and its intermediate metabolite aminomethylphosphonic acid (AMPA) in microbes. The microbial degradation of glyphosate has been reported, but the underlying molecular mechanism is still unclear. Therefore, in this study, the interaction mechanism of GOX and C-P lyase with glyphosate and AMPA were investigated by using molecular docking and molecular dynamics (MD) simulations. The results indicate that glyphosate contacts with the active site of GOX and C-P lyase by hydrogen bonds as well as hydrophobic and van der Waals interactions in aqueous solution to maintain its stability. The presence of glyphosate and AMPA in the active site significantly changes the conformation of GOX and C-P lyase. The results of the MD simulations confirm that GOX and C-P lyase complexes are stable during the catalytic reaction. This study offers a molecular level of understanding of the expression and function of GOX and C-P lyase for the bioremediation of glyphosate.

Exploring the genetic potential of a fosmid metagenomic library from an oil-impacted mangrove sediment for metabolism of aromatic compounds

Aromatic hydrocarbons (AH) are widely distributed in nature, and many of them have been reported as relevant environmental pollutants and valuable carbon sources for different microorganisms. In this work, high-throughput sequencing of a metagenomic fosmid library was carried out to evaluate the functional and taxonomic diversity of genes involved in aromatic compounds degradation in oil-impacted mangrove sediments. In addition, activity-based approach and gas chromatography were used to assess the degradation potential of fosmid clones. Results indicated that AH degradation genes, such as monooxygenases and dioxygenases, were grouped into the following categories: anaerobic degradation of aromatic compounds (20.34%), metabolism of central aromatic intermediates (35.40%) and peripheral pathways for catabolism of aromatic compounds (22.56%). Taxonomic affiliation of genes related to aromatic compounds metabolism revealed the prevalence of the classes Alphaproteobacteria, Actinobacteria, Betaproteobacteria, Gammaproteobacteria and Deltaproteobacteria. Aromatic hydrocarbons (phenol, naphthalene, phenanthrene, pyrene and benzopyrene) were used as the only carbon source to screen clones with degradation potential. Of the 2500 clones tested, 48 showed some respiratory activity in at least one of the five carbon sources used. The hydrocarbon degradation ability of the top ten fosmid clones was confirmed by GC-MS. Further, annotation of assembled metagenomic fragments revealed ORFs corresponding to proteins and functional domains directly or indirectly involved in the aromatic compound metabolism, such as catechol 2,3-dioxygenase and ferredoxin oxidoreductase. Finally, these data suggest that the indigenous mangrove sediment microbiota developed essential mechanisms towards ecosystem remediation of petroleum hydrocarbon impact.

Exploration of the reduction mechanism of Cr(VI) in anaerobic hydrogen fermenter

The bio-reduction of hexavalent chromium (Cr(VI)) by anaerobic fermentation is considered as a promising, low-cost and environment-friendly way. However, it is unclear for the reduction mechanisms of Cr(VI) in an anaerobic hydrogen fermenter, such as reduction kinetics, related electron donors, migration and transformation, reduction site and key components, and related microorganisms. To clarify these issues, a hydrogen fermenter was designed to reduce Cr(VI) at 55 °C with glucose as initial substrate. Results show that 100 mg/L Cr(VI) can be completely reduced (99.5%) to trivalent chromium (Cr(III) through chemical and biological reactions. Bio-reduction dominates Cr(VI) removal in a first-order exponential decay mode with both glucose and its metabolites (volatile fatty acids) as electron donors. Moreover, volatile fatty acids are more suitable as electron donors for Cr(VI) bio-reduction than glucose. Bacilli, Clostridia and Thermotogae in the fermenter dominated the reduction of Cr(VI) by regulating the production and composition of extracellular polymers (EPSs), in which carboxyl and hydroxyl groups play an important role for Cr(VI) reduction by coordination. The results can guide us to regulate the bio-reduction of Cr(VI), and provide reference for the development of bio-reduction technology of Cr(VI).

Effect of surfactants on the interaction of phenol with laccase: Molecular docking and molecular dynamics simulation studies

Some surfactants can enhance the removal of phenol by laccase (Lac) in various industrial effluents. Their behavior and function in the biodegradation of phenolic wastewater have been experimentally reported by many researchers, but the underlying molecular mechanism is still unclear. Therefore, the interaction mechanisms of phenol with Lac from Trametes versicolor were investigated in the presence or absence of Triton X-100 (TX100) or rhamnolipid (RL) by molecular docking and molecular dynamics (MD) simulations. The results indicate that phenol contacts with an active site of Lac by hydrogen bonds (HBs) and van der Waals (vdW) interactions in aqueous solution for maintaining its stability. The presence of TX100 or RL results in the significant changes of enzymatic conformations. Meanwhile, the hydrophobic parts of surfactants contact with the outside surface of Lac. These changes lead to the decrease of binding energy between phenol and Lac. The migration behavior of water molecules within hydration shell is also inevitably affected. Therefore, the amphipathic TX100 or RL may influence the phenol degradation ability of Lac by modulating their interactions and water environment. This study offers molecular level of understanding on the function of surfactants in biosystem.

Family-Four Aldehyde Dehydrogenases Play an Indispensable Role in the Pathogenesis of Magnaporthe oryzae

The oxidative degradation of lipids through lipid peroxidation processes results in the generation of free fatty acid radicals. These free radicals including reactive oxygen species (ROS) serve as a substrate for generating reactive aldehydes. The accumulation of free fatty acid radicals, ROS, and reactive aldehydes in cell compartments beyond physiological threshold levels tends to exert a damaging effect on proximal membranes and distal tissues. Living organisms deploy a wide array of efficient enzymes including superoxide dismutase (SOD), catalase (CAT), peroxidase (POD), and aldehyde dehydrogenases (ALDHs) for scavenging reactive molecules and intermediates produced from membrane lipid peroxidation events. Although the contributions of SOD, CAT, and POD to the pathogenesis of microbial plant pathogens are well known, the influence of ALDH genes on the morphological and infectious development of plant pathogenic microbes is not well understood. In this study, we deployed RNA interference (RNAi) techniques and successfully silenced two putative family-four aldehyde dehydrogenase genes potassium-activated aldehyde dehydrogenase (MoKDCDH) and delta-1-pyrrorine-5-carboxylate dehydrogenase (MoP5CDH) in the rice blast pathogen Magnaporthe oryzae. The results obtained from the phenotypic analysis of individual knock-down strains showed that the RNAi-mediated inactivation of MoKDCDH and MoP5CDH triggered a significant reduction in conidiogenesis and vegetative growth of ΔMokdcdh and ΔMop5cdh strains. We further observed that downregulating the expression of MoKDCDH and MoP5CDH severely compromised the pathogenesis of the rice blast fungus. Also, the disruption of MoKDCDH and MoP5CDH M. oryzae undermined membrane integrity and rendered the mutant strains highly sensitive to membrane stress inducing osmolytes. However, the MoKDCDH and MoP5CDH knock-down strains generated in this study displayed unaltered cell wall integrity and thus suggested that family-four ALDHs play a dispensable role in enforcing cell wall-directed stress tolerance in M. oryzae. From these results, we deduced that family-four ALDHs play a conserved role in fostering membrane integrity in M. oryzae possibly by scavenging reactive aldehydes, fatty acid radicals, and other alcohol derivatives. The observation that downregulating the expression activities of MoKDCDH had a lethal effect on potential mutants further emphasized the need for comprehensive and holistic evaluation of the numerous ALDHs amassed by the rice blast fungus for their possible engagement as suitable targets as antiblast agents.

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

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

AGeNNT: annotation of enzyme families by means of refined neighborhood networks

BACKGROUND

Large enzyme families may contain functionally diverse members that give rise to clusters in a sequence similarity network (SSN). In prokaryotes, the genome neighborhood of a gene-product is indicative of its function and thus, a genome neighborhood network (GNN) deduced for an SSN provides strong clues to the specific function of enzymes constituting the different clusters. The Enzyme Function Initiative ( http://enzymefunction.org/ ) offers services that compute SSNs and GNNs.

RESULTS

We have implemented AGeNNT that utilizes these services, albeit with datasets purged with respect to unspecific protein functions and overrepresented species. AGeNNT generates refined GNNs (rGNNs) that consist of cluster-nodes representing the sequences under study and Pfam-nodes representing enzyme functions encoded in the respective neighborhoods. For cluster-nodes, AGeNNT summarizes the phylogenetic relationships of the contributing species and a statistic indicates how unique nodes and GNs are within this rGNN. Pfam-nodes are annotated with additional features like GO terms describing protein function. For edges, the coverage is given, which is the relative number of neighborhoods containing the considered enzyme function (Pfam-node). AGeNNT is available at https://github.com/kandlinf/agennt .

CONCLUSIONS

An rGNN is easier to interpret than a conventional GNN, which commonly contains proteins without enzymatic function and overly specific neighborhoods due to phylogenetic bias. The implemented filter routines and the statistic allow the user to identify those neighborhoods that are most indicative of a specific metabolic capacity. Thus, AGeNNT facilitates to distinguish and annotate functionally different members of enzyme families.

Integrative View of the Diversity and Evolution of SWEET and SemiSWEET Sugar Transporters

Sugars Will Eventually be Exported Transporter (SWEET) and SemiSWEET are recently characterized families of sugar transporters in eukaryotes and prokaryotes, respectively. SemiSWEETs contain 3 transmembrane helices (TMHs), while SWEETs contain 7. Here, we performed sequence-based comprehensive analyses for SWEETs and SemiSWEETs across the biosphere. In total, 3,249 proteins were identified and ≈60% proteins were found in green plants and Oomycota, which include a number of important plant pathogens. Protein sequence similarity networks indicate that proteins from different organisms are significantly clustered. Of note, SemiSWEETs with 3 or 4 TMHs that may fuse to SWEET were identified in plant genomes. 7-TMH SWEETs were found in bacteria, implying that SemiSWEET can be fused directly in prokaryote. 15-TMH extraSWEET and 25-TMH superSWEET were also observed in wild rice and oomycetes, respectively. The transporters can be classified into 4, 2, 2, and 2 clades in plants, Metazoa, unicellular eukaryotes, and prokaryotes, respectively. The consensus and coevolution of amino acids in SWEETs were identified by multiple sequence alignments. The functions of the highly conserved residues were analyzed by molecular dynamics analysis. The 19 most highly conserved residues in the SWEETs were further confirmed by point mutagenesis using SWEET1 from Arabidopsis thaliana. The results proved that the conserved residues located in the extrafacial gate (Y57, G58, G131, and P191), the substrate binding pocket (N73, N192, and W176), and the intrafacial gate (P43, Y83, F87, P145, M161, P162, and Q202) play important roles for substrate recognition and transport processes. Taken together, our analyses provide a foundation for understanding the diversity, classification, and evolution of SWEETs and SemiSWEETs using large-scale sequence analysis and further show that gene duplication and gene fusion are important factors driving the evolution of SWEETs.