The non-mevalonate pathway of isoprenoids: genes, enzymes and intermediates

Inhibition of nitrogen assimilation promotes carbon-based secondary metabolism in callus of Cyclocarya paliurus

The biosynthesis and accumulation of secondary metabolites are critical important to quality formation of medicinal plants, which are usually give way to primary processes and growth. Here, methionine sulfoximine (MSO) was used to inhibit the nitrogen assimilation in callus of Cyclocarya paliurus. The newly assimilated nitrogen characterized by ¹⁵N atom percentage excess, and the levels of amino acid and protein were reduced. The other primary processes such as carbohydrate metabolism and lipid metabolism were also repressed. In addition, the expression of the growth-related target of rapamycin (TOR) signaling was repressed, indicating nitrogen assimilation inhibition led to a systematic down-regulated primary metabolisms and resulted in a disruption of growth. In contrast, the biosynthesis of flavonoids and triterpenoids, antioxidase system, and the SnRK2-mediated abscisic acid (ABA) and jasmonic acid (JA) signaling were induced, which can improve plant stress resistance and defense. Nitrogen assimilation inhibition led to the carbon metabolic flux redirection from primary processes to secondary pathways, and facilitated the biosynthesis of flavonoids and triterpenoids in calluses of C. paliurus. Our results provide a comprehensive understanding of metabolic flux redirection between primary and secondary metabolic pathways and a potential means to improve the quality of medicinal plants.

Integrated genomics and proteomics analysis of Paenibacillus peoriae IBSD35 and insights into its antimicrobial characteristics

Antimicrobial resistance has been developing fast and incurring a loss of human life, and there is a need for new antimicrobial agents. Naturally occurring antimicrobial peptides offer the characteristics to counter AMR because the resistance development is low or no resistance. Antimicrobial peptides from Paenibacillus peoriae IBSD35 cell-free supernatant were salted out and purified using chromatography and characterized with liquid chromatography-tandem-mass spectrometry. The extract has shown a high and broad spectrum of antimicrobial activity. Combining the strain IBSD35 genome sequence with its proteomic data enabled the prediction of biosynthetic gene clusters by connecting the peptide from LC-MS/MS data to the gene that encode. Antimicrobial peptide databases offered a platform for the effective search, prediction, and design of AMPs and expanded the studies on their isolation, structure elucidation, biological evaluation, and pathway engineering. The genome-based taxonomy and comparisons have shown that P. peoriae IBSD35 is closely related to Paenibacillus peoriae FSL J3-0120. P. peoriae IBSD35 harbored endophytic trait genes and nonribosomal peptide synthases biosynthetic gene clusters. The comparative genomics revealed evolutionary insights and facilitated the discovery of novel SMs using proteomics from the extract of P. peoriae IBSD35. It will increase the potential to find novel bio-molecules to counter AMR.

The convergence of bacterial natural products from evolutionarily distinct pathways

As bacteria readily convert simple starting materials into a diverse array of complex molecules with useful bioactivities, these microorganisms and their biosynthetic machinery represent attractive alternatives to traditional chemical syntheses. While the well-documented divergent evolution of biosynthesis has allowed bacteria to explore wide swaths of natural product chemical space, the convergent evolution of these pathways remains a comparably rare phenomenon. The emergence of similar phenotypes within disparate genetic contexts provides a unique opportunity to probe the limitations of natural selection and the predictability and reproducibility of evolution under different constraints. Here, we report several recent examples of functional and structural convergence of bacterial natural products, as well as intra- and inter-domain convergence of bacterial biosynthetic machinery. While the genetic underpinnings of biosynthetic pathway evolution are of fundamental interest, the evolutionary constraints exemplified by phenotypic convergence also have immediate implications for efforts to engineer microorganisms for therapeutic small molecule production.

Functional Characterization of the 1-Deoxy-D-Xylulose 5-Phosphate Synthase Genes in Morus notabilis

Terpenoids are considered to be the largest group of secondary metabolites and natural products. Studies have revealed 1-deoxy-D-xylulose 5-phosphate synthase (DXS) is the first and rate-limiting enzyme in the plastidial methylerythritol phosphate pathway, which produces isopentenyl diphosphate and its isoform dimethylallyl diphosphate as terpenoid biosynthesis precursors. Mulberry (Morus L.) is an economically and ecologically important perennial tree with diverse secondary metabolites, including terpenoids that protect plants against bacteria and insects and may be useful for treating human diseases. However, there has been relatively little research regarding DXS genes in mulberry and other woody plant species. In this study, we cloned and functionally characterized three Morus notabilis DXS genes (MnDXS1, MnDXS2A, and MnDXS2B). Bioinformatics analyses indicated MnDXS1 belongs to clade 1, whereas MnDXS2A and MnDXS2B are in clade 2. The three encoded MnDXS proteins are localized to chloroplasts. Additionally, substantial differences in MnDXS expression patterns were observed in diverse tissues and in response to insect feeding and methyl jasmonate treatment. Moreover, overexpression of MnDXS1 in Arabidopsis thaliana increased the gibberellic acid content and resulted in early flowering, whereas overexpression of MnDXS2A enhanced root growth and increased the chlorophyll and carotenoid content. Our findings indicate that MnDXS functions vary among the clades, which may be useful for further elucidation of the functions of the DXS genes in mulberry.

Fungal Pigments: Potential Coloring Compounds for Wide Ranging Applications in Textile Dyeing

Synthetic pigments/non-renewable coloring sources used normally in the textile industry release toxic substances into the environment, causing perilous ecological challenges. To be safer from such challenges of synthetic colorants, academia and industries have explored the use of natural colorants such as microbial pigments. Such explorations have created a fervent interest among textile stakeholders to undertake the dyeing of textile fabrics, especially with fungal pigments. The biodegradable and sustainable production of natural colorants from fungal sources stand as being comparatively advantageous to synthetic dyes. The prospective scope of fungal pigments has emerged in the opening of many new avenues in textile colorants for wide ranging applications. Applying the biotechnological processes, fungal pigments like carotenoids, melanins, flavins, phenazines, quinones, monascins, violacein, indigo, etc. could be extracted on an industrial scale. This review appraises the studies and applications of various fungal pigments in dyeing textile fabrics and is furthermore shedding light on the importance of toxicity testing, genetic manipulations of fungal pigments, and their future perspectives under biotechnological approaches.

Bacterial Semiochemicals and Transkingdom Interactions with Insects and Plants

