Rule-Based Modelling of Cellular Signalling

Modelling is becoming a necessity in studying biological signalling pathways, because the combinatorial complexity of such systems rapidly overwhelms intuitive and qualitative forms of reasoning. Yet, this same combinatorial explosion makes the traditional modelling paradigm based on systems of differential equations impractical. In contrast, agentbased or concurrent languages, such as ? [1,2,3] or the closely related BioNetGen language [4,5,6,7,8,9,10], describe biological interactions in terms of rules, thereby avoiding the combinatorial explosion besetting differential equations. Rules are expressed in an intuitive graphical form that transparently represents biological knowledge. In this way, rules become a natural unit of model building, modification, and discussion. We illustrate this with a sizeable example obtained from refactoring two models of EGF receptor signalling that are based on differential equations [11,12]. An exciting aspect of the agent-based approach is that it naturally lends itself to the identification and analysis of the causal structures that deeply shape the dynamical, and perhaps even evolutionary, characteristics of complex distributed biological systems. In particular, one can adapt the notions of causality and conflict, familiar from concurrency theory, to ?, our representation language of choice. Using the EGF receptor model as an example, we show how causality enables the formalization of the colloquial concept of pathway and, perhaps more surprisingly, how conflict can be used to dissect the signalling dynamics to obtain a qualitative handle on the range of system behaviours. By taming the combinatorial explosion, and exposing the causal structures and key kinetic junctures in a model, agent- and rule-based representations hold promise for making modelling more powerful, more perspicuous, and of appeal to a wider audience.

[1]  Ugo Montanari,et al.  Unfolding and Event Structure Semantics for Graph Grammars , 1999, FoSSaCS.

[2]  Chi-Ying F. Huang,et al.  Ultrasensitivity in the mitogen-activated protein kinase cascade. , 1996, Proceedings of the National Academy of Sciences of the United States of America.

[3]  William S. Hlavacek,et al.  BioNetGen: software for rule-based modeling of signal transduction based on the interactions of molecular domains , 2004, Bioinform..

[4]  H. Kitano,et al.  A comprehensive pathway map of epidermal growth factor receptor signaling , 2005, Molecular systems biology.

[5]  Corrado Priami,et al.  Application of a stochastic name-passing calculus to representation and simulation of molecular processes , 2001, Inf. Process. Lett..

[6]  Vincent Danos,et al.  Computational self-assembly , 2008, Theor. Comput. Sci..

[7]  Vincent Danos,et al.  Formal Molecular Biology Done in CCS-R , 2007, Electron. Notes Theor. Comput. Sci..

[8]  Pierpaolo Degano,et al.  Causal pi-Calculus for Biochemical Modelling , 2003, CMSB.

[9]  H. Lane,et al.  ERBB Receptors and Cancer: The Complexity of Targeted Inhibitors , 2005, Nature Reviews Cancer.

[10]  R. Milner,et al.  Bigraphical Reactive Systems , 2001, CONCUR.

[11]  R. Weinberg,et al.  The Biology of Cancer , 2006 .

[12]  Luca Cardelli,et al.  Brane Calculi , 2004, CMSB.

[13]  T. Pawson,et al.  Assembly of Cell Regulatory Systems Through Protein Interaction Domains , 2003, Science.

[14]  Cosimo Laneve,et al.  Core Formal Molecular Biology , 2003, ESOP.

[15]  D. Koshland,et al.  An amplified sensitivity arising from covalent modification in biological systems. , 1981, Proceedings of the National Academy of Sciences of the United States of America.

[16]  R. Bresalier,et al.  The Biology of Cancer. Garland Science, Oxford, UK (2006), 864 pp. $104.95, ISBN: 0-8153-4076-1. Web site for ordering: www.garandscience.com , 2007 .

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

[18]  Michael Hucka,et al.  A Correction to the Review Titled "Rules for Modeling Signal-Transduction Systems" by W. S. Hlavacek et al. , 2006, Science's STKE.

[19]  Ehud Shapiro,et al.  Cells as Computation , 2003, CMSB.

[20]  Cosimo Laneve,et al.  Formal molecular biology , 2004, Theor. Comput. Sci..

[21]  Aviv Regev,et al.  Representation and Simulation of Biochemical Processes Using the pi-Calculus Process Algebra , 2000, Pacific Symposium on Biocomputing.

[22]  D. Gillespie Exact Stochastic Simulation of Coupled Chemical Reactions , 1977 .

[23]  Vincent Danos,et al.  Formal Molecular Biology done in CCS , 2003 .

[24]  Jin Yang,et al.  Graph Theory for Rule-Based Modeling of Biochemical Networks , 2006, Trans. Comp. Sys. Biology.

[25]  E. Gilles,et al.  Computational modeling of the dynamics of the MAP kinase cascade activated by surface and internalized EGF receptors , 2002, Nature Biotechnology.

[26]  Corrado Priami,et al.  Beta Binders for Biological Interactions , 2004, CMSB.

[27]  D. Fell,et al.  Differential feedback regulation of the MAPK cascade underlies the quantitative differences in EGF and NGF signalling in PC12 cells , 2000, FEBS letters.

[28]  William S. Hlavacek,et al.  Graphical rule-based representation of signal-transduction networks , 2005, SAC '05.

[29]  B. Kholodenko,et al.  Quantification of Short Term Signaling by the Epidermal Growth Factor Receptor* , 1999, The Journal of Biological Chemistry.

[30]  Glynn Winskel,et al.  An introduction to event structures , 1988, REX Workshop.

[31]  E. Shapiro,et al.  Cellular abstractions: Cells as computation , 2002, Nature.

[32]  José Meseguer,et al.  Pathway Logic: Symbolic Analysis of Biological Signaling , 2001, Pacific Symposium on Biocomputing.

[33]  Corrado Priami,et al.  Modelling biochemical pathways through enhanced pi-calculus , 2004, Theor. Comput. Sci..

[34]  Glynn Winskel,et al.  Petri Nets, Event Structures and Domains , 1979, Semantics of Concurrent Computation.

[35]  Glynn Winskel,et al.  Petri Nets, Event Structures and Domains, Part I , 1981, Theor. Comput. Sci..

[36]  W. S. Hlavacek,et al.  A network model of early events in epidermal growth factor receptor signaling that accounts for combinatorial complexity. , 2006, Bio Systems.

[37]  Luca Cardelli,et al.  BioAmbients: an abstraction for biological compartments , 2004, Theor. Comput. Sci..

[38]  B. Kholodenko,et al.  Negative feedback and ultrasensitivity can bring about oscillations in the mitogen-activated protein kinase cascades. , 2000, European journal of biochemistry.

[39]  Mogens Nielsen,et al.  Models for Concurrency , 1992 .

[40]  Davide Sangiorgi,et al.  Communicating and Mobile Systems: the π-calculus, , 2000 .

[41]  Jin Yang,et al.  Depicting signaling cascades , 2006, Nature Biotechnology.

[42]  Robin Milner,et al.  Communicating and mobile systems - the Pi-calculus , 1999 .

[43]  M L Blinov,et al.  Combinatorial complexity and dynamical restriction of network flows in signal transduction. , 2004, Systems biology.