Team:NJAU China/Modeling

From 2014.igem.org

(Difference between revisions)
 
(17 intermediate revisions not shown)
Line 20: Line 20:
padding:0;
padding:0;
text-decoration:none;
text-decoration:none;
 +
                font family:Andalus;
margin:0;}
margin:0;}
         body {
         body {
Line 25: Line 26:
             background:#a2dcf3;
             background:#a2dcf3;
             color: #5a5a5a;
             color: #5a5a5a;
-
             font-size: 12px;
+
             font-size: 14px;
 +
            font-family:Andalus;
         }
         }
pre{
pre{
font-size:14px;
font-size:14px;
-
font-family:}
+
font-family:Andalus;
 +
}
         #myheader{
         #myheader{
             width:1366px;
             width:1366px;
Line 48: Line 51:
             padding-top: 100px;
             padding-top: 100px;
             margin-bottom: 20px;
             margin-bottom: 20px;
-
             font-family: Georgia, Times, "Times New Roman", serif;
+
             font-family: Georgia, Times, "Andalus", serif;
         }
         }
         h1, h2, h3, h4, h5 {
         h1, h2, h3, h4, h5 {
-
             font-family: Georgia, Times, "Times New Roman", serif;
+
             font-family: Georgia, Times, "Andalus", serif;
         }
         }
Line 114: Line 117:
         .container {
         .container {
-
             font-size: 14px;
+
             font-size: 16px;
             width:700px;
             width:700px;
             float:right;
             float:right;
Line 145: Line 148:
margin:0 0 20px 0;
margin:0 0 20px 0;
padding:0;
padding:0;
 +
}
 +
#center{
 +
text-align:center;
}
}
#my_abstract p{
#my_abstract p{
Line 165: Line 171:
margin:20px 0;}
margin:20px 0;}
#my_con p{
#my_con p{
-
font-size:14px;
+
font-size:16px;
 +
        font-family:Andalus;
text-indent:30px;
text-indent:30px;
margin:15px 0;}
margin:15px 0;}
-
.center{
 
-
margin:0 auto;
 
-
}
 
#main{width:1000px; height:auto; overflow:hidden; margin: auto; padding-bottom:10px;}
#main{width:1000px; height:auto; overflow:hidden; margin: auto; padding-bottom:10px;}
Line 236: Line 240:
           <li id="ackonwledgement_li"><a href="#">Acknowledgement</a>
           <li id="ackonwledgement_li"><a href="#">Acknowledgement</a>
                 <dl id="ackonwledgement">
                 <dl id="ackonwledgement">
-
                    <dd><a href="https://2014.igem.org/Team:NJAU_China/Inside Campus">Inside Campus</a></dd>
+
                      <dd><a href="https://2014.igem.org/Team:NJAU_China/Inside Campus">Inside Campus&Arrtibutions</a></dd>
                     <dd><a href="https://2014.igem.org/Team:NJAU_China/Outside Campus">Outside Campus</a></dd>
                     <dd><a href="https://2014.igem.org/Team:NJAU_China/Outside Campus">Outside Campus</a></dd>
                 </dl>
                 </dl>
Line 277: Line 281:
<b>Abstract</b> The process which this equipment “Copper Terminator” uptaking Cu<sup>2+</sup> is modeled by ordinary differential equations. Parameters in this models are both inferred from our experiments or obtained from literature. Second order Runge-Kutta method is utilized in simulation and the result shows that Cu<sup>2+</sup> is almost entirely  eliminated within 300 minutes from 1mg/L. Finally, Morris sensitivity analysis is carried out and shows that four controllable parameters are most important to the effect of elimination and this demonstrate the possibility of its application. </p>
<b>Abstract</b> The process which this equipment “Copper Terminator” uptaking Cu<sup>2+</sup> is modeled by ordinary differential equations. Parameters in this models are both inferred from our experiments or obtained from literature. Second order Runge-Kutta method is utilized in simulation and the result shows that Cu<sup>2+</sup> is almost entirely  eliminated within 300 minutes from 1mg/L. Finally, Morris sensitivity analysis is carried out and shows that four controllable parameters are most important to the effect of elimination and this demonstrate the possibility of its application. </p>
     <h1> 1.Introduction</h1>
     <h1> 1.Introduction</h1>
