Team:Tianjin/Modeling

From 2013.igem.org

(Difference between revisions)
(Analysis & Discussion)
 
(35 intermediate revisions not shown)
Line 1: Line 1:
<html>
<html>
-
</div>
+
 
-
</div>
+
</div>
</div>
<style type="text/css">
<style type="text/css">
-
*{margin:0;padding:0;font-family:"微软雅黑","Arial";}
+
*{margin:0;padding:0;font-family:"微软雅黑","Times New Roman";}
body{
body{
Line 17: Line 16:
}
}
-
h1,h2,h3,h4,h5,h6{
+
h2,h3,h4,h5,h6{
border:none;}
border:none;}
p {margin:0.5em 0 !important;
p {margin:0.5em 0 !important;
text-decoration:none;
text-decoration:none;
-
font-family:Arial, Helvetica, sans-serif;
+
font-family:Times New Roman, Helvetica, sans-serif;
font-size:12px;
font-size:12px;
}
}
Line 36: Line 35:
float:left;
float:left;
margin-top:5px;
margin-top:5px;
-
background:url(https://static.igem.org/mediawiki/2013/d/d8/Modeling_banner.jpg) no-repeat;
+
background:url(https://static.igem.org/mediawiki/igem.org/1/16/Tju-Project-banner.jpg) no-repeat;
}
}
.note{
.note{
Line 62: Line 61:
padding:30px 10px 10px 30px;
padding:30px 10px 10px 30px;
}
}
 +
#footer{
#footer{
margin:0 auto;
margin:0 auto;
Line 70: Line 70:
float:left;
float:left;
}
}
-
#goTopBtn {POSITION: fixed;TEXT-ALIGN: center;LINE-HEIGHT: 30px;WIDTH: 100px;BOTTOM: 35px;HEIGHT: 100px;FONT-SIZE: 12px;RIGHT: 30px;}
+
 
 +
#goTopBtn  
 +
{POSITION: fixed;TEXT-ALIGN: center;LINE-HEIGHT: 30px;WIDTH: 100px;BOTTOM: 35px;HEIGHT: 100px;FONT-SIZE: 12px;RIGHT: 30px;}
 +
 
 +
#globalWrapper{
 +
width:960px;
 +
height:auto;
 +
background-color:#fff;
 +
}
Line 151: Line 159:
.main img{border:hidden;margin-bottom:5px;}
.main img{border:hidden;margin-bottom:5px;}
.main div,li,p{font-family:Arial;line-height:150%;word-spacing:0.1em;}
.main div,li,p{font-family:Arial;line-height:150%;word-spacing:0.1em;}
-
 
