Phylogenomics reveals subfamilies of fungal nonribosomal peptide synthetases and their evolutionary relationships

Background  

Nonribosomal peptide synthetases (NRPSs) are multimodular enzymes, found in fungi and bacteria, which biosynthesize peptides without the aid of ribosomes. Although their metabolite products have been the subject of intense investigation due to their life-saving roles as medicinals and injurious roles as mycotoxins and virulence factors, little is known of the phylogenetic relationships of the corresponding NRPSs or whether they can be ranked into subgroups of common function. We identified genes (NPS) encoding NRPS and NRPS-like proteins in 38 fungal genomes and undertook phylogenomic analyses in order to identify fungal NRPS subfamilies, assess taxonomic distribution, evaluate levels of conservation across subfamilies, and address mechanisms of evolution of multimodular NRPSs. We also characterized relationships of fungal NRPSs, a representative sampling of bacterial NRPSs, and related adenylating enzymes, including α-aminoadipate reductases (AARs) involved in lysine biosynthesis in fungi.     

In silico analysis of the adenylation domains of the freestanding enzymes belonging to the eucaryotic nonribosomal peptide synthetase-like family

This work presents a computational analysis of the molecular characteristics shared by the adenylation domains from traditional nonribosomal peptide synthetases (NRPSs) and the group of the freestanding homologous enzymes: alpha-aminoadipate semialdehyde dehydrogenase, alpha-aminoadipate reductase and the protein Ebony. The results of systematic sequence comparisons allow us to conclude that a specificity-conferring code, similar to that described for the NRPSs, can be recognized in such enzymes. The structural and functional roles of the residues involved in the substrate selection and binding are proposed through the analysis of the predicted interactions of the model active sites and their respective substrates. The indications deriving from this study can be useful for the programming of experiments aimed at a better characterization and at the engineering of this emerging group of single NRPS modules that are responsible for amino acid selection, activation and modification in the absence of other NRPS assembly line components.     

Evolution of Chemical Diversity in Echinocandin Lipopeptide Antifungal Metabolites

The echinocandins are a class of antifungal drugs that includes caspofungin, micafungin, and anidulafungin. Gene clusters encoding most of the structural complexity of the echinocandins provided a framework for hypotheses about the evolutionary history and chemical logic of echinocandin biosynthesis. Gene orthologs among echinocandin-producing fungi were identified. Pathway genes, including the nonribosomal peptide synthetases (NRPSs), were analyzed phylogenetically to address the hypothesis that these pathways represent descent from a common ancestor. The clusters share cooperative gene contents and linkages among the different strains. Individual pathway genes analyzed in the context of similar genes formed unique echinocandin-exclusive phylogenetic lineages. The echinocandin NRPSs, along with the NRPS from the inp gene cluster in Aspergillus nidulans and its orthologs, comprise a novel lineage among fungal NRPSs. NRPS adenylation domains from different species exhibited a one-to-one correspondence between modules and amino acid specificity that is consistent with models of tandem duplication and subfunctionalization. Pathway gene trees and Ascomycota phylogenies are congruent and consistent with the hypothesis that the echinocandin gene clusters have a common origin. The disjunct Eurotiomycete-Leotiomycete distribution appears to be consistent with a scenario of vertical descent accompanied by incomplete lineage sorting and loss of the clusters from most lineages of the Ascomycota. We present evidence for a single evolutionary origin of the echinocandin family of gene clusters and a progression of structural diversification in two fungal classes that diverged approximately 290 to 390 million years ago. Lineage-specific gene cluster evolution driven by selection of new chemotypes contributed to diversification of the molecular functionalities.     

The polynucleotide ligase and RNA capping enzyme superfamily of covalent nucleotidyltransferases

ATP- and NAD(+)-dependent DNA ligases, ATP-dependent RNA ligases and GTP-dependent mRNA capping enzymes comprise a superfamily of proteins that catalyze nucleotidyl transfer to polynucleotide 5' ends via covalent enzyme-(lysyl-N)-NMP intermediates. The superfamily is defined by five peptide motifs that line the nucleotide-binding pocket and contribute amino acid sidechains essential for catalysis. Early crystal structures revealed a shared core tertiary structure for DNA ligases and capping enzymes, which are composed minimally of a nucleotidyltransferase domain fused to a distal OB-fold domain. Recent structures of viral and bacterial DNA ligases, and a fungal mRNA capping enzyme underscore how the substrate-binding and chemical steps of the ligation and capping pathways are coordinated with large rearrangements of the component protein domains and with remodeling of the atomic contacts between the enzyme and the nucleotide at the active site. The first crystal structure of an RNA ligase suggests that contemporary DNA ligases, RNA ligases and RNA capping enzymes evolved by fusion of ancillary effector domains to an ancestral catalytic module involved in RNA repair.     

