Programming biological models in Python using PySB

Molecular Systems Biology

Volumn 9, Issue 1, 2013, Pages

  Lopez, Carlos F ^a   Muhlich, Jeremy L ^b   Bachman, John A ^b   Sorger, Peter K ^b  

^a VANDERBILT UNIVERSITY SCHOOL OF MEDICINE   (United States)

^b HARVARD MEDICAL SCHOOL   (United States)

Author keywords

apoptosis; modeling; rule based; software engineering

Indexed keywords

APOPTOSIS; ARTICLE; BIOLOGICAL MODEL; CATALYSIS; CHEMICAL ANALYSIS; CHEMICAL BINDING; COMPUTER PROGRAM; ECOSYSTEM; LANGUAGE; LIBRARY; POLYMERIZATION; PRIORITY JOURNAL; PYSB MODEL; PYTHON MODEL; COMPUTER LANGUAGE; COMPUTER SIMULATION; MITOCHONDRION; PHYSIOLOGY;

PROTEIN BCL 2;

APOPTOSIS; COMPUTER SIMULATION; MITOCHONDRIA; MODELS, BIOLOGICAL; PROGRAMMING LANGUAGES; PROTO-ONCOGENE PROTEINS C-BCL-2; SOFTWARE;

EID: 84882758675     PISSN: None     EISSN: 17444292     Source Type: Journal    
DOI: 10.1038/msb.2013.1     Document Type: Article

References (84)

