Difference between revisions of "Team:HUST-China/Modeling on Ecosystem Level"

 
(46 intermediate revisions by 4 users not shown)
Line 79: Line 79:
 
  .header{position: fixed;  
 
  .header{position: fixed;  
 
      width: 1348px;
 
      width: 1348px;
−
                            text-align:center;
+
                text-align:center;
 
    height:200px;
 
    height:200px;
 
    margin-top: 0px;
 
    margin-top: 0px;
Line 88: Line 88:
 
    z-index:1;}
 
    z-index:1;}
  
−
.first-menu{float: left;
+
.first-menu{float: left;
−
margin-left: 1%;
+
 
  font-size:19px;  
 
  font-size:19px;  
−
      max-width: 22%;
+
      max-width: 20%;
 
                 padding-top:5px;}
 
                 padding-top:5px;}
  
 
  .other-menu{float: left;
 
  .other-menu{float: left;
−
  margin-left:1.2%;
+
  margin-left:1%;
 
  font-size:19px;
 
  font-size:19px;
−
      max-width: 22%;
+
      max-width: 20%;
 
                     padding-top:5px;}
 
                     padding-top:5px;}
 +
  
 
  .navbar .nav > li .dropdown-menu {margin: 0;}
 
  .navbar .nav > li .dropdown-menu {margin: 0;}
Line 141: Line 141:
 
     }
 
     }
 
     .pic_a{float: left;
 
     .pic_a{float: left;
−
             margin-top:100px;
+
             margin-top:35px;
 
             width: 1349px;
 
             width: 1349px;
 
             text-align:center;
 
             text-align:center;
Line 230: Line 230:
 
 
 
<!--悬浮菜单-->
 
<!--悬浮菜单-->
 +
<!--悬浮菜单-->
 
<div class="navbar">
 
<div class="navbar">
 
    <div class="navbar-inner">
 
    <div class="navbar-inner">
Line 236: Line 237:
 
              <li class="first-menu"><a href="https://2015.igem.org/Team:HUST-China">HOME</a>
 
              <li class="first-menu"><a href="https://2015.igem.org/Team:HUST-China">HOME</a>
 
                </li>      
 
                </li>      
−
                <li class="dropdown other-menu" id="accountmenu"><a href="https://2015.igem.org/Team:HUST-China/Background">PROJECT<b class="caret"></b></a>
+
                <li class="dropdown other-menu" id="accountmenu"><a href="https://2015.igem.org/Team:HUST-China/Project">PROJECT<b class="caret"></b></a>
 
                <ul class="dropdown-menu">
 
                <ul class="dropdown-menu">
 
                                 <li><a href="https://2015.igem.org/Team:HUST-China/Background">Background</a></li>
 
                                 <li><a href="https://2015.igem.org/Team:HUST-China/Background">Background</a></li>
Line 246: Line 247:
 
                         </li>
 
                         </li>
  
−
                         <li class="dropdown other-menu" id="accountmenu"><a href="https://2015.igem.org/Team:HUST-China/Results">WETLAB<b class="caret"></b></a>
+
                         <li class="dropdown other-menu" id="accountmenu"><a href="https://2015.igem.org/Team:HUST-China/Wetlab">WETLAB<b class="caret"></b></a>
 
                <ul class="dropdown-menu">
 
                <ul class="dropdown-menu">
 
<li><a href="https://2015.igem.org/Team:HUST-China/Results">Results</a></li>
 
<li><a href="https://2015.igem.org/Team:HUST-China/Results">Results</a></li>
Line 258: Line 259:
 
                         </li>
 
                         </li>
  
−
                         <li class="dropdown other-menu" id="accountmenu"><a href="https://2015.igem.org/Team:HUST-China/Basic_part">PARTS<b class="caret"></b></a>
+
                         <li class="dropdown other-menu" id="accountmenu"><a href="https://2015.igem.org/Team:HUST-China/Parts">PARTS<b class="caret"></b></a>
 
                <ul class="dropdown-menu">                          
 
                <ul class="dropdown-menu">                          
 
                        <li><a href="https://2015.igem.org/Team:HUST-China/Basic_part">Basic Parts</a></li>
 
                        <li><a href="https://2015.igem.org/Team:HUST-China/Basic_part">Basic Parts</a></li>
 
                        <li class="divider"></li>
 
                        <li class="divider"></li>
−
                        <li><a href="https://2015.igem.org/Team:HUST-China/Basic_part#2">Composite Parts</a></li>   
+
                        <li><a href="https://2015.igem.org/Team:HUST-China/Basic_part#2">Composite Parts</a></li>                          <li class="divider"></li>
 +
                        <li><a href="https://2015.igem.org/Team:HUST-China/Part_Collection">Part Collection</a></li>   
 
                    </ul>
 
                    </ul>
 
                         </li>
 
                         </li>
Line 268: Line 270:
 
                <li class="dropdown other-menu" id="accountmenu"><a href="https://2015.igem.org/Team:HUST-China/Modeling">MODELING<b class="caret"></b></a>  
 
                <li class="dropdown other-menu" id="accountmenu"><a href="https://2015.igem.org/Team:HUST-China/Modeling">MODELING<b class="caret"></b></a>  
 
                             <ul class="dropdown-menu">
 
                             <ul class="dropdown-menu">
 +
                                        <li><a href="https://2015.igem.org/Team:HUST-China/Modeling">Overiew</a></li>
 +
                                        <li class="divider"></li>
 
                          <li><a href="https://2015.igem.org/Team:HUST-China/Modeling on Cellular Level">Modeling on Cellular Level</a></li>
 
                          <li><a href="https://2015.igem.org/Team:HUST-China/Modeling on Cellular Level">Modeling on Cellular Level</a></li>
 
                          <li class="divider"></li>
 
                          <li class="divider"></li>
Line 275: Line 279:
  
 
                                 
 
                                 
−
                 
+
                <li class="dropdown other-menu" id="accountmenu"><a href="https://2015.igem.org/Team:HUST-China/Practices">HUMAN PRACTICES</a>
−
                <li class="dropdown other-menu" id="accountmenu">
+
                                       
−
                    <a class="dropdown-toggle" data-toggle="dropdown" href="https://2015.igem.org/Team:HUST-China/Practices">HUMAN PRACTICES</a>
+
−
                                         
+
 
                </li>
 
                </li>
  
−
              <li class="dropdown other-menu" id="accountmenu">
+
              <li class="dropdown other-menu" id="accountmenu"><a href="https://2015.igem.org/Team:HUST-China/Safety">OTHERS<b class="caret"></b></a>
−
                    <a class="dropdown-toggle" data-toggle="dropdown" href="https://2015.igem.org/Team:HUST-China/Safety">OTHERS<b class="caret"></b></a>
+
 
                    <ul class="dropdown-menu">
 
                    <ul class="dropdown-menu">
 
                        <li><a href="https://2015.igem.org/Team:HUST-China/Safety">Safety</a></li>
 
                        <li><a href="https://2015.igem.org/Team:HUST-China/Safety">Safety</a></li>
Line 303: Line 304:
 
           $('.dropdown-toggle').dropdown(); }); </script>
 
           $('.dropdown-toggle').dropdown(); }); </script>
 
