Team:KULeuven/Model/Filter

From 2008.igem.org

(Difference between revisions)
(3. Proof of filtering capacities)
m
 
(27 intermediate revisions not shown)
Line 1: Line 1:
 +
{{:Team:KULeuven/Tools/Styling}}
 +
{{:Team:KULeuven/Tools/Scripting}}
{{:Team:KULeuven/Tools/Header}}
{{:Team:KULeuven/Tools/Header}}
-
<div style="float: right;">[[Image:pictogram_filter.png|120px]]</div>
+
[[Image:pictogram_filter.png|120px|right]]
== Filter ==
== Filter ==
Line 19: Line 21:
<body>
<body>
<p>
<p>
-
<a href="https://static.igem.org/mediawiki/2008/a/a4/Filter.pdf">
+
<a href="https://static.igem.org/mediawiki/2008/2/2a/FilterODE.pdff">
<img border="0" src="https://2008.igem.org/wiki/skins/common/images/icons/fileicon-pdf.png" width="65" height="60">
<img border="0" src="https://2008.igem.org/wiki/skins/common/images/icons/fileicon-pdf.png" width="65" height="60">
</a>
</a>
Line 40: Line 42:
| 0.00155 s<sup>-1</sup>
| 0.00155 s<sup>-1</sup>
| UmuD tag added to speed up degradation of otherwise too stable T7 polymerase
| UmuD tag added to speed up degradation of otherwise too stable T7 polymerase
-
| [http://www.openwetware.org/wiki/IGEM:Tsinghua/2007/Projects/RAP#Model_and_simulation link]
+
| [https://2008.igem.org/Team:KULeuven/Model/Filter#References [4<html>]</html> [5<html>]</html> [8<html>]</html>]
|-
|-
| d<sub>mRNA RIBOKEY</sub>
| d<sub>mRNA RIBOKEY</sub>
| 0.00462 s<sup>-1</sup>
| 0.00462 s<sup>-1</sup>
| estimate: because this RNA isn't translated, it degrades faster
| estimate: because this RNA isn't translated, it degrades faster
-
|[http://www.pubmedcentral.nih.gov/picrender.fcgi?artid=124983&blobtype=pdf link]
+
|[https://2008.igem.org/Team:KULeuven/Model/Filter#References [3<html>]</html>]
|-
|-
| d<sub>closed mRNA T7</sub>
| d<sub>closed mRNA T7</sub>
| 0.00462 s<sup>-1</sup>
| 0.00462 s<sup>-1</sup>
| estimate: because this mRNA isn't translated, it degrades faster
| estimate: because this mRNA isn't translated, it degrades faster
-
|[http://www.pubmedcentral.nih.gov/picrender.fcgi?artid=124983&blobtype=pdf link]
+
|[https://2008.igem.org/Team:KULeuven/Model/Filter#References [3<html>]</html>]
|-
|-
| d<sub>open mRNA T7</sub>
| d<sub>open mRNA T7</sub>
| 0.00231 s<sup>-1</sup>
| 0.00231 s<sup>-1</sup>
|  
|  
-
|[http://www.pubmedcentral.nih.gov/picrender.fcgi?artid=124983&blobtype=pdf link]
+
|[https://2008.igem.org/Team:KULeuven/Model/Filter#References [3<html>]</html>]
|-
|-
| d<sub>open mRNA T7 complex</sub>
| d<sub>open mRNA T7 complex</sub>
| 0.00231 s<sup>-1</sup>
| 0.00231 s<sup>-1</sup>
|  
|  
-
|[http://www.pubmedcentral.nih.gov/picrender.fcgi?artid=124983&blobtype=pdf link]
+
|[https://2008.igem.org/Team:KULeuven/Model/Filter#References [3<html>]</html>]
|-
|-
! colspan="4" style="border-bottom: 1px solid #003E81;" | Equilibrium constants
! colspan="4" style="border-bottom: 1px solid #003E81;" | Equilibrium constants
Line 66: Line 68:
| K<sub>eq 1</sub>
| K<sub>eq 1</sub>
| 0,015 [M]
| 0,015 [M]
-
| between closed and open T7 mRNA, experimental
+
| between closed and open T7 mRNA, models competition, experimental
-
| [http://parts2.mit.edu/wiki/index.php/Berkeley2006-RiboregulatorsMain link]
+
| [https://2008.igem.org/Team:KULeuven/Model/Filter#References [1<html>]</html> [2<html>]</html>]
|-
|-
| K<sub>eq 2</sub>
| K<sub>eq 2</sub>
| 0.0212 [M]
| 0.0212 [M]
-
| between closed T7 mRNA and key unlocked mRNA complex, experimental
+
| between closed T7 mRNA and key unlocked mRNA complex, models competition, experimental
-
| [http://parts2.mit.edu/wiki/index.php/Berkeley2006-RiboregulatorsMain link]
+
| [https://2008.igem.org/Team:KULeuven/Model/Filter#References [1<html>]</html> [2<html>]</html>]
|-
|-
! colspan="4" style="border-bottom: 1px solid #003E81;" | Rate constants
! colspan="4" style="border-bottom: 1px solid #003E81;" | Rate constants
Line 78: Line 80:
| k<sub>dis</sub>
| k<sub>dis</sub>
| 0.00416 s<sup>-1</sup>
| 0.00416 s<sup>-1</sup>
-
| derived from experimental values
+
| estimate derived from experimental values
-
| [http://parts2.mit.edu/wiki/index.php/Berkeley2006-RiboregulatorsMain link]
+
| [https://2008.igem.org/Team:KULeuven/Model/Filter#References [1<html>]</html>]
|-
|-
| k<sub>complex</sub>
| k<sub>complex</sub>
| 0.00237 s<sup>-1</sup>
| 0.00237 s<sup>-1</sup>
-
| derived from experimental values
+
| estimate derived from experimental values
-
| [http://parts2.mit.edu/wiki/index.php/Berkeley2006-RiboregulatorsMain link]
+
| [https://2008.igem.org/Team:KULeuven/Model/Filter#References [1<html>]</html>]
|-
|-
| k<sub>closed</sub>
| k<sub>closed</sub>
| 500 s<sup>-1</sup>
| 500 s<sup>-1</sup>
-
| derived from experimental values
+
| estimate derived from experimental values
-
| [http://parts2.mit.edu/wiki/index.php/Berkeley2006-RiboregulatorsMain link]
+
| [https://2008.igem.org/Team:KULeuven/Model/Filter#References [1<html>]</html>]
|-
|-
| k<sub>open</sub>
| k<sub>open</sub>
| 7.5 s<sup>-1</sup>
| 7.5 s<sup>-1</sup>
-
| derived from experimental values
+
| estimate derived from experimental values
-
| [http://parts2.mit.edu/wiki/index.php/Berkeley2006-RiboregulatorsMain link]
+
| [https://2008.igem.org/Team:KULeuven/Model/Filter#References [1<html>]</html>]
|-
|-
| k<sub>translation</sub>
| k<sub>translation</sub>
| 0.167 s<sup>-1</sup>
| 0.167 s<sup>-1</sup>
-
| lock defined translation rate for T7 RNA pol
+
| estimate: lock defined translation rate for T7 RNA pol
-
| [http://partsregistry.org/Part:BBa_B0032 link]
+
| [https://2008.igem.org/Team:KULeuven/Model/Filter#References [7<html>]</html>]
|-
|-
! colspan="4" style="border-bottom: 1px solid #003E81;" | Transcription rates
! colspan="4" style="border-bottom: 1px solid #003E81;" | Transcription rates
Line 106: Line 108:
| p(TetR) dependent  
| p(TetR) dependent  
| (RiboKey) between 5E-5 and 0.0125 s<sup>-1</sup> ~ [aTc]
| (RiboKey) between 5E-5 and 0.0125 s<sup>-1</sup> ~ [aTc]
-
|
+
| [https://2008.igem.org/Team:KULeuven/Model/Filter#References [6<html>]</html>]
|-
|-
| k<sub>mRNA T7</sub>
| k<sub>mRNA T7</sub>
| 0,0011 s<sup>-1</sup>
| 0,0011 s<sup>-1</sup>
| weak constitutive promoter J23109
| weak constitutive promoter J23109
-
| [http://partsregistry.org/Part:BBa_J23109 link]
+
| [https://2008.igem.org/Team:KULeuven/Model/Filter#References [7<html>]</html>]
|}
|}
Line 129: Line 131:
[[Image:Sim_filter_1.png|800px|center|filter]]
[[Image:Sim_filter_1.png|800px|center|filter]]
-
 
-
[[Image:Filter 1.PNG|right]]
 
====1. AND gate of the filter====
====1. AND gate of the filter====
In the simulation we can clearly see this series of events:
In the simulation we can clearly see this series of events:
 +
[[Image:Filter 1.PNG|right]]
* when dark blue(ribokey) starts to increase, red (T7) also starts to increase, giving rise to an increasing amount of lactonase (blue) = AND-GATE.
* when dark blue(ribokey) starts to increase, red (T7) also starts to increase, giving rise to an increasing amount of lactonase (blue) = AND-GATE.
* when dark blue(ribokey) starts to decrease, red (T7) also starts to decrease, but much slower. The lactonase also starts to decrease, as it should be.
* when dark blue(ribokey) starts to decrease, red (T7) also starts to decrease, but much slower. The lactonase also starts to decrease, as it should be.
Line 157: Line 158:
====3. Proof of filtering capacities====
====3. Proof of filtering capacities====
-
As proof of the filtering capacities, we'll compare the maximum amount of GFP with the maximum amount of Lactonase as a function of the pulse length of TetR. We've chosen Lactonase and GFP because they only differ in the absence or presence of the filter. The used model for the proof consists only of the filter, the output and the lactonase productions subsystems.
+
As proof of the filtering capacities, we'll compare the maximum amount of GFP with the maximum amount of Lactonase as a function of the pulse length of TetR. We've chosen Lactonase and GFP because they differ only in the absence or presence of the filter, all other parameters are identical (eg. degradation rate). The model used here consists only of the filter, the output and the lactonase productions subsystems.
-
The left figures shows us a clear difference between the amount of GFP and lactonase for almost every pulse length. The explanation for the (almost constant) big difference for pulses longer than 1000 seconds is the following: For such long pulses, the filter is fully active such that the only differences between the output and the lactonase production are the limitations by the ribokey and the difference in transcription factors.  
+
The left figures shows us a clear difference between the amount of GFP and lactonase for almost every pulse length. The explanation for the (almost constant) big difference for pulses longer than 1000 seconds is the following. For such long pulses, the filter can be considered as fully active, such that the only differences between the output and the lactonase production are characteristics of the ribokey/lock and the difference in transcription rates between E. coli's RNA polymerase II and phage T7's polymerase.  
-
For a pulse of less than 1000 seconds the filter is not fully active. This results in a filtering effect. If we zoom for only the first 1000 seconds (right figures), we see that the amount of Lactonase remains more or less the same while the amount of GFP increases very fast.
+
For a pulse shorter than 1000 seconds the filter is not fully active. This creates a filtering effect. If we just look at the first 1000 seconds (right figures), we see that the amount of Lactonase remains more or less the same while the amount of GFP increases very fast.
<html>
<html>
<div class="center">
<div class="center">
Line 169: Line 170:
</div></div></div></html>
</div></div></div></html>
-
Another way to see the filtering effect is to look at the ratio of GFP with Lactonase. This figure clearly shows the disproportionality between GFP and lactonase for the pulses shorter than 1000 seconds.
+
<br>
-
[[Image:GFP_Lactonase.png|center|550px]]
+
-
=== Sensitivity Analysis ===
+
Another way to see the filtering effect is to look at the ratio of GFP to Lactonase. This figure clearly shows the disproportionality between GFP and lactonase for the pulses shorter than 1000 seconds. The best filtering effect occurs for pulses around 100 seconds. This figure also shows that placing the filter in front of a signal wil lower this signal with a factor of at least 25 (full throughput in the filter). This base factor should be multiplied by ~8 for pules around 100 seconds. This is the filtering effect quantitated: pulses of only 100 seconds will produce almost 25 * 8 = 200 times more GFP than Lactonase. A factor 25 characteristic of the system and a variable extra filtering factor that is a function of pulse length and is maximal (8) for pulses of 100 seconds.
-
[[Image:Sens Filter.png|center]]
+
<html>
 +
<div class="center">
 +
<div class="noborder" style="overflow: auto; width: 900px; height: 370px;">
 +
<div class="noborder" style="width: 850px;">
 +
<img src="https://static.igem.org/mediawiki/2008/0/0c/GFP_Lactonase.png" style="float: left; width: 400px; height: 350px; margin: 0 5px;" />
 +
<img src="https://static.igem.org/mediawiki/2008/1/13/GFP_lactonase_norm.png" style="float: left; width: 400px; height: 350px; margin: 0 5px;" />
 +
</div></div></div></html>
 +
 
 +
<br>
 +
 
 +
Another nice way to visualazise the same filtering effect can be seen in the following figures. In the left figure, 4 input pulses are given in sequence. The first lasts 100 seconds, the second 500, the third 2500 while the last one remains indefinateley.
 +
The maximal output (result of the 4th pulse), both non-filtered and filtered are set to 1, normalizing the signals. It is now very easy to see the filtering in effect on the zoomed-in right figure and the ratio's can be compared with the y-values of the previous figure.
 +
 
 +
<html>
 +
<div class="center">
 +
<div class="noborder" style="overflow: auto; width: 900px; height: 320px;">
 +
<div class="noborder" style="width: 900px;">
 +
<img src="https://static.igem.org/mediawiki/2008/1/18/Filter_extra.png" style="float: left; width: 440px; height: 300px; margin: 0 5px;" />
 +
<img src="https://static.igem.org/mediawiki/2008/0/0b/Filter_extra_zoom.png" style="float: left; width: 440px; height: 300px; margin: 0 5px;" />
 +
</div></div></div></html>
 +
 
 +
<br>
 +
 
 +
=== References ===
 +
<html xmlns="http://www.w3.org/1999/xhtml" xml:lang="en">
 +
<head>
 +
<meta http-equiv="Content-Type" content="text/html; charset=utf-8"/>
 +
<title>Bibliography</title>
 +
</head>
 +
<body>
 +
<table style="border-collapse:collapse;line-height:1.1em;">
 +
<tr style="vertical-align:top;"><td>[1]</td><td style="padding-left:4pt;">“Berkeley2006-RiboregulatorsMain - IGEM”; http://parts2.mit.edu/wiki/index.php/Berkeley2006-RiboregulatorsMain.</td></tr>
 +
<tr><td colspan="2">&nbsp;</td></tr>
 +
<tr style="vertical-align:top;"><td>[2]</td><td style="padding-left:4pt;">F.J. Isaacs et al., “Engineered riboregulators enable post-transcriptional control of gene expression,” <span style="font-style:italic;">Nat Biotech</span>,  vol. 22, Jul. 2004, pp. 841-847. <span class="Z3988" title="url_ver=Z39.88-2004&amp;ctx_ver=Z39.88-2004&amp;rft_id=info%3Adoi/10.1038/nbt986&amp;rft_val_fmt=info%3Aofi%2Ffmt%3Akev%3Amtx%3Ajournal&amp;rft.genre=article&amp;rft.atitle=Engineered%20riboregulators%20enable%20post-transcriptional%20control%20of%20gene%20expression&amp;rft.jtitle=Nat%20Biotech&amp;rft.stitle=Nat%20Biotech&amp;rft.volume=22&amp;rft.issue=7&amp;rft.aufirst=Farren%20J&amp;rft.aulast=Isaacs&amp;rft.au=Farren%20J%20Isaacs&amp;rft.au=Daniel%20J%20Dwyer&amp;rft.au=Chunming%20Ding&amp;rft.au=Dmitri%20D%20Pervouchine&amp;rft.au=Charles%20R%20Cantor&amp;rft.au=James%20J%20Collins&amp;rft.date=2004-07&amp;rft.pages=841-847&amp;rft.issn=1087-0156"></span></td></tr>
 +
 
 +
<tr><td colspan="2">&nbsp;</td></tr>
 +
<tr style="vertical-align:top;"><td>[3]</td><td style="padding-left:4pt;">J.A. Bernstein et al., “Global analysis of mRNA decay and abundance in Escherichia coli at single-gene resolution using two-color fluorescent DNA microarrays,” <span style="font-style:italic;">Proceedings of the National Academy of Sciences of the United States of America</span>,  vol. 99, Jul. 2002, pp. 9697–9702. <span class="Z3988" title="url_ver=Z39.88-2004&amp;ctx_ver=Z39.88-2004&amp;rft_id=info%3Adoi/10.1073/pnas.112318199&amp;rft_val_fmt=info%3Aofi%2Ffmt%3Akev%3Amtx%3Ajournal&amp;rft.genre=article&amp;rft.atitle=Global%20analysis%20of%20mRNA%20decay%20and%20abundance%20in%20Escherichia%20coli%20at%20single-gene%20resolution%20using%20two-color%20fluorescent%20DNA%20microarrays&amp;rft.jtitle=Proceedings%20of%20the%20National%20Academy%20of%20Sciences%20of%20the%20United%20States%20of%20America&amp;rft.stitle=Proc%20Natl%20Acad%20Sci%20U%20S%20A.%20&amp;rft.volume=99&amp;rft.issue=15&amp;rft.aufirst=Jonathan%20A.&amp;rft.aulast=Bernstein&amp;rft.au=Jonathan%20A.%20Bernstein&amp;rft.au=Arkady%20B.%20Khodursky&amp;rft.au=Pei-Hsun%20Lin&amp;rft.au=Sue%20Lin-Chao&amp;rft.au=Stanley%20N.%20Cohen&amp;rft.date=2002-07-23&amp;rft.pages=9697%E2%80%939702"></span></td></tr>
 +
<tr><td colspan="2">&nbsp;</td></tr>
 +
<tr style="vertical-align:top;"><td>[4]</td><td style="padding-left:4pt;">“IGEM:Tsinghua/2007/Projects/RAP - OpenWetWare”; http://www.openwetware.org/wiki/IGEM:Tsinghua/2007/Projects/RAP#Model_and_simulation.</td></tr>
 +
<tr><td colspan="2">&nbsp;</td></tr>
 +
<tr style="vertical-align:top;"><td>[5]</td><td style="padding-left:4pt;">M. Gonzalez et al., “Lon-mediated proteolysis of the Escherichia coli UmuD mutagenesis protein: in vitro degradation and identification of residues required for proteolysis,” <span style="font-style:italic;">Genes Dev.</span>,  vol. 12, Dec. 1998, pp. 3889-3899. <span class="Z3988" title="url_ver=Z39.88-2004&amp;ctx_ver=Z39.88-2004&amp;rft_id=info%3Adoi/10.1101/gad.12.24.3889&amp;rft_val_fmt=info%3Aofi%2Ffmt%3Akev%3Amtx%3Ajournal&amp;rft.genre=article&amp;rft.atitle=Lon-mediated%20proteolysis%20of%20the%20Escherichia%20coli%20UmuD%20mutagenesis%20protein%3A%20in%20vitro%20degradation%20and%20identification%20of%20residues%20required%20for%20proteolysis&amp;rft.jtitle=Genes%20Dev.&amp;rft.volume=12&amp;rft.issue=24&amp;rft.aufirst=Martin&amp;rft.aulast=Gonzalez&amp;rft.au=Martin%20Gonzalez&amp;rft.au=Ekaterina%20G.%20Frank&amp;rft.au=Arthur%20S.%20Levine&amp;rft.au=Roger%20Woodgate&amp;rft.date=1998-12-15&amp;rft.pages=3889-3899"></span></td></tr>
 +
<tr><td colspan="2">&nbsp;</td></tr>
 +
 
 +
<tr style="vertical-align:top;"><td>[6]</td><td style="padding-left:4pt;">M. Bon, S.J. McGowan, and P.R. Cook, “Many expressed genes in bacteria and yeast are transcribed only once per cell cycle,” <span style="font-style:italic;">FASEB J.</span>,  vol. 20, Aug. 2006, pp. 1721-1723. <span class="Z3988" title="url_ver=Z39.88-2004&amp;ctx_ver=Z39.88-2004&amp;rft_id=info%3Adoi/10.1096/fj.06-6087fje&amp;rft_val_fmt=info%3Aofi%2Ffmt%3Akev%3Amtx%3Ajournal&amp;rft.genre=article&amp;rft.atitle=Many%20expressed%20genes%20in%20bacteria%20and%20yeast%20are%20transcribed%20only%20once%20per%20cell%20cycle&amp;rft.jtitle=FASEB%20J.&amp;rft.volume=20&amp;rft.issue=10&amp;rft.aufirst=Michael&amp;rft.aulast=Bon&amp;rft.au=Michael%20Bon&amp;rft.au=Simon%20J.%20McGowan&amp;rft.au=Peter%20R.%20Cook&amp;rft.date=2006-08-01&amp;rft.pages=1721-1723"></span></td></tr>
 +
<tr><td colspan="2">&nbsp;</td></tr>
 +
<tr style="vertical-align:top;"><td>[7]</td><td style="padding-left:4pt;">“Part:BBa J23109 - partsregistry.org”; http://partsregistry.org/Part:BBa_J23109.</td></tr>
 +
<tr><td colspan="2">&nbsp;</td></tr>
 +
<tr style="vertical-align:top;"><td>[8]</td><td style="padding-left:4pt;">J.P. McDonald et al., “Regulation of UmuD cleavage: role of the amino-terminal tail,” <span style="font-style:italic;">Journal of Molecular Biology</span>,  vol. 282, Oct. 1998, pp. 721-730. <span class="Z3988" title="url_ver=Z39.88-2004&amp;ctx_ver=Z39.88-2004&amp;rft_id=info%3Adoi/10.1006/jmbi.1998.2044&amp;rft_val_fmt=info%3Aofi%2Ffmt%3Akev%3Amtx%3Ajournal&amp;rft.genre=article&amp;rft.atitle=Regulation%20of%20UmuD%20cleavage%3A%20role%20of%20the%20amino-terminal%20tail&amp;rft.jtitle=Journal%20of%20Molecular%20Biology&amp;rft.volume=282&amp;rft.issue=4&amp;rft.aufirst=John%20P&amp;rft.aulast=McDonald&amp;rft.au=John%20P%20McDonald&amp;rft.au=Erinn%20E%20Maury&amp;rft.au=Arthur%20S%20Levine&amp;rft.au=Roger%20Woodgate&amp;rft.date=1998-10-02&amp;rft.pages=721-730"></span></td></tr>
 +
<tr><td colspan="2">&nbsp;</td></tr>
 +
</table></body>
 +
 
 +
</html>

Latest revision as of 23:12, 29 October 2008

  dock/undock dropdown  

Pictogram filter.png

Contents

Filter

Position in the system

The filter is positioned immediately after the input, because its job is to filter out possible noise signals or background signals that aren't caused by the "desease". It is the starting piece of the whole system, situated before the invertimer- and the reset-subsystem.

Describing the system

see also: Project:Filter

Filter BioBrick.jpg

ODE's

Parameters

Parameter values Filter
Name Value Comments Reference
Degradation rates
dpT7 tag 0.00155 s-1 UmuD tag added to speed up degradation of otherwise too stable T7 polymerase [4] [5] [8]
dmRNA RIBOKEY 0.00462 s-1 estimate: because this RNA isn't translated, it degrades faster [3]
dclosed mRNA T7 0.00462 s-1 estimate: because this mRNA isn't translated, it degrades faster [3]
dopen mRNA T7 0.00231 s-1 [3]
dopen mRNA T7 complex 0.00231 s-1 [3]
Equilibrium constants
Keq 1 0,015 [M] between closed and open T7 mRNA, models competition, experimental [1] [2]
Keq 2 0.0212 [M] between closed T7 mRNA and key unlocked mRNA complex, models competition, experimental [1] [2]
Rate constants
kdis 0.00416 s-1 estimate derived from experimental values [1]
kcomplex 0.00237 s-1 estimate derived from experimental values [1]
kclosed 500 s-1 estimate derived from experimental values [1]
kopen 7.5 s-1 estimate derived from experimental values [1]
ktranslation 0.167 s-1 estimate: lock defined translation rate for T7 RNA pol [7]
Transcription rates
TetR_var_transcr_rate p(TetR) dependent (RiboKey) between 5E-5 and 0.0125 s-1 ~ [aTc] [6]
kmRNA T7 0,0011 s-1 weak constitutive promoter J23109 [7]

Remark: The key-lock system has been enhanced to 0.3%-14% (new parameters have been added)

Models

CellDesigner (SBML file)

filter

Matlab (SBML file)

Filter Matlab.jpg

Simulations

filter

1. AND gate of the filter

In the simulation we can clearly see this series of events:

Filter 1.PNG
  • when dark blue(ribokey) starts to increase, red (T7) also starts to increase, giving rise to an increasing amount of lactonase (blue) = AND-GATE.
  • when dark blue(ribokey) starts to decrease, red (T7) also starts to decrease, but much slower. The lactonase also starts to decrease, as it should be.

The short lifetime of the ribokey compared to the lifetime of the T7-protein, guarantees that the AND-GATE always works perfectly fine: when there's no more input, the ribokey will rapidly decrease (and disappear) and makes sure that the AND-GATE is not activated anymore, even when the T7-protein is slowly decreasing.

2. Filtering in practice

In the simulation three kinds of inputpulses have been used:

  • First pulse: 300s
    • small peak of lactonase
    • no influence on the timer
  • Second pulse: 1000s
    • medium peak of lactonase
    • influences the timer by levelling the timing capabilities, but it doesn't reset the timer
  • Third pulse: 5000s
    • huge peak of lactonase
    • reset of the timer: amount of complex goes to zero

3. Proof of filtering capacities

As proof of the filtering capacities, we'll compare the maximum amount of GFP with the maximum amount of Lactonase as a function of the pulse length of TetR. We've chosen Lactonase and GFP because they differ only in the absence or presence of the filter, all other parameters are identical (eg. degradation rate). The model used here consists only of the filter, the output and the lactonase productions subsystems.

The left figures shows us a clear difference between the amount of GFP and lactonase for almost every pulse length. The explanation for the (almost constant) big difference for pulses longer than 1000 seconds is the following. For such long pulses, the filter can be considered as fully active, such that the only differences between the output and the lactonase production are characteristics of the ribokey/lock and the difference in transcription rates between E. coli's RNA polymerase II and phage T7's polymerase. For a pulse shorter than 1000 seconds the filter is not fully active. This creates a filtering effect. If we just look at the first 1000 seconds (right figures), we see that the amount of Lactonase remains more or less the same while the amount of GFP increases very fast.


Another way to see the filtering effect is to look at the ratio of GFP to Lactonase. This figure clearly shows the disproportionality between GFP and lactonase for the pulses shorter than 1000 seconds. The best filtering effect occurs for pulses around 100 seconds. This figure also shows that placing the filter in front of a signal wil lower this signal with a factor of at least 25 (full throughput in the filter). This base factor should be multiplied by ~8 for pules around 100 seconds. This is the filtering effect quantitated: pulses of only 100 seconds will produce almost 25 * 8 = 200 times more GFP than Lactonase. A factor 25 characteristic of the system and a variable extra filtering factor that is a function of pulse length and is maximal (8) for pulses of 100 seconds.


Another nice way to visualazise the same filtering effect can be seen in the following figures. In the left figure, 4 input pulses are given in sequence. The first lasts 100 seconds, the second 500, the third 2500 while the last one remains indefinateley. The maximal output (result of the 4th pulse), both non-filtered and filtered are set to 1, normalizing the signals. It is now very easy to see the filtering in effect on the zoomed-in right figure and the ratio's can be compared with the y-values of the previous figure.


References

Bibliography

[1]“Berkeley2006-RiboregulatorsMain - IGEM”; http://parts2.mit.edu/wiki/index.php/Berkeley2006-RiboregulatorsMain.
 
[2]F.J. Isaacs et al., “Engineered riboregulators enable post-transcriptional control of gene expression,” Nat Biotech, vol. 22, Jul. 2004, pp. 841-847.
 
[3]J.A. Bernstein et al., “Global analysis of mRNA decay and abundance in Escherichia coli at single-gene resolution using two-color fluorescent DNA microarrays,” Proceedings of the National Academy of Sciences of the United States of America, vol. 99, Jul. 2002, pp. 9697–9702.
 
[4]“IGEM:Tsinghua/2007/Projects/RAP - OpenWetWare”; http://www.openwetware.org/wiki/IGEM:Tsinghua/2007/Projects/RAP#Model_and_simulation.
 
[5]M. Gonzalez et al., “Lon-mediated proteolysis of the Escherichia coli UmuD mutagenesis protein: in vitro degradation and identification of residues required for proteolysis,” Genes Dev., vol. 12, Dec. 1998, pp. 3889-3899.
 
[6]M. Bon, S.J. McGowan, and P.R. Cook, “Many expressed genes in bacteria and yeast are transcribed only once per cell cycle,” FASEB J., vol. 20, Aug. 2006, pp. 1721-1723.
 
[7]“Part:BBa J23109 - partsregistry.org”; http://partsregistry.org/Part:BBa_J23109.
 
[8]J.P. McDonald et al., “Regulation of UmuD cleavage: role of the amino-terminal tail,” Journal of Molecular Biology, vol. 282, Oct. 1998, pp. 721-730.