A peculiar feature of all living beings is their capability to communicate. With the discovery of the quorum sensing phenomenon in bioluminescent bacteria in the late 1960s, it became clear that intraspecies and interspecies communications and social behaviors also occur in simple microorganisms such as bacteria. However, at that time, it was difficult to imagine how such small organisms-invisible to the naked eye-could influence the behavior and wellbeing of the larger, more complex and visible organisms they colonize. Now that we know this information, the challenge is to identify the myriad of bacterial chemical signals and communication networks that regulate the life of what can be defined, in a whole, as a meta-organism. In this review, we described the transkingdom crosstalk between bacteria, insects, and plants from an ecological perspective, providing some paradigmatic examples. Second, we reviewed what is known about the genetic and biochemical bases of the bacterial chemical communication with other organisms and how explore the semiochemical potential of a bacterium can be explored. Finally, we illustrated how bacterial semiochemicals managing the transkingdom communication may be exploited from a biotechnological point of view.

An ImmunoPEGliposome for Targeted Antimalarial Combination Therapy at the Nanoscale

Combination therapies, where two drugs acting through different mechanisms are administered simultaneously, are one of the most efficient approaches currently used to treat malaria infections. However, the different pharmacokinetic profiles often exhibited by the combined drugs tend to decrease treatment efficacy as the compounds are usually eliminated from the circulation at different rates. To circumvent this obstacle, we have engineered an immunoliposomal nanovector encapsulating hydrophilic and lipophilic compounds in its lumen and lipid bilayer, respectively. The antimalarial domiphen bromide has been encapsulated in the liposome membrane with good efficiency, although its high IC₅₀ of ca. 1 µM for living parasites complicates its use as immunoliposomal therapy due to erythrocyte agglutination. The conjugation of antibodies against glycophorin A targeted the nanocarriers to Plasmodium-infected red blood cells and to gametocytes, the sole malaria parasite stage responsible for the transmission from the human to the mosquito vector. The antimalarials pyronaridine and atovaquone, which block the development of gametocytes, have been co-encapsulated in glycophorin A-targeted immunoliposomes. The co-immunoliposomized drugs have activities significantly higher than their free forms when tested in in vitro Plasmodium falciparum cultures: Pyronaridine and atovaquone concentrations that, when encapsulated in immunoliposomes, resulted in a 50% inhibition of parasite growth had no effect on the viability of the pathogen when used as free drugs.

Transcriptomic analysis in Anemone flaccida rhizomes reveals ancillary pathway for triterpene saponins biosynthesis and differential responsiveness to phytohormones

Anemone flaccida Fr. Schmidt is a perennial medicinal herb that contains pentacyclic triterpenoid saponins as the major bioactive constituents. In China, the rhizomes are used as treatments for a variety of ailments including arthritis. However, yields of the saponins are low, and little is known about the plant's genetic background or phytohormonal responsiveness. Using one-quarter of the 454 pyrosequencing information from the Roche GS FLX Titanium platform, we performed a transcriptomic analysis to identify 157 genes putatively encoding 26 enzymes involved in the synthesis of the bioactive compounds. It was revealed that there are two biosynthetic pathways of triterpene saponins in A. flaccida. One pathway depends on β-amyrin synthase and is similar to that found in other plants. The second, subsidiary ("backburner") pathway is catalyzed by camelliol C synthase and yields β-amyrin as minor byproduct. Both pathways used cytochrome P450-dependent monooxygenases (CYPs) and family 1 uridine diphosphate glycosyltransferases (UGTs) to modify the triterpenoid backbone. The expression of CYPs and UGTs were quite different in roots treated with the phytohormones methyl jasmonate, salicylic acid and indole-3-acetic acid. This study provides the first large-scale transcriptional dataset for the biosynthetic pathways of triterpene saponins and their phytohormonal responsiveness in the genus Anemone.

Structural characterization of 1-deoxy-D-xylulose 5-phosphate Reductoisomerase from Vibrio vulnificus

Vibrio vulnificus, a gram-negative bacterium, is the leading cause of seafood-borne illnesses and mortality in the United States. Previous studies have identified metabolites 2-C-methylerythritol 4-phosphate (MEP) as being essential for V. vulnificus growth and function. It was shown that 1-deoxy-D-xylulose-5-phosphate reductoisomerase (Dxr) is a critical enzyme in the viability of V. vulnificus, and many other bacteria, as it catalyzes the rearrangement of 1-deoxy-D-xylulose-5-phosphate (Dxp) to 2-C-methylerythritol 4-phosphate (MEP) within the MEP pathway, found in plants and bacteria. The MEP pathway produces the isoprenoids, isopentenyl diphosphate and dimethylallyl pyrophosphate. In this study, we produced and structurally characterized V. vulnificus Dxr. The enzyme forms a dimeric assembly and contains a metal ion in the active site. Protein produced in Escherichia coli co-purifies with Mg²⁺ ions, however the Mg²⁺ cations may be substituted with Mn²⁺, as both of these metals may be utilized by Dxrs. These findings will provide a basis for the design of Dxr inhibitors that may find application as antimicrobial compounds.

Structural modeling identifies Plasmodium vivax 4-diphosphocytidyl-2C-methyl-d-erythritol kinase (IspE) as a plausible new antimalarial drug target

Malaria parasites utilize Methylerythritol phosphate (MEP) pathway for synthesis of isoprenoid precursors which are essential for maturation and survival of parasites during erythrocytic and gametocytic stages. The absence of MEP pathway in the human host establishes MEP pathway enzymes as a repertoire of essential drug targets. The fourth enzyme, 4-diphosphocytidyl-2C-methyl-d-erythritol kinase (IspE) has been proved essential in pathogenic bacteria, however; it has not yet been studied in any Plasmodium species. This study was undertaken to investigate genetic polymorphism and concomitant structural implications of the Plasmodium vivax IspE (PvIspE) by employing sequencing, modeling and bioinformatics approach. We report that PvIspE gene displayed six non-synonymous mutations which were restricted to non-conserved regions within the gene from seven topographically distinct malaria-endemic regions of India. Phylogenetic studies reflected that PvIspE occupies unique status within Plasmodia genus and reflects that Plasmodium vivax IspE gene has a distant and non-conserved relation with human ortholog Mevalonate Kinase (MAVK). Structural modeling analysis revealed that all PvIspE Indian isolates have critically conserved canonical galacto-homoserine-mevalonate-phosphomevalonate kinase (GHMP) domain within the active site lying in a deep cleft sandwiched between ATP and CDPME-binding domains. The active core region was highly conserved among all clinical isolates, may be due to >60% β-pleated rigid architecture. The mapped structural analysis revealed the critically conserved active site of PvIspE, both sequence, and spacially among all Indian isolates; showing no significant changes in the active site. Our study strengthens the candidature of Plasmodium vivax IspE enzyme as a future target for novel antimalarials.

