Team:NYMU-Taipei/Modeling/Linear epidemic model

From 2013.igem.org

Revision as of 23:56, 27 September 2013 by Mastershot (Talk | contribs)

National Yang Ming University


Contents

Ethanol model

Background:

trxC promoter is an oxyR-activated promoter. Behind trxC promoter are several terminals to inhibit T7 polymerase-producing gene from being opened easily ; namely, it will not open until the concentration of oxyR is high enough, which means the bees are facing disastrous infection. Once oxyR concentration overcomes the threshold and conquers the terminal obstacles, T7 polymerase is produced and will bind to T7 promoter, which is a specific promoter binding only to T7 polymerase. Behind the T7promoter are enzyme PDC and ADH-producing genes, which will convert pyruvate to ethanol, and thus kill spores of Nosema Ceranae.

Pathway of ethanol:

PDC= pyruvate decarboxylase; ADH= Acetaldehyde

The purpose of this modeling is:


1.To simulate how many terminal do we need as a threshold to have the T7 polymerase-producing gene open at a proper time.
2.To simulate the threshold concentration of oxyR to conquer the terminal.
3.To determine the time of opening to see if the circuit could open in time (useful).


We assume that the equilibrium tendency of pyruvate to acetaldehyde is overwhelmingly rightward, while acetaldehyde to ethanol is bidirectional. We also assume that mRNA degrading rate is neglectable. As for how we measure promoter strengths, we choose PoPS, which is the rate of RNA polymerase binding to the DNA and trigger the transcription of the certain gene. It is assumed that the dose of pyruvate is sufficient in bees’ body through enough proliferation/copy number of Beecoli; the possibility of RNApolymerase skipping terminators is assumed to be proportional to time span; equilibrium between pyruvate and acetaldehyde is dominantly rightward, while equilibrium between acetaldehyde and ethanol is bidirectional; concentration of ethanol will not easily decrease and will sustain for a period of time.

Equation1:


1.KdROSoxyR = dissociation constant of ROSoxyR
2.n ROSoxyR = Hill coefficient of ROSoxyR
3.PoPStrxC = promoter strength of trxC
4.KdT7 = dissociation constant of T7
5.NT7 = Hill coefficient of T7
6.PoPST7 = promoter strength of T7 promoter
7.kdegmRNA = degrading constant of sensor promoter mRNA
8.N = number of plasmid in a single cell
9.V = volume of a cell


The aim of the equation is to know mRNA of T7 polymerase production rate and when it can reach the level to translate enough T7 poltmerase.


Equation2:

\frac{d[protein of T7 polymerase]}{dt}={RBS}\times{PoPStrxC}\times{\frac {[ROSoxyR]^{n1} }{[Kd]}^{n1}+[ROSoxyR]^{n1}}\times{\frac{N}{V}}\times{a}+PoPST7\times{T7effect}\times{\frac{N}{V}}-KdegT7\times{T7}

1.PoPStrxC = promotor strength of trxC promoter
2.PoPST7 = promoter strength of T7 promoter
3.RBS = binding site strength
4.kdegT7 = degrading constant of T7 polymerase
5.a = posibility of conquering the threshold concentration
6.n1 = Hill coefficient of ROSoxyR (complex of ROS+oxyR)
7.N = number of plasmid in a single cell
8.V = volume of a cell

The aim of the equation is to know the threshold concentration of oxyR to conquer the terminal and when T7 polymerase can reach the required concentration to activate T7 promoter.


Equation3:

\frac{d[mRNAPDC]}{dt}=\frac{ [T7]^{nT7} }{ KdT7^{nT7}+[T7]^{nY7} }\timesPoPST7\times\frac{N}{V}-kdegmRNA\times[mRNAPDC]

1.KdT7 = dissociation constant of T7
2.NT7 = Hill coefficient of T7
3.PoPST7 = promoter strength of T7 promoter
4.kdegmRNA = degrading constant of sensor promoter mRNA
5.N = number of plasmid in a single cell
6.V = volume of a cell


The aim of the equation is to know mRNA of enzyme PDC production rate and when it can reach the level to translate enough PDC.