Team:NCTU Formosa/RNA discussion

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.

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.)

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.)

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.