</div>
 
</div>
 +
 
<!--标题栏-->
 
<!--标题栏-->
 
  <div id="pic" >
 
  <div id="pic" >
−
  <img class="title" src="https://static.igem.org/mediawiki/2015/c/c4/HUST_MODELING.png"/>
+
  <img class="title" src="https://static.igem.org/mediawiki/2015/3/31/HUST_MODELING_7.PNG"/>
 
  <br>
 
  <br>
−
  <div class="pic_a" >
+
  <div class="pic_a" >
 +
<h4 align="center" style="color:white"><b>click it~</b></h4>
 
      <img style="cursor:pointer;" id="to_des" src="https://static.igem.org/mediawiki/2015/8/80/White.png"/>
 
      <img style="cursor:pointer;" id="to_des" src="https://static.igem.org/mediawiki/2015/8/80/White.png"/>
 
  </div>
 
  </div>
Line 319: Line 322:
 
  <li class="li"><a href="#2" class="btn btn-default btn-lg">The “wake-up” problem</a></li>
 
  <li class="li"><a href="#2" class="btn btn-default btn-lg">The “wake-up” problem</a></li>
 
  <li class="li"><a href="#3" class="btn btn-default btn-lg">The permeation problem</a></li>
 
  <li class="li"><a href="#3" class="btn btn-default btn-lg">The permeation problem</a></li>
−
  <li class="li"><a href="#4" class="btn btn-default btn-lg">Robustness Analysis</a></li>
+
  <li class="li"><a href="#4" class="btn btn-default btn-lg">R&S Analysis
 +
</a></li>
 
  </ul>
 
  </ul>
 
  </div>
 
  </div>
Line 363: Line 367:
 
  <h3 style="color:black" align="left"><b>Parameters:</b></h3><br>
 
  <h3 style="color:black" align="left"><b>Parameters:</b></h3><br>
 
             <p>The description of parameters, their values and the references involved in this model are listed in a table.  
 
             <p>The description of parameters, their values and the references involved in this model are listed in a table.  
−
            The parameters table (2) can be downloaded<a href="#"><b> here.</b></a>
+
<h4 style="width:800px;margin-top:30px;text-align:center" align:center;>Brief Parameter Table of Diffusion Model</h4>
 +
<table border="1">
 +
            <tr>
 +
<th> Parameter </th>
 +
<th> Description </th>
 +
</tr>
 +
<tr>
 +
<td> a(x,y,z,t) </td>
 +
<td> the concentration of Mcfp-3 </td>
 +
</tr>
 +
<tr>
 +
<td> A(x,y,z,t) </td>
 +
<td> the concentration of fixed bacterial </td>
 +
</tr>
 +