RNA-sequencing-based transcriptome and biochemical analyses of steroidal saponin pathway in a complete set of Allium fistulosum-A. cepa monosomic addition lines

The genus Allium is a rich source of steroidal saponins, and its medicinal properties have been attributed to these bioactive compounds. The saponin compounds with diverse structures play a pivotal role in Allium’s defense mechanism. Despite numerous studies on the occurrence and chemical structure of steroidal saponins, their biosynthetic pathway in Allium species is poorly understood. The monosomic addition lines (MALs) of the Japanese bunching onion (A. fistulosum, FF) with an extra chromosome from the shallot (A. cepa Aggregatum group, AA) are powerful genetic resources that enable us to understand many physiological traits of Allium. In the present study, we were able to isolate and identify Alliospiroside A saponin compound in A. fistulosum with extra chromosome 2A from shallot (FF2A) and its role in the defense mechanism against Fusarium pathogens. Furthermore, to gain molecular insight into the Allium saponin biosynthesis pathway, high-throughput RNA-Seq of the root, bulb, and leaf of AA, MALs, and FF was carried out using Illumina's HiSeq 2500 platform. An open access Allium Transcript Database (Allium TDB, http://alliumtdb.kazusa.or.jp) was generated based on RNA-Seq data. The resulting assembled transcripts were functionally annotated, revealing 50 unigenes involved in saponin biosynthesis. Differential gene expression (DGE) analyses of AA and MALs as compared with FF (as a control) revealed a strong up-regulation of the saponin downstream pathway, including cytochrome P450, glycosyltransferase, and beta-glucosidase in chromosome 2A. An understanding of the saponin compounds and biosynthesis-related genes would facilitate the development of plants with unique saponin content and, subsequently, improved disease resistance.

Carotenoid Biosynthesis in Fusarium

Many fungi of the genus Fusarium stand out for the complexity of their secondary metabolism. Individual species may differ in their metabolic capacities, but they usually share the ability to synthesize carotenoids, a family of hydrophobic terpenoid pigments widely distributed in nature. Early studies on carotenoid biosynthesis in Fusariumaquaeductuum have been recently extended in Fusarium fujikuroi and Fusarium oxysporum, well-known biotechnological and phytopathogenic models, respectively. The major Fusarium carotenoid is neurosporaxanthin, a carboxylic xanthophyll synthesized from geranylgeranyl pyrophosphate through the activity of four enzymes, encoded by the genes carRA, carB, carT and carD. These fungi produce also minor amounts of β-carotene, which may be cleaved by the CarX oxygenase to produce retinal, the rhodopsin's chromophore. The genes needed to produce retinal are organized in a gene cluster with a rhodopsin gene, while other carotenoid genes are not linked. In the investigated Fusarium species, the synthesis of carotenoids is induced by light through the transcriptional induction of the structural genes. In some species, deep-pigmented mutants with up-regulated expression of these genes are affected in the regulatory gene carS. The molecular mechanisms underlying the control by light and by the CarS protein are currently under investigation.

Cloning and Expression Analysis of MEP Pathway Enzyme-encoding Genes in Osmanthus fragrans

The 2-C-methyl-d-erythritol 4-phosphate (MEP) pathway is responsible for the biosynthesis of many crucial secondary metabolites, such as carotenoids, monoterpenes, plastoquinone, and tocopherols. In this study, we isolated and identified 10 MEP pathway genes in the important aromatic plant sweet osmanthus (Osmanthus fragrans). Multiple sequence alignments revealed that 10 MEP pathway genes shared high identities with other reported proteins. The genes showed distinctive expression profiles in various tissues, or at different flower stages and diel time points. The qRT-PCR results demonstrated that these genes were highly expressed in inflorescences, which suggested a tissue-specific transcript pattern. Our results also showed that OfDXS1, OfDXS2, and OfHDR1 had a clear diurnal oscillation pattern. The isolation and expression analysis provides a strong foundation for further research on the MEP pathway involved in gene function and molecular evolution, and improves our understanding of the molecular mechanism underlying this pathway in plants.

Isoprenoid-Based Biofuels: Homologous Expression and Heterologous Expression in Prokaryotes

Enthusiasm for mining advanced biofuels from microbial hosts has increased remarkably in recent years. Isoprenoids are one of the highly diverse groups of secondary metabolites and are foreseen as an alternative to petroleum-based fuels. Most of the prokaryotes synthesize their isoprenoid backbone via the deoxyxylulose-5-phosphate pathway from glyceraldehyde-3-phosphate and pyruvate, whereas eukaryotes synthesize isoprenoids via the mevalonate pathway from acetyl coenzyme A (acetyl-CoA). Microorganisms do not accumulate isoprenoids in large quantities naturally, which restricts their application for fuel purposes. Various metabolic engineering efforts have been utilized to overcome the limitations associated with their natural and nonnatural production. The introduction of heterologous pathways/genes and overexpression of endogenous/homologous genes have shown a remarkable increase in isoprenoid yield and substrate utilization in microbial hosts. Such modifications in the hosts' genomes have enabled researchers to develop commercially competent microbial strains for isoprenoid-based biofuel production utilizing a vast array of substrates. The present minireview briefly discusses the recent advancement in metabolic engineering efforts in prokaryotic hosts for the production of isoprenoid-based biofuels, with an emphasis on endogenous, homologous, and heterologous expression strategies.

RNA-seq analysis of Paris polyphylla var. yunnanensis roots identified candidate genes for saponin synthesis

