Team:UCSF/Modeling Results

Though we don’t know the exact parameters in the experiments, and sometimes even cannot tune every parameter in reality, the more robust a dynamic system is, the easier the experiments could succeed.

To test the robustness of this dynamic system, we sampled the parameters using Latin-cubic sampling, and searched the suitable groups of parameters which could make the population ratio of the two strains (x=n1/n2) get to a constant number, and at the same time could also enable the final steady ratio to be neither too large nor too small, say, not larger than 10, and not smaller than 0.1.

Finally we found 261 out of 10000 groups of parameters that meet our requirement. This number gives an evaluation of the robustness.

What’s more, to give guidance on how to tune parameters to make experiments work out, we should find out which parameter is sensitive, say, which one to tune and how to tune it to make an experiment successful.

Due to the reasons above, we did sensitivity analysis using C++, the results are as following:

The result of sensitivity analysis indicates that k1 or k2 should be tuned if we fail to get a steady population ratio of the two strains. But this result leaves us a dilemma in tuning the parameters in reality in experiments, because k_1 and k_2, the maximum growth rates of Escherichia coli (gDM/(Ls)), are the natural characteristics of the cells, which makes this situation intractable. We could tune α and β to get the final steady population ratio we want, but if the system couldn’t even go to a steady ratio, that doesn’t make sense in reality. That’s why a new approach should be proposed.

Result of Simulation

To test our analytical result and the robustness of this system, we simulated the system in Matlab with different groups of parameters:

k1=4, k2=3.9, K1=1, K2=1, α1=1, α2=4, β1=1, β2=1
Setting the initial values to [0.1, 0.2, 0, 0], [0.1, 0.5, 0, 0], [0.1, 1.0, 0, 0], [0.3, 0.2, 0, 0], [0.3, 0.5, 0, 0], [0.3, 1.0, 0, 0], [0.7, 0.2, 0, 0], [0.7, 0.5, 0, 0], [0.7, 1.0, 0, 0] respectively, we follow the values of x of each line as time goes on, and a phase plane is shown below

However, when the population ratio of the two strains gets to a steady state, y (c_1/c_2 ) doesn’t necessarily be constant. As an example, with another group of parameters

k1=101,k2=100,K1=1e^(-5),K2=1e^(-5),α1=1e^(-4),α2=1e^(-4),β1=1e^(-6),β2=1e^(-6),
the behavior of this dynamic system changes. Though the ratio of the two strains and c_1 come to a steady state finally, c_2 increases continuously as time goes on. Setting the initial values to [0.1, 0.2, 0, 0], [0.1, 0.5, 0, 0], [0.1, 1.0, 0, 0], [0.3, 0.2, 0, 0], [0.3, 0.5, 0, 0], [0.3, 1.0, 0, 0], [0.7, 0.2, 0, 0], [0.7, 0.5, 0, 0], [0.7, 1.0, 0, 0] respectively ( The initial values are the same as the former situation.), we follow the values of x of each line as time goes on, and a phase plane is shown below.

Retrieved from "http://2012.igem.org/Team:UCSF/Modeling_Results"