Deep sequencing of non-ribosomal peptide synthetases and polyketide synthases from the microbiomes of Australian marine sponges

The biosynthesis of non-ribosomal peptide and polyketide natural products is facilitated by multimodular enzymes that contain domains responsible for the sequential condensation of amino and carboxylic subunits. These conserved domains provide molecular targets for the discovery of natural products from microbial metagenomes. This study demonstrates the application of tag-encoded FLX amplicon pyrosequencing (TEFAP) targeting non-ribosomal peptide synthetase (NRPS) and polyketide synthase (PKS) genes as a method for determining the identity and diversity of natural product biosynthesis genes. To validate this approach, we assessed the diversity of NRPS and PKS genes within the microbiomes of six Australian marine sponge species using both TEFAP and metagenomic whole-genome shotgun sequencing approaches. The TEFAP approach identified 100 novel ketosynthase (KS) domain sequences and 400 novel condensation domain sequences within the microbiomes of the six sponges. The diversity of KS domains within the microbiome of a single sponge species Scopalina sp. exceeded that of any previously surveyed marine sponge. Furthermore, this study represented the first to target the condensation domain from NRPS biosynthesis and resulted in the identification of a novel condensation domain lineage. This study highlights the untapped potential of Australian marine sponges for the isolation of novel bioactive natural products. Furthermore, this study demonstrates that TEFAP approaches can be applied to functional genes, involved in natural product biosynthesis, as a tool to aid natural product discovery. It is envisaged that this approach will be used across multiple environments, offering an insight into the biological processes that influence the production of secondary metabolites.     

Polyketide synthases of bacterial symbionts in sponges – Evolution-based applications in natural products research

Marine sponges are an unusually rich source of bioactive natural products with clinical potential. They also often harbor rich communities of symbiotic bacteria that have often been suspected as the true producers of sponge-derived compounds. To date, these bacteria can in most cases not be cultivated, but culture-independent methods, such as isolating and analyzing biosynthetic gene clusters using metagenomic strategies, have recently provided first insights into their chemical potential. This review summarizes recent work of our laboratory on the study of polyketide synthases (PKSs). These studies revealed two evolutionarily distinct, unusual PKS types that are commonly found in sponge metagenomes and were shown to be of bacterial origin. One, the sup PKS, dominates sponge metagenomic DNA libraries, occurs widespread in bacteriosponges and is to date exclusively known from such animals. Data suggest that it is a type of synthase that generates methyl-branched fatty acids, which are commonly present in sponges. The other PKS type, termed trans-acyltransferase (AT) PKS, is responsible for the biosynthesis of complex, bioactive polyketides, such as the onnamides, and also occurs in free-living bacteria. The diversity of PKS genes present in a single sponge metagenome can be enormous. However, the phylogenetic approaches outlined in this review can provide valuable insights into the PKS function and structures of polyketides and can assist in the targeted isolation of gene clusters.     

Structural insights into biosynthesis of resorcinolic lipids by a type III polyketide synthase in Neurospora crassa

Microbial type III polyketide synthases (PKSs) have revealed remarkable mechanistic as well as functional versatility. Recently, a type III PKS homolog from Azotobacter has been implicated in the biosynthesis of resorcinolic lipids, thus adding a new functional significance to this class of proteins. Here, we report the structural and mutational investigations of a novel type III PKS protein from Neurospora crassa involved in the biosynthesis of resorcinolic metabolites by utilizing long chain fatty acyl-CoAs. The structure revealed a long hydrophobic tunnel responsible for its fatty acyl chain length specificity resembling that of PKS18, a mycobacterial type III PKS. Structure-based mutational studies to block the tunnel not only altered the fatty acyl chain specificity but also resulted in change of cyclization pattern affecting the product profile. This first structural characterization of a resorcinolic lipid synthase provides insights into the coordinated functioning of cyclization and a substrate-binding pocket, which shows mechanistic intricacy underlying type III PKS catalysis.     

Towards Prediction of Metabolic Products of Polyketide Synthases: An In Silico Analysis

