Skip to main content
  • Research article
  • Open access
  • Published:

Simulation of a Petri net-based Model of the Terpenoid Biosynthesis Pathway



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 [1]. 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 [4]. The most popular way of describing the behaviour of such a model is the Ordinary Differential Equations (ODEs) approach [57]. 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 [8] and Gepasi [9] 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, 1214] 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 [16]. 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) [1720]. 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 [17].

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, 2125]. 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, 2529] motivated us to use this technique in our modelling. HFPNe architecture has been successfully applied to the modelling of processes such as vulval development [30] and the operation of microRNA regulatory networks [31] in Caenorhabditis elegans.

Figure 1
figure 1

HFPNe entities. Places, transitions and arcs are basic entities in the HFPNe architecture. Places and transitions can be subdivided into continuous, discrete and generic.

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 [3436]. 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, 3942].

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, 4345].

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

Model overview

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).

Figure 2
figure 2

Overall sketch of the terpenoid biosynthetic pathways. The entire network is subdivided into MEV and MEP subnetworks, which are placed in different compartments. A rough sketch of the network is shown in A, with the MEV subnetwork on the left and the MEP subnetwork isolated in another compartment on the right. A clearer sketch of the MEV subnetwork is shown in B, and the MEP subnetwork is shown in C. The identified enzymatic and crosstalk mechanisms are shown in the sketches.

Figure 3
figure 3

The overall layout of the terpenoid biosynthetic pathway model. All metabolites are represented using continous places whereas the regulatory components such as the on/off switches and parameter modulators are represented using generic places. Each transition entity represents a process such as phosphorylation, synthesis and degradation. Each substrate or intermediate is connected to the transition entity by a normal arc whereas an enzyme is connected to its transition entity using a test arc. The MEV subnetwork is on the left and the MEP subnetwork is secluded in a different compartment to depict that the pathway is operative in plastids of eubacterias or chloroplasts or plants.

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.

Table 1 Reactions involved in the MEV pathway (in order of occurrence)
Table 2 Reactions involved in the MEP pathway (in order of occurrence)
Figure 4
figure 4

The MEV and MEP pathway constructed using Cell Illustrator. The pathways are numbered according to sequence, as shown in Tables 1 and 2.

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:

  1. 1)

    normal production of sesquiterpenes through the MEV subnetwork, with the MEP subnetwork and all inhibitory processes deactivated;

  2. 2)

    normal production of monoterpenes through the MEP subnetwork, with the MEV subnetwork and all inhibitory processes deactivated;

  3. 3)

    overexpression of ATP in the fifth reaction of the MEV subnetwork, with the MEP subnetwork and other inhibitory processes deactivated;

  4. 4)

    overproduction of FPP in the ninth reaction of the MEV subnetwork, with the MEP subnetwork and other inhibitory processes deactivated;

  5. 5)

    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 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 catalyses the production of mevalonate diphosphate (MEVPP) and ADP from mevalonate-5-phosphate (MEVP) and ATP (Table 1) [4648]. 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 [46]. 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) [46].

Figure 5
figure 5

Concentration (unit) changes in metabolites involved in the feedback loop mechanism over time (pt). Oscillations between the concentrations describe how MEVPP production halts when high concentrations of ATP are present and resumes when the concentration of MEVP increases.

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; EC2.7.1.36). 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).

Figure 6
figure 6

Concentration (unit) changes in metabolites involved in the feedback loop mechanism caused by overproduction of FPP over time (pt). The simulation successfully describes how high concentrations of FPP affect the production of MEVP. The oscillatory behaviour also affects the production of sesquiterpenes.

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).

c. Crosstalk

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; EC1.1.1.267). 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) [51]. 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.'

Figure 7
figure 7

Concentration (unit) changes in metabolites affected by the crosstalk mechanism over time (pt). When the concentration of fosmidomycin increases over time, monoterpene production decelerates at 100 pt; however, a slight increase in sesquiterpene concentration is observed.

Figure 8
figure 8

Concentration (unit) changes in levels of the main precursors involved in the crosstalk mechanism between both pathways over time (pt). IPP_2 and DMAPP_2 are precursors from the MEP subnetworks, whereas GPP and FPP serve as precursors in the MEV subnetwork. The transfer of fluxes from IPP_2 and DMAPP_2 to the MEV subnetwork enables GPP and FPP to be synthesised and thus produce sesquiterpenes.

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, 5256].

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, 2225].

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.



(adenosine diphosphate)


(adenosine triphosphate)


(dimethylallyl diphosphate)


(deoxyxylulose-5-phosphate reductoisomerase)


(enzyme commission)


(geranyl diphosphate)