Paris polyphylla Smith var. yunnanensis (Franch.) Hand.-Mazz. is a rhizomatous, herbaceous, perennial plant that has been used for more than a thousand years in traditional Chinese medicine. It is facing extinction due to overharvesting. Steroids are the major therapeutic components in Paris roots, the commercial value of which increases with age. To date, no genomic data on the species have been available. In this study, transcriptome analysis of an 8-year-old root and a 4-year-old root provided insight into the metabolic pathways that generate the steroids. Using Illumina sequencing technology, we generated a high-quality sequence and demonstrated de novo assembly and annotation of genes in the absence of prior genome information. Approximately 87,577 unique sequences, with an average length of 614 bases, were obtained from the root cells. Using bioinformatics methods, we annotated approximately 65.51% of the unique sequences by conducting a similarity search with known genes in the National Center for Biotechnology Information's non-redundant database. The unique transcripts were functionally classified using the Gene Ontology hierarchy and the Kyoto Encyclopedia of Genes and Genomes database. Of 3082 genes that were identified as significantly differentially expressed between roots of different ages, 1518 (49.25%) were upregulated and 1564 (50.75%) were downregulated in the older root. Metabolic pathway analysis predicted that 25 unigenes were responsible for the biosynthesis of the saponins steroids. These data represent a valuable resource for future genomic studies on this endangered species and will be valuable for efforts to genetically engineer P. polyphylla and facilitate saponin-rich plant development.

A Terpene Synthase Is Involved in the Synthesis of the Volatile Organic Compound Sodorifen of Serratia plymuthica 4Rx13

Bacteria release a plethora of volatile organic compounds, including compounds with extraordinary structures. Sodorifen (IUPAC name: 1,2,4,5,6,7,8-heptamethyl-3-methylenebicyclo[3.2.1]oct-6-ene) is a recently identified and unusual volatile hydrocarbon that is emitted by the rhizobacterium Serratia plymuthica 4R×13. Sodorifen comprises a bicyclic ring structure solely consisting of carbon and hydrogen atoms, where every carbon atom of the skeleton is substituted with either a methyl or a methylene group. This unusual feature of sodorifen made a prediction of its biosynthetic origin very difficult and so far its biosynthesis is unknown. To unravel the biosynthetic pathway we performed genome and transcriptome analyses to identify candidate genes. One knockout mutant (SOD_c20750) showed the desired negative sodorifen phenotype. Here it was shown for the first time that this gene is indispensable for the synthesis of sodorifen and strongly supports the hypothesis that sodorifen descends from the terpene metabolism. SOD_c20750 is the first bacterial terpene cyclase isolated from Serratia spp. and Enterobacteriales. Homology modeling revealed a 3D structure, which exhibits a functional role of amino acids for intermediate cation stabilization (W325) and putative proton acception (Y332). Moreover, the size and hydrophobicity of the active site strongly indicates that indeed the enzyme may catalyze the unusual compound sodorifen.

Likelihood-based gene annotations for gap filling and quality assessment in genome-scale metabolic models

Genome-scale metabolic models provide a powerful means to harness information from genomes to deepen biological insights. With exponentially increasing sequencing capacity, there is an enormous need for automated reconstruction techniques that can provide more accurate models in a short time frame. Current methods for automated metabolic network reconstruction rely on gene and reaction annotations to build draft metabolic networks and algorithms to fill gaps in these networks. However, automated reconstruction is hampered by database inconsistencies, incorrect annotations, and gap filling largely without considering genomic information. Here we develop an approach for applying genomic information to predict alternative functions for genes and estimate their likelihoods from sequence homology. We show that computed likelihood values were significantly higher for annotations found in manually curated metabolic networks than those that were not. We then apply these alternative functional predictions to estimate reaction likelihoods, which are used in a new gap filling approach called likelihood-based gap filling to predict more genomically consistent solutions. To validate the likelihood-based gap filling approach, we applied it to models where essential pathways were removed, finding that likelihood-based gap filling identified more biologically relevant solutions than parsimony-based gap filling approaches. We also demonstrate that models gap filled using likelihood-based gap filling provide greater coverage and genomic consistency with metabolic gene functions compared to parsimony-based approaches. Interestingly, despite these findings, we found that likelihoods did not significantly affect consistency of gap filled models with Biolog and knockout lethality data. This indicates that the phenotype data alone cannot necessarily be used to discriminate between alternative solutions for gap filling and therefore, that the use of other information is necessary to obtain a more accurate network. All described workflows are implemented as part of the DOE Systems Biology Knowledgebase (KBase) and are publicly available via API or command-line web interface.

Synthetic biology and metabolic engineering for marine carotenoids: new opportunities and future prospects

Carotenoids are a class of diverse pigments with important biological roles such as light capture and antioxidative activities. Many novel carotenoids have been isolated from marine organisms to date and have shown various utilizations as nutraceuticals and pharmaceuticals. In this review, we summarize the pathways and enzymes of carotenoid synthesis and discuss various modifications of marine carotenoids. The advances in metabolic engineering and synthetic biology for carotenoid production are also reviewed, in hopes that this review will promote the exploration of marine carotenoid for their utilizations.

β-carotene-producing bacteria residing in the intestine provide vitamin A to mouse tissues in vivo

Vitamin A deficiency (VAD) is an overwhelming public health problem that affects hundreds of millions of people worldwide. A definitive solution to VAD has yet to be identified. Because it is an essential nutrient, vitamin A or its carotenoid precursor β-carotene can only be obtained from food or supplements. In this study, we wanted to establish whether β-carotene produced in the mouse intestine by bacteria synthesizing the provitamin A carotenoid could be delivered to various tissues within the body. To achieve this, we took advantage of the Escherichia coli MG1655*, an intestine-adapted spontaneous mutant of E. coli MG1655, and the plasmid pAC-BETA, containing the genes coding for the 4 key enzymes of the β-carotene biosynthetic pathway (geranylgeranyl pyrophosphate synthase, phytoene synthase, phytoene desaturase, and lycopene cyclase) from Erwinia herbicola. We engineered the E. coli MG1655* to produce β-carotene during transformation with pAC-BETA (MG1655*-βC) and gavaged wild-type and knockout mice for the enzyme β-carotene 15,15'-oxygenase with this recombinant strain. Various regimens of bacteria administration were tested (single vs. multiple and low vs. high doses). β-Carotene concentration was measured by HPLC in mouse serum, liver, intestine, and feces. Enumeration of MG1655*-βC cells in the feces was performed to assess efficiency of intestinal colonization. We demonstrated in vivo that probiotic bacteria could be used to deliver vitamin A to the tissues of a mammalian host. These results have the potential to pave the road for future investigations aimed at identifying alternative, novel approaches to treat VAD.

Crystal structures of IspF from Plasmodium falciparum and Burkholderia cenocepacia: comparisons inform antimicrobial drug target assessment

Background

