Imperial College/Courses/Spring2008/Synthetic Biology/Computer Modelling Practicals/Practical 2
<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>
Practical 2
Objectives:
- To explore computationally some simple genetic motifs:
- Constitutive gene expression.
- Activated and Repressed gene expression.
- Positive and Negative feedback gene expression.
Deliverables
- A report is expected.
- 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: Constitutive Gene Expression
In this section, we explore a model to represent constitutive gene expression. The model is based on a simple interpretation of the central dogma: gene --> mRNA --> Protein. Both the mRNA and the protein are also naturally degraded.
Model | CellDesigner Instructions |
---|---|
|
Define the topology of the reaction network:
|
Following the law of mass action, we can write: |
|
Questions:
- From the ODE system given, write down the steady-state expression of [mRNA] and [Protein], with regards to k_1, k_2, d_1, d_2.
- Knowing that average mRNA per gene is 2.5 in E.Coli, what is the average transcription rate ?
- Knowing that the average number of proteins per gene is 1000 in E.Coli, what is the average translation rate ?
- Define all the necessary kinetics laws on the model, and create all the appropriate parameters.
- Run a simulation, and comment on the simulation ouputs.
Next, we want to explore the quasi-steady-state assumption on the mRNA expression. From the previous simulations, you must have noticed that the concentration of mRNA is reaching steady-state very quickly, compared to the Protein concentration. In that case, we can approximate the previous model by considering that [math]\displaystyle{ \frac{d[mRNA]}{dt}=0 }[/math].
Model | CellDesigner Instructions |
---|---|
|
|
Following the law of mass action, we can write: |
|
Questions:
- Taking into account the quasi-steady-state assumption on the [mRNA], work-out the value of 'SynthesisRate', and 'DegradationRate' with regards to K_1, K_2, d_1, d_2.
- Simulate the full model and the quasi-steady state approximation. Comment on how good this approximation seems to be.
- Additional Resources:
Part II: Activated and Repressed Gene Expression
Very few genes are known to have a purely constitutive expression, most of genes have their expression controlled by outside signals. We will particularly be interested in the study of the transcription factors, these proteins can either have an activation effect on the gene (positive control), or a repression effect (negative control). In this part of the tutorial, we investigate the most common model used to describe this type of interactions.
Let's first consider the case of a transcription factor acting as a repressor. A repressor will bind to the DNA so that it prevents the initiation of transcription. Typically, we expect the transcription rate to decrease as the concentration of repressor increases. A very useful family of functions to describe this effect is the Hill function: [math]\displaystyle{ ... }[/math] This function has 3 parameters: Vmax, n, Km
Model | CellDesigner Instructions |
---|---|
|
Define the topology of the reaction network:
|
Following the law of mass action, we can write: |
Then, let's consider the case of a transcription factor acting as an activator. A activator will bind to the DNA so that it promotes the initiation of transcription. Typically, we expect the transcription rate to increase as the concentration of repressor increases. Once again, q very useful family of functions to describe this effect is the Hill function. It is slightly different from the previous one:
[math]\displaystyle{
...
}[/math]
This function has 3 parameters: Vmax, n, Km
Model | CellDesigner Instructions |
---|---|
|
Define the topology of the reaction network:
|
Following the law of mass action, we can write: |
Questions:
- Draw the transfer function between the Activator concentration and the steady-state protein concentration. And illustrate the effect of the parameters 'Km' and 'n'.
- Imagine an application where this behaviour might be useful.
- Draw the transfer function between the Repressor concentration and the steady-state protein concentration. And illustrate the effect of the parameters 'Km' and 'n'.
- Imagine an application where this behaviour might be useful.
Additional Resources:
Part III: Positive and Negative Feedbacks
In this section, ....
Model | CellDesigner Instructions |
---|---|
|
|
Following the law of mass action, we can write: |
|
Questions:
- ...