Revision as of 08:28, 12 September 2008

<html> <style type="text/css"> .firstHeading {display: none;} </style> </html> <html> <style type="text/css">

   table.calendar          { margin:0; padding:2px; }

table.calendar td { margin:0; padding:1px; vertical-align:top; } table.month .heading td { padding:1px; background-color:#FFFFFF; text-align:center; font-size:120%; font-weight:bold; } table.month .dow td { text-align:center; font-size:110%; } table.month td.today { background-color:#3366FF } table.month td {

   border:2px;
   margin:0;
   padding:0pt 1.5pt;
   font-size:8pt;
   text-align:right;
   background-color:#FFFFFF;
   }

bodyContent table.month a { background:none; padding:0 }

.day-active { font-weight:bold; } .day-empty { color:black; } </style> </html>

<html><a href=http://openwetware.org/wiki/IGEM:IMPERIAL/2008/Prototype><img width=50px src=http://openwetware.org/images/f/f2/Imperial_2008_Logo.png></img</a></html>	Home	The Project	B.subtilis Chassis	Wet Lab	Dry Lab	Notebook

<html> <style type="text/css"> div.Section { font:11pt/16pt Calibri, Verdana, Arial, Geneva, sans-serif; }

/* Text (paragraphs) */ div.Section p { font:11pt/16pt Calibri, Verdana, Arial, Geneva, sans-serif; text-align:justify; margin-top:0px; margin-left:30px; margin-right:30px; }