<tr>
 +
<td> B </td>
 +
<td> the vector of prerequisite coefficient </td>
 +
</tr>
 +
<tr>
 +
<td> C(x,y,z,t) </td>
 +
<td> the concentration of free bacterial </td>
 +
</tr>
 +
<tr>
 +
<td> b<sub>i</sub> </td>
 +
<td> the permeation coefficient of on the i direction </td>
 +
</tr>
 +
<tr>
 +
<td> eb </td>
 +
<td> the rate of Mcfp-3 expression </td>
 +
</tr>
 +
<tr>
 +
<td> ed </td>
 +
<td> the rate of Mcfp-3 degradation </td>
 +
</tr>
 +
<tr>
 +
<td> l </td>
 +
<td> the edge length of each lattice </td>
 +
</tr>
 +
<tr>
 +
<td> T<sub>0</sub> </td>
 +
<td> the time unit of prerequisite </td>
 +
</tr>
 +
<tr>
 +
<td> α </td>
 +
<td> the possibility of effective collision </td>
 +
</tr>
 +
<tr>
 +
<td> λ </td>
 +
<td> average length of steps of Brownian motion </td>
 +
</tr>
 +
<tr>
 +
<td> τ </td>
 +
<td> average free time of Brownian motion </td>
 +
</tr>
 +
<tr>
 +
<td> N<sub>A</sub> </td>
 +
<td> the Avogadro constant </td>
 +
</tr>
 +
<tr>
 +
<td> T </td>
 +
<td> the absolute temperature </td>
 +
</tr>
 +
<tr>
 +
<td> R </td>
 +
<td> the Rydberg constant </td>
 +
</tr>
 +
<tr>
 +
<td> r & d </td>
 +
<td> the radius and diameter of our yeast </td>
 +
</tr>
 +
<tr>
 +
<td> ν </td>
 +
<td> the speed of the prerequisite </td>
 +
</tr>
 +
<tr>
 +
<td> ρ </td>
 +
<td> the density of our yeast </td>
 +
</tr>
 +
<tr>
 +
<td> η </td>
 +
<td> the dynamic viscosity coefficient </td>
 +
</tr>
 +
<tr>
 +
<td> C<sub>i</sub> </td>
 +
<td> the total number of bacterial of the 26 neighboring lattices of the lattice (x,y,z) </td>
 +
</tr>
 +
</table>
 +
 
 +
            More detailed parameters table of diffusion model can be downloaded<a href="https://static.igem.org/mediawiki/2015/d/df/Parameters_table.zip"><b> here.</b></a>
 
</p>
 
</p>
  
 
    <h3 style="color:black" align="left"><b>Formulary:</b></h3><br>
 
    <h3 style="color:black" align="left"><b>Formulary:</b></h3><br>
−
            <img class="gongshi" src="https://static.igem.org/mediawiki/2015/4/4f/Pm1.png">
+
          <div >
−
            <h4 style="float:right;">(2.1)</h4><br>
+
             <img style="width:850px" class="gongshi" src="https://static.igem.org/mediawiki/2015/d/dd/HUST_gongshi2.jpg">
−
             <img class="gongshi" src="https://static.igem.org/mediawiki/2015/4/4c/Pm2.png">
+
          </div>
−
            <h4 style="float:right;">(2.2)</h4><br>
+
−
            <img class="gongshi" src="https://static.igem.org/mediawiki/2015/f/f9/Pm3.png">
+
−
            <h4 style="float:right;">(2.3)</h4><br>
+
−
            <img class="gongshi" src="https://static.igem.org/mediawiki/2015/f/f8/Pm4.png">
+
−
            <h4 style="float:right;">(2.4)</h4><br>
+
−
            <img class="gongshi" src="https://static.igem.org/mediawiki/2015/9/94/Pm5.png">
+
−
            <h4 style="float:right;">(2.5)</h4><br>
+
  
 
         <h3 style="color:black" align="left"><b>Monte Carlo Simulation:</b></h3><br>
 
         <h3 style="color:black" align="left"><b>Monte Carlo Simulation:</b></h3><br>
Line 396: Line 481:
 
             <p>Comparing the effect of sedimentation due to gravity and the Brownian motion, we reach a conclusion that the latter one outweigh the former one (about 100 times). And therefore the most effective factor of the motion of Euk.Cement in the sands is Brownian motion.</p>
 
             <p>Comparing the effect of sedimentation due to gravity and the Brownian motion, we reach a conclusion that the latter one outweigh the former one (about 100 times). And therefore the most effective factor of the motion of Euk.Cement in the sands is Brownian motion.</p>
 
              
 
              
−
             <p><a href=""><b>The code of the Monte Carlo Simulation can be downloaded here.</b></a>Its flow chart is</p>
+
             <p>The code of the Monte Carlo Simulation can be downloaded<a href=""><b> here.</b></a> Its flow chart is</p>
 
             <img class="picture" src="https://static.igem.org/mediawiki/2015/d/d3/Q.png"><br>
 
             <img class="picture" src="https://static.igem.org/mediawiki/2015/d/d3/Q.png"><br>
 
             <p>After running the Monte Carlo Simulation with 106 yeasts for multiple times, we finally determine the values of the following parameters</p>
 
             <p>After running the Monte Carlo Simulation with 106 yeasts for multiple times, we finally determine the values of the following parameters</p>
Line 402: Line 487:
 
   
 
   
 
         <h3 style="color:black" align="left"><b>Lattices Method:</b></h3><br>
 
         <h3 style="color:black" align="left"><b>Lattices Method:</b></h3><br>
−
             <p>With sufficient parameters, we could finally simulate the permeation process. The code of our permeation model could be downloaded<a href=""><b>  here.</b></a></p>
+
             <p>With sufficient parameters, we could finally simulate the permeation process. The code of our permeation model could be downloaded <a href=""><b>  here.</b></a></p>
 
             <p>We chose to use numerical solution to solve our Permeation model. We divided the space into the same lattices, and each lattice has only one phase. Hence, each variable in the same coordinate is considered the same. The timeline was divided into the minimum and discrete values. By calculating the variables in each lattice through the discretized timeline, we could attain the values of a(x,y,z,t),A(x,y,z,t)and C(x,y,z,t). If the lattices were small enough, the PDEs above could be replaced by Lattices Method. The equations of Lattices Method (Cellular Automaton) are listed below.</p>
 
             <p>We chose to use numerical solution to solve our Permeation model. We divided the space into the same lattices, and each lattice has only one phase. Hence, each variable in the same coordinate is considered the same. The timeline was divided into the minimum and discrete values. By calculating the variables in each lattice through the discretized timeline, we could attain the values of a(x,y,z,t),A(x,y,z,t)and C(x,y,z,t). If the lattices were small enough, the PDEs above could be replaced by Lattices Method. The equations of Lattices Method (Cellular Automaton) are listed below.</p>
−
             <img class="gongshi" src="https://static.igem.org/mediawiki/2015/1/1d/Pm14.png">
+
            <div >
−
             <h4 style="float:right;">(2.9)</h4><br>
+
             <img style="width:850px" class="gongshi" src="https://static.igem.org/mediawiki/2015/5/50/HUST_gongshi.jpg">
−
            <img class="gongshi" src="https://static.igem.org/mediawiki/2015/e/ed/Pm15.png">
+
          </div>
−
            <h4 style="float:right;">(2.10)</h4><br>
+
             <p>With the equations of Cellular Automaton as well as Boundary equations and Initial state equations, we could attain the values of the concentration of Mcfp-3, the number of fixed Euk.Cement and free Euk.Cement at any time and at any coordinate.</p>
−
            <img class="gongshi" src="https://static.igem.org/mediawiki/2015/8/81/Pm16.png">
+
            <p>We ran the simulation with a server and printed a heat map</p>
−
            <h4 style="float:right;">(2.11)</h4><br>
+
              <img class="picture" style="width:650px;height:auto" src="https://static.igem.org/mediawiki/2015/8/89/HUST_heat_map3.gif">
−
            <p>With the equations of Cellular Automaton as well as Boundary equations and Initial state equations, we could attain the values of the concentration of Mcfp-3, the number of fixed Euk.Cement and free Euk.Cement at any time and at any coordinate.<br>The heat map for diffusion of Euk.Cement and Mcfp-3 were generated based on the equations above as well as Lattice Method. Results were shown as in Figure XXX.</p>
+
    <div style="font-size:14px;margin-bottom:30px" align="center">Figure 8: Diffusion Model of Free Euk.Cement Cells</div>
−
             <img class="picture" src=""><br>
+
            <img class="picture" style="width:650px;height:auto" src="https://static.igem.org/mediawiki/2015/7/7d/HUST_heat_map2.gif">
 +
    <div style="font-size:14px;margin-bottom:30px" align="center">Figure 8: Diffusion Model of Bound Euk.Cement Cells</div>
 +
             <img class="picture" style="width:650px;height:auto" src="https://static.igem.org/mediawiki/2015/6/6f/HUST_heat_map1.gif">
 +
    <div style="font-size:14px;margin-bottom:30px" align="center">Figure 7: Diffusion Model of Secreted Mcfp-3 Protein</div>
 +
<br>
  
−
         <h3 style="color:black" align="left"><b>Results:</b></h3><br>
+
         <h4 style="color:black" align="left"><b>Results:</b></h4><br>
 
             <p>From the permeation process shown in the heat map, we can derive the following conclusion:<br>
 
             <p>From the permeation process shown in the heat map, we can derive the following conclusion:<br>
 
- The Euk.Cement permeates very efficiently over a large space. Therefore, perhaps we have to enclose it.<br>
 
- The Euk.Cement permeates very efficiently over a large space. Therefore, perhaps we have to enclose it.<br>
Line 430: Line 519:
 
The total number of Euk.Cement bound with sands in 30s = 28224;<br>
 
