Team:Wageningen UR/Project/ModelingProj1

From 2011.igem.org

Revision as of 20:17, 21 September 2011 by Suzette (Talk | contribs)

Building a Synchronized Oscillatory System

Modeling synchronized oscillations


Mathematical model of the Hasty construct

Since our BioBricked oscillatory system is based on the circuit published by Danino et al. in the paper “A synchronized quorum of genetic clocks” [1], our first model of the system is a reproduction of the mathematical model in the supplementary information [2] accompanying the publication mentioned. In their simulations, Danino et al. used a set of four delay differential equations, which we also used as starting point for our modeling work.

The steps (transcription, translation, maturation etc.) from the luxI and aiiA genes to the corresponding proteins are not modeled separately. Instead, the delay of the correlation between the internal AHL concentration, which triggers the expression of the genes, and the corresponding AiiA and LuxI concentrations is simulated by a Hill function. This Hill function takes the history of the system into account, i.e. the concentration of AHL at the time it binds to LuxR to form the activation complex.



back to top


Writing a modeling tool in Matlab

For the first simulations, the same parameters were used as in the cited paper, the only variables thus being the cell density and flow rate. To get an idea about the different conditions in which oscillations could occur, our team created a script for a Matlab modeling tool which uses nested for-loops to vary the flow rate and cell density over a range of values. The resulting tool allows the user to enter the range in which the variables should be varied. The script then iterates over the values and plots graphs of all combinations possible for that range of values. Figure 2 shows an example output of the script.


Example output graphs WUR.png

Fig.2: Variation of output graphs depending on the different starting conditions


back to top


Preliminary conclusions for our system

The first observation from the model was that, for oscillations to occur, the flow rates may not be too fast, especially at lower cell densities (see output graphs above). Since the device used for our system has larger dimensions than the microfluidic device used by Danino et al. the flow rates required could not be achieved by varying height differences alone.

Furthermore, an interesting observation was that applying a flow rate over the cells was not essential to obtain oscillations. Figure 3 shows how, according to the model, oscillations can potentially occur at 0 flow rate. Interesting is that this happens only at high cell densities.

Model0flow WUR.png

Fig.3: Graphs showing oscillatory behavior increasing with cell density even when the flow rate is set to 0

Disregarding the strictly mathematical model, some more conclusions were drawn from modeling work as to whether to use a microdish or microsieve as platform for measuring oscillations. This is described in more detail in the section about the setup of the device.


Expansion of the model for the double tunable oscillatory construct

To model the double tunable oscillator, the original set of delay differential equations [2] was expanded. This involved adding the effect of the repressor and activator in the form of a Hill function to the differential equations for LuxI and AiiA. The equations for the internal and external AHL concentrations were not changed. The new set of equations is seen below.

Equations enhanced WUR.png

Top: Adjusted set of differential equations for the double tunable oscillator

Since the lacI repressed hybrid promoter used for this system was designed by the Tokyo iGEM 2007 team [3], the additional terms were derived according to their model. The tetR repressor was then modeled using the same template. The set of equations seen above therefore integrates the features of both models.

This set of equations was then used to write a new matlab modeling tool, of which the GUI is depicted below.


GUI 2 synchroscillator WUR.png

As can be seen in the GUI, the new model allows the construct to be controlled by changing the values for IPTG and aTc. The parameters used are still in arbitrary units and the values not yet verified with wet lab experiments. However, some example conditions were simulated, of which the output graphs are seen below.


back to top

Links and references:

[1][http://www.nature.com/nature/journal/v463/n7279/abs/nature08753.html Danino et al. 2010]

[2][http://www.nature.com/nature/journal/v463/n7279/suppinfo/nature08753.html Supplementary information]

[3][http://parts.mit.edu/igem07/index.php/Tokyo/Works Tokyo iGEM 2007]