-
       <img src="https://static.igem.org/mediawiki/2014/c/cc/Qq1.png" class="center"/>
+
       <div style="text-align:center"><img src="https://static.igem.org/mediawiki/2014/c/cc/Qq1.png"/></div>
<p> In our “Copper Terminator” (Figure 3.2.1), <i>E.coli</i> reproduce in a relatively stable environment in which the temperature is set constant by a temperature sensor , a heater and a microcontroller. Polymers are covered on a framework to keep <i>E.coli</i> form spreading into outside space but let Cu<sup>2+</sup> pass through. Cu<sup>2+</sup> which have entered this equipment is uptaken by <i>E.coli</i>. </p>
<p> In our “Copper Terminator” (Figure 3.2.1), <i>E.coli</i> reproduce in a relatively stable environment in which the temperature is set constant by a temperature sensor , a heater and a microcontroller. Polymers are covered on a framework to keep <i>E.coli</i> form spreading into outside space but let Cu<sup>2+</sup> pass through. Cu<sup>2+</sup> which have entered this equipment is uptaken by <i>E.coli</i>. </p>
     <h1> 2.Assumptions</h1>
     <h1> 2.Assumptions</h1>
Line 288: Line 292:
<b> a)Variables and Parameters</b> </p>
<b> a)Variables and Parameters</b> </p>
<p> i.Variables</p>
<p> i.Variables</p>
-
     <img src=" https://static.igem.org/mediawiki/2014/2/2c/Qq2.png" />
+
     <div style="text-align:center"><img src=" https://static.igem.org/mediawiki/2014/2/2c/Qq2.png" /></div>
<p> ii.Parameters</p>
<p> ii.Parameters</p>
-
     <img src=" https://static.igem.org/mediawiki/2014/5/56/Qq3.1.png" />
+
     <div style="text-align:center"><img src=" https://static.igem.org/mediawiki/2014/5/56/Qq3.1.png" /></div>
-
<img src=" https://static.igem.org/mediawiki/2014/5/57/Qq3.2.png" />
+
<div style="text-align:center"><img src=" https://static.igem.org/mediawiki/2014/5/57/Qq3.2.png" /></div>
<p>
<p>
<b> b)Model Development</b> </p>
<b> b)Model Development</b> </p>
Line 297: Line 301:
<img src=" https://static.igem.org/mediawiki/2014/c/c6/Qqt1.png" />
<img src=" https://static.igem.org/mediawiki/2014/c/c6/Qqt1.png" />
<p> Similarly, a back diffusion from j to i is described as:</p>
<p> Similarly, a back diffusion from j to i is described as:</p>
-
<img src=" https://static.igem.org/mediawiki/2014/8/86/Qqt2.png" />
+
<div style="text-align:center"><img src=" https://static.igem.org/mediawiki/2014/8/86/Qqt2.png" /></div>
<p> The change in the substance of each compartment can be derivated as :</p>
<p> The change in the substance of each compartment can be derivated as :</p>
-
<img src=" https://static.igem.org/mediawiki/2014/0/0e/Qqt3.png" />
+
<div style="text-align:center"><img src=" https://static.igem.org/mediawiki/2014/0/0e/Qqt3.png" /></div>
<p> In our model, taking the effective surface area of each “Copper Terminator” into consideration, we derivate the change of Cu<sup>2+</sup> concentration outside the equipment as:</p>
<p> In our model, taking the effective surface area of each “Copper Terminator” into consideration, we derivate the change of Cu<sup>2+</sup> concentration outside the equipment as:</p>
-
<img src=" https://static.igem.org/mediawiki/2014/1/17/Qq4.png" />
+
<div style="text-align:center"><img src=" https://static.igem.org/mediawiki/2014/1/17/Qq4.png" /></div>
<p> Parameter <i>k</i> here has a different unit from previous ones.</p>
<p> Parameter <i>k</i> here has a different unit from previous ones.</p>
<p> The concentration of Cu<sup>2+</sup> inside is increased by diffusion while decreased by E.coli uptaking. So that the change of Cu<sup>2+</sup> concentration inside can be described as :</p>
<p> The concentration of Cu<sup>2+</sup> inside is increased by diffusion while decreased by E.coli uptaking. So that the change of Cu<sup>2+</sup> concentration inside can be described as :</p>
-
<img src=" https://static.igem.org/mediawiki/2014/6/65/Qqt4.png" />
+
<div style="text-align:center"><img src=" https://static.igem.org/mediawiki/2014/6/65/Qqt4.png" /></div>
<p> The growth of E.coli follow logistic equation :</p>
<p> The growth of E.coli follow logistic equation :</p>
-
<img src=" https://static.igem.org/mediawiki/2014/e/e4/Qqt5.png" />
+
<div style="text-align:center"><img src=" https://static.igem.org/mediawiki/2014/e/e4/Qqt5.png" /></div>
<p>
<p>
<b> c)Model Simulation</b> </p>
<b> c)Model Simulation</b> </p>
  <p> Simulation is carried out by second order Runge-Kutta method with a step of 0.1 min and 3000 times iteration. The parameters are both from our experiment and some literatures.</p>
  <p> Simulation is carried out by second order Runge-Kutta method with a step of 0.1 min and 3000 times iteration. The parameters are both from our experiment and some literatures.</p>