(Hybrid Dynamic Petri net)


(Hybrid Functional Petri net with extension)


(Hybrid Petri net)


(isopentenyl diphosphate)


(Functional Petri net)


(farnesyl diphosphate)


(Petri net time)








(mevalonate diphosphate)


(mevalonate kinase)


(ordinary differential equation)


(phosphomevalonate kinase).


  1. Kitano H: Systems Biology: A Brief Overview. Science 2002, 295: 1662–1664. 10.1126/science.1069492

    Article  CAS  PubMed  Google Scholar 

  2. 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

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  3. 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

    Article  CAS  PubMed  Google Scholar 

  4. 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

    Article  CAS  PubMed  Google Scholar 

  5. 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

    Article  CAS  PubMed  Google Scholar 

  6. 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

    Article  CAS  Google Scholar 

  7. 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

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  8. 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.

    Article  CAS  Google Scholar 

  9. 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.

    CAS  PubMed  Google Scholar 

  10. 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.

    CAS  PubMed  Google Scholar 

  11. 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

    Article  CAS  PubMed  Google Scholar 

  12. 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

    Article  CAS  PubMed  Google Scholar 

  13. 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

    Chapter  Google Scholar 

  14. 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

    Article  CAS  PubMed  Google Scholar 

  15. 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

    Article  PubMed  PubMed Central  Google Scholar 

  16. 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

    Google Scholar 

  17. Nagasaki M: A platform for biopathway modeling/simulation and recreating biopathway databases towards simulation. In PhD Thesis. University of Tokyo; 2004.

    Google Scholar 

  18. 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.

    CAS  PubMed  Google Scholar 

  19. 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.

    Google Scholar 

  20. Nagasaki M, Saito A, Doi A, Matsuno H, Miyano S: Foundations of Systems Biology. London: Springer; 2009.

    Chapter  Google Scholar 

  21. Matsuno H, Doi A, Miyano S: Hybrid Petri net Representation of Gene Regulatory Network. Pacific Symposium on Biocomputing: 2000 2000, 338–349.

    Google Scholar 

  22. 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):

  23. 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.

    CAS  PubMed  Google Scholar 

  24. 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

    Article  Google Scholar 

  25. 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.

    Google Scholar 

  26. 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.

    Google Scholar 

  27. 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

    Article  CAS  PubMed  Google Scholar 

  28. 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

    Article  PubMed  PubMed Central  Google Scholar 

  29. 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

    Article  CAS  PubMed  Google Scholar 

  30. 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):

  31. 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.

    CAS  PubMed  Google Scholar 

  32. 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

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  33. 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):

  34. 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

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  35. 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

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  36. 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

    Article  CAS  PubMed  Google Scholar 

  37. 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

    Article  Google Scholar 

  38. Hunter WN: The non-mevalonate pathway of isoprenoid precursor biosynthesis. Journal of Biological Chemistry 2007, 282(30):21573–21577. 10.1074/jbc.R700005200

    Article  CAS  PubMed  Google Scholar 

  39. 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

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  40. 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

    Article  CAS  Google Scholar 

  41. 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

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  42. 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

    Article  CAS  PubMed  Google Scholar 

  43. 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

    Article  CAS  PubMed  Google Scholar 

  44. 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

    Article  CAS  PubMed  Google Scholar 

  45. 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

    Article  CAS  Google Scholar 

  46. Eyzaguirre J, Valdebenito D, Cardemil E: Pig Liver Phosphomevalonate kinase: Kinetic Mechanism. Archives of Biochemistry and Biophysics 2006, 454: 188–195. 10.1016/

    Article  Google Scholar 

  47. Doun SS, Burgner JW II, Briggs SD, Rodwell VW: Enterococcus faecalis phosphomevalonate kinase. Protein Science 2005, 14: 1134–1139. 10.1110/ps.041210405

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  48. 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

    Article  CAS  PubMed  Google Scholar 

  49. 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

    Article  CAS  Google Scholar 

  50. Gray JC, Kekwick RGO: Mevalonate Kinase in Green Leaves and Etiolated Cotyledons of the French Bean Phaseolus vulgaris . Biochemical Journal 1973, 133: 335–347.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  51. 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

    Article  CAS  PubMed  Google Scholar 

  52. 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.

    Google Scholar 

  53. 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

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  54. 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

    Google Scholar 

  55. 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

    Google Scholar 

  56. 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.

    Google Scholar 

Download references


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.

Author information

Authors and Affiliations


Corresponding author

Correspondence to Zeti-Azura Mohamed-Hussein.

Additional information

Authors' contributions

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

Rights and permissions

This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Reprints and permissions

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).

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: