Team:Colombia Uniandes/Parameters

From 2013.igem.org

Parameters

The second step in the deterministic model is to find the parameters that define the system. There are 3 possible ways to do this:

  1. Literature: Many of the genes that we use in iGEM have been characterised before, and the constants that model their behaviour have been found. Sometimes an exact value is not found, so we have to find at least the order of magnitude or a range where the parameter can be found
  2. Experiments: This is the most reliable way of knowing the parameters. In our case, we need them to perform the deterministic model first, in order to see if the system works as expected in the lab. If there is an estimated order of magnitude the model can be run while the exact value is found experimentally
  3. Screening: If the parameters we are looking of are not found in literature we proceed to perform a screening of parameters as described in Colombia 2012 Team Wiki

Below is shown how we found the parameters for both of the projects:

Glucocorticoid sensor

After looking for parameters in literature, the following values were found:

ParamGC.png

  • The values for these parameters varied a lot depending on the glucocorticoid detected. As we plan to test the circuit with dexamethasone, all the parameters were adjusted for this hormone.
  • The values in red are the parameters that could not be found. An estimated order of magnitude was assessed supported in the literature read.
  • As the part of the chimera that is going to activate transcription is Gal4p, the values for the activation of the gene are assumed to be the same as those of the pure protein.
  • Some of the destruction rates were not found so we assumed that the concentration of the protein is significantly reduced only when the cell divides.
Sensitivity Analysis

In order to see which of the parameters had a significant effect in the signal response, we performed a sensitivity analysis. For each parameter, we took the value found in the literature (or guessed value) and evaluated the model in a range that varied from one order of magnitude below to one order of magnitude above. Depending on how much the response varied to each parameter the screening was planned.

After perfoming the analysis we found some parameters that had little effect on the desired response of the system. For example, the destruction rate of glucocorticoids inside the cell and the reverse kinetic constant for the formation of the complex:

Nosens.png

The parameters related to the Hill expression had a significant effect on the desired response. For this reason, and given the fact that two out of three parameters had an estimated order of magnitude, the screening of parameters was focused on this part of the model.

Hillsens.png
Kmsens.png


You can see the complete sensitivity analysis here.

The parameters that need to be found with the model are those involved in the activation of the promoter by GAL4. Before this promoter exists an Upstream Activation Sequence (UAS) that enhances the production of M-CHERRY. This sequence can be as many times as wanted upstream de promoter, and depending on the number of times it exists the strength of the production is affected. Little is found in the literature about the parameters that concern to the activation by GAL4 and less is found in the difference of the parameters with the different arrangements of the UAS.

Without the parameter information of at least one arrangement of UAS the screening of parameters can not be performed. In order to solve this problem we designed an experiment, in which we tested the M-cherry's expression trough time with different UAS arrangements (1,3,6 and 9 UAS before de promoter). With the results of this procedure we aim to find an expression for the activation of the promoter pGAL1 with different UAS and its parameters. THIS EXPERIMENT CURRENTLY IS BEING PERFORMED AT OUR LABS!


Nickel removal System

After a research for parameters in the literature research was done the following values were found:

ParamNic.png

  • The values in red are the parameters that could not be found. An estimated order of magnitude was assessed supported in the literature read.
  • Some of the destruction rates are not found so it is assumed the concentration of the protein is significantly reduced only when the cell divides
Sensitivity Analysis

A complete sensitivity analysis was performed in order to see which parameters have a significant importance in the model. Each parameter was evaluated in various orders of magnitude and the behaviour of the nickel intake was analysed.

Most of the parameters had a significant effect. Only the rate constant for the formation of the tetramer and the repressor destruction rate were the parameters that compared with the others had little effect on the nickel intake. The first one needed to be found with the model, but after seeing that it has not a significant effect on the desired response an estimated order of magnitude is enough to run the deterministic model.

Repressordestructionses.png

The repressor basal production rate is one of the parameters that needs to be found with the model. As can be seen below this parameter has a strong influence on the response of the system in a small range (0.8-1.3) . After 1.4 the response of the system is not influenced by the change of the basal production. The parameter can be chosen to be higher than 1.4 because it has no significant influence in the response and the intake of nickel is high enough.

Basalprodrepanalysis.png

As said before most of the parameters showed to have influence in the model. The parameters related to the Hill's equation: The Maximum expression rate and the cooperation constant had the highest effect on the model and its value needs to be screened.

EffectNickel.png


Here is the complete document with the results for all the parameters