Chemoenzymatic and template-directed synthesis of bioactive macrocyclic peptides.

Non-ribosomally synthesized peptides have compelling biological activities ranging from antimicrobial to immunosuppressive and from cytostatic to antitumor. The broad spectrum of applications in modern medicine is reflected in the great structural diversity of these natural products. They contain unique building blocks, such as d-amino acids, fatty acids, sugar moieties, and heterocyclic elements, as well as halogenated, methylated, and formylated residues. In the past decades, significant progress has been made toward the understanding of the biosynthesis of these secondary metabolites by nonribosomal peptide synthetases (NRPSs) and their associated tailoring enzymes. Guided by this knowledge, researchers genetically redesigned the NRPS template to synthesize new peptide products. Moreover, chemoenzymatic strategies were developed to rationally engineer nonribosomal peptides products in order to increase or alter their bioactivities. Specifically, chemical synthesis combined with peptide cyclization mediated by nonribosomal thioesterase domains enabled the synthesis of glycosylated cyclopeptides, inhibitors of integrin receptors, peptide/polyketide hybrids, lipopeptide antibiotics, and streptogramin B antibiotics. In addition to the synthetic potential of these cyclization catalysts, which is the main focus of this review, different enzymes for tailoring of peptide scaffolds as well as the manipulation of carrier proteins with reporter-labeled coenzyme A analogs are discussed.     

A proteomic survey of nonribosomal peptide and polyketide biosynthesis in actinobacteria

Actinobacteria such as streptomycetes are renowned for their ability to produce bioactive natural products including nonribosomal peptides (NRPs) and polyketides (PKs). The advent of genome sequencing has revealed an even larger genetic repertoire for secondary metabolism with most of the small molecule products of these gene clusters still unknown. Here, we employed a "protein-first" method called PrISM (Proteomic Investigation of Secondary Metabolism) to screen 26 unsequenced actinomycetes using mass spectrometry-based proteomics for the targeted detection of expressed nonribosomal peptide synthetases or polyketide synthases. Improvements to the original PrISM screening approach (Nat. Biotechnol. 2009, 27, 951-956), for example, improved de novo peptide sequencing, have enabled the discovery of 10 NRPS/PKS gene clusters from 6 strains. Taking advantage of the concurrence of biosynthetic enzymes and the secondary metabolites they generate, two natural products were associated with their previously "orphan" gene clusters. This work has demonstrated the feasibility of a proteomics-based strategy for use in screening for NRP/PK production in actinomycetes (often >8 Mbp, high GC genomes) versus the bacilli (2-4 Mbp genomes) used previously.     

Synthetic biology,genome mining,and combinatorial biosynthesis of NRPS-derived antibiotics: a perspective

Combinatorial biosynthesis of novel secondary metabolites derived from nonribosomal peptide synthetases (NRPSs) has been in slow development for about a quarter of a century. Progress has been hampered by the complexity of the giant multimodular multienzymes. More recently, advances have been made on understanding the chemical and structural biology of these complex megaenzymes, and on learning the design rules for engineering functional hybrid enzymes. In this perspective, I address what has been learned about successful engineering of complex lipopeptides related to daptomycin, and discuss how synthetic biology and microbial genome mining can converge to broaden the scope and enhance the speed and robustness of combinatorial biosynthesis of NRPS-derived natural products for drug discovery.     

Exploring the biosynthesis of natural products and their inherent suitability for the rational design of desirable compounds through genetic engineering

