Biomod/2011/TeamJapan/Tokyo/Project/Simulations

From OpenWetWare
Revision as of 09:51, 2 November 2011 by Takafumi Hayashi (talk | contribs)
Jump to navigationJump to search


<html> <!--<head>--> <!-- We refer to the wiki of team iGEM Heidlberg 2009. -->

      <style type="text/css">

/* ====================

  主に全体に関わるCSS

   ==================== */

      body {behavior: url(http://www.xs4all.nl/~peterned/htc/csshover3-source.htc);}
      .clear {clear:both;}

body.mediawiki {

                       font-size: 14px;

background-color:#707070; background-position: center center; background-attachment: fixed; background-repeat: no-repeat; font-family: Calibri, Verdana, helvetica, sans-serif; }

               h1 {
                       padding: 0px 20px 5px 20px;
                       font-size: 34px;
                       font-weight: bold;
               }
               h2 {
                       padding: 20px 20px 5px 20px;
                       font-size: 25px;
                       color: #0083eb;
                       text-decoration: none;
                       font-weight: bold;}
                       border: none;
               h2 a {
                       color: #eb8300;
               }
               h3 {
                       padding: 20px 20px 5px 20px;
                       font-size: 20px;
                       color: #000;
                       font-decoration: none;
                       font-weight: bold;
               }
               h1.firstHeading {

display: none;

               }

p { text-align: justify; } a:link { color: #00a5ea; text-decoration: none } a:visited { color:#00a5ea; text-decoration: none } a:hover {

                       color: #eb8300;

text-decoration: none } a:active { color:#f29400; text-decoration: none } #bodyContent {

                       width: 970px;
                       margin: 0px 0px 0px;
                       background-color:#ffffff;
                       border-width: 0px 1px 0px 1px;
                       border-color: #000000;
               }
               #content {
                       padding-left: 0px; width: 970px;
               }

table#team_members {

                       text-align: justify;
                       border: 0;
               }

table#team_members h2, table#team_members h3 {

                       clear: both;
               }

#content * a:hover {

                       text-decoration:none;
               }

#main_wrapper {

                       position:absolute;
                       left:0px;
                       top:20px;

margin-top: 0; width: 969px; height: 221px; align: center; border-style: solid;

                       border-color: white;

} /* ====================

  メニューの画像を変更できる部分
  ==================== */

#header {

                       position:relative;
                       left:0px;
                       top:0px;

margin-top: 0; width: 969px; height: 221px; align: left; background-color: #FFFFFF;

                       background-image: url(http://openwetware.org/images/e/e0/NEW_header.jpg); 

} /* ====================

  以下、特殊なclassに適用される
  ==================== */

#navigation { position:absolute;

                       left:18px; top:155px; width:1200px; height:69px;

z-index:100;

                       background-color: transparent;
                       float: left; 

color: #0000FF; } #super_main_wrapper {

                       position:absolute;
                       left:0px;
                       top:227px;

width: 975px; align: center; background-color: #ede8e2;

                       heigth: auto;
               }
               #SubWrapper {
                       width: 645px;
                       padding: 0px;
                       border-left:4px solid #ede8e2;
                       float: left;
                       margin-top: 0px;
                       background-color: #ede8e2;
               }
               #SubWrapper * p, #SubWrapper p  {
                       padding: 0 20px;
                       text-align: justify;
                       font-size: 12px;
               }
               #SubWrapper * h3, #SubWrapper h3 {
                       padding-top: 10px;
                       font-size: 18px;
               }
               #news {
                       width: 322px;
                       margin-top: 0px;
                       float: left;
                       background-color: #d8d5d0;
                       border-right:4px solid #ede8e2;
               }
               #news p {
                       padding: 0 20px 20px 20px;
                       text-align: justify;
                       font-size: 12px;
               }
               #news h3 {
                       padding: 10px 20px;
                       font-size: 18px;
               }
               #mission_box {
                       width:650px;
                       float: left;
               }
               #team_box, #heartbeat_box, #notebook_box, #parts_box, #gallery_box, #sponsors_box, .boxy {
                       width:215px;
                       float: left;
                       padding: 10px 0 0 0;
               }
               div.tleft {
                       border-width: 0px;
                       margin:0;
                       padding:0;
                       border-color:transparent;
               }

