Difference between revisions of "Team:UFAM-UEA Brazil/MathModeling"

Line 72: Line 72:
 
https://static.igem.org/mediawiki/2016/7/76/UFAM_UEA_Team_math_%2813%29.png" width=80 height=40/> then <img src="
 
https://static.igem.org/mediawiki/2016/7/76/UFAM_UEA_Team_math_%2813%29.png" width=80 height=40/> then <img src="
 
https://static.igem.org/mediawiki/2016/c/c5/UFAM_UEA_Team_math_%2814%29.png" width=150 height=60/> so we have that the variation of the interior mercury is</p>
 
https://static.igem.org/mediawiki/2016/c/c5/UFAM_UEA_Team_math_%2814%29.png" width=150 height=60/> so we have that the variation of the interior mercury is</p>
 +
 +
<p>(7)<img src="
 +
https://static.igem.org/mediawiki/2016/5/53/UFAM_UEA_Team_math_%2815%29.png" width=350 height=80/></p>
 +
<p>And assume that mRNA is at Quasi-Steady state, ie. <img src="
 +
https://static.igem.org/mediawiki/2016/5/52/UFAM_UEA_Team_math_%2816%29.png" width=80 height=40/>, then we have
 +
that the amount of mRNA is
 +
</p>
 +
 +
<p>(8)<img src="https://static.igem.org/mediawiki/2016/4/4d/UFAM_UEA_Team_math_%2817%29.png" width=350 height=80/></p>
 +
 +
<p>Then we have that the velocity of production of RFP is</p>
 +
 +
<p>(9)<img src="https://static.igem.org/mediawiki/2016/0/0f/UFAM_UEA_Team_math_%2818%29.png" width=350 height=80/></p>
 +
 +
<p>So we will also assume that the degradation of RFP is very slow, δ2 = 0, and also the stochastic production of mRNA is very slow, α = 0. Then finaly we have that the production of RFP is</p>
 +
<p>(10)<img src="
 +
https://static.igem.org/mediawiki/2016/c/c7/UFAM_UEA_Team_math_%2819%29.png" width=350 height=80/></p>
 +
 +
<p>Where <img src="
 +
https://static.igem.org/mediawiki/2016/e/e2/UFAM_UEA_Team_math_%2820%29.png" width=80 height=40/></p>
 +
  
  

Revision as of 20:41, 18 October 2016


MODEL OF RFP PRODUCTION

1. Modelling

We want to characterize the promoter’s velocity of expression due to presence of mercury, so we will attach an RFP gene to it to produce.

First we will model the production of RFP due to Hg2+ at stationary time phase with this configuration. Then we add the initial phase consideration

1.1 Stationary time phase.

The GMO we have is:

Where MerR is a repressor and releases when H2+ is presented in the interior of the cell, we will note this amount as Hgin and the exterior as Hgout. And Promoter nP is the promoter which velocity we want to characterize. The activation of the gene in the pressence of Hgin is represented like:

where G is the inactive form of gene repressed by merR, and X is its active form. Through the law of mass action we derive the diferential equation

(1)

At the equillibrium state, ie. , we have that the proportion of genes in the activated gene is

(2)

Where

So this is the average production rate of a typical gene, so the average mRNA production will be

(3)

Where α is the production of mRNA due to the stochastic nature of the binding of merR for repressing the RNA Polymerase and δ1 is the degradation factor of mRNA. The reaction for production of RFP from mRNA is

so we derive the diferential equation

(4)

where δ2 is the degradation constant of [RFP]. The variation of exterior mercury is

from here we derive the differential equation

(5)

As the interior mercury is used by the inactive form of the gene, the variation of Hgin2+ is

(6)

To solve this system we will assume that the permeability of the cell memebrane to mercury is instantaneus then we have that then so we have that the variation of the interior mercury is

(7)

And assume that mRNA is at Quasi-Steady state, ie. , then we have that the amount of mRNA is

(8)

Then we have that the velocity of production of RFP is

(9)

So we will also assume that the degradation of RFP is very slow, δ2 = 0, and also the stochastic production of mRNA is very slow, α = 0. Then finaly we have that the production of RFP is

(10)

Where