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

BackgroundThe 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).ResultsThe 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.ConclusionsAlthough 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.

[1]  Monika Heiner,et al.  Application of Petri net based analysis techniques to signal transduction pathways , 2006, BMC Bioinformatics.

[2]  Ina Koch,et al.  Petri net modelling of gene regulation of the Duchenne muscular dystrophy , 2008, Biosyst..

[3]  Luhua Lai,et al.  Dynamic Simulations on the Arachidonic Acid Metabolic Network , 2007, PLoS Comput. Biol..

[4]  Pedro Mendes,et al.  GEPASI: a software package for modelling the dynamics, steady states and control of biochemical and other systems , 1993, Comput. Appl. Biosci..

[5]  Kiyoko F. Aoki-Kinoshita,et al.  From genomics to chemical genomics: new developments in KEGG , 2005, Nucleic Acids Res..

[6]  R. Croteau,et al.  Isopentenyl diphosphate is the terminal product of the deoxyxylulose-5-phosphate pathway for terpenoid biosynthesis in plants , 1999 .

[7]  William N. Hunter,et al.  The Non-mevalonate Pathway of Isoprenoid Precursor Biosynthesis* , 2007, Journal of Biological Chemistry.

[8]  Masao Nagasaki,et al.  A versatile petri net based architecture for modeling and simulation of complex biological processes. , 2004, Genome informatics. International Conference on Genome Informatics.

[9]  M. Chye,et al.  Brassica juncea 3-hydroxy-3-methylglutaryl (HMG)-CoA synthase 1: expression and characterization of recombinant wild-type and mutant enzymes. , 2004, The Biochemical journal.

[10]  Russ B. Altman,et al.  Modelling biological processes using workflow and Petri Net models , 2002, Bioinform..

[11]  Oliver E. Sturm,et al.  Computational modelling of the receptor-tyrosine-kinase-activated MAPK pathway. , 2005, The Biochemical journal.

[12]  H. Miziorko,et al.  Staphylococcus aureus Mevalonate Kinase: Isolation and Characterization of an Enzyme of the Isoprenoid Biosynthetic Pathway , 2004, Journal of bacteriology.

[13]  Olaf Wolkenhauer,et al.  SBML export interface for the systems biology toolbox for MATLAB , 2007, Bioinform..

[14]  Hiroshi Matsuno,et al.  How to Model and Simulate Biological Pathways with Petri Nets - A New Challenge for Systems Biology - ∗ , 2004 .

[15]  V. Dubey,et al.  An overview of the non-mevalonate pathway for terpenoid biosynthesis in plants , 2003, Journal of Biosciences.

[16]  J. Schwender,et al.  Cloning and heterologous expression of a cDNA encoding 1‐deoxy‐D‐xylulose‐5‐phosphate reductoisomerase of Arabidopsis thaliana 1 , 1999, FEBS letters.

[17]  B. M. Lange,et al.  A family of transketolases that directs isoprenoid biosynthesis via a mevalonate-independent pathway. , 1998, Proceedings of the National Academy of Sciences of the United States of America.

[18]  Masao Nagasaki,et al.  Cell fate simulation model of gustatory neurons with MicroRNAs double-negative feedback loop by hybrid functional Petri net with extension. , 2006, Genome informatics. International Conference on Genome Informatics.

[19]  Masao Nagasaki,et al.  Constructing biological pathway models with hybrid functional petri nets. , 2011, Studies in health technology and informatics.

[20]  Ronald J A Wanders,et al.  Phosphomevalonate kinase is a cytosolic protein in humans Published, JLR Papers in Press, January 16, 2004. DOI 10.1194/jlr.M300373-JLR200 , 2004, Journal of Lipid Research.

[21]  R. Verpoorte,et al.  Purification and characterization of phosphomevalonate kinase from Catharanthus roseus , 1999 .

[22]  Lawrence Hunter,et al.  Pacific symposium on biocomputing 2006 , 2005, PSB 2016.

[23]  Russ B. Altman,et al.  Research Paper: Using Petri Net Tools to Study Properties and Dynamics of Biological Systems , 2004, J. Am. Medical Informatics Assoc..

[24]  Atsushi Doi,et al.  Biopathways representation and simulation on hybrid functional Petri net , 2003, Silico Biol..

[25]  Antje Chang,et al.  BRENDA, AMENDA and FRENDA the enzyme information system: new content and tools in 2009 , 2008, Nucleic Acids Res..

[26]  W. Eisenreich,et al.  The non-mevalonate pathway of isoprenoids: genes, enzymes and intermediates. , 2001, Current opinion in chemical biology.

