Simulation of a Petri net-based Model of the Terpenoid Biosynthesis Pathway
BMC Bioinformatics volume 11, Article number: 83 (2010)
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).
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.
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.
Biological networks are complex and may consist of hundreds of reactions that directly and indirectly affect each other. In order to understand the relationship between these reactions in detail, it is necessary to view the network as a whole . Computational models of metabolic networks have been developed in order to provide an overview of the biosynthetic processes involved in multiple pathways [2, 3]. Using such models, changes in the concentration of each metabolite corresponding to those that occur during normal and perturbed conditions can be simulated and used to determine how such changes may contribute to the entire biosynthetic process. The use of computational models of this type gives researchers an in-depth view of the problems that need to be solved and points to new strategies and alternatives.
Various approaches can be taken to building and simulating a biological pathway model . The most popular way of describing the behaviour of such a model is the Ordinary Differential Equations (ODEs) approach [5–7]. ODEs are typically derived from the Michaelis-Menten equation. They are then embedded into the model through script writing using programs such as SimBiology in Matlab  and Gepasi  that were developed for building and editing such models. However, choosing the right approach to building a particular biological model relies on the type of biological pathway it is to represent. Typically, there are three categories of biological pathways: gene regulatory networks [10, 11], metabolic pathways [6, 12–14] and signalling pathways [4, 7, 15]. The model presented here, which describes the biosynthesis of terpenoids via two independent pathways, one involving mevalonate (MEV) and one involving methylerythritol phosphate (MEP), falls into the metabolic pathways category.
We chose to model the terpenoid biosynthetic network using Hybrid Functional Petri net with extension (HFPNe). This method allows easy modelling of complex biological information using a graphical approach and also allows various types of entities to be modelled together. The technique is derived from the traditional Petri net theory that was first described by Carl Adam Petri in his 1962 PhD dissertation . HFPNe is an extension of the HFPN architecture with extended features resulting from fusion of three extended Petri net techniques, namely Hybrid Petri net (HPN), Dynamic Petri net (HDN) and Functional Petri net (FPN) [17–20]. This technique is suitable for our modelling purposes for several reasons: (a) it is capable of utilising primitive types (Boolean and string); (b) it is able to treat more than one value in an entity; and (c) it is capable of treating object type for complicated bioprocesses. Finally, HFPNe introduces the new notion of a generic entity, as well as processing and enhancing simulation processes for this notion .
HFPNe-based biological models are built based on three basic entities called places, transitions and arcs; these are shown in Figure 1[21, 22]. In addition to possessing the ability to model continuous and discrete entities (places and transitions) concurrently as described in HDNs, HFPNe introduces generic entities that can hold values in the form of both real number and Boolean strings. This feature enables on/off switches and parameter modulators to be created, which can then function as regulatory components in the model. These HFPNe entities form a network and each is assigned with specific parameters that describe its behaviour, making this technique distinct from other Petri net techniques [10, 21–25]. The assignment of specific parameters to the HFPNe entities allows more options to be customised in order to control the behaviour of the model. The parameters of the entities and the topology of their network are generated from biological facts that can be obtained through experimental procedures and by literature and database searches [10, 23]. HFPNe and other types of Petri net techniques allow an intuitive approach to modelling; their ability to provide a visual modelling framework for development and simulation [13, 14, 21, 25–29] motivated us to use this technique in our modelling. HFPNe architecture has been successfully applied to the modelling of processes such as vulval development  and the operation of microRNA regulatory networks  in Caenorhabditis elegans.
The biosynthesis of the universal precursors isopentenyl diphosphate (IPP) and dimethylallyl diphosphate (DMAPP) occurs by two independent pathways, the mevalonate (MEV) pathway and the methylerythritol phosphate (MEP) pathway. These two pathways are distinct in terms of organism, cellular localisation and products as well as in the metabolites and precursors involved. The classical pathway, the MEV pathway, was thought to be the only pathway involved in terpenoid biosynthesis until the MEP pathway was discovered in the early 90s [32, 33]. The MEV pathway usually operates in the cytoplasm of eukaryotic organisms such as Gallus gallus and Rattus norvegicus, whereas the MEP pathway takes place in the plastids of eubacteria such as Escherichia coli and Streptomyces cerevisiae, as well as in the chloroplasts of plants [34–36]. The MEP pathway is not found in animals or fungi, making it attractive for drug design studies because the enzymes in this pathway are using a specific substrate or generate a specific product and their function cannot be compensated for by another enzyme [37, 38]. However, both pathways are operational in higher-level plants such as Arabidopsis thaliana and Helianthus annuu s, as well as in selected eubacteria such as Streptomyces sp. [34, 39–42].
The several classes of terpenoids (e.g., monoterpenes, sesquiterpenes, triterpenes and diterpenes) are organised according to the number of isoprene units in their backbone. Isoprene units are five-carbon compounds that act as the building block for terpenoids. Isoprenes (C5) consist of one isoprene unit, monoterpenes (C10) consist of two isoprene units, sesquiterpenes (C15) consist of three isoprene units, diterpenes (C20) consist of four isoprene units and the addition of more units will make up bigger terpene classes. The type of terpene produced also depends on the route of biosynthesis. The MEV pathway is more likely to produce sesqui- and triterpenes, whereas the MEP pathway produces mono- and diterpenes [34, 37, 43–45].
The topological network of the pathways in the model was constructed based on reaction stoichiometry and enzymatic mechanisms of known enzymatic reactions in the terpene biosynthetic pathway. Underlying biological information, including production and consumption rates, was used to generate the parameters of each HFPNe element that govern the behaviour of the model during simulations. Simulation and validation processes were then performed iteratively to produce an optimum system.
In this study, the possibility of discovering new alternatives for amplifying the production of terpenoids was the major goal. Therefore, a complete understanding of how the compounds are synthesised is vital.
Results and Discussion
This model consists of two subnetworks, the MEV and MEP subnetworks, which represent the respective terpene biosynthetic pathways. The study focuses on the synthesis of two classes of terpenes, sesquiterpenoids (the end product of the MEV pathway) and monoterpenoids (the end product of the MEP pathway). These two pathways are placed in different compartments in the model to reflect their different subcellular locations (Figure 2). All of the metabolites involved in the pathways are represented using HFPNe elements and are interconnected to each other based on their reaction stoichiometry and enzymatic mechanisms. The network consists of 101 continuous transitions representing various reactions and production/degradation processes; 66 places representing 61 metabolites; three on/off switches; and two parameter modulators. All metabolites were represented using continuous places, whereas generic places were used to create the regulatory components (on/off switches and parameter modulators). All metabolites were represented using continuous places whereas generic places were used to create the regulatory components (on/off switches and parameter modulators) (Figure 3, see also Additional file 1).
The MEV subnetwork consists of ten main enzymatic reactions and two feedback loops, whereas the MEP subnetwork consists of nine main enzymatic reactions and a crosstalk mechanism that also involves the MEV subnetwork. The 19 main enzymatic reactions participating in both subnetworks are shown in Tables 1 and 2 and Figure 4. The behaviour of each reaction in the pathway is described using HFPNe element parameters such as speed of transition formulae, initial place values and threshold arc values [18, 21, 22, 25]. To simulate the behaviour of the model in normal and/or perturbed conditions, the entities representing on/off switches are adjusted. The initial values of these entities can be manipulated to create different conditions during the simulation, leading to a better understanding of the pathway and pointing to reactions to be further manipulated. All simulations carried out in this study are based on known information retrieved from literature and database searches.
Based on the introduction of inhibitors and increments of metabolite concentration in both pathways, together with the regulation of the three switch entities, we here summarise five simulation conditions:
normal production of sesquiterpenes through the MEV subnetwork, with the MEP subnetwork and all inhibitory processes deactivated;
normal production of monoterpenes through the MEP subnetwork, with the MEV subnetwork and all inhibitory processes deactivated;
overexpression of ATP in the fifth reaction of the MEV subnetwork, with the MEP subnetwork and other inhibitory processes deactivated;
overproduction of FPP in the ninth reaction of the MEV subnetwork, with the MEP subnetwork and other inhibitory processes deactivated;
establishment of a crosstalk mechanism with the MEP subnetwork activated but with the MEV subnetwork as well as other inhibitory processes deactivated.
Simulation and Validation
The simulation and validation processes, as well as the developmental stages of the model performed in this study, were carried out using Cell Illustrator https://cionline.hgc.jp. The model was validated against known data and information obtained from extensive literature and database searches. Simulation results are returned as concentrations (unit) versus Petri net time (pt) graphs. Here, we discuss three of the five simulated conditions carried out in the study.
a. Overproduction of ATP
In the fifth reaction of the MEV pathway, phosphomevalonate kinase (PMK; EC 188.8.131.52) catalyses the production of mevalonate diphosphate (MEVPP) and ADP from mevalonate-5-phosphate (MEVP) and ATP (Table 1) [46–48]. This reaction acts in a bi bi (also known as two-substrate/two-product) sequential enzymatic mechanism with respect to the two substrates; MEVP first binds to the free enzyme, followed by ATP; the enzyme then releases the main product and finally ADP . Studies have shown that high concentration of ATP halts the production of MEVPP and thus blocks the entire pathway from continuing to synthesise other intermediates. ATP acts as a competitive inhibitor, binding to both sites of the free enzyme and forming a dead-end complex. This obstructs the binding of MEVP to the enzyme and impedes the downstream reactions in the pathway. However, it was also shown that by increasing the concentration of MEVP, the inhibitory effects of ATP could be reversed. Higher concentrations of MEVP inhibit the synthesis of ATP, lowering the concentration of ATP to a manageable level. This allows the synthesis of MEVPP, as well as downstream reactions, to resume (Figure 5) .
To simulate this mutually antagonist inhibitory feedback loop, the 'MEV_switch' component is set to 'true' and the other two switches are set to 'false.' The parameter that controls the production of ATP is adjusted to a higher value of 30 unit/pt from 10 unit/pt, resulting in overproduction of ATP by 3-fold and activates the inhibitory arc. The threshold value for the inhibitory arc that suppresses the production of ATP is set to 80 units to allow the inhibitory effect to take place.
b. Overproduction of FPP
The second feedback loop identified in the MEV pathway is a positive feedback mechanism. It is initiated when FPP levels become high enough to inhibit the catalytic activity of mevalonate kinase (MK; EC184.108.40.206). MK catalyses the fourth reaction in the MEV pathway, phosphorylation of mevalonate to produce mevalonate-5-phosphate (MEVP) (Table 1). The inhibitory activity of FPP halts the phosphorylation of mevalonate, thus slowing the downstream processes, including the production of FPP itself [49, 50]. This will return FPP concentration to its normal level, releasing its inhibitory effect on MK (Figure 6).
To simulate this behaviour, the 'MEV_switch' component is set to 'true' and the two other switches are set to 'false.' No other inhibitory processes except the inhibitory arc connecting FPP to the catalytic activity of MK are activated. The threshold value of the arc is set to 10 units from a higher value of 100 units to allow the inhibitory effect to take place. The simulation result shows an oscillatory behaviour of the concentrations of sesquiterpene, FPP and MEVPP over time (Figure 6).
A crosstalk mechanism between the MEP and MEV subnetworks is established with the introduction of the inhibitor fosmidomycin. This inhibitor suppresses the catalytic activity of deoxyxylulose-5-phosphate reductoisomerase (DXR; EC220.127.116.117). Suppression of DXR activity reduces the rates of the downstream reactions in the MEP pathway, affecting the production of monoterpenes (Figure 7). However, the inhibition also triggers a flux of precursors (IPP and DMAPP) from the MEP pathway to the MEV pathway, resulting in increased sesquiterpene concentration despite the inactivation of the MEV pathway (Figure 8) . This event is described in the model by assigning function values to the transition elements of the components that are involved in this mechanism. These functions are then regulated by a switch component called 'crossSwitch.'
To simulate the crosstalk mechanism, the 'MEV_switch' component is set to 'false.' The 'MEP_switch' component is set to 'true' while the 'crossSwitch' component is given a script value. The script value enables the activation of the 'fosmidomycin' component after 150 pt, thereby establishing the crosstalk mechanism. The simulation result successfully described oscillatory changes in the concentrations of the metabolites involved over time.
The simulation and validation processes performed using the model are consistent with known biological information and data. The intuitive approach introduced by the HFPNe technique enables intricate modelling tasks to be viewed and solved in a graphical perspective. The model serves as a tool for better understanding the reactions involved in both pathways and how they affect each other. Perturbations performed on the model provide insight into the effects of manipulation of a single reaction on the whole network and thus should facilitate both industrial and biomedical terpenoid bioengineering.
Information and data acquisition
The biological information and data needed to construct this model included information on the substrates, products and enzymes involved in the MEV and MEP pathways, as well as detailed enzymatic reaction information such as reaction stoichiometry and enzyme mechanisms. Consumption and production rates were also taken into consideration in developing the dynamic model. This information was used to design the model layout and parameter assignments for each of the HFPNe elements. The data were obtained through the use of extensive literature and biological databases including KEGG, BRENDA, ENZYME, IUBMB, MetaCyc and PATHWAY DATABASE [23, 52–56].
Model template design and parameter assignment
The HFPN elements (place, transition and arc) represent the metabolites and processes that comprise the MEV and MEP biosynthetic pathways. Continuous places represent the substrates, products, inhibitors and enzymes in both subnetworks, whereas continuous transitions represent biological processes such as catalysis, phosphorylation and degradation. Generic places represent the on/off switches and parameter modulators. The normal arcs then connect places (enzyme, substrate, product) to transitions while inhibitory arcs connect inhibitors to their prospective inhibited processes. All these components were arranged according to the reaction stoichiometry and enzymatic mechanism in their order of occurrence in the respective pathways (Tables 1 and 2).
Each of the HFPNe elements possesses different attributes that characterise its function in the model. Parameter values are assigned to each of these attributes to control the behaviour of the model during simulation. These values are generated based on the biological information obtained and are translated into virtual concentration values (initial values of places), ordinary differential equations (speed formulas of transitions), numeric values (threshold of normal and inhibitory arcs) and scripts (on/off switches) [10, 18, 22–25].
Iterative model development
The test-driven models were simulated and validated against known data obtained through literature searches. Simulations were carried out under five different conditions as explained in the Results and Discussion section. Simulation results were displayed as concentration (unit) versus time (pt) graphs. Petri net time (pt) indicates virtual Petri net time units that do not correspond to real time; concentration also is given in general concentration units (unit) that do not specifically correspond to actual concentration units such as mM and μM. The changes in metabolite concentrations (unit) over time predicted by each simulation were validated against known biological data to identify gaps and inconsistencies. Uncorrelated biological behaviours were re-examined and previous developmental steps (topological design and parameter assignment) were repeated. The simulation process was then re-executed and revalidated against known biological information. The entire process was carried out iteratively in order to rule out inconsistencies and obtain an optimum system.
(Hybrid Dynamic Petri net)
(Hybrid Functional Petri net with extension)
(Hybrid Petri net)
(Functional Petri net)
(Petri net time)
(ordinary differential equation)
Kitano H: Systems Biology: A Brief Overview. Science 2002, 295: 1662–1664. 10.1126/science.1069492
Sweetlove LJ, Last RL, Fernie AR: Predictive Metabolic Engineering: A Goal for Systems Biology. Plant Physiology 2003, 132: 420–425. 10.1104/pp.103.022004
Sweetlove LJ, Fernie AR: Regulation of Metabolic Networks: Understanding Metabolic Complexity in the Systems Biology Era. New Phytologist 2005, 168: 9–24. 10.1111/j.1469-8137.2005.01513.x
Gilbert D, FuB H, Gu X, Orton R, Robinson S, Vyshemirsky V, Kurth MJ, Downes CS, Dubitzky W: Computational Methodologies for Modelling, Analysis and Simulation of Signalling Networks. Briefings in Bioinformatics 2006, 7(4):339–353. 10.1093/bib/bbl043
Goryanin I, Hodgman TC, Selkov E: Mathematical Simulation and Analysis of Celular Metabolism and Regulation. Bioinformatics 1999, 15(9):749–758. 10.1093/bioinformatics/15.9.749
Yang K, Ma W, Liang H, Ouyang Q, Tang C, Lai L: Dynamic Simulations on the Arachidonic Acid Metabolic Network. PLoS Computational Biology 2007, 3(3):523–530. 10.1371/journal.pcbi.0030055
Orton RJ, Sturm OE, Vyshemirsky V, Calder M, Gilbert DR, Kolch W: Computational modelling of the receptor-tyrosine-kinase-activated MAPK pathway. Biochemical Journal 2005, 392: 249–261. 10.1042/BJ20050908
Schmidt H, Drews G, Vera J, Wolkenhauer O: SBML export interface for the systems biology toolbox for MATLAB. Bioinformatics Applications Note 2007, 23(10):1297–1298.
Mendes P: GEPASI a software for modeling the dynamics, steady states and control of biochemical and other systems. Comput Appl Biosci 1993, 9: 563–571.
Doi A, Nagasaki M, Matsuno H, Miyano S: Simulation-based validation of the p53 transcriptional activity with Hybrid Functional Petri Net. In Silico Biology 2006, 6(1):1–13.
Grunwald S, Speer A, Ackermann J, Koch I: Petri net modelling of gene regulation of the Duchenne muscular dystrophy. Biosystems 2008, 92(2):189–205. 10.1016/j.biosystems.2008.02.005
Oliveira JS, Bailey CG, Jones-Oliveira JB, Dixon DA, Gull DW, Chandler ML: A computational model for the identification of biochemical pathways in the Krebs cycle. Journal of Computational Biology 2003, 10(1):57–82. 10.1089/106652703763255679
Koch I, Heiner M: Petri nets in analysis of biological network. In Biological Network Analysis Edited by: Junker BH, Schreiber F. 2008, 139–180. full_text
Koch I, Junker BH, Heiner M: Application of Petri net Theory for Modelling and Validation of the Sucrose Breakdown Pathway in the Potato Tuber. Bioinformatics 2005, 21(7):1219–1226. 10.1093/bioinformatics/bti145
Sackmann A, Heiner M, Koch I: Application of Petri net based analysis techniques to signal transduction pathways. BMC Bioinformatics 2006, 7: 482–498. 10.1186/1471-2105-7-482
Petri CA, Bonn : Petri net. In Schriften des IIM Nr. 2 1962. Volume 1. 2nd edition. New York: Griffiss Air Force Base, Technical Report RADC-TR-65–377; 1966. English Translation English Translation
Nagasaki M: A platform for biopathway modeling/simulation and recreating biopathway databases towards simulation. In PhD Thesis. University of Tokyo; 2004.
Nagasaki M, Doi A, Matsuno H, Miyano S: A Versatile Petri net Based Architecture for Modeling and Simulation of Complex Biological Processes. Genome Informatics 2004, 15(1):180–197.
Nagasaki M, Doi A, Matsuno H, Miyano S: Petri net based description and modeling of biological pathways. In Proceedings of the first international conference on algebraic biology ---computer algebra in biology---: 2005. Tokyo: Universal Academy Press; 2005:19–31.
Nagasaki M, Saito A, Doi A, Matsuno H, Miyano S: Foundations of Systems Biology. London: Springer; 2009.
Matsuno H, Doi A, Miyano S: Hybrid Petri net Representation of Gene Regulatory Network. Pacific Symposium on Biocomputing: 2000 2000, 338–349.
Matsuno H, Tanaka Y, Aoshima H, Doi A, Matsui M, Miyano S: Biopathways representation and simulation on hybrid functional Petri net. In Silico Biology 2003., 3(32):
Doi A, Fujita S, Matsuno H, Nagasaki M, Miyano S: Constructing biological pathway models with hybrid functional Petri nets. In Silico Biology 2004, 4(3):271–291.
Matsuno H, Li C, Miyano S: Petri net based descriptions for systematic understanding of biological pathways. IECICE Trans Fundamentals 2006, E89(11):3166–3174. 10.1093/ietfec/e89-a.11.3166
Miyano S: How to Model and Simulate Biological Pathways with Petri nets - A new challenge for systems Biology-. 25th International Conference on Application and Theory of Petri nets Bologna, Italy 2004.
Rodriguez EM, del Rosario RCH, Rudy A, Vollmar A, Mendoza E: Cephalostatin 1-induced Apoptosis in Leukemic Cells: From Petri net Model to Kinetic Model. International Conference on Molecular Systems Biology: 2008; Manila, Filipina 2008, 114–117.
Steggles LJ, Banks R, Shaw O, Wipat A: Qualitatively modelling and analysing genetic regulatory networks: a Petri net approach. Bioinformatics 2007, 23(3):336–343. 10.1093/bioinformatics/btl596
Peleg M, Rubin D, Altman RB: Using Petri net tools to study properties and dynamics of biological systems. Journal of the American Medical Informatics Association 2005, 12(2):181–199. 10.1197/jamia.M1637
Peleg M, Yeh I, Altman RB: Modelling biological processes using workflow and Petri net models. Bioinformatics 2002, 18(6):825–837. 10.1093/bioinformatics/18.6.825
Li C, Nagasaki M, Ueno K, Miyano S: Simulation-based model checking approach to cell fate specification during Caenorhabditis elegans vulval development by hybrid functional Petri netPetri net with extension. BMC Systems Biology 2009., 3(32):
Saito A, Nagasaki M, Doi A, Ueno K, Miyano S: Cell fate simulation model of gustatory neurons with microRNAs double-negative feedback loop by Hybrid Functional Petri net with extension. Genome Informatics 2006, 17(1):100–111.
Hahn FM, Eubanks LM, Testa CA, Blagg BSJ, Baker JA, Poulter CD: 1-Deoxy-D-xylulose 5-Phosphate synthase, the gene product of open reading frame (ORF) 2816 and ORF 2895 in Rhodobacter capsulatus . Journal of Bacteriology 2001, 183(1):1–11. 10.1128/JB.183.1.1-11.2001
Lange BM, Wildung MR, McCaskill D, Croteau R: A family of transketolases that directs isoprenoid biosynthesis via a mevalonate-independent pathway. Proceedings of the National Academy of Sciences of the United States of America 1998., 95(2100–2104):
Rohdich F, Hecht S, Gartner K, Adam P, Krieger C, Amslinger S, Arigoni D, Bacher A, Eisenrich W: Studies on the nonmevalonate terpene biosynthetic pathway: Metabolic role of IspH (LytB) protein. Proceedings of the National Academy of Sciences of the United States of America 2002, 99(3):1158–1163. 10.1073/pnas.032658999
Voynova NE, Rios SE, Miziorko HM: Staphylococcus aureus mevalonate kinase: Isolation and characterization of an enzyme of the isoprenoid biosynthetic pathway. Journal of Bacteriology 2004, 186(1):61–67. 10.1128/JB.186.1.61-67.2004
Hoeffler J, Tritsch D, Grosdemange-Biliard C, Rohmer M: Isoprenoid biosynthesis via the methylerythritol phosphate pathway-Mechanistic investigations of the 1-deoxy-D-xylulose 5-phosphate reductoisomerase. European Journal of Biochemistry 2002, 269: 4446–4457. 10.1046/j.1432-1033.2002.03150.x
Eisenrich W, Rohdich F, Bacher A: Deoxyxylulose phosphate pathway to terpenoids. TRENDS in Plant Science 2001, 6(2):78–84. 10.1016/S1360-1385(00)01812-4
Hunter WN: The non-mevalonate pathway of isoprenoid precursor biosynthesis. Journal of Biological Chemistry 2007, 282(30):21573–21577. 10.1074/jbc.R700005200
Nagegowda DA, Bach TJ, Chye M: Brassica juncea 3-hydroxy-3-methylglutaryl (HMG)-CoA synthase 1: expression and characterization of recombinant wild-type and mutant enzymes. Biochemical Journal 2004, 383: 517–527. 10.1042/BJ20040721
Schulte AE, Heijden R, Verpoorte R: Purification and characterization of phosphomevalonate kinase from Catharanthus roseus . Phytochemistry 1999, 52(6):975–983. 10.1016/S0031-9422(99)00382-9
Kuzuyama T, Takagi M, Takahashi S, Seto H: Cloning and characterization of 1-Deoxy-D-Xylulose 5-Phosphate Synthase from Streptomyces sp. Strain CL190, which uses both the mevalonate and nonmevalonate pathways for isopentenyl diphosphate biosynthesis. Journal of Bacteriology 2000, 182(4):891–897. 10.1128/JB.182.4.891-897.2000
Kuzuyama T, Takahashi S, Takagi M, Seto H: Characterization of 1-Deoxy-D-xylulose 5-phosphate reductoisomerase, an enzyme involved in Isopentenyl diphosphate Biosynthesis and identification of its catalytic amino acid residues. The Journal of Biological Chemistry 2000, 275(26):19928–19932. 10.1074/jbc.M001820200
Dubey VS, Bhalla R, Luthra R: An overview of the non-mevalonate pathway for terpenoid biosynthesis in plants. Journal of Biosciences 2003, 28(5):637–646. 10.1007/BF02703339
Rohdich F, Kis K, Bacher A, Eisenrich W: The non-mevalonate pathway of isoprenoids: genes, enzymes and intermediates. Current Opinion in Chemical Biology 2001, 5: 535–540. 10.1016/S1367-5931(00)00240-4
McCaskill D, Croteau R: Isopentenyl Diphosphate is the terminal product of the deoxyxylulose-5-phosphate pathway for terpenoid biosynthesis in plants. Tetrahedron Letters 1999, 40: 653–656. 10.1016/S0040-4039(98)02496-4
Eyzaguirre J, Valdebenito D, Cardemil E: Pig Liver Phosphomevalonate kinase: Kinetic Mechanism. Archives of Biochemistry and Biophysics 2006, 454: 188–195. 10.1016/j.abb.2006.08.002
Doun SS, Burgner JW II, Briggs SD, Rodwell VW: Enterococcus faecalis phosphomevalonate kinase. Protein Science 2005, 14: 1134–1139. 10.1110/ps.041210405
Hogenboom S, Tuyp JJM, Espeel M, Koster J, Wanders RJA, Waterham HR: Phosphomevalonate kinase is a cytosolic protein in humans. Journal of Lipid Research 2004, 45: 697–705. 10.1194/jlr.M300373-JLR200
Schulte AE, Llamas Duran EM, Heijden R, Verpoorte R: Mevalonate kinase activity in Catharanthus roseus plants and suspension cultured cells. Plant Science 2000, 150: 59–69. 10.1016/S0168-9452(99)00164-8
Gray JC, Kekwick RGO: Mevalonate Kinase in Green Leaves and Etiolated Cotyledons of the French Bean Phaseolus vulgaris . Biochemical Journal 1973, 133: 335–347.
Schwender J, Muller C, Zeidler J, Lichtenthanler HK: Cloning and heterologous expression of a cDNA encoding 1-deoxy-D-xylulose-5-phosphate reductoisomerase of Arabidopsis thaliana . FEBS letters 1999, 455: 140–144. 10.1016/S0014-5793(99)00849-2
Kanehisa M, Araki M, Goto S, Hattori M, Hirakawa M, itoh M, Katayama T, Kawashima S, Okuda S, Tokimatsu T, Yamanishi Y: KEGG for linking genomes to life and the environment. Nucleic Acids Research 2008, (36 Database):D480-D484.
Kanehisa M, Goto S: KEGG: Kyoto Encyclopedia of Genes and Genomes. Nucleic Acids Research 2000, 28(1):27–30. 10.1093/nar/28.1.27
Kanehisa M, Goto S, Hattori M, Aoki-Kinoshita KF, Itoh M, Kawashima S, Katayama T, Araki M, Hirakawa M: From genomics to chemical genomics: new developements in KEGG. Nucleic Acids Research 2006, (34 Database):D354-D357. 10.1093/nar/gkj102
Chang A, Scheer M, Grote A, Schomburg I, Schomburg D: BRENDA, AMENDA and FRENDA the enzyme information system: new content and tools in 2009. Nucleic Acids Research 2009, (37 Database):D588-D592. 10.1093/nar/gkn820
Caspi R, Foerster H, Fulcher CA, Kaipa P, Krummenacker M, Latendresse M, Paley S, Rhee SY, Shearer AG, Tissier C, et al.: The MetaCyc Database of Metabolic Pathways and Enzymes and the BioCyc Collection of Pathway/Genome Databases. Nucleic Acids Research 2008, (36 Database):D623-D631.
AHH is a National Science Fellowship scholar under the Ministry of Science, Technology and Innovation (MOSTI), Malaysia. This study was supported by the Universiti Kebangsaan Malaysia (UKM) and Malaysia Genome Institute (MGI) initiative research project (grant UKM-MGI-NBD0008-2007 and UKM-MGI-NBD0005-2007). This work was carried out at the Bioinformatics Lab, Institute of Systems Biology (INBIOSIS) UKM. The facility and financial aid are duly acknowledged.
AHH performed the research and drafted the manuscript. ZAMH formulated the study, gave valuable insights, participated in its design and helped to draft and refine the manuscript. All authors read and approved the final manuscript.
Aliah Hazmah Hawari and Zeti-Azura Mohamed-Hussein contributed equally to this work.
Electronic supplementary material
Authors’ original submitted files for images
Below are the links to the authors’ original submitted files for images.
About this article
Cite this article
Hawari, A.H., Mohamed-Hussein, ZA. Simulation of a Petri net-based Model of the Terpenoid Biosynthesis Pathway. BMC Bioinformatics 11, 83 (2010). https://doi.org/10.1186/1471-2105-11-83
- Isoprene Unit
- Terpenoid Biosynthesis
- Methylerythritol Phosphate