2C-methyl-D-erythritol-2,4-cyclodiphosphate synthase (IspF) catalyzes the conversion of 4-diphosphocytidyl-2C-methyl-D-erythritol-2-phosphate to 2C-methyl-D-erythritol-2,4-cyclodiphosphate and cytidine monophosphate in production of isoprenoid-precursors via the methylerythritol phosphate biosynthetic pathway. IspF is found in the protozoan Plasmodium falciparum, a parasite that causes cerebral malaria, as well as in many Gram-negative bacteria such as Burkholderia cenocepacia. IspF represents a potential target for development of broad-spectrum antimicrobial drugs since it is proven or inferred as essential in these pathogens and absent from mammals. Structural studies of IspF from these two important yet distinct pathogens, and comparisons with orthologues have been carried out to generate reagents, to support and inform a structure-based approach to early stage drug discovery.

Results

Efficient recombinant protein production and crystallization protocols were developed, and high-resolution crystal structures of IspF from P. falciparum (Emphasis/Emphasis>IspF) and B. cenocepacia (BcIspF) in complex with cytidine nucleotides determined. Comparisons with orthologues, indicate a high degree of order and conservation in parts of the active site where Zn²⁺ is bound and where recognition of the cytidine moiety of substrate occurs. However, conformational flexibility is noted in that area of the active site responsible for binding the methylerythritol component of substrate. Unexpectedly, one structure of BcIspF revealed two molecules of cytidine monophosphate in the active site, and another identified citrate coordinating to the catalytic Zn²⁺. In both cases interactions with ligands appear to help order a flexible loop at one side of the active site. Difficulties were encountered when attempting to derive complex structures with other ligands.

Conclusions

High-resolution crystal structures of IspF from two important human pathogens have been obtained and compared to orthologues. The studies reveal new data on ligand binding, with citrate coordinating to the active site Zn²⁺ and when present in high concentrations cytidine monophosphate displays two binding modes in the active site. Ligand binding appears to order a part of the active site involved in substrate recognition. The high degree of structural conservation in and around the IspF active site suggests that any structural model might be suitable to support a program of structure-based drug discovery.

Methylerythritol 4-phosphate (MEP) pathway metabolic regulation

The methylerythritol 4-phosphate pathway provides precursors for isoprenoids in bacteria, some eukaryotic parasites, and chloroplasts of plants. Metabolic regulatory mechanisms control flux through the pathway and the concentration of a central intermediate, methylerythritol cyclodiphosphate.

De Novo transcriptome sequencing reveals important molecular networks and metabolic pathways of the plant, Chlorophytum borivilianum

Chlorophytum borivilianum, an endangered medicinal plant species is highly recognized for its aphrodisiac properties provided by saponins present in the plant. The transcriptome information of this species is limited and only few hundred expressed sequence tags (ESTs) are available in the public databases. To gain molecular insight of this plant, high throughput transcriptome sequencing of leaf RNA was carried out using Illumina's HiSeq 2000 sequencing platform. A total of 22,161,444 single end reads were retrieved after quality filtering. Available (e.g., De-Bruijn/Eulerian graph) and in-house developed bioinformatics tools were used for assembly and annotation of transcriptome. A total of 101,141 assembled transcripts were obtained, with coverage size of 22.42 Mb and average length of 221 bp. Guanine-cytosine (GC) content was found to be 44%. Bioinformatics analysis, using non-redundant proteins, gene ontology (GO), enzyme commission (EC) and kyoto encyclopedia of genes and genomes (KEGG) databases, extracted all the known enzymes involved in saponin and flavonoid biosynthesis. Few genes of the alkaloid biosynthesis, along with anticancer and plant defense genes, were also discovered. Additionally, several cytochrome P450 (CYP450) and glycosyltransferase unique sequences were also found. We identified simple sequence repeat motifs in transcripts with an abundance of di-nucleotide simple sequence repeat (SSR; 43.1%) markers. Large scale expression profiling through Reads per Kilobase per Million mapped reads (RPKM) showed major genes involved in different metabolic pathways of the plant. Genes, expressed sequence tags (ESTs) and unique sequences from this study provide an important resource for the scientific community, interested in the molecular genetics and functional genomics of C. borivilianum.

Transcriptome analysis of Salicornia europaea under saline conditions revealed the adaptive primary metabolic pathways as early events to facilitate salt adaptation

BACKGROUND

Halophytes such as Salicornia europaea have evolved to exhibit unique mechanisms controlled by complex networks and regulated by numerous genes and interactions to adapt to habitats with high salinity. However, these mechanisms remain unknown.

METHODS

To investigate the mechanism by which halophytes tolerate salt based on changes in the whole transcriptome, we performed transcriptome sequencing and functional annotation by database search. Using the unigene database, we conducted digital gene expression analysis of S. europaea at various time points after these materials were treated with NaCl. We also quantified ion uptakes. Gene functional enrichment analysis was performed to determine the important pathways involved in this process.

RESULTS

A total of 57,151 unigenes with lengths of >300 bp were assembled, in which 57.5% of these unigenes were functionally annotated. Differentially expressed genes indicated that cell wall metabolism and lignin biosynthetic pathways were significantly enriched in S. europaea to promote the development of the xylem under saline conditions. This result is consistent with the increase in sodium uptake as ions pass through the xylem. Given that PSII efficiency remained unaltered, salt treatment activated the expression of electron transfer-related genes encoded by the chloroplast chromosome. Chlorophyll biosynthesis was also inhibited, indicating the energy-efficient state of the electron transfer system of S. europaea.

CONCLUSIONS

The key function of adjusting important primary metabolic pathways in salt adaption was identified by analyzing the changes in the transcriptome of S. europaea. These pathways could involve unique salt tolerance mechanisms in halophytes. This study also provided information as the basis of future investigations on salt response genes in S. europaea. Ample gene resources were also provided to improve the genes responsible for the salt tolerance ability of crops.

Formal Synthesis of 4-diphosphocytidyl-2-C-methyl D-erythritol From D-(+)-Arabitol

2-C-methyl-D-erythritol-4-phosphate (MEP) is a key chemical intermediate of the non-mevalonate pathway for isoprenoid biosynthesis employed by many pathogenic microbes. MEP is also the precursor for the synthesis of 4-diphosphocytidyl-2-C-methyl D-erythritol (CDP-ME), another key intermediate of the non-mevalonate pathway. As this pathway is non-existent in higher animals, including humans, it represents great opportunities for novel antimicrobial development. To facilitate the in-depth studies of this pathway, we reported here a formal synthesis of CDP-ME through a new synthesis of 2-C-Methyl-D-erythritol-4-phosphoric acid from D-(+)-arabitol.

Tanshinones: sources, pharmacokinetics and anti-cancer activities