1. Albeck JG, Burke JM, Aldridge BB, Zhang M, Lauffenburger DA, Sorger PK, (2008a) Quantitative analysis of pathways controlling extrinsic apoptosis in single cells. Mol Cell 30: 11-25
2. Albeck JG, Burke JM, Spencer SL, Lauffenburger DA, Sorger PK, (2008b) Modeling a snap-action, variable-delay switch controlling extrinsic cell death. PLoS Biol 6: 2831-2852
3. Aldridge BB, Gaudet S, Lauffenburger DA, Sorger PK, (2011) Lyapunov exponents and phase diagrams reveal multi-factorial control over TRAIL-induced apoptosis. Mol Syst Biol 7: 553
4. Bachman JA, Sorger P, (2011) New approaches to modeling complex biochemistry. Nat Meth 8: 130-131
5. Bagci E, Vodovotz Y, Billiar T, Ermentrout G, Bahar I, (2006) Bistability in apoptosis: roles of bax, bcl-2, and mitochondrial permeability transition pores. Biophys J 90: 1546-1559
6. Batchelor E, Mock CS, Bhan I, Loewer A, Lahav G, (2008) Recurrent initiation: a mechanism for triggering p53 pulses in response to DNA damage. Mol Cell 30: 277-289
7. Bentele M, Lavrik I, Ulrich M, Stösser S, Heermann DW, Kalthoff H, Krammer PH, Eils R, (2004) Mathematical modeling reveals threshold mechanism in CD95-induced apoptosis. J Cell Biol 166: 839-851
8. Billen LP, Kokoski CL, Lovell JF, Leber B, Andrews DW, (2008) Bcl-XL inhibits membrane permeabilization by competing with Bax. PLoS Biol 6: e147
9. Blinov ML, Faeder JR, Goldstein B, Hlavacek WS, (2006) A network model of early events in epidermal growth factor receptor signaling that accounts for combinatorial complexity. BioSystems 83: 136-151
10. Brown P, Hindmarsh A, Byrne G, (1989) VODE: a variable coefficient ODE solver. SIAM J Sci Stat Comput 10: 1038-1051
11. Certo M, Del Gaizo Moore V, Nishino M, Wei G, Korsmeyer S, Armstrong SA, Letai A, (2006) Mitochondria primed by death signals determine cellular addiction to antiapoptotic BCL-2 family members. Cancer Cell 9: 351-365
12. Chen C, Cui J, Lu H, Wang R, Zhang S, Shen P, (2007a) Modeling of the role of a Bax-activation switch in the mitochondrial apoptosis decision. Biophys J 92: 4304-4315
13. Chen C, Cui J, Zhang W, Shen P, (2007b) Robustness analysis identifies the plausible model of the Bcl-2 apoptotic switch. FEBS Lett 581: 5143-5150
14. Chen L, Willis SN, Wei A, Smith BJ, Fletcher JI, Hinds MG, Colman PM, Day CL, Adams JM, Huang DCS, (2005) Differential targeting of prosurvival Bcl-2 proteins by their BH3-only ligands allows complementary apoptotic function. Mol Cell 17: 393-403
15. Chen WW, Schoeberl B, Jasper PJ, Niepel M, Nielsen UB, Lauffenburger DA, Sorger PK, (2009) Input-output behavior of ErbB signaling pathways as revealed by a mass action model trained against dynamic data. Mol Syst Biol 5: 239
16. Chylek LA, Hu B, Blinov ML, Emonet T, Faeder JR, Goldstein B, Gutenkunst RN, Haugh JM, Lipniacki T, Posner RG, Yang J, Hlavacek WS, (2011) Guidelines for visualizing and annotating rule-based models. Mol Biosyst 7: 2779-2795
17. Cui J, Chen C, Lu H, Sun T, Shen P, (2008) Two independent positive feedbacks and bistability in the Bcl-2 apoptotic switch. PLoS ONE 3: e1469
18. Danos V, (2007) Agile modelling of cellular signalling. AIP Conference Proceedings 963: 611-614
19. Danos V, Feret J, Fontana W, Krivine J, (2007a) Scalable simulation of cellular signaling networks. Lect Notes Comput Sci 4807: 139-157
20. Danos V, Feret J, Fontana W, Krivine J, (2008) Abstract interpretation of cellular signalling networks. Lect Notes Comput Sci 4905: 83-97
21. Danos V, Feret J, Fontana W, Harmer R, Krivine J, (2007b) Rule-based modelling of cellular signalling. Lect Notes Comput Sci 4703: 17-41
22. Danos V, Feret J, Fontana W, Harmer R, Krivine J, (2009) Rule-based modelling and model perturbation. Lect Notes Comput Sci 5750: 116-137
23. Debartolo J, Dutta S, Reich L, Keating AE, (2012) Predictive bcl-2 family binding models rooted in experiment or structure. J Mol Biol 422: 124-144
24. Deeds EJ, Krivine J, Feret J, Danos V, Fontana W, (2012) Combinatorial complexity and compositional drift in protein interaction networks. PLoS ONE 7: e32032
25. Del Vecchio D, Ninfa AJ, Sontag ED, (2008) Modular cell biology: retroactivity and insulation. Mol Syst Biol 4: 161
26. Düssmann H, Rehm M, Concannon CG, Anguissola S, Würstle M, Kacmar S, Völler P, Huber HJ, Prehn JHM, (2010) Single-cell quantification of Bax activation and mathematical modelling suggest pore formation on minimal mitochondrial Bax accumulation. Cell Death Differ 17: 278-290
27. Faeder JR, Blinov ML, Hlavacek WS, (2009) Rule-based modeling of biochemical systems with BioNetGen. Methods Mol Biol 500: 113-167
28. Feret J, Danos V, Krivine J, Harmer R, Fontana W, (2009) Internal coarse-graining of molecular systems. Proc Natl Acad Sci USA 106: 6453-6458
29. Fire E, Gullá SV, Grant RA, Keating AE, (2010) Mcl-1-Bim complexes accommodate surprising point mutations via minor structural changes. Prot Sci 19: 507-519
30. Fricker N, Beaudouin J, Richter P, Eils R, Krammer PH, Lavrik IN, (2010) Model-based dissection of CD95 signaling dynamics reveals both a pro-and antiapoptotic role of c-FLIPL. J Cell Biol 190: 377-389
31. Gansner ER, North SC, (2000) An open graph visualization system and its applications to software engineering. Software: Practice and Experience 30: 1203-1233
32. Gaudet S, Spencer SL, Chen WW, Sorger PK, (2012) Exploring the contextual sensitivity of factors that determine cell-to-cell variability in receptor-mediated apoptosis. PLoS Comput Biol 8: e1002482
33. Gavathiotis E, Reyna DE, Davis ML, Bird GH, Walensky LD, (2010) BH3-triggered structural reorganization drives the activation of proapoptotic BAX. Mol Cell 40: 481-492
34. Gnad F, Estrada J, Gunawardena J, (2012) Proteus: a web-based, context-specific modelling tool for molecular networks. Bioinformatics 28: 1284-1286
35. Goldstein JC, Kluck RM, Green DR, (2000) A single cell analysis of apoptosis. Ordering the apoptotic phenotype. Ann N Y Acad Sci 926: 132-141
36. Hlavacek W, Faeder J, Blinov M, Posner R, (2006) Rules for modeling signal-transduction systems. Sci STKE 2006: re6
37. Hlavacek WS, (2009) How to deal with large models? Mol Syst Biol 5: 240
38. Hoops S, Sahle S, Gauges R, Lee C, Pahle J, Simus N, Singhal M, Xu L, Mendes P, Kummer U, (2006) COPASI-a COmplex PAthway SImulator. Bioinformatics 22: 3067-3074
39. Howells CC, Baumann WT, Samuels DC, Finkielstein CV, (2010) The Bcl-2-associated death promoter (BAD) lowers the threshold at which the Bcl-2-interacting domain death agonist (BID) triggers mitochondria disintegration. J Theor Biol 271: 114-123
40. Huber HJ, Dussmann H, Kilbride SM, Rehm M, Prehn JHM, (2011) Glucose metabolism determines resistance of cancer cells to bioenergetic crisis after cytochrome-c release. Mol Syst Biol 7: 470
41. Hucka M, Finney A, Sauro HM, Bolouri H, Doyle JC, Kitano H, Arkin AP, Bornstein BJ, Bray D, Cornish-Bowden A, Cuellar AA, Dronov S, Gilles ED, Ginkel M, Gor V, Goryanin II, Hedley WJ, Hodgman TC, Hofmeyr JH, Hunter PJ, et al, (2003) The systems biology markup language (SBML): a medium for representation and exchange of biochemical network models. Bioinformatics 19: 524-531
42. Hunter JD, (2007) Matplotlib: A 2D graphics environment. Comput Sci Eng 9: 90-95
43. Kim H, Rafiuddin-Shah M, Tu H-C, Jeffers JR, Zambetti GP, Hsieh JJ-D, Cheng EH-Y, (2006) Hierarchical regulation of mitochondrion-dependent apoptosis by BCL-2 subfamilies. Nat Cell Biol 8: 1348-1358
44. Kim H, Tu H-C, Ren D, Takeuchi O, Jeffers JR, Zambetti GP, Hsieh JJ-D, Cheng EH-Y, (2009) Stepwise activation of BAX and BAK by tBID, BIM, and PUMA initiates mitochondrial apoptosis. Mol Cell 36: 487-499
45. Kitano H, Funahashi A, Matsuoka Y, Oda K, (2005) Using process diagrams for the graphical representation of biological networks. Nat Biotechnol 23: 961-966
46. Kleiman LB, Maiwald T, Conzelmann H, Lauffenburger DA, Sorger PK, (2011) Rapid phospho-turnover by receptor tyrosine kinases impacts downstream signaling and drug binding. Mol Cell 43: 723-737
47. Krakauer DC, Collins JP, Erwin D, Flack JC, Fontana W, Laubichler MD, Prohaska SJ, West GB, Stadler PF, (2011) The challenges and scope of theoretical biology. J Theor Biol 276: 269-276
48. Kutumova EO, Kiselev IN, Sharipov RN, Lavrik IN, Kolpakov FA, (2012) A modular model of the apoptosis machinery. Adv Exp Med Biol 736: 235-245
49. Kuwana T, Bouchier-Hayes L, Chipuk JE, Bonzon C, Sullivan BA, Green DR, Newmeyer DD, (2005) BH3 domains of BH3-only proteins differentially regulate Bax-mediated mitochondrial membrane permeabilization both directly and indirectly. Mol Cell 17: 525-535
50. Le Novère N, Finney A, Hucka M, Bhalla US, Campagne F, Collado-Vides J, Crampin EJ, Halstead M, Klipp E, Mendes P, Nielsen P, Sauro H, Shapiro B, Snoep JL, Spence HD, Wanner BL, (2005) Minimum information requested in the annotation of biochemical models (MIRIAM). Nat Biotechnol 23: 1509-1515
51. Leber B, Lin J, Andrews DW, (2010) Still embedded together binding to membranes regulates Bcl-2 protein interactions. Oncogene 29: 5221-5230
52. Legewie S, Blüthgen N, Herzel H, (2006) Mathematical modeling identifies inhibitors of apoptosis as mediators of positive feedback and bistability. PLoS Comput Biol 2: e120
53. Letai A, Bassik MC, Walensky LD, Sorcinelli MD, Weiler S, Korsmeyer SJ, (2002) Distinct BH3 domains either sensitize or activate mitochondrial apoptosis, serving as prototype cancer therapeutics. Cancer Cell 2: 183-192
54. Llambi F, Moldoveanu T, Tait SWG, Bouchier-Hayes L, Temirov J, Mccormick LL, Dillon CP, Green DR, (2011) A unified model of mammalian BCL-2 protein family interactions at the mitochondria. Mol Cell 44: 517-531
55. Lovell JF, Billen LP, Bindner S, Shamas-Din A, Fradin C, Leber B, Andrews DW, (2008) Membrane binding by tBid initiates an ordered series of events culminating in membrane permeabilization by Bax. Cell 135: 1074-1084
56. Maiwald T, Timmer J, (2008) Dynamical modeling and multi-experiment fitting with PottersWheel. Bioinformatics 24: 2037-2043
57. Mallavarapu A, Thomson M, Ullian B, Gunawardena J, (2008) Programming with models: modularity and abstraction provide powerful capabilities for systems biology. J R Soc Interface 6: 257-270
58. Martinez-Caballero S, Dejean LM, Kinnally MS, Oh KJ, Mannella CA, Kinnally KW, (2009) Assembly of the mitochondrial apoptosis-induced channel, MAC. J Biol Chem 284: 12235-12245
59. Mathworks, Inc. (2012) MATLAB release 2012b The Mathworks, Inc.: Natick, Massachusetts
60. Mérino D, Giam M, Hughes PD, Siggs OM, Heger K, O'Reilly LA, Adams JM, Strasser A, Lee EF, Fairlie WD, Bouillet P, (2009) The role of BH3-only protein Bim extends beyond inhibiting Bcl-2-like prosurvival proteins. J Cell Biol 186: 355-362
61. Mirschel S, Steinmetz K, Rempel M, Ginkel M, Gilles ED, (2009) PROMOT: modular modeling for systems biology. Bioinformatics 25: 687-689
62. Moraru II, Schaff JC, Slepchenko BM, Loew LM, (2002) The virtual cell: an integrated modeling environment for experimental and computational cell biology. Ann N Y Acad Sci 971: 595-596
63. Muzzey D, Gõmez-Uribe CA, Mettetal JT, van Oudenaarden A, (2009) A systems-level analysis of perfect adaptation in yeast osmoregulation. Cell 138: 160-171
64. Neumann L, Pforr C, Beaudouin J, Pappa A, Fricker N, Krammer PH, Lavrik IN, Eils R, (2010) Dynamics within the CD95 death-inducing signaling complex decide life and death of cells. Mol Syst Biol 6: 352
65. O'Connor CL, Anguissola S, Huber HJ, Dussmann H, Prehn JHM, Rehm M, (2008) Intracellular signaling dynamics during apoptosis execution in the presence or absence of X-linked-inhibitor-of-apoptosis-protein. Biochim Biophys Acta 1783: 1903-1913
66. Oliphant TE, (2007) Python for Scientific Computing. Comput. Sci. Eng 9: 10-20
67. Pedersen M, Plotkin G, (2008) A language for biochemical systems. Lecture Notes in Computer Science 5307: 63-82
68. Rehm M, Dussmann H, Janicke RU, Tavare JM, Kogel D, Prehn JHM, (2002) Single-cell fluorescence resonance energy transfer analysis demonstrates that caspase activation during apoptosis is a rapid process. Role of caspase-3. J Biol Chem 277: 24506-24514
69. Rehm M, Huber HJ, Dussmann H, Prehn JHM, (2006) Systems analysis of effector caspase activation and its control by X-linked inhibitor of apoptosis protein. EMBO J 25: 4338-4349
70. Rehm M, Huber HJ, Hellwig CT, Anguissola S, Dussmann H, Prehn JHM, (2009) Dynamics of outer mitochondrial membrane permeabilization during apoptosis. Cell Death Differ 16: 613-623
71. Ren D, Tu H-C, Kim H, Wang GX, Bean GR, Takeuchi O, Jeffers JR, Zambetti GP, Hsieh JJ-D, Cheng EH-Y, (2010) BID, BIM, and PUMA are essential for activation of the BAX-and BAK-dependent cell death program. Science 330: 1390-1393
72. Schleich K, Warnken U, Fricker N, Öztürk S, Richter P, Kammerer K, Schnölzer M, Krammer PH, Lavrik IN, (2012) Stoichiometry of the CD95 death-inducing signaling complex: experimental and modeling evidence for a death effector domain chain model. Mol Cell 47: 306-319
73. Smith AM, Xu W, Sun Y, Faeder JR, Marai GE, (2012) RuleBender: integrated modeling, simulation and visualization for rule-based intracellular biochemistry. BMC Bioinformatics 13 (Suppl 8): S3
74. Sneddon MW, Faeder JR, Emonet T, (2010) Efficient modeling, simulation and coarse-graining of biological complexity with NFsim. Nat Meth 8: 177-183
75. Spencer SL, Gaudet S, Albeck JG, Burke JM, Sorger PK, (2009) Non-genetic origins of cell-to-cell variability in TRAIL-induced apoptosis. Nature 459: 428-432
76. SymPy Development Team. (2012) SymPy: Python library for symbolic mathematics. Available at http://www.sympy.org
77. Tan C, Dlugosz PJ, Peng J, Zhang Z, Lapolla SM, Plafker SM, Andrews DW, Lin J, (2006) Auto-activation of the apoptosis protein Bax increases mitochondrial membrane permeability and is inhibited by Bcl-2. J Biol Chem 281: 14764-14775
78. Tiger C-F, Krause F, Cedersund G, Palmér R, Klipp E, Hohmann S, Kitano H, Krantz M, (2012) A framework for mapping, visualisation and automatic model creation of signal-transduction networks. Mol Syst Biol 8: 578
79. Waltemath D, Adams R, Bergmann FT, Hucka M, Kolpakov F, Miller AK, Moraru II, Nickerson D, Sahle S, Snoep JL, Le Novère N, (2011) Reproducible computational biology experiments with SED-ML-the simulation experiment description markup language. BMC Syst Biol 5: 198
80. Willis SN, Chen L, Dewson G, Wei A, Naik E, Fletcher JI, Adams JM, Huang DCS, (2005) Proapoptotic Bak is sequestered by Mcl-1 and Bcl-xL, but not Bcl-2, until displaced by BH3-only proteins. Genes Dev 19: 1294-1305
81. Willis SN, Fletcher JI, Kaufmann T, van Delft MF, Chen L, Czabotar PE, Ierino H, Lee EF, Fairlie WD, Bouillet P, Strasser A, Kluck RM, Adams JM, Huang DCS, (2007) Apoptosis initiated when BH3 ligands engage multiple Bcl-2 homologs, not Bax or Bak. Science 315: 856-859
82. Wolfram Research, Inc. (2010) Mathematica 8th edn. Wolfram Research, Inc.: Champaign, Illinois
83. Xu T-R, Vyshemirsky V, Gormand A, Kriegsheim von A, Girolami M, Baillie GS, Ketley D, Dunlop AJ, Milligan G, Houslay MD, Kolch W, (2010) Inferring signaling pathway topologies from multiple perturbation measurements of specific biochemical species. Sci Signal 3: ra20
84. Youle RJ, Strasser A, (2008) The BCL-2 protein family: opposing activities that mediate cell death. Nat Rev Mol Cell Biol 9: 47-59