非核糖体肽合成酶(NRPSs)作用机理与应用的研究进展

许多微生物能利用非核糖体肽合成酶(NRPSs)合成结构复杂、种类繁多的的生物活性肽。非核糖体肽因其独特的理化特性和药理学特性已被广泛关注,极具商业开发潜力。NRPSs由多个模块组成,模块的不同空间排列顺序决定其多肽产物的氨基酸序列特异性。NRPSs以多载体巯基化模板机理进行多肽合成,其底物特异性由腺苷酰化结构域和缩合结构域共同实现。目前,人们已经利用天然的NRPSs、某些特定结构域、将已知NRPSs的模块或特定结构域进行组合甚至杂合组合而构建成的新的NRPSs来合成目的多肽。     

Four homologous domains in the primary structure of GrsB are related to domains in a superfamily of adenylate-forming enzymes

The entire nucleotide sequence of the Bacillus brevis grsB gene encoding the gramicidin S synthetase 2, which activates and condenses the four amino acids proline, valine, ornithine and leucine has been determined. The gene contains an open reading frame of 13,359 bp which encodes a protein of 4453 amino acids with a predicted Mr of 510,287. The gene is located within the gramicidin S biosynthetic operon, also containing the genes grsT and grsA, whose nucleotide sequences have been determined previously. Within the GrsB amino acid sequence four conserved and repeated domains of about 600 amino acids (45-50% identity) have been identified. The four domains are separated by non-homologous sequences of about 500 amino acids. The domains also share a high degree of similarity (20-70%) with eight peptide synthetases of bacterial and fungal origin as well as with conserved sequences of nine other adenylate-forming enzymes of diverse origin. On the basis of sequence homology and functional similarities, we infer that those enzymes share a common evolutionary origin and present a phylogenetic tree for this superfamily of domain-bearing enzymes.     

Molecular function prediction for a family exhibiting evolutionary tendencies toward substrate specificity swapping: Recurrence of tyrosine aminotransferase activity in the Iα subfamily

The subfamily Iα aminotransferases are typically categorized as having narrow specificity toward carboxylic amino acids (AATases), or broad specificity that includes aromatic amino acid substrates (TATases). Because of their general role in central metabolism and, more specifically, their association with liver‐related diseases in humans, this subfamily is biologically interesting. The substrate specificities for only a few members of this subfamily have been reported, and the reliable prediction of substrate specificity from protein sequence has remained elusive. In this study, a diverse set of aminotransferases was chosen for characterization based on a scoring system that measures the sequence divergence of the active site. The enzymes that were experimentally characterized include both narrow‐specificity AATases and broad‐specificity TATases, as well as AATases with broader‐specificity and TATases with narrower‐specificity than the previously known family members. Molecular function and phylogenetic analyses underscored the complexity of this family's evolution as the TATase function does not follow a single evolutionary thread, but rather appears independently multiple times during the evolution of the subfamily. The additional functional characterizations described in this article, alongside a detailed sequence and phylogenetic analysis, provide some novel clues to understanding the evolutionary mechanisms at work in this family. Proteins 2013. © 2013 Wiley Periodicals, Inc.     

Phylogenomic Relationships between Amylolytic Enzymes from 85 Strains of Fungi

Fungal amylolytic enzymes, including α-amylase, gluocoamylase and α-glucosidase, have been extensively exploited in diverse industrial applications such as high fructose syrup production, paper making, food processing and ethanol production. In this paper, amylolytic genes of 85 strains of fungi from the phyla Ascomycota, Basidiomycota, Chytridiomycota and Zygomycota were annotated on the genomic scale according to the classification of glycoside hydrolase (GH) from the Carbohydrate-Active enZymes (CAZy) Database. Comparisons of gene abundance in the fungi suggested that the repertoire of amylolytic genes adapted to their respective lifestyles. Amylolytic enzymes in family GH13 were divided into four distinct clades identified as heterologous α- amylases, eukaryotic α-amylases, bacterial and fungal α-amylases and GH13 α-glucosidases. Family GH15 had two branches, one for gluocoamylases, and the other with currently unknown function. GH31 α-glucosidases showed diverse branches consisting of neutral α-glucosidases, lysosomal acid α-glucosidases and a new clade phylogenetically related to the bacterial counterparts. Distribution of starch-binding domains in above fungal amylolytic enzymes was related to the enzyme source and phylogeny. Finally, likely scenarios for the evolution of amylolytic enzymes in fungi based on phylogenetic analyses were proposed. Our results provide new insights into evolutionary relationships among subgroups of fungal amylolytic enzymes and fungal evolutionary adaptation to ecological conditions.     