Tanshinones are a class of abietane diterpene compound isolated from Salvia miltiorrhiza (Danshen or Tanshen in Chinese), a well-known herb in Traditional Chinese Medicine (TCM). Since they were first identified in the 1930s, more than 40 lipophilic tanshinones and structurally related compounds have been isolated from Danshen. In recent decades, numerous studies have been conducted to investigate the isolation, identification, synthesis and pharmacology of tanshinones. In addition to the well-studied cardiovascular activities, tanshinones have been investigated more recently for their anti-cancer activities in vitro and in vivo. In this review, we update the herbal and alternative sources of tanshinones, and the pharmacokinetics of selected tanshinones. We discuss anti-cancer properties and identify critical issues for future research. Whereas previous studies have suggested anti-cancer potential of tanshinones affecting multiple cellular processes and molecular targets in cell culture models, data from in vivo potency assessment experiments in preclinical models vary greatly due to lack of uniformity of solvent vehicles and routes of administration. Chemical modifications and novel formulations had been made to address the poor oral bioavailability of tanshinones. So far, human clinical trials have been far from ideal in their design and execution for the purpose of supporting an anti-cancer indication of tanshinones.

Differentially expressed transcripts from leaf and root tissue of Chlorophytum borivilianum: a plant with high medicinal value

Chlorophytum borivilianum is one of the important medicinal plants used for treating different health problems such as diabetes, arthritis, physical weakness, etc. Saponins present in C. borivilianum are the primary source of its significant medicinal properties and are synthesized by mevalonate and non-mevalonate pathways in plants. However, the biosynthesis of these compounds at molecular level is not studied in C. borivilianum. Cloning and sequencing of genes involved in metabolic processes are prerequisite to study the gene expression, their regulation and genetic engineering experiments. Expressed sequence tags (ESTs) provide a quick insight into various genes and their tissue specific expression. Suppression subtractive hybridization (SSH) libraries were constructed using mRNA from leaf and root tissues of C. borivilianum. High quality non-redundant 506 and 303 ESTs were generated from leaf and root specific libraries respectively. These sequences were analyzed using bioinformatics tools and grouped into different categories based on their similarity and cellular functions such as photosynthesis, metabolism, transcription factors, cell signaling, defense, stress response etc. ESTs also showed similarity with genes involved in saponins biosynthesis such as squalene synthase, squalene epoxidase, cytochrome p450, glycosyltransferase, etc. Semi-quantitative analysis of some of the ESTs involved in saponins biosynthesis confirmed their differential regulation in leaves and roots. These ESTs will provide an efficient resource to accelerate gene discovery in C. borivilianum and will help in determining promising targets for genetic engineering of saponins pathway.

Pre-steady-state kinetic analysis of 1-deoxy-D-xylulose-5-phosphate reductoisomerase from Mycobacterium tuberculosis reveals partially rate-limiting product release by parallel pathways

As part of the non-mevalonate pathway for the biosynthesis of the isoprenoid precursor isopentenyl pyrophosphate, 1-deoxy-D-xylulose-5-phosphate (DXP) reductoisomerase (DXR) catalyzes the conversion of DXP into 2-C-methyl-D-erythritol 4-phosphate (MEP) by consecutive isomerization and NADPH-dependent reduction reactions. Because this pathway is essential to many infectious organisms but is absent in humans, DXR is a target for drug discovery. In an attempt to characterize its kinetic mechanism and identify rate-limiting steps, we present the first complete transient kinetic investigation of DXR. Stopped-flow fluorescence measurements with Mycobacterium tuberculosis DXR (MtDXR) revealed that NADPH and MEP bind to the free enzyme and that the two bind together to generate a nonproductive ternary complex. Unlike the Escherichia coli orthologue, MtDXR exhibited a burst in the oxidation of NADPH during pre-steady-state reactions, indicating a partially rate-limiting step follows chemistry. By monitoring NADPH fluorescence during these experiments, the transient generation of MtDXR·NADPH·MEP was observed. Global kinetic analysis supports a model involving random substrate binding and ordered release of NADP(+) followed by MEP. The partially rate-limiting release of MEP occurs via two pathways--directly from the binary complex and indirectly via the MtDXR·NADPH·MEP complex--the partitioning being dependent on NADPH concentration. Previous mechanistic studies, including kinetic isotope effects and product inhibition, are discussed in light of this kinetic mechanism.

Identification and characterization of xiamycin A and oxiamycin gene cluster reveals an oxidative cyclization strategy tailoring indolosesquiterpene biosynthesis

Xiamycin A (XMA) and oxiamycin (OXM) are bacterial indolosesquiterpenes featuring rare pentacyclic ring systems and are isolated from a marine-derived Streptomyces sp. SCSIO 02999. The putative biosynthetic gene cluster for XMA/OXM was identified by a partial genome sequencing approach. Eighteen genes were proposed to be involved in XMA/OXM biosynthesis, including five genes for terpene synthesis via a non-mevalonate pathway, eight genes encoding oxidoreductases, and five genes for regulation and resistance. Targeted disruptions of 13 genes within the xia gene cluster were carried out to probe their encoded functions in XMA/OXM biosynthesis. The disruption of xiaK, encoding an aromatic ring hydroxylase, led to a mutant producing indosespene and a minor amount of XMA. Feeding of indosespene to XMA/OXM nonproducing mutants revealed indosespene as a common precursor for XMA/OXM biosynthesis. Most notably, the flavin dependent oxygenase XiaI was biochemically characterized in vitro to convert indosespene to XMA, revealing an unusual oxidative cyclization strategy tailoring indolosesquiterpene biosynthesis.

IspE inhibitors identified by a combination of in silico and in vitro high-throughput screening

CDP-ME kinase (IspE) contributes to the non-mevalonate or deoxy-xylulose phosphate (DOXP) pathway for isoprenoid precursor biosynthesis found in many species of bacteria and apicomplexan parasites. IspE has been shown to be essential by genetic methods and since it is absent from humans it constitutes a promising target for antimicrobial drug development. Using in silico screening directed against the substrate binding site and in vitro high-throughput screening directed against both, the substrate and co-factor binding sites, non-substrate-like IspE inhibitors have been discovered and structure-activity relationships were derived. The best inhibitors in each series have high ligand efficiencies and favourable physico-chemical properties rendering them promising starting points for drug discovery. Putative binding modes of the ligands were suggested which are consistent with established structure-activity relationships. The applied screening methods were complementary in discovering hit compounds, and a comparison of both approaches highlights their strengths and weaknesses. It is noteworthy that compounds identified by virtual screening methods provided the controls for the biochemical screens.