Sequence data arising from an increasing number of partial and complete genome projects is revealing the presence of the polyketide synthase (PKS) family of genes not only in microbes and fungi but also in plants and other eukaryotes. PKSs are huge multifunctional megasynthases that use a variety of biosynthetic paradigms to generate enormously diverse arrays of polyketide products that posses several pharmaceutically important properties. The remarkable conservation of these gene clusters across organisms offers abundant scope for obtaining novel insights into PKS biosynthetic code by computational analysis. We have carried out a comprehensive in silico analysis of modular and iterative gene clusters to test whether chemical structures of the secondary metabolites can be predicted from PKS protein sequences. Here, we report the success of our method and demonstrate the feasibility of deciphering the putative metabolic products of uncharacterized PKS clusters found in newly sequenced genomes. Profile Hidden Markov Model analysis has revealed distinct sequence features that can distinguish modular PKS proteins from their iterative counterparts. For iterative PKS proteins, structural models of iterative ketosynthase (KS) domains have revealed novel correlations between the size of the polyketide products and volume of the active site pocket. Furthermore, we have identified key residues in the substrate binding pocket that control the number of chain extensions in iterative PKSs. For modular PKS proteins, we describe for the first time an automated method based on crucial intermolecular contacts that can distinguish the correct biosynthetic order of substrate channeling from a large number of non-cognate combinatorial possibilities. Taken together, our in silico analysis provides valuable clues for formulating rules for predicting polyketide products of iterative as well as modular PKS clusters. These results have promising potential for discovery of novel natural products by genome mining and rational design of novel natural products.     

A Single Module Type I Polyketide Synthase Directs de Novo Macrolactone Biogenesis during Galbonolide Biosynthesis in Streptomyces galbus

Galbonolide (GAL) A and B are antifungal macrolactone polyketides produced by Streptomyces galbus. During their polyketide chain assembly, GAL-A and -B incorporate methoxymalonate and methylmalonate, respectively, in the fourth chain extension step. The methoxymalonyl-acyl carrier protein biosynthesis locus (galG to K) is specifically involved in GAL-A biosynthesis, and this locus is neighbored by a gene cluster composed of galA-E. GalA-C constitute a single module, highly reducing type I polyketide synthase (PKS). GalD and GalE are cytochrome P450 and Rieske domain protein, respectively. Gene knock-out experiments verified that galB, -C, and -D are essential for GAL biosynthesis. A galD mutant accumulated a GAL-C that lacked two hydroxyl groups and a double bond when compared with GAL-B. A [U-¹³C]propionate feeding experiment indicated that no rare precursor other than methoxymalonate was incorporated during GAL biogenesis. A search of the S. galbus genome for a modular type I PKS system, the type that was expected to direct GAL biosynthesis, resulted in the identification of only one modular type I PKS gene cluster. Homology analysis indicated that this PKS gene cluster is the locus for vicenistatin biosynthesis. This cluster was previously reported in Streptomyces halstedii. A gene deletion of the vinP2 ortholog clearly demonstrated that this modular type I PKS system is not involved in GAL biosynthesis. Therefore, we propose that GalA-C direct macrolactone polyketide formation for GAL. Our studies provide a glimpse into a novel biochemical strategy used for polyketide synthesis; that is, the iterative assembly of propionates with highly programmed β-keto group modifications.     

Diversity and evolution of polyketide biosynthesis gene clusters in the Ceratocystidaceae

Polyketides are secondary metabolites with diverse biological activities. Polyketide synthases (PKS) are often encoded from genes clustered in the same genomic region. Functional analyses and genomic studies show that most fungi are capable of producing a repertoire of polyketides. We considered the potential of Ceratocystidaceae for producing polyketides using a comparative genomics approach. Our aims were to identify the putative polyketide biosynthesis gene clusters, to characterize them and predict the types of polyketide compounds they might produce. We used sequences from nineteen species in the genera, Ceratocystis, Endoconidiophora, Davidsoniella, Huntiella, Thielaviopsis and Bretziella, to identify and characterize PKS gene clusters, by employing a range of bioinformatics and phylogenetic tools. We showed that the genomes contained putative clusters containing a non-reducing type I PKS and a type III PKS. Phylogenetic analyses suggested that these genes were already present in the ancestor of the Ceratocystidaceae. By contrast, the various reducing type I PKS-containing clusters identified in these genomes appeared to have distinct evolutionary origins. Although one of the identified clusters potentially allows for the production of melanin, their functional characterization will undoubtedly reveal many novel and important compounds implicated in the biology of the Ceratocystidaceae.     

