Difference between revisions of "Team:Tokyo Tech/Model"

Line 398: Line 398:
 
</div> <!-- population_contents-->
 
</div> <!-- population_contents-->
 
</div> <!-- population-->
 
</div> <!-- population-->
 +
 +
<div onclick="obj=document.getElementById('open3').style; obj.display=(obj.display=='none')?'block':'none';">
 +
<a style="cursor:pointer;">▼Read more</a>
 +
</div>
 +
<!--// 折り畳み展開ポインタ -->
 +
<!-- 折り畳まれ部分 -->
 +
<div id="open3" style="display:none;clear:both;">
  
 
<div id="toxin">
 
<div id="toxin">
Line 507: Line 514:
 
    </div> <!--more_contents-->
 
    </div> <!--more_contents-->
 
    </div><!--more-->
 
    </div><!--more-->
 
+
</div>
 
</div><!-- /fitting_contents -->
 
</div><!-- /fitting_contents -->
 
</div><!-- /fitting -->
 
</div><!-- /fitting -->

Revision as of 17:01, 19 October 2016

1. Overview

To recreate the story of ”Snow White”, we have designed a cell-cell communication system with improved or characterized parts and collected data from comprehensive experiments. Furthermore, we constructed the mathematical model to simulate the behavior of the whole system and to confirm the feasibility of our story. This simulation successfully contributed to give the suggestions to wet lab experiments. In addition, in order to help us utilize our Toxin-Antitoxin (TA) system, we developed a new software in Java for adjusting the number of ACA sequences, which MazF dimer recognizes and cleaves in mRNAs.

2. Story simulation

2-1. Mathematical model

In order to simulate our gene circuits, we developed an ordinary differential equation model.

[Model development]


2-2. Results

We obtained and confirmed the desirable behavior of the whole system by modifying and improving parts. As described below, our simulation showed an appropriate transition of concentration of RFP and GFP for the story.


Fig.5-2. Time-dependent change of the concentrations of fluorescence proteins


In the blue area of Fig.5-2-2, the concentration of fluorescence proteins starts to increase. The concentration of RFP of Snow White coliexceeds the concentration of GFP of the Queen coli.
It is as if Snow White got fairer more and more.

In the pink area of Fig.5-2-2, C12 is being synthesized and increased thanks to the appearance of C4. As a result, the MazF inside Snow White coli and the Queen coli, suppressing the increment of fluorescent proteins.
It is as if the Queen, influenced by the Mirror's answer, transforming into a Witch in order to give Snow White a poisoned apple.

In the green area of Fig.5-2-2, C12 increases and the MazF inside Snow White coli induced by it increases even more. So the concentration GFP exceeds that of RFP.
It looks as if Snow White bit the apple, sinking into unconsciousness promptly.

In the yellow area of Fig.5-2-2, the AmiE synthesized by the introduced the Prince coli decomposes C12 so the MazF inside Snow White coli diminishes and C4 increases.
It looks as if the Prince lifted Snow White and she opened her eyes.

3. Fitting

3-1. Population growth

First, we tried to model the growth curve of the system. When the number of E. coli approaches a certain value, the growth will stop. We defined this value in the culture as Pmax. Then the population growth equation for our system is described as follows:
$$ \frac{dP}{dt} = g\left(1 - \frac{P}{P_{max}}\right)P$$
where g is the population growth rate.
This equation can be analytically solved as:
$$ P = \frac{P_{0} P_{max} e^{gt}}{P_{max} - P_{0} + P_{0} e^{gt}}$$
where P 0 is the population at t = 0. We used this equation to fit the experimental data.


Fig.5-3-1. Modeled growth curve of E. coli fitted to experiment data

Using the experimental data from the Toxin assay for this fitting, we estimated the following parameters:

g = 0.0123
and
Pmax =3.3

respectively.
These parameters can be used for Snow White coli, the Queen coli and the Prince coli in the same way.

4. Sensitivity analysis

We performed the sensitivity analysis descried in this section in order to examine which parameter dominates the story.


4-1. Requirements

We defined these requirements as the “successful Snow White story.”
1) At t = 150
     concentration of RFP > concentration of GFP
