"
Team:Valencia Biocampus/modeling1
From 2012.igem.org
Alejovigno (Talk | contribs) (→Modeling Bacteria) |
Alejovigno (Talk | contribs) |
||
| Line 5: | Line 5: | ||
==='''Are you hungry?'''=== | ==='''Are you hungry?'''=== | ||
| + | |||
| + | The first models we did is for our glucose-sensitive construction. We carried out three different experiments changing the carbon source, to see its influence. First we used glucose as the only carbon source, then we added galactose as an extra carbon source and then we put sodium acetate as an extra carbon source instead of galactose. | ||
| + | This is because we checked in previous tests that low concentrations of glucose compromised cell growth. | ||
| + | With the experimental data we fit the following repressible hill function. | ||
| + | [[Image:ecu_glu.png]] | ||
| + | |||
| + | And we estimated the following parameters for the 3 different conditions. | ||
| + | [[Image:table_glu.png]] | ||
Revision as of 19:57, 26 September 2012
Contents |
Modeling Bacteria
Here we modeled several of the Bacteria constructions.
Are you hungry?
The first models we did is for our glucose-sensitive construction. We carried out three different experiments changing the carbon source, to see its influence. First we used glucose as the only carbon source, then we added galactose as an extra carbon source and then we put sodium acetate as an extra carbon source instead of galactose.
This is because we checked in previous tests that low concentrations of glucose compromised cell growth.
With the experimental data we fit the following repressible hill function.
And we estimated the following parameters for the 3 different conditions.
Are you hot?
For the case of the Heat Shock promoter, we decided to model the static transfer function of the Fluorescent Intensity as a result of the Temperature with the following equation:
In order to estimate the parameters of this model we performed a heat shock experiment.
Analysis corresponding to a heat shock applied during 10 minutes at different temperatures ranging from 20ºC up to 51ºC. For each temperature three samples were tested. Fluorescence Intensity (FI) and optical density (OD) were measured at t1=60 minutes and t2=120 minutes after application of the heat pulse had finished.
The first practical problem we encountered was the lack of enough heating coats. Because of these we took the culture, and we splitted the experiments in two batches. FI and OD were measured at t1=60, and t2=120 for batch number 1, while for batch number 2 the heat shock was delayed by 40 minutes with respect to the first batch. Therefore, measu-rements for batch no.2 were taken at t1=100, and t2=160 minutes respectively.
During this 40 minutes delay, the microorganism was growing. Therefore, the heat shock in batch no. 2 was applied to a larger population than in the case of batch no.1. Consequently, the measurements of FI and DO gave much higher values. This affects the ratio FI/DO, were the two batches can be clearly distinguished (see figure 1).
 |
| Figure 1. Raw data from the heat shock experiment. |
Our first modeling work was to perform data reconciliation. To this end, we assume that during this 40 minutes delay the only thing that happended was that cells grew. Growth was approximately linear in the range analyzed, as seen in figure 1 (blue data). We fitted a simple model DO=p1+p2*t. The result can be seen in figure 2, were: p1 = 0.161573 p2 = 0.000866314 Therefore, the ratio FI/DO above can be compensated by multiplying data by the frac-tions (p1+p2*t100)/(p1+p2*t), and (p1+p2*t160)/(p1+p2*t). Actually, what we did was to compensate all measurements to the same time base. The reconcialited result can be seen in figure 3.
 |
| Figure 2. Model fitting for OD w.r.t. time. |
 |
| Figure 3. Recoinciliated data for t=60min and t=120mins. |
Once we had all data reconciliated we fitted it to a static model relating the ratio be-ween fluorescence intensity (FI), and the culture optical density (OD) with the temperature applied during the heat shock.
Since the model is nonlinear with respect to parameters K1, and n, we applied a global op-timization algorithm based on stochastic search. Actually we applied it to fit the data at t=60 minutes, and t=120 minutes independently. The estimated parameters were:
The results can be seen in figures 4 (t60), and 5 (t120) respectively.
 |
| Figure 4. Model and experimental data for t=60min. |
 |
| Figure 5. Model and experimental data for t=120min. |
Compare this results with Are you hot? results
Modeling Yeast
For our yeast construct: How long have you been fermenting? we propose a phenomenological model to explain the relationship between glucose concentration and fluorescent protein production. The following model was assumed:
where
The graph below shows a simulated theoretical time evolution for the three variables.
As observed in the figure, once the glucose is depleted, the concentration of fluorescent protein increases from a basal value. Compare this results with Yeast fluorescence results
