Debasisa Mohanty

NRPS-PKS: a knowledge-based resource for analysis of NRPS/PKS megasynthases

NRPS-PKS is web-based software for analysing large multi-enzymatic, multi-domain megasynthases that are involved in the biosynthesis of pharmaceutically important natural products such as cyclosporin, rifamycin and erythromycin. NRPS-PKS has been developed based on a comprehensive analysis of the sequence and structural features of several experimentally characterized biosynthetic gene clusters. The results of these analyses have been organized as four integrated searchable databases for elucidating domain organization and substrate specificity of nonribosomal peptide synthetases and three types of polyketide synthases. These databases work as the backend of NRPS-PKS and provide the knowledge base for predicting domain organization and substrate specificity of uncharacterized NRPS/PKS clusters. Benchmarking on a large set of biosynthetic gene clusters has demonstrated that, apart from correct identification of NRPS and PKS domains, NRPS-PKS can also predict specificities of adenylation and acyltransferase domains with reasonably high accuracy. These features of NRPS-PKS make it a valuable resource for identification of natural products biosynthesized by NRPS/PKS gene clusters found in newly sequenced genomes. The training and test sets of gene clusters included in NRPS-PKS correlate information on 307 open reading frames, 2223 functional protein domains, 68 starter/extender precursors and their specific recognition motifs, and also the chemical structure of 101 natural products from four different families. NRPS-PKS is a unique resource which provides a user-friendly interface for correlating chemical structures of natural products with the domains and modules in the corresponding nonribosomal peptide synthetases or polyketide synthases. It also provides guidelines for domain/module swapping as well as site-directed mutagenesis experiments to engineer biosynthesis of novel natural products. NRPS-PKS can be accessed at http://www.nii.res.in/nrps-pks.html.

Enzymic activation and transfer of fatty acids as acyl-adenylates in mycobacteria

The metabolic repertoire in nature is augmented by generating hybrid metabolites from a limited set of gene products. In mycobacteria, several unique complex lipids are produced by the combined action of fatty acid synthases and polyketide synthases (PKSs), although it is not clear how the covalently sequestered biosynthetic intermediates are transferred from one enzymatic complex to another. Here we show that some of the 36 annotated fadD genes, located adjacent to the PKS genes in the Mycobacterium tuberculosis genome, constitute a new class of long-chain fatty acyl-AMP ligases (FAALs). These proteins activate long-chain fatty acids as acyl-adenylates, which are then transferred to the multifunctional PKSs for further chain extension. This mode of activation and transfer of fatty acids is contrary to the previously described universal mechanism involving the formation of acyl-coenzyme A thioesters. Similar mechanisms may operate in the biosynthesis of other lipid-containing metabolites and could have implications in engineering novel hybrid products.

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.

SBSPKS: structure based sequence analysis of polyketide synthases

Polyketide synthases (PKSs) catalyze biosynthesis of a diverse family of pharmaceutically important secondary metabolites. Bioinformatics analysis of sequence and structural features of PKS proteins plays a crucial role in discovery of new natural products by genome mining, as well as in design of novel secondary metabolites by biosynthetic engineering. The availability of the crystal structures of various PKS catalytic and docking domains, and mammalian fatty acid synthase module prompted us to develop SBSPKS software which consists of three major components. Model_3D_PKS can be used for modeling, visualization and analysis of 3D structure of individual PKS catalytic domains, dimeric structures for complete PKS modules and prediction of substrate specificity. Dock_Dom_Anal identifies the key interacting residue pairs in inter-subunit interfaces based on alignment of inter-polypeptide linker sequences to the docking domain structure. In case of modular PKS with multiple open reading frames (ORFs), it can predict the cognate order of substrate channeling based on combinatorial evaluation of all possible interface contacts. NRPS–PKS provides user friendly tools for identifying various catalytic domains in the sequence of a Type I PKS protein and comparing them with experimentally characterized PKS/NRPS clusters cataloged in the backend databases of SBSPKS. SBSPKS is available at http://www.nii.ac.in/sbspks.html.

SEARCHPKS: a program for detection and analysis of polyketide synthase domains

SEARCHPKS is a software for detection and analysis of polyketide synthase (PKS) domains in a polypeptide sequence. Modular polyketide synthases are unusually large multi-enzymatic multi-domain megasynthases, which are involved in the biosynthesis of pharmaceutically important natural products using an assembly-line mechanism. This program facilitates easy identification of various PKS domains and modules from a given polypeptide sequence. In addition, it also predicts the specificity of the potential acyltransferase domains for various starter and extender precursor units. SEARCHPKS is a user-friendly tool for correlating polyketide chemical structures with the organization of domains and modules in the corresponding modular polyketide synthases. This program also allows the user to extensively analyze and assess the sequence homology of various polyketide synthase domains, thus providing guidelines for carrying out domain and module swapping experiments. SEARCHPKS can also aid in identification of polyketide products made by PKS clusters found in newly sequenced genomes. The computational approach used in SEARCHPKS is based on a comprehensive analysis of various characterized clusters of modular polyketide synthases compiled in PKSDB, a database of modular polyketide synthases. SEARCHPKS can be accessed at http://www.nii.res.in/searchpks.html.

A new family of type III polyketide synthases in Mycobacterium tuberculosis.

