Recombinatorial biosynthesis of polyketides

Modular polyketide synthases (PKSs) from Streptomyces and related genera of bacteria produce many important pharmaceuticals. A program called CompGen was developed to carry out in silico homologous recombination between gene clusters encoding PKSs and determine whether recombinants have cluster architectures compatible with the production of polyketides. The chemical structure of recombinant polyketides was also predicted. In silico recombination was carried out for 47 well-characterised clusters. The predicted recombinants would produce 11,796 different polyketide structures. The molecular weights and average degree of reduction of the chemical structures are dispersed around the parental structures indicating that they are likely to include pharmaceutically interesting compounds. The details of the recombinants and the chemical structures were entered in a database called r-CSDB. The virtual compound library is a useful resource for computer-aided drug design and chemoinformatics strategies for finding pharmaceutically relevant chemical entities. A strategy to construct recombinant Streptomyces strains to produce these polyketides is described and the critical steps of mobilizing large biosynthetic clusters and producing new linear cloning vectors are illustrated by experimental data.     

A novel docking domain interface model predicting recombination between homoeologous modular biosynthetic gene clusters

An in silico model for homoeologous recombination between gene clusters encoding modular polyketide synthases (PKS) or non-ribosomal peptide synthetases (NRPS) was developed. This model was used to analyze recombination between 12 PKS clusters from Streptomyces species and related genera to predict if new clusters might give rise to new products. In many cases, there were only a limited number of recombination sites (about 13 per cluster pair), suggesting that recombination may pose constraints on the evolution of PKS clusters. Most recombination events occurred between pairs of ketosynthase (KS) domains, allowing the biosynthetic outcome of the recombinant modules to be predicted. About 30% of recombinants were predicted to produce polyketides. Four NRPS clusters from Streptomyces strains were also used for in silico recombination. They yielded a comparable number of recombinants to PKS clusters, but the adenylation (A) domains contained the largest proportion of recombination events; this might be a mechanism for producing new substrate specificities. The extreme G + C-content, the presence of linear chromosomes and plasmids, as well as the lack of a mutSL-mismatch repair system should favor production of recombinants in Streptomyces species.     

Characterization of the enzymatic domains in the modular polyketide synthase involved in rifamycin B biosynthesis by Amycolatopsis mediterranei

Five clustered polyketide synthase (PKS) genes, rifA–rifE, involved in rifamycin (Rf) biosynthesis in Amycolatopsis mediterranei S699 have been cloned and sequenced (August, P.R. et al., 1998. Chem. Biol. 5, 69–79). The five multifunctional polypeptides constitute a type I modular PKS that contains ten modules, each responsible for a specific round of polyketide chain elongation. Sequence comparisons of the Rf PKS proteins with other prokaryotic modular PKSs elucidated the regions that have an important role in enzyme activity and specificity. The β-ketoacyl:acyl carrier protein synthase (KS) domains show the highest degree of similarity between themselves (86–90%) and to other PKSs (78–85%) among all the constituent domains. Both malonyl-coenzyme A (MCoA) and methylmalonyl-coenzyme A (mMCoA) are substrates for chain elongation steps carried out by the Rf PKS. Since acyltransferase (AT) domains of modular PKSs can distinguish between these two substrates, comparison of the sequence of all ten AT domains of the Rf PKS with those found in the erythromycin (Er) (Donadio, S. and Katz, L., 1992. Gene 111, 51–60) and rapamycin (Rp) (Haydock, S. et al., 1995. FEBS Lett. 374, 246–248) PKSs revealed that the AT domains in module 2 of RifA and module 9 of RifE are specific for MCoA, whereas the other eight modules specify mMCoA. Dehydration of the β-hydroxyacylthioester intermediates should occur during the reactions catalysed by module 4 of RifB and modules 9 and 10 of RifE, yet only the active site region of module 4 conforms closely to the dehydratase (DH) motifs in the Er and Rp PKSs. The DH domains of modules 9 and 10 diverge significantly from the consensus sequence defined by the Er and Rp PKSs, except for the active site His residues. Deletions in the DH active sites of module 1 in RifA and module 5 in RifB and in the N- and C-terminal regions of module 8 of RifD should inactivate these domains, and module 2 of RifA lacks a DH domain, all of which are consistent with the proposed biosynthesis of Rf. In contrast, module 6 of RifB and module 7 of RifC appear to contain intact DH domains even though DH activity is not apparently required in these modules. Module 2 of RifA lacks a β-ketoacyl:acyl carrier protein reductase (KR) domain and the one in module 3 has an apparently inactive NADPH binding motif, similar to one found in the Er PKS, while the other eight KR domains of the Rf PKS should be functional. These observations are consistent with biosynthetic predictions. All the acyl carrier protein (ACP) domains, while clearly functional, nevertheless have active site signature sequences distinctive from those of the Er and Rp PKSs. Module 2 of RifA has only the core domains (KS, AT and ACP). The starter unit ligase (SUL) and ACP domains present in the N-terminus of RifA direct the selection and loading of the starter unit, 3-amino-5-hydroxybenzoic acid (AHBA), onto the PKS. AHBA is made by the products of several other genes in the Rf cluster through a variant of the shikimate pathway (August, P.R. et al., inter alia). RifF, produced by the gene immediately downstream of rifE, is thought to catalyse the intramolecular cyclization of the PKS product, thereby forming the ansamacrolide precursor of Rf B.     

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.     

