Imperial College/Courses/Spring2008/Synthetic Biology/Computer Modelling Practicals/Practical 1
<html> <body> <!-- Start of StatCounter Code --> <script type="text/javascript"> var sc_project=3315875; var sc_invisible=0; var sc_partition=36; var sc_security="779debd0"; </script> <script type="text/javascript" src="http://www.statcounter.com/counter/counter_xhtml.js"></script><noscript><div class="statcounter"><a class="statcounter" href="http://www.statcounter.com/"><img class="statcounter" src="http://c37.statcounter.com/3315875/0/779debd0/0/" alt="web metrics" /></a></div></noscript> <!-- End of StatCounter Code --> </body> </html>
...under development...
Practical 1
Objectives:
- Learn how to use a computational modelling tool for biochemical reaction simulations.
- Building biochemical networks
- Simulating the time evolution of the reactions
- Explore the properties of simple biochemical reactions.
- A --> B --> C model
- Synthesis-Degradation model
- Michaelis-Menten model
Deliverables
- A report is expected by ... (Word or PDF format, sent to XXX@XXX)
- 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.
Part I: Introduction to Computer Modelling
- Presentation Slides: "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:
- Describe the time evolution of A, B and C, taking into account the default parameters.
- Using the 'Parameter Scan' function, investigate how parameters 'k1' and 'k2' influence the production of 'C'.
- Find the set of parameters (k1, k2), within a 10% range of their initial value, so that B is maximal at some point in time.
- Additional Resources:
Part IV: Synthesis-Degradation Model
In this section, we are going to investigate a very common motif in biochemistry. It models the synthesis of a compound and its natural degradation. From a Mathematical point of view, the model is described as a first-order linear ordinary differential equation.
Model | CellDesigner Instructions |
---|---|
Build the topology of the reaction network
| |
Dynamical system, Law of Mass action:
|
Define the kinetic driving the reaction network
|
Simulate the dynamical behaviour |
|
- Questions:
- Run a simulation over t=1000s, 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 [math]\displaystyle{ \frac{d[A]}{dt}=0 }[/math]
- Using the 'Parameter Scan' feature, illustrate the influence of both parameters, on the steady state level of 'A' (illustration needed).
- now consider that k_1=0, and [math]\displaystyle{ [A]_{t=0}=A_{0} \gt 0 }[/math]. Keep k_2=0.01. Illustrate the concept of half-life for the compound 'A'.
- Additional Resources:
Part V: Michaelis-Menten Model
Model | CellDesigner Instructions |
---|---|
| |
Dynamical system, Law of Mass action: [math]\displaystyle{ a }[/math]
|
|
Simulate the dynamical behaviour |
|
- Questions:
- What is the influence of each parameter on how quickly the Product is formed ?
- Consider the Michaelis-Menten assumption, and derive the simplified system.
- Update the reaction network at the bottom so that the two reactions are equivalent.
- Simulate the full set of reactions and comment on the difference.