The Mycobacterium tuberculosis genome has revealed a remarkable array of polyketide synthases (PKSs); however, no polyketide product has been isolated thus far. Most of the PKS genes have been implicated in the biosynthesis of complex lipids. We report here the characterization of two novel type III PKSs from M. tuberculosis that are involved in the biosynthesis of long-chain α-pyrones. Measurement of steady-state kinetic parameters demonstrated that the catalytic efficiency of PKS18 protein was severalfold higher for long-chain acyl-coenzyme A substrates as compared with the small-chain precursors. The specificity of PKS18 and PKS11 proteins toward long-chain aliphatic acyl-coenzyme A (C12 to C20) substrates is unprecedented in the chalcone synthase (CHS) family of condensing enzymes. Based on comparative modeling studies, we propose that these proteins might have evolved by fusing the catalytic machinery of CHS and β-ketoacyl synthases, the two evolutionarily related members with conserved thiolase fold. The mechanistic and structural importance of several active site residues, as predicted by our structural model, was investigated by performing site-directed mutagenesis. The functional identification of diverse catalytic activity in mycobacterial type III PKSs provide a fascinating example of metabolite divergence in CHS-like proteins.

Mechanistic and functional insights into fatty acid activation in Mycobacterium tuberculosis

The recent discovery of fatty acyl-AMP ligases (FAALs) in Mycobacterium tuberculosis (Mtb) provided a new perspective to fatty acid activation dogma. These proteins convert fatty acids to corresponding adenylates, which is an intermediate of acyl-CoA-synthesizing fatty acyl-CoA ligases (FACLs). Presently, it is not evident how obligate pathogens like Mtb have evolved such new themes of functional versatility and whether the activation of fatty acids to acyl-adenylates could indeed be a general mechanism. Here, based on elucidation of the first structure of a FAAL protein and by generating loss- as well as gain-of-function mutants that interconvert FAAL and FACL activities, we demonstrate that an insertion motif dictates formation of acyl-adenylate. Since FAALs in Mtb are crucial nodes in biosynthetic network of virulent lipids, inhibitors directed against these proteins provide a unique multi-pronged approach of simultaneously disrupting several pathways.

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.

In silico analysis of methyltransferase domains involved in biosynthesis of secondary metabolites

BackgroundSecondary metabolites biosynthesized by polyketide synthase (PKS) and nonribosomal peptide synthetase (NRPS) family of enzymes constitute several classes of therapeutically important natural products like erythromycin, rapamycin, cyclosporine etc. In view of their relevance for natural product based drug discovery, identification of novel secondary metabolite natural products by genome mining has been an area of active research. A number of different tailoring enzymes catalyze a variety of chemical modifications to the polyketide or nonribosomal peptide backbone of these secondary metabolites to enhance their structural diversity. Therefore, development of powerful bioinformatics methods for identification of these tailoring enzymes and assignment of their substrate specificity is crucial for deciphering novel secondary metabolites by genome mining.ResultsIn this work, we have carried out a comprehensive bioinformatics analysis of methyltransferase (MT) domains present in multi functional type I PKS and NRPS proteins encoded by PKS/NRPS gene clusters having known secondary metabolite products. Based on the results of this analysis, we have developed a novel knowledge based computational approach for detecting MT domains present in PKS and NRPS megasynthases, delineating their correct boundaries and classifying them as N-MT, C-MT and O-MT using profile HMMs. Analysis of proteins in nr database of NCBI using these class specific profiles has revealed several interesting examples, namely, C-MT domains in NRPS modules, N-MT domains with significant homology to C-MT proteins, and presence of NRPS/PKS MTs in association with other catalytic domains. Our analysis of the chemical structures of the secondary metabolites and their site of methylation suggested that a possible evolutionary basis for the presence of a novel class of N-MT domains with significant homology to C-MT proteins could be the close resemblance of the chemical structures of the acceptor substrates, as in the case of pyochelin and yersiniabactin. These two classes of MTs recognize similar acceptor substrates, but transfer methyl groups to N and C positions on these substrates.ConclusionWe have developed a novel knowledge based computational approach for identifying MT domains present in type I PKS and NRPS multifunctional enzymes and predicting their site of methylation. Analysis of nr database using this approach has revealed presence of several novel MT domains. Our analysis has also given interesting insight into the evolutionary basis of the novel substrate specificities of these MT proteins.

SEARCHGTr: a program for analysis of glycosyltransferases involved in glycosylation of secondary metabolites

SEARCHGTr is a web-based software for the analysis of glycosyltransferases (GTrs) involved in the biosynthesis of a variety of pharmaceutically important compounds like adriamycin, erythromycin, vancomycin etc. This software has been developed based on a comprehensive analysis of sequence/structural features of 102 GTrs of known specificity from 52 natural product biosynthetic gene clusters. SEARCHGTr is a powerful tool that correlates sequences of GTrs to the chemical structures of their corresponding substrates. This software indicates the donor/acceptor specificity and also identifies putative substrate binding residues. In addition, it provides interfaces to other public databases like GENBANK, SWISS-PROT, CAZY, PDB, PDBSum and PUBMED for extracting various information on GTrs homologous to the query sequence. SEARCHGTr would provide new dimension to our previously developed bioinformatics tool NRPS-PKS. Together, these tools facilitate comprehensive computational analysis of proteins involved in biosynthesis of aglycone core and its downstream glycosylations. Apart from presenting opportunities for rational design of novel natural products, these tools would assist in the identification of biosynthetic products of secondary metabolite gene clusters found in newly sequenced genomes. SEARCHGTr can be accessed at .