The total number of Euk.Cement bound with sands in 30s = 28224;<br>
 
The total number of Euk.Cement drifted away in 30s = 74778;<br></p>
 
The total number of Euk.Cement drifted away in 30s = 74778;<br></p>
−
<img class="picture2" src="https://static.igem.org/mediawiki/2015/a/a6/HUST_model2.3.png"><br><br>
+
<table border="0" style="width:300px" align:center>
 +
            <tr>
 +
<th> val(;:;:1)= </th>
 +
</tr>
 +
 
 +
<tr>
 +
<td>  </td>
 +
<td> 0 </td>
 +
                                <td> 8 </td>
 +
<td> 0 </td>
 +
</tr>
 +
<tr>
 +
<td>  </td>
 +
<td> 4 </td>
 +
<td> 12439 </td>
 +
<td> 7 </td>
 +
 
 +
</tr>
 +
<tr>    <td>  </td>
 +
<td> 0 </td>
 +
<td> 4 </td>
 +
<td> 0 </td>
 +
</tr>
 +
<tr>
 +
<th> val(;:;:2)= </th>
 +
</tr>
 +
            <tr>
 +
<td>  </td>
 +
<td> 5 </td>
 +
                                <td> 12474 </td>
 +
<td> 4 </td>
 +
</tr>
 +
            <tr>
 +
<td>  </td>
 +
<td> 12456 </td>
 +
                                <td> 1000000 </td>
 +
<td> 12659 </td>
 +
</tr>
 +
            <tr>
 +
<td>  </td>
 +
<td> 2 </td>
 +
                                <td> 12327 </td>
 +
<td> 9 </td>
 +
</tr>
 +
<th> val(;:;:3)= </th>
 +
</tr>
 +
            <tr>
 +
<td>  </td>
 +
<td> 0 </td>
 +
                                <td> 7 </td>
 +
<td> 0 </td>
 +
</tr>
 +
            <tr>
 +
<td>  </td>
 +
<td> 11 </td>
 +
                                <td> 12423 </td>
 +
<td> 4 </td>
 +
</tr>
 +
            <tr>
 +
<td>  </td>
 +
<td> 0 </td>
 +
                                <td> 4 </td>
 +
<td> 0 </td>
 +
</tr>
 +
</table>
 
   
 
   
 
  <p>When the α of MC = 0.0013, the α of LM = 0.0365:<br>
 
  <p>When the α of MC = 0.0013, the α of LM = 0.0365:<br>
 
The total number of Euk.Cement bound with sands in 30s = 36464;<br>
 
The total number of Euk.Cement bound with sands in 30s = 36464;<br>
 
The total number of Euk.Cement drifted away in 30s = 74925;<br></p>
 
The total number of Euk.Cement drifted away in 30s = 74925;<br></p>
−
<img class="picture2" src="https://static.igem.org/mediawiki/2015/a/ad/HUST_model2.4.png"><br><br>
+
<table border="0" style="width:300px" align:center>
 +
            <tr>
 +
<th> val(;:;:1)= </th>
 +
</tr>
 +
 
 +
<tr>
 +
<td>  </td>
 +
<td> 0 </td>
 +
                                <td> 8 </td>
 +
<td> 0 </td>
 +
</tr>
 +
<tr>
 +
<td>  </td>
 +
<td> 3 </td>
 +
<td> 12538 </td>
 +
<td> 9 </td>
 +
 
 +
</tr>
 +
<tr>    <td>  </td>
 +
<td> 0 </td>
 +
<td> 7 </td>
 +
<td> 0 </td>
 +
</tr>
 +
<tr>
 +
<th> val(;:;:2)= </th>
 +
</tr>
 +
            <tr>
 +
<td>  </td>
 +
<td> 5 </td>
 +
                                <td> 12553 </td>
 +
<td> 3 </td>
 +
</tr>
 +
            <tr>
 +
<td>  </td>
 +
<td> 12549 </td>
 +
                                <td> 1000000 </td>
 +
<td> 12550 </td>
 +
</tr>
 +
            <tr>
 +
<td>  </td>
 +
<td> 2 </td>
 +
                                <td> 12398 </td>
 +
<td> 5 </td>
 +
</tr>
 +
<th> val(;:;:3)= </th>
 +
</tr>
 +
            <tr>
 +
<td>  </td>
 +
<td> 0 </td>
 +
                                <td> 5 </td>
 +
<td> 0 </td>
 +
</tr>
 +
            <tr>
 +
<td>  </td>
 +
<td> 9 </td>
 +
                                <td> 12337 </td>
 +
<td> 4 </td>
 +
</tr>
 +
            <tr>
 +
<td>  </td>
 +
<td> 0 </td>
 +
                                <td> 10 </td>
 +
<td> 0 </td>
 +
</tr>
 +
</table>
 
 
 
 
 
  <p>When the α of MC = 0.0015, the α of LM = 0.0421:<br>
 
  <p>When the α of MC = 0.0015, the α of LM = 0.0421:<br>
 