2) At t = 700
      concentration of RFP< concentration of GFP
3) At t = 1500
     concentration of RFP > concentration of GFP

We the began to analyze the graph that satisfies these requirements that we could say that recreates the story correctly. In the table bellow we show the range in which we modified each parameter. They were modified one step size at a time.

Parameter Range Step size
D $$ 0.0001 < D < 0.001 $$ 0.0001
pC4 $$ 0.0001 < p_{C4} < 1 $$ 0.001
pC12 $$ 0.0001 < p_{C12} < 1 $$ 0.001
α $$ 0.01 < α < 0.2 $$ 0.01
dAmiE $$ $$ 0.001 < d_{AmiE} < 1 $$ 0.001

As a result, we obtained the following parameter ranges.

Parameter Value
D $$ 0.0056 < D < 0.001 $$
pC4 $$ 0.0029 < p_{C4} < 0.778 $$
pC12 $$ 0.001 < p_{C12} < 0.217 $$
α $$ 0.01 < α < 0.16 $$
dAmiE $$ 0.001 < d_{AmiE} < 1 $$


4-2. the Prince coli should be put in during the process

We run simulations in order to determine whether we would get a better behavior let we introduce the Prince coli at the beginning or halfway of the story.

Fig.5-4-2-1. Number of individuals when the Prince coli is introduced from the beginning

Fig.5-4-2-2. AHL concentrations when the Prince coli is introduced from the beginning

As a result, if we introduce the Prince coli from the beginning, the number of Prince coli increases too much (Fig.5-4-2-1) so the AmiE the Prince coli produces augments and the decomposition of C12 also occurs overly. So C12 is almost inexistent in the medium (Fig.5-4-2-2).


Fig.5-4-2-3. Number of individuals when the Prince coli is introduced at t = 700


Fig.5-4-2-4. AHL concentrations when the Prince coli is introduced at t = 700

On the other hand, if we introduce the Prince coli at t = 600, the number of Prince coli does not increment much (Fig.5-4-2-3), so C12 can exist until t = 600 and then decreases thanks to the augment of AmiE (Fig.5-4-2-4).
In conclusion, if we introduce the Prince coli at t = 700, the circuit will behave accordingly.


4-3. Prhl should be changed

In order to confirm the feasibility of the story with our gene circuit by the combination of the existing promoters, we performed some simulations based on the results of our assays.


Fig.5-4-3. The intensity of Plux and Prhl promoters

The diagram above shows that the intensity of the two promoters should be in the red region of the figure.The combination of promoters which we were originally going to use is shown in the graph by the green point. To move this point into the red region, we had to improve Prhl to raise its expression level.


4-4. Production rate of C4HSL and 3OC12HSL by RhlI and LasI

The signaling molecule production rates by RhlI and LasI can be changed by modifying RhlI and LasI to make more or less signaling molecule in silico.


Fig.5-4-4-1. Concentrations of GFP and RFP dependencies of production rate of C4 by RhlI

Each line corresponds to the transition of the concentration of RFP and GFP with a certain production rate of C12. Red and green lines correspond to RFP and GFP, respectively. Blighter one indicates higher production rate of C12.
If the production rate of C4 is between 0.029 and 0.778, our system can recreate the story. If this parameter is too small, the production of LasI by the Queen coli is insufficiently inhibited by MazF so C12 increases greatly. As a result, The concentration of GFP overcomes the concentration of RFP and the story does not develop correctly. And if this parameter is too big, the production of LasI by the Queen coli is overly inhibited by MazF so C12 does not increase.
As a result, the concentration of RFP is always greater than the GFP concentration and the story does not develop either.


Fig.5-4-4-2. Concentrations of GFP and RFP dependencies of production rate of C12 by LasI

Each line corresponds to the transition of the concentration of RFP/GFP with a certain production rate of C4. Red and green lines correspond to RFP and GFP, respectively. Blighter one indicates higher production rate of C4.
If production rate of C12 is between 0.001 and 0.217, our system can recreate the story.


4-5. Translation rate of protein

Translation rate affects the production of protein.


Fig.5-4-5. Concentrations of GFP and RFP dependencies of translation rate of proteins

