Team:TU Delft/Project/rbs-characterization/characterization

From 2010.igem.org

(Difference between revisions)
(Characterization)
(Data normalization method)
 
(44 intermediate revisions not shown)
Line 3: Line 3:
<html><center><img src="https://static.igem.org/mediawiki/2010/0/00/TU_Delft_project_navigation.jpg" usemap="#projectnavigation" border="0" /></center><map id="projectnavigation" name="projectnavigation"><area shape="rect" alt="Characterization" title="" coords="309,3,591,45" href="https://2010.igem.org/Team:TU_Delft#page=Project/rbs-characterization/characterization" target="" /><area shape="rect" alt="Results" title="" coords="609,3,891,44" href="https://2010.igem.org/Team:TU_Delft#page=Project/rbs-characterization/results" target="" /><area shape="rect" alt="Parts" title="" coords="9,3,290,44" href="https://2010.igem.org/Team:TU_Delft#page=Project/rbs-characterization/parts" target="" /></map></html>
<html><center><img src="https://static.igem.org/mediawiki/2010/0/00/TU_Delft_project_navigation.jpg" usemap="#projectnavigation" border="0" /></center><map id="projectnavigation" name="projectnavigation"><area shape="rect" alt="Characterization" title="" coords="309,3,591,45" href="https://2010.igem.org/Team:TU_Delft#page=Project/rbs-characterization/characterization" target="" /><area shape="rect" alt="Results" title="" coords="609,3,891,44" href="https://2010.igem.org/Team:TU_Delft#page=Project/rbs-characterization/results" target="" /><area shape="rect" alt="Parts" title="" coords="9,3,290,44" href="https://2010.igem.org/Team:TU_Delft#page=Project/rbs-characterization/parts" target="" /></map></html>
==Characterization==
==Characterization==
-
<i>Summary</i>
+
<i>By varying only the RBS sequence in a protein generator construct the dependence of the protein expression level on the RBS sequence could be determined by simple fluorescence and biomass curve analysis. The methods are further explained below.</i>
-
By varying only the RBS sequence in a protein generator construct the dependence of the expression level of the encoded protein on the RBS sequence could be determined by simple fluorescence and biomass curve analysis. The methods are further explained below.
+
-
=== What is this 'strength' we're attempting to characterize? ===
+
=== What is RBS strength, really ? ===
In the field of synthetic biology the 'strength' of an RBS is defined as the rate at which translation is initiated by the RBS sequence on any given mRNA molecule. Going back a step, we see that translation initiation occurs preferentially at certain mRNA sequences, which show a similarity to the consensus -Shine-Delgarno- sequence. This is due to the optimal binding of the 16S rRNA at these regions. In other words, the RBS 'strength' may loosely be defined as the rate of ribosome binding to any given mRNA molecule. (Click [http://partsregistry.org/Help:Ribosome_Binding_Sites/Mechanism here] for more information.)
In the field of synthetic biology the 'strength' of an RBS is defined as the rate at which translation is initiated by the RBS sequence on any given mRNA molecule. Going back a step, we see that translation initiation occurs preferentially at certain mRNA sequences, which show a similarity to the consensus -Shine-Delgarno- sequence. This is due to the optimal binding of the 16S rRNA at these regions. In other words, the RBS 'strength' may loosely be defined as the rate of ribosome binding to any given mRNA molecule. (Click [http://partsregistry.org/Help:Ribosome_Binding_Sites/Mechanism here] for more information.)
For our project it is not directly the binding equilibria of ribosome to mRNA that is of interest to us, but rather the netto rate of protein production. Thus, in our experimentation we define the RBS strength as the production rate of a given protein downstream of an RBS.
For our project it is not directly the binding equilibria of ribosome to mRNA that is of interest to us, but rather the netto rate of protein production. Thus, in our experimentation we define the RBS strength as the production rate of a given protein downstream of an RBS.
-
===Experimental Setup===
+
==The experiment==
-
We took five different RBS sequences from the [http://partsregistry.org/Ribosome_Binding_Sites/Prokaryotic/Constitutive/Anderson Anderson RBS family] (J61100, J61101, J61107, J61117, J61127). The standard RBS [http://partsregistry.org/Part:BBa_B0032 B0032] was used as a reference for our characterization.  
+
The <i>E.coli strains</i> carrying the parts to be tested were cultured for 16 hours on 96-well plates while measuring fluorescence of GFP and biomass absorbance at 10 minute intervals. Click [https://2010.igem.org/Team:TU_Delft/Protocols#RBS_characterization_experiment here] for the more detailed protocol.
-
All these RBS sequences were placed in front of the standard [http://partsregistry.org/Part:BBa_E0040 GFP coding sequence], so expression could be measured.
 
-
The Biobricks generated in order to perform the experiments were: [http://partsregistry.org/Part:BBa_K398500 K398500], [http://partsregistry.org/Part:BBa_K398501 K398501], [http://partsregistry.org/Part:BBa_K398502 K398502], [http://partsregistry.org/Part:BBa_K398503 K398503], [http://partsregistry.org/Part:BBa_K398504 K398504]. The general map of the construction is shown below, where the RBS is displayed in fucsia.
 
-
 
-
[[Image:RBS1.jpg|left|650px]]
 
-
 
-
{| style="color:black; background-color:white;" cellpadding="5" cellspacing="0" border="1"
 
-
|'''Feature'''
 
-
|'''Function'''
 
-
|-
 
-
|[http://partsregistry.org/Part:BBa_B0032 AmpR]
 
-
|Ampicillin resistance
 
-
|-
 
-
|[http://partsregistry.org/Part:BBa_B0015 B0015]
 
-
|Transcriptional (double) terminator
 
-
|-
 
-
|[http://partsregistry.org/Part:BBa_B0062 B0062]
 
-
|Transcriptional terminator
 
-
|-
 
-
|[http://partsregistry.org/Part:BBa_E0040 E0040]
 
-
|GFP
 
-
|-
 
-
|[http://partsregistry.org/Part:BBa_G00000 G00000]
 
-
|Standard prefix
 
-
|-
 
-
|[http://partsregistry.org/Part:BBa_G00001 G00001]
 
-
|Standard suffix
 
-
|-
 
-
|[http://partsregistry.org/Part:BBa_G00100 G00100]
 
-
|VF2 primer binding site
 
-
|-
 
-
|[http://partsregistry.org/Part:BBa_G00101 G00101]
 
-
|VR primer binding site
 
-
|-
 
-
|[http://partsregistry.org/Part:BBa_J61100 J61100]
 
-
|RBS Anderson family
 
-
|-
 
-
|[http://partsregistry.org/Part:BBa_J23100 J23100]
 
-
|Promoter
 
-
|}
 
-
 
-
<html><div style='clear:both'>
 
-
</html>
 
-
 
-
=== Known methods ===
 
-
We've seen that previous teams have performed their RBS characterization measurements only include current protein level (GFP measurements) and current biomass concentration (Absorption measurements). During measurements, the protein concentration is also influenced by other factors:
 
-
* The rapid growth of the bacteria will dilute the protein considerably
 
-
* Proteins degrade over time.
 
-
 
-
The conventional method for characterization of RBS or promotor sequences is to assume that the derivative of GFP concentration is the production rate, as shown in the formula from Kelly et al (shown below).
 
-
[[Image:TUDelft2010_RBS_1754-1611-3-4-i9.gif]]
 
-
 
-
This method only works if cells are subcultured a few times in the exponential phase. In that case, steady state GFP concentration can be assumed and dilution can be ignored.
 
-
In our experiment however, this was not the case. Specific GFP is actually decreasing because of dilution, so it needs to be modeled in order to make any use of the data.
 
-
 
-
=== Growth curve fitting ===
 
-
From the biomass, or optical density (OD), graph below, it can be seen that the fastest growth occurs from 1 hour until about 3 hours into the experiment. Within this timespan, it is assumed that growth is exponential, and growth rate can be calculated. The graph below shows 72 growth curves. 12 seperate wells were used for every RBS sequence.
 
-
 
-
[[Image:Tud2010_RBS_OD.png]]
 
-
 
-
The exponential part of this graph is curve-fitted to an exponential growth model:
 
-
 
-
[[Image:TUD2010_RBS_Growth_model.PNG]]
 
-
 
-
Which results in a growth rate ('mu') that can be used to add dilution to further calculation.
 
-
 
-
=== Biomass growth will dilute the GFP ===
 
-
[[Image:TUD2010_Gfpod_all.png|right|400px]]
 
-
At the start of our measurements, the GFP protein has reached a steady state high concentration.
 
-
The specific GFP will decrease during growth because of dilution, and this needs to be accounted for when calculating the RBS strength.
 
-
 
-
<br style="clear: both" />
 
-
 
-
== Protein production model ==
 
-
As said, the protein production model needs to take dilution into account. The GFP protein that is used is very stable (half life 41 hours), so the degradation term is negligible compared to dilution.
 
-
 
-
[[Image:RBS_Expression_model.PNG]]
 
-
 
-
In above formula:
 
-
* Y is the GFP concentration
 
-
* Beta is the production rate (The constant that we are interested in)
 
-
* Alpha is the dilution+degradation rate
 
-
 
-
Converting this equation to an explicit form results in:
 
-
 
-
[[Image:TU2010_Gfp_explicit.PNG]]
 
-
 
-
The 2 unknown parameters (beta and Y1) where curve-fitted against the GFP/OD data, using the matlab lsqcurvefit command.
 
-
The beta constant for each of the 72 experiments were then to calculate RBS strengths, listed in the  [[Team:TU_Delft/Project/rbs-characterization/results|results]] page.
 
=== Data normalization method ===
=== Data normalization method ===
-
To go from original measurements to RBS strenghts, we ran the following steps using a matlab script:
+
To go from the raw measurements to the RBS strengths, we ran the following steps using a matlab script:
[[Image:TUDelft2010_Rbs-control-gfp.png|right|thumb|300px]]
[[Image:TUDelft2010_Rbs-control-gfp.png|right|thumb|300px]]
-
Our plate reader data consists of 12 columns of
+
the plate reader data consists of 8 rows, each containing multiples of 12:
-
* 5 strains with Anderson RBS'es to be characterized
+
* 5 strains with Anderson RBSs to be characterized
* B0032 strain to compare to
* B0032 strain to compare to
-
* I13401 strain as control (GFP without promotor and RBS)
+
* I13401 strain as control (GFP-dT without promotor and RBS)
* LB+Amp as blank
* LB+Amp as blank
-
To prepare the data for further calculation, some things have to be done:
+
To prepare the data for further calculation, the following was done:
-
* Subtract the blank GFP and OD values from the data with E coli in it. This takes care of the platereader error independent of biomass
+
* Subtract the blank GFP and OD values from the data of wells containing E coli. This takes care of any background noise caused by the LB medium.
-
* Biomass without GFP also produces some nonzero measurements (See right). From the control strain measurements we can calculate how much 'fake' GFP is reported due to biomass.  
+
* The strain containing I13401 produces some nonzero measurements (See right). From this control strain measurements we can calculate how much autofluorescence is reported due to biomass or due to leaky expression.
-
* For each well, we calculate the real GFP that is not a result of biomass:
+
* For each well, we calculate the GFP that is not a result of biomass or leaky expression:
-
real_gfp = gfp - fake_gfp_per_od * od
+
[[Image:TUDelft_2010_PPM_Autofluorescence.png]]
-
* These values are curve fitted into our protein production model, resulting in a (exponential phase) protein production rate for each well.
+
* These values are curve fitted into our [[Team:TU_Delft/Modeling/protein-production-model|protein production model]], described in the In Silico section. This results in a (exponential phase) protein production rate for each plate reader well. These production rates are then averaged and compared to the B0032 reference RBS to produce our [[Team:TU_Delft/Project/rbs-characterization/results|resulting RBS strengths]].
-
 
+
-
=== Source code and data ===
+
-
https://static.igem.org/mediawiki/2010/8/88/TUDelft_2010_RBS_characterization.zip
+
-
 
+
-
=== References ===
+
-
;Kelly 2009
+
-
:Jason R. Kelly, Adam J. Rubin, Joseph H. Davis, ''et al'' (March 2009). "[http://dx.doi.org/10.1186/1754-1611-3-4 Measuring the activity of BioBrick promoters using an in vivo reference standard]". ''Journal of Biological Engineering'' '''3''': 4
+
<html><center><img src="https://static.igem.org/mediawiki/2010/0/00/TU_Delft_project_navigation.jpg" usemap="#projectnavigation" border="0" /></center><map id="projectnavigation" name="projectnavigation"><area shape="rect" alt="Characterization" title="" coords="309,3,591,45" href="https://2010.igem.org/Team:TU_Delft#page=Project/rbs-characterization/characterization" target="" /><area shape="rect" alt="Results" title="" coords="609,3,891,44" href="https://2010.igem.org/Team:TU_Delft#page=Project/rbs-characterization/results" target="" /><area shape="rect" alt="Parts" title="" coords="9,3,290,44" href="https://2010.igem.org/Team:TU_Delft#page=Project/rbs-characterization/parts" target="" /></map></html>
<html><center><img src="https://static.igem.org/mediawiki/2010/0/00/TU_Delft_project_navigation.jpg" usemap="#projectnavigation" border="0" /></center><map id="projectnavigation" name="projectnavigation"><area shape="rect" alt="Characterization" title="" coords="309,3,591,45" href="https://2010.igem.org/Team:TU_Delft#page=Project/rbs-characterization/characterization" target="" /><area shape="rect" alt="Results" title="" coords="609,3,891,44" href="https://2010.igem.org/Team:TU_Delft#page=Project/rbs-characterization/results" target="" /><area shape="rect" alt="Parts" title="" coords="9,3,290,44" href="https://2010.igem.org/Team:TU_Delft#page=Project/rbs-characterization/parts" target="" /></map></html>

Latest revision as of 19:47, 27 October 2010

CharacterizationResultsParts

Characterization

By varying only the RBS sequence in a protein generator construct the dependence of the protein expression level on the RBS sequence could be determined by simple fluorescence and biomass curve analysis. The methods are further explained below.

What is RBS strength, really ?

In the field of synthetic biology the 'strength' of an RBS is defined as the rate at which translation is initiated by the RBS sequence on any given mRNA molecule. Going back a step, we see that translation initiation occurs preferentially at certain mRNA sequences, which show a similarity to the consensus -Shine-Delgarno- sequence. This is due to the optimal binding of the 16S rRNA at these regions. In other words, the RBS 'strength' may loosely be defined as the rate of ribosome binding to any given mRNA molecule. (Click here for more information.)

For our project it is not directly the binding equilibria of ribosome to mRNA that is of interest to us, but rather the netto rate of protein production. Thus, in our experimentation we define the RBS strength as the production rate of a given protein downstream of an RBS.

The experiment

The E.coli strains carrying the parts to be tested were cultured for 16 hours on 96-well plates while measuring fluorescence of GFP and biomass absorbance at 10 minute intervals. Click here for the more detailed protocol.


Data normalization method

To go from the raw measurements to the RBS strengths, we ran the following steps using a matlab script:

TUDelft2010 Rbs-control-gfp.png

the plate reader data consists of 8 rows, each containing multiples of 12:

  • 5 strains with Anderson RBSs to be characterized
  • B0032 strain to compare to
  • I13401 strain as control (GFP-dT without promotor and RBS)
  • LB+Amp as blank

To prepare the data for further calculation, the following was done:

  • Subtract the blank GFP and OD values from the data of wells containing E coli. This takes care of any background noise caused by the LB medium.
  • The strain containing I13401 produces some nonzero measurements (See right). From this control strain measurements we can calculate how much autofluorescence is reported due to biomass or due to leaky expression.
  • For each well, we calculate the GFP that is not a result of biomass or leaky expression:

TUDelft 2010 PPM Autofluorescence.png

  • These values are curve fitted into our protein production model, described in the In Silico section. This results in a (exponential phase) protein production rate for each plate reader well. These production rates are then averaged and compared to the B0032 reference RBS to produce our resulting RBS strengths.

CharacterizationResultsParts