User:Geoff S. Baldwin/Sandbox
From OpenWetWare
Contents 
My Sandbox
Section 1
SubSection
 bullet point
 bullet point
 bullet point
 bullet point
Section 2
 topic 1
 topic 2
 topic 3
Practical 1
Objectives:
 To learn how to use a computational modelling tool for biochemical reaction simulations.
 To build biochemical networks
 To simulate the time evolution of the reactions
 To explore the properties of simple biochemical reactions.
 A > B > C model
 SynthesisDegradation model
 MichaelisMenten model
Deliverables
 A brief report is expected ( see structure). Report due Monday 18th February 2008.
 When you find in the text (illustration needed), it means that you will have to provide an image export of your simulation results in your report.
 Instructors: Vincent Rouilly, Geoff Baldwin.
Part I: Introduction to Computer Modelling
 "All models are wrong, but some of them are useful", George Box.
Part II: Getting to know CellDesigner
 Read through the tutorial example, and get familiar with CellDesigner features. Official CellDesigner Tutorial
 Open a sample file: File > Open > Samples/...
 Select items, move them around, delete, undo...
Part III: Building Your First Model: A > B > C
In this section, you will build your first model from scratch with CellDesigner, and you will learn to run a simulation.
The model explored describe a system where a compound 'A' is transformed into a compound 'B', which is consequently transformed into a compound 'C'.
To start, launch the CellDesigner Application: Double Click on the Icon found on your Desktop. Then follow the instructions below to build the model.
Model  CellDesigner Instructions 

 
Following the Law of Mass action, the dynamic of the system is described as: 

Simulate the dynamical behaviour 

 Questions:(see report structure)
 How do A, B and C, change with time using these default parameters?
 Now swap the values of k1 and k2 (k1=10 and k2=1)under the parameters tab
 How does this alter the formation of C?
 How does B change?
 Explain these results
 If you had real life data showing the accumulation of C for an ABC reaction you could fit the data using this model and two rate constants would be returned. Could you assign these rate constants to k1 or k2 (yes or no)?
 What additional data would you need to assign k1 and k2?
 Additional Resources:
Part IV: SynthesisDegradation Model
In this section, we investigate a very common motif in biochemistry. It models the continuous and constant synthesis of a compound, and its natural degradation.
From a Mathematical point of view, the model is described as a firstorder linear ordinary differential equation.
Model  CellDesigner Instructions 


Build the topology of the reaction network

From the law of mass action, we can write:

Define the kinetics driving the reaction network

Simulate the dynamical behaviour 

 Questions:(see report structure)
 Run a simulation over t=1000s, nb points=1000. Comment on the time evolution of 'A'. (illustration needed).
 Using the dynamical system definition, what is the steady state level of 'A' with regards to the parameters k1 and k2 ? (Steady state means that
 Using the 'Parameter Scan' feature, illustrate the influence of both parameters (k_1 and k_2), on the steady state level of 'A' (illustration needed).
 Bonus: Give the analytical solution of the ODE system.
 Now, consider that k_1=0, and [A]_{t = 0} = A_{0} > 0. Keep k_2=0.01. Illustrate the concept of halflife for the compound 'A'.
 Bonus: Derive the analytical expression of the halflife of 'A', with regards to k_1 and k_2.
 Additional Resources:
Part V: MichaelisMenten Model
An enzyme converts a substrate into a product, this is usually an irreversible reaction and is treated as such in the MichaelisMenten model. An enzyme reaction constitutes a dynamic process and can be studied as such. One may look at the time courses of the reactants, or look at the steadystates and their stability properties.
This part of the tutorial deals with wellknown MichaelisMenten formula.
Here, we will focus on comparing the MichaelisMenten approximation to the full enzymatic reaction network.
Model  CellDesigner Instructions 



Following law of mass action, we can write:

Questions:(see report structure)
 From the ODE system description, create all the necessary kinetics reactions in the network provided. We will be considering K_{1} = 10^{5}M^{ − 1}s^{ − 1},K_{2} = 1000s^{ − 1},K_{3} = 10^{ − 1},K_{4} = 2M^{ − 1}s^{ − 1},[E]_{t = 0} = 0.01M,[S]_{t = 0} = 0.01M,[P]_{t = 0} = 0
 Open the Simulation Panel, set Time=150, NbPoints=1000.
 Run a simulation, and comment on the different phases during the product formation. Pay special attention to the formation and decay of the [ES] complex. Note that this is a full simulation of the reaction scheme and so does not rely on any assumptions.
Model  CellDesigner Instructions 


 We want now to investigate the MichaelisMenten approximation. Show that under the assumption that the complex [ES] is at steadystate (), we can write: . (Note that [E]_{t = 0} = [E]_{t} + [ES]_{t}). Also, make sure that the concentration of the substrate is at least 10 fold greater than the concentration of the enzyme.
 Express (Km and Vmax) with regards to K_1, K_2, K_3 and [E]_{0}; see the links below on how the MichaelisMenten equation is derived if you are not sure.
 Now create a new reaction in CellDesigner(as shown above)with an Enzyme that acts on the reaction. Define the maths for this reaction based on the above form of the Michaelis Menten equation. Make sure that both models are equivalent with regards to their parameters.
 Run simulations, and comment on the differences observed between the full model, and the MichaelisMenten approximation.