Each line corresponds to the transition of the concentration of RFP and GFP with a certain translation rate. Red and green lines correspond to RFP and GFP, respectively. Blighter one indicates higher translation rate.
If translation rate of protein is between 0.01 and 0.16, our system can recreate the story.


4-6. Decomposition rate of C12 by AmiE

AmiE decomposes C12. The concentration of C12 after input of the Prince coli is changed by AmiE. If the decomposition rate of C12 is too small, C12 does not decrease enough so the MazF inside Snow White coli continues being expressed and the concentration of RFP decreases.


Fig.5-4-6. Concentrations of GFP and RFP dependencies of decomposition rate of C12 by AmiE

Each line corresponds to the transition of the concentration of RFP and GFP with a certain decomposition rate of C12. Red and green lines correspond to RFP and GFP, respectively. Blighter one indicates higher value of decomposition rate of C12.
If degradation rate of C12 is higher than 0.00056, our system can recreate the story.


4-7. Degradation rate of AmiE


Fig.5-4-7. Concentration of GFP and RFP dependencies of degradation rate of AmiE

Each line corresponds to the transition of the concentration of RFP and GFP with a certain degradation rate of AmiE. Red and green lines correspond to RFP and GFP, respectively. Blighter one indicates higher degradation rate of AmiE.
Even if we modify the values of the parameters inside the defined range, the concentration of RFP overcomes the concentration of GFP. And even if the degradation rate of AmiE is small, the decomposition rate of C12 by AmiE is high enough so C12 decreases sufficiently.


4-8. Degradation rate of RFP and GFP

The degradation rate of RFP and GFP is key to the success the story of ‘"Snow White".
These parameters are closely related to the concentrations of GFP and RFP, so we conjectured that if they do not take appropriate values the story can not be correctly recreated.


Fig.5-4-8 Relation of degradation rate of GFP and RFP

The story is recreated only if the degradation rate of RFP and GFP are the same.

5. Software


5-1. Abstract


We developed a new software named ACADwarfs. This software helps to control the sensitivity of the protein to MazF by regulating the number of ACA sequences in the mRNA sequence. ACADwarfs can increase or decrease the number of ACA sequences on mRNA without changing the amino acid sequences that the mRNA specifies or frameshifts resulted from insertion of bases without considering.
Then we improve the practicality of the characteristic of the mazEF system. For example you can let protein A express constantly by eliminating ACA sequences of the sequence, while letting protein B stop being expressed, at the desirable timing, by expression of MazF.
This software also evades the use of rare codons, so you don’t have to worry about them.


5-2. Key achievements

・Provided the tool regulating the number of ACA sequences

・Released under open-source license so everyone can use it

・Able to correspond to any base arrangements

・Rare codons are evaded

・Extend the application field of mazEF system


5-3. Work flow

  1. Start ACADwarfs

    Download the zip file and open it.

  2. Put in the part sequence

    Put the part sequence you want to regulate the effect of the mazEF system on in the upper window.

  3. Choose operation

    If you want to magnify the effect of mazEF system, choose “increase ACA”.

    If you want to lessen the effect of mazEF system, choose “decrease ACA”.

  4. Get sequences

    You can get the modified part sequences in the lower window. Modified codons will be written in capital letters and while the rest will be in small letters so you can easily locate the modified parts.
    And under this window you can see the number of ACA sequences on the sequence before the modification of the following "previous number" and the number of ACA sequences on the sequence after the modification of the following "new number".


5-4. Demonstration

We created a demo to present the features of this software. Using this, we regulated the number of ACA sequences and control the sensitivity of the protein to MazF.

Gene Pre number of ACA sequences Post number of ACA sequences
RFP 10 30
GFP 23 39
MazF 2 1
MazE 2 1

We increased the numbers of ACA sequences of RFP and GFP decreased the numbers of ACA sequences of MazF and MazE .


Fig.5-4. Comparison between the results of simulations using original sequences and modified sequences

We can see that the concentration of expressed MazF reacts more keenly after adjusting the ACA sequences than before doing so.


5-5. Download

To download click here.

The code is available on github.