Much effort has been invested in studying how natural products are biosynthesized, and great advances have been made in understanding how these compounds acquire their structural complexity and biological activities. In recent years, significant progress has been made due to the devoted efforts of scientists in this field and to technological advancements. Numerous details, applications, and innovative findings have been elucidated by scientists using biochemical, genetic, and molecular biological approaches. Here I present a comprehensive overview of highly valued biosynthetic proteins, polyketide synthase and nonribosomal peptide synthetase. I begin with "Diels-Alderase" a captivating enzyme that has the ability to catalyze a Diels-Alder reaction valued by chemists for its usefulness in chemical synthesis. A handful of these enzymes have been characterized and chemically authenticated. The most well understood enzyme of this category is macrophomate synthase. Secondly, I focus on the polyketide and nonribosomal peptide biosynthetic pathways and the enzyme assembly for producing its metabolite. Many important natural products are produced by this biosynthetic pathway as secondary metabolites, such as erythromycin, rifamycin, and FK520, as antibiotics and an immunosuppressive, respectively. I conclude with a discussion of nonribosomal peptides and their mechanistic pathways. Special attention will be devoted to de novo production of echinomycin in a heterologous manner, the earliest example of totally engineered biosynthesis of the biologically active form of a nonribosomal peptide host in Escherichia coli.     

Molecular biology of peptide and polyketide biosynthesis in cyanobacteria

Cyanobacteria produce numerous and structurally diverse secondary metabolites, in particular nonribosomal peptide and polyketide structures. Various bioactivities could be assigned to these compounds, and some may prove useful either for development into commercial drugs or as biochemical research tools. Microcystin, a worldwide common cyanobacterial hepatotoxin, was the first metabolite whose nonribosomal biosynthesis could be confirmed by knock-out mutagenesis. The microcystin synthetase complex consists of peptide synthetases, polyketide synthases, and hybrid enzymes, and reveals a number of novel enzymatic features, signifying the potential of cyanobacterial biosynthetic systems for combinatorial biochemistry. Recent studies have shown the presence of peptide synthetase genes and polyketide synthase genes within a number of cyanobacterial genomes. This knowledge may be very valuable for future screening projects aimed at the detection of new bioactive compounds.     

Genomic and transcriptomic analysis of the endophytic fungus Pestalotiopsis fici reveals its lifestyle and high potential for synthesis of natural products

Background

In recent years, the genus Pestalotiopsis is receiving increasing attention, not only because of its economic impact as a plant pathogen but also as a commonly isolated endophyte which is an important source of bioactive natural products. Pestalotiopsis fici Steyaert W106-1/CGMCC3.15140 as an endophyte of tea produces numerous novel secondary metabolites, including chloropupukeananin, a derivative of chlorinated pupukeanane that is first discovered in fungi. Some of them might be important as the drug leads for future pharmaceutics.

Results

Here, we report the genome sequence of the endophytic fungus of tea Pestalotiopsis fici W106-1/CGMCC3.15140. The abundant carbohydrate-active enzymes especially significantly expanding pectinases allow the fungus to utilize the limited intercellular nutrients within the host plants, suggesting adaptation of the fungus to endophytic lifestyle. The P. fici genome encodes a rich set of secondary metabolite synthesis genes, including 27 polyketide synthases (PKSs), 12 non-ribosomal peptide synthases (NRPSs), five dimethylallyl tryptophan synthases, four putative PKS-like enzymes, 15 putative NRPS-like enzymes, 15 terpenoid synthases, seven terpenoid cyclases, seven fatty-acid synthases, and five hybrids of PKS-NRPS. The majority of these core enzymes distributed into 74 secondary metabolite clusters. The putative Diels-Alderase genes have undergone expansion.

Conclusion

The significant expansion of pectinase encoding genes provides essential insight in the life strategy of endophytes, and richness of gene clusters for secondary metabolites reveals high potential of natural products of endophytic fungi.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-014-1190-9) contains supplementary material, which is available to authorized users.     

Single cell genome amplification accelerates identification of the apratoxin biosynthetic pathway from a complex microbial assemblage

