"
Team:NCTU Formosa/modeling
From 2011.igem.org
Measurement
Modeling and simulations of high temperature induced device BBa_K098995 – cI promoter & cI repressor
In order to characterize this high temperature induced device BBa_K098995, the fluorescence intensity ofBBa_K098988 is measured by the flow cytometry (Figure. 1).
Figure1. Part BBa_K098988 Design. The heat induced device BBa_K098995uses gene BBa_K098997 coding for cI repressor to inhibit the cI promoter BBa_R0051. The activity of cI repressor is decreased by elevating temperature from 30 ℃ to 42 ℃. A differential equation is used to calculate protein expression activity of BBa_K098995 as follows.
This equation describes the concentration of GFP in BBa_K098988change with time (Figure. 1). Alpha-Temp is the protein expression rates corresponding to BBa_K098995which is a temperature sensitive expression device. To describe transition during log phase and stationary phase, the alpha-Temp is assumed to zero in stationary phase. Gamma-GFP are decay rates of the GFP 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 ratio of medium (Figure. 2). 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 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 out model.
Figure 3. The behavior of high temperature induced device BBa_K098988 at 25°C, 37 °C and 42°C. Experimental data (dot) and simulated results (line) of the model suggest this temperature-dependent device 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 protein expression activity of BBa_K098988during log phase and stationary phase.
Using least squares estimation from experimental data, the relative the protein expression activity of BBa_K098988 at 25°C, 37 °C and 42°C were estimated (Figure. 4).
Figure 4. The relative the protein expression activity of BBa_K098988at 25°C, 37 °C and 42°C estimated using least squares estimation from experimental data. The protein expression activity at 42°C is higher than 25°C, 37 °C
According to the fitting results (Figure. 3), the dynamic model successfully approximated the behavior of our high-temperature induced system. The model equation presents 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.
References
Alon, U. (2007) An Introduction to Systems Biology: Design Principles of Biological Circuits. Chapman & Hall/CRC.
