ETHZ/Simulation

From 2007.igem.org

(Difference between revisions)
(Sensitivity Analysis)
(Sensitivity Analysis)
Line 41: Line 41:
There for we simulated the timing behavior while sweeping the parameters by 5%. In detail we did the following for e.g. test case 1:
There for we simulated the timing behavior while sweeping the parameters by 5%. In detail we did the following for e.g. test case 1:
-
For each paramter we did the timing simulation like shown in Figure 1 once while increasing the paramter by 5% and once while decreasing the parameter by 5%. Since in test case 1 we are interested in having a high concentration of GFP during the recognition phase with respect to the base line we divide the steady state GFP concentration at 4000 minutes by the baseline concentration at 1000 minutes. Finally as shown in the following equation we define the sensitivity by subtracting the resulting factor that we get with an increase parameter from the resulting factor that we get by decreasing the paramter.
+
For each paramter we did the timing simulation like shown in Figure 1 once while increasing the paramter by 5% and once while decreasing the parameter by 5%. Since in test case 1 we are interested in having a high concentration of GFP during the recognition phase with respect to the base line we divide the steady state GFP concentration at 4000 minutes by the baseline concentration at 1000 minutes. Finally as shown in the following equation we define the sensitivity by subtracting the resulting factor that we get with an increase parameter from the resulting factor that we get by decreasing the parameter.
[[Image:ETHZSensititivityequation.png|center|320px]]
[[Image:ETHZSensititivityequation.png|center|320px]]

Revision as of 12:34, 15 October 2007

Contents

Protocol

Fig. 1: Flow diagram

Simulation of Test Cases

To test for all conditions we run timing simulations for all 4 different branches in the flow diagram of our protocol and plotted the concentrations of the inducers, memory proteins and reporter proteins.

  • During the first 1000 minutes of the simulation we did not add any inducers to the system to check for the baseline production of proteins.
  • Between 1000 and 2000 minutes we added the inducer that should be learned and let the system reach steady state.
  • After the 2000th minute we added the inducer AHL for memory formation and run the simulation for additional 1000 minutes to reach again steady state.
  • Between 3000 and 4000 minutes we finally tested if the system reacts in the desired way by showing the right color.

For the simulation parameters that have been used please see ADD LINK TO PARAMETER TABLE

Test 1: Learn and recognize IPTG

Fig 1: Test 1: learn and recognize IPTG

Figure 1 shows the simulated behavior of the system when IPTG is presented both during the learing phase and the recognition phase. In both the learning and recognition phase the system reports by producing green florescence proteins which matches the desired behavior.

Test 2: Learn IPTG and get tested for aTc

Fig 2: Test 2: learn IPTG and get tested for aTc

Figure 2 shows the simulated behavior of the system when IPTG is presented only during the learing phase but aTc during the recognition phase. In contrast to Figure 1 the system reports by producing red florescence proteins during the recognition phase which matches the desired behavior.

Test 3: Learn and recognize aTc

Fig 3: Test 3: learn and recognize aTc

Figure 3 shows the simulated behavior of the system when aTc is presented both during the learing phase and the recognition phase. In both the learning and recognition phase the system reports by producing yellow florescence proteins which matches the desired behavior.

Test 4: Learn aTc and get tested for IPTG

Fig 4: Test 4: learn aTc and get tested for IPTG

Figure 4 shows the simulated behavior of the system when aTc is presented only during the learing phase but IPTG during the recognition phase. In contrast to Figure 3 the system reports by producing cyan florescence proteins during the recognition phase which matches the desired behavior.

Sensitivity Analysis

We were not only interested in finding a set of parameters for which our system is working. To provide a stable performance also in the presence of noise and parameter missmatch we were interested in analysing the sensitivity of the system to find the most sensitive parameters.

There for we simulated the timing behavior while sweeping the parameters by 5%. In detail we did the following for e.g. test case 1: For each paramter we did the timing simulation like shown in Figure 1 once while increasing the paramter by 5% and once while decreasing the parameter by 5%. Since in test case 1 we are interested in having a high concentration of GFP during the recognition phase with respect to the base line we divide the steady state GFP concentration at 4000 minutes by the baseline concentration at 1000 minutes. Finally as shown in the following equation we define the sensitivity by subtracting the resulting factor that we get with an increase parameter from the resulting factor that we get by decreasing the parameter.

ETHZSensititivityequation.png

So if there is no difference in the concentration of the desired output during the recognition phase with respect to baseline when sweeping a parameter by +/-5% the sensitivity of that parameter is equal to zero. Otherwise the sensitivity gives us a measure of how much the concentration changes.

We quantified the sensitivity for all parameters and test cases. The results are plotted in Figure 5.

Fig. 5: Sensitivity analysis

We found the following results:

  • Our system is most sensitve to the parameter a which descibes the base production. This is not a surpising result since we depend on having a good switch and our signal to noise or with other words output to baseline production level highly depends on the fact that we can reliably switch off the production of the memory proteins cI and p22cII. For a more detailed analysis of the switching behavior please see ADD LINK TO ANALYSIS!!!
  • Other sensitive parameters are those that are directly related to the production and decay of the memory proteins cI and p22cII. This can be explained by the fact that our system prefers symmetry between the parameter sets of cI and p22cII to support good switching behavior where the same concentration of p22cII leads to the same repression of cI than vice versa. We can especially see that those parameters seem to be more sensitive for the test cases 1 and 2 where cI has to be stored inside the memory. The reason for this behavior is that due to the missing symmetry between the parameters for cI and p22cII we currently have a bias inside the system towards the production of p22cII.