Filamentous marine cyanobacteria are extraordinarily rich sources of structurally novel, biomedically relevant natural products. To understand their biosynthetic origins as well as produce increased supplies and analog molecules, access to the clustered biosynthetic genes that encode for the assembly enzymes is necessary. Complicating these efforts is the universal presence of heterotrophic bacteria in the cell wall and sheath material of cyanobacteria obtained from the environment and those grown in uni-cyanobacterial culture. Moreover, the high similarity in genetic elements across disparate secondary metabolite biosynthetic pathways renders imprecise current gene cluster targeting strategies and contributes sequence complexity resulting in partial genome coverage. Thus, it was necessary to use a dual-method approach of single-cell genomic sequencing based on multiple displacement amplification (MDA) and metagenomic library screening. Here, we report the identification of the putative apratoxin. A biosynthetic gene cluster, a potent cancer cell cytotoxin with promise for medicinal applications. The roughly 58 kb biosynthetic gene cluster is composed of 12 open reading frames and has a type I modular mixed polyketide synthase/nonribosomal peptide synthetase (PKS/NRPS) organization and features loading and off-loading domain architecture never previously described. Moreover, this work represents the first successful isolation of a complete biosynthetic gene cluster from Lyngbya bouillonii, a tropical marine cyanobacterium renowned for its production of diverse bioactive secondary metabolites.     

ASMPKS: an analysis system for modular polyketide synthases

Background  

Polyketides are secondary metabolites of microorganisms with diverse biological activities, including pharmacological functions such as antibiotic, antitumor and agrochemical properties. Polyketides are synthesized by serialized reactions of a set of enzymes called polyketide synthase(PKS)s, which coordinate the elongation of carbon skeletons by the stepwise condensation of short carbon precursors. Due to their importance as drugs, the volume of data on polyketides is rapidly increasing and creating a need for computational analysis methods for efficient polyketide research. Moreover, the increasing use of genetic engineering to research new kinds of polyketides requires genome wide analysis.     

Identification of a biosynthetic gene cluster and the six associated lipopeptides involved in swarming motility of Pseudomonas syringae pv. tomato DC3000

Pseudomonas species are known to be prolific producers of secondary metabolites that are synthesized wholly or in part by nonribosomal peptide synthetases. In an effort to identify additional nonribosomal peptides produced by these bacteria, a bioinformatics approach was used to "mine" the genome of Pseudomonas syringae pv. tomato DC3000 for the metabolic potential to biosynthesize previously unknown nonribosomal peptides. Herein we describe the identification of a nonribosomal peptide biosynthetic gene cluster that codes for proteins involved in the production of six structurally related linear lipopeptides. Structures for each of these lipopeptides were proposed based on amino acid analysis and mass spectrometry analyses. Mutations in this cluster resulted in the loss of swarming motility of P. syringae pv. tomato DC3000 on medium containing a low percentage of agar. This phenotype is consistent with the loss of the ability to produce a lipopeptide that functions as a biosurfactant. This work gives additional evidence that mining the genomes of microorganisms followed by metabolite and phenotypic analyses leads to the identification of previously unknown secondary metabolites.     

Practical 4'-phosphopantetheine active site discovery from proteomic samples

Polyketide and nonribosomal peptides constitute important classes of small molecule natural products. Due to the proven biological activities of these compounds, novel methods for discovery and study of the polyketide synthase (PKS) and nonribosomal peptide synthetase (NRPS) enzymes responsible for their production remains an area of intense interest, and proteomic approaches represent a relatively unexplored avenue. While these enzymes may be distinguished from the proteomic milieu by their use of the 4'-phosphopantetheine (PPant) post-translational modification, proteomic detection of PPant peptides is hindered by their low abundance and labile nature which leaves them unassigned using traditional database searching. Here we address key experimental and computational challenges to facilitate practical discovery of this important post-translational modification during shotgun proteomics analysis using low-resolution ion-trap mass spectrometers. Activity-based enrichment maximizes MS input of PKS/NRPS peptides, while targeted fragmentation detects putative PPant active sites. An improved data analysis pipeline allows experimental identification and validation of these PPant peptides directly from MS2 data. Finally, a machine learning approach is developed to directly detect PPant peptides from only MS2 fragmentation data. By providing new methods for analysis of an often cryptic post-translational modification, these methods represent a first step toward the study of natural product biosynthesis in proteomic settings.     

Expansion of RiPP biosynthetic space through integration of pan-genomics and machine learning uncovers a novel class of lanthipeptides