Functional and Phylogenetic Divergence of Fungal Adenylate-Forming Reductases

A key step in fungal l-lysine biosynthesis is catalyzed by adenylate-forming l-α-aminoadipic acid reductases, organized in domains for adenylation, thiolation, and the reduction step. However, the genomes of numerous ascomycetes and basidiomycetes contain an unexpectedly large number of additional genes encoding similar but functionally distinct enzymes. Here, we describe the functional in vitro characterization of four reductases which were heterologously produced in Escherichia coli. The Ceriporiopsis subvermispora serine reductase Nps1 features a terminal ferredoxin-NADP⁺ reductase (FNR) domain and thus belongs to a hitherto undescribed class of fungal multidomain enzymes. The second major class is characterized by the canonical terminal short-chain dehydrogenase/reductase domain and represented by Ceriporiopsis subvermispora Nps3 as the first biochemically characterized l-α-aminoadipic acid reductase of basidiomycete origin. Aspergillus flavus l-tyrosine reductases LnaA and LnbA are members of a distinct phylogenetic clade. Phylogenetic analysis supports the view that fungal adenylate-forming reductases are more diverse than previously recognized and belong to four distinct classes.     

Purification and characterization of benzoate-coenzyme A ligase and 2-aminobenzoate-coenzyme A ligases from a denitrifying Pseudomonas sp.

The enzymes catalyzing the formation of coenzyme A (CoA) thioesters of benzoate and 2-aminobenzoate were studied in a denitrifying Pseudomonas sp. anaerobically grown with these aromatic acids and nitrate as sole carbon and energy sources. Three different rather specific aromatic acyl-CoA ligases, E1, E2, and E3, were found which catalyze the formation of CoA thioesters of benzoate, fluorobenzoates, and 2-aminobenzoate. ATP is cleaved into AMP and pyrophosphate. The enzymes were purified, their N-terminal amino acid sequences were determined, and their catalytic and molecular properties were studied. Cells anaerobically grown on benzoate and nitrate contain one CoA ligase (AMP forming) for benzoic acid (E1). It is a homodimer of Mr 120,000 which prefers benzoate as a substrate but shows some activity also with 2-aminobenzoate and fluorobenzoates, although with lower Km. Cells anaerobically grown on 2-aminobenzoate and nitrate contain three different CoA ligases for aromatic acids. The first one is identical with benzoate-CoA ligase (E1). The second enzyme is a 2-aminobenzoate-CoA ligase (E2). It is a monomer of Mr 60,000 which prefers 2-aminobenzoate but also activates benzoate, fluorobenzoates and, less effectively, 2-methylbenzoate, with lower affinities to the latter substrates. The enzymes E1 and E2 have similar activity levels; a third minor CoA ligase activity is due to a different 2-aminobenzoate-CoA ligase. The enzyme (E3) is a monomer of Mr, 65,000 which 2-aminobenzoate pathway (U. Altenschmidt, C. Eckerskorn, and G. Fuchs, Eur. J. Biochem. 194:647-653, 1990); apparently, it is not completely repressed under anaerobic conditions and therefore also is induced to a small extent by 2-aminobenzoate under anaerobic growth conditions.     

Interkingdom gene transfer of a hybrid NPS/PKS from bacteria to filamentous Ascomycota

Nonribosomal peptides (NRPs) and polyketides (PKs) are ecologically important secondary metabolites produced by bacteria and fungi using multidomain enzymes called nonribosomal peptide synthetases (NRPSs) and polyketide synthases (PKSs), respectively. Previous phylogenetic analyses of fungal NRPSs and PKSs have suggested that a few of these genes were acquired by fungi via horizontal gene transfer (HGT) from bacteria, including a hybrid NPS/PKS found in Cochliobolus heterostrophus (Dothideomycetes, Ascomycota). Here, we identify this hybrid gene in fungi representing two additional classes of Ascomycota (Aspergillus spp., Microsporum canis, Arthroderma spp., and Trichophyton spp., Eurotiomycetes; Chaetomium spp. and Metarhizium spp., Sordariomycetes) and use phylogenetic analyses of the most highly conserved domains from NRPSs (adenylation (A) domain) and PKSs (ketoacyl synthase (KS) domain) to examine the hypothesis that the hybrid NPS7/PKS24 was acquired by fungi from bacteria via HGT relatively early in the evolution of the Pezizomycotina. Our results reveal a unique ancestry of the A domain and KS domain in the hybrid gene relative to known fungal NRPSs and PKSs, provide strong evidence for HGT of the hybrid gene from a putative bacterial donor in the Burkholderiales, and suggest the HGT event occurred early in the evolution of the filamentous Ascomycota.     