The total number of Euk.Cement bound with sands in 30s =42164;<br>
 
The total number of Euk.Cement bound with sands in 30s =42164;<br>
 
The total number of Euk.Cement drifted away in 30s = 74661;<br></p>
 
The total number of Euk.Cement drifted away in 30s = 74661;<br></p>
−
<img class="picture2" src="https://static.igem.org/mediawiki/2015/0/0e/HUST_model2.5.png"><br><br>
+
<table border="0" style="width:300px" align:center>
 +
            <tr>
 +
<th> val(;:;:1)= </th>
 +
</tr>
 +
 
 +
<tr>
 +
<td>  </td>
 +
<td> 0 </td>
 +
                                <td> 2 </td>
 +
<td> 0 </td>
 +
</tr>
 +
<tr>
 +
<td>  </td>
 +
<td> 2 </td>
 +
<td> 12619 </td>
 +
<td> 2 </td>
 +
 
 +
</tr>
 +
<tr>    <td>  </td>
 +
<td> 0 </td>
 +
<td> 7 </td>
 +
<td> 0 </td>
 +
</tr>
 +
<tr>
 +
<th> val(;:;:2)= </th>
 +
</tr>
 +
            <tr>
 +
<td>  </td>
 +
<td> 7 </td>
 +
                                <td> 12638 </td>
 +
<td> 8 </td>
 +
</tr>
 +
            <tr>
 +
<td>  </td>
 +
<td> 12307 </td>
 +
                                <td> 1000000 </td>
 +
<td> 12394 </td>
 +
</tr>
 +
            <tr>
 +
<td>  </td>
 +
<td> 5 </td>
 +
                                <td> 12351 </td>
 +
<td> 6 </td>
 +
</tr>
 +
<th> val(;:;:3)= </th>
 +
</tr>
 +
            <tr>
 +
<td>  </td>
 +
<td> 0 </td>
 +
                                <td> 7 </td>
 +
<td> 0 </td>
 +
</tr>
 +
            <tr>
 +
<td>  </td>
 +
<td> 8 </td>
 +
                                <td> 12352 </td>
 +
<td> 9 </td>
 +
</tr>
 +
            <tr>
 +
<td>  </td>
 +
<td> 0 </td>
 +
                                <td> 4 </td>
 +
<td> 1 </td>
 +
</tr>
 +
</table>
  
−
             <h3 style="color:black" align="left"><b>Results:</b></h3><br>
+
             <h4 style="color:black" align="left"><b>Results:</b></h4><br>
 
             <p>From the result as shown above we could figure out that with 81.79% addition of the α, the bound Euk.Cement had an increase of 81.94%. This result indicates that the substance inside each lattice could be exchanged sufficiently to attain an even concentration of Euk.Cement as well as its secretions.<br>Therefore, we could change theα by selecting different Si-tag to achieve different targets, such as the intensity of the sands and the scope of the permeation. For example, we could change the scope Euk.Cement permeates and the intensity of sand by choosing different domain of Si-tag.<br><br>
 
             <p>From the result as shown above we could figure out that with 81.79% addition of the α, the bound Euk.Cement had an increase of 81.94%. This result indicates that the substance inside each lattice could be exchanged sufficiently to attain an even concentration of Euk.Cement as well as its secretions.<br>Therefore, we could change theα by selecting different Si-tag to achieve different targets, such as the intensity of the sands and the scope of the permeation. For example, we could change the scope Euk.Cement permeates and the intensity of sand by choosing different domain of Si-tag.<br><br>
 
<h3 style="color:black" align="left"><b>References:</b></h3><br>[1] Chen, W.W.; Neipel, M.; Sorger, P.K. (2010). "Classic and contemporary approaches to modeling biochemical reactions". Genes Dev 24 (17): 1861–1875.doi:10.1101/gad.1945410. PMC 2932968. PMID 20810646.<br>[2] Kenneth Connors, Chemical Kinetics, 1990, VCH Publishers, pg. 14<br>[3] Wolf-Gladrow D A. Lattice-gas cellular automata and lattice Boltzmann models: An Introduction[M]. Springer Science & Business Media, 2000.<br>[4] Wolf-Gladrow D A. Lattice-gas cellular automata and lattice Boltzmann models: An Introduction[M]. Springer Science & Business Media, 2000.<br>[5] Sorokina O, Kapus A, Terecskei K, et al. A switchable light-input, light-output system modelled and constructed in yeast[J]. J Biol Eng, 2009, 3: 15.<br>
 
