Team:HUST-China/Modelling

From 2013.igem.org

The genetic pathway we used in modelling work is described as below.

Fig 1.The pathway of genetic oscillator used in our project.

Approach and Result

Our modeling work includes 3 parts: (1)Feasibility, stability and sensitivity of single genetic oscillator and its parameters;
(2)Feasibility of multi cells’ oscillator made up by a group of single genetic oscillators and comparison with experiment data;
(3)Statistically analyzing our wet-lab result.

To finish part (1), we construct DDEs model, which superb model using ODEs because it is capable of describing time interval between expression of regulating protein and their binding to promoters; in other words, our model is much more precise than ODEs model. Moreover, DDEs model is more flexible than ODEs model, as we can degrade it into an ODEs one by setting lag equal to 0. Result of our model is satisfying, as our oscillator is stable while circumstance(Arabinose and IPTG) varies, and sensitive when itself($\tau$) changes. That is to say, it is easy for manipulate, which is good for further application.

To accomplish part (2), we use data from part (1) and simulate using Agent-Based Model. We then compare simulation outcome with our wet-lab data. There comes a promising result, as florescent intensity linearly related to expression of AraC of these two groups of data match each other well, namely our model is supported by our experiment. This model is extendable because it has capability simulating other problem in large population. Also method that used to compar wetlab and drylab data, which is Q-Q plot, provides a quick way for analyzing groups of data correspondent time interval remained unknown.

To fulfill part (3), we use Q-Q plot and K-S test to confirm that single transformed cells and co-transformed cells are different. We then use S-W test to demonstrate that the oscillatory behavior of our oscillator is not a Brownian motion, which is often neglected regardless of its importance.

Notes

Throughout the whole modelling work, we choose to study AraC instead of mRFP since they are in the same plasmid and we assumed that the expression rate of both protein is similar. By doing this, we can reduce the number of equations.