Biochemistry of the non-mevalonate isoprenoid pathway

The non-mevalonate pathway of isoprenoid (terpenoid) biosynthesis is essential in many eubacteria including the major human pathogen, Mycobacterium tuberculosis, in apicomplexan protozoa including the Plasmodium spp. causing malaria, and in the plastids of plants. The metabolic route is absent in humans and is therefore qualified as a promising target for new anti-infective drugs and herbicides. Biochemical and structural knowledge about all enzymes involved in the pathway established the basis for discovery and development of inhibitors by high-throughput screening of compound libraries and/or structure-based rational design.

Integrative genome-scale metabolic analysis of Vibrio vulnificus for drug targeting and discovery

Chromosome 1 of Vibrio vulnificus tends to contain larger portion of essential or housekeeping genes on the basis of the genomic analysis and gene knockout experiments performed in this study, while its chromosome 2 seems to have originated and evolved from a plasmid.

The genome-scale metabolic network model of V. vulnificus was reconstructed based on databases and literature, and was used to identify 193 essential metabolites.

Five essential metabolites finally selected after the filtering process are 2-amino-4-hydroxy-6-hydroxymethyl-7,8-dihydropteridine (AHHMP), D-glutamate (DGLU), 2,3-dihydrodipicolinate (DHDP), 1-deoxy-D-xylulose 5-phosphate (DX5P), and 4-aminobenzoate (PABA), which were predicted to be essential in V. vulnificus, absent in human, and are consumed by multiple reactions.

Chemical analogs of the five essential metabolites were screened and a hit compound showing the minimal inhibitory concentration (MIC) of 2 μg/ml and the minimal bactericidal concentration (MBC) of 4 μg/ml against V. vulnificus was identified.

Discovering new antimicrobial targets and consequently new antimicrobials is important as drug resistance of pathogenic microorganisms is becoming an increasingly serious problem in human healthcare management (Fischbach and Walsh, 2009). There clearly exists a gap between genomic studies and drug discovery as the accumulation of knowledge on pathogens at genome level has not successfully transformed into the development of effective drugs (Mills, 2006; Payne et al, 2007). In this study, we dissected the genome of a microbial pathogen in detail, and subsequently developed a systems biological strategy of employing genome-scale metabolic modeling and simulation together with metabolite essentiality analysis for effective drug targeting and discovery. This strategy was used for identifying new drug targets in an opportunistic pathogen Vibrio vulnificus CMCP6 as a model.

V. vulnificus is a Gram-negative halophilic bacterium that is found in estuarine waters, brackish ponds, or coastal areas, and its Biotype 1 is an opportunistic human pathogen that can attack immune-compromised patients, and causes primary septicemia, necrotized wound infections, and gastroenteritis. We previously found that many metabolic genes were specifically induced in vivo, suggesting that specific metabolic pathways are essential for in vivo survival and virulence of this pathogen (Kim et al, 2003; Lee et al, 2007). These results motivated us to carry out systems biological analysis of the genome and the metabolic network for new drug target discovery.

V. vulnificus CMCP6 has two chromosomes. We first re-sequenced genomic regions assembled in low quality and low depth, and subsequently re-annotated the whole genome of V. vulnificus. Horizontal gene transfer was suspected to be responsible for the diversification of each chromosome of V. vulnificus, and the presence of metabolic genes was more biased to chromosome 1 than chromosome 2. Further studies on V. vulnificus genome revealed that chromosome 2 is more prone to diversification for better adaptation to the environment than its chromosome 1, while chromosome 1 tends to expand their genetic repertoire while maintaining the core genes at a constant level.

Next, a genome-scale metabolic network VvuMBEL943 was reconstructed based on literature, databases and experiments for systematic studies on the metabolism of this pathogen and prediction of drug targets. The VvuMBEL943 model is composed of 943 reactions and 765 metabolites, and covers 673 genes. The model was validated by comparing its simulated cell growth phenotype obtained by constraints-based flux analysis with the V. vulnificus-specific experimental data previously reported in the literature. In this study, constraints-based flux analysis is an optimization-based simulation method that calculates intracellular fluxes under the specific genetic and environmental condition (Kim et al, 2008). As a result, 17 growth phenotypes were correctly predicted out of 18 cases, which demonstrate the validity of VvuMBEL943.

The main objective of constructing VvuMBEL943 in this study is to predict potential drug targets by system-wide analysis of the metabolic network for the effective treatment of V. vulnificus. To achieve this goal, a set of drug target candidates was predicted by taking a metabolite-centric approach. Metabolite essentiality analysis is a concept recently introduced for the study of cellular robustness to complement conventional reaction or gene-centric approach (Kim et al, 2007b). Metabolite essentiality analysis observes changes in flux distribution by removing each metabolite from the in silico metabolic network. Hence, metabolite essentiality predicts essential metabolites whose absence causes cell death. By selecting essential metabolites, it is possible to directly screen only their structural analogs, which substantially reduces the number of chemical compounds to screen from the chemical compound library. As a result of implementing this approach, 193 metabolites were initially identified to be essential to the cell. These essential metabolites were then further filtered based on the predetermined criteria, mainly organism specificity and multiple connectivity associated with each metabolite, in order to reduce the number of initial target candidates towards identifying the most effective ones.

Five essential metabolites finally selected are 2-amino-4-hydroxy-6-hydroxymethyl-7,8-dihydropteridine (AHHMP), D-glutamate (DGLU), 2,3-dihydrodipicolinate (DHDP), 1-deoxy-D-xylulose 5-phosphate (DX5P), and 4-aminobenzoate (PABA). Enzymes that consume these essential metabolites were experimentally verified to be essential, which indeed demonstrates the essentiality of these five metabolites. On the basis of the structural information of these five essential metabolites, whole-cell screening assay was performed using their analogs for possible antibacterial discovery. We screened 352 chemical analogs of the essential metabolites selected from the chemical compound library, and found a hit compound 24837, which shows the minimal inhibitory concentration (MIC) of 2 μg/ml and minimal bactericidal concentration (MBC) of 4 μg/ml, showing good antibacterial activity without further structural modification. Although this study demonstrates a proof-of-concept, the approaches and their rationale taken here should serve as a general strategy for discovering novel antibiotics and drugs based on systems-level analysis of metabolic networks.