Microbial natural products constitute a wide variety of chemical compounds, many which can have antibiotic, antiviral, or anticancer properties that make them interesting for clinical purposes. Natural product classes include polyketides (PKs), nonribosomal peptides (NRPs), and ribosomally synthesized and post-translationally modified peptides (RiPPs). While variants of biosynthetic gene clusters (BGCs) for known classes of natural products are easy to identify in genome sequences, BGCs for new compound classes escape attention. In particular, evidence is accumulating that for RiPPs, subclasses known thus far may only represent the tip of an iceberg. Here, we present decRiPPter (Data-driven Exploratory Class-independent RiPP TrackER), a RiPP genome mining algorithm aimed at the discovery of novel RiPP classes. DecRiPPter combines a Support Vector Machine (SVM) that identifies candidate RiPP precursors with pan-genomic analyses to identify which of these are encoded within operon-like structures that are part of the accessory genome of a genus. Subsequently, it prioritizes such regions based on the presence of new enzymology and based on patterns of gene cluster and precursor peptide conservation across species. We then applied decRiPPter to mine 1,295 Streptomyces genomes, which led to the identification of 42 new candidate RiPP families that could not be found by existing programs. One of these was studied further and elucidated as a representative of a novel subfamily of lanthipeptides, which we designate class V. The 2D structure of the new RiPP, which we name pristinin A3 (1), was solved using nuclear magnetic resonance (NMR), tandem mass spectrometry (MS/MS) data, and chemical labeling. Two previously unidentified modifying enzymes are proposed to create the hallmark lanthionine bridges. Taken together, our work highlights how novel natural product families can be discovered by methods going beyond sequence similarity searches to integrate multiple pathway discovery criteria.

This study shows that decRiPPter, an innovative algorithmic approach using pan-genomics and machine learning, can discover novel types of ribosomally synthesized peptide (RIPP) natural products, including a new class of lanthipeptides.     

Chemoenzymatic synthesis of cryptophycin/arenastatin natural products

Microbially derived modular polyketide synthase and nonribosomal peptide synthetase biosynthetic pathways are a rich source of novel natural products. Development of these systems for the engineered biosynthesis of diverse secondary metabolites continues to progress as a robust source of chemical diversity. Recent efforts that employ individual enzymes and catalytic domains for the production or modification of small molecules have met with growing success. In this study, the thioesterase domain from the cryptophycin biosynthetic pathway was isolated and its function evaluated with a series of linear chain elongation intermediates in developing a novel chemoenzymatic synthesis of the cryptophycin/arenastatin class of antitumor agents. The results show the high efficiency of the thioesterase in generating the 16-membered depsipeptide ring of this important natural product system. Moreover, analysis of selected substrates revealed considerable tolerance for structural variation within the seco-cryptophycin unit C beta-alanine residue, but strict structural requirements at the phenyl group position of the unit A delta-hydroxy octadienoate chain elongation intermediates.     

Function of MbtH homologs in nonribosomal peptide biosynthesis and applications in secondary metabolite discovery

Mycobacterium tuberculosis encodes mycobactin, a peptide siderophore that is biosynthesized by a nonribosomal peptide synthetase (NRPS) mechanism. Within the mycobactin biosynthetic gene cluster is a gene that encodes a 71-amino-acid protein MbtH. Many other NRPS gene clusters harbor mbtH homologs, and recent genetic, biochemical, and structural studies have begun to shed light on the function(s) of these proteins. In some cases, MbtH-like proteins are required for biosynthesis of their cognate peptides, and non-cognate MbtH-like proteins have been shown to be partially complementary. Biochemical studies revealed that certain MbtH-like proteins participate in tight binding to NRPS proteins containing adenylation (A) domains where they stimulate adenylation reactions. Expression of MbtH-like proteins is important for a number of applications, including optimal production of native and genetically engineered secondary metabolites produced by mechanisms that employ NRPS enzymes. They also may serve as beacons to identify gifted actinomycetes and possibly other bacteria that encode multiple functional NRPS pathways for discovery of novel secondary metabolites by genome mining.     

Genome mining for natural product biosynthetic gene clusters in the Subsection V cyanobacteria

Background

