Team:NCTU Formosa/RNA discussion

From 2011.igem.org



RNA Thermometer

Modeling

A low-temperature release system was designed and constructed in the Pathway Commander to control target protein expression. A specific ribosome binding site (RBS) BBa_K115002 with high translation activity at high temperature(> 37°C) and low translation activity at room temperature was used to design the temperature-dependent genetic circuit in E. coli, with a green fluorescent protein (GFP) used as the reporter protein(Figure1.) . We analyzed fluorescence intensity during E. coli growth at log phase and stationary phase at temperatures 25°C, 30°C, 37°C and 40°C. Our experimental results indicate that high temperature decreased the translation rate of the target protein, and this temperature-dependent genetic circuit can control the expression level of the target protein by the host cell's incubation temperature. However, the translational activity of the RBS BBa_K115002 at different temperatures cannot be quantified directly from experimental data. To overcome this problem, we provided a dynamic model which can quantitatively assess the translation strength of the RBS BBa_K115002 at temperatures 25°C, 30°C, 37°C and 40°C.


Circuit A: BBa_J23101 is a constitutive promoter, and the 37℃induced RBS we use is BBa_K115002. Ptet BBa _R0040 is a constitutive repressible promoter, and the RBS we use is BBa_B0034. The mGFP is the gene coding for “Green Fluorescence Protein”. The terminators are all BBa_J61048.

Low-temperature release system consists of two genes: tetR and GFP. The expression level of TetR protein is the input to the system, and the concentration of GFP is the output.
This system can be modeled by differential equations as follows.

The first equation describes the temperature control in strand A (Figure 1.). Alpha-Temp is the production rates corresponding to transcriptional rate of constitutive promoter and the translation rate of the RBS BBa_K115002which is a temperature sensitive post-transcriptional regulator. The second equation describes the concentration of GFP change with time. Alpha-B is production rates of the GFP, which are assumed to be given constants. To describe transition during log phase and stationary phase, the alpha-Temp and alpha-B and is assumed to zero when the Terminator in stationary phase. Gamma-TetR, and gamma-GFP are decay rates of the corresponding proteins. When bacteria divide, the molecular in a bacterium will be dilute. Because bacteria grow faster, the dilution rate d(t) is included in this model and can be calculated from OD(optical density) ratio of medium (Figure 2.). For an inhibition of TetR protein, Hill function is an S-shaped curve which can be described in the form 1 / (1 +x^n) (Alon, 2007). The values of the kinetic parameters used in the simulation were initially obtained from the literature and experimental data. Data computations were performed with Matlab software. A program was written and used as a subroutine in Matlab for parameter optimization using nonlinear regression (Figure 3.)



Figure 2: The reporter circuit:
The OD ratio is increased faster in log phase than it in stationary phase. The dilution rate d(t) can be calculated from OD ratio and used in our model.



Figure 3: The behavior of low temperature release circuit at 25°C, 30°C, 37 °C and 40°C. Experimental data (dot) and simulated results (line) of the model suggest this temperature-dependent genetic circuit can control the expression level of the target protein by the host cell’s incubation. The fitting results indicate our dynamic model can quantitatively assess the relative translational activity of RBS during log phase and stationary phase.

Using least squares estimation from experimental data, the relative translational activity of this RBS BBa_K115002at 25°C, 30°C, 37 °C and 40°C were estimated. (Figure 4.)



Figure 4: The relative translation activity of this RBS (BBa_K115002) at 25°C, 30°C, 37°C and 40°C estimated using least squares estimation from experimental data. This means our dynamic model can accurately quantify the translational activity of the RBS from experimental data.

According to the fitting results (Figure 3.), the dynamic model successfully approximated the behavior of our low-temperature release system. The model equations present interesting mathematical properties that can be used to explore how qualitative features of the genetic circuit depend on reaction parameters. This method of dynamic modeling can be used to guide the choice of genetic ‘parts’ for implementation in circuit design in the future.