Starter unit specificity directs genome mining of polyketide synthase pathways in fungi

Search of the protein database with the aflatoxin pathway polyketide synthase (PKS) revealed putative PKSs in the pathogenic fungi Coccidioides immitis and Coccidioides posadasii that could require partnerships with a pair of fatty acid synthase (FAS) subunits for the biosynthesis of fatty acid-polyketide hybrid metabolites. A starter unit:acyl-carrier protein transacylase (SAT) domain was discovered in the nonreducing PKS. This domain is thought to accept the fatty acid product from the FAS to initiate polyketide synthesis. We expressed the C. immitis SAT domain in Escherichia coli and showed that this domain, unlike that from the aflatoxin pathway PKS, transferred octanoyl-CoA four times faster than hexanoyl-CoA. The SAT domain also formed a covalent octanoyl intermediate and transferred this group to a free-standing ACP domain. Our results suggest that C. immitis/posadasii, both human fungal pathogens, contain a FAS/PKS cluster with functional similarity to the aflatoxin cluster found in Aspergillus species. Dissection of the PKS and determination of in vitro SAT domain specificity provides a tool to uncover the growing number of similar sequenced pathways in fungi, and to guide elucidation of the fatty acid-polyketide hybrid metabolites that they produce.     

Approaches to stabilization of inter-domain recombination in polyketide synthase gene expression plasmids

Regions of extremely high sequence identity are recurrent in modular polyketide synthase (PKS) genes. Such sequences are potentially detrimental to the stability of PKS expression plasmids used in the combinatorial biosynthesis of polyketide metabolites. We present two different solutions for circumventing intra-plasmid recombination within the megalomicin PKS genes in Streptomyces coelicolor. In one example, a synthetic gene was used in which the codon usage was reengineered without affecting the primary amino acid sequence. The other approach utilized a heterologous subunit complementation strategy to replace one of the problematic regions. Both methods resulted in PKS complexes capable of 6-deoxyerythronolide B analogue biosynthesis in S. coelicolor CH999, permitting reproducible scale-up to at least 5-l stirred-tank fermentation and a comparison of diketide precursor incorporation efficiencies between the erythromycin and megalomicin PKSs. Electronic Publication     

Advances in cloning, functional analysis and heterologous expression of fungal polyketide synthase genes

Fungal polyketides comprise a diverse group of secondary metabolites that play an important role for drug discovery, as pigments, and as mycotoxins. Their biosynthesis is governed by multidomain enzymes, so-called fungal type I polyketide synthases (PKS). Investigating the molecular basis of polyketide biosynthesis in fungi is of great importance for ecological and pharmacological reasons. In addition, cloning, functional analysis and expression of fungal PKS genes also set the basis for engineering the yet largely untapped biosynthetic potential.     

Characterization of two novel non-reducing polyketide synthase genes from the lichen-forming fungus Hypogymnia physodes

Lichens are known to produce a variety of secondary metabolites including polyketides, which have valuable biological activities. Some polyketides are produced solely by lichens. The biosynthesis of these compounds is primarily governed by iterative type I polyketide synthases. Hypogymnia physodes synthesize polyketides such as physodic, physodalic and hydroxyphysodic acid and atranorin, which are non-reducing polyketides. Two novel non-reducing polyketide synthase (PKS) genes were isolated from a fosmid genomic library of a mycobiont of H. physodes using a 409bp fragment corresponding to part of the reductase (R) domain as a probe. H. physodes PKS1 (Hyopks1) and PKS2 (Hypopks2) contain keto synthase (KS), acyl transferase (AT), acyl carrier protein (ACP), methyl transferase (ME) and R domains. Classification based on phylogeny analysis using the translated KS and AT domains demonstrated that Hypopks1 and Hypopks2 are members of the fungal non-reducing PKSs clade III. This is the first report of non-reducing PKSs containing the R domain-mediated release mechanisms in lichens, which are also rare fungal type I PKS in non-lichenized filamentous fungi.     