[27]  Monika Heiner,et al.  Application of Petri net theory for modelling and validation of the sucrose breakdown pathway in the potato tuber , 2005, Bioinform..

[28]  Alisdair R Fernie,et al.  Regulation of metabolic networks: understanding metabolic complexity in the systems biology era. , 2005, The New phytologist.

[29]  W. Eisenreich,et al.  Studies on the nonmevalonate terpene biosynthetic pathway: Metabolic role of IspH (LytB) protein , 2002, Proceedings of the National Academy of Sciences of the United States of America.

[30]  Yoshihiro Yamanishi,et al.  KEGG for linking genomes to life and the environment , 2007, Nucleic Acids Res..

[31]  J. Eyzaguirre,et al.  Pig liver phosphomevalonate kinase: kinetic mechanism. , 2006, Archives of biochemistry and biophysics.

[32]  H Matsuno,et al.  Hybrid Petri net representation of gene regulatory network. , 1999, Pacific Symposium on Biocomputing. Pacific Symposium on Biocomputing.

[33]  Masao Nagasaki,et al.  Simulation-Based Validation of the p53 Transcriptional Activity with Hybrid Functional Petri Net , 2006, Silico Biol..

[34]  Masao Nagasaki,et al.  Simulation-based model checking approach to cell fate specification during Caenorhabditis elegans vulval development by hybrid functional Petri net with extension , 2009, BMC Systems Biology.

[35]  Alisdair R Fernie,et al.  Predictive Metabolic Engineering: A Goal for Systems Biology1 , 2003, Plant Physiology.

[36]  M. Rohmer,et al.  Isoprenoid biosynthesis via the methylerythritol phosphate pathway. Mechanistic investigations of the 1-deoxy-D-xylulose 5-phosphate reductoisomerase. , 2002, European journal of biochemistry.

[37]  Hiroshi Matsuno,et al.  Petri Net Based Descriptions for Systematic Understanding of Biological Pathways , 2006, IEICE Trans. Fundam. Electron. Commun. Comput. Sci..

[38]  H. Kitano Systems Biology: A Brief Overview , 2002, Science.

[39]  W. Eisenreich,et al.  Deoxyxylulose phosphate pathway to terpenoids. , 2001, Trends in plant science.

[40]  J. Gray,et al.  Mevalonate kinase in green leaves and etiolated cotyledons of the french bean Phaseolus vulgaris. , 1973, The Biochemical journal.

[41]  Janet B. Jones-Oliveira,et al.  A Computational Model for the Identification of Biochemical Pathways in the Krebs Cycle , 2003, J. Comput. Biol..

[42]  Shunji Takahashi,et al.  Cloning and Characterization of 1-Deoxy-d-Xylulose 5-Phosphate Synthase fromStreptomyces sp. Strain CL190, Which Uses both the Mevalonate and Nonmevalonate Pathways for Isopentenyl Diphosphate Biosynthesis , 2000, Journal of bacteriology.

[43]  Hiroyuki Ogata,et al.  KEGG: Kyoto Encyclopedia of Genes and Genomes , 1999, Nucleic Acids Res..

[44]  長崎 正朗,et al.  A platform for biopathway modeling/simulation and recreating biopathway databases towards simulation , 2004 .

[45]  Michael K. Molloy,et al.  Petri net , 2003 .

[46]  B. Blagg,et al.  1-Deoxy-d-Xylulose 5-Phosphate Synthase, the Gene Product of Open Reading Frame (ORF) 2816 and ORF 2895 inRhodobacter capsulatus , 2001, Journal of bacteriology.

[47]  Shunji Takahashi,et al.  Characterization of 1-Deoxy-d-xylulose 5-Phosphate Reductoisomerase, an Enzyme Involved in Isopentenyl Diphosphate Biosynthesis, and Identification of Its Catalytic Amino Acid Residues* , 2000, The Journal of Biological Chemistry.

[48]  S. Briggs,et al.  Enterococcus faecalis phosphomevalonate kinase , 2005, Protein science : a publication of the Protein Society.

[49]  Richard Banks,et al.  Qualitatively modelling and analysing genetic regulatory networks: a Petri net approach , 2007, Bioinform..

[50]  R. Verpoorte,et al.  Mevalonate kinase activity in Catharanthus roseus plants and suspension cultured cells , 2000 .

[51]  Igor Goryanin,et al.  Mathematical simulation and analysis of cellular metabolism and regulation , 1999, Bioinform..

[52]  Masao Nagasaki,et al.  Foundations of Systems Biology - Using Cell Illustrator® and Pathway Databases , 2009, Computational Biology.

[53]  David R. Gilbert,et al.  Computational methodologies for modelling, analysis and simulation of signalling networks , 2006, Briefings Bioinform..