390 likes | 508 Views
CS 572 Fall 2013 Lecture 6 Tool for setting up kinetic simulations. Example 2: second order and enzymatic reactions. t02.py. Cell Designer is tool for creating and sharing models. Early version. When used in conjunction with Systems Biology Workbench, Cell designer will
E N D
CS 572 Fall 2013 Lecture 6 Tool for setting up kinetic simulations
Cell Designer is tool for creating and sharing models Early version
When used in conjunction with Systems Biology Workbench, Cell designer will (1) prepare MATLAB files to solve ODEs (2) save your model as an XML file in systems biology markup language (Hucka et al, Bioinformatics 19 (2004) 524 SourceForge.net
http://celldesigner.org/help/Install/ With Cell Designer, you draw the network
Boris network drawn in Cell Designer Cell Designer sets up the tables of reactants, reactions, and parameters
id defined in order that species is added to graph To get a particular order for spices id (e.g. Boris table I), put all species on graph before reactions
Directions of arrows correspond to rate laws in Boris Table II Reactions added to the model in numerical order
reactions listed in the same order as Boris Table II Species id used to define reactants and products. User decides names of rate parameters User has option to define each rate law (i.e. “math”)
Put in the value of rate constants. These values are different from Boris Table II to make number per cell the dynamic variable
Choose MATLAB from options offered by Systems Biology Workbench
Supplement MATLAB m-files as needed to make graphs in Figure 4A of Boris paper.
MATLAB main program assigns initial values, sets time span, calls solver, and plots results
More recent versions of Cell Designer have more features http://celldesigner.org/help/Install/
MAPK cascade (with feedback) in SBGN Denotes enzyme MKK = MAPK kinase; example MEK MKKK= MAPK kinase kinase: example RAF ERK is an example of MAPK
Model of MAPK cascade with feedback can oscillate
CellDesigner 4.1 supports “stochastic” simulations Concentrations used in ODE method are average quantities. When number of molecules is small, fluctuations become important. Stochastic method: discrete time steps reaction rates probabilities choose next reaction by random numbers update number of molecules after each reaction
Questions on lecture 6 What is Systems Biology Markup Language (SBML)? What is Systems Biology Graphis Notation (SBGN)? Contrast Cell Designer and Reactor with respect the following: model creation, rate law specification parameter input, ODE generation and solution model sharing.
More questions on lecture 6 Contrast ODE-based and stochastic simulations of signaling kinetics. Under what condition is stochastic simulation needed? Give examples of enzymes MAPK, MKK, and MKKK. In SBNG, how do one show that E is an enzyme in reaction A B (draw) What type of enzymes reverse phosphorylation cascades?
Cell Designer: You draw the network Dr. Miller: Here's the link for CellDesigner: http://celldesigner.org/ It contains everything you need about it. If you need any help about installation and how to use it, please let me know. Best wishes, Shuangshuang
Down load the Systems Biology Workbench before installing Cell Designer Can be accomplished using folder “windows install for students” on class web page • Extract sbw-1.0.5-src • Run sbw-1.0.5-win-full (you will be ask to install java VM) • Run CDInstall-2-2-win (choose option install with SBW 2.0)