Team:TU Delft/Modeling/HC regulation/Sensitivity

From 2010.igem.org

Revision as of 14:09, 27 October 2010 by Lbergwerff (Talk | contribs)

(diff) ← Older revision | Latest revision (diff) | Newer revision → (diff)

Sensitivity analysis

For all beta parameters (lumped variables for transcriptions speed, translations speed, promoter strength and rbs strength) and all K parameters (equilibrium constant for binding of AlkS with alkanes and coefficients for the activity of promoters) a sensitivity analysis has been performed. This was done by varying the parameters from 50% of their original values to 200% of their original values.

For this range of parameter values it is checked how much the steady state concentration of AlkB is and how long it takes for 95% of the steady state concentration of AlkB to be reached from steady state values for no alkanes to steady state values for 1 µM of alkanes.

In the figures below are the results for this sensitivity analysis.

Figure 1 – Sensitivity analysis of parameter β^max_S1

Figure 2 – Sensitivity analysis of parameter β^max_S2

Figure 3 – Sensitivity analysis of parameter β^max_B

Figure 4 – Sensitivity analysis of parameter β^max_AlkS

Figure 5 – Sensitivity analysis of parameter β^max_AlkB

Figure 6 – Sensitivity analysis of parameter K_HC

Figure 7 – Sensitivity analysis of parameter K_S1

Figure 8 – Sensitivity analysis of parameter K_S2

Figure 9 – Sensitivity analysis of parameter K_B

Conclusion

Most parameters have a linear or slightly exponential effect on the steady state concentration. The time it takes to reach steady state is more complexly correlated to the parameters. The wave functions for β^max_B and β^max_AlkB can be explained by the numerical integration steps. The number of time steps of the integration is finite, through interpolation the value the concentration reaches 95% of the steady state concentration is approximated. This can result however in numerical artifacts as displayed by these two functions.

These results are summarized in table 1. A '0' means insensitive, a '+' means somewhat sensitive and a '++' is very sensitive.

Table 1; results for sensitivity analysis for the gene regulation model

Parameter	Sensitive for steady state concentration	Sensitive for reaching steady state concentration
β^max_S1	0	+
β^max_S2	+	+
β^max_B	++	0
β^max_AlkS	+	+
β^max_AlkB	++	0
K_HC	+	+
K_S1	0	+
K_S2	+	+
K_B	+	+

β^max_S1 and β^max_S1 have the strongest effect on the steady state concentration of AlkB. This is logical as they directly influence the concentrations of AlkB and the mRNA of AlkB in the ODE. Interestingly they have no significant influence on how fast the system reaches steady state.

The second observation that can be made from the results is that β^max_S1 and K_S1 have no significant influence on the steady state concentration of AlkB. This is logical as AlkB is expressed with higher alkane concentrations. In this situation pAlkS1 is more repressed and does not contribute much at the total AlkS and AlkB concentrations.

Navigation

Go back to the model explanation here

Go back to the main modeling page here