IGEM:IMPERIAL/2008/New/Genetic Circuit

From OpenWetWare

(Difference between revisions)
Jump to: navigation, search
(Genetic Circuit)
(Genetic Circuit)
Line 19: Line 19:
[[Image:Phase 1.PNG|thumb|300px]]
[[Image:Phase 1.PNG|thumb|300px]]
A simple [[Imperial_College/Courses/Spring2008/Synthetic_Biology/Computer_Modelling_Practicals/Practical_2 | synthesis-degradation model]] is assumed for the modelling of the expression of a protein under the control of a constitutive promoter, with the same model assumed for all four [[IGEM:IMPERIAL/2008/Prototype/Wetlab/test_constructs | promoter-RBS constructs]]. The synthesis-degradation model assumes a steady state level of mRNA.
A simple [[Imperial_College/Courses/Spring2008/Synthetic_Biology/Computer_Modelling_Practicals/Practical_2 | synthesis-degradation model]] is assumed for the modelling of the expression of a protein under the control of a constitutive promoter, with the same model assumed for all four [[IGEM:IMPERIAL/2008/Prototype/Wetlab/test_constructs | promoter-RBS constructs]]. The synthesis-degradation model assumes a steady state level of mRNA.
 +
<math>\frac{d[protein]}{dt} = k_{1} - d_{1}[protein]</math><p>
<math>\frac{d[protein]}{dt} = k_{1} - d_{1}[protein]</math><p>
 +
In this case, <math>[protein]</math> represents the concentration of GFP, <math>k_{1}</math> represents the rate of sythesis and <math>d_{1}</math> represents the degradation rate.
In this case, <math>[protein]</math> represents the concentration of GFP, <math>k_{1}</math> represents the rate of sythesis and <math>d_{1}</math> represents the degradation rate.
We can easily simulate this synthesis-degradation model using matlab:<br>
We can easily simulate this synthesis-degradation model using matlab:<br>

Revision as of 06:58, 10 September 2008




Contents

Genetic Circuit

Authors: Erika

Editors: Erika

Erika's to do list: Upload F2620 citation. Change internal links to redirect to new wiki

Why model the genetic circuit?

An accurate mathematical description of genetic circult behaviour is one of the foundations of synthetic biology. Such descriptions are an integral component of part submission to the registry, as exemplified by the canonical characterised part F2620. <citation needed>. The ability to capture part behaviour as a mathematical relationship between input and output is useful for future re-use of the part modification of integration into novel genetic circuits.

Modelling Constitutive Gene Expression

A simple synthesis-degradation model is assumed for the modelling of the expression of a protein under the control of a constitutive promoter, with the same model assumed for all four promoter-RBS constructs. The synthesis-degradation model assumes a steady state level of mRNA.

\frac{d[protein]}{dt} = k_{1} - d_{1}[protein]

In this case, [protein] represents the concentration of GFP, k1 represents the rate of sythesis and d1 represents the degradation rate. We can easily simulate this synthesis-degradation model using matlab:
ODE
Simulation File

We can also solve this ODE analytically.
At \frac{d[protein]}{dt}=0, [protein]=\frac{k_1}{d_1} and you can see this relationship in the parameter scan graphs.
From the wetlab experiments it is likely that we will obtain steady-state data for each of the four promoter-RBS constructs. If we assume the same rate of degradation of GFP in each case, we can have some measure of the relative rate of transcription through each promoter which will help us with the selection of the most appropriate promoter to use for Phase 2. In order to obtain an absolute measure of transcription (as opposed to a relative measure of transcriptional strength) we require constitutive expression in terms of molecules per cell (as opposed to fluorescene in arbitrary units).
Note from the parameter scan graphs:

  • In the case where k1 = 0, no GFP is sythesised.
  • In the case where d1 = 0, the concentration of protein does not reach a steady state.

Modelling Inducible Gene Expression

The repressor is constitutively expressed. Hence we can assume the constitutive expression model from the previous characterisation step.

\frac{d[R]}{dt} = k_{1} - d_{1}[R]

When the inducer is added it binds reversibly to the repressor.

 R + I \Leftrightarrow RI

Free repressor only binds to the promoter, we think this will show co-operative binding as there are two repressor binding sites on the promoter sequence. Then transcription will be a function of free repressor concentration.

 Transcription = \frac{\beta.{[R]}^n}{{K_m}^n+[R]^n}

And overall protein expression can be described as

\frac{d[protein]}{dt} = Transcription - d_2[protein]



    [[IGEM:IMPERIAL/2008/New/{{{1}}}|< Previous]]
    [[IGEM:IMPERIAL/2008/New/{{{2}}}|Next >]]

</div></td></tr></table>

Personal tools