From genetic diversity to metabolic unity: studies on the biosynthesis of aurafurones and aurafuron-like structures in myxobacteria and streptomycetes

The myxobacterial polyketide secondary metabolites aurafuron A and B were identified by genome mining in the myxobacterial strain Stigmatella aurantiaca DW4/3-1. The compounds contain an unusual furanone moiety and resemble metabolites isolated from soil-dwelling and marine actinobacteria, a fungus and mollusks. We describe here the cloning and functional analysis of the aurafuron biosynthetic gene cluster, including site-directed mutagenesis and feeding studies using labeled precursors. The polyketide core of the aurafurones is assembled by a modular polyketide synthase (PKS). As with many such systems described from myxobacteria, the aurafuron PKS exhibits a number of unusual features, including the apparent iterative use of a module, redundant modules and domains, a trans acting dehydratase and the absence of a terminal thioesterase domain. Four oxidoreductases are encoded within the gene locus, some of which likely participate in formation of the furanone moiety via a Baeyer-Villiger type oxidation. Indeed, inactivation of a gene encoding a cytochrome P₄₅₀ monooxygenase completely abolished production of both compounds. We also compare the complete gene locus to biosynthetic gene clusters from two Streptomyces sp., which produce close structural analogues of the aurafurones. A portion of the post-PKS biosynthetic machinery is strikingly similar in all three cases, in contrast to the PKS genes, which are highly divergent. Phylogenetic analysis of the ketosynthase domains further indicates that the PKSs have developed independently (polyphyletically) during evolution. These findings point to a currently unknown but important biological function of aurafuron-like compounds for the producing organisms.     

Cloning of modular type I polyketide synthase genes from salinomycin producing strain of Streptomyces albus

Cloning of polyether polyketide synthase (PKS) genes for salinomycin biosynthesis was attempted from Streptomyces albus. Seven beta-ketoacyl synthase (KS) core regions were obtained by PCR amplification using primers designed based on the conserved KS domains of type I PKSs. Using the KS fragment as a probe, screening of an S. albus genomic DNA library was carried out by colony hybridization. From the positive cosmid clone isolated, a 4.5-kbBamHI fragment was subcloned and sequenced. It showed high homology with bacterial type I PKSs and was deduced to code for KS, malonyl transferase, and ketoreductase motifs. By gene disruption with this 4.5-kb BamHI fragment, the cloned gene was shown to be a part of the salinomycin biosynthetic gene cluster of S. albus.     

Oxazolomycin Biosynthesis in Streptomyces albus JA3453 Featuring an “Acyltransferase-less” Type I Polyketide Synthase That Incorporates Two Distinct Extender Units