Hybrid peptide-polyketide natural products: biosynthesis and prospects toward engineering novel molecules

The structural and catalytic similarities between modular nonribosomal peptide synthetase (NRPS) and polyketide synthase (PKS) inspired us to search for hybrid NRPS-PKS systems. By examining the biochemical and genetic data known to date for the biosynthesis of hybrid peptide-polyketide natural products, we show (1) that the same catalytic sites are conserved between the hybrid NRPS-PKS and normal NRPS or PKS systems, although the ketoacyl synthase domain in NRPS/PKS hybrids is unique, and (2) that specific interpolypeptide linkers exist at both the C- and N-termini of the NRPS and PKS proteins, which presumably play a critical role in facilitating the transfer of the growing peptide or polyketide intermediate between NRPS and PKS modules in hybrid NRPS-PKS systems. These findings provide new insights for intermodular communications in hybrid NRPS-PKS systems and should now be taken into consideration in engineering hybrid peptide-polyketide biosynthetic pathways for making novel "unnatural" natural products.     

Engineered biosynthesis of regioselectively modified aromatic polyketides using bimodular polyketide synthases

Bacterial aromatic polyketides such as tetracycline and doxorubicin are a medicinally important class of natural products produced as secondary metabolites by actinomyces bacteria. Their backbones are derived from malonyl-CoA units by polyketide synthases (PKSs). The nascent polyketide chain is synthesized by the minimal PKS, a module consisting of four dissociated enzymes. Although the biosynthesis of most aromatic polyketide backbones is initiated through decarboxylation of a malonyl building block (which results in an acetate group), some polyketides, such as the estrogen receptor antagonist R1128, are derived from nonacetate primers. Understanding the mechanism of nonacetate priming can lead to biosynthesis of novel polyketides that have improved pharmacological properties. Recent biochemical analysis has shown that nonacetate priming is the result of stepwise activity of two dissociated PKS modules with orthogonal molecular recognition features. In these PKSs, an initiation module that synthesizes a starter unit is present in addition to the minimal PKS module. Here we describe a general method for the engineered biosynthesis of regioselectively modified aromatic polyketides. When coexpressed with the R1128 initiation module, the actinorhodin minimal PKS produced novel hexaketides with propionyl and isobutyryl primer units. Analogous octaketides could be synthesized by combining the tetracenomycin minimal PKS with the R1128 initiation module. Tailoring enzymes such as ketoreductases and cyclases were able to process the unnatural polyketides efficiently. Based upon these findings, hybrid PKSs were engineered to synthesize new anthraquinone antibiotics with predictable functional group modifications. Our results demonstrate that (i) bimodular aromatic PKSs present a general mechanism for priming aromatic polyketide backbones with nonacetate precursors; (ii) the minimal PKS controls polyketide chain length by counting the number of atoms incorporated into the backbone rather than the number of elongation cycles; and (iii) in contrast, auxiliary PKS enzymes such as ketoreductases, aromatases, and cyclases recognize specific functional groups in the backbone rather than overall chain length. Among the anthracyclines engineered in this study were compounds with (i) more superior activity than R1128 against the breast cancer cell line MCF-7 and (ii) inhibitory activity against glucose-6-phosphate translocase, an attractive target for the treatment of Type II diabetes.     

Characterization of Two Polyketide Synthase Genes Involved in Zearalenone Biosynthesis in Gibberella zeae

