Programming with models: modularity and abstraction provide powerful capabilities for systems biology

Mathematical models are increasingly used to understand how phenotypes emerge from systems of molecular interactions. However, their current construction as monolithic sets of equations presents a fundamental barrier to progress. Overcoming this requires modularity, enabling sub-systems to be specified independently and combined incrementally, and abstraction, enabling generic properties of biological processes to be specified independently of specific instances. These, in turn, require models to be represented as programs rather than as datatypes. Programmable modularity and abstraction enables libraries of modules to be created, which can be instantiated and reused repeatedly in different contexts with different components. We have developed a computational infrastructure that accomplishes this. We show here why such capabilities are needed, what is required to implement them and what can be accomplished with them that could not be done previously.

[1]  Leslie Greengard,et al.  A mathematical tool for exploring the dynamics of biological networks , 2007, Proceedings of the National Academy of Sciences.

[2]  David A. Moon,et al.  The philosophy of Lisp , 1991, CACM.

[3]  Andreas Kremling,et al.  Modular Modeling of Cellular Systems with ProMoT/Diva , 2003, Bioinform..

[4]  P. Lawrence The making of a fly , 1992 .

[5]  P. Roach,et al.  Multisite and hierarchal protein phosphorylation. , 1991, The Journal of biological chemistry.

[6]  James M. Bower,et al.  The Book of GENESIS , 1994, Springer New York.

[7]  Kevin N. Dalby,et al.  The Kinetic Mechanism of the Dual Phosphorylation of the ATF2 Transcription Factor by p38 Mitogen-activated Protein (MAP) Kinase α , 2001, The Journal of Biological Chemistry.

[8]  金明桓,et al.  Multiprocessor 컴퓨터 시스템 , 1984 .

[9]  Hiroaki Kitano,et al.  The systems biology markup language (SBML): a medium for representation and exchange of biochemical network models , 2003, Bioinform..

[10]  D. Sachs,et al.  Reconstitution with syngeneic plus allogeneic or xenogeneic bone marrow leads to specific acceptance of allografts or xenografts , 1984, Nature.

[11]  J. Gunawardena Multisite protein phosphorylation makes a good threshold but can be a poor switch. , 2005, Proceedings of the National Academy of Sciences of the United States of America.

[12]  Michael Hucka,et al.  The Systems Biology Markup Language (SBML) Level 2 Version 2 , 2007 .

[13]  Lea Sistonen,et al.  Multisite phosphorylation provides sophisticated regulation of transcription factors. , 2002, Trends in biochemical sciences.

[14]  William J. R. Longabaugh,et al.  Computational representation of developmental genetic regulatory networks. , 2005, Developmental biology.

[15]  Gaudenz Danuser,et al.  Linking data to models: data regression , 2006, Nature Reviews Molecular Cell Biology.

[16]  Lewis Wolpert,et al.  Principles of Development , 1997 .

[17]  Hugh D. Spence,et al.  Minimum information requested in the annotation of biochemical models (MIRIAM) , 2005, Nature Biotechnology.

[18]  D. Rand,et al.  Uncovering the design principles of circadian clocks: mathematical analysis of flexibility and evolutionary goals. , 2006, Journal of theoretical biology.

[19]  D. Lauffenburger,et al.  Physicochemical modelling of cell signalling pathways , 2006, Nature Cell Biology.

[20]  James M. Bower,et al.  Neural Modeling with GENESIS , 1998 .

[21]  H. Kitano,et al.  Computational systems biology , 2002, Nature.

[22]  Eduardo Sontag,et al.  Untangling the wires: A strategy to trace functional interactions in signaling and gene networks , 2002, Proceedings of the National Academy of Sciences of the United States of America.

[23]  Yu Zhao,et al.  The Mechanism of Dephosphorylation of Extracellular Signal-regulated Kinase 2 by Mitogen-activated Protein Kinase Phosphatase 3* , 2001, The Journal of Biological Chemistry.

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

[25]  G. Jicha,et al.  Hierarchical Phosphorylation of Recombinant Tau by the Paired‐Helical Filament‐Associated Protein Kinase Is Dependent on Cyclic AMP‐Dependent Protein Kinase , 1999, Journal of neurochemistry.

[26]  Mark A. Miller,et al.  A Dynamic Model of Guard Cell Abscisic Acid Signaling , 2006 .

