Abstract
Hybrid Petri net (HPN) is an extension of the Petri net formalism, which enables us to handle continuous information in addition to discrete information. Firstly, this paper demonstrates how biological pathways can be modeled by the integration of discrete and continuous elements, with an example of the λ phage genetic switch system including induction and retroregulation mechanisms. Although HPN allows intuitive modeling of biological pathways, some fundamental biological processes such as complex formation cannot be represented with HPN. Thus, this paper next provides the formal definition of hybrid functional Petri net with extension (HFPNe), which has high potential for modeling various kinds of biological processes. Cell Illustrator is a software tool developed on the basis of the definition of HFPNe. Hypothesis creation by Cell Illustrator is demonstrated with the example of the cyanobacterial circadian gene clock system. Finally, our ongoing tasks, which include the development of a computational platform for systems biology, are presented.
Similar content being viewed by others
Notes
The HPN model in (Matsuno et al. 2000) employs the firing rule “the firing rate of a continuous transition can be expressed as a function of the input places of the continuous transition,” which was defined in the notion of hybrid dynamic net in (Drath et al. 1999). The incorporation of this firing rule in the HPN allows us to represent the rate of a reaction in a biological pathway as a function of the substances of the reaction. This representation has been commonly used in the dynamic modeling of biological pathways [see e.g. Gilbert and Heiner (2006)].
For example, discrete places P L and P R in Fig. 5 are used to represent the promoter status of operons starting from N and Cro genes, respectively. The status of the promoter, i.e., the presence of a token in the place P L or P R in this case, is determined by the token provision from its adjacent transition which is further controlled by the content(s) of its upstream place(s).
Figure 5 involves two operons starting from discrete places P L and P R , which are comprised of the sequences of N-CIII-Xis-Int B and Cro-CII-O-P-Q, respectively.
For the actual values of thresholds in Fig. 5 including this th 1, refer to the CSML file of Cell Illustrator (see subsection 5.2) that can be download from the website (http://www.genome.ib.sci.yamaguchi-u.ac.jp/~matsuno/lambda_CI.zip).
References
Alla H, David R (1998) Continuous and hybrid Petri nets. J Circuits Syst Comput 8:159–188
Drath R, Engmann U, Schwuchow S (1999) Hybrid aspects of modelling manufacturing systems using modified Petri nets. Proceedings of the 5th workshop on intelligent manufacturing systems, Granado, Brasil
Gilbert D, Heiner M (2006) From Petri nets to differential equations—an integrative approach for biochemical network analysis. 27th International conference on application and theory of Petri nets (LNCS4024), pp 181–200
Hardy S, Robillard PN (2008) Petri net-based method for the analysis of the dynamics of signal propagation in signaling pathways. Bioinformatics 24(2):209–217
Ishiura M, Kutsuna S, Aoki S, Iwasaki H, Andersson CR, Tanabe A, Golden SS, Johnson CH, Kondo T (1998) Expression of a gene cluster kaiABC as a circadian feedback process in cyanobacteria. Science 281:1519–1523
Iwasaki H, Williams SB, Kitayama Y, Ishiura M, Golden SS, Kondo T (2000) A KaiC-interacting sensory histidine kinase, SasA, necessary to sustain robust circadian oscillation in cyanobacteria. Cell 101:223–233
Iwasaki H, Nishiwaki T, Kitayama Y, Nakajima M, Kondo T (2002) KaiA-stimulated KaiC phosphorylation in circadian timing loops in cyanobacteria. Proc Natl Acad Sci USA 99(24):15788–15793
Jeong E, Nagasai M, Saito A, Miyano S (2007a) Cell System Ontology: representation for modeling, visualizing, and simulating biological pathways. In Silico Biol 7:0055
Jeong E, Nagasaki M, Miyano S (2007b) Conversion from BioPAX to CSO for system dynamics and visualization of biological pathway. Genome Inform 18:225–236
Kitayama Y, Iwasaki H, Nishiwaki T, Kondo T (2003) KaiB functions as an attenuator of KaiC phosphorylation in the cyanobacterial circadian clock system. EMBO J 22(9):2127–2134
Koh G, Teong HFC, Clément MV, Hsu D, Thiagarajan PS (2006) A decompositional approach to parameter estimation in pathway modeling: a case study of the Akt and MAPK pathways and their crosstalk. Bioinformatics 22(14):e271–e280
Kojima K, Nagasaki M, Jeong E, Kato M, Miyano S (2007) An efficient grid layout algorithm for biological networks utilizing various biological attributes. BMC Bioinform 8:76
Matsuno H, Doi A, Nagasaki M, Miyano S (2000) Hybrid Petri net representation of gene regulatory network. Pac Symp Biocomput 5:338–349 http://www.psb.stanford.edu/psb-online/proceedings/psb00/matsuno.pdf
Matsuno H, Tanaka Y, Aoshima H, Doi A, Matsui M, Miyano S (2003a) Biopathways representation and simulation on hybrid functional Petri net. In Silico Biol 3(3):389–404
Matsuno H, Fujita S, Doi A, Nagasaki M, Miyano S (2003b) Towards biopathway modeling and simulation, vol 2679. In: Proceedings of the 24th international conference on applications and theory of Petri nets (ICATPN2003), Lecture Notes in Computer Science, pp 3–22
Nagasaki M, Doi A, Matsuno H, Miyano S (2003) Genomic object net: I. A platform for modeling and simulating biopathways. Appl Bioinform 2(3):181–184
Nagasaki M, Doi A, Matsuno H, Miyano S (2004) A versatile Petri net based architecture for modeling and simulation of complex biological processes. Genome Inform 15(1):180–197
Nagasaki M, Doi A, Matsuno H, Miyano S (2005) Computational modeling of biological processes with Petri net based architecture. In: Chen YP (ed) Bioinformatics technologies. Springer Press, Berlin, pp 179–243
Nagasaki M, Yamaguchi R, Yoshida R, Imoto S, Doi A, Tamada Y, Matsuno H, Miyano S, Higuchi T (2006) Genomic data assimilation for estimating hybrid functional Petri net from time-course gene expression data. Genome Inform 17(1):46–61
Nagasaki M, Saito A, Li C, Jeong E, Miyano S (2008) Systematic reconstruction of TRANSPATH data into cell system markup language. BMC System Biol 2:53
Ptashne M (1992) A genetic switch, 2nd edn. Blackwell Science & Cell Press, MA
Tasaki S, Nagasaki M, Oyama M, Hata H, Ueno K, Yoshida R, Higuchi T, Sugano S, Miyano S (2006) Modeling and estimation of dynamic EGFR pathway by data assimilation approach using time series proteomic data. Genome Inform 17(2):226–228
Troncale S, Tahi F, Campard D, Vannier JP, Guespin J (2006) Modeling and simulation with hybrid functional Petri nets of the role of Interleukin-6 in human early haematopoiesis. Pacific Symp Biocomput 11:427–438
Acknowledgements
Hiroshi Matsuno, one of the authors of this paper, would like to express his gratitude to Dr. Atsushi Doi at Cell Innovator Inc., and Ms. Mao Yamaji at CORE Corporation who contributed to the constructions of the hybrid functional Petri net models of λ phage genetic switch and cyanobacteria circadian rhythm, respectively. This work is partially supported by Grant-in-Aid for Scientific Research on Priority Areas “Systems Genomics”and Grant-in-Aid for Scientific Research (B) from the Ministry of Education, Culture, Sports, Science, and Technology of Japan.
Author information
Authors and Affiliations
Corresponding author
Rights and permissions
About this article
Cite this article
Matsuno, H., Nagasaki, M. & Miyano, S. Hybrid Petri net based modeling for biological pathway simulation. Nat Comput 10, 1099–1120 (2011). https://doi.org/10.1007/s11047-009-9164-6
Published:
Issue Date:
DOI: https://doi.org/10.1007/s11047-009-9164-6