Functional domains of assimilatory nitrate reductases and nitrite reductases

Biochemical investigation of nitrate assimilation enzymes spans the past four decades. With the molecular cloning of genes for nitrate reductases and nitrite reductases, exciting new prospects are developing for the study of these enzymes. As large, complex enzymes with multiple redox centers, these two types of reductases should help us gain understanding of structural, functional and evolutionary relationships among the diverse group of multicenter redox enzymes.     

Development of the genetic code: Insights from a fungal codon reassignment

The high conservation of the genetic code and its fundamental role in genome decoding suggest that its evolution is highly restricted or even frozen. However, various prokaryotic and eukaryotic genetic code alterations, several alternative tRNA-dependent amino acid biosynthesis pathways, regulation of tRNA decoding by diverse nucleoside modifications and recent in vivo incorporation of non-natural amino acids into prokaryotic and eukaryotic proteins, show that the code evolves and is surprisingly flexible. The cellular mechanisms and the proteome buffering capacity that support such evolutionary processes remain unclear. Here we explore the hypothesis that codon misreading and reassignment played fundamental roles in the development of the genetic code and we show how a fungal codon reassignment is enlightening its evolution.     

Phosphopantetheinyl transferase CfwA/NpgA is required for Aspergillus nidulans secondary metabolism and asexual development

Polyketide synthases (PKSs) and/or nonribosomal peptide synthetases (NRPSs) are central components of secondary metabolism in bacteria, plants, and fungi. In filamentous fungi, diverse PKSs and NRPSs participate in the biosynthesis of secondary metabolites such as pigments, antibiotics, siderophores, and mycotoxins. However, many secondary metabolites as well as the enzymes involved in their production are yet to be discovered. Both PKSs and NRPSs require activation by enzyme members of the 4'-phosphopantetheinyl transferase (PPTase) family. Here, we report the isolation and characterization of Aspergillus nidulans strains carrying conditional (cfwA2) and null (DeltacfwA) mutant alleles of the cfwA gene, encoding an essential PPTase. We identify the polyketides shamixanthone, emericellin, and dehydroaustinol as well as the sterols ergosterol, peroxiergosterol, and cerevisterol in extracts from A. nidulans large-scale cultures. The PPTase CfwA/NpgA was required for the production of these polyketide compounds but dispensable for ergosterol and cerevisterol and for fatty acid biosynthesis. The asexual sporulation defects of cfwA, DeltafluG, and DeltatmpA mutants were not rescued by the cfwA-dependent compounds identified here. However, a cfwA2 mutation enhanced the sporulation defects of both DeltatmpA and DeltafluG single mutants, suggesting that unidentified CfwA-dependent PKSs and/or NRPSs are involved in the production of hitherto-unknown compounds required for sporulation. Our results expand the number of known and predicted secondary metabolites requiring CfwA/NpgA for their biosynthesis and, together with the phylogenetic analysis of fungal PPTases, suggest that a single PPTase is responsible for the activation of all PKSs and NRPSs in A. nidulans.     

Structure of PA1221, a nonribosomal peptide synthetase containing adenylation and peptidyl carrier protein domains

Many bacteria use large modular enzymes for the synthesis of polyketide and peptide natural products. These multidomain enzymes contain integrated carrier domains that deliver bound substrates to multiple catalytic domains, requiring coordination of these chemical steps. Nonribosomal peptide synthetases (NRPSs) load amino acids onto carrier domains through the activity of an upstream adenylation domain. Our lab recently determined the structure of an engineered two-domain NRPS containing fused adenylation and carrier domains. This structure adopted a domain-swapped dimer that illustrated the interface between these two domains. To continue our investigation, we now examine PA1221, a natural two-domain protein from Pseudomonas aeruginosa. We have determined the amino acid specificity of this new enzyme and used domain specific mutations to demonstrate that loading the downstream carrier domain within a single protein molecule occurs more quickly than loading of a nonfused carrier domain intermolecularly. Finally, we have determined crystal structures of both apo- and holo-PA1221 proteins, the latter using a valine-adenosine vinylsulfonamide inhibitor that traps the adenylation domain-carrier domain interaction. The protein adopts an interface similar to that seen with the prior adenylation domain-carrier protein construct. A comparison of these structures with previous structures of multidomain NRPSs suggests that a large conformational change within the NRPS adenylation domains guides the carrier domain into the active site for thioester formation.     