[27]  P. Cohen The role of protein phosphorylation in human health and disease. The Sir Hans Krebs Medal Lecture. , 2001, European journal of biochemistry.

[28]  Xiang-Dong Fu,et al.  Mass Spectrometric and Kinetic Analysis of ASF/SF2 Phosphorylation by SRPK1 and Clk/Sty* , 2005, Journal of Biological Chemistry.

[29]  Charles L. Forgy,et al.  Rete: a fast algorithm for the many pattern/many object pattern match problem , 1991 .

[30]  Cristina M Furdui,et al.  Autophosphorylation of FGFR1 kinase is mediated by a sequential and precisely ordered reaction. , 2006, Molecular cell.

[31]  G. Odell,et al.  Design and constraints of the Drosophila segment polarity module: robust spatial patterning emerges from intertwined cell state switches. , 2002, The Journal of experimental zoology.

[32]  J. Weinstein,et al.  Molecular interaction maps of bioregulatory networks: a general rubric for systems biology. , 2005, Molecular biology of the cell.

[33]  J. Hopfield,et al.  From molecular to modular cell biology , 1999, Nature.

[34]  Hamid Bolouri,et al.  Dizzy: Stochastic Simulation of Large-scale Genetic Regulatory Networks , 2005, J. Bioinform. Comput. Biol..

[35]  Bard Ermentrout,et al.  Simulating, analyzing, and animating dynamical systems - a guide to XPPAUT for researchers and students , 2002, Software, environments, tools.

[36]  Vincent Danos,et al.  Scalable Simulation of Cellular Signaling Networks , 2007, APLAS.

[37]  Melanie I. Stefan,et al.  Molecules for memory: modelling CaMKII , 2007, BMC Systems Biology.

[38]  Vincent Danos,et al.  Abstract Interpretation of Cellular Signalling Networks , 2008, VMCAI.

[39]  Catherine M Lloyd,et al.  CellML: its future, present and past. , 2004, Progress in biophysics and molecular biology.

[40]  Kit-Yi Leung,et al.  Novel Phosphorylation Sites in Tau from Alzheimer Brain Support a Role for Casein Kinase 1 in Disease Pathogenesis* , 2007, Journal of Biological Chemistry.

[41]  Jehoshua Bruck,et al.  Scaffold proteins may biphasically affect the levels of mitogen-activated protein kinase signaling and reduce its threshold properties. , 2000, Proceedings of the National Academy of Sciences of the United States of America.

[42]  Bruce E. Shapiro,et al.  Cellerator: extending a computer algebra system to include biochemical arrows for signal transduction simulations , 2003, Bioinform..

[43]  U. Bhalla,et al.  Emergent properties of networks of biological signaling pathways. , 1999, Science.

[44]  Paul Graham ANSI Common Lisp , 1995 .

[45]  Roger Brent,et al.  Automatic generation of cellular reaction networks with Moleculizer 1.0 , 2005, Nature Biotechnology.

[46]  G. Odell,et al.  The segment polarity network is a robust developmental module , 2000, Nature.

[47]  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.

[48]  Peter D. Karp,et al.  Querying and computing with BioCyc databases , 2005, Bioinform..

[49]  H. Redkey,et al.  A new approach. , 1967, Rehabilitation record.

[50]  Stewart Shuman,et al.  Schizosaccharomyces pombe Carboxyl-terminal Domain (CTD) Phosphatase Fcp1 , 2004, Journal of Biological Chemistry.

[51]  Ronan M. T. Fleming,et al.  Quantitative prediction of cellular metabolism with constraint-based models: the COBRA Toolbox v2.0 , 2007, Nature Protocols.

[52]  Nicholas T Ingolia,et al.  Topology and Robustness in the Drosophila Segment Polarity Network , 2004, PLoS biology.

[53]  A. Bangs,et al.  Finding value in in silico biology , 2003 .

[54]  A. Arkin Synthetic cell biology. , 2001, Current opinion in biotechnology.

[55]  Yukiko Matsuoka,et al.  Using process diagrams for the graphical representation of biological networks , 2005, Nature Biotechnology.

[56]  Rainer Breitling,et al.  What is Systems Biology? , 2010, Front. Physiology.

[57]  T. Ideker,et al.  A new approach to decoding life: systems biology. , 2001, Annual review of genomics and human genetics.