Cyanobacteria are well known for the production of a range of secondary metabolites. Whilst recent genome sequencing projects has led to an increase in the number of publically available cyanobacterial genomes, the secondary metabolite potential of many of these organisms remains elusive. Our study focused on the 11 publically available Subsection V cyanobacterial genomes, together with the draft genomes of Westiella intricata UH strain HT-29-1 and Hapalosiphon welwitschii UH strain IC-52-3, for their genetic potential to produce secondary metabolites. The Subsection V cyanobacterial genomes analysed in this study are reported to produce a diverse range of natural products, including the hapalindole-family of compounds, microcystin, hapalosin, mycosporine-like amino acids and hydrocarbons.

Results

A putative gene cluster for the cyclic depsipeptide hapalosin, known to reverse P-glycoprotein multiple drug resistance, was identified within three Subsection V cyanobacterial genomes, including the producing cyanobacterium H. welwitschii UH strain IC-52-3. A number of orphan NRPS/PKS gene clusters and ribosomally-synthesised and post translationally-modified peptide gene clusters (including cyanobactin, microviridin and bacteriocin gene clusters) were identified. Furthermore, gene clusters encoding the biosynthesis of mycosporine-like amino acids, scytonemin, hydrocarbons and terpenes were also identified and compared.

Conclusions

Genome mining has revealed the diversity, abundance and complex nature of the secondary metabolite potential of the Subsection V cyanobacteria. This bioinformatic study has identified novel biosynthetic enzymes which have not been associated with gene clusters of known classes of natural products, suggesting that these cyanobacteria potentially produce structurally novel secondary metabolites.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-015-1855-z) contains supplementary material, which is available to authorized users.     

A 3-hydroxy-3-methylglutaryl-CoA synthase-based probe for the discovery of the acyltransferase-less type I polyketide synthases

Acyltransferase (AT)-less type I polyketide synthases (PKSs) produce complex natural products due to the presence of many unique tailoring enzymes. The 3-hydroxy-3-methylglutaryl coenzyme A synthases (HCSs) are responsible for β-alkylation of the growing polyketide intermediates in AT-less type I PKSs. In this study, we discovered a large group of HCSs, closely associated with the characterized and orphan AT-less type I PKSs through in silico genome mining, sequence and genome neighbourhood network analyses. Using HCS-based probes, the survey of 1207 in-house strains and 18 soil samples from different geographic locations revealed the vast diversity of HCS-containing AT-less type I PKSs. The presence of HCSs in many AT-less type I PKSs suggests their co-evolutionary relationship. This study provides a new probe to study the abundance and diversity of AT-less type I PKSs in the environment and microbial strain collections. Our study should inspire future efforts to discover new polyketide natural products from AT-less type I PKSs.     

Non-ribosomal peptide synthetases: Identifying the cryptic gene clusters and decoding the natural product

Non-ribosomal peptide synthetases (NRPSs) and polyketide synthases (PKSs) present in bacteria and fungi are the major multi-modular enzyme complexes which synthesize secondary metabolites like the pharmacologically important antibiotics and siderophores. Each of the multiple modules of an NRPS activates a different amino or aryl acid, followed by their condensation to synthesize a linear or cyclic natural product. The studies on NRPS domains, the knowledge of their gene cluster architecture and tailoring enzymes have helped in the in silico genetic screening of the ever-expanding sequenced microbial genomic data for the identification of novel NRPS/PKS clusters and thus deciphering novel non-ribosomal peptides (NRPs). Adenylation domain is an integral part of the NRPSs and is the substrate selecting unit for the final assembled NRP. In some cases, it also requires a small protein, the MbtH homolog, for its optimum activity. The presence of putative adenylation domain and MbtH homologs in a sequenced genome can help identify the novel secondary metabolite producers. The role of the adenylation domain in the NRPS gene clusters and its characterization as a tool for the discovery of novel cryptic NRPS gene clusters are discussed.     

Proteome analysis of Myxococcus xanthus by off-line two-dimensional chromatographic separation using monolithic poly-(styrene-divinylbenzene) columns combined with ion-trap tandem mass spectrometry