Molecular Evolution Within the L-Malate and L-Lactate Dehydrogenase Super-Family

The NAD(P)-dependent malate (L-MalDH) and NAD-dependent lactate (L-LDH) form a large super-family that has been characterized in organisms belonging to the three domains of life. In the first part of this study, the group of [LDH-like] L-MalDH, which are malate dehydrogenases resembling lactate dehydrogenase, were analyzed and clearly defined with respect to the other enzymes. In the second part, the phylogenetic relationships of the whole super-family were presented by taking into account the [LDH-like] L-MalDH. The inferred tree unambiguously shows that two ancestral genes duplications, and not one as generally thought, are needed to explain both the distribution into two enzymatic functions and the observation of three main groups within the super-family: L-LDH, [LDH-like] L-MalDH, and dimeric L-MalDH. In addition, various cases of functional changes within each group were observed and analyzed. The direction of evolution was found to always be polarized: from enzymes with a high stringency of substrate recognition to enzymes with a broad substrate specificity. A specific phyletic distribution of the L-LDH, [LDH-like] L-MalDH, and dimeric L-MalDH over the Archaeal, Bacterial, and Eukaryal domains was observed. This was analyzed in the light of biochemical, structural, and genomic data available for the L-LDH, [LDH-like] L-MalDH, and dimeric L-MalDH. This analysis led to the elaboration of a refined evolutionary scenario of the super-family, in which the selection of L-LDH and the fate of L-MalDH during mitochrondrial genesis are presented.     

Development of a radiolabeled ATP assay for carboxylic acid:CoA Ligases and its use in the characterization of the xenobiotic carboxylic acid:CoA ligases of bovine liver mitochondria

A radiolabeled ATP assay was developed for measuring carboxylic acid:CoA ligase activity. The assay was designed to measure the formation of [γ-³³P]pyrophosphate from [γ-³³P]ATP in the course of the reaction. The assay was linear with protein concentration, and rates as low as 1 pmol/min were measurable. Rates determined with this assay were in agreement with rates determined with [14C]carboxylic acids. The assay was used to characterize the substrate specificity of the XL-I, XL-II, and XL-III ligases from bovine liver mitochondria. Forty carboxylic acids were tested for activity. The enzymes differed in their substrate specificities with XL-I and XL-II being the most similar and XL-III having the broadest specificity. This study has uncovered 19 new carboxylic acids that are substrates for these enzymes. © 1998 John Wiley & Sons, Inc. J Biochem Toxicol 12: 151–155, 1998     

Phylogenetic analysis of the UDP-glycosyltransferase multigene family of Arabidopsis thaliana

A class of UDP-glycosyltransferases (UGTs) defined by the presence of a C-terminal consensus sequence is found throughout the plant and animal kingdoms. Whereas mammalian enzymes use UDP-glucuronic acid, the plant enzymes typically use UDP-glucose in the transfer reactions. A diverse array of aglycones can be glucosylated by these UGTs. In plants, the aglycones include plant hormones, secondary metabolites involved in stress and defense responses, and xenobiotics such as herbicides. Glycosylation is known to regulate many properties of the aglycones such as their bioactivity, their solubility, and their transport properties within the cell and throughout the plant. As a means of providing a framework to start to understand the substrate specificities and structure-function relationships of plant UGTs, we have now applied a molecular phylogenetic analysis to the multigene family of 99 UGT sequences in Arabidopsis. We have determined the overall organization and evolutionary relationships among individual members with a surprisingly high degree of confidence. Through constructing a composite phylogenetic tree that also includes all of the additional plant UGTs with known catalytic activities, we can start to predict both the evolutionary history and substrate specificities of new sequences as they are identified. The tree already suggests that while the activities of some subgroups of the UGT family are highly conserved among different plant species, others subgroups shift substrate specificity with relative ease.     

Module evolution and substrate specificity of fungal nonribosomal peptide synthetases involved in siderophore biosynthesis

Background  

Most filamentous ascomycete fungi produce high affinity iron chelators called siderophores, biosynthesized nonribosomally by multimodular adenylating enzymes called nonribosomal peptide synthetases (NRPSs). While genes encoding the majority of NRPSs are intermittently distributed across the fungal kingdom, those encoding ferrichrome synthetase NRPSs, responsible for biosynthesis of ferrichrome siderophores, are conserved, which offers an opportunity to trace their evolution and the genesis of their multimodular domain architecture. Furthermore, since the chemistry of many ferrichromes is known, the biochemical and structural 'rules' guiding NRPS substrate choice can be addressed using protein structural modeling and evolutionary approaches.