Zearalenone, a mycotoxin produced by several Fusarium spp., is most commonly found as a contaminant in stored grain and has chronic estrogenic effects on mammals. Zearalenone is a polyketide derived from the sequential condensation of multiple acetate units by a polyketide synthase (PKS), but the genetics of its biosynthesis are not understood. We cloned two genes, designated ZEA1 and ZEA2, which encode polyketide synthases that participate in the biosynthesis of zearalenone by Gibberella zeae (anamorph Fusarium graminearum). Disruption of either gene resulted in the loss of zearalenone production under inducing conditions. ZEA1 and ZEA2 are transcribed divergently from a common promoter region. Quantitative PCR analysis of both PKS genes and six flanking genes supports the view that the two polyketide synthases make up the core biosynthetic unit for zearalenone biosynthesis. An appreciation of the genetics of zearalenone biosynthesis is needed to understand how zearalenone is synthesized under field conditions that result in the contamination of grain.     

Novel intermolecular iterative mechanism for biosynthesis of mycoketide catalyzed by a bimodular polyketide synthase

In recent years, remarkable versatility of polyketide synthases (PKSs) has been recognized; both in terms of their structural and functional organization as well as their ability to produce compounds other than typical secondary metabolites. Multifunctional Type I PKSs catalyze the biosynthesis of polyketide products by either using the same active sites repetitively (iterative) or by using these catalytic domains only once (modular) during the entire biosynthetic process. The largest open reading frame in Mycobacterium tuberculosis, pks12, was recently proposed to be involved in the biosynthesis of mannosyl-beta-1-phosphomycoketide (MPM). The PKS12 protein contains two complete sets of modules and has been suggested to synthesize mycoketide by five alternating condensations of methylmalonyl and malonyl units by using an iterative mode of catalysis. The bimodular iterative catalysis would require transfer of intermediate chains from acyl carrier protein domain of module 2 to ketosynthase domain of module 1. Such bimodular iterations during PKS biosynthesis have not been characterized and appear unlikely based on recent understanding of the three-dimensional organization of these proteins. Moreover, all known examples of iterative PKSs so far characterized involve unimodular iterations. Based on cell-free reconstitution of PKS12 enzymatic machinery, in this study, we provide the first evidence for a novel "modularly iterative" mechanism of biosynthesis. By combination of biochemical, computational, mutagenic, analytical ultracentrifugation and atomic force microscopy studies, we propose that PKS12 protein is organized as a large supramolecular assembly mediated through specific interactions between the C- and N-terminus linkers. PKS12 protein thus forms a modular assembly to perform repetitive condensations analogous to iterative proteins. This novel intermolecular iterative biosynthetic mechanism provides new perspective to our understanding of polyketide biosynthetic machinery and also suggests new ways to engineer polyketide metabolites. The characterization of novel molecular mechanisms involved in biosynthesis of mycobacterial virulent lipids has opened new avenues for drug discovery.     

Evolution of ketosynthase domains of polyketide synthase genes in the Cladonia chlorophaea species complex (Cladoniaceae)

Lichen-forming fungi synthesize a diversity of polyketides, but only a few non-reducing polyketide synthase (PKS) genes from a lichen-forming fungus have been linked with a specific polyketide. While it is a challenge to link the large number of PKS paralogs in fungi with specific products, it might be expected that the PKS paralogs from closely related species would be similar because of recent evolutionary divergence. The objectives of this study were to reconstruct a PKS gene phylogeny of the Cladonia chlorophaea species complex based on the ketosynthase domain, a species phylogeny of the complex, and to explore the presence of PKS gene paralogs among members of the species complex. DNA was isolated from 51 individuals of C. chlorophaea and allies to screen for the presence of 13 PKS paralogs. A 128 sequence PKS gene phylogeny using deduced amino acid sequences estimated from the 13 PKS paralogs and sequences subjected to BLASTx comparisons showed losses of each of two PKS domains (reducing and methylation). This research provided insight into the evolution of PKS genes in the C. chlorophaea group, species evolution in the group, and it identified potential directions for further investigation of polyketide synthesis in the C. chlorophaea species complex.     