Although the genomes of many microbial pathogens have been studied to help identify effective drug targets and novel drugs, such efforts have not yet reached full fruition. In this study, we report a systems biological approach that efficiently utilizes genomic information for drug targeting and discovery, and apply this approach to the opportunistic pathogen Vibrio vulnificus CMCP6. First, we partially re-sequenced and fully re-annotated the V. vulnificus CMCP6 genome, and accordingly reconstructed its genome-scale metabolic network, VvuMBEL943. The validated network model was employed to systematically predict drug targets using the concept of metabolite essentiality, along with additional filtering criteria. Target genes encoding enzymes that interact with the five essential metabolites finally selected were experimentally validated. These five essential metabolites are critical to the survival of the cell, and hence were used to guide the cost-effective selection of chemical analogs, which were then screened for antimicrobial activity in a whole-cell assay. This approach is expected to help fill the existing gap between genomics and drug discovery.

From Zn to Mn: the study of novel manganese-binding groups in the search for new drugs against tuberculosis

In most eubacteria, apicomplexans, and most plants, including the causal agents for diseases such as malaria, leprosy, and tuberculosis, the methylerythritol phosphate pathway is the route for the biosynthesis of the C₅ precursors to the essential isoprenoid class of compounds. Owing to their absence in humans, the enzymes of the methylerythritol phosphate pathway have become attractive targets for drug discovery. This work investigates a new class of inhibitors against the second enzyme of the pathway, 1-deoxy-d-xylulose 5-phosphate reductoisomerase. Inhibition of this enzyme may involve the chelation of a crucial active site Mn ion, and the metal-chelating moieties studied here have previously been shown to be successful in application to the zinc-dependent metalloproteinases. Quantum mechanics and docking calculations presented in this work suggest the transferability of these metal-chelating compounds to Mn-containing 1-deoxy-d-xylulose 5-phosphate reductoisomerase enzyme, as a promising starting point to the development of potent inhibitors.

IspG-catalyzed positional isotopic exchange in methylerythritol cyclodiphosphate of the deoxyxylulose phosphate pathway: mechanistic implications

H(2)(18)O under the bridge: Recently, the deoxyxylulose phosphate (DXP) pathway was discovered to be a second pathway supplying isoprenoid biosynthetic precursors. One of steps is an IspG-catalyzed reductive deoxygenation of methylerythritol cyclodiphosphate (MEcPP) to 4-hydroxyl-3-methyl-2-(E)-1-diphosphate (HMBPP). Using [2-(13) C,(18) O]-MEcPP, we detected the positional isotopic exchange for the bridging oxygen in MEcPP.

Methylerythritol cyclodiphosphate (MEcPP) in deoxyxylulose phosphate pathway: synthesis from an epoxide and mechanisms

In this communication, we reported another unique IspG-catalyzed transformation, the production of its substrate, MEcPP, from (2R,3R)-4-hydroxy-3-methyl-2,3-epoxybutanyl diphosphate (Epoxy-HMBPP) when reductants are excluded from the reaction mixture.

A single-molecule force spectroscopy nanosensor for the identification of new antibiotics and antimalarials

An important goal of nanotechnology is the application of individual molecule handling techniques to the discovery of potential new therapeutic agents. Of particular interest is the search for new inhibitors of metabolic routes exclusive of human pathogens, such as the 2-C-methyl-D-erythritol-4-phosphate (MEP) pathway essential for the viability of most human pathogenic bacteria and of the malaria parasite. Using atomic force microscopy single-molecule force spectroscopy (SMFS), we have probed at the single-molecule level the interaction of 1-deoxy-D-xylulose 5-phosphate synthase (DXS), which catalyzes the first step of the MEP pathway, with its two substrates, pyruvate and glyceraldehyde-3-phosphate. The data obtained in this pioneering SMFS analysis of a bisubstrate enzymatic reaction illustrate the substrate sequentiality in DXS activity and allow for the calculation of catalytic parameters with single-molecule resolution. The DXS inhibitor fluoropyruvate has been detected in our SMFS competition experiments at a concentration of 10 μM, improving by 2 orders of magnitude the sensitivity of conventional enzyme activity assays. The binding of DXS to pyruvate is a 2-step process with dissociation constants of k(off) = 6.1 × 10(-4) ± 7.5 × 10(-3) and 1.3 × 10(-2) ± 1.0 × 10(-2) s(-1), and reaction lengths of x(β) = 3.98 ± 0.33 and 0.52 ± 0.23 Å. These results constitute the first quantitative report on the use of nanotechnology for the biodiscovery of new antimalarial enzyme inhibitors and open the field for the identification of compounds represented only by a few dozens of molecules in the sensor chamber.

IspG converts an epoxide substrate analogue to (E)-4-hydroxy-3-methylbut-2-enyl diphosphate: implications for IspG catalysis in isoprenoid biosynthesis

IspG is an intriguing enzyme in bacteria, parasite, and plant isoprenoid biosynthesis, and its catalytic mechanism remains elusive. We report here the synthesis of (2R,3R)-4-hydroxy-3-methyl-2,3-epoxybutanyl diphosphate (Epoxy-HMBPP), a proposed intermediate in one of the frequently cited mechanistic models. We have also demonstrated that this epoxide analogue is a catalytically competent IspG substrate. This study represents the first mechanistic study of this important enzyme.

Simulation of a Petri net-based model of the terpenoid biosynthesis pathway

Background

The development and simulation of dynamic models of terpenoid biosynthesis has yielded a systems perspective that provides new insights into how the structure of this biochemical pathway affects compound synthesis. These insights may eventually help identify reactions that could be experimentally manipulated to amplify terpenoid production. In this study, a dynamic model of the terpenoid biosynthesis pathway was constructed based on the Hybrid Functional Petri Net (HFPN) technique. This technique is a fusion of three other extended Petri net techniques, namely Hybrid Petri Net (HPN), Dynamic Petri Net (HDN) and Functional Petri Net (FPN).

Results

The biological data needed to construct the terpenoid metabolic model were gathered from the literature and from biological databases. These data were used as building blocks to create an HFPNe model and to generate parameters that govern the global behaviour of the model. The dynamic model was simulated and validated against known experimental data obtained from extensive literature searches. The model successfully simulated metabolite concentration changes over time (pt) and the observations correlated with known data. Interactions between the intermediates that affect the production of terpenes could be observed through the introduction of inhibitors that established feedback loops within and crosstalk between the pathways.

Conclusions

Although this metabolic model is only preliminary, it will provide a platform for analysing various high-throughput data, and it should lead to a more holistic understanding of terpenoid biosynthesis.