350 likes | 454 Views
Cellerator is a system that simulates complex biochemical reaction networks using ODEs. It translates biochemical pictures into ODEs, providing a quantitative understanding. With a GUI and modeling meta-language, Cellerator conducts simulations of cellular activities like cell division. It enables the translation of biochemistry into canonical forms and offers predictions for various biological processes. With applications in diverse fields like myogenesis, plant growth, and leukemia, Cellerator facilitates collaborative research in systems biology and biomedicine.
E N D
Cellerator: A System for Simulating Biochemical Reaction Networks Bruce E Shapiro Jet Propulsion Laboratory California Institute of Technology bshapiro@jpl.nasa.gov
Part of a Biochemical Network From: Kohn (1999) Molecular interaction map of the mammalian cell cycle control and DNA repair systems. Mol Biol Cell 10:2703-2734
Biochemical Networks Are... • Complex • Mutually interacting • Large • Number of reactions grows exponentially with number of states • Best understood pictorially • Best described quantitatively by a large system of differential equations (ODEs) Need to translate pictures to ODEs
Online network databases exist ... http://www.genome.ad.jp/kegg/
... but mathematical simulations of these networks are hopelessly naive...
A B C Input Canonical Form Biochemical Notation Output Canonical Form System of ODEs Activity (e.g., Cell Division) Solver Concentrations vs. Time
Caltech ERATO* Simulator Architecture A GUI and Modeling meta-language C B Text Transfer Protocol XML based protocol Application Application Application Application *Exploratory Research for Advanced Technology (Japan Science & Technology Corporation) http://www.systems-biology.org Application
A simpler network for cell division C=Cyclin: enzyme that gets things going M=MPF promoting factor. M>Threshold induces cell division X=Cyclin Protease: enzyme that breaks down C Goldbeter, A (1991) A minimal cascade model for the mitotic oscillator involving cyclin and cdc2 kinase. PNAS 88:9107-9111
Cellerator canonical form for input Reactions are input with a biochemical based notation Prints out ODES STN = {{reaction, rate-constants}, {reaction, rate-constants},…}; interpret[STN]; Simulation = predictTimeCourse[STN, options]; Returns tables of values as a function of time, with optional plots
Simple Cooperative Conversion Creation, Degradation Enzymatic Reversible Enzymatic Transcription (Gene RNA) Post-transcriptional Processing Translation (RNA Protein) Diffusion and more ... The Basis of Cellerator: Chemical Reactions
Translation of Biochemical Formula to ODE rate constant Concentrations • Law of Mass Action • Two-way Reaction • Complex reactions built from simple reactions is described by Similar ODE’s can be written for B and C is described by
Enzyme Kinetic (Catalytic) Reaction • Enzyme Ecatalyzes the production of product P from substrate (source) S • Write more compactly as 3 Reactions written two different ways Explicit Cellerator syntax for this set of reactions Hidden Rate constants
Two-way catalytic reaction • A second enzyme F catalyzes the reverse reaction • Total of Six Elementary Reactions • Write more compactly as Explicit Rate constants Hidden Cellerator syntax for this set of reactions
Canonical Forms for Translation: Chemical reactions • Input Canonical Form for Chemical Reaction • Output Canonical Form: Terms in an ODE
INPUT OUTPUT MAP Kinase Cascade
Object Oriented Implementation:“Domains” and “Fields” • Domain: object • Field: function that maps domains to R • Field of Domains: maps domain elements to domains • Example • graphDomain: represents tissue • node Domains: cells • neighbors[g,n] returns a list of nodeDomains that are neighbors of node n n in graph g
Myogenesis: Collaboration with Laboratory Dr. Barbara Wold (Chris Hart), Caltech
Plant Growth: Collaboration with Laboratory Dr. Elliot Meyerowitz, Caltech
Secondary Leukemia: Collaboration with City of Hope National Medical Center (NASA/BSRP) Focus: Pathogenesis of myelodysplasia & acute myeloid leukemia following high-dose chemo/radiotherapy and autologous peripheral blood stem cell transplantation for treatment of Hodgkin’s disease and non-Hodgkin’s lymphoma
JPL Collaborations using Cellerator • Effects of microgravity during space flight on bone and muscle development (Caltech, JSC, and UCI) • Development of childhood leukemias (Caltech, Children’s Hospital of LA, and UC, Irvine) • Description of “core” signal transduction units (Johns Hopikins) • Improving algorithms for micro-array data analysis (Caltech, Harvey Mudd) • Systems Biology Workbench (Caltech, JST/Erato)
Acknowledgements • Eric Mjolsness* - UC, Irivine • Andre Levchenko* - Johns Hopkins University • Barbara Wold - Caltech • Elliot Meyerowitz - Caltech * Original Developers