[58]  Nicholas T. Carnevale,et al.  The NEURON Book: Epilogue , 2006 .

[59]  Philip R. Cohen,et al.  The development and therapeutic potential of protein kinase inhibitors. , 1999, Current opinion in chemical biology.

[60]  U. Alon,et al.  Robustness in bacterial chemotaxis , 2022 .

[61]  A. Harwood,et al.  Regulation of GSK-3 A Cellular Multiprocessor , 2001, Cell.

[62]  N. Perrimon,et al.  The emergence of geometric order in proliferating metazoan epithelia , 2006, Nature.

[63]  M. Feinberg The existence and uniqueness of steady states for a class of chemical reaction networks , 1995 .

[64]  Christopher R. Myers,et al.  Universally Sloppy Parameter Sensitivities in Systems Biology Models , 2007, PLoS Comput. Biol..

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

[66]  J. Gunawardena,et al.  Unlimited multistability in multisite phosphorylation systems , 2009, Nature.

[67]  Tian Jin,et al.  Key Role of Local Regulation in Chemosensing Revealed by a New Molecular Interaction-Based Modeling Method , 2006, PLoS Comput. Biol..

[68]  Jacky L. Snoep,et al.  BioModels Database: a free, centralized database of curated, published, quantitative kinetic models of biochemical and cellular systems , 2005, Nucleic Acids Res..

[69]  W. R. Burack,et al.  The activating dual phosphorylation of MAPK by MEK is nonprocessive. , 1997, Biochemistry.

[70]  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.

[71]  P. Pellicena,et al.  Processive Phosphorylation of p130Cas by Src Depends on SH3-Polyproline Interactions* , 2001, The Journal of Biological Chemistry.

[72]  Xiang-Dong Fu,et al.  Processive phosphorylation of alternative splicing factor/splicing factor 2 , 2003, Proceedings of the National Academy of Sciences of the United States of America.

[73]  John J. Tyson,et al.  Modeling Networks of Coupled Enzymatic Reactions Using the Total Quasi-Steady State Approximation , 2007, PLoS Comput. Biol..

[74]  T. Henzinger,et al.  Executable cell biology , 2007, Nature Biotechnology.

[75]  E. O’Shea,et al.  Multi-site phosphorylation of Pho4 by the cyclin-CDK Pho80-Pho85 is semi-processive with site preference. , 2001, Journal of molecular biology.

[76]  G J Griffiths,et al.  Decreased internalisation of erbB1 mutants in lung cancer is linked with a mechanism conferring sensitivity to gefitinib. , 2006, Systems biology.

[77]  Michael Travers,et al.  BioLingua: a programmable knowledge environment for biologists , 2005, Bioinform..

[78]  Leslie M Loew,et al.  Computational cell biology: spatiotemporal simulation of cellular events. , 2002, Annual review of biophysics and biomolecular structure.

[79]  S. V. Aksenov,et al.  A systems biology dynamical model of mammalian G1 cell cycle progression , 2007, Molecular systems biology.

[80]  Kevin N. Dalby,et al.  The kinetic mechanism of the dual phosphorylation of the ATF2 transcription factor by p38 MAP kinaseα: Implications for signal/response Profiles of MAP Kinase Pathways , 2000 .

[81]  G. Odell,et al.  Ingeneue: a versatile tool for reconstituting genetic networks, with examples from the segment polarity network. , 2002, The Journal of experimental zoology.

[82]  Tony White,et al.  UML as a cell and biochemistry modeling language. , 2005, Bio Systems.

[83]  James E. Ferrell,et al.  Mechanistic Studies of the Dual Phosphorylation of Mitogen-activated Protein Kinase* , 1997, The Journal of Biological Chemistry.

[84]  D. Morrison,et al.  Regulation of MAP kinase signaling modules by scaffold proteins in mammals. , 2003, Annual review of cell and developmental biology.

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

[86]  T. Bartol,et al.  Monte Carlo Methods for Simulating Realistic Synaptic Microphysiology Using MCell , 2000 .

[87]  Marc W Kirschner,et al.  The Meaning of Systems Biology , 2005, Cell.

[88]  Charles L. Forgy,et al.  Rete: A Fast Algorithm for the Many Patterns/Many Objects Match Problem , 1982, Artif. Intell..