<h3 style="color:black" align="left"><b>References:</b></h3><br>[1] Chen, W.W.; Neipel, M.; Sorger, P.K. (2010). "Classic and contemporary approaches to modeling biochemical reactions". Genes Dev 24 (17): 1861–1875.doi:10.1101/gad.1945410. PMC 2932968. PMID 20810646.<br>[2] Kenneth Connors, Chemical Kinetics, 1990, VCH Publishers, pg. 14<br>[3] Wolf-Gladrow D A. Lattice-gas cellular automata and lattice Boltzmann models: An Introduction[M]. Springer Science & Business Media, 2000.<br>[4] Wolf-Gladrow D A. Lattice-gas cellular automata and lattice Boltzmann models: An Introduction[M]. Springer Science & Business Media, 2000.<br>[5] Sorokina O, Kapus A, Terecskei K, et al. A switchable light-input, light-output system modelled and constructed in yeast[J]. J Biol Eng, 2009, 3: 15.<br>

Latest revision as of 15:50, 18 September 2015

Team:HUST-China:Modeling


click it~


Modeling on Ecosystem Level


With the benefit of the DDEs model we built, we could acquire necessary data to move on to the Ecosystem level modeling based on Cellular Automaton. In this modeling, we intended to foresee the important points of the application. For example, we could determine how should we release Euk.Cement, how would Euk.Cement spread and how much Euk.Cement do we need to produce the best result. The ecosystem could derive the best strategy to make full use of Euk.Cement.


Part one: the “wake-up” problem:


Since it takes some time for the darkness induction system to shut down, we have to treat our strains in darkness to make sure they recover from the shut-down state. However, the problem is: How long should we put the strains into darkness?

As we can see from the DDEs model, our system has bi-stable state: the stable state in darkness and another stable state in light. The longer our strains are exposed to light, the longer it takes to convert them to the stable-state of darkness. Besides, our strains are more sensitive to light than darkness. In another words, if our strains are exposed to light for 10 minutes, it has to take 60 to 70 minutes for the strains to return to the original state.

For storage and transportation, we have to keep our strains exposed to light all the time to ensure they would not produce Si-tag or Mcfp-3. When we need them to work, we have to turn off the light and put them into darkness and wait for them to “wake up”. With the DDEs model, we can simulate the process of how do the strains “wake up”.

Figure 6: Simulation of Darkness Induction System of Euk.Cement

Results:


The figure shows that it takes about 3600 minutes for Euk.Cement to “wake up” from the stable-state of light. Therefore, before releasing our kit into the sands, we should turn off the light and treat them for 60 hours. If some Euk.Cement drift away from the sands, they will be exposed to sunlight and their downstream systems would be shut down by darkness induction system.


Part two: the permeation problem:


After releasing Euk.Cement, it would permeate the sands and we intended to find out what’s the final outcome of its permeation and what strategy should we take to attain the best result. To achieve this goal, we have built a permeation model with Lattice Method (Cellular Automaton)[3] .

Hypothesis:


In order to make the model rigorous, we hypothesized the following prerequisites.
- There is no current in the sands.
- The total number of our yeast is constant. No yeast will be born or dead during the permeation process.
- The variability of temperature has no effect on our Euk.Cement or it could be ignored.
- The size of the sands and our Euk.Cement is constant.
- Every unit of sands has the same amount of silica.
The Cellular Automation simulated the process that our Euk.Cement expresses Si-tag and Mcfp-3 while it permeates in the sands. The PDEs (Partial Differential Equations) of this model are listed below.

Parameters:


The description of parameters, their values and the references involved in this model are listed in a table.

Brief Parameter Table of Diffusion Model

Parameter Description
a(x,y,z,t) the concentration of Mcfp-3
A(x,y,z,t) the concentration of fixed bacterial
B the vector of prerequisite coefficient
C(x,y,z,t) the concentration of free bacterial
bi the permeation coefficient of on the i direction
eb the rate of Mcfp-3 expression
ed the rate of Mcfp-3 degradation
l the edge length of each lattice
T0 the time unit of prerequisite
α the possibility of effective collision
λ average length of steps of Brownian motion
τ average free time of Brownian motion
NA the Avogadro constant
T the absolute temperature
R the Rydberg constant
r & d the radius and diameter of our yeast
ν the speed of the prerequisite
ρ the density of our yeast
η the dynamic viscosity coefficient
Ci the total number of bacterial of the 26 neighboring lattices of the lattice (x,y,z)
More detailed parameters table of diffusion model can be downloaded here.

Formulary:


Monte Carlo Simulation:


To determine the values of some necessary parameters in the permeation model, we run the Monte Carlo simulation for a single lattice.
The Euk.Cement moves in the sands and has a certain possibility to bind the silica. Besides, the motion of Euk.Cement includes the Brownian motion, the sedimentation due to gravity and the motion due to the current. We hypothesized that there is no water flow in the sands, and therefore the motion of Euk.Cement is largely depended on the former two.

Brownian motion:


(2.6)


We know that the speed of Brownian motion: is relevant to the value of free time. The average length of steps of Brownian motion has the same scale with the cracks of the sands. Hence, we can set the unit time of Monte Carlo Simulation as 1s (It’s not the unit time of permeation). And with this setting, we could derive that every yeast has one collision with the sands.

Sedimentation due to Gravity:


(2.7)


and therefore

(2.8)