+
.main li{list-style:disc;}
.box p{color:#000;font-family:Arial, Helvetica, sans-serif;font-size:14px;line-height:150%;text-align:left; clear:both;}
.box p{color:#000;font-family:Arial, Helvetica, sans-serif;font-size:14px;line-height:150%;text-align:left; clear:both;}
Line 169: Line 177:
$(window).scroll(function(e){
$(window).scroll(function(e){
s = $(document).scrollTop();
s = $(document).scrollTop();
-
if(s > t - 100){
+
if(s > t - 70){
$('.fixed').css('position','fixed');
$('.fixed').css('position','fixed');
if(s + fh - 2000 > mh){
if(s + fh - 2000 > mh){
Line 215: Line 223:
         <ul>  
         <ul>  
             <li class="hmain" style="margin-top:20px;">  
             <li class="hmain" style="margin-top:20px;">  
-
                    <a href="#anchor-m-01"> Model Objective</a>  
+
                <a href="#Model_Objective">Model Objective</a>  
             </li>  
             </li>  
             <li class="hmain">  
             <li class="hmain">  
-
                 <a href="#anchor-m-02">Problem Description</a>
+
                 <a href="#Problem_Description">Problem Description</a>
             </li>  
             </li>  
             <li class="hmain">  
             <li class="hmain">  
-
                 <a href="#anchor-m-03">Problem Abstraction
+
                 <a href="#Problem_Abstraction">Problem Abstraction</a>  
-
</a>  
+
             </li>  
             </li>  
             <li class="hmain">  
             <li class="hmain">  
-
                 <a href="#anchor-m-04">Assumption</a>  
+
                 <a href="#Assumption">Assumption</a>  
             </li>
             </li>
             <li class="hmain">  
             <li class="hmain">  
-
                 <a href="#anchor-m-05">Model Development</a>  
+
                 <a href="#Model_Development">Model Development</a>  
             </li>
             </li>
             <li class="hmain">  
             <li class="hmain">  
-
                 <a href="#anchor-m-06">Analysis & Discussion</a>  
+
                 <a href="#Analysis_.26_Discussion">Analysis & Discussion</a>  
             </li>
             </li>
             <li class="hmain">  
             <li class="hmain">  
-
                 <a href="#anchor-m-07">Achievement</a>  
+
                 <a href="#Achievement">Achievement</a>  
             </li>
             </li>
Line 250: Line 257:
<div class="main">
<div class="main">
-
<center><span style="font-family:Arial;font-size:46px;color:#000;"> Mathematic Analysis on AlkSensor</span></center>
 
-
 
-
<a name="anchor-m-01" id="anchor-m-01"></a>
 
<br/>
<br/>
-
 
+
<center><span style="font-family:Arial;font-size:45px;color:#000;"> Mathematic Analysis on AlkSensor </span></center>
 +
<br/>
</html>
</html>
 +
= Model Objective =
 +
<br>
 +
To perform a mathematic analysis on AlkSensor, to find out the relationship between AlkSensor’s input and output, to find ways to specifically regulate the function of AlkSensor, to find the factors affecting sensor’s leakage.
-
= Model Objective=
+
<br>
-
<html>
+
=Problem Description=
-
<hr /><br />
+
<br>
-
<p>We perform this mathematical analysis on AlkSensor</p>
+
AlkSensor is composed of protein ALKR and promoter alkM. Genes of ALKR and alkM were synthesized and constructed into a plasmid, as shown in figure 1. As mentioned in the introduction, protein ALKR is a transcription factor which can recognize certain alkanes. ALKR is under a constitutive promoter and is constitutively expressed. Alkane molecules is recognized by ALKR and a dimerized ALKR-alkane complex is formed. Subsequently the reporter’s promoter alkM is induced by the complex and the genes in the downstream of alkM are expressed. However, promoter alkM could also be slightly induced when inducers are absent from the system. There are possibilities that dierized ALKR binds on promoter alkM and activates the downstream genes.
-
<ul class="modeling-list">
+
[[Image:Modeling-01.png|thumb|600px|center|<b>Figure 1.</b> Scheme of AlkSensor’s mechanism]]
-
<li> to perform a mathematic analysis on AlkSensor,</li>
+
-
<li> to find out the relationship between AlkSensor’s input and output,</li>
+
-
<li> to find ways to specifically regulate the function of AlkSensor.</li>
+
-
</ul>
+
 +
<br>
-
<a name="anchor-m-02" id="anchor-m-02"></a>
+
=Problem Abstraction=
-
<br/>
+
<br>
 +
The mechanism of AlkSensor can be represented by a set of chemical reactions, as shown below.
-
</html>
+
[[Image:TJU-F-11.png|center]]
-
= Problem Description=
+
<center><b>Figure 2.</b> Chemical reactions of AlkSensor</center>
-
<html>
+
*Reaction 1: the generation of protein ALKR.
 +
*Reaction 2: the degradation of protein ALKR
-
<hr /><br />
+
*Reaction 3: the dimerization of protein ALKR(AR2-A denotes the dimerized ALKR-alkane complex)
-
<p> AlkSensor is composed of protein ALKR and promoter alkM. Genes of ALKR and alkM were synthesized and constructed into a plasmid, as shown in figure 1. As mentioned in the introduction, protein ALKR is a transcription factor which can recognize certain alkanes. ALKR is under a constitutive promoter and is constitutively expressed. Alkane molecules is recognized by ALKR and a dimerized ALKR-alkane complex is formed. Subsequently the reporter’s promoter alkM is induced by the complex. </p>
+
*Reaction 4: the combination of dimerized ALKR with alkane(Pm denotes promoter alkM)
-
<div style="margin:0 50px 15px 50px;padding:5px 5px 5px 5px;background-color:#fafafa;border:thin solid #999;width:600px;">
+
-
<div>
+
-
<a href="https://static.igem.org/mediawiki/2013/e/ee/Modeling-01.png" target="_blank" ><img src="https://static.igem.org/mediawiki/2013/e/ee/Modeling-01.png" width="600px" height="338px" border="thin solid #999";/></a>
+
-
<div id="Enlarge">
+
-
<a href="https://static.igem.org/mediawiki/2013/e/ee/Modeling-01.png"  title="Enlarge" target="_blank"><img src=" https://static.igem.org/mediawiki/2013/9/90/Enlarge.jpg " width="20" height="20" alt="" align="right"/></a>
+
-
<b>Figure 1.</b> Scheme of AlkSensor’s mechanism </div>
+
-
</div></div>
+
-
<a name="anchor-m-03" id="anchor-m-03"></a>
+
*Reaction 5: the combination of dimerized ALKR-alkane complex with promoter alkM(Pm’ denotes the promoter alkM binding with inducer)
-
<br/>
+
-
</html>
+
-
= Problem Abstraction=
+
*Reaction 6: the combination of dimerized ALKR with promoter alkM(Pm’’ denotes the promoter alkM binding with dimerized ALKR)
-
<html>
+
*Reaction 7: the generation of reporter (RNAP denote the RNA polymerase, Rp denotes reporter )
-
<hr /><br />
+
*Reaction 8: the leakage of AlkSensor
-
<p> The mechanism of AlkSensor can be represented by a set of chemical reactions, as shown below.</p>
+
*Reaction 9: the the degradation of reporter
-
<div style="text-align:center;vertical-align:middle;"><a href="https://2013.igem.org/File:Formula-01.gif" target="_blank" ><img src="https://static.igem.org/mediawiki/2013/f/ff/Formula-01.gif"/></a></div>
+
<br>
-
<style type="text/css">
+
= Assumption =
-
.modeling-list{font-size:14px;line-height:150%;}
+
<br>
-
.modeling-list li{list-style-type:disc}
+
-
</style>
+
-
<ul class="modeling-list">
+
There are several assumptions in the abstraction of AlkSensor.
-
<li> Reaction 1: the generation of protein ALKR.</li>
+
-
<li> Reaction 2: the degradation of protein ALKR</li>
+
-
<li> Reaction 3: the dimerization of protein ALKR(AR2-A denotes the dimerized ALKR-alkane complex)</li>
+
-
<li> Reaction 4: the combination of dimerized ALKR with alkane(Pm denotes promoter alkM)</li>
+
-
<li>Reaction 5: the combination of dimerized ALKR-alkane complex with promoter alkM(Pm’ denotes the promoter alkM binding with inducer)</li>
+
-
<li> Reaction 6: the generation of reporter(RNAP denote the RNA polymerase, Rp denotes reporter )</li>
+
-
<li> Reaction 7: the the degradation of reporter</li>
+
-
</ul>
+
 +
*Because the genes of ALKR is under a constitutive promoter, we assume that the generation rate of ALKR is constant.
-
<a name="anchor-m-04" id="anchor-m-04"></a>
+
*The four binding reactions, i.e, reaction 3, 4, 5 and 6 are fast reversible reactions since DNA binding and unbinding of the repressors dimers and the dimerization itself occur within seconds, whereas the synthesis (transcription, translation, folding) and degradation of monomers takes minutes to sometimes an hour.
-
<br/>
+
-
</html>
+
-
= Assumption=
+
Reaction 3, 4, 5 and 6 can be integrated into two reversible questions. Therefore, the 9 reactions can be simplified as 7 reactions.
-
<html>
+
[[Image:TJU-F-12.png|center|]]
-
<hr /><br />
+
-
<p> There are several assumptions in the abstraction of AlkSensor.</p>
+
-
<ul class="modeling-list">
+
<center><b>Figure 3.</b> Simplified form of AlkSensor’s chemical reactions</center>
-
<li> Because the the genes of ALKR is under a constitutive promoter, we assume that the generation rate of ALKR is constant.</li>
+
-
<li> The three binding reactions, i.e, reaction 3, 4 and 5 are fast reversible reactions since DNA binding and unbinding of the repressors dimers and the dimerization itself occur within seconds, whereas the synthesis (transcription, translation, folding) and degradation of monomers takes minutes to sometimes an hour. </li>
+
-
<li> For simplification, we do not take into consideration the combination of dimerization ALKR with promoter alkM.</li>
+
-
</ul>
+
 +
<br>
-
<a name="anchor-m-05" id="anchor-m-05"></a>
+
=Model Development=
-
<br/>
+
<br>
-
</html>
+
-
= Model Development=
+
Reaction 3 and 4 are fast reversible reactions because DNA binding and unbinding of the transcription factor dimers and the dimerization itself occur within seconds. However the synthesis and degradation of monomers takes minutes to sometimes an hour. Therefore, the three reactions are in equilibrium and the steady state concentrations are given in terms of the equilibrium constants. The accumulation rates of ALKR and reporter are subject to their generation rate and degradation rate, which is described by two ordinary differential equations (ODE). So the mechanism of AlkSensor can be described by a mathematical model consisting of 2 ordinary differential equations (ODE) and three equilibrium equation. The input of AlkSensor can be defined as the concentration alkane, the output can be defined as the concentration of reporter protein, as shown in figure 2.
-
<html>
+
[[Image:TJU-model-4.png|thumb|600px|center|<b>Figure 4.</b> Mathematical representation of AlkSensor]]
-
<hr /><br />
+
In the first equation, there are four terms that determine the accumulation rate of reporter, the generation rate of reporter, the leakage rate of AlkSensor, the basal expression rate of reporter, the degradation rate of reporter. In the second equation, there are four terms that determine the accumulation rate of ALKR, the generation rate of ALKR, the degradation rate of ALKR, the dilution rate of ALKR. The last 3 equation describe the quantitative relation in reaction 3 and 4.
-
<p> In the set of reactions reaction 3, 4, and 5 are fast reversible reactions because DNA binding and unbinding of the transcription factor dimers and the dimerization itself occur within seconds. However the synthesis and degradation of monomers takes minutes to sometimes an hour. Therefore, the three reactions are in equilibrium and the steady state concentrations are given in terms of the equilibrium constants. The accumulation rates of ALKR and reporter are subject to their generation rate and degradation rate, which is described by two ordinary differential equations(ODE). So the mechanism of AlkSensor can be described by a mathematical model consisting of 2 ordinary differential equations(ODE) and three equilibrium equation. The input of AlkSensor can be defined as the concentration alkane, the output can be defined as the concentration of reporter protein, as shown in figure 2.
+
The [Pm’] and [Pm’’] are still intermediate variables, we next deduce the value of [Pm’] and [Pm’’].
-
</p>
+
-
<div style="padding:5px 5px 5px 5px;background-color:#fafafa;border:thin solid #999;width:700px;">
+
[[Image:TJU-F-1.png|center]]
-
<div>
+
-
<a href="https://static.igem.org/mediawiki/2013/b/b7/Modeling-02.png" target="_blank" ><img src="https://static.igem.org/mediawiki/2013/b/b7/Modeling-02.png" width="700px" height="310px";/></a>
+
-
<div id="Enlarge">
+
-
<a href="https://static.igem.org/mediawiki/2013/b/b7/Modeling-02.png"  title="Enlarge" target="_blank"><img src="https://static.igem.org/mediawiki/2013/9/90/Enlarge.jpg" width="20" height="20" alt="" align="right"/></a>
+
-
<b>Figure 2.</b> Mathematical representation of AlkSensor </div>
+
-
</div></div>
+
-
<p> In the first equation, there are three terms that determine the accumulation rate of ALKR, the generation rate of ALKR, the degradation rate of ALKR, the dilution rate of ALKR.
+
[[Image:TJU-F-2.png|center]]
-
In the second equation, there are four terms that determine the accumulation rate of reporter, the generation rate of reporter, the basal expression rate of reporter, the degradation rate of reporter, the dilution rate of reporter. The last 3 equation describe the quantitative relation in reaction 3,4 and 5.</p>
+
-
<p> The [Pm’] is still a intermediate variable, we next deduce the value of [Pm’].</p>
+
The total concentration of promoter alkM is thought to be a constant, which is proportional to the copy number of alkM.
-
<div style="text-align:center;vertical-align:middle;"><a href="https://2013.igem.org/File:Formula-02-1.gif" target="_blank" ><img src="https://static.igem.org/mediawiki/2013/7/78/Formula-02-1.gif" /></a></div>
+
[[Image:TJU-F-3.png|center]]
-
<p> The total concentration of promoter alkM is though to be a constant, which is proportional to the copy number of alkM.</p>
+
So we can get the expression of [Pm], [Pm’] and [Pm’’]
-
<div style="text-align:center;vertical-align:middle;"><a href="https://2013.igem.org/File:Formula-03.gif" target="_blank" ><img src="https://static.igem.org/mediawiki/2013/a/a8/Formula-03.gif" /></a></div>
+
[[Image:TJU-F-4.png|center]]
-
<p> So we can get the the expression of [Pm] and [Pm’]</p>
+
[[Image:TJU-F-5.png|center]]
-
<div style="text-align:center;vertical-align:middle;"><a href="https://2013.igem.org/File:Formula-04.gif" target="_blank" ><img src="https://static.igem.org/mediawiki/2013/f/f2/Formula-04.gif" /></a></div>
+
[[Image:TJU-F-6.png|center]]
-
<p> Then we substitute [Pm’] with the terms on the right side of the equation.</p>
+
Then we substitute [Pm’] and [Pm’] with the terms on the right side of the equation.
-
<div style="text-align:center;vertical-align:middle;"><a href="https://2013.igem.org/File:Formula-05.gif" target="_blank" ><img src="https://static.igem.org/mediawiki/2013/3/37/Formula-05.gif" /></a></div>
+
[[Image:TJU-F-7.png|center]]
-
<p> For simplification, we substitute K<sub>4</sub>K<sub>3</sub>K<sub>2</sub> with K<sub>ap</sub>. K<sub>ap</sub>, equals to K<sub>4</sub>K<sub>3</sub>K<sub>2</sub>, is the apparent equilibrium constant of reaction 8, .</p>
+
When the leakage of ALKR’s promoter is little, the terms r can be ignored. When the system reach a steady state, the accumulation rate of Rp is zero. In such a condition, the response function of AlkSensor can be described as the equation below.
-
<div style="text-align:center;vertical-align:middle;"><a href="https://2013.igem.org/File:Formula-06.gif" target="_blank" ><img src="https://static.igem.org/mediawiki/2013/7/71/Formula-06.gif" /></a></div>
+
[[Image:TJU-F-8.png|center]]
-
<p> Finally we get the ODE about the reporter.</p>
+
We can simplify the response function into an elegant form, in which x denotes the input, y denotes the output, a, b and c are three parameters.
-
<div style="text-align:center;vertical-align:middle;"><a href="https://2013.igem.org/File:Formula-07.gif" target="_blank" ><img src="https://static.igem.org/mediawiki/2013/a/a1/Formula-07.gif" /></a></div>
+
[[Image:TJU-F-9.png|center]]
-
<p> When the leakage of ALKR’s promoter is little, the terms r can be ignored. When the system reach a steady state, the accumulation rate of Rp is zero. </p>
+
*k’: the expression rate constant of promoter alkM;
 +
*k<sub>d</sub>: the degradation constant of protein ALKR, thought to be a constant;
 +
*P<sub>0</sub>: the concentration of RNAP, constant;
 +
*[P<sub>t</sub>]: the total concentration of promoter alkM, which is proportional to the copy number. This is easily changeable. We can change it through construct alkM on different plasmid with varying copy number.
-
<div style="text-align:center;vertical-align:middle;"><a href="https://2013.igem.org/File:Formula-08.gif" target="_blank" ><img src="https://static.igem.org/mediawiki/2013/d/d0/Formula-08.gif" /></a></div>
+
<br>
-
 
+
-
<p> It is easy to get the equation that describe the relation between AlkSensor’s input and output.</p>
+
-
 
+
-
<div style="text-align:center;vertical-align:middle;"><a href="https://2013.igem.org/File:Formula-09.gif" target="_blank" ><img src="https://static.igem.org/mediawiki/2013/e/eb/Formula-09.gif" /></a></div>
+
-
 
+
-
<p> The equation can be simplified into a very elegant form. Let x represent [Alkane], y represent [Reporter], the equation can be simplified as <a href="https://2013.igem.org/File:Formula-10.gif" target="_blank" ><img src="https://static.igem.org/mediawiki/2013/c/cd/Formula-10.gif" /></a>, in which <a href="https://2013.igem.org/File:Formula-10-2.gif" target="_blank" ><img src="https://static.igem.org/mediawiki/2013/4/40/Formula-10-2.gif"/></a>,<a href="https://2013.igem.org/File:Formula-11.gif" target="_blank" ><img src="https://static.igem.org/mediawiki/2013/4/41/Formula-11.gif" /></a>. </p>
+
-
 
+
-
 
+
-
<a name="anchor-m-06" id="anchor-m-06"></a>
+
-
<br/>
+
-
</html>
+
= Analysis & Discussion =
= Analysis & Discussion =
 +
<br>
-
<html>
+
The output has a monotone increasing relation with the input. When the input is large enough, the output reaches its maximum value, a, which is proportional to the amount of promoter alkM. When there is no input, the output has a smallest value, i.e. the value of leakage. We can easily obtain the value of dynamic range by dividing the maximum response into the value of leakage.
-
<hr /><br />
+
-
<p><b> What is the response curve of AlkSensor like?</b></p>
+
-
<p> Figure 3 shows twelve different response curves of AlkSensor. Although they are in different shapes, they are in the same pattern——<a href="https://2013.igem.org/File:Formula-12.gif" target="_blank" ><img src="https://static.igem.org/mediawiki/2013/f/fc/Formula-12.gif"/></a>.</p>
+
-
<p> When the input is very low, the response curve is nearly liner and has a maximum value of slope, aK. With the increase of input, the value of slope decrease. Finally the response curve approach a horizontal line y=a.</p>
+
-
<br />
+
[[Image:TJU-F-10.png|center]]
-
<p><b> How to evaluate the response curve of AlkSensor?</b></p>
+
-
<p> The response curve describes how the output is determined by the input. The slope of the curve represent output’s sensitivity to input. The bigger the slope is, the more sensitive output is to input. If the slope is zero, the output does’t change with input. This is when the AlkSensor loses its function.</p>
+
-
<p> We can derive two critical parameters from the response curve to evaluate AlkSensor, <b>sensitive range of input(SRI)</b> and <b>response range of output(RRO)</b>. SRI is the range in which the output is sensitive enough to the input. RRO is the range of AlkSensor’s signal strength. The crossing zone of AlkSensor’s SRI and RRO is the effective working zone of AlkSensor.</p>
+
-
<br />
+
There are two parameters that can be rationally regulated, [P<sub>t</sub>] and [ALKR]. We can control [P<sub>t</sub>] by adjust the copy number of promoter alkM and control [ALKR] by adjusting the strength of its promoter. Therefore, we have the ability to regulate the value of a and b, that is to say, the maximum response, leakage and dynamic range of AlkSensor can all be rationally regulated.
-
<p><b> How does <i>a</i> and <i>K</i> influence the shape of the response curve?</b></p>
+
-
<p> We list a set of response of curves with different values of <b>a</b> and <b>K</b>, as shown in figure 3. It is easy to get that the RRO is <b>a</b>. The SRI is subject to <b>K</b>, the lower <b>K</b> is, the wider the SRI is.</p>
+
-
<div style="padding:5px 5px 5px 5px;background-color:#fafafa;border:thin solid #999;width:700px;">
+
Figure 5 shows a set of response curves of AlkSensor with different expression levels of ALKR. When there is no ALKR expressed, the response curve is a horizontal line, suggesting that AlkSensor has no ability to distinguish different inputs. With the increase of [ALKR], the leakage increase, however, the maximum response remains the same value, which makes the dynamic range decrease. When [ALKR] is too large, the response curve turns to be a horizontal line again. AlkSensor loses its function in this case.
-
<div>
+
-
<a href="https://static.igem.org/mediawiki/2013/8/85/Modeling-03.png" target="_blank" ><img src="https://static.igem.org/mediawiki/2013/8/85/Modeling-03.png" width="700px" height="400px";/></a>
+
-
<div id="Enlarge">
+
-
<a href="https://static.igem.org/mediawiki/2013/8/85/Modeling-03.png"  title="Enlarge" target="_blank"><img src=" https://static.igem.org/mediawiki/2013/9/90/Enlarge.jpg " width="20" height="20" alt="" align="right"/></a>
+
-
<b>Figure 3.</b> AlkSensor’s response curve with different values of a and K </div>
+
-
</div></div>
+
 +
[[Image:TJU-model-5.png|thumb|600px|center|<b>Figure 5.</b> The influence of [ALKR] on the response curve of AlkSensor]]
-
<br />
+
<br>
-
<p><b> How to regulate the value of <i>a</i> and <i>K</i> of AlkSensor?</b></p>
+
-
<p> To achieve the goal of rationally regulating AlkSensor, we need to know how to adjust the value of <b>a</b> and <b>K</b>. We already have the equation of <b>a</b> and <b>K</b>.</p>
+
-
 
+
-
<div style="text-align:center;vertical-align:middle;"><a href="https://2013.igem.org/File:Formula-13.gif" target="_blank" ><img src="https://static.igem.org/mediawiki/2013/9/9f/Formula-13.gif" align="middle"/></a></div>
+
-
 
+
-
<p> The value of <b>a</b> is determined by five values. They are:</p>
+
-
<ul class="modeling-list">
+
-
<li> k<sub>t</sub>: the combination constant of RNAP with Promoter alkM. It is hardly changeable, so we can consider it to be a constant;</li>
+
-
<li> k<sub>d</sub>: the degradation constant of protein ALKR, also thought to be a constant;</li>
+
-
<li> P<sub>0</sub>: the concentration of RNAP, constant;</li>
+
-
<li>[P<sub>t</sub>]: the total concentration of  promoter alkM, which is proportional to the copy number. This is easily changeable. We can change it through construct alkM on different plasmid with varying copy number.</li>
+
-
<li>μ:the dilution constant caused by the cell cycle, which is hardly changeable.</li>
+
-
</ul>
+
-
<p> As mention before, the <b>a</b> represent the response range of the sensor’s out put, so through changing the copy number, we can rationally control AlkSensor’s RRO.</p>
+
-
<p> Next comes the <b>K</b>, which influence the sensitive range of the sensor’s input. <b>K</b> equals the product of K<sub>ap</sub>, and the square of ALKR’s concentration. K<sub>ap</sub> is the product of three reactions’ equilibrium constant: the dimerization of ALKR, the formation of inducer, the activation of promoter alkM. The structure of ALKR certainly influence the value of <b>K</b>. The value of <b>K</b> also reflect AlkSensor’s specificity towards input. Through altering the structure of ALKR, we can change AlkSensor specificity towards certain alkane. ALKR’s concentration is determined by the generation rate of ALKR, which can easily be regulated by changing the promoter of ALKR. As mentioned before, ALKR is under a constitutive promoter, the strength and copy number of the promoter have a positive relationship with ALKR’s concentration. So, through changing the strength and copy number of ALKR’s promoter, we can rationally regulate the sensitive range of AlkSensor.</p>
+
-
 
+
-
 
+
-
 
+
-
<a name="anchor-m-07" id="anchor-m-07"></a>
+
-
<br/>
+
-
</html>
+
= Achievement =
= Achievement =
<html>
<html>
-
<hr /><br />
+
<br />
<p> In this model, we have</p>
<p> In this model, we have</p>
Line 467: Line 398:
<li style="list-style-image: url(https://static.igem.org/mediawiki/2013/e/ed/Tick.gif)"> developed strategies to rationally regulate AlkSensor.</li>
<li style="list-style-image: url(https://static.igem.org/mediawiki/2013/e/ed/Tick.gif)"> developed strategies to rationally regulate AlkSensor.</li>
</ul>
</ul>
 +
Line 505: Line 437:
<!--top button section-->
<!--top button section-->
<div style="display: block" id="goTopBtn">
<div style="display: block" id="goTopBtn">
-
<a href="#top" title="Top"><img border=0 src="https://static.igem.org/mediawiki/2013/6/61/Top.jpg"></a>
+
<a href="#top" title="Top"><img border=0 src="https://static.igem.org/mediawiki/2013/0/0b/Top-1.png"></a>
</div>
</div>
</html>
</html>

Latest revision as of 00:50, 29 October 2013


Mathematic Analysis on AlkSensor

Contents

Model Objective


To perform a mathematic analysis on AlkSensor, to find out the relationship between AlkSensor’s input and output, to find ways to specifically regulate the function of AlkSensor, to find the factors affecting sensor’s leakage.


Problem Description


AlkSensor is composed of protein ALKR and promoter alkM. Genes of ALKR and alkM were synthesized and constructed into a plasmid, as shown in figure 1. As mentioned in the introduction, protein ALKR is a transcription factor which can recognize certain alkanes. ALKR is under a constitutive promoter and is constitutively expressed. Alkane molecules is recognized by ALKR and a dimerized ALKR-alkane complex is formed. Subsequently the reporter’s promoter alkM is induced by the complex and the genes in the downstream of alkM are expressed. However, promoter alkM could also be slightly induced when inducers are absent from the system. There are possibilities that dierized ALKR binds on promoter alkM and activates the downstream genes.

Figure 1. Scheme of AlkSensor’s mechanism


Problem Abstraction


The mechanism of AlkSensor can be represented by a set of chemical reactions, as shown below.

TJU-F-11.png
Figure 2. Chemical reactions of AlkSensor
  • Reaction 1: the generation of protein ALKR.
  • Reaction 2: the degradation of protein ALKR
  • Reaction 3: the dimerization of protein ALKR(AR2-A denotes the dimerized ALKR-alkane complex)
  • Reaction 4: the combination of dimerized ALKR with alkane(Pm denotes promoter alkM)
  • Reaction 5: the combination of dimerized ALKR-alkane complex with promoter alkM(Pm’ denotes the promoter alkM binding with inducer)
  • Reaction 6: the combination of dimerized ALKR with promoter alkM(Pm’’ denotes the promoter alkM binding with dimerized ALKR)
  • Reaction 7: the generation of reporter (RNAP denote the RNA polymerase, Rp denotes reporter )
  • Reaction 8: the leakage of AlkSensor
  • Reaction 9: the the degradation of reporter


Assumption


There are several assumptions in the abstraction of AlkSensor.

  • Because the genes of ALKR is under a constitutive promoter, we assume that the generation rate of ALKR is constant.
  • The four binding reactions, i.e, reaction 3, 4, 5 and 6 are fast reversible reactions since DNA binding and unbinding of the repressors dimers and the dimerization itself occur within seconds, whereas the synthesis (transcription, translation, folding) and degradation of monomers takes minutes to sometimes an hour.

Reaction 3, 4, 5 and 6 can be integrated into two reversible questions. Therefore, the 9 reactions can be simplified as 7 reactions.

TJU-F-12.png
Figure 3. Simplified form of AlkSensor’s chemical reactions


Model Development


Reaction 3 and 4 are fast reversible reactions because DNA binding and unbinding of the transcription factor dimers and the dimerization itself occur within seconds. However the synthesis and degradation of monomers takes minutes to sometimes an hour. Therefore, the three reactions are in equilibrium and the steady state concentrations are given in terms of the equilibrium constants. The accumulation rates of ALKR and reporter are subject to their generation rate and degradation rate, which is described by two ordinary differential equations (ODE). So the mechanism of AlkSensor can be described by a mathematical model consisting of 2 ordinary differential equations (ODE) and three equilibrium equation. The input of AlkSensor can be defined as the concentration alkane, the output can be defined as the concentration of reporter protein, as shown in figure 2.

Figure 4. Mathematical representation of AlkSensor

In the first equation, there are four terms that determine the accumulation rate of reporter, the generation rate of reporter, the leakage rate of AlkSensor, the basal expression rate of reporter, the degradation rate of reporter. In the second equation, there are four terms that determine the accumulation rate of ALKR, the generation rate of ALKR, the degradation rate of ALKR, the dilution rate of ALKR. The last 3 equation describe the quantitative relation in reaction 3 and 4.

The [Pm’] and [Pm’’] are still intermediate variables, we next deduce the value of [Pm’] and [Pm’’].

TJU-F-1.png
TJU-F-2.png

The total concentration of promoter alkM is thought to be a constant, which is proportional to the copy number of alkM.

TJU-F-3.png

So we can get the expression of [Pm], [Pm’] and [Pm’’]

TJU-F-4.png
TJU-F-5.png
TJU-F-6.png

Then we substitute [Pm’] and [Pm’] with the terms on the right side of the equation.

TJU-F-7.png

When the leakage of ALKR’s promoter is little, the terms r can be ignored. When the system reach a steady state, the accumulation rate of Rp is zero. In such a condition, the response function of AlkSensor can be described as the equation below.

TJU-F-8.png

We can simplify the response function into an elegant form, in which x denotes the input, y denotes the output, a, b and c are three parameters.

TJU-F-9.png
  • k’: the expression rate constant of promoter alkM;
  • kd: the degradation constant of protein ALKR, thought to be a constant;
  • P0: the concentration of RNAP, constant;
  • [Pt]: the total concentration of promoter alkM, which is proportional to the copy number. This is easily changeable. We can change it through construct alkM on different plasmid with varying copy number.


Analysis & Discussion


The output has a monotone increasing relation with the input. When the input is large enough, the output reaches its maximum value, a, which is proportional to the amount of promoter alkM. When there is no input, the output has a smallest value, i.e. the value of leakage. We can easily obtain the value of dynamic range by dividing the maximum response into the value of leakage.

TJU-F-10.png

There are two parameters that can be rationally regulated, [Pt] and [ALKR]. We can control [Pt] by adjust the copy number of promoter alkM and control [ALKR] by adjusting the strength of its promoter. Therefore, we have the ability to regulate the value of a and b, that is to say, the maximum response, leakage and dynamic range of AlkSensor can all be rationally regulated.

Figure 5 shows a set of response curves of AlkSensor with different expression levels of ALKR. When there is no ALKR expressed, the response curve is a horizontal line, suggesting that AlkSensor has no ability to distinguish different inputs. With the increase of [ALKR], the leakage increase, however, the maximum response remains the same value, which makes the dynamic range decrease. When [ALKR] is too large, the response curve turns to be a horizontal line again. AlkSensor loses its function in this case.

Figure 5. The influence of [ALKR] on the response curve of AlkSensor


Achievement


In this model, we have

  • identified the pattern of AlkSensor’s response curve,
  • proposed a template for experimental data of AlkSensor,
  • found ways to evaluate AlkSensor,
  • developed strategies to rationally regulate AlkSensor.

Retrieved from "http://2013.igem.org/Team:Tianjin/Modeling"