-
<img src=" https://static.igem.org/mediawiki/2014/7/75/Qq5.png" />
+
<div style="text-align:center"><img src=" https://static.igem.org/mediawiki/2014/a/a6/Qq5.jpg" /></div>
<p> Figure 3.2.2 shows the Cu<sup>2+</sup> concentration change with time both inside and outside.</p>
<p> Figure 3.2.2 shows the Cu<sup>2+</sup> concentration change with time both inside and outside.</p>
-
  <img src=" https://static.igem.org/mediawiki/2014/6/62/Qq6.png" />
+
  <div style="text-align:center"><img src=" https://static.igem.org/mediawiki/2014/1/1a/Qq6.jpg" /></div>
<p> Figure 3.2.3 shows the growth curve of <i>E.coli</i>.</p>
<p> Figure 3.2.3 shows the growth curve of <i>E.coli</i>.</p>
<p>
<p>
Line 321: Line 325:
         <p> Standard deviation of the effect, σ, which estimates the higher-order characteristics of the parameter (such as curvatures and interactions).</p>
         <p> Standard deviation of the effect, σ, which estimates the higher-order characteristics of the parameter (such as curvatures and interactions).</p>
         <p> We define t<sub>0.5</sub> as the time when Cu<sup>2+</sup> concentration outside reaches 0.5mg/L, which is the maximum value allowed in city water. Monte Carlo method is applied to generate 80 groups of 7 parameters randomly with certain distributions. t<sub>0.5</sub> of each sample can be calculated and then they are undergoing Morris sensitivity analysis.</p>
         <p> We define t<sub>0.5</sub> as the time when Cu<sup>2+</sup> concentration outside reaches 0.5mg/L, which is the maximum value allowed in city water. Monte Carlo method is applied to generate 80 groups of 7 parameters randomly with certain distributions. t<sub>0.5</sub> of each sample can be calculated and then they are undergoing Morris sensitivity analysis.</p>
-
<img src=" https://static.igem.org/mediawiki/2014/d/d3/Qq7.png" />
+
<div style="text-align:center"><img src="https://static.igem.org/mediawiki/2014/f/f7/Qq7.jpg" /></div>
     <p> According to figure 3.2.4, the most significant parameters which have most effect on the result in this model are R, h, D, and k. Those parameters are all controllable parameters. The result of Morris analysis gives us a direction on design and application of these “Copper Terminator” equipment: As long as we set proper R, h, D, and k, we will make an outstanding improvement on the performance.</p>
     <p> According to figure 3.2.4, the most significant parameters which have most effect on the result in this model are R, h, D, and k. Those parameters are all controllable parameters. The result of Morris analysis gives us a direction on design and application of these “Copper Terminator” equipment: As long as we set proper R, h, D, and k, we will make an outstanding improvement on the performance.</p>
       <h3>Refrences:</h3>
       <h3>Refrences:</h3>

Latest revision as of 23:06, 17 October 2014

Home