Comparing the effect of sedimentation due to gravity and the Brownian motion, we reach a conclusion that the latter one outweigh the former one (about 100 times). And therefore the most effective factor of the motion of Euk.Cement in the sands is Brownian motion.

The code of the Monte Carlo Simulation can be downloaded here. Its flow chart is


After running the Monte Carlo Simulation with 106 yeasts for multiple times, we finally determine the values of the following parameters


Lattices Method:


With sufficient parameters, we could finally simulate the permeation process. The code of our permeation model could be downloaded here.

We chose to use numerical solution to solve our Permeation model. We divided the space into the same lattices, and each lattice has only one phase. Hence, each variable in the same coordinate is considered the same. The timeline was divided into the minimum and discrete values. By calculating the variables in each lattice through the discretized timeline, we could attain the values of a(x,y,z,t),A(x,y,z,t)and C(x,y,z,t). If the lattices were small enough, the PDEs above could be replaced by Lattices Method. The equations of Lattices Method (Cellular Automaton) are listed below.

With the equations of Cellular Automaton as well as Boundary equations and Initial state equations, we could attain the values of the concentration of Mcfp-3, the number of fixed Euk.Cement and free Euk.Cement at any time and at any coordinate.

We ran the simulation with a server and printed a heat map

Figure 8: Diffusion Model of Free Euk.Cement Cells
Figure 8: Diffusion Model of Bound Euk.Cement Cells
Figure 7: Diffusion Model of Secreted Mcfp-3 Protein

Results:


From the permeation process shown in the heat map, we can derive the following conclusion:
- The Euk.Cement permeates very efficiently over a large space. Therefore, perhaps we have to enclose it.
- A rock (1000 cm3) is added to the sands. However, it’s not even noticeable from the outcome. Hence, the rocks or stone inside the sands have little effect on the permeation.
- The final shape of the permeation is approximately global and not sensitive to the initial input.
- The Euk.Cement distributed quite evenly in the whole space they spread but not densely at the surface of the sands.
- The Euk.Cement expressed sufficient amount of Mcfp-3.


Robustness and Parameter Sensitivity Analysis:


The influential factors include noise of gene expression, variability of temperature, variability of the size of the sands and Euk.Cement, etc. However, our permeation model is most sensitive to the parameter α (the possibility of effective collision). Therefore, we focused on the analysis of the sensitivity of the parameter α.

When the α of MC = 0.001, the α of LM = 0.0282:
The total number of Euk.Cement bound with sands in 30s = 28224;
The total number of Euk.Cement drifted away in 30s = 74778;

val(;:;:1)=
0 8 0
4 12439 7
0 4 0
val(;:;:2)=
5 12474 4
12456 1000000 12659
2 12327 9
val(;:;:3)=
0 7 0
11 12423 4
0 4 0

When the α of MC = 0.0013, the α of LM = 0.0365:
The total number of Euk.Cement bound with sands in 30s = 36464;
The total number of Euk.Cement drifted away in 30s = 74925;

val(;:;:1)=
0 8 0
3 12538 9
0 7 0
val(;:;:2)=
5 12553 3
12549 1000000 12550
2 12398 5
val(;:;:3)=
0 5 0
9 12337 4
0 10 0

When the α of MC = 0.0015, the α of LM = 0.0421:
The total number of Euk.Cement bound with sands in 30s =42164;
The total number of Euk.Cement drifted away in 30s = 74661;

val(;:;:1)=
0 2 0
2 12619 2
0 7 0
val(;:;:2)=
7 12638 8
12307 1000000 12394
5 12351 6
val(;:;:3)=
0 7 0
8 12352 9
0 4 1

Results:


From the result as shown above we could figure out that with 81.79% addition of the α, the bound Euk.Cement had an increase of 81.94%. This result indicates that the substance inside each lattice could be exchanged sufficiently to attain an even concentration of Euk.Cement as well as its secretions.
Therefore, we could change theα by selecting different Si-tag to achieve different targets, such as the intensity of the sands and the scope of the permeation. For example, we could change the scope Euk.Cement permeates and the intensity of sand by choosing different domain of Si-tag.

References:


[1] Chen, W.W.; Neipel, M.; Sorger, P.K. (2010). "Classic and contemporary approaches to modeling biochemical reactions". Genes Dev 24 (17): 1861–1875.doi:10.1101/gad.1945410. PMC 2932968. PMID 20810646.
[2] Kenneth Connors, Chemical Kinetics, 1990, VCH Publishers, pg. 14
[3] Wolf-Gladrow D A. Lattice-gas cellular automata and lattice Boltzmann models: An Introduction[M]. Springer Science & Business Media, 2000.
[4] Wolf-Gladrow D A. Lattice-gas cellular automata and lattice Boltzmann models: An Introduction[M]. Springer Science & Business Media, 2000.
[5] Sorokina O, Kapus A, Terecskei K, et al. A switchable light-input, light-output system modelled and constructed in yeast[J]. J Biol Eng, 2009, 3: 15.