Natural biocombinatorics in the polyketide synthase genes of the actinobacterium Streptomyces avermitilis

Modular polyketide synthases (PKSs) of bacteria provide an enormous reservoir of natural chemical diversity. Studying natural biocombinatorics may aid in the development of concepts for experimental design of genes for the biosynthesis of new bioactive compounds. Here we address the question of how the modularity of biosynthetic enzymes and the prevalence of multiple gene clusters in Streptomyces drive the evolution of metabolic diversity. The phylogeny of ketosynthase (KS) domains of Streptomyces PKSs revealed that the majority of modules involved in the biosynthesis of a single compound evolved by duplication of a single ancestor module. Using Streptomyces avermitilis as a model organism, we have reconstructed the evolutionary relationships of different domain types. This analysis suggests that 65% of the modules were altered by recombinational replacements that occurred within and between biosynthetic gene clusters. The natural reprogramming of the biosynthetic pathways was unambiguously confined to domains that account for the structural diversity of the polyketide products and never observed for the KS domains. We provide examples for natural acyltransferase (AT), ketoreductase (KR), and dehydratase (DH)–KR domain replacements. Potential sites of homologous recombination could be identified in interdomain regions and within domains. Our results indicate that homologous recombination facilitated by the modularity of PKS architecture is the most important mechanism underlying polyketide diversity in bacteria.     

A gene cluster involved in nogalamycin biosynthesis fromStreptomyces nogalater: sequence analysis and complementation of early-block mutations in the anthracycline pathway

We have analyzed an anthracycline biosynthesis gene cluster fromStreptomyces nogalater. Based on sequence analysis, a contiguous region of 11 kb is deduced to include genes for the early steps in anthracycline biosynthesis, a regulatory gene (snoA) promoting the expression of the biosynthetic genes, and at least one gene whose product might have a role in modification of the glycoside moiety. The three ORFs encoding a minimal polyketide synthase (PKS) are separated from the regulatory gene (snoA) by a comparatively AT-rich region (GC content 60%). Subfragments of the DNA region were transferred toStreptomyces galilaeus mutants blocked in aclacinomycin biosynthesis, and to a regulatory mutant ofS. nogalater. TheS. galilaeus mutants carrying theS. nogalater minimal PKS genes produced auramycinone glycosides, demonstrating replacement of the starter unit for polyketide biosynthesis. The product ofsnoA seems to be needed for expression of at least the genes for the minimal PKS.     

A novel tunnel in mycobacterial type III polyketide synthase reveals the structural basis for generating diverse metabolites

The superfamily of plant and bacterial type III polyketide synthases (PKSs) produces diverse metabolites with distinct biological functions. PKS18, a type III PKS from Mycobacterium tuberculosis, displays an unusual broad specificity for aliphatic long-chain acyl-coenzyme A (acyl-CoA) starter units (C(6)-C(20)) to produce tri- and tetraketide pyrones. The crystal structure of PKS18 reveals a 20 A substrate binding tunnel, hitherto unidentified in this superfamily of enzymes. This remarkable tunnel extends from the active site to the surface of the protein and is primarily generated by subtle changes of backbone dihedral angles in the core of the protein. Mutagenic studies combined with structure determination provide molecular insights into the structural elements that contribute to the chain length specificity of the enzyme. This first bacterial type III PKS structure underlines a fascinating example of the way in which subtle changes in protein architecture can generate metabolite diversity in nature.     

Non-modular polyketide synthases in myxobacteria

Myxobacteria are prolific producers of a wide variety of secondary metabolites. The vast majority of these compounds are complex polyketides which are biosynthesised by multimodular polyketide synthases (PKSs). In contrast, few myxobacterial metabolites isolated to date are derived from non-modular PKSs, in particular type III PKSs. This review reports our progress on the characterisation of type III PKSs in myxobacteria. We also summarize current knowledge on bacterial type III PKSs, with a special focus on the evolutionary relationship between plant and bacterial enzymes. The biosynthesis of a quinoline alkaloid in Stigmatella aurantiaca by a non-modular PKS is also discussed.