/* Headings */ div.Section h1 { font:22pt Calibri, Verdana, Arial, Geneva, sans-serif; text-align:left; color:#3366FF; }

/* Subheadings */ div.Section h2 { font:18pt Calibri, Verdana, Arial, Geneva, sans-serif; color:#3366FF; margin-left:5px; }

/* Subsubheadings */ div.Section h3 { font:16pt Calibri, Verdana, Arial, sans-serif; font-weight:bold; color:#3366FF; margin-left:10px; }

/* Subsubsubheadings */ div.Section h4 { font:12pt Calibri, Verdana, Arial, sans-serif; color:#3366FF; margin-left:15px; }

/* Subsubsubsubheadings */ div.Section h5 { font:12pt Calibri, Verdana, Arial, sans-serif; color:#3366FF; margin-left:20px; }

/* References */ div.Section h6 { font:12pt Calibri, Verdana, Arial, sans-serif; font-weight:bold; font-style:italic; color:#3366FF; margin-left:25px; }

/* Hyperlinks */ div.Section a {

}

div.Section a:hover {

}

/* Tables */ div.Section td { font:11pt/16pt Calibri, Verdana, Arial, Geneva, sans-serif; text-align:justify; vertical-align:top; padding:2px 4px 2px 4px; }

/* Lists */ div.Section li { font:11pt/16pt Calibri, Verdana, Arial, Geneva, sans-serif; text-align:left; margin-top:0px; margin-left:30px; margin-right:0px; }

/* TOC stuff */ table.toc { margin-left:10px; }

table.toc li { font: 11pt/16pt Calibri, Verdana, Arial, Geneva, sans-serif; text-align: justify; margin-top: 0px; margin-left:2px; margin-right:2px; }

/* [edit] links */ span.editsection { color:#BBBBBB; font-size:10pt; font-weight:normal; font-style:normal; vertical-align:bottom; } span.editsection a { color:#BBBBBB; font-size:10pt; font-weight:normal; font-style:normal; vertical-align:bottom; } span.editsection a:hover { color:#3366FF; font-size:10pt; font-weight:normal; font-style:normal; vertical-align:bottom; }

sddm {

margin: 0; padding: 0; z-index: 30 }

sddm li {

margin: 0; padding: 0; list-style: none; float: center; font: bold 12pt Calibri, Verdana, Arial, Geneva, sans-serif; border: 0px }

sddm li a {

display: block; margin: 0px 0px 0px 0px; padding: 0 0 12px 0; background: #33bbff; color: #FFFFFF; text-align: center; text-decoration: none; }

sddm li a:hover {

border: 0px }

sddm div {

position: absolute; visibility: hidden; margin: 0; padding: 0; background: #33bbff; border: 1px solid #33bbff } #sddm div a { position: relative; display: block; margin: 0; padding: 5px 10px; width: auto; white-space: nowrap; text-align: left; text-decoration: none; background: #FFFFFF; color: #2875DE; font: 11pt Calibri, Verdana, Arial, Geneva, sans-serif } #sddm div a:hover { background: #33bbff; color: #FFFFFF } </style></html>

ODE models of test constructs

The genetic circuitry involved in each test construct can be individually modelled.

These individual models can be compared with experimental results at each stage, and thus the individual models can be refined.

Individual models can be combined to generate a model of the complete circuit.

This model can be compared to the experimental results to evaluate how well it represents the behaviour of the system.

[1]

Characterisation of Constitutive promoter and RBS

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. [math]\displaystyle{ \frac{d[A]}{dt} = k_{1} - d_{1}[A] }[/math] Where [math]\displaystyle{ [A] }[/math] represents the concentration of protein (in this case GFP), [math]\displaystyle{ k_{1} }[/math] represents the rate of sythesis of the protein and [math]\displaystyle{ d_{1} }[/math] represents the degradation rate. We can easily simulate this synthesis-degradation model using matlab:
ODE
Simulation File

$Parameter scan on [math]\displaystyle{ k_1 }[/math]$

Parameter scan on [math]\displaystyle{ k_1 }[/math]
$Parameter scan on [math]\displaystyle{ d_1 }[/math]$

Parameter scan on [math]\displaystyle{ d_1 }[/math]

We can also solve this ODE analytically.
At [math]\displaystyle{ \frac{d[A]}{dt}=0 }[/math], [math]\displaystyle{ [A]=\frac{k_1}{d_1} }[/math] 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 [math]\displaystyle{ k_1 = 0 }[/math], no GFP is sythesised.
In the case where [math]\displaystyle{ d_1 = 0 }[/math], the concentration of protein does not reach a steady state.

Characterisation of inducible promoter and RBS

Characterising the inducible promoter is a bit less straightforward. There are four inducible promoter + RBS combinations.

pHyper-spank-gsiB,
pHyper-spank-SpoVG,
pxyl-gsiB,
pxyl-SpoVG

For each of the above four inducible promoter + RBS combinations, we want to experimentally determine the maximum protein expression level. (In these test constructs the reporter protein is GFP). By varying the concentration of the inducer we want to characterise the transfer function between the inducer concentration and the steady-state GFP expression level. In addition, we want to characterise the time taken for induction, and how this is affected by concentrations of inducer.

The image on the right shows the genetic construct for the inducible promoter. A transcriptional regulator is constitutively expressed: In the case of the promoter pHyper-spank, the transcriptional regulator is LacI, and in the case of the promoter pxyl, the transcriptional regulator is XylR. Constitutive expression levels of these transcriptional repressors can be taken from the above characterisation using GFP as a reporter.

Factors affecting the steady state protein level downstream of the inducible promoter:

Basal ("leaky") expression through the inducible promoter in the absence of inducer.
Expression level of repressor.
Concentration of inducer.
Functional relationship between concentration of inducer, concentration of repressor and rate of transcription.
Degradation rate of protein.

We describe repressor-inducer interactions by using LacI and IPTG in an example:

LacI reversibly binds IPTG.

[math]\displaystyle{ [LacI] + [IPTG] \longleftrightarrow [IPTG-LacI] }[/math]

With a [math]\displaystyle{ k_{on} }[/math] of [math]\displaystyle{ k_2 }[/math] and a [math]\displaystyle{ k_{off} }[/math] of [math]\displaystyle{ k_3 }[/math]

LacI is constitutively expressed and enters into an equilibrium with the repressor, IPTG:

[math]\displaystyle{ \frac{d[LacI]}{dt}=k_1-d_1[LacI]-k_2[IPTG][LacI]+k_3[IPTG-LacI] }[/math]

The concentration of free IPTG is dependent on the concentration of free LacI:

[math]\displaystyle{ \frac{d[IPTG]}{dt}=-k_2[IPTG][LacI]+k_3[IPTG-LacI] }[/math]

And the concentration of [IPTG-LacI] must also be accounted for in the ODE model.

[math]\displaystyle{ \frac{d[IPTG-LacI]}{dt}=k_2[LacI][IPTG]-k_3[IPTG-LacI] }[/math]

The difficulty is that the interaction between the inducer, the repressor protein and the operator DNA is not simple. It depends on the relative concentrations and the state of the repressor protein (free or bound to the operator) [2].

First Model

In the first model, the rate of synthesis GFP is affected by [LacI]. LacI binds cooperatively to the promoter and this relationship if described using a Hill function. Cooperative binding is a suitable model [3] and the promoter pHyper-spank has two LacI binding sites.

[math]\displaystyle{ \frac{d[GFP]}{dt}=\frac{k_4k_{m}^n}{k_m^n+[LacI]^n}-d_2[GFP] }[/math]

The trouble with this model is that [LacI] evolves to a steady state regardless of [math]\displaystyle{ [IPTG]_{t=0} }[/math]. As [math]\displaystyle{ [GFP]_{steady-state} }[/math] depends soley on [math]\displaystyle{ [LacI] }[/math], the paramters of the Hill function and the degradation rate of GFP and is not dependent on [math]\displaystyle{ [IPTG]_{t=0} }[/math], and [math]\displaystyle{ [LacI]_{steady-state} }[/math] is the same for all [math]\displaystyle{ [IPTG]_{t=0} }[/math], changing [math]\displaystyle{ [IPTG]_{t=0} }[/math] affects only the dynamic and not the steady-state behaviour of [math]\displaystyle{ [GFP] }[/math].

Here are the simulations in matlab:

ODEs for first model

Simulation file for first model

And here are four simulations at different [math]\displaystyle{ [IPTG]_{t=0} }[/math]:

$[math]\displaystyle{ [IPTG]_{t=0}=0.1 }[/math]$

[math]\displaystyle{ [IPTG]_{t=0}=0.1 }[/math]
$[math]\displaystyle{ [IPTG]_{t=0}=1 }[/math]$

[math]\displaystyle{ [IPTG]_{t=0}=1 }[/math]
$[math]\displaystyle{ [IPTG]_{t=0}=10 }[/math]$

[math]\displaystyle{ [IPTG]_{t=0}=10 }[/math]
$[math]\displaystyle{ [IPTG]_{t=0}=100 }[/math]$

[math]\displaystyle{ [IPTG]_{t=0}=100 }[/math]

Shown below:

[math]\displaystyle{ [LacI]_{steady-state} }[/math] and [math]\displaystyle{ [GFP]_{steady-state} }[/math] as a function of [math]\displaystyle{ [IPTG]_{t=0} }[/math].

By inspecting the simulations at different [math]\displaystyle{ [IPTG]_{t=0} }[/math] you can see that the maximum [math]\displaystyle{ [GFP] }[/math] reached before [math]\displaystyle{ [GFP] }[/math] attains steady-state is affected by [math]\displaystyle{ [IPTG]_{t=0} }[/math]. In order to characterise the dynamic behaviour of the system these values are also shown. If the wetlab team is able to measure these maximum values this will aid us in model selection.

$[math]\displaystyle{ [LacI]_{steady-state} }[/math] at different [math]\displaystyle{ [IPTG]_{t=0} }[/math]$

[math]\displaystyle{ [LacI]_{steady-state} }[/math] at different [math]\displaystyle{ [IPTG]_{t=0} }[/math]
$[math]\displaystyle{ [GFP]_{steady-state} }[/math] at different [math]\displaystyle{ [IPTG]_{t=0} }[/math]$

[math]\displaystyle{ [GFP]_{steady-state} }[/math] at different [math]\displaystyle{ [IPTG]_{t=0} }[/math]
$Max [math]\displaystyle{ [GFP] }[/math] at different [math]\displaystyle{ [IPTG]_{t=0} }[/math]$

Max [math]\displaystyle{ [GFP] }[/math] at different [math]\displaystyle{ [IPTG]_{t=0} }[/math]

Second Model

In this model, the interactions between IPTG and LacI are as in the first model, but the rate of expression of GFP is related to the concentration of free IPTG and not the concentration of LacI.

We use the "activator" form of the Hill function for this model:

[math]\displaystyle{ \frac{d[GFP]}{dt}=\frac{k_4[IPTG]^n}{k_m^n+[IPTG]^n}-d_2[GFP] }[/math]

ODE

Simulation File

$[math]\displaystyle{ [IPTG]_{t=0}=0.1 }[/math]$

[math]\displaystyle{ [IPTG]_{t=0}=0.1 }[/math]
$[math]\displaystyle{ [IPTG]_{t=0}=1.0 }[/math]$

[math]\displaystyle{ [IPTG]_{t=0}=1.0 }[/math]
$[math]\displaystyle{ [IPTG]_{t=0}=10 }[/math]$

[math]\displaystyle{ [IPTG]_{t=0}=10 }[/math]
$[math]\displaystyle{ [IPTG]_{t=0}=100 }[/math]$

[math]\displaystyle{ [IPTG]_{t=0}=100 }[/math]

$[math]\displaystyle{ [GFP]_{steady-state} }[/math] at different [math]\displaystyle{ [IPTG]_{t=0} }[/math]$

[math]\displaystyle{ [GFP]_{steady-state} }[/math] at different [math]\displaystyle{ [IPTG]_{t=0} }[/math]
$[math]\displaystyle{ [LacI]_{steady-state} }[/math] at different initial [math]\displaystyle{ [IPTG] }[/math]$

[math]\displaystyle{ [LacI]_{steady-state} }[/math] at different initial [math]\displaystyle{ [IPTG] }[/math]
$[math]\displaystyle{ [GFP]_{max} }[/math] at different initial [math]\displaystyle{ [IPTG] }[/math]$

[math]\displaystyle{ [GFP]_{max} }[/math] at different initial [math]\displaystyle{ [IPTG] }[/math]

Characterisation of light-inducible promoter and RBS

Having established the transfer function between inducer concentration and GFP concentration (which will be useful for part characterisation and establishing the most appropriate RBS even if these inducible promoters won't be used in our final constructs) we need to establish:

Whether our model describing the constitutive expression of GFP can be applied to the constitutive expression of YtvA? If we can establish the concentration of GFP in molecules per cell this will be a better measure than arbitrary units/fluorescence but even so, YtvA will likely have a different degradation rate (can we characterise this)?
The relationship between YtvA (over)expression and expression through the sigma-B sensitive promoters. What is a good model for this relationship?

i.e. can we assume
[math]\displaystyle{ \frac{d[ytvA]}{dt}=k_1-d_1[ytvA] }[/math]
making use of the same values of [math]\displaystyle{ k_1 }[/math] and [math]\displaystyle{ d_1 }[/math] that were found from the characterisation of the constitutive promoters?
Would it be more relevant to use
[math]\displaystyle{ \frac{d[ytvA]}{dt}=k_1-d_2[ytvA] }[/math]
i.e. same expression level, different degradation rate, and how to characterise this degradation rate?
We could try the model
[math]\displaystyle{ \frac{d[ytvA]}{dt}=f(\frac{d[GFP]}{dt}) }[/math]
Where [math]\displaystyle{ f }[/math] represents some function to be determined.
Once we have characterised this relationship, we need to characterise the relationship between ytvA expression, light intensity/duration of light exposure and protein expression downstream of the [math]\displaystyle{ \sigma_{B} }[/math] promoter.

Characterisation of EpsE

Characterisation of EpsE has two considerations:

Succeessful expression
Desired functionality

It would be desirable to verify expression of EpsE by using anti-EpsE antibody and an immunohistochemical technique.

However we can functionally characterise EpsE by measuring motility as a function of inducer concentration.

Characterisation of Biomaterials

Ideally we wish to quantitatively characterise Biomaterial expression and secretion to enable us to select the optimal signal peptide - biomaterial combination. We can use our data about GFP expression to estimate expression through the light-inducible system but secretion is more tricky. Hmm.

Light induced expression of EpsE

Light induced expression of biomaterials

All Medline abstracts: PubMed | HubMed

<html><center><table style="color:#ffffff;background-color:#66aadd;" cellpadding="3" cellspacing="1" border="0" bordercolor="#ffffff" align="center"> <tr><td><ul id="sddm"></html>[[IGEM:IMPERIAL/2008/New/{{{1}}}|< Previous]]<html></ul> </td><td><ul id="sddm"><a href="#">Back to top</a></ul> </td><td><ul id="sddm"></html>[[IGEM:IMPERIAL/2008/New/{{{2}}}|Next >]]<html></ul> </td></tr></table> </center></html>

@@ Line 169: / Line 169: @@
 <biblio>
 #1 pmid=18612302
+#3 pmid=7000772
 #2 pmid=2182324
-#3 pmid=7000772
 </biblio>
 {{Imperial/EndPage}}

IGEM:IMPERIAL/2008/Prototype/Drylab/ODE models of test constructs: Difference between revisions

Revision as of 08:28, 12 September 2008

Contents

ODE models of test constructs

Characterisation of Constitutive promoter and RBS

Characterisation of inducible promoter and RBS

First Model

Second Model

Characterisation of light-inducible promoter and RBS

Characterisation of EpsE

Characterisation of Biomaterials

Light induced expression of EpsE

Light induced expression of biomaterials

Navigation menu

Page actions

Page actions

Personal tools

Navigation

Search

research

Tools