The oxazolomycins (OZMs) are a growing family of antibiotics produced by several Streptomyces species that show diverse and important antibacterial, antitumor, and anti-human immunodeficiency virus activity. Oxazolomycin A is a peptide-polyketide hybrid compound containing a unique spiro-linked β-lactone/γ-lactam, a 5-substituted oxazole ring. The oxazolomycin biosynthetic gene cluster (ozm) was identified from Streptomyces albus JA3453 and localized to 79.5-kb DNA, consisting of 20 open reading frames that encode non-ribosomal peptide synthases, polyketide synthases (PKSs), hybrid non-ribosomal peptide synthase-PKS, trans-acyltransferases (trans-ATs), enzymes for methoxymalonyl-acyl carrier protein (ACP) synthesis, putative resistance genes, and hypothetical regulation genes. In contrast to classical type I polyketide or fatty acid biosynthases, all 10 PKS modules in the gene cluster lack cognate ATs. Instead, discrete ATs OzmM (with tandem domains OzmM-AT1 and OzmM-AT2) and OzmC were equipped to carry out all of the loading functions of both malonyl-CoA and methoxymalonyl-ACP extender units. Strikingly, only OzmM-AT2 is required for OzmM activity for OZM biosynthesis, whereas OzmM-AT1 seemed to be a cryptic AT domain. The above findings, together with previous results using isotope-labeled precursor feeding assays, are assembled for the OZM biosynthesis model to be proposed. The incorporation of both malonyl-CoA (by OzmM-AT2) and methoxymalonyl-ACP (by OzmC) extender units seemed to be unprecedented for this class of trans-AT type I PKSs, which might be fruitfully manipulated to create structurally diverse novel compounds.     

Alternative method for site-directed mutagenesis of complex polyketide synthase in Streptomyces albus JA3453

Sequence analysis of oxazolomycin (OZM) biosynthetic gene cluster from Streptomyces albus JA3453 revealed a gene, ozmH , encoding a hybrid polyketide and non-ribosomal pep-tide enzyme. Tandem ketosynthase (KS) domains (KS_10–1 and KS_10–2) were characterized and they show significant homol-ogy with known KSs. Using an alternative method that involves RecA-mediated homologous recombination, the negative selection marker sacB gene, and temperature-sensitive replications, site-directed mutagenesis of the catalytic triad amino acid cysteines were carried out in each of the tandem KS domains totest the function they play in OZM biosynthesis. HPLC-mass spectrometry analysis of the resulting mutant strains showed that KS_10–2 is essential for OZM biosynthesis but KS_10–1 is not indispensable and might serve as a "redundant" domain. These results confirmed the existence of an "extra domain" in complex polyketide synthase.     

Metagenomic analysis reveals diverse polyketide synthase gene clusters in microorganisms associated with the marine sponge Discodermia dissoluta

Sponge-associated bacteria are thought to produce many novel bioactive compounds, including polyketides. PCR amplification of ketosynthase domains of type I modular polyketide synthases (PKS) from the microbial community of the marine sponge Discodermia dissoluta revealed great diversity and a novel group of sponge-specific PKS ketosynthase domains. Metagenomic libraries totaling more than four gigabases of bacterial genomes associated with this sponge were screened for type I modular PKS gene clusters. More than 90% of the clones in total sponge DNA libraries represented bacterial DNA inserts, and 0.7% harbored PKS genes. The majority of the PKS hybridizing clones carried small PKS clusters of one to three modules, although some clones encoded large multimodular PKSs (more than five modules). The most abundant large modular PKS appeared to be encoded by a bacterial symbiont that made up < 1% of the sponge community. Sequencing of this PKS revealed 14 modules that, if expressed and active, is predicted to produce a multimethyl-branched fatty acid reminiscent of mycobacterial lipid components. Metagenomic libraries made from fractions enriched for unicellular or filamentous bacteria differed significantly, with the latter containing numerous nonribosomal peptide synthetase (NRPS) and mixed NRPS-PKS gene clusters. The filamentous bacterial community of D. dissoluta consists mainly of Entotheonella spp., an unculturable sponge-specific taxon previously implicated in the biosynthesis of bioactive peptides.     

Developing tools for engineering hybrid polyketide synthetic pathways