Myxobacteria are potent producers of secondary metabolites exhibiting diverse biological activities and pharmacological potential. The proteome of Myxococcus xanthus DK1622 was characterized by two-dimensional chromatographic separation of tryptic peptides from a lysate followed by tandem mass spectrometric identification. The high degree of orthogonality of the separation system employing polymer-based strong cation-exchange and monolithic reversed-phase stationary phases was clearly demonstrated. Upon automated database searching, 1312 unique peptides were identified, which were associated with 631 unique proteins. High-molecular polyketide synthetases and nonribosomal peptide synthetases, known to be involved in the biosynthesis of various secondary metabolites, were readily detected. Besides the identification of gene products associated with the production of known secondary metabolites, proteins could also be identified for six gene clusters, for which no biosynthetic product has been known so far.     

Characterization of the Nodularin Synthetase Gene Cluster and Proposed Theory of the Evolution of Cyanobacterial Hepatotoxins

Nodularia spumigena is a bloom-forming cyanobacterium which produces the hepatotoxin nodularin. The complete gene cluster encoding the enzymatic machinery required for the biosynthesis of nodularin in N. spumigena strain NSOR10 was sequenced and characterized. The 48-kb gene cluster consists of nine open reading frames (ORFs), ndaA to ndaI, which are transcribed from a bidirectional regulatory promoter region and encode nonribosomal peptide synthetase modules, polyketide synthase modules, and tailoring enzymes. The ORFs flanking the nda gene cluster in the genome of N. spumigena strain NSOR10 were identified, and one of them was found to encode a protein with homology to previously characterized transposases. Putative transposases are also associated with the structurally related microcystin synthetase (mcy) gene clusters derived from three cyanobacterial strains, indicating a possible mechanism for the distribution of these biosynthetic gene clusters between various cyanobacterial genera. We propose an alternative hypothesis for hepatotoxin evolution in cyanobacteria based on the results of comparative and phylogenetic analyses of the nda and mcy gene clusters. These analyses suggested that nodularin synthetase evolved from a microcystin synthetase progenitor. The identification of the nodularin biosynthetic gene cluster and evolution of hepatotoxicity in cyanobacteria reported in this study may be valuable for future studies on toxic cyanobacterial bloom formation. In addition, an appreciation of the natural evolution of nonribosomal biosynthetic pathways will be vital for future combinatorial engineering and rational design of novel metabolites and pharmaceuticals.     

Characterization of the nodularin synthetase gene cluster and proposed theory of the evolution of cyanobacterial hepatotoxins

Nodularia spumigena is a bloom-forming cyanobacterium which produces the hepatotoxin nodularin. The complete gene cluster encoding the enzymatic machinery required for the biosynthesis of nodularin in N. spumigena strain NSOR10 was sequenced and characterized. The 48-kb gene cluster consists of nine open reading frames (ORFs), ndaA to ndaI, which are transcribed from a bidirectional regulatory promoter region and encode nonribosomal peptide synthetase modules, polyketide synthase modules, and tailoring enzymes. The ORFs flanking the nda gene cluster in the genome of N. spumigena strain NSOR10 were identified, and one of them was found to encode a protein with homology to previously characterized transposases. Putative transposases are also associated with the structurally related microcystin synthetase (mcy) gene clusters derived from three cyanobacterial strains, indicating a possible mechanism for the distribution of these biosynthetic gene clusters between various cyanobacterial genera. We propose an alternative hypothesis for hepatotoxin evolution in cyanobacteria based on the results of comparative and phylogenetic analyses of the nda and mcy gene clusters. These analyses suggested that nodularin synthetase evolved from a microcystin synthetase progenitor. The identification of the nodularin biosynthetic gene cluster and evolution of hepatotoxicity in cyanobacteria reported in this study may be valuable for future studies on toxic cyanobacterial bloom formation. In addition, an appreciation of the natural evolution of nonribosomal biosynthetic pathways will be vital for future combinatorial engineering and rational design of novel metabolites and pharmaceuticals.     

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.