Team:Tokyo Tech/Modeling/RPS-game/RPS-game

From 2011.igem.org

(Difference between revisions)
 
(32 intermediate revisions not shown)
Line 35: Line 35:
list-style: none;
list-style: none;
float: left;
float: left;
-
height: 20px;
+
height: 30px;
width: 945px
width: 945px
}
}
Line 53: Line 53:
display: none;
display: none;
position: absolute;
position: absolute;
-
font-size: 14px;
+
font-size: 16px;
opacity: 0.8;
opacity: 0.8;
list-style: none;
list-style: none;
Line 77: Line 77:
background-color: #EFEBEC;
background-color: #EFEBEC;
color: #000000;
color: #000000;
 +
font-size:13px;
}
}
Line 84: Line 85:
left: 0;
left: 0;
clear: both;
clear: both;
-
height: 180px;
+
height: 190px;
-
width: 752px;
+
}
}
.main
.main
Line 107: Line 107:
background-color: #00508D;
background-color: #00508D;
color:#FFFFFF;
color:#FFFFFF;
-
text-align: center;
+
text-align: left;
padding-left: 5px;
padding-left: 5px;
clear: both;
clear: both;
 +
}
 +
 +
h2
 +
{
 +
text-indent: 1em;
}
}
 +
h3
 +
{
 +
text-indent: 2em;
 +
}
 +
 +
h4
 +
{
 +
text-indent: 3em;
 +
}
 +
 +
h5
 +
{
 +
text-indent: 4em;
 +
}
 +
 +
h6
 +
{
 +
text-indent: 5em;
 +
}
 +
.footer
.footer
{
{
Line 145: Line 170:
}
}
-
        img.center