Abstract The process which this equipment “Copper Terminator” uptaking Cu2+ is modeled by ordinary differential equations. Parameters in this models are both inferred from our experiments or obtained from literature. Second order Runge-Kutta method is utilized in simulation and the result shows that Cu2+ is almost entirely eliminated within 300 minutes from 1mg/L. Finally, Morris sensitivity analysis is carried out and shows that four controllable parameters are most important to the effect of elimination and this demonstrate the possibility of its application.

1.Introduction

In our “Copper Terminator” (Figure 3.2.1), E.coli reproduce in a relatively stable environment in which the temperature is set constant by a temperature sensor , a heater and a microcontroller. Polymers are covered on a framework to keep E.coli form spreading into outside space but let Cu2+ pass through. Cu2+ which have entered this equipment is uptaken by E.coli.

2.Assumptions

Due to our experiment, Cu2+ concentration have little influence on the growth of E.coli. The critical values of tolerance of Cu2+ for DH5α and BL21 is 165mg/L and 160mg/L which are all beyond the normal concentration in water. We assume that the growth of E.coli is not influenced by Cu2+ concentration.

Since we have polymers covering the surface of this equipment, the environment inside is relatively independent from the environment outside. The growth of E.coli can hardly be disturbed by environment outside. The assumption is that the growth of E.coli is subject to logistic model, in other words, the environmental resistance only comes from the competition on nutrition within the population.

How Cu2+ is transported into E.coli is still undiscovered. So that it is difficult to come up with an exact description by mathematical equations for this process. To simplify this model, we assume that the rate of Cu2+ uptaken is not related to the physiological status of E.coli. They are equal under different stage of a single E.coli.

For this model here, we assume uniformity of Cu2+ and E.coli. in the solution.

3.Model Details

a)Variables and Parameters

i.Variables

ii.Parameters

b)Model Development

The simple model which describes forward diffusion from compartment i to j :

Similarly, a back diffusion from j to i is described as:

The change in the substance of each compartment can be derivated as :

In our model, taking the effective surface area of each “Copper Terminator” into consideration, we derivate the change of Cu2+ concentration outside the equipment as:

Parameter k here has a different unit from previous ones.

The concentration of Cu2+ inside is increased by diffusion while decreased by E.coli uptaking. So that the change of Cu2+ concentration inside can be described as :

The growth of E.coli follow logistic equation :

c)Model Simulation

Simulation is carried out by second order Runge-Kutta method with a step of 0.1 min and 3000 times iteration. The parameters are both from our experiment and some literatures.

Figure 3.2.2 shows the Cu2+ concentration change with time both inside and outside.

Figure 3.2.3 shows the growth curve of E.coli.

d)Global Sensitivity Analysis

Global sensitivity analysis is a method to analyse all the parameters at one time to find out the influence on the result for each parameter and the interaction between those parameters.

Morris (1991) proposed conducting individually randomized experiments that evaluate the effect of changing one parameter at a time. Each input may assume a discrete number of values, called levels, that are selected from within an allocated range of variation for the parameter.

For each parameter, two sensitivity measures are proposed by Morris (1991): (1) the mean, μ, which estimates the overall effect of the parameter on a given output; and (2) the.

Standard deviation of the effect, σ, which estimates the higher-order characteristics of the parameter (such as curvatures and interactions).

We define t0.5 as the time when Cu2+ concentration outside reaches 0.5mg/L, which is the maximum value allowed in city water. Monte Carlo method is applied to generate 80 groups of 7 parameters randomly with certain distributions. t0.5 of each sample can be calculated and then they are undergoing Morris sensitivity analysis.

According to figure 3.2.4, the most significant parameters which have most effect on the result in this model are R, h, D, and k. Those parameters are all controllable parameters. The result of Morris analysis gives us a direction on design and application of these “Copper Terminator” equipment: As long as we set proper R, h, D, and k, we will make an outstanding improvement on the performance.

Refrences:

[1] Keen, Robert E., and James D. Spain. Computer simulation in biology: a BASIC introduction. John Wiley & Sons, Inc., 1994.

[2] Morris, Max D. "Factorial sampling plans for preliminary computational experiments." Technometrics 33.2 (1991): 161-174.