Biomod/2012/TeamSendai/Simulation

From OpenWetWare

< Biomod | 2012(Difference between revisions)
Jump to: navigation, search
Current revision (23:06, 27 October 2012) (view source)
 
(47 intermediate revisions not shown.)
Line 1: Line 1:
<html>
<html>
-
<head>
 
-
    <title>Team Sendai Top</title>
 
-
    <script type="text/javascript" src="http://code.jquery.com/jquery-latest.min.js"></script>
 
-
    <script type="text/javascript">
 
-
        $(function() {
 
-
          if ($.browser.msie && $.browser.version.substr(0,1)<7)
 
-
          {
 
-
$('li').has('ul').mouseover(function(){
 
-
$(this).children('ul').show();
 
-
}).mouseout(function(){
 
-
$(this).children('ul').hide();
 
-
})
 
-
          }
 
-
        });       
 
-
    </script>
 
-
 
-
<style>
 
-
#Container
 
-
{
 
-
width: 960px;
 
-
margin: 40px auto;
 
-
}
 
-
/* Main menu */
+
<head><style>#column-one { display:none; width:0px;}
-
#menu
+
.container{background-color: #f5f5f5; margin-top:50px}
-
{
+
.OWWNBcpCurrentDateFilled {display: none;}
-
width: 100%;
+
#content {width: 0px; margin: 0 auto auto 0; padding: 1em 1em 1em 1em; align: center;}
-
margin: 0;
+
#column-content {width: 0px; float: left; margin: 0 0 0 0;padding: 0;}
-
padding: 10px 0 0 0;
+
.firstHeading {display:none; width:0px;}
-
list-style: none;
+
#globalWrapper{width:1280px; margin:auto}
-
background: #111;
+
body {background: #F0F0F0 !important;}
-
background: -moz-linear-gradient(#444, #111);  
+
#column-one {display:none; width:0px;background-color: #f0f0f0;}
-
    background: -webkit-gradient(linear,left bottom,left top,color-stop(0, #111),color-stop(1, #444));
+
#content{border:none;margin: 0 0 0 0; padding: 1em 1em 1em 1em; position: center; width: 800px;background-color: #f0f0f0; }
-
background: -webkit-linear-gradient(#444, #111);
+
.container{ width:1280px; margin: auto; background-color: #f0f0f0; text-align:justify; font-family: helvetica, arial, sans-serif; color:#f0f0f0; margin-top:25px; }
-
background: -o-linear-gradient(#444, #111);
+
#bodyContent{ width: 1267px; align: center; background-color: #f0f0f0;}
-
background: -ms-linear-gradient(#444, #111);
+
#column-content{width: 1280px;background-color: #f0f0f0;}
-
background: linear-gradient(#444, #111);
+
.firstHeading { display:none;width:0px;background-color: #f0f0f0;}
-
-moz-border-radius: 50px;
+
#header{position: center; width: 800px;background-color: #f0f0f0;}
-
border-radius: 50px;
+
#footer{position: center; width:1280px;}
-
-moz-box-shadow: 0 2px 1px #9c9c9c;
+
</style><header>
-
-webkit-box-shadow: 0 2px 1px #9c9c9c;
+
-
box-shadow: 0 2px 1px #9c9c9c;
+
-
}
+
-
#menu li
+
</html>
-
{
+
-
float: left;
+
-
padding: 0 0 10px 0;
+
-
position: relative;
+
-
line-height: 0;
+
-
}
+
-
#menu a
 
-
{
 
-
float: left;
 
-
height: 25px;
 
-
padding: 0 25px;
 
-
color: #999;
 
-
text-transform: uppercase;
 
-
font: bold 12px/25px Arial, Helvetica;
 
-
text-decoration: none;
 
-
text-shadow: 0 1px 0 #000;
 
-
}
 
-
#menu li:hover > a
 
-
{
 
-
color: #fafafa;
 
-
}
 
-
*html #menu li a:hover /* IE6 */
 
-
{
 
-
color: #fafafa;
 
-
}
 
-
#menu li:hover > ul
 
-
{
 
-
display: block;
 
-
}
 
-
/* Sub-menu */
+
{{:Biomod/2012/Tohoku/Team_Sendai/header}}
-
#menu ul
 
-
{
 
-
    list-style: none;
 
-
    margin: 0;
 
-
    padding: 0;   
 
-
    display: none;
 
-
    position: absolute;
 
-
    top: 35px;
 
-
    left: 0;
 
-
    z-index: 99999;   
 
-
    background: #444;
 
-
    background: -moz-linear-gradient(#444, #111);
 
-
    background: -webkit-gradient(linear,left bottom,left top,color-stop(0, #111),color-stop(1, #444));
 
-
    background: -webkit-linear-gradient(#444, #111);   
 
-
    background: -o-linear-gradient(#444, #111);
 
-
    background: -ms-linear-gradient(#444, #111);
 
-
    background: linear-gradient(#444, #111);
 
-
    -moz-box-shadow: 0 0 2px rgba(255,255,255,.5);
 
-
    -webkit-box-shadow: 0 0 2px rgba(255,255,255,.5);
 
-
    box-shadow: 0 0 2px rgba(255,255,255,.5);
 
-
    -moz-border-radius: 5px;
 
-
    border-radius: 5px;
 
-
}
 
-
#menu ul ul
+
<div id="Container">
-
{
+
<h1>Aim of Simulation</h1>
-
  top: 0;
+
Molecular robots works in meso-scale (from nm to sub-micro meter). In the meso-scale, characters of each molecule and cooperative movement of molecules could not be ignored. Without calculation of potentials or molecular dynamics simulations, it is difficult to predict the behaviors of the molecules in this scale. Thus, we calculated potentials to confirm the GATE effect, and simulated molecular dynamics of DNAs to strength reliability of our design.
-
  left: 150px;
+
-
}
+
-
#menu ul li
+
<h1>Numerical Analysis of Electrostatic Potential</h1>
-
{
+
Electric potential of the GATE was calculated inside and outside the GATE. {{-}}The phosphate groups in the backbone of DNA are negatively-charged. Because the GATE is made by DNA origami, the GATE has the Coulomb force.
-
    float: none;
+
-
    margin: 0;
+
-
    padding: 0;
+
-
    display: block; 
+
-
    -moz-box-shadow: 0 1px 0 #111111, 0 2px 0 #777777;
+
-
    -webkit-box-shadow: 0 1px 0 #111111, 0 2px 0 #777777;
+
-
    box-shadow: 0 1px 0 #111111, 0 2px 0 #777777;
+
-
}
+
-
#menu ul li:last-child
+
<h2>Model of DNA</h2>
-
{  
+
A point-charge model is applied to DNA model. Each phosphate of a nucleotide was converted into one electric beads, which has elementary electoric charge (1.602176565 x 10^−19 [C]). In DNA helix, negative charges appear along the axis of the double helix every 10.4 base pairs. {{-}}Debye–Hückel equation was used to calculate electric potential. Electric potential of each point around the GATE was deduced from summation of all the Electric Potential for every phosphate presented on DNA.{{-}}
-
    -moz-box-shadow: none;
+
-
    -webkit-box-shadow: none;
+
-
    box-shadow: none;   
+
-
}
+
-
#menu ul a
 
-
{   
 
-
    padding: 10px;
 
-
height: 10px;
 
-
width: 130px;
 
-
height: auto;
 
-
    line-height: 1;
 
-
    display: block;
 
-
    white-space: nowrap;
 
-
    float: none;
 
-
text-transform: none;
 
-
}
 
-
*html #menu ul a /* IE6 */
+
Debye–Hückel equation:
-
{  
+
{{-}}
-
height: 10px;
+
[[Image: Potential_fomula.png |left|300px]]
-
}
+
{{-}}
 +
reference: Debye length
 +
{{-}}
 +
[[Image: Debyelen.png |left|300px]]
 +
{{-}}
-
*:first-child+html #menu ul a /* IE7 */
 
-
{   
 
-
height: 10px;
 
-
}
 
-
#menu ul a:hover
+
Following GIF animation shows the scheme of a point-charge model and potential calculations.{{-}}
-
{
+
(It takes about 7 seconds to watch one cycle of the GIF animation){{-}}
-
    background: #0186ba;
+
{{-}}
-
background: -moz-linear-gradient(#04acec,  #0186ba);
+
{{-}}
-
background: -webkit-gradient(linear, left top, left bottom, from(#04acec), to(#0186ba));
+
[[Image: Helix.gif |left|500px|thumb]]
-
background: -webkit-linear-gradient(#04acec,  #0186ba);
+
{{-}}
-
background: -o-linear-gradient(#04acec,  #0186ba);
+
-
background: -ms-linear-gradient(#04acec,  #0186ba);
+
-
background: linear-gradient(#04acec,  #0186ba);
+
-
}
+
-
#menu ul li:first-child > a
 
-
{
 
-
    -moz-border-radius: 5px 5px 0 0;
 
-
    border-radius: 5px 5px 0 0;
 
-
}
 
-
#menu ul li:first-child > a:after
+
<h2>Results</h2>
-
{
+
Electric potential was calculated by using C language.{{-}}
-
    content: '';
+
Under our conditions, the target DNA has 25 electric beads.{{-}}
-
    position: absolute;
+
All calculation were performed under these condition: Temperature 298[K], Na+ 50mM.{{-}}
-
    left: 30px;
+
-
    top: -8px;
+
-
    width: 0;
+
-
    height: 0;
+
-
    border-left: 5px solid transparent;
+
-
    border-right: 5px solid transparent;
+
-
    border-bottom: 8px solid #444;
+
-
}
+
-
#menu ul ul li:first-child a:after
+
<h3>Electric potential from the top of view of the GATE</h3>
-
{
+
The following figure shows electric potentials between the target DNA and the GATE from the top of view of the GATE.{{-}}
-
    left: -8px;
+
The electric potentials were standardized by thermal energy, <I>K<sub>B</sub>T</I> (K<sub>B</sub>; Boltzmann constant).{{-}}
-
    top: 12px;
+
Potential energies were shown by a heat map.{{-}}
-
    width: 0;
+
-
    height: 0;
+
-
    border-left: 0;
+
-
    border-bottom: 5px solid transparent;
+
-
    border-top: 5px solid transparent;
+
-
    border-right: 8px solid #444;
+
-
}
+
-
#menu ul li:first-child a:hover:after
+
[[Image: Potential2.png |center|400px| Electric potential from the top of view of the GATE ]]
-
{
+
{{-}}
-
    border-bottom-color: #04acec;
+
{{-}}
-
}
+
-
#menu ul ul li:first-child a:hover:after
+
Potential energy along axis indicated by the arrows are shown as the following figure.{{-}}
-
{
+
Together, electric potential is too high for the target DNA to enter inside the GATE. {{-}}
-
    border-right-color: #04acec;
+
-
    border-bottom-color: transparent;
+
-
}
+
 +
{{-}}
 +
[[Image: Potential3.png |center|600px| Electric potential from the top of view of the GATE along an axis]]
 +
{{-}}
 +
{{-}}
-
#menu ul li:last-child > a
+
<h3>Electric potential at the center of hole in the GATE</h3>
-
{
+
-
    -moz-border-radius: 0 0 5px 5px;
+
-
    border-radius: 0 0 5px 5px;
+
-
}
+
-
/* Clear floated elements */
+
The difficulty to go through inside hole of the GATE was also indicated by following figures.{{-}}
-
#menu:after
+
The next figures shows electric potential at the center of hole in the GATE.{{-}}
-
{
+
Potentials were calculated along the axis indicated by the red arrow.{{-}}
-
visibility: hidden;
+
Blue lines shows electric potential in the case of 1.5 fold radius of the GATE.{{-}}
-
display: block;
+
-
font-size: 0;
+
-
content: " ";
+
-
clear: both;
+
-
height: 0;
+
-
}
+
-
* html #menu            { zoom: 1; } /* IE6 */
+
{{-}}
-
*:first-child+html #menu { zoom: 1; } /* IE7 */
+
[[Image: Potential_K.png |center|800px| Electric potential at the center of hole in the GATE ]]
 +
{{-}}
 +
{{-}}
 +
By the figure, we concluded {{-}}
 +
i ) Our designed size is suitable size for the GATE function{{-}}
 +
ii) Enlarging radius of the GATE decreases electric repulsion effect{{-}}
-
/*目次*/
 
-
div#mokuji {width: 800px;
 
-
margin-left: auto;
 
-
margin-right: auto;
 
-
background-color: #f5f5dc}
 
-
div#mokuji h2 { background-color: # f5f5dc;
+
<h3>Hybridization energy of Porters overcomes electric potential of the GATE</h3>
-
font-size: 1.50em;
+
-
color: #000000;
+
-
line-height: 45px;
+
-
padding-left: 12px;
+
-
margin-bottom: 0}
+
-
ol#mokuji {font-size: 1.00em;
+
We expected the hybridization energy of the Porters overcomes the inhibition of entrance of the GATE by electric repulsion.{{-}}
-
margin-left: 0;
+
The Gibbs energies of hybridization was calculated by neighbor-joining methods.{{-}}{{-}}
-
padding-left: 0}
+
The result of Gibbs energy of the fixed Porter were defined as the maximum length of each form of Porter. For examples, Porter 1 has three forms: Stretched form, 1 loop form, and 2 loop form. In this Porter 1 case, we defined Porter 1 has three Gibbs energy according to the hybridization state. Since the maximum lengths of each form were different, the minimum energy of possible length was selected.{{-}} Gibbs energies were calculated at the center of hole in the GATE.{{-}}
 +
The following figure shows the calculated energy of Porters and Toehold DNA used as controls in our project.
 +
{{-}}
 +
[[Image: Porter1.png |center|500px| Three different form ]]
 +
{{-}}
 +
{{-}}
 +
[[Image: Gibbs kai.png |center|800px| Gibbs energies by the Porters in the GATE ]]
 +
{{-}}
 +
{{-}}
 +
The Gibbs energies by the Porters shown above were summed with electric potential of the GATE.{{-}}
 +
The summation of energies were calculated at the center of hole in the GATE.{{-}}
 +
The result was shown in the next figure.
 +
{{-}}
 +
{{-}}
 +
[[Image: Gibbs2.png |center|800px| Electric potential from the top of view of the GATE ]]
 +
{{-}}
-
/* コンテンツ */
+
The results of calculation indicate that the outside Porter can catch the target DNA, and inner Porter can pull in the GATE step by step. On the other hand, short DNA like the toehold DNA we called n this study cannot catch the target DNA. {{-}}
-
div#Content {width: 800px;
+
-
margin-left: auto;
+
-
margin-right: auto}
+
-
div#Content h2 { background-color: # f5f5dc;
 
-
font-size: 1.50em;
 
-
color: #000000;
 
-
line-height: 10px;
 
-
padding-left: 12px;
 
-
margin-top: 20px;
 
-
margin-bottom: 0}
 
-
div#Content h3 {background-color: #ffffff;
 
-
background-image: url(maru-skyblue.png);
 
-
background-repeat: no-repeat;
 
-
background-position: 0px 2px;
 
-
font-size: 0.875em;
 
-
line-height: 22px;
 
-
padding-left: 26px;
 
-
margin-top: 30px;
 
-
margin-bottom: 0;
 
-
margin-left: 12px;
 
-
margin-right: 12px}
 
-
div#Content p {font-size: 1.25em;
+
<h1>MD Simulation</h1>
-
line-height: 1.6;
+
To further strengthen feasibility of our design, we examined the capturing mechanism and
-
margin-top: 10px;
+
the effectiveness of our structure “CELL GATE" by molecular dynamics (MD) simulation.
-
margin-left: 12px;
+
-
margin-right: 12px}
+
-
</style>
+
-
</head>
+
-
<body>
 
-
<div id="Container">
 
-
   
 
-
<!-- Menu -->
 
-
<ul id="menu">
 
-
<li><a href="http://openwetware.org/wiki/Biomod/2012/Tohoku/Team_Sendai ">Top</a></li>
 
-
<li><a href=" http://openwetware.org/wiki/Biomod/2012/TeamSendai/Idea ">Project</a></li>
 
-
<li><a href=" http://openwetware.org/wiki/Biomod/2012/TeamSendai/Simulation">Simulation</a>
 
-
</li>
 
-
<li><a href=" http://openwetware.org/wiki/Biomod/2012/TeamSendai/Design">Design</a>
 
-
</li>
 
-
<li>
 
-
<a href=" http://openwetware.org/wiki/Biomod/2012/TeamSendai/Experiment ">Experiment</a>
 
-
<ul>
 
-
<li><a href=" http://openwetware.org/wiki/Biomod/2012/TeamSendai/Method">Method</a>
 
-
<ul>
 
-
<li>
 
-
<a href=" http://openwetware.org/wiki/Biomod/2012/TeamSendai/Result#Porter">Porter</a>
 
-
<li>
 
-
<a href=" http://openwetware.org/wiki/Biomod/2012/TeamSendai/Result#Cylinder">Cylinder</a> </li>
 
-
<li> <a href=" http://openwetware.org/wiki/Biomod/2012/TeamSendai/Result# Vesicle">Vesicle</a>
 
-
</li>
 
-
</ul>
 
-
</li>
 
-
<li>
 
-
<a href=" http://openwetware.org/wiki/Biomod/2012/TeamSendai/Result">Result</a>
 
-
<ul>
 
-
<li>
 
-
<a href=" http://openwetware.org/wiki/Biomod/2012/TeamSendai/Result#Porter">Porter</a>
 
-
<li>
 
-
<a href=" http://openwetware.org/wiki/Biomod/2012/TeamSendai/Result#Cylinder">Cylinder</a> </li>
 
-
<li> <a href=" http://openwetware.org/wiki/Biomod/2012/TeamSendai/Result# Vesicle">Vesicle</a>
 
-
</li>
 
-
</ul>
 
-
</li>
 
-
</ul>
 
-
</li>
 
-
<li>
 
-
<a href=" http://openwetware.org/wiki/Biomod/2012/TeamSendai/Achievement">Achievement</a>
 
-
</li>
 
-
<li>
 
-
<a href=" http://openwetware.org/wiki/Biomod/2012/TeamSendai/Diary">Diary</a>
 
-
</li>
 
-
<li>
 
-
<a href=" http://openwetware.org/wiki/Biomod/2012/TeamSendai/Team ">Team</a>
 
-
</li>
 
-
<li>
 
-
<a href=" http://openwetware.org/wiki/Biomod/2012/TeamSendai/FAQ">FAQ</a>
 
-
</li>
 
-
</ul>
 
-
<!--目次 -->
+
<h2>DNA Model</h2>
-
<div id="mokuji">
+
[[Image: Format_DNA_rasen.jpg |center|638px]]
-
<h2>Contents</h2>
+
For simplicity, a course-grained DNA model was used in our simulation. One DNA nucleotide was represented by one bead in the model and each bead can be hybridized with a complementary bead.
-
<ol>
+
-
<li><a href="#Numerical Calculation for Electric Potential">Electric Potential Numerical Calculation</a></li>
+
-
<ol>
+
-
<li><a href="#Model">Model</a></li>
+
-
<li><a href="#">Results</a></li>
+
-
</ol>
+
-
<li><a href="#MD Simulation">MD Simulation</a></li>
+
The potential energy of the system includes 5 distinct contributions.
-
<ol>
+
-
<li><a href="#DNA Model">DNA Model</a></li>
+
-
<ol>
+
-
<li><a href="#Results">Results</a></li>
+
-
</ol>
+
-
<li><a href="#Comparison of capture ability">Comparison of capture ability</a></li>
+
[[Image: Format_suusiki1.jpg |left|711px]]
-
<ol>
+
{{-}}
-
<li><a href="#Results">Results</a></li>
+
-
</ol>
+
-
<li><a href="#Reference">Reference</a></li>
+
-
</ol>
+
-
</ol>
+
The first three terms are intramolecular interactions, bonds, bond angles, and dihedral angles. These terms are important for maintaining structure. In order to express the “tether like structure”, only bond interactions are considered in our DNA model.
-
</div>
+
-
<p>
+
And the latter two terms are non-bonded interactions. Coulomb interactions are taken into account using the Debye-Huckel approximation which enables to internalize  counter-ions contribution.
-
<br><br>
+
 +
Parameters of these potentials were fit to the reference literature ; Thomas A. Knotts ''et al''.  A coarse grain model of DNA .
 +
The force on bead i is given by a Langevin equation
-
</p>
+
[[Image: Format_suusiki2.jpg |left|357px]]
-
<a name="Numerical Calculation"></a><h2>Numerical Calculation</h2>
+
{{-}}
-
<p>
+
-
A phosphodiester bond make up the backbone of each helical strand of DNA. <br>
+
The first term donates a conservative force derived from the potential Vtot and the second is a viscosity dependent friction.
-
The phosphate groups in the phosphodiester bond are negatively-charged.<br>
+
-
Because gate is produced by DNA, we can not ignore the influence of the Coulomb force.<br>
+
-
So we calculate the electric potential near the gate.
+
-
</p>
+
The third term is a white Gaussian noise and effects of collision with solvent molecules which causes brownian motion are internalized in this term.
-
<a name="Model"></a><h2>Model</h2>
+
Langevin equation is integrated using a Velocity-Verlet method.
-
<p>
+
-
<br>
+
-
Sets the coordinates as follows.<br>
+
-
<img src="http://openwetware.org/images/9/90/Cy.png" width="350px" height="300px">
 
-
<img src="http://openwetware.org/images/6/66/Lin.jpg" width="420px" height="300px"><br>
 
-
<br><br><br>
 
-
Point-charge model is used.<br>
+
<h2> Verification of the Coarse-grained model: Toehold displacement of dsDNA</h2>
-
Assumesd the phosphate groups negative charge,and<br>
+
[[Image: Format_toea_toeb.jpg |right|400px]]
-
negative charge circles the axis of the double helix once every 10.4 base pairs like DNA.<br>
+
In order to test the model, here we carried out a simulation of Toehold displacement between two strands.
-
And we use follow fomula to calculate electric potential.<br><br>
+
Toehold strands are modeled according to designs of [http://openwetware.org/wiki/Biomod/2012/TeamSendai/Experiment#Comparison_of_Porter_and_Toehold the experiment section].
-
Debye–Hückel equation<br>
+
-
<img src="http://openwetware.org/images/e/ec/Potential_fomula.png" width="300px" height="90px"><br>
+
-
<br>
+
-
Debye length<br>
+
-
<img src="http://openwetware.org/images/f/f9/Debyelen.png" width="400px" height="250px"><br><br>
+
 +
Length of strands and simulation condition were as follows:
 +
Target strand / Toehold A / Toehold B : 16nt / 9nt / 13nt <br>
 +
Temperature : 300K <br>
 +
Time-step / simulation length : 0.01ps / 100ns<br>
 +
Ion concentration : 50mM Na+<br>
-
Add all potential by negative charge DNA which compose gate have.<br>
+
<html><div style="clear:both;"></div></html>
-
(used C language to output the numbers)<br><br>
+
-
<img src="http://openwetware.org/images/4/45/Helix.gif" width="500px" height="290px"><br>
+
 +
<h3>Results</h3>
 +
<table border="0" align="left" vertical-align: middle;>
 +
<tr>
 +
<td>
 +
<youtube width="450" align="left">FGpemnMLUMk&border=1&color1=0x6699&color2=0x54abd6</youtube>
 +
<html><div style="clear:both;"></div></html>
 +
  </td>
 +
</tr>
 +
<tr>
 +
<td align="center">
 +
Movie1 : MD simulation of toehold strand displacement
 +
</td>
 +
</tr>
 +
</table>
 +
<html><div style="clear:both;"></div></html>
-
Condition<br>
 
-
Temperature 298[K]<br>
 
-
Na+ 50mM<br>
 
-
<img src="http://openwetware.org/images/f/fc/Add2.png" width="500px" height="180px"><br>
+
Movie 1 shows the trajectory of each strand from the simulation. The target strand moves from Toehold A strand to Toehold B strand which are fixed on the field. This result agrees with the energy gradient. The result supported our coarse-grained DNA model is feasible.
 +
<h2> Comparison of capture ability</h2>
 +
One of constructional features of our structure ”CELL-GATE”is the use of a novel strand displacement method, called porter system.
-
</p>
+
By comparing our porter strand to toehold strand, the most popular method for strand displacement, we looked at the effectiveness of our structure in terms of it's strand capturing ability.
-
<a name="Results"></a><h2>Results</h2>
+
<h3> Model and Method</h3>
-
<p>
+
According to the design of [http://openwetware.org/wiki/Biomod/2012/TeamSendai/Experiment#Comparison_of_Porter_and_Toehold the experiment section], we modeled the porter strand and the toehold strand as shown below. <br><br>
-
<br>
+
[[Image: Format_selector.jpg |left|680px]]
-
Electric potential changing z-axis at x-axis and y-axis is 0.<br>
+
{{-}}
-
<img src="http://openwetware.org/images/0/09/1014x0y0potential.png" width="620px" height="450px"><br>
+
<html><div style="clear:both;"></div></html>
-
the length of the gate is 88bp, 30nm.
+
-
Target base pair 25 を点電荷と仮定する
+
Hex-cylinder is represented as the assembly of fixed electrically-charged mass points.<br><br>
-
</p>
+
[[Image: Format_Hexagon.jpg |left|750px]]
 +
<html><div style="clear:both;"></div></html>
-
<a name="MD Simulation"></a><h2>MD Simulation</h2>
+
Simulation was carried out at the following condition:
-
<p>
+
-
We carried out molecular dynamics simulation to examine the the mechanism and
+
-
the effectiveness of our structure “Cell Gate”.
+
-
</p>
+
Temperature : 300K<br>
 +
Ion concentration : Na+ 50mM<br>
 +
Box size : 20nm×20nm×20nm (periodic boundary condition)<br>
 +
Time-step / simulation length : 0.01ps / 10ns<br>
 +
{{-}}
-
<a name="DNA Model"></a><h2>DNA Model</h2>
 
-
<p>
 
-
For simplicity, course-grained DNA model is used in our simulation. <br>
 
-
One DNA nucleotide is represented by one bead in the model and each bead can be<br>
 
-
hybridized with complementary bead.<br>
 
-
  <<モデル載せる>><br><br>
 
-
The potential energy of the system includes 5 distinct contributions.<br>
 
-
  <<ポテンシャル載せる>><br><br><br>
 
-
The first three terms are intramolecular interactions , bonds , bond angles, and<br>
 
-
dihedral angles. In order to express “tether like structure”, only bond interactions<br>
 
-
are active in our DNA model.<br>
 
-
And the latter two are non-bonded interactions. Coulomb interactions are taken into<br>
 
-
account using the Debye-Huckel approximation which enables to internalize<br>
 
-
counterions contribution.<br>
 
-
Constants of these potentials are achieved from references.<br>
 
-
The force on bead i is given by a Langevin equation<br><br><br>
 
-
Langevin equation<br><br>
+
<h3>Results and Conclusion</h3>
-
+
-
<img src="http://openwetware.org/images/1/11/Langevin.png" width="220px" height="80px"><br><br>
+
-
<img src="http://openwetware.org/images/2/23/F%3D.png" width="150px" height="80px"><br>
+
-
The first term donates a conservative force derived from the potential U and the<br>
+
<table border="0" align="left" vertical-align: middle;>
-
second is a viscosity dependent friction.<br>
+
<tr>
-
The third term is a white Gaussian noise and effects of solvent molecules are<br>
+
<td>
-
internalized in this term.<br>
+
<youtube width="450" align="left">NXo4cYKkrF0&border=1&color1=0x6699&color2=0x54abd6</youtube>
-
Langevin equation is integrated using a Velocity-Verlet method.<br><br><br>
+
<html><div style="clear:both;"></div></html>
-
Toehold displacement of dsDNA<br>
+
</td>
-
In order to test predictive capability of the model, here we carried out a simulation<br>
+
<td width="20">
-
of Toehold displacement between two strands.<br>
+
</td>
-
Length of strands and simulation situation was as follows.<br><br>
+
<td>
-
Target strand/Toehold A/Toehold B : 25nt / 9nt (+10nt spacer) / 13nt (+10nt
+
<youtube width="450" align="left">qidZS1pI0lc&border=1&color1=0x6699&color2=0x54abd6</youtube>
-
spacer)<br>
+
<html><div style="clear:both;"></div></html>
-
Temperature : 300K<br>
+
</td>
-
Time-step size / simulation length : 0.01ps / 100ns<br>
+
</tr>
-
Ion concentration : 50mM Na+<br><br>
+
<tr>
-
results<br>
+
<td align="center">
-
<<後ほど>>
+
Movie2 : MD simulation of porter strand
 +
</td>
 +
<td width="20">
 +
</td>
 +
<td align="center">
 +
Movie3 : MD simulation of toehold strand
 +
</td>
 +
</tr>
 +
</table>
 +
<html><div style="clear:both;"></div></html>
-
</p>
 
-
<a name="Comparison of capture ability"></a><h2>Comparison of capture ability</h2>
+
Movie2 and 3 shows the result of each simulation, porter-target and toehold-target.
-
<p>
+
-
One of constructional features of our structure ”Cell-Gate” is the use of new strand
+
-
displacement method.<br>
+
-
By comparing our selector strand and a toehold strand, the most popular method for<br>
+
-
strand displacement, we show the effectiveness our structure in terms of capture
+
-
ability.<br><br><br>
+
-
Model and Method<br>
+
-
According to the design of experiment section, we designed models as below of the<br>
+
-
selector strand and the toehold strand.<br>
+
-
<<モデル載せる>><br><br><br>
+
-
Hex-cylinder is represented as the assembly of electrically-charged mass points<br>
+
-
fixed on the field.<br>
+
-
<<モデル載せる>><br><br><br>
+
-
Simulation was carried out at the following condition.<br>
+
-
Temperature : 300K<br>
+
-
Ion concentration : Na+ 50mM<br>
+
-
Box size : 20nm×20nm×20nm (periodic boundary condition)<br>
+
-
Time-step size / simulation length : 0.01ps / 10ns<br>
+
-
Results<br>
+
-
<<後ほど>><br>
+
-
</p>
+
We note that this simulation was carried out under a periodic boundary condition, where the size of the simulation box is 20nm×20nm×20nm. Then, the distance between the target strand and the Hex-cylinder is maintained virtually constant.
-
<a name="Reference"></a><h2>Reference</h2>
+
One of the advantages of the porter strand is shrinking ability. The porter strand hybridizes to the target making a loop which makes it possible to extend the strand length without changing the final structure's length.
-
<p>
+
-
1. Thomas A. Knotts et al. A coarse grain model of DNA , J.Chem.Phys
+
-
126,084901(2007)<br>
+
-
2. Carsten Svaneborg et al. DNA Self-Assembly and Computation Studied with a
+
-
Coarse-Grained Dynamic Bonded Model, DNA 18,LNCS 7433, pp.123-134,
+
-
2012<br>
+
-
3. Xhuysn Guo & D.Thirumalai, Kinetics of Protein Folding: Nucleation
+
-
Mechanism, Time Scales, and Pathways, Biopolymars, Vol.36, 83-102 (1995)<br>
+
-
4. GROMACS manual ()<br>
+
-
5. Cafemol manual ( http://www.cafemol.org/ )<br>
+
-
</p>
 
-
</div>
+
Results obtained from this simulation show that the porter strand can catch the target strand exists outside of the Hex-cylinder and hybridize completely while the toehold strand never hybridize to the target strand in simulation time.
-
</body>
+
We run additional 5 simulations for each capturing mechanism under the same conditions (NOTE: random seeds were different) and very similar results were obtained. Thus, the simulation results are very feasible under our model.
-
</html>
+
 
 +
 
 +
By considering results of electrostatic potential calculation around the hex-cylinder and MD simulation, it is clear that the electrostatic field prevents the entrance of DNA strands into the Hex-cylinder and the porter strand helps it to get into the cylinder.
 +
 
 +
Therefore, we concluded that '''our novel porter strand provides a high capture ability to our system “CELL-GATE”.'''
 +
 
 +
 
 +
 
 +
<h2>References</h2>
 +
1. Thomas A. Knotts ''et al''.  A coarse grain model of DNA , J.Chem.Phys 126,084901(2007)
 +
 
 +
2. Carsten Svaneborg ''et al''. DNA Self-Assembly and Computation Studied with a Coarse-Grained Dynamic Bonded Model, DNA 18,LNCS 7433, pp.123-134, (2012)
 +
 
 +
3. Xhuysn Guo & D.Thirumalai, Kinetics of Protein Folding: Nucleation Mechanism, Time Scales, and Pathways, Biopolymars, Vol.36, 83-102 (1995)
 +
 
 +
4. GROMACS manual ( http://www.gromacs.org/ )
 +
 
 +
5. Cafemol manual  ( http://www.cafemol.org/ )
 +
 
 +
 
 +
 
 +
 
 +
</div>

Current revision



Team Sendai Top


Contents

Aim of Simulation

Molecular robots works in meso-scale (from nm to sub-micro meter). In the meso-scale, characters of each molecule and cooperative movement of molecules could not be ignored. Without calculation of potentials or molecular dynamics simulations, it is difficult to predict the behaviors of the molecules in this scale. Thus, we calculated potentials to confirm the GATE effect, and simulated molecular dynamics of DNAs to strength reliability of our design.

Numerical Analysis of Electrostatic Potential

Electric potential of the GATE was calculated inside and outside the GATE.
The phosphate groups in the backbone of DNA are negatively-charged. Because the GATE is made by DNA origami, the GATE has the Coulomb force.

Model of DNA

A point-charge model is applied to DNA model. Each phosphate of a nucleotide was converted into one electric beads, which has elementary electoric charge (1.602176565 x 10^−19 [C]). In DNA helix, negative charges appear along the axis of the double helix every 10.4 base pairs.
Debye–Hückel equation was used to calculate electric potential. Electric potential of each point around the GATE was deduced from summation of all the Electric Potential for every phosphate presented on DNA.


Debye–Hückel equation:


reference: Debye length



Following GIF animation shows the scheme of a point-charge model and potential calculations.
(It takes about 7 seconds to watch one cycle of the GIF animation)




Results

Electric potential was calculated by using C language.
Under our conditions, the target DNA has 25 electric beads.
All calculation were performed under these condition: Temperature 298[K], Na+ 50mM.

Electric potential from the top of view of the GATE

The following figure shows electric potentials between the target DNA and the GATE from the top of view of the GATE.
The electric potentials were standardized by thermal energy, KBT (KB; Boltzmann constant).
Potential energies were shown by a heat map.

Electric potential from the top of view of the GATE



Potential energy along axis indicated by the arrows are shown as the following figure.
Together, electric potential is too high for the target DNA to enter inside the GATE.


Electric potential from the top of view of the GATE along an axis



Electric potential at the center of hole in the GATE

The difficulty to go through inside hole of the GATE was also indicated by following figures.
The next figures shows electric potential at the center of hole in the GATE.
Potentials were calculated along the axis indicated by the red arrow.
Blue lines shows electric potential in the case of 1.5 fold radius of the GATE.


Electric potential at the center of hole in the GATE



By the figure, we concluded
i ) Our designed size is suitable size for the GATE function
ii) Enlarging radius of the GATE decreases electric repulsion effect


Hybridization energy of Porters overcomes electric potential of the GATE

We expected the hybridization energy of the Porters overcomes the inhibition of entrance of the GATE by electric repulsion.
The Gibbs energies of hybridization was calculated by neighbor-joining methods.

The result of Gibbs energy of the fixed Porter were defined as the maximum length of each form of Porter. For examples, Porter 1 has three forms: Stretched form, 1 loop form, and 2 loop form. In this Porter 1 case, we defined Porter 1 has three Gibbs energy according to the hybridization state. Since the maximum lengths of each form were different, the minimum energy of possible length was selected.
Gibbs energies were calculated at the center of hole in the GATE.
The following figure shows the calculated energy of Porters and Toehold DNA used as controls in our project.


Three different form



Gibbs energies by the Porters in the GATE



The Gibbs energies by the Porters shown above were summed with electric potential of the GATE.
The summation of energies were calculated at the center of hole in the GATE.
The result was shown in the next figure.

Electric potential from the top of view of the GATE


The results of calculation indicate that the outside Porter can catch the target DNA, and inner Porter can pull in the GATE step by step. On the other hand, short DNA like the toehold DNA we called n this study cannot catch the target DNA.


MD Simulation

To further strengthen feasibility of our design, we examined the capturing mechanism and the effectiveness of our structure “CELL GATE" by molecular dynamics (MD) simulation.


DNA Model

For simplicity, a course-grained DNA model was used in our simulation. One DNA nucleotide was represented by one bead in the model and each bead can be hybridized with a complementary bead.

The potential energy of the system includes 5 distinct contributions.


The first three terms are intramolecular interactions, bonds, bond angles, and dihedral angles. These terms are important for maintaining structure. In order to express the “tether like structure”, only bond interactions are considered in our DNA model.

And the latter two terms are non-bonded interactions. Coulomb interactions are taken into account using the Debye-Huckel approximation which enables to internalize counter-ions contribution.

Parameters of these potentials were fit to the reference literature ; Thomas A. Knotts et al. A coarse grain model of DNA .

The force on bead i is given by a Langevin equation


The first term donates a conservative force derived from the potential Vtot and the second is a viscosity dependent friction.

The third term is a white Gaussian noise and effects of collision with solvent molecules which causes brownian motion are internalized in this term.

Langevin equation is integrated using a Velocity-Verlet method.


Verification of the Coarse-grained model: Toehold displacement of dsDNA

In order to test the model, here we carried out a simulation of Toehold displacement between two strands.

Toehold strands are modeled according to designs of the experiment section.


Length of strands and simulation condition were as follows:

Target strand / Toehold A / Toehold B : 16nt / 9nt / 13nt
Temperature : 300K
Time-step / simulation length : 0.01ps / 100ns
Ion concentration : 50mM Na+

Results

Movie1 : MD simulation of toehold strand displacement


Movie 1 shows the trajectory of each strand from the simulation. The target strand moves from Toehold A strand to Toehold B strand which are fixed on the field. This result agrees with the energy gradient. The result supported our coarse-grained DNA model is feasible.


Comparison of capture ability

One of constructional features of our structure ”CELL-GATE”is the use of a novel strand displacement method, called porter system.

By comparing our porter strand to toehold strand, the most popular method for strand displacement, we looked at the effectiveness of our structure in terms of it's strand capturing ability.

Model and Method

According to the design of the experiment section, we modeled the porter strand and the toehold strand as shown below.


Hex-cylinder is represented as the assembly of fixed electrically-charged mass points.

Simulation was carried out at the following condition:

Temperature : 300K
Ion concentration : Na+ 50mM
Box size : 20nm×20nm×20nm (periodic boundary condition)
Time-step / simulation length : 0.01ps / 10ns


Results and Conclusion

Movie2 : MD simulation of porter strand

Movie3 : MD simulation of toehold strand


Movie2 and 3 shows the result of each simulation, porter-target and toehold-target.

We note that this simulation was carried out under a periodic boundary condition, where the size of the simulation box is 20nm×20nm×20nm. Then, the distance between the target strand and the Hex-cylinder is maintained virtually constant.

One of the advantages of the porter strand is shrinking ability. The porter strand hybridizes to the target making a loop which makes it possible to extend the strand length without changing the final structure's length.


Results obtained from this simulation show that the porter strand can catch the target strand exists outside of the Hex-cylinder and hybridize completely while the toehold strand never hybridize to the target strand in simulation time.

We run additional 5 simulations for each capturing mechanism under the same conditions (NOTE: random seeds were different) and very similar results were obtained. Thus, the simulation results are very feasible under our model.


By considering results of electrostatic potential calculation around the hex-cylinder and MD simulation, it is clear that the electrostatic field prevents the entrance of DNA strands into the Hex-cylinder and the porter strand helps it to get into the cylinder.

Therefore, we concluded that our novel porter strand provides a high capture ability to our system “CELL-GATE”.


References

1. Thomas A. Knotts et al. A coarse grain model of DNA , J.Chem.Phys 126,084901(2007)

2. Carsten Svaneborg et al. DNA Self-Assembly and Computation Studied with a Coarse-Grained Dynamic Bonded Model, DNA 18,LNCS 7433, pp.123-134, (2012)

3. Xhuysn Guo & D.Thirumalai, Kinetics of Protein Folding: Nucleation Mechanism, Time Scales, and Pathways, Biopolymars, Vol.36, 83-102 (1995)

4. GROMACS manual ( http://www.gromacs.org/ )

5. Cafemol manual ( http://www.cafemol.org/ )



Personal tools