Bacterial type I polyketide synthases (PKSs) are complex, multifunctional enzymes that synthesize structurally diverse and medicinally important natural products. Given their modular organization, the manipulation of type I PKSs holds tremendous promise for the generation of novel compounds that are not easily accessible by standard synthetic chemical approaches. In theory, hybrid polyketide synthetic pathways can be constructed through the rational recombination of catalytic domains or modules from a variety of PKS systems; however, the general success of this strategy has been elusive, largely due to a poor understanding of the interactions between catalytic domains, as well as PKS modules. Over the past several years, a fundamental knowledge of these issues, and others, has begun to emerge, offering refined strategies for the facile engineering of hybrid polyketide pathways.     

Computational approach for prediction of domain organization and substrate specificity of modular polyketide synthases

Modular polyketide synthases (PKSs) are large multi-enzymatic, multi-domain megasynthases, which are involved in the biosynthesis of a class of pharmaceutically important natural products, namely polyketides. These enzymes harbor a set of repetitive active sites termed modules and the domains present in each module dictate the chemical moiety that would add to a growing polyketide chain. This modular logic of biosynthesis has been exploited with reasonable success to produce several novel compounds by genetic manipulation. However, for harnessing their vast potential of combinatorial biosynthesis, it is essential to develop knowledge based in silico approaches for correlating the sequence and domain organization of PKSs to their polyketide products. In this work, we have carried out extensive sequence analysis of experimentally characterized PKS clusters to develop an automated computational protocol for unambiguous identification of various PKS domains in a polypeptide sequence. A structure based approach has been used to identify the putative active site residues of acyltransferase (AT) domains, which control the specificities for various starter and extender units during polyketide biosynthesis. On the basis of the analysis of the active site residues and molecular modelling of substrates in the active site of representative AT domains, we have identified a crucial residue that is likely to play a major role in discriminating between malonate and methylmalonate during selection of extender groups by this domain. Structural modelling has also explained the experimentally observed chiral preference of AT domain in substrate selection. This computational protocol has been used to predict the domain organization and substrate specificity for PKS clusters from various microbial genomes. The results of our analysis as well as the computational tools for prediction of domain organization and substrate specificity have been organized in the form of a searchable computerized database (PKSDB). PKSDB would serve as a valuable tool for identification of polyketide products biosynthesized by uncharacterized PKS clusters. This database can also provide guidelines for rational design of experiments to engineer novel polyketides.     

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.     

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.     

Unraveling polyketide synthesis in members of the genus Aspergillus

Aspergillus species have the ability to produce a wide range of secondary metabolites including polyketides that are generated by multi-domain polyketide synthases (PKSs). Recent biochemical studies using dissected single or multiple domains from PKSs have provided deep insight into how these PKSs control the structural outcome. Moreover, the recent genome sequencing of several species has greatly facilitated the understanding of the biosynthetic pathways for these secondary metabolites. In this review, we will highlight the current knowledge regarding polyketide biosynthesis in Aspergillus based on the domain architecture of non-reducing, highly reducing, and partially reducing PKSs, and PKS-non-ribosomal peptide synthetases.     

A new reducing polyketide synthase gene from the lichen-forming fungus Cladonia metacorallifera

Lichens produce unique polyketide secondary metabolites including depsides, depsidones, dibenzofurans and depsones. The biosynthesis of these compounds is governed by polyketide synthase (PKS), but the mechanism via which they are produced has remained unclear until now. We reported the 6-methylsalicylic acid synthase (6-MSAS) type of PKS gene, which is a member of the fungal-reducing PKSs. A cultured mycobiont of Cladonia metacorallifera was employed in the isolation and characterization of a polyketide synthase gene (CmPKS1). The complete sequence information for CmPKS1 was acquired via the screening of a Fosmid genomic library with a 456 bp fragment corresponding to part of the acyl transferase (AT) domain as a probe. CmPKS1 contains β-ketoacyl synthase (KS), AT, dehydratase (DH), ketoreductase (KR) and phosphopantetheine attachment site (PP) domains.: The domain organization of CmPKS1 (KS-AT-DH-KR-PP) is a typical 6-MSAS-type PKS, and the results of phylogenetic analysis showed that CmPKS1 grouped with other fungal-reducing PKSs. Quantitative real time PCR analyses showed that CmPKS1 was expressed preferentially in the early growth stage of the axenically cultured mycobiont. Furthermore CmPKS1 expression was found to be dependent on the carbon sources and concentrations in the medium.