+
img.center
{
{
float: center;
float: center;
Line 190: Line 215:
}
}
-
window.onload = function() {
+
function getElementsByClass() {
 +
    var classElements = new Array();
 +
    var allElements = document.getElementsByTagName("*");
 +
    for (i = 0 ; i < allElements.length; i++) {
 +
if (allElements[i].className == 'firstHeading') {
 +
    allElements[i].style.display = 'none';
 +
  // window.alert('発見。');
 +
}
 +
    }
 +
}
 +
 +
window.onload = function() {
 +
getElementsByClass();
var menu = window.document.getElementById('TopMenu');
var menu = window.document.getElementById('TopMenu');
if(menu==null)
if(menu==null)
Line 223: Line 260:
<!-- top box -->
<!-- top box -->
<div class="top">
<div class="top">
-
<OBJECT classid="clsid:D27CDB6E-AE6D-11cf-96B8-444553540000" codebase="http://download.macromedia.com/pub/shockwave/cabs/flash/swflash.cab#version=6,0,0,0" WIDTH="965" HEIGHT="150" id="Yourfilename" ALIGN="">
+
<object classid="clsid:D27CDB6E-AE6D-11cf-96B8-444553540000" codebase="http://download.macromedia.com/pub/shockwave/cabs/flash/swflash.cab#version=6,0,0,0" width="965" height="150" id="Yourfilename" align="">
-
<PARAM NAME=movie VALUE="https://static.igem.org/mediawiki/2011/4/4d/Header.swf"> <PARAM NAME=quality VALUE=high> <PARAM NAME=bgcolor VALUE=#FFFFFF> <EMBED src="https://static.igem.org/mediawiki/2011/4/4d/Header.swf" quality=high bgcolor=#FFFFFF WIDTH="965" HEIGHT="150" NAME="Yourfilename" ALIGN="" TYPE="application/x-shockwave-flash" PLUGINSPAGE="http://www.macromedia.com/go/getflashplayer"></EMBED> </OBJECT>  
+
<param name="movie" value="https://static.igem.org/mediawiki/2011/4/4d/Header.swf">
 +
<param name="quality" value="high">
 +
<param name="bgcolor" value="#FFFFFF">
 +
<embed src="https://static.igem.org/mediawiki/2011/4/4d/Header.swf" quality="high" bgcolor="#FFFFFF" width="965" height="150" name="Yourfilename" align="" type="application/x-shockwave-flash" pluginspage="http://www.macromedia.com/go/getflashplayer">
 +
</embed>
 +
</object>  
<!-- list of top menu -->
<!-- list of top menu -->
<div id="navigation">
<div id="navigation">
Line 255: Line 297:
<li id="menu_Extra">
<li id="menu_Extra">
More
More
-
<ul>
+
<ul style="width:210px;">
<li><a href="https://2011.igem.org/Team:Tokyo_Tech/Safety">Safety</a></li>
<li><a href="https://2011.igem.org/Team:Tokyo_Tech/Safety">Safety</a></li>
<li><a href="https://2011.igem.org/Team:Tokyo_Tech/Attribution_and_Contributions.htm">Attribution and Contributions</a></li>
<li><a href="https://2011.igem.org/Team:Tokyo_Tech/Attribution_and_Contributions.htm">Attribution and Contributions</a></li>
Line 274: Line 316:
<!-- left menu list -->
<!-- left menu list -->
-
<div style="min-height: 6500px; float: left;">
+
<div style="min-height: 6100px; float: left;">
<div id="LeftMenu">
<div id="LeftMenu">
<!--list of page menu: DO NOT WRITE LINKS NOT WRITTEN IN THIS PAGE -->
<!--list of page menu: DO NOT WRITE LINKS NOT WRITTEN IN THIS PAGE -->
<ul>
<ul>
-
<li>
+
<li><a href="#1.1">1. Introduction</a></li>
-
<a href="#1.abst">1. Abstract</a>
+
<li><a href="#1.2">2. Adjusting the Model to create a True Randomizer</a></li>
-
</li>
+
<li><a href="#1.3">3. How the Three Types of Bacteria Compete for Survival</a></li>
 +
<li><a href="#1.4">4. The Old Model</a></li>
 +
<li><a href="#1.5">5. Our New Model</a></li>
 +
<li><a href="#1.6">6. The Biological Meaning of our Model</a></li>
 +
<li><a href="#1.7">7. Making it Obvious</a></li>
</ul>
</ul>
</div>
</div>
Line 294: Line 340:
<p>
<p>
-
                <img src="https://static.igem.org/mediawiki/2011/2/28/Image042.png" />
+
<div align="center">
-
             
+
<img src="https://static.igem.org/mediawiki/2011/2/28/Image042.png" width="500px" />
-
        <h2 id="1.1">Introduction: Minimal differences determine who will survive</h2>
+
</div>
-
                <p>
+
             
-
In this section we will show a shocking scenario of evolution: the future of each of three different rival strains (whether the strain will die or survive) is marked by minimal differences between the initial population densities of the strains. Furthermore, we will also show that we can apply this very interesting result to create a randomizer that can be used in our Rock-Paper-Scissors game, due to the fact that only one of the rival strains will survive. More specifically, we assign to each of the three rival strains either of Rock, Paper or Scissors, make them compete for survival and take the surviving strain to represent the bacteria’s choice for the RPS game.
+
<h2 id="1.1">1. Introduction: Very small differences determine who will survive</h2>
-
                </p>
+
-
             
+
<p>
-
                <h2 id="12"> Adjusting the Model to create a True Randomizer</h2>
+
In this section we will show a shocking scenario of evolution:  
-
                <p>
+
the future of each of three different rival strains (whether the strain will die or survive)  
-
The idea for creating this randomizer was born from a paper written in 1996 by Durret and Levin. In it, the authors described a system of three types of bacteria that competed for survival in dynamic that resembled a Rock-Paper-Scissors (RPS) game. However, the model proposed in this paper is not fully appropriate for our RPS randomizer, since one of the three types of bacteria cannot ultimately survive (although it can dominate the system, i.e. have the highest population density, for definite periods of time). We will discuss more on the limitations we found in this model to be adopted as a randomizer and the modifications we made to create a true randomizer.
+
is marked by very small differences between the initial population densities of the strains(a phenomenon also known as the "butterfly effect").  
-
                </p>
+
Furthermore, we will also show that we can apply this very interesting result to create a  
-
               
+
randomizer that can be used in our Rock-Paper-Scissors game, due to the fact that only one  
-
                <h2 id="1.3"> How the Three Types of Bacteria Compete for Survival</h2>
+
of the rival strains will survive. More specifically, we assign to each of the three  
 +
rival strains either of Rock, Paper or Scissors, make them compete for survival and take  
 +
the surviving strain to represent the bacteria’s choice for the RPS game.
 +
</p>
 +
             
 +
<h2 id="1.2">2. Adjusting the Model to create a True Randomizer</h2>
 +
 +
<p>
 +
The idea for creating this randomizer was born from a paper written in 1996 by Durrett and Levin.  
 +
In it, the authors described a system of three types of bacteria that competed for  
 +
survival in dynamic that resembled a Rock-Paper-Scissors (RPS) game. However,  
 +
the model proposed in this paper is not fully appropriate for our RPS randomizer,  
 +
since one of the three types of bacteria cannot ultimately survive (although it  
 +
can dominate the system, i.e. have the highest population density, for definite periods  
 +
of time). We will discuss more on the limitations we found in this model to be adopted as  
 +
a randomizer and the modifications we made to create a true randomizer.
 +
</p>
 +
               
 +
<h2 id="1.3">3. How the Three Types of Bacteria Compete for Survival</h2>
 +
 +
<p>
 +
The three types of bacteria that compete for survival use three tactics to
 +
outcompete their rivals: the production of a toxin (a bacteriocin called colicin)
 +
that is toxic to other strains, resistance to the toxin produced by other strains,
 +
and a higher birth rate than their rival strains. Namely, the three types
 +
of bacteria are: colicin-producing <span class="name">E. coli</span> (R),
 +
colicin-resistant <span class="name">E. coli</span> (P) and colicin-sensitive
 +
<span class="name">E. coli</span> (S). The colicin-producer outcompetes the
 +
colicin-sensitive by producing the colicin. The colicin-sensitive bacteria
 +
outcompetes the colicin-resistant because its birth rate is higher
 +
than that of the colicin-resistant. The colicin-resistant outcompetes the
 +
colicin producer because its birth rate is higher than that of the colicin producer.
 +
The colicin resistant bacteria are also able to produce colicin,
 +
but at a lower energetic cost, which allows them to have a higher birth rate.
 +
</p>
 +
 +
<div align="center">
 +
<img src="https://static.igem.org/mediawiki/2011/b/bf/Image044.png" width="500px" />
 +
</div>
 +
 +
<p>
 +
The system was described by the following general differential equations
 +
</p>
 +
 +
<div align="center">
 +
<img src="https://static.igem.org/mediawiki/2011/2/28/Colisin1.png" />
 +
</div>
-
                <p>
+
<p>
-
The three types of bacteria that compete for survival use three tactics to outcompete their rivals: the production of a toxin (a bacteriocin called colicin) that is toxic to other strains, resistance to the toxin produced by other strains, and a higher birth rate than their rival strains. Namely, the three types of bacteria are: colicin-producing E. coli (R), colicin-resistant E. coli (P) and colicin-sensitive E. coli (S). The colicin-producer outcompetes the colicin-sensitive by producing the colicin. The colicin-sensitive bacteria outcompetes the colicin-resistant because its birth rate is higher than that of the colicin-resistant. The colicin-resistant outcompetes the colicin producer because it’ birth rate is higher than that of the colicin producer. The colicin resistant bacteria are also able to produce colicin, but at a lower energetic cost, which allows them to have a higher birth rate.
+
Where
-
                </p>
+
</p>
-
<div align="center">
+
<div align="center">
 +
<img src="https://static.igem.org/mediawiki/2011/c/c0/Colisin2.png" />
 +
</div>
-
                <img src="https://static.igem.org/mediawiki/2011/b/bf/Image044.png" />
 
-
</div>
 
-
                <p>
 
-
The system was described by the following general differential equations
 
-
                </p>
 
-
 
-
<div align="center">
 
-
                <img src="https://static.igem.org/mediawiki/2011/2/28/Colisin1.png" />
 
-
</div>
 
-
                <p>
 
-
Where
 
-
                </p>
 
-
<div align="center">
+
<h2 id="1.4">4. The Old Model</h2>
-
                <img src="https://static.igem.org/mediawiki/2011/c/c0/Colisin2.png" />
+
-
</div>
+
<p>
-
                <h2 id="1.4">The Old Model</h2>
+
In the model described by Durrett and Levin’s paper the equations were as follows:
-
                <p>
+
</p>
-
In the model described by Durret and Levin’s paper the equations were as follows:
+
-
                </p>
+
-
               
+
-
<div align="center">
+
-
                <p>Producer</p>
+
-
                <img src="https://static.igem.org/mediawiki/2011/1/1d/Colisin3.png" /><br />
+
-
                <p>Resistant</p>
+
-
                <img src="https://static.igem.org/mediawiki/2011/0/04/Colisin4.png" /><br />
+
-
                <p>Sensitive</p>
+
-
                <img src="https://static.igem.org/mediawiki/2011/3/36/Colisin5.png" /><br />
+
-
</div>
+
-
                <p>
+
-
These equations show that the colicin-resistant bacteria are completely immune to colicin (there is not death factor associated to colicin in the equation for du2/dt). However, as will be explained afterwards, this results in a loss of balance that does not allow building a true randomizing system.
+
-
                </p>
+
-
                <p>
+
-
Now, setting the parameters as follows, the graph below was created by Durret and Levin.
+
-
                </p>
+
-
<div align="center">
+
<div align="center">
 +
Producer<br />
 +
<img src="https://static.igem.org/mediawiki/2011/1/1d/Colisin3.png" /><br />
 +
Resistant<br />
 +
<img src="https://static.igem.org/mediawiki/2011/0/04/Colisin4.png" /><br />
 +
Sensitive<br />
 +
<img src="https://static.igem.org/mediawiki/2011/3/36/Colisin5.png" /><br />
 +
</div>
-
                <img src="https://static.igem.org/mediawiki/2011/4/41/Colisin6.png" /><br />
+
<p>
-
                <img src="https://static.igem.org/mediawiki/2011/e/e3/Colisin7.png" />
+
These equations show that the colicin-resistant bacteria are completely immune to colicin
-
                <img src="https://static.igem.org/mediawiki/2011/thumb/e/e5/Image086.png/669px-Image086.png" width="354px"/><br />
+
(there is not death factor associated to colicin in the equation for du<sub>2</sub>/dt). However, as will be explained afterwards, this results in a loss of balance that does not allow building a true randomizing system.
-
                <img src="https://static.igem.org/mediawiki/2011/thumb/b/b1/Colisin8.png/800px-Colisin8.png" width="800px"/><br />
+
</p>
-
                <img src="https://static.igem.org/mediawiki/2011/b/b8/Image094.png" width="354px"/><br />
+
-
</div>
+
-
                <h2 id="1.5"> Our New Model</h2>
+
-
                <p>
+
-
As mentioned before, the model proposed by Durret and Levin has critical limitations as a randomizer for the RPS game. To be able to create a true randomizer, we modified the differential equations of the model taking care to give it a biological meaning. With our new differential equations, any of the three types of bacteria can ultimately survive by outcompeting the other two strains, which will die. More specifically, we limited the resistance of the colicin-resistant bacteria in the sense that it would produce a type of bacteriocin that is only toxic to itself and to the sensitive strain, and additionally the resistant strain would also be vulnerable to the colicin produced by the colicin-producer. Since which strain will be the one that survives is determined by very small differences in the initial concentrations of the three different populations of bacteria, in practice this systems becomes a randomizer because of the imprecisions in the measurements that result, for example, when using micropipettes. This randomizer describes a new competition dynamic that could not be reproduced in the previous model proposed by Durret and Levin due to the instability along the
+
-
<img src="https://static.igem.org/mediawiki/2011/e/e0/Image088.png" alt="u1">
+
<p>
 +
Now, setting the parameters as follows, the graph below was created by Durrett and Levin.
 +
</p>
-
axis.
+
<div align="center">
-
                </p>
+
<img src="https://static.igem.org/mediawiki/2011/4/41/Colisin6.png" /><br />
-
<div align="center">
+
<img src="https://static.igem.org/mediawiki/2011/e/e3/Colisin7.png" />
-
               
+
<img src="http://partsregistry.org/wiki/images/b/bd/Modeling1.png" width="354px"/><br />
-
                <img src="https://static.igem.org/mediawiki/2011/c/c1/Colimodel1.png" alt="our new model" />
+
<img src="https://static.igem.org/mediawiki/2011/thumb/b/b1/Colisin8.png/800px-Colisin8.png" width="800px"/><br />
-
</div>
+
<img src="https://static.igem.org/mediawiki/2011/5/55/Modeling2.png" width="354px"/><br />
-
                <p>If we set the parameters as follows</p>
+
</div>
-
<div align="center">
+
-
                <img src="https://static.igem.org/mediawiki/2011/0/00/Colimodel2.png" alt="new model's coefficient" />
+
-
</div>
+
-
                <p>
+
<h2 id="1.5">5. Our New Model</h2>
-
and we graph this equations using a Matlab program, we get a graph which clearly shows there are stable points on each of the three axes (Figure 1, Up).
+
<p>
-
                </p>
+
As mentioned before, the model proposed by Durrett and Levin has critical limitations as a randomizer for the RPS game.
 +
To be able to create a true randomizer, we modified the differential equations of the model taking care to give it a  
 +
biological meaning. With our new differential equations, any of the three types of bacteria can ultimately survive by
 +
outcompeting the other two strains, which will die. More specifically, we limited the resistance of the colicin-resistant
 +
bacteria in the sense that it would produce a type of bacteriocin that is only toxic to itself and to the sensitive strain,
 +
and additionally the resistant strain would also be vulnerable to the colicin produced by the colicin-producer. 
 +
Since which strain will be the one that survives is determined by very small differences in the initial concentrations of  
 +
the three different populations of bacteria, in practice this systems becomes a randomizer because of the imprecisions in
 +
the measurements that result, for example, when using micropipettes. This randomizer describes a new competition dynamic
 +
that could not be reproduced in the previous model proposed by Durrett and Levin due to the instability along the
 +
<img src="https://static.igem.org/mediawiki/2011/e/e0/Image088.png" alt="u1">axis.
 +
</p>
 +
 +
<div align="center">
 +
<img src="https://static.igem.org/mediawiki/2011/c/c1/Colimodel1.png" alt="our new model" />
 +
</div>
 +
<p>If we set the parameters as follows</p>
 +
<div align="center">
 +
<img src="https://static.igem.org/mediawiki/2011/0/00/Colimodel2.png" alt="new model's coefficient" />
 +
</div>
 +
<p>
 +
and we graph this equations using a Matlab program, we get a graph which clearly shows there are stable points on each of the three axes (Figure 1, Up).
 +
</p>
 +
 +
<img src="https://static.igem.org/mediawiki/2011/b/b8/Model2.png" width="815px" style="float:center;"/><br />
 +
<img src="https://static.igem.org/mediawiki/2011/0/05/Image106.png" width="815px" style="float:center;"/>
 +
 +
<div align="center">
 +
Figure 1. Up: Our New model. Down: The Old Model
 +
</div>
-
<img src="https://static.igem.org/mediawiki/2011/b/b8/Model2.png" width="100%" style="float:center;"/>
+
<p>
-
<img src="https://static.igem.org/mediawiki/2011/0/05/Image106.png" width="100%" style="float:center;"/>
+
These stable points (u<sub>1</sub>,0,0), (0,u<sub>2</sub>,0) and (0,0,u<sub>3</sub>) indicate that for the equations we have set
 +
all of the three strains may ultimately survive for infinite peiriods of time. The differences between
 +
our model and the model of Durrett and Levin can be seen graphically in Figure 1. These graphs were plotted using Matlab.
 +
</p>
-
                <p>
+
<p>
-
                <center>Figure 1. Up: Our New model. Down: The Old Model</center>
+
Note that the parameters we have set for our equations satisfy the initail
-
                </p>
+
conditions of the model proposed by Durrett and Levin (indicated in black font)<br />
-
                <p>
+
<img src="https://static.igem.org/mediawiki/2011/a/a4/Colimodel4.png" alt="new terms" width="800px"/>
 +
</p>
-
These stable points (u1,0,0), (0,u2,0) and (0,0,u3) indicate that for the equations we have set all of the three strains may ultimately survive for infinite peiriods of time. The differences between our model and the model of Durret and Levin can be seen graphically in Figure 1. These graphs were plotted using Matlab.
+
<h2 id="1.6">6. The Biological Meaning of our Model</h2>
-
</p>
+
-
                <p>
+
<p>
-
Note that the parameters we have set for our equations satisfy the initail conditions of the model proposed by Durret and Levin (indicated in black font)
+
From a biological perspective, our model describes the existence of two strains of
-
                </p>
+
bacteria that produce two different types of bacteriocins. One of these strains is not
 +
completely resistant to its own bacteriocin nor to the bacteriocin produced by its
 +
rival strain. This can be justified as the consequence of insufficient/ineffective
 +
resistance protein production by the &ldquo;resistant&rdquo; strain. This limitation in  
 +
the production of resistance protein could be thought of as a consequence of the
 +
&ldquo;resistant&rdquo; strain being a mutant of a colicin-sensitive strain.
 +
</p>
-
                <img src="https://static.igem.org/mediawiki/2011/a/a4/Colimodel4.png" alt="new terms" width="800px"/>
+
<h2 id="1.7">7. Making it Obvious</h2>
-
 
+
-
                <h4 id="3.3.6">3.3.6 The Biological Meaning of our Model</h4>
+
-
 
+
-
                <p>
+
-
From a biological perspective, our model describes the existence of two strains of bacteria that produce two different types of bacteriocins. One of these strains is not completely resistant to its own bacteriocin nor to the bacteriocin produced by its rival strain. This can be justified as the consequence of insufficient/ineffective resistance protein production by the “resistant” strain. This limitation in the production of resistance protein could be thought of as a consequence of the “resistant” strain being a mutant of a colicin-sensitive strain.
+
-
                </p>
+
-
 
+
-
                <h2 id="1.7"> Making it Obvious</h2>
+
-
 
+
-
<p>
+
-
      From the graph of our new model (Figure 1, left) it can be deduced that there are paths that converge at stable points (u1,0,0), (0,u2,0) and (0,0,u3), and that this paths all have an approximately common origin. In this section we would like to show that the origin of these paths is practically the same, and that in that sense we have designed a true randomizer (since, as mentioned before, the imprecisions that result in the experimental measurements will make it impossible to make the initial population density of the three strains the same).
+
-
</p>
+
-
 
+
-
                <p>
+
-
In the following set of graphs we will make it obvious that each of the three different strains of E. coli to survive in a random fashion by minimal differences on the initial concentrations of each strain.
+
-
                </p>
+
-
                <p>
+
-
We modeled our results using Matlab.  As can be seen in the graphs below, each of the strains can survive if their initial density in only tree hundredths (a.u.) greater than the other two strains’ initial concentrations.
+
-
                </p>
+
-
<div align= "center">
+
-
                <img src="https://static.igem.org/mediawiki/2011/thumb/1/19/Image130.png/800px-Image130.png" alt="output1" Width="600px" />
+
-
                <img src="https://static.igem.org/mediawiki/2011/1/11/Tt-modeling1.png" alt="coefficient of output1" Width="600px"/>
+
-
                <img src="https://static.igem.org/mediawiki/2011/thumb/a/a1/Image132.png/800px-Image132.png" alt="output2" Width="600px" />
+
-
                <img src="https://static.igem.org/mediawiki/2011/9/99/Tt-modeling2.png" alt="coefficient of output2" Width="600px"/>
+
-
                <img src="https://static.igem.org/mediawiki/2011/thumb/a/a2/Image134.png/800px-Image134.png" alt="output3" Width="600px" />
+
-
                <img src="https://static.igem.org/mediawiki/2011/b/b7/Tt-modeling3.png" alt="coefficient of output3" Width="600px"/>
+
-
</div>
+
-
 
+
-
                <p>
+
-
We this graphs it becomes clear that the imprecisions in experimental measurements (i.e. pipetting) are enough to cause the outcome of Rock, Paper or Scissors signaling molecule to be random. Consequently, we can conclude that this randomizer is not only feasible but also practical and effective (let alone interesting).
+
-
                </p>
+
 +
<p>
 +
From the graph of our new model (Figure 1, left) it can be deduced that there are
 +
paths that converge at stable points (u<sub>1</sub>,0,0), (0,u<sub>2</sub>,0) and (0,0,u<sub>3</sub>), and that this
 +
paths all have an approximately common origin. In this section we would like to show
 +
that the origin of these paths is practically the same, and that in that sense we have
 +
designed a true randomizer (since, as mentioned before, the imprecisions that result
 +
in the experimental measurements will make it impossible to make the initial
 +
population density of the three strains the same).
 +
</p>
 +
<p>
 +
In the following set of graphs we will make it obvious that each of the three different
 +
strains of <span class="name">E. coli</span> to survive in a random fashion by very small
 +
differences on the initial concentrations of each strain.
 +
</p>
 +
<p>
 +
We modeled our results using Matlab. As can be seen in the graphs below, each of the strains
 +
can survive if their initial density in only tree hundredths (a.u.) greater than
 +
the other two strains' initial concentrations.
 +
</p>
 +
 +
<table>
 +
<tr>
 +
<td>
 +
<a href="https://static.igem.org/mediawiki/2011/thumb/1/19/Image130.png/800px-Image130.png">
 +
<img src="https://static.igem.org/mediawiki/2011/thumb/1/19/Image130.png/800px-Image130.png" alt="output1" Width="400px" />
 +
</a>
 +
</td>
 +
<td>
 +
<a href="https://static.igem.org/mediawiki/2011/thumb/a/a1/Image132.png/800px-Image132.png">
 +
<img src="https://static.igem.org/mediawiki/2011/thumb/a/a1/Image132.png/800px-Image132.png" alt="output2" Width="400px" />
 +
</a>
 +
</td>
 +
</tr>
 +
<tr>
 +
<td>
 +
<img src="https://static.igem.org/mediawiki/2011/1/11/Tt-modeling1.png" alt="coefficient of output1" Width="400px"/>
 +
</td>
 +
<td>
 +
<img src="https://static.igem.org/mediawiki/2011/9/99/Tt-modeling2.png" alt="coefficient of output2" Width="400px"/>
 +
</td>
 +
</tr>
 +
<tr>
 +
<td>
 +
<a href="https://static.igem.org/mediawiki/2011/thumb/a/a2/Image134.png/800px-Image134.png">
 +
<img src="https://static.igem.org/mediawiki/2011/thumb/a/a2/Image134.png/800px-Image134.png" alt="output3" Width="400px" />
 +
</a>
 +
</td>
 +
<td rowspan="2">
 +
With these graphs it becomes clear that the imprecisions in experimental measurements
 +
(i.e. pipetting) are enough to cause the outcome of Rock, Paper or Scissors signaling
 +
molecule to be random. Consequently, we can conclude that this randomizer is not
 +
only feasible but also practical and effective (let alone interesting).
 +
</td>
 +
</tr>
 +
<tr>
 +
<td><img src="https://static.igem.org/mediawiki/2011/b/b7/Tt-modeling3.png" alt="coefficient of output3" Width="400px"/></td>
 +
</tr>
 +
</table>
 +
 +
</p>
</p>
 +
<br/>
 +
<br/>
 +
<h2>Reference </h2>
 +
<p>Durrett, R., Levin, S. Allelopathy in Spatially Distributed Populations (1997). Journal of Theoretical Biology, 185, 165-171.</p>
<!-- ############ End of main contents ############ -->
<!-- ############ End of main contents ############ -->
Line 448: Line 585:
</p>
</p>
<!-- end of under contents -->
<!-- end of under contents -->
 +
<div style="bottom: 0;left: 30; position:fixed; background-color: #000000; opacity:0.8;">
 +
<a onclick="scrollTo(0,0);return false;" onfocus="this.blur();" style="color:#FFFFFF;">Return to Page Top</a>
 +
</div>
<!-- DO NOT WRITE UNDER HERE -->
<!-- DO NOT WRITE UNDER HERE -->
</body>
</body>
</html>
</html>

Latest revision as of 14:09, 28 October 2011

Tokyo Tech 2011

Survival of One Strain

1. Introduction: Very small differences determine who will survive

In this section we will show a shocking scenario of evolution: the future of each of three different rival strains (whether the strain will die or survive) is marked by very small differences between the initial population densities of the strains(a phenomenon also known as the "butterfly effect"). Furthermore, we will also show that we can apply this very interesting result to create a randomizer that can be used in our Rock-Paper-Scissors game, due to the fact that only one of the rival strains will survive. More specifically, we assign to each of the three rival strains either of Rock, Paper or Scissors, make them compete for survival and take the surviving strain to represent the bacteria’s choice for the RPS game.

2. Adjusting the Model to create a True Randomizer

The idea for creating this randomizer was born from a paper written in 1996 by Durrett and Levin. In it, the authors described a system of three types of bacteria that competed for survival in dynamic that resembled a Rock-Paper-Scissors (RPS) game. However, the model proposed in this paper is not fully appropriate for our RPS randomizer, since one of the three types of bacteria cannot ultimately survive (although it can dominate the system, i.e. have the highest population density, for definite periods of time). We will discuss more on the limitations we found in this model to be adopted as a randomizer and the modifications we made to create a true randomizer.

3. How the Three Types of Bacteria Compete for Survival

The three types of bacteria that compete for survival use three tactics to outcompete their rivals: the production of a toxin (a bacteriocin called colicin) that is toxic to other strains, resistance to the toxin produced by other strains, and a higher birth rate than their rival strains. Namely, the three types of bacteria are: colicin-producing E. coli (R), colicin-resistant E. coli (P) and colicin-sensitive E. coli (S). The colicin-producer outcompetes the colicin-sensitive by producing the colicin. The colicin-sensitive bacteria outcompetes the colicin-resistant because its birth rate is higher than that of the colicin-resistant. The colicin-resistant outcompetes the colicin producer because its birth rate is higher than that of the colicin producer. The colicin resistant bacteria are also able to produce colicin, but at a lower energetic cost, which allows them to have a higher birth rate.

The system was described by the following general differential equations

Where

4. The Old Model

In the model described by Durrett and Levin’s paper the equations were as follows:

Producer

Resistant

Sensitive

These equations show that the colicin-resistant bacteria are completely immune to colicin (there is not death factor associated to colicin in the equation for du2/dt). However, as will be explained afterwards, this results in a loss of balance that does not allow building a true randomizing system.

Now, setting the parameters as follows, the graph below was created by Durrett and Levin.





5. Our New Model

As mentioned before, the model proposed by Durrett and Levin has critical limitations as a randomizer for the RPS game. To be able to create a true randomizer, we modified the differential equations of the model taking care to give it a biological meaning. With our new differential equations, any of the three types of bacteria can ultimately survive by outcompeting the other two strains, which will die. More specifically, we limited the resistance of the colicin-resistant bacteria in the sense that it would produce a type of bacteriocin that is only toxic to itself and to the sensitive strain, and additionally the resistant strain would also be vulnerable to the colicin produced by the colicin-producer. Since which strain will be the one that survives is determined by very small differences in the initial concentrations of the three different populations of bacteria, in practice this systems becomes a randomizer because of the imprecisions in the measurements that result, for example, when using micropipettes. This randomizer describes a new competition dynamic that could not be reproduced in the previous model proposed by Durrett and Levin due to the instability along the u1axis.

our new model

If we set the parameters as follows

new model's coefficient

and we graph this equations using a Matlab program, we get a graph which clearly shows there are stable points on each of the three axes (Figure 1, Up).


Figure 1. Up: Our New model. Down: The Old Model

These stable points (u1,0,0), (0,u2,0) and (0,0,u3) indicate that for the equations we have set all of the three strains may ultimately survive for infinite peiriods of time. The differences between our model and the model of Durrett and Levin can be seen graphically in Figure 1. These graphs were plotted using Matlab.

Note that the parameters we have set for our equations satisfy the initail conditions of the model proposed by Durrett and Levin (indicated in black font)
new terms

6. The Biological Meaning of our Model

From a biological perspective, our model describes the existence of two strains of bacteria that produce two different types of bacteriocins. One of these strains is not completely resistant to its own bacteriocin nor to the bacteriocin produced by its rival strain. This can be justified as the consequence of insufficient/ineffective resistance protein production by the “resistant” strain. This limitation in the production of resistance protein could be thought of as a consequence of the “resistant” strain being a mutant of a colicin-sensitive strain.

7. Making it Obvious

From the graph of our new model (Figure 1, left) it can be deduced that there are paths that converge at stable points (u1,0,0), (0,u2,0) and (0,0,u3), and that this paths all have an approximately common origin. In this section we would like to show that the origin of these paths is practically the same, and that in that sense we have designed a true randomizer (since, as mentioned before, the imprecisions that result in the experimental measurements will make it impossible to make the initial population density of the three strains the same).

In the following set of graphs we will make it obvious that each of the three different strains of E. coli to survive in a random fashion by very small differences on the initial concentrations of each strain.

We modeled our results using Matlab. As can be seen in the graphs below, each of the strains can survive if their initial density in only tree hundredths (a.u.) greater than the other two strains' initial concentrations.

output1 output2
coefficient of output1 coefficient of output2
output3 With these graphs it becomes clear that the imprecisions in experimental measurements (i.e. pipetting) are enough to cause the outcome of Rock, Paper or Scissors signaling molecule to be random. Consequently, we can conclude that this randomizer is not only feasible but also practical and effective (let alone interesting).
coefficient of output3



Reference

Durrett, R., Levin, S. Allelopathy in Spatially Distributed Populations (1997). Journal of Theoretical Biology, 185, 165-171.

Return to Page Top