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

From 2010.igem.org

(Difference between revisions)
(The Experiment)
(The Model)
Line 13: Line 13:
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 [here] for the more detailed protocol.
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 [here] for the more detailed protocol.
-
===The Model===
+
==The Model==
-
====Previous models====
+
===Previous models===
We've seen that [http://2009.igem.org/Team:Groningen/Promoters previous teams] have performed their protein expression rate analysis to only include current protein level (fluorescence measurements) and current biomass concentration (absorption measurements). However, we've found that during measurements, the protein concentration is also influenced by other factors:
We've seen that [http://2009.igem.org/Team:Groningen/Promoters previous teams] have performed their protein expression rate analysis to only include current protein level (fluorescence measurements) and current biomass concentration (absorption measurements). However, we've found that during measurements, the protein concentration is also influenced by other factors:
* The rapid growth of the bacteria dilutes the protein considerably
* The rapid growth of the bacteria dilutes the protein considerably
Line 21: Line 21:
See below for more on how we approached these issues.
See below for more on how we approached these issues.
-
====Growth curve fitting====
+
===Growth curve fitting===
From the biomass, or optical density (OD), graph below, it could be concluded that the fastest growth occurs from 1 hour to 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. All in all 12 wells were inoculated per strain carrying a different RBS sequence or 'empty' plasmid (I13401).
From the biomass, or optical density (OD), graph below, it could be concluded that the fastest growth occurs from 1 hour to 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. All in all 12 wells were inoculated per strain carrying a different RBS sequence or 'empty' plasmid (I13401).
Line 32: Line 32:
Which resulted in a growth rate ('mu') which was implemented in the dilution factor in the further calculations.
Which resulted in a growth rate ('mu') which was implemented in the dilution factor in the further calculations.
-
====Biomass growth will dilute GFP====
+
===Biomass growth will dilute GFP===
[[Image:TUD2010_Gfpod_all.png|right|400px]]
[[Image:TUD2010_Gfpod_all.png|right|400px]]
The conventional method for characterization of RBS or promotor sequences is to assume that the derivative of the GFP fluorescence is equal to the production rate, as shown in the formula from Kelly et al [1].
The conventional method for characterization of RBS or promotor sequences is to assume that the derivative of the GFP fluorescence is equal to the production rate, as shown in the formula from Kelly et al [1].

Revision as of 16:32, 23 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 this 'strength' we are attempting to characterize?

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.

The Model

Previous models

We've seen that previous teams have performed their protein expression rate analysis to only include current protein level (fluorescence measurements) and current biomass concentration (absorption measurements). However, we've found that during measurements, the protein concentration is also influenced by other factors:

  • The rapid growth of the bacteria dilutes the protein considerably
  • Proteins degrade over time.

See below for more on how we approached these issues.

Growth curve fitting

From the biomass, or optical density (OD), graph below, it could be concluded that the fastest growth occurs from 1 hour to 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. All in all 12 wells were inoculated per strain carrying a different RBS sequence or 'empty' plasmid (I13401).

Tud2010 RBS OD.png

The exponential section of this graph was thus curve-fitted to an exponential growth model:

TUD2010 RBS Growth model.PNG

Which resulted in a growth rate ('mu') which was implemented in the dilution factor in the further calculations.

Biomass growth will dilute GFP

TUD2010 Gfpod all.png

The conventional method for characterization of RBS or promotor sequences is to assume that the derivative of the GFP fluorescence is equal to the production rate, as shown in the formula from Kelly et al [1]. TUDelft2010 RBS 1754-1611-3-4-i9.gif

However, this method would only be representable if the strains are subcultured multiple times from the exponential phase. In that case, a steady state GFP concentration (fluorescence) may be assumed and dilution is ignored. Due to the decrease in the biomass specific GFP fluorescence at the start of the measurements it was assumed that dilution effects played an important role in our data and thus was taken into account for our model.


Protein production model

As discussed earlier, the protein production model will take dilution into account. Furthermore, the GFP protein that is used is very stable (half life 41 hours), thus the degradation term is negligible compared to dilution.

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:

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 results page.

Data normalization method

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

TUDelft2010 Rbs-control-gfp.png

Our plate reader data consists of 12 columns of

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

To prepare the data for further calculation, some things have to be 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
  • 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.
  • For each well, we calculate the real GFP that is not a result of biomass:
real_gfp = gfp - fake_gfp_per_od * od
  • These values are curve fitted into our protein production model, resulting in a (exponential phase) protein production rate for each well.

Source code and data

http://2010.igem.org/wiki/images/8/88/TUDelft_2010_RBS_characterization.zip

References

[1] Jason R. Kelly, Adam J. Rubin, Joseph H. Davis, et al (March 2009). "Measuring the activity of BioBrick promoters using an in vivo reference standard". Journal of Biological Engineering 3: 4

CharacterizationResultsParts