/* ====================

  ここからプルダウン周辺のデザイン
  ==================== */
  1. menu * {

margin: 0; padding: 0; }

  1. menu {

behavior: url(http://www.xs4all.nl/~peterned/htc/csshover3-source.htc); <!-- fuer ie6 --> font-family: calibri, verdana, sans-serif;

       font-color: #ffffff;

font-size: 19px; background-color: transparent; float:left; padding: 12px 0 0 0; }

  1. menu ul {

float: left; list-style: none; }

  1. menu ul li {

background-color:transparent;

       position:relative;

float:left; list-style: none; padding: 10px 20px 0 0;

       font-weight: bold;
       width: auto;

}

  1. menu a {

color: #FFFFFF; display: inline; text-decoration: none; }

  1. menu a:visited {

color:#FFFFFF; text-decoration: none }

  1. menu a:hover {

color: #00a5ea; }

  1. menu ul li ul {

display: none; position: absolute; left: 0px;

       width: 155px;
       heigth: 1%;

font-size: 19px; opacity: 0.8; list-style: none;

       top: 30px;
       padding-top: 20px;
       z-index:500;

}

  1. menu ul li:hover ul {

display: inline;

       background-position: bottom;

}

  1. menu ul li ul li {

width: 100%; list-style: none;

       background-color: #000;
       margin: -1px;
       padding: 0px 0 0 5px;
       display: inline;

}

</style> <!-- </head> -->


<!-- ############################# BODY!!!!! #################################### -->

<!-- <body> --> <!--[if IE]> <style type="text/css">

  1. SubWrapper {width: 645px;}
  2. news {width: 317px;}
  3. super_main_wrapper {position:static;}
  4. navigation {left: 15px;}
  5. menu ul li ul {left: -2px; top: 22px;}

</style> <![endif]--> <div id="main_wrapper">

 <div id="header">

<div id="navigation"> <div id="menu" style="position:static"> <ul> <li><a class="aMain" href="http://openwetware.org/wiki/Biomod/2011/TeamJapan/Tokyo">Home</a></li> <li><a class="aTeam" href="http://openwetware.org/wiki/Biomod/2011/TeamJapan/Tokyo/Team/Students">Team</a></li> <li><a class="aProject" href="http://openwetware.org/wiki/Biomod/2011/TeamJapan/Tokyo/Project">Project</a> <!-- <ul> <li><a href="http://openwetware.org/wiki/Biomod/2011/TeamJapan/Tokyo/Project">Overview</a></li> <li><a href="http://openwetware.org/wiki/Biomod/2011/TeamJapan/Tokyo/Project/introduction">Introduction</a></li> <li><a href="http://openwetware.org/wiki/Biomod/2011/TeamJapan/Tokyo/Project/Model">Model</a></li> <li><a href="http://openwetware.org/wiki/Biomod/2011/TeamJapan/Tokyo/Project/Devices">Devices</a></li> <li><a href="http://openwetware.org/wiki/Biomod/2011/TeamJapan/Tokyo/Project/Modes">Modes</a></li> <li><a href="http://openwetware.org/wiki/Biomod/2011/TeamJapan/Tokyo/Project/Results">Results</a></li> <li><a href="http://openwetware.org/wiki/Biomod/2011/TeamJapan/Tokyo/Project/Achievements">Achievements</a></li> <li><a href="http://openwetware.org/wiki/Biomod/2011/TeamJapan/Tokyo/Project/Future_works">Future works</a></li> </ul> --> <li><font color="#ffffff">Results</font> <ul> <li><a href="http://openwetware.org/wiki/Biomod/2011/TeamJapan/Tokyo/Project/Results">Experiments</a></li> <li><a href="http://openwetware.org/wiki/Biomod/2011/TeamJapan/Tokyo/Project/Simulations">Simulations</a></li> <li><a href="http://openwetware.org/wiki/Biomod/2011/TeamJapan/Tokyo/Achievements/DNA_Devices">DNA Design</a></li> </ul></li> <!-- <li><a class="Simulation" href="http://openwetware.org/wiki/Biomod/2011/TeamJapan/Tokyo/Project/Simulations">Simulations</a></li> <li><a class="DNA design" href="http://openwetware.org/wiki/Biomod/2011/TeamJapan/Tokyo/Achievements/DNA_Devices">DNA Designs</a></li> --> <li><a href="http://openwetware.org/wiki/Biomod/2011/TeamJapan/Tokyo/Project/Achievements">Achievements</a></li> <li><a href="http://openwetware.org/wiki/Biomod/2011/TeamJapan/Tokyo/Project/Future_works">Future works</a></li> <li><a href="http://openwetware.org/wiki/Biomod/2011/TeamJapan/Tokyo/Notebook/Protocols">Protocols</a></li> <li><a href="http://openwetware.org/wiki/Biomod/2011/TeamJapan/Tokyo/Notebook/Lab.notebook">Notes</a></li> <li><a class="aNotebook" href="http://openwetware.org/wiki/Biomod/2011/TeamJapan/Tokyo/Sponsors/">Sponsors</a></li> <li><a class="aSitemap" href="http://openwetware.org/wiki/Biomod/2011/TeamJapan/Tokyo/Sitemap">Sitemap</a></li>

</ul> </div><!-- end drop menu --> </div> <!-- end navigation --> </div> </div> <!-- For IE 6 --> <script type="text/javascript">

/***********************************************

  • CSS Vertical List Menu- by JavaScript Kit (www.javascriptkit.com)
  • Menu interface credits: http://www.dynamicdrive.com/style/csslibrary/item/glossy-vertical-menu/
  • This notice must stay intact for usage
  • Visit JavaScript Kit at http://www.javascriptkit.com/ for this script and 100s more
                                                                                              • /

var menuids=new Array("verticalmenu") //Enter id(s) of UL menus, separated by commas var submenuoffset=-2 //Offset of submenus from main menu. Default is -2 pixels.

function createcssmenu(){ for (var i=0; i<menuids.length; i++){

 var ultags=document.getElementById(menuids[i]).getElementsByTagName("ul")
   for (var t=0; t<ultags.length; t++){
   var spanref=document.createElement("span")

spanref.className="arrowdiv" spanref.innerHTML="&nbsp;&nbsp;" ultags[t].parentNode.getElementsByTagName("a")[0].appendChild(spanref)

   ultags[t].parentNode.onmouseover=function(){
   this.getElementsByTagName("ul")[0].style.left=this.parentNode.offsetWidth+submenuoffset+"px"
   this.getElementsByTagName("ul")[0].style.display="block"
   }
   ultags[t].parentNode.onmouseout=function(){
   this.getElementsByTagName("ul")[0].style.display="none"
   }
   }
 }

} if (window.addEventListener) window.addEventListener("load", createcssmenu, false) else if (window.attachEvent) window.attachEvent("onload", createcssmenu) </script> <!-- </body> --> </html>













Simulations for track walking mode

Principles and methods of simulations

To confirm this mode, we simulated whether DNA ciliate moves intended direction on substrate DNAs track by cellular automaton framework. Cellular automaton framework is discrete model. Its field is assumed that the width is infinite and cells are square-block type. Next time step state is calculated by present own value and neighbor’s value. In this simulation, probability for movement is determined by each cell’s free energy. Probability for moving to substrate is larger than moving to cleaved substrate because substrate’s energy is smaller than cleaved substrate’s energy.
This program is the simulation of the behavior of DNA ciliate, which works base on cellular-automaton framework.

Method

In our cellular-automaton simulations, we defined cells as followings.
  • Length of side of a cell is the same to the diameter of DNA ciliate body.
  • A yellow cell represents DNA ciliate.
  • Blue cells represent oligodeoxynucleotide substrates with a single ribose moiety.
  • Red cells represent oligodeoxynucleotide substrates without a ribose moiety.
  • Green cells represent UV-switching DNAs which are spotted UV.
  • A Yellow cell makes a move depend on its located cell. The rate of moving on each cell is derived as following steps.
  • Yellow cells don’t move to the yellow cell. (This means DNA ciliate )


We calculated the average time that DNA ciliate moves one micrometer, the rate that yellow cells don’t move at the present cell, and the rate that yellow cells moves to each cell. The detail of the calculations is following.


Firstly, we calculated the average time that DNA ciliate moves one micrometer by Brownian movement as following formula.
R: Gas constant. 8.3145[P/mol*K]
T: Temperature. 298[K].
Na=Avogadro number. 6.022141×1023
η: Viscosity. 0.00089[Pa/s]
a: Radius of DNA ciliate body.[m]
x: Length of side of a cell. [m]


Secondly, we calculated the rate that yellow cells don’t move based on two assumptions. One assumption is that a DNA ciliate moves to another cell with cleaving substrate in 120 seconds at a rate of 0.5. In the other assumption, a DNA ciliate moves to another cell without cleaving substrate in 1.2 seconds at a rate of 0.5. This assumption is based on this thesis[1]. In this thesis, the time which DNA releases substrate by using enzyme activity is slower about 100 times than the time which DNA releases substrate by not using enzyme activity. We defined the “step” as an indicator of time. The “step” is the unit time. The step of DNA ciliate which is 1um in diameter is about 1.02 seconds. The step of DNA ciliate which is 200nm in diameter is about 0.00815 seconds.
t: Average time to move one cell. [s/μm]
y: Steps to move at rate 0.5. [s]
Then, each rate is following.
The diameter of DNA ciliate body is 1 micro-meter.
  On the blue cells, yellow cells do not move at a rate of 0.9941.
  On the red cells, yellow cells do not move at a rate of 0.555.
  On the white cells, yellow cells do not move at a rate of 0 (means to move absolutely).
 
The diameter of DNA ciliate body is 200 nanometer.
  On the blue cells, yellow cells do not move at a rate of 0.9999529.
  On the red cells, yellow cells do not move at a rate of 0.9953006.
  On the white cells, yellow cells do not move at a rate of 0 (means to move absolutely).


Thirdly, we calculated the rate that yellow cells move to each cell.
  • Yellow cells can move to adjacent cells at one step and cannot move to distant cells.
  • The direction of yellow cells' movement is stochastically dependent on free-energy of the adjacent cells. The free-energy of each cell is followings. Free-energy of blue cells is -19.02[KJ/mol]. Free-energy of red cells is -16.03[KJ/mol]. Free-energy of green cells is -49.90[KJ/mol]. For examples, Yellow cells are about 211 times more likely to move to blue cells than to red cells.


ΔGi: An free energy of cell i which is adjacent;
R: Gas constant. 8.3145[P/mol*K]
T: Temperature. 298[K].
Based on these three calculations, we simulated the movement of DNA ciliates.

Results

We simulated movement of DNA ciliate

Simulation movie of track walking mode (a single DNA ciliate)

We simulated the movement of DNA ciliate on DNA track. In the simulation, the width of the DNA track was 4um, a DNA ciliate was 1um in diameter, and the length of the DNA track was 50um.

<html><body><center><iframe width="450" height="259" src="http://www.youtube.com/embed/jbS0hgjK7q8?rel=0" frameborder="0" allowfullscreen></iframe></center></body></html>

Execute a Simulation
In this simulation, we confirmed that DNA ciliate moved from left to right directly.


To confirm “tracks walking mode” works correctly, we made line graphs which means the migration length of DNA ciliate. This graph is the average of the movement by several simulations. The vertical axis shows the distance from origin. The horizontal axis shows the number of steps. Three lines are different in the DNA types of DNA track. The red line's DNA track issubstrate, so DNA ciliate moves with cleaving substrate. The blue line's DNA track is cleaved substrate, so DNA ciliate moves without cleaving in this situation. The green line's DNA track is complementary strands for deoxyribozyme, so DNA ciliate can hybridize with, but cannot cleave the DNA. Figure1 is the graph of the movement of DNA ciliate which is 1um, and figure2 is the graph of the movement of DNA ciliate which is 200nm.
Figure1
Figure1:The result of the simulation. (The DNA ciliate is 1um in diameter)
Figure2
Figure2:The result of the simulation. (The DNA ciliate is 200nm in diameter)
From these graphs, we found that DNA ciliate moving on substrate moves faster and more directly than the others.


To confirm DNA ciliate is forced motion, we make mean square displacement graph. In mean square displacement, Brownian motion is presented by linear function and forced motion is presented by quadratic function.
Figure3
Figure3:The result of the simulation.
In this simulation, the moving distance on substrate field is presented by quadratic function. We confirmed DNA ciliate is forced directly and moves intended direction on substrate in this simulation. Furthermore, DNA ciliate on DNA track is faster than on cleaved substrate.


Simulation movie of DNA ciliates on DNA track (several DNA ciliates)

Actually, it is difficult to put a single DNA ciliate on DNA track, so it is necessary to confirm DNA ciliates can moves directly on DNA track if several DNA ciliates are put on at a time. We simulated the movement of several DNA ciliates on DNA track. The number of DNA ciliates is 12, the width of DNA track is 30 um, and the length of DNA track is 100 um.

<html><body><center><iframe width="450" height="259" src="http://www.youtube.com/embed/EQoU17FJeLM?hl=ja&fs=1" frameborder="0" allowfullscreen></iframe></center></body></html>

From the result, most of DNA ciliates move to right directly. It is confirmed several DNA ciliates moves directly on a DNA track if several DNA ciliates are put on at a time.


Simulation of light-irradiate mode on DNA track

We simulated the movement of several DNA ciliates on DNA track which contains UV-switching DNA. The UV-switching DNA is spotted UV, so DNA ciliates are trapped by the DNA. The left video is not attached UV-switching DNA, on the other hand, the right video is attached UV-switching DNA. The number of DNA ciliates is 12, the width of DNA track is 30 um, and the length of DNA track is 100 um.

<html><body><center><iframe width="450" height="259" src="http://www.youtube.com/embed/4BE5f5qPobo?hl=ja&fs=1" frameborder="0" allowfullscreen></iframe></center></body></html>

From the result, most of DNA ciliates gathered at the area of UV-switching DNA, so it is confirmed DNA ciliate on DNA track gather at UV-switching DNA.
[Program Source]

Discussion

We simulated the movement of a DNA ciliate on a long DNA track. The following graph is the mean square displacement graph for movement of a DNA ciliate. This graph is same to figure3, except in the length of the x axis.
Figure4
Figure4:The result of the simulation.
In this graph, the red line (on DNA track) changes quadratic function to linear function little by little. It is thought that the movement of a DNA ciliate which walk on DNA track change track walking mode to free moving mode.
To consider about this incomprehensible, we made line graphs which mean the movement a DNA ciliate. Each line of these graph means a single simulation, NOT the average of several simulations. The left graph is mean square displacement of a DNA ciliate on DNA track. The right graph is the movement of a DNA ciliate on cleaved substrate.
Figure5
Figure5:The result of the simulation.
In left graph, all lines are very rough and approximated linear functions, so it is thought DNA ciliate moves at ramdom.
In right graph, most of lines are approximated quadratic functions at first, but some of lines become rough and are approximated linear functions. It is thought DNA ciliate on DNA track walk on track at first, but a short time later, some DNA ciliates become moving on cleaved substrate. It is thought that the number of the DNA ciliates which move on cleaved substrate become large little by little because the probability of moving to cleaved substrate and starting free moving become large according to the area of cleaved substrate spreads.
To reduce the probability of moving to cleaved substrate and starting free moving, it is thought to reduce the free-energy of substrate. However, it makes the time that DNA ciliate separate from substrate or cleaved substrate become long, so the speed of moving will be slow. The speed and the distance of keep walking on DNA track are trade-off.


Conclusion

In simulations, we confirmed that DNA ciliate moves by directional force, not by Brownian motion in “track walking mode”. Furthermore, we found the speed of movement on DNA track is faster than on cleaved substrate.


Program sources

Simulation for no light-irradiated DNA and single DNA ciliate

[Program Source]

Simulation for no light-irradiated DNA and several DNA ciliates

[Program Source]

Simulation for light-irradiated DNA and several DNA ciliates

[Program Source]


References

[1] Multivalent random walkers: a model for deoxyribozyme walkers