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

From 2011.igem.org

(Difference between revisions)
 
(30 intermediate revisions not shown)
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 235: 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" heightT="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="movie" value="https://static.igem.org/mediawiki/2011/4/4d/Header.swf">
<param name="quality" value="high">
<param name="quality" value="high">
Line 291: 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><a href="#1.abst">1. Abstract</a></li>
+
<li><a href="#1.1">1. Introduction</a></li>
 +
<li><a href="#1.2">2. Adjusting the Model to create a True Randomizer</a></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 310: Line 341:
<div align="center">
<div align="center">
-
<img src="https://static.igem.org/mediawiki/2011/2/28/Image042.png" />
+
<img src="https://static.igem.org/mediawiki/2011/2/28/Image042.png" width="500px" />
</div>
</div>
               
               
-
<h2 id="1.1">Introduction: Minimal differences determine who will survive</h2>
+
<h2 id="1.1">1. Introduction: Very small differences determine who will survive</h2>
<p>
<p>
In this section we will show a shocking scenario of evolution:  
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)  
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.  
+
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  
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  
randomizer that can be used in our Rock-Paper-Scissors game, due to the fact that only one  
Line 326: Line 357:
</p>
</p>
               
               
-
<h2 id="1.2">Adjusting the Model to create a True Randomizer</h2>
+
<h2 id="1.2">2. Adjusting the Model to create a True Randomizer</h2>
<p>
<p>
-
The idea for creating this randomizer was born from a paper written in 1996 by Durret and Levin.  
+
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  
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,  
survival in dynamic that resembled a Rock-Paper-Scissors (RPS) game. However,  
Line 339: Line 370:
</p>
</p>
                 
                 
-
<h2 id="1.3"> How the Three Types of Bacteria Compete for Survival</h2>
+
<h2 id="1.3">3. How the Three Types of Bacteria Compete for Survival</h2>
<p>
<p>
Line 358: Line 389:
<div align="center">
<div align="center">
-
<img src="https://static.igem.org/mediawiki/2011/b/bf/Image044.png" />
+
<img src="https://static.igem.org/mediawiki/2011/b/bf/Image044.png" width="500px" />
</div>
</div>
Line 378: Line 409:
-
<h2 id="1.4">The Old Model</h2>
+
<h2 id="1.4">4. The Old Model</h2>
<p>
<p>
-
In the model described by Durret and Levin’s paper the equations were as follows:
+
In the model described by Durrett and Levin’s paper the equations were as follows:
</p>
</p>
<div align="center">
<div align="center">
-
<p>Producer</p>
+
Producer<br />
<img src="https://static.igem.org/mediawiki/2011/1/1d/Colisin3.png" /><br />
<img src="https://static.igem.org/mediawiki/2011/1/1d/Colisin3.png" /><br />
-
<p>Resistant</p>
+
Resistant<br />
<img src="https://static.igem.org/mediawiki/2011/0/04/Colisin4.png" /><br />
<img src="https://static.igem.org/mediawiki/2011/0/04/Colisin4.png" /><br />
-
<p>Sensitive</p>
+
Sensitive<br />
<img src="https://static.igem.org/mediawiki/2011/3/36/Colisin5.png" /><br />
<img src="https://static.igem.org/mediawiki/2011/3/36/Colisin5.png" /><br />
</div>
</div>
Line 399: Line 430:
<p>
<p>
-
Now, setting the parameters as follows, the graph below was created by Durret and Levin.
+
Now, setting the parameters as follows, the graph below was created by Durrett and Levin.
</p>
</p>
Line 405: Line 436:
<img src="https://static.igem.org/mediawiki/2011/4/41/Colisin6.png" /><br />
<img src="https://static.igem.org/mediawiki/2011/4/41/Colisin6.png" /><br />
<img src="https://static.igem.org/mediawiki/2011/e/e3/Colisin7.png" />
<img src="https://static.igem.org/mediawiki/2011/e/e3/Colisin7.png" />
-
<img src="https://static.igem.org/mediawiki/2011/thumb/e/e5/Image086.png/669px-Image086.png" width="354px"/><br />
+
<img src="http://partsregistry.org/wiki/images/b/bd/Modeling1.png" width="354px"/><br />
<img src="https://static.igem.org/mediawiki/2011/thumb/b/b1/Colisin8.png/800px-Colisin8.png" width="800px"/><br />
<img src="https://static.igem.org/mediawiki/2011/thumb/b/b1/Colisin8.png/800px-Colisin8.png" width="800px"/><br />
-
<img src="https://static.igem.org/mediawiki/2011/b/b8/Image094.png" width="354px"/><br />
+
<img src="https://static.igem.org/mediawiki/2011/5/55/Modeling2.png" width="354px"/><br />
</div>
</div>
-
<h2 id="1.5"> Our New Model</h2>
+
<h2 id="1.5">5. Our New Model</h2>
<p>
<p>
-
As mentioned before, the model proposed by Durret and Levin has critical limitations as a randomizer for the RPS game.  
+
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  
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  
biological meaning. With our new differential equations, any of the three types of bacteria can ultimately survive by  
Line 421: Line 452:
the three different populations of bacteria, in practice this systems becomes a randomizer because of the imprecisions in  
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  
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
+
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.
<img src="https://static.igem.org/mediawiki/2011/e/e0/Image088.png" alt="u1">axis.
</p>
</p>
Line 445: Line 476:
<p>
<p>
-
These stable points (u1,0,0), (0,u2,0) and (0,0,u3) indicate that for the equations we have set  
+
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  
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.  
+
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>
<p>
Note that the parameters we have set for our equations satisfy the initail  
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)<br />
+
conditions of the model proposed by Durrett and Levin (indicated in black font)<br />
<img src="https://static.igem.org/mediawiki/2011/a/a4/Colimodel4.png" alt="new terms" width="800px"/>
<img src="https://static.igem.org/mediawiki/2011/a/a4/Colimodel4.png" alt="new terms" width="800px"/>
</p>
</p>
-
<h4 id="1.6"> The Biological Meaning of our Model</h4>
+
<h2 id="1.6">6. The Biological Meaning of our Model</h2>
<p>
<p>
Line 468: Line 499:
</p>
</p>
-
<h2 id="1.7">Making it Obvious</h2>
+
<h2 id="1.7">7. Making it Obvious</h2>
<p>
<p>
From the graph of our new model (Figure 1, left) it can be deduced that there are  
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 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  
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  
that the origin of these paths is practically the same, and that in that sense we have  
Line 482: Line 513:
<p>
<p>
In the following set of graphs we will make it obvious that each of the three different  
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 minimal
+
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.
differences on the initial concentrations of each strain.
</p>
</p>
Line 492: Line 523:
</p>
</p>
-
<div align= "center">
+
<table>
-
                <img src="https://static.igem.org/mediawiki/2011/thumb/1/19/Image130.png/800px-Image130.png" alt="output1" Width="600px" />
+
<tr>
-
                <img src="https://static.igem.org/mediawiki/2011/1/11/Tt-modeling1.png" alt="coefficient of output1" Width="600px"/>
+
<td>
-
                <img src="https://static.igem.org/mediawiki/2011/thumb/a/a1/Image132.png/800px-Image132.png" alt="output2" Width="600px" />
+
<a href="https://static.igem.org/mediawiki/2011/thumb/1/19/Image130.png/800px-Image130.png">
-
                <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/1/19/Image130.png/800px-Image130.png" alt="output1" Width="400px" />
-
                <img src="https://static.igem.org/mediawiki/2011/thumb/a/a2/Image134.png/800px-Image134.png" alt="output3" Width="600px" />
+
</a>
-
                <img src="https://static.igem.org/mediawiki/2011/b/b7/Tt-modeling3.png" alt="coefficient of output3" Width="600px"/>
+
</td>
-
</div>
+
<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>
 
-
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>
</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 524: 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