Line 54: | Line 54: | ||
// padding: 18px; | // padding: 18px; | ||
font-family: 'Century', serif; | font-family: 'Century', serif; | ||
+ | } | ||
+ | #contents_menu{ | ||
+ | line-height: 1em; | ||
} | } | ||
Line 68: | Line 71: | ||
} | } | ||
− | |||
− | |||
− | |||
--></style> | --></style> | ||
<script src="https://2015.igem.org/common/MathJax-2.5-latest/MathJax.js?config=TeX-AMS-MML_HTMLorMML"></script> | <script src="https://2015.igem.org/common/MathJax-2.5-latest/MathJax.js?config=TeX-AMS-MML_HTMLorMML"></script> | ||
Line 84: | Line 84: | ||
<div id="main_contents"> | <div id="main_contents"> | ||
<div id="page_header" class="container container_top"> | <div id="page_header" class="container container_top"> | ||
− | <h1 align="center"> | + | <h1 align="center">Model</h1> |
</div><!-- /page_header --> | </div><!-- /page_header --> | ||
<div id="contents" class="container"> | <div id="contents" class="container"> | ||
Line 92: | Line 92: | ||
<div id="contents_contents" class="container_contents"> | <div id="contents_contents" class="container_contents"> | ||
<div id="contents_menu"> | <div id="contents_menu"> | ||
− | <h3 class="link"><a href="# | + | <h3 class="link"><a href="#overview">1. Overview</a></h3> |
− | <h3 class="link"><a href="# | + | <h3 class="link"><a href="#story_simulation">2. Story Simulation</a></h3> |
− | <h3 class="link"><a href="# | + | <h3 class="link"><a href="#mathematical_model"><font size="2.7"> 2-1. Mathematical Model</font></a></h3> |
− | <h3 class="link"><a href="# | + | <h3 class="link"><a href="#results"><font size="2.7"> 2-2. Results</font></a></h3> |
− | <h3 class="link"><a href="# | + | <h3 class="link"><a href="#fitting">3. Fitting</a></h3> |
− | <h3 class="link"><a href="# | + | <h3 class="link"><a href="#analysis">4. Sensitivity Analysis</a></h3> |
− | <h3 class="link"><a href="# | + | <h3 class="link"><a href="#requirements"><font size="2.7"> 4-1. Requirements</font></a></h3> |
− | <h3 class="link"><a href=" | + | <h3 class="link"><a href="#prince_coli"><font size="2.7"> 4-2. the Prince <span style ="font-style : italic">coli</span> should be put in in the middle</font></a></h3> |
− | <h3 class="link"><a href=" | + | <h3 class="link"><a href="prhl"><font size="2.7"> 4-3. Prhl should be changed</font></a></h3> |
− | <h3 class="link"><a href=" | + | <h3 class="link"><a href="production_ahl"><font size="2.7"> 4-4. Production rate of C4HSL and 3OC12HSL by RhlI and LasI</font></a></h3> |
− | <h3 class="link"><a href=" | + | <h3 class="link"><a href="translation"><font size="2.7"> 4-5. Translation rate of protein</font></a></h3> |
− | <h3 class="link"><a href=" | + | <h3 class="link"><a href="decomposition"><font size="2.7"> 4-6. Decomposition rate of O3C12HSL by AmiE</font></a></h3> |
− | <h3 class="link"><a href=" | + | <h3 class="link"><a href="degradation_amie"><font size="2.7"> 4-6. Degradation rate of AmiE</font></a></h3> |
− | <h3 class="link"><a href="# | + | <h3 class="link"><a href="degradation_protein"><font size="2.7"> 4-7. Degradation rate of RFP and GFP</font></a></h3> |
− | <h3 class="link"><a href="# | + | <h3 class="link"><a href="#software">5. Software</a></h3> |
− | <h3 class="link"><a href="# | + | <h3 class="link"><a href="#abstract"><font size="2.7"> 5-1. Abstract</font></a></h3> |
− | <h3 class="link"><a href="# | + | <h3 class="link"><a href="#key_achievement"><font size="2.7"> 5-2. Key achievement</font></a></h3> |
− | <h3 class="link"><a href="# | + | <h3 class="link"><a href="#work_flow"><font size="2.7"> 5-3. Work flow</font></a></h3> |
+ | <h3 class="link"><a href="#demo"><font size="2.7"> 5-4. Demonstration</font></a></h3> | ||
+ | <h3 class="link"><a href="#download"><font size="2.7"> 5-5. Download</font></a></h3> | ||
</div><!-- /contents_menu --> | </div><!-- /contents_menu --> | ||
</div><!-- /contents_contents --> | </div><!-- /contents_contents --> | ||
</div><!-- /contents --> | </div><!-- /contents --> | ||
+ | <div id="overview" class="container"> | ||
+ | <div id="overview_header" class="container_header"> | ||
+ | <h2><span>1. Overview</span></h2> | ||
+ | </div><!-- /overview_header --> | ||
+ | <div id="overview_contents" class="container_contents"> | ||
+ | <p class="normal_text">To reproduce the story of ”Snow White”, we have designed the cell-cell communication system with improved or characterized parts and collected data from comprehensive experiments. Furthermore, we constructed the mathematical model in order 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 and confirm the feasibility of our system. In addition, in order to utilize TA (Toxin-Antitoxin) system, we developed a new software in Java for adjusting the number of ACA sequences, which MazF dimer recognizes and cleaves in mRNAs.</p> | ||
+ | </div><!-- /overview_contents --> | ||
+ | </div><!-- /overview --> | ||
− | <div id=" | + | <div id="story_simulation" class="container"> |
− | <div id=" | + | <div id="story_simulation_header" class="container_header"> |
− | <h2><span> | + | <h2><span>2. Story Simulation</span></h2> |
− | </div><!-- / | + | </div><!-- /story_simulation_header --> |
− | <div id=" | + | <div id="story_simulation_contents" class="container_contents"> |
− | <p class="normal_text"> | + | <div id="mathematical_model"> |
− | + | <div id="mathematical_model_header"> | |
− | + | <h3><span>2-1. Mathematical Model</span></h3> | |
− | + | </div><!-- /_header --> | |
− | + | <p class="normal_text">To simulate our gene circuits, we developed an ordinary differential equation model.</p> | |
− | + | <p style="text-align:center;">[<a href="https://2016.igem.org/Team:Tokyo_Tech/Modeling_Details">Detailed Description for Modeling</a>]</p> | |
+ | <p class="normal_text" style="text-align:center;"><a href="javascript:void(0);" onClick="show('modeling_detail');" class="showHidden">Read More</a></p> | ||
+ | <div id="modeling_detail" class="off"> | ||
+ | <div id="modeling_detail_wrapper"> | ||
+ | <div id="modeling_detail_expressions"> | ||
+ | <h2>Differencial Equations</h2> | ||
+ | <h3>Snow White</h3> | ||
+ | \begin{equation} | ||
+ | \frac{d[mRNA_{RFP}]}{dt} = k - d[mRNA_{RFP}] - F_{DiMazF}(1-(1-f)^{f_{mRNA_{RFP}}})[mRNA_{RFP}][DiMazF] | ||
+ | \end{equation} | ||
+ | \begin{equation} | ||
+ | \frac{d[mRNA_{RhlI}]}{dt} = leak_{P_{lux}} + \frac{\kappa_{Lux}[C12]^{n_{Lux}}}{K_{mLux}^{n_{Lux}} + [C12]^{n_{Lux}}} - d[mRNA_{RhlI}] - F_{DiMazF}f[mRNA_{RhlI}][DiMazF] | ||
+ | \end{equation} | ||
+ | \begin{equation} | ||
+ | \frac{d[RFP]}{dt} = \alpha [mRNA_{RFP}] - d_{RFP}[RFP] | ||
+ | \end{equation} | ||
+ | \begin{equation} | ||
+ | \frac{d[RhlI]}{dt} = \alpha [mRNA_{RhlI}] - d_{RhlI}[RhlI] | ||
+ | \end{equation} | ||
+ | \begin{equation} | ||
+ | \frac{d[C4]}{dt} = p_{C4}[RhlI]P_{Snow White} - d_{C4}[C4] | ||
+ | \end{equation} | ||
+ | \begin{equation} | ||
+ | \frac{d[mRNA_{MazF}]}{dt} = leak_{P_{lux}} + \frac{\kappa_{Lux}[C12]^{n_{Lux}}}{K_{mLux}^{n_{Lux}}+ [C12]^{n_{Lux}}} \\ | ||
+ | - d[mRNA_{MazF}] - F_{DiMazF}(1-(1-f)^{f_{mRNA_{MazF}}})[mRNA_{MazF}][DiMazF] | ||
+ | \end{equation} | ||
+ | \begin{equation} | ||
+ | \frac{d[mRNA_{MazE}]}{dt} = k - d[mRNA_{MazE}] - F_{DiMazF}(1-(1-f)^{f_{mRNA_{MazE}}})[mRNA_{MazE}][DiMazF] | ||
+ | \end{equation} | ||
+ | \begin{equation} | ||
+ | \frac{d[MazF]}{dt} = \alpha [mRNA_{MazF}] - 2k_{Di_{MazF}}[MazF] + 2k_{-Di_{MazF}}[DiMazF] - d_{MazF}[MazF] | ||
+ | \end{equation} | ||
+ | \begin{equation} | ||
+ | \frac{d[DiMazF]}{dt} = k_{Di_{MazF}}[MazF] - k_{-Di_{MazF}}[DiMazF] - 2k_{Hexa}[DiMazE][DiMazF]^2 \\ | ||
+ | + 2k_{-Hexa}[MazHexamer] - d_{DiMazF}[DiMazF] \end{equation} | ||
+ | \begin{equation} | ||
+ | \frac{d[MazE]}{dt} = \alpha [mRNA_{MazE}] - 2k_{Di_{MazE}}[MazE] + 2k_{-Di_{MazE}}[DiMazE] - d_{MazE}[MazE] | ||
+ | \end{equation} | ||
+ | \begin{equation} | ||
+ | \frac{d[DiMazE]}{dt} = k_{Di_{MazE}}[MazE] - k_{-Di_{MazE}}[DiMazE] - k_{Hexa}[DiMazE][DiMazF]^2 \\ | ||
+ | + k_{-Hexa}[MazHexamer] - d_{DiMazE}[DiMazE] | ||
+ | \end{equation} | ||
+ | \begin{equation} | ||
+ | \frac{d[MazHexa]}{dt} = k_{Hexa}[DiMazE][DiMazF]^2 - k_{-Hexa}[MazHexa] - d_{Hexa}[MazHexa] | ||
+ | \end{equation} | ||
+ | \begin{equation} | ||
+ | \frac{dP_{Snow White}}{dt} = g \frac{E_{DiMazF}}{E_{DiMazF}+[DiMazF]}\left(1- \frac{P_{Snow White}+P_{Queen}+P_{Prince}}{P_{max}} \right) P_{Snow White} | ||
+ | \end{equation} | ||
+ | <h3>Queen</h3> | ||
+ | \begin{equation} | ||
+ | \frac{d[mRNA_{GFP}]}{dt} = k - d[mRNA_{GFP}] - F_{DiMazF}(1-(1-f)^{f_{mRNA_{GFP}}})[mRNA_{GFP}][DiMazF] | ||
+ | \end{equation} | ||
+ | \begin{equation} | ||
+ | \frac{d[mRNA_{LasI}]}{dt} = leak_{P_{rhl}} + \frac{\kappa_{Rhl}[C4]^{n_{Rhl}}}{K_{mRhl}^{n_{Rhl}} + [C4]^{n_{Rhl}}} \\ | ||
+ | - d[mRNA_{LasI}] - F_{DiMazF}(1-(1-f)^{f_{mRNA_{LasI}}})[mRNA_{LasI}][DiMazF] | ||
+ | \end{equation} | ||
+ | \begin{equation} | ||
+ | \frac{d[GFP]}{dt} = \alpha [mRNA_{GFP}] - d_{GFP}[GFP] | ||
+ | \end{equation} | ||
+ | \begin{equation} | ||
+ | \frac{d[LasI]}{dt} = \alpha [mRNA_{LasI}] - d_{LasI}[LasI] | ||
+ | \end{equation} | ||
+ | \begin{equation} | ||
+ | \frac{d[C12]}{dt} = p_{C12}[LasI]P_{Queen} - d_{C12}[C12] - D[C12][AmiE] | ||
+ | \end{equation} | ||
+ | \begin{equation} | ||
+ | \frac{d[mRNA_{MazF}]}{dt} = leak_{P_{lux}} + \frac{\kappa_{Rhl}[C4]^{n_{Rhl}}}{K_{mRhl}^{n_{Rhl}} + [C4]^{n_{Rhl}}} \\ | ||
+ | - d[mRNA_{MazF}] - F_{DiMazF}(1-(1-f)^{f_{mRNA_{MazF}}})[mRNA_{MazF}][DiMazF] | ||
+ | \end{equation} | ||
+ | \begin{equation} | ||
+ | \frac{d[mRNA_{MazE}]}{dt} = k - d[mRNA_{MazE}] - F_{DiMazF}(1-(1-f)^{f_{mRNA_{MazE}}})[mRNA_{MazE}][DiMazF] | ||
+ | \end{equation} | ||
+ | \begin{equation} | ||
+ | \frac{d[MazF]}{dt} = \alpha [mRNA_{MazF}] - 2k_{Di_{MazF}}[MazF] + 2k_{-Di_{MazF}}[DiMazF] - d_{MazF}[MazF] | ||
+ | \end{equation} | ||
+ | \begin{equation} | ||
+ | \frac{d[DiMazF]}{dt} = k_{Di_{MazF}}[MazF] - k_{-Di_{MazF}}[DiMazF] - 2k_{Hexa}[DiMazE][DiMazF]^2 \\ | ||
+ | + 2k_{-Hexa}[MazHexamer] - d_{DiMazF}[DiMazF] | ||
+ | \end{equation} | ||
+ | \begin{equation} | ||
+ | \frac{d[MazE]}{dt} = \alpha [mRNA_{MazE}] - 2k_{Di_{MazE}}[MazE] + 2k_{-Di_{MazE}}[DiMazE] - d_{MazE}[MazE] | ||
+ | \end{equation} | ||
+ | \begin{equation} | ||
+ | \frac{d[DiMazE]}{dt} = k_{Di_{MazE}}[MazE] - k_{-Di_{MazE}}[DiMazE] - k_{Hexa}[DiMazE][DiMazF]^2 \\ | ||
+ | + k_{-Hexa}[MazHexamer] - d_{DiMazE}[DiMazE] | ||
+ | \end{equation} | ||
+ | \begin{equation} | ||
+ | \frac{d[MazHexa]}{dt} = k_{Hexa}[DiMazE][DiMazF]^2 - k_{-Hexa}[MazHexa] - d_{Hexa}[MazHexa] | ||
+ | \end{equation} | ||
+ | \begin{equation} | ||
+ | \frac{dP_{Queen}}{dt} = g \frac{E_{DiMazF}}{E_{DiMazF}+[DiMazF]}\left(1- \frac{P_{Snow White}+P_{Queen}+P_{Prince}}{P_{max}}\right) P_{Queen}\\ | ||
+ | \end{equation} | ||
+ | <h3>Prince</h3> | ||
+ | \begin{equation} | ||
+ | \frac{d[mRNA_{AmiE}]}{dt} = leak_{P_{lux}} + \frac{\kappa_{Lux}[C12]^n}{K_{mLux}^n + [C12]^n} - d[mRNA_{AmiE}] | ||
+ | \end{equation} | ||
+ | \begin{equation} | ||
+ | \frac{d[AmiE]}{dt} = \alpha [mRNA_{AmiE}]P_{Prince} - d_{AmiE}[AmiE] | ||
+ | \end{equation} | ||
+ | \begin{equation} | ||
+ | \frac{dP_{Prince}}{dt} = g\left(1- \frac{P_{Snow White}+P_{Queen}+P_{Prince}}{P_{max}}\right) P_{Prince} | ||
+ | \end{equation} | ||
+ | </div><!-- /modeling_detail_expressions --> | ||
+ | <div id="modeling_detail_parameter"> | ||
+ | <h2>Explanation about Parameters</h2> | ||
+ | <table border="1" style="margin: auto;"> | ||
+ | <tbody> | ||
+ | <tr><td>Parameter </td><td> Description </td></tr> | ||
+ | <tr><td>$$g$$ </td><td> Growth rate of each cells</td></tr> | ||
+ | <tr><td>$$P_{max}$$ </td><td> Carrying capacity </td></tr> | ||
+ | <tr><td>$$E_{DiMazF}$$ </td><td> Effect of MazF dimer on growth rate</td></tr> | ||
+ | <tr><td>$$k$$ </td><td> Transcription rate of mRNA under \(P_{tet}\) </td></tr> | ||
+ | <tr><td>$$leak_{P_{lux}}$$ </td><td> Leakage of \(P_{lux}\) </td></tr> | ||
+ | <tr><td>$$leak_{P_{rhl}}$$ </td><td> Leakage of \(P_{rhl}\) </td></tr> | ||
+ | <tr><td>$$\kappa_{Lux}$$ </td><td> Maximum transcription rate of mRNA under \(P_{lux}\)</td></tr> | ||
+ | <tr><td>$$\kappa_{Rhl}$$ </td><td> Maximum transcription rate of mRNA under \(P_{rhl}\) </td></tr> | ||
+ | <tr><td>$$n_{Lux}$$ </td><td> Hill coefficient for \(P_{lux}\)</td></tr> | ||
+ | <tr><td>$$n_{Rhl}$$ </td><td> Hill coefficient for \(P_{rhl}\)</td></tr> | ||
+ | <tr><td>$$K_{mLux}$$ </td><td> Lumped paremeter for the Lux System</td></tr> | ||
+ | <tr><td>$$K_{mRhl}$$ </td><td> Lumped paremeter for the Rhl System</td></tr> | ||
+ | <tr><td>$$F_{DiMazF}$$ </td><td> Cutting rate at ACA sequences on mRNA by MazF dimer</td></tr> | ||
+ | <tr><td>$$f$$ </td><td> The probability of distinction of ACA sequencess in each mRNA</td></tr> | ||
+ | <tr><td>$$f_{mRNA_{RFP}}$$ </td><td> The number of ACA sequences in \(mRNA_{RFP}\)</td></tr> | ||
+ | <tr><td>$$f_{mRNA_{GFP}}$$ </td><td> The number of ACA sequences in \(mRNA_{GFP}\)</td></tr> | ||
+ | <tr><td>$$f_{mRNA_{RhlI}}$$ </td><td> The number of ACA sequences in \(mRNA_{RhlI}\) </td></tr> | ||
+ | <tr><td>$$f_{mRNA_{LasI}}$$ </td><td> The number of ACA sequences in \(mRNA_{LasI}\)</td></tr> | ||
+ | <tr><td>$$f_{mRNA_{MazF}}$$ </td><td> The number of ACA sequences in \(mRNA_{MazF}\) </td></tr> | ||
+ | <tr><td>$$f_{mRNA_{MazE}}$$ </td><td> The number of ACA sequences in \(mRNA_{MazE}\) </td></tr> | ||
+ | <tr><td>$$\alpha$$ </td><td> Translation rate of Protein </td></tr> | ||
+ | <tr><td>$$k_{Di_{MazF}}$$ </td><td> Formation rate of MazF dimer </td></tr> | ||
− | + | <tr><td>$$k_{-Di_{MazF}}$$ </td><td> Dissociation rate of MazF dimer </td></tr> | |
+ | <tr><td>$$k_{Di_{MazE}}$$ </td><td> Formation rate of MazE dimer </td></tr> | ||
+ | <tr><td>$$k_{-Di_{MazE}}$$ </td><td> Dissociation rate of MazE dimer </td></tr> | ||
+ | <tr><td>$$k_{Hexa}$$ </td><td> Formation rate of Maz hexamer </td></tr> | ||
+ | <tr><td>$$k_{-Hexa}$$ </td><td> Dissociation rate of Maz hexamer</td></tr> | ||
+ | <tr><td>$$p_{C4}$$ </td><td> Production rate of C4HSL by RhlI</td></tr> | ||
+ | <tr><td>$$p_{C12}$$ </td><td> Production rate of 3OC12HSL by LuxI </td></tr> | ||
+ | <tr><td>$$D$$ </td><td> Decomposition rate of 3OC12HSL by AmiE </td></tr> | ||
+ | <tr><td>$$d$$ </td><td> Degradation rate of mRNA </td></tr> | ||
+ | <tr><td>$$d_{RFP}$$ </td><td> Degradation rate of RFP</td></tr> | ||
+ | <tr><td>$$d_{GFP}$$ </td><td> Degradation rate of GFP </td></tr> | ||
+ | <tr><td>$$d_{RhlI}$$ </td><td> Degradation rate of RhlI</td></tr> | ||
+ | <tr><td>$$d_{LasI}$$ </td><td> Degradation rate of LasI</td></tr> | ||
+ | <tr><td>$$d_{MazF}$$ </td><td> Degradation rate of MazF</td></tr> | ||
+ | <tr><td>$$d_{DiMazF}$$ </td><td> Degradation rate of MazF dimer</td></tr> | ||
+ | <tr><td>$$d_{MazE}$$ </td><td> Degradation rate of MazE </td></tr> | ||
+ | <tr><td>$$d_{DiMazE}$$ </td><td> Degradation rate of MazE dimer </td></tr> | ||
+ | <tr><td>$$d_{Hexa}$$ </td><td> Degradation rate of Maz Hexamer </td></tr> | ||
+ | <tr><td>$$d_{C4}$$ </td><td> Degradation rate of C4HSL </td></tr> | ||
+ | <tr><td>$$d_{C12}$$ </td><td> Degradation rate of 3OC12HSL </td></tr> | ||
+ | <tr><td>$$d_{AmiE}$$ </td><td> Degradation rate of AmiE </td></tr> | ||
− | + | </tbody> | |
− | + | </table> | |
+ | </div><!-- /modeling_detail_parameter --> | ||
+ | </div><!-- /modeling_detail_wrapper --> | ||
+ | </div><!-- /modeling_detail --> | ||
+ | </div><!-- /story_simalation_contents --> | ||
− | + | <br> | |
− | + | <div id="results"> | |
− | + | <div id="results_header"> | |
− | + | <h3><span>2-2. Results</span></h3> | |
− | + | </div><!-- /_header --> | |
− | + | <p class="normal_text">As a result, we obtained and confirmed the desirable behavior of the whole system by modifying and improving parts. As described below, our simulation showed appropriate transition of fluorescence for the story. | |
− | + | </p> | |
− | + | <div align="center"><img src="https://static.igem.org/mediawiki/2016/7/7f/T--Tokyo_Tech--2-1-1.png" height ="300"><br></div> | |
− | + | <div align="center"><p class="caption" style="font-size: 16px; text-align: center;"><span style="font-weight: bold;">Fig. 5-1.Time-dependent change of the concentrations of fluorescence proteins</span> | |
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | <div align="center"><p class="caption" style="font-size: 16px; text-align: center;"><span style="font-weight: bold;">Fig. | + | |
</p></div><br> | </p></div><br> | ||
− | + | <div align="center"><img src="https://static.igem.org/mediawiki/2016/7/79/T--Tokyo_Tech--4koma.png" width ="800"class="align_right" /><br></div> | |
− | <p class="normal_text"> | + | <div align="center"><p class="caption" style="font-size: 16px; text-align: center;"><span style="font-weight: bold;">Fig. 5-2.The story of "Snow White"</span> |
− | + | <p class="normal_text"> In the blue area of Fig, the fluorescence intensity starts to increase. Snow White coli’s fluorescence intensity exceeds that of the Queen coli’s. It represents Snow White got fairer more and more. <br> | |
− | + | In the pink area of Fig, C12 is being synthesized thanks to the C4’s appearance. And the concentration of C12 increases. As a result, the toxin MazF augments inside of Snow White coli and the Queen coli, suppressing the increment of fluorescent proteins. It is as if the Mirror’s answer transformed the Queen into a Witch, so she can give Snow White a poisoned apple.<br> | |
− | + | In the green area of Fig, C12 increases and the MazF inside Snow White coli induced by it increases even more. So the GFP exceeds the RFP. It looks as if Snow White bit the apple, sinking into unconsciousness promptly.<br> | |
− | + | In the yellow area of Fig, the AmiE synthesized by the introduced Prince coli decomposes the C12 molecules so the amount of MazF inside Snow White diminishes and the amount of C4 increases. It looks as if Prince lifted Snow White and she opened her eyes.<br> | |
− | + | ||
</p> | </p> | ||
− | </div><!-- / | + | </div><!-- /results_contents --> |
− | </div><!-- / | + | </div><!-- /results --> |
+ | </div><!-- /story_simulation_contents --> | ||
+ | </div><!-- /story_smulation --> | ||
+ | |||
+ | <div id="fitting" class="container"> | ||
+ | <div id="fitting_header" class="container_header"> | ||
+ | <h2><span>3. Fitting</span></h2> | ||
+ | </div><!-- /fitting_header --> | ||
+ | <div id="fitting_contents" class="container_contents"> | ||
+ | <p class="normal_text">TEST</p> | ||
+ | </div><!-- /fitting_contents --> | ||
+ | </div><!-- /fitting --> | ||
− | <div id=" | + | <div id="analysis" class="container"> |
− | <div id=" | + | <div id="analysis_header" class="container_header"> |
− | <h2><span> | + | <h2><span>4. Sensitivity Analysis</span></h2> |
− | </div><!-- / | + | </div><!-- /analysis_header --> |
− | <div id=" | + | <div id="analysis_contents" class="container_contents"> |
− | <p class="normal_text"> | + | <p class="normal_text">We performed the sensitivity analysis descried in this section in order to examine which parameter dominates the story.</p> |
− | + | <br> | |
− | + | <div id="requirements"> | |
− | + | <div id="requirements_header"> | |
− | + | <h3><span>4-1. Requirements</span></h3> | |
− | + | </div><!-- /_header --> | |
− | + | <div id="requirements_contents"> | |
− | <div id=" | + | <p class="normal_text">We defined these requirements as the “successful Snow White story.”<br> |
− | <h3><span> | + | 1) At t = 150<br> |
− | </div><!-- /_header --> | + | concentration of RFP > concentration of GFP<br> |
− | <div id=" | + | 2) At t = 700<br> |
− | + | concentration of RFP< concentration of GFP<br> | |
− | + | 3) At t = 1500<br> | |
− | <p class="normal_text"> | + | concentration of RFP > concentration of GFP<br> |
− | + | We the began to analyze the graph that satisfies these requirements that we could say that reproduces 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.<br></p> | |
− | + | <table border="1" style="margin: auto;"> | |
− | + | <tbody> | |
− | + | <tr><td>Parameter </td><td> Range </td><td> Step size </td></tr> | |
− | + | <tr><td>$$D$$ </td><td> $$0.0001<D<0.001$$</td><td>$$0.0001$$</td></tr> | |
− | < | + | <tr><td>$$p_{C4}$$ </td><td> $$0.0001<p_{C4}<1$$</td><td>$$0.001$$</td></tr> |
− | + | <tr><td>$$p_{C12}$$ </td><td> $$0.0001<p_{C12}<1$$</td><td>$$0.001$$</td></tr> | |
− | + | <tr><td>$$α$$ </td><td> $$0.01<α<0.2$$</td><td>$$0.01$$</td></tr> | |
− | + | <tr><td>$$d_{AmiE}$$ </td><td> $$0.001<d_{AmiE}<1$$</td><td>$$0.001$$</td></tr> | |
− | + | </tbody> | |
− | + | </table> | |
− | + | <!--ここにレンジの表が入る--> | |
− | + | ||
− | + | <p class="normal_text">As a result, we obtained the following parameter ranges.</p> | |
− | + | <table border="1" style="margin: auto;"> | |
− | + | <tbody> | |
− | + | <tr><td>Parameter </td><td> Value </td></tr> | |
− | + | <tr><td>$$D$$ </td><td> $$0.0056<D<0.001$$</td></tr> | |
− | + | <tr><td>$$p_{C4}$$ </td><td> $$0.0029<p_{C4}<0.778$$</td></tr> | |
− | + | <tr><td>$$p_{C12}$$ </td><td> $$0.001<p_{C12}<0.217$$</td></tr> | |
− | + | <tr><td>$$α$$ </td><td> $$0.01<α<0.16$$</td></tr> | |
− | + | <tr><td>$$d_{AmiE}$$ </td><td> $$0.001<d_{AmiE}<1$$</td></tr> | |
− | + | </tbody> | |
− | + | </table> | |
− | + | <!--ここにレンジの表がもいっちょ入る--> | |
− | + | </p> | |
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | 4 | + | <br> |
+ | <div id="prince_coli"> | ||
+ | <div id="prince_coli_header"> | ||
+ | <h3><span>4-2. the Prince <span style ="font-style : italic">coli</span> should be put in in the middle</span></h3> | ||
+ | </div><!-- /_header --> | ||
+ | <div id="prince_coli_contents"> | ||
+ | <p class="normal_text">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.</p> | ||
+ | |||
+ | <div align="left"><img src="https://static.igem.org/mediawiki/2016/1/1f/T--Tokyo_Tech--Population-miss.png" width ="500"class="align_left" /><br></div> | ||
+ | <div align="left"><p class="caption" style="font-size: 16px; text-align: center;"><span style="font-weight: bold;">Fig. ??.Number of individuals when the Prince <span style ="font-style : italic">coli</span> is introduced from the beginning</span> | ||
+ | <div align="right"><img src="https://static.igem.org/mediawiki/2016/f/f1/T--Tokyo_Tech--AHL-miss.png" width ="500"class="align_right" /><br></div> | ||
+ | <div align="right"><p class="caption" style="font-size: 16px; text-align: center;"><span style="font-weight: bold;">Fig. ??.AHL concentrations when the Prince <span style ="font-style : italic">coli</span> is introduced from the beginning</span> | ||
+ | |||
+ | <p class="normal_text">As a result, if we introduce the Prince coli from the beginning, the number of Prince coli increases too much (Fig) so the AmiE the Prince coli produces augments and the decomposition of C12 also overly occurs. So C12 is almost inexistent in the medium (Fig).</p> | ||
+ | |||
+ | <div align="left"><img src="https://static.igem.org/mediawiki/2016/e/e7/T--Tokyo_Tech--Population.png" width ="500"class="align_left" /><br></div> | ||
+ | <div align="left"><p class="caption" style="font-size: 16px; text-align: center;"><span style="font-weight: bold;">Fig. ??.Number of individuals when the Prince <span style ="font-style : italic">coli</span> is introduced at t = 700</span> | ||
+ | <div align="right"><img src="https://static.igem.org/mediawiki/2016/4/4f/T--Tokyo_Tech--AHL.png" width ="500"class="align_right" /><br></div> | ||
+ | <div align="right"><p class="caption" style="font-size: 16px; text-align: center;"><span style="font-weight: bold;">Fig. ??.AHL concentrations when the Prince <span style ="font-style : italic">coli</span> is introduced at t = 700</span> | ||
+ | <p class="normal_text">On the other hand, if we introduce the Prince coli at t=700, the number of Prince coli does not increment much (Fig ), so C12 can exist until t=700 and then decreases thanks to the augment of AmiE (Fig ).<br>In conclusion, if we introduce the Prince coli at t=700, the circuit will behave accordingly.</p> | ||
+ | </p> | ||
− | + | <br> | |
− | + | <div id="prhl"> | |
+ | <div id="prhl_header"> | ||
+ | <h3><span>4-3. Prhl should be changed</span></h3> | ||
+ | </div><!-- /_header --> | ||
+ | <div id="prhl_contents"> | ||
+ | <p class="normal_text">In order to confirm whether by using a combination of the current promoters we would be able to reproduce the story with our gene circuit, we performed some simulations based on the results of our assays.</p> | ||
<div align="center"><img src="https://static.igem.org/mediawiki/2016/6/68/T--Tokyo_Tech--projct_model1.png" height ="300"><br></div> | <div align="center"><img src="https://static.igem.org/mediawiki/2016/6/68/T--Tokyo_Tech--projct_model1.png" height ="300"><br></div> | ||
<div align="center"><p class="caption" style="font-size: 16px; text-align: center;"><span style="font-weight: bold;">Fig. </span> The intensity of Plux and Prhl promoters | <div align="center"><p class="caption" style="font-size: 16px; text-align: center;"><span style="font-weight: bold;">Fig. </span> The intensity of Plux and Prhl promoters | ||
− | + | <p class="normal_text">The diagram above shows that the intensity of the two promoters should be in the red region of the figure.The combination of promoters 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 on Prhl so its expression level increases.</p> | |
+ | </p> | ||
− | <p class="normal_text">The | + | <br> |
− | + | <div id="production_ahl"> | |
− | + | <div id="production_ahl_header"> | |
+ | <h3><span>4-4. Production rate of C4HSL and 3OC12HSL by RhlI and LasI</span></h3> | ||
+ | </div><!-- /_header --> | ||
+ | <div id="production_ahl_contents"> | ||
+ | <p class="normal_text">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.</p> | ||
+ | |||
+ | <div align="center"><img src="https://static.igem.org/mediawiki/2016/9/9f/T--Tokyo_Tech--Production-AHL.png" height ="300"><br></div> | ||
+ | <div align="center"><p class="caption" style="font-size: 16px; text-align: center;"><span style="font-weight: bold;">Fig. </span> Concentrations of GFP and RFP dependencies of production rate of C4 by RhlI | ||
+ | |||
+ | <p class="normal_text">Each line corresponds to the transition of the concentration of RFP/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.<br> | ||
+ | if the production rate of C4 is between 0.029 and 0.778, our system can reproduce the story.<br> | ||
+ | If this parameter is too small, The production of LasI by the Queen <span style ="font-style : italic">coli</span> 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 <span style ="font-style : italic">coli</span> 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.</p> | ||
+ | |||
+ | <div align="center"><img src="https://static.igem.org/mediawiki/2016/0/05/T--Tokyo_Tech--Production-AHL2.png" height ="300"><br></div> | ||
+ | <div align="center"><p class="caption" style="font-size: 16px; text-align: center;"><span style="font-weight: bold;">Fig. </span> Concentrations of GFP and RFP dependencies of production rate of C12 by LasI | ||
+ | |||
+ | <p class="normal_text">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.<br> | ||
+ | If production rate of C12 is between 0.001 and 0.217, our system can reproduce the story.</p> | ||
+ | </p> | ||
− | + | <br> | |
− | + | <div id="translation"> | |
− | + | <div id="translation_header"> | |
− | + | <h3><span>4-5. Translation rate of protein</span></h3> | |
− | + | </div><!-- /_header --> | |
+ | <div id="translation_contents"> | ||
+ | <p class="normal_text">Translation rate affects the production of protein.</p> | ||
+ | |||
+ | <div align="center"><img src="https://static.igem.org/mediawiki/2016/0/03/T-Tokyo_Tech--Translation.png" height ="300"><br></div> | ||
+ | <div align="center"><p class="caption" style="font-size: 16px; text-align: center;"><span style="font-weight: bold;">Fig. </span> Concentrations of GFP and RFP dependencies of translation rate of proteins | ||
+ | |||
+ | <p class="normal_text">Each line corresponds to the transition of the concentration of RFP/GFP with a certain translation rate. Red and green lines correspond to RFP and GFP, respectively. Blighter one indicates higher translation rate.<br> | ||
+ | If translation rate of protein is between 0.01 and 0.16, our system can reproduce the story. | ||
+ | </p> | ||
− | <div align="center"><img src="https://static.igem.org/mediawiki/2016/ | + | <br> |
− | <div align="center"><p class="caption" style="font-size: 16px; text-align: center;"><span style="font-weight: bold;">Fig. </span> | + | <div id="decomposition"> |
− | < | + | <div id="decomposition_header"> |
+ | <h3><span>4-6. Decomposition rate of C12 by AmiE</span></h3> | ||
+ | </div><!-- /_header --> | ||
+ | <div id="decomposition_contents"> | ||
+ | <p class="normal_text">AmiE decomposes C12. The concentration of C12 after input of the Prince <span style ="font-style : italic">coli</span> is changed by AmiE. If the decomposition rate of C12 is too small, C12 does not decrease enough so the MazF inside Snow White <span style ="font-style : italic">coli</span> continues being expressed and the concentration of RFP decreases.</p> | ||
+ | |||
+ | <div align="center"><img src="https://static.igem.org/mediawiki/2016/0/08/T--Tokyo_Tech--Decomposition.png" height ="300"><br></div> | ||
+ | <div align="center"><p class="caption" style="font-size: 16px; text-align: center;"><span style="font-weight: bold;">Fig. </span> Concentrations of GFP and RFP dependencies of decomposition rate of C12 by AmiE | ||
+ | |||
+ | <p class="normal_text">Each line corresponds to the transition of the concentration of RFP/GFP with a certain value of D. Red and green lines correspond to RFP and GFP, respectively. Blighter one indicates higher value of decomposition rate of C12.<br> | ||
+ | If degradation rate of C12 is higher than 0.00056, our system can reproduce the story.</p> | ||
+ | </p> | ||
− | + | <br> | |
+ | <div id="degradation _amie"> | ||
+ | <div id="degradation_amie_header"> | ||
+ | <h3><span>4-7. Degradation rate of AmiE</span></h3> | ||
+ | </div><!-- /_header --> | ||
+ | <div id="degradation_amie_contents"> | ||
+ | |||
+ | <div align="center"><img src="https://static.igem.org/mediawiki/2016/d/d4/T--Tokyo_Tech--DegradationAmiE.png" height ="300"><br></div> | ||
+ | <div align="center"><p class="caption" style="font-size: 16px; text-align: center;"><span style="font-weight: bold;">Fig. </span> Concentration of GFP and RFP dependencies of degradation rate of AmiE | ||
+ | |||
+ | <p class="normal_text">Each line corresponds to the transition of the concentration of RFP/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.<br> | ||
+ | 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.</p> | ||
+ | </p> | ||
+ | <br> | ||
+ | <div id="degradation_protein"> | ||
+ | <div id="degradation_protein"> | ||
+ | <h3><span>4-8. Degradation rate of RFP and GFP</span></h3> | ||
+ | </div><!-- /_header --> | ||
+ | <div id="degradation_protein_contents"> | ||
+ | <p class="normal_text">The degradation rate of RFP and GFP is key to the success the story of ‘Snow White’.<br> | ||
+ | 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 reproduced.</p> | ||
+ | |||
+ | <div align="center"><img src="https://static.igem.org/mediawiki/2016/a/a0/T--Tokyo_Tech--DegradationGFP.png" height ="300"><br></div> | ||
+ | <div align="center"><p class="caption" style="font-size: 16px; text-align: center;"><span style="font-weight: bold;">Fig. </span> Relation of degradation rate of GFP and RFP | ||
+ | |||
+ | <p class="normal_text">The story is reproduced only if the degradation rate of RFP and GFP are the same.</p> | ||
+ | </p> | ||
+ | </div><!-- /analysis_contents --> | ||
+ | </div><!-- /analysis --> | ||
+ | <br> | ||
+ | <div id="software" class="container"> | ||
+ | <div id="software_header" class="container_header"> | ||
+ | <h2><span>5. Software</span></h2> | ||
+ | </div><!-- /software_header --> | ||
+ | <div id="software_contents" class="container_contents"> | ||
+ | <br> | ||
+ | <div id="abstract"> | ||
+ | <div id="abstract_header"> | ||
+ | <h3><span>5-1. Abstract</span></h3> | ||
+ | </div><!-- /_header --> | ||
+ | <div id="abstract_contents"> | ||
+ | <img src="https://static.igem.org/mediawiki/2016/8/8b/T--Tokyo_Tech--ACA_Dwarfs.jpg" height ="300"><br></div> | ||
+ | <p class="normal_text">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. <br>Therewith we improve the practicality of the characteristic of the <span style ="font-style : italic"> mazEF </span> system. For example you can let protein A express constantly by eliminating ACA sequences of the sequence, while letting protein B suspend to express, at the desirable timing, by expression of MazF.<br>This software also evades using rare codons, so you don’t have to worry about them.</p> | ||
+ | <br> | ||
+ | <div id="key_achievement"> | ||
+ | <div id="key_achievement_header"> | ||
+ | <h3><span>5-2. Key achievement</span></h3> | ||
+ | </div><!-- /_header --> | ||
+ | <p class="normal_text"> | ||
+ | ・Provided the tool regulating the number of ACA sequences<br> | ||
+ | ・Released under open-source license so everyone can use it<br> | ||
+ | ・Able to correspond to any base arrangements<br> | ||
+ | ・Rare codons are evaded<br> | ||
+ | ・Extend the application field of <span style ="font-style : italic">mazE</span> system<br> | ||
+ | </p> | ||
− | + | <br> | |
− | < | + | <div id="work_flow"> |
− | </div> | + | <div id="work_flow_header"> |
− | + | <h3><span>5-3. Work flow</span></h3> | |
− | <div id=" | + | </div><!-- /_header --> |
− | <h3><span>4. | + | <p class="normal_text">TEST |
+ | </p> | ||
+ | <br> | ||
+ | <div id="demo"> | ||
+ | <div id="demo_header"> | ||
+ | <h3><span>5-4. Demonstration</span></h3> | ||
+ | </div><!-- /_header --> | ||
+ | <p class="normal_text">TEST | ||
+ | </p> | ||
+ | <br> | ||
+ | <div id="download"> | ||
+ | <div id="download_header"> | ||
+ | <h3><span>5-5. Download</span></h3> | ||
</div><!-- /_header --> | </div><!-- /_header --> | ||
− | + | <p class="normal_text">To download click <a href="https://static.igem.org/mediawiki/2016/3/32/T--Tokyo_Tech--ACA_Dwarfs.zip">here</a>.</p> | |
− | + | <p class="normal_text">The code is available on <a href="https://github.com/Ryuta339/ACADwarfs">github</a>.</p> | |
− | + | </div><!-- /software_contents --> | |
− | + | </div><!-- /software --> | |
− | + | </div><!-- /main_contents --> | |
− | + | <script type="text/javascript"> | |
− | + | function show(idName){ | |
− | + | var elem = document.getElementById(idName); | |
− | + | if(elem.className === 'on'){ | |
− | + | elem.className = 'off'; | |
− | + | var elem2 = document.getElementsByClassName('showHidden'); | |
− | + | elem2[0].textContent = 'Read More'; | |
− | + | } else { | |
− | + | elem.className = 'on'; | |
− | + | var elem2 = document.getElementsByClassName('showHidden'); | |
− | + | elem2[0].textContent = 'Hide'; | |
− | + | } | |
− | + | } | |
− | + | </script> | |
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
</body> | </body> | ||
</html> | </html> |
Revision as of 06:15, 18 October 2016
![](https://static.igem.org/mediawiki/2016/f/f6/T--Tokyo_Tech--goTop.png)
Model
Contents
1. Overview
To reproduce the story of ”Snow White”, we have designed the cell-cell communication system with improved or characterized parts and collected data from comprehensive experiments. Furthermore, we constructed the mathematical model in order 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 and confirm the feasibility of our system. In addition, in order to utilize TA (Toxin-Antitoxin) 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
To simulate our gene circuits, we developed an ordinary differential equation model.
[Detailed Description for Modeling]
Differencial Equations
Snow White
\begin{equation} \frac{d[mRNA_{RFP}]}{dt} = k - d[mRNA_{RFP}] - F_{DiMazF}(1-(1-f)^{f_{mRNA_{RFP}}})[mRNA_{RFP}][DiMazF] \end{equation} \begin{equation} \frac{d[mRNA_{RhlI}]}{dt} = leak_{P_{lux}} + \frac{\kappa_{Lux}[C12]^{n_{Lux}}}{K_{mLux}^{n_{Lux}} + [C12]^{n_{Lux}}} - d[mRNA_{RhlI}] - F_{DiMazF}f[mRNA_{RhlI}][DiMazF] \end{equation} \begin{equation} \frac{d[RFP]}{dt} = \alpha [mRNA_{RFP}] - d_{RFP}[RFP] \end{equation} \begin{equation} \frac{d[RhlI]}{dt} = \alpha [mRNA_{RhlI}] - d_{RhlI}[RhlI] \end{equation} \begin{equation} \frac{d[C4]}{dt} = p_{C4}[RhlI]P_{Snow White} - d_{C4}[C4] \end{equation} \begin{equation} \frac{d[mRNA_{MazF}]}{dt} = leak_{P_{lux}} + \frac{\kappa_{Lux}[C12]^{n_{Lux}}}{K_{mLux}^{n_{Lux}}+ [C12]^{n_{Lux}}} \\ - d[mRNA_{MazF}] - F_{DiMazF}(1-(1-f)^{f_{mRNA_{MazF}}})[mRNA_{MazF}][DiMazF] \end{equation} \begin{equation} \frac{d[mRNA_{MazE}]}{dt} = k - d[mRNA_{MazE}] - F_{DiMazF}(1-(1-f)^{f_{mRNA_{MazE}}})[mRNA_{MazE}][DiMazF] \end{equation} \begin{equation} \frac{d[MazF]}{dt} = \alpha [mRNA_{MazF}] - 2k_{Di_{MazF}}[MazF] + 2k_{-Di_{MazF}}[DiMazF] - d_{MazF}[MazF] \end{equation} \begin{equation} \frac{d[DiMazF]}{dt} = k_{Di_{MazF}}[MazF] - k_{-Di_{MazF}}[DiMazF] - 2k_{Hexa}[DiMazE][DiMazF]^2 \\ + 2k_{-Hexa}[MazHexamer] - d_{DiMazF}[DiMazF] \end{equation} \begin{equation} \frac{d[MazE]}{dt} = \alpha [mRNA_{MazE}] - 2k_{Di_{MazE}}[MazE] + 2k_{-Di_{MazE}}[DiMazE] - d_{MazE}[MazE] \end{equation} \begin{equation} \frac{d[DiMazE]}{dt} = k_{Di_{MazE}}[MazE] - k_{-Di_{MazE}}[DiMazE] - k_{Hexa}[DiMazE][DiMazF]^2 \\ + k_{-Hexa}[MazHexamer] - d_{DiMazE}[DiMazE] \end{equation} \begin{equation} \frac{d[MazHexa]}{dt} = k_{Hexa}[DiMazE][DiMazF]^2 - k_{-Hexa}[MazHexa] - d_{Hexa}[MazHexa] \end{equation} \begin{equation} \frac{dP_{Snow White}}{dt} = g \frac{E_{DiMazF}}{E_{DiMazF}+[DiMazF]}\left(1- \frac{P_{Snow White}+P_{Queen}+P_{Prince}}{P_{max}} \right) P_{Snow White} \end{equation}Queen
\begin{equation} \frac{d[mRNA_{GFP}]}{dt} = k - d[mRNA_{GFP}] - F_{DiMazF}(1-(1-f)^{f_{mRNA_{GFP}}})[mRNA_{GFP}][DiMazF] \end{equation} \begin{equation} \frac{d[mRNA_{LasI}]}{dt} = leak_{P_{rhl}} + \frac{\kappa_{Rhl}[C4]^{n_{Rhl}}}{K_{mRhl}^{n_{Rhl}} + [C4]^{n_{Rhl}}} \\ - d[mRNA_{LasI}] - F_{DiMazF}(1-(1-f)^{f_{mRNA_{LasI}}})[mRNA_{LasI}][DiMazF] \end{equation} \begin{equation} \frac{d[GFP]}{dt} = \alpha [mRNA_{GFP}] - d_{GFP}[GFP] \end{equation} \begin{equation} \frac{d[LasI]}{dt} = \alpha [mRNA_{LasI}] - d_{LasI}[LasI] \end{equation} \begin{equation} \frac{d[C12]}{dt} = p_{C12}[LasI]P_{Queen} - d_{C12}[C12] - D[C12][AmiE] \end{equation} \begin{equation} \frac{d[mRNA_{MazF}]}{dt} = leak_{P_{lux}} + \frac{\kappa_{Rhl}[C4]^{n_{Rhl}}}{K_{mRhl}^{n_{Rhl}} + [C4]^{n_{Rhl}}} \\ - d[mRNA_{MazF}] - F_{DiMazF}(1-(1-f)^{f_{mRNA_{MazF}}})[mRNA_{MazF}][DiMazF] \end{equation} \begin{equation} \frac{d[mRNA_{MazE}]}{dt} = k - d[mRNA_{MazE}] - F_{DiMazF}(1-(1-f)^{f_{mRNA_{MazE}}})[mRNA_{MazE}][DiMazF] \end{equation} \begin{equation} \frac{d[MazF]}{dt} = \alpha [mRNA_{MazF}] - 2k_{Di_{MazF}}[MazF] + 2k_{-Di_{MazF}}[DiMazF] - d_{MazF}[MazF] \end{equation} \begin{equation} \frac{d[DiMazF]}{dt} = k_{Di_{MazF}}[MazF] - k_{-Di_{MazF}}[DiMazF] - 2k_{Hexa}[DiMazE][DiMazF]^2 \\ + 2k_{-Hexa}[MazHexamer] - d_{DiMazF}[DiMazF] \end{equation} \begin{equation} \frac{d[MazE]}{dt} = \alpha [mRNA_{MazE}] - 2k_{Di_{MazE}}[MazE] + 2k_{-Di_{MazE}}[DiMazE] - d_{MazE}[MazE] \end{equation} \begin{equation} \frac{d[DiMazE]}{dt} = k_{Di_{MazE}}[MazE] - k_{-Di_{MazE}}[DiMazE] - k_{Hexa}[DiMazE][DiMazF]^2 \\ + k_{-Hexa}[MazHexamer] - d_{DiMazE}[DiMazE] \end{equation} \begin{equation} \frac{d[MazHexa]}{dt} = k_{Hexa}[DiMazE][DiMazF]^2 - k_{-Hexa}[MazHexa] - d_{Hexa}[MazHexa] \end{equation} \begin{equation} \frac{dP_{Queen}}{dt} = g \frac{E_{DiMazF}}{E_{DiMazF}+[DiMazF]}\left(1- \frac{P_{Snow White}+P_{Queen}+P_{Prince}}{P_{max}}\right) P_{Queen}\\ \end{equation}Prince
\begin{equation} \frac{d[mRNA_{AmiE}]}{dt} = leak_{P_{lux}} + \frac{\kappa_{Lux}[C12]^n}{K_{mLux}^n + [C12]^n} - d[mRNA_{AmiE}] \end{equation} \begin{equation} \frac{d[AmiE]}{dt} = \alpha [mRNA_{AmiE}]P_{Prince} - d_{AmiE}[AmiE] \end{equation} \begin{equation} \frac{dP_{Prince}}{dt} = g\left(1- \frac{P_{Snow White}+P_{Queen}+P_{Prince}}{P_{max}}\right) P_{Prince} \end{equation}Explanation about Parameters
Parameter | Description |
$$g$$ | Growth rate of each cells |
$$P_{max}$$ | Carrying capacity |
$$E_{DiMazF}$$ | Effect of MazF dimer on growth rate |
$$k$$ | Transcription rate of mRNA under \(P_{tet}\) |
$$leak_{P_{lux}}$$ | Leakage of \(P_{lux}\) |
$$leak_{P_{rhl}}$$ | Leakage of \(P_{rhl}\) |
$$\kappa_{Lux}$$ | Maximum transcription rate of mRNA under \(P_{lux}\) |
$$\kappa_{Rhl}$$ | Maximum transcription rate of mRNA under \(P_{rhl}\) |
$$n_{Lux}$$ | Hill coefficient for \(P_{lux}\) |
$$n_{Rhl}$$ | Hill coefficient for \(P_{rhl}\) |
$$K_{mLux}$$ | Lumped paremeter for the Lux System |
$$K_{mRhl}$$ | Lumped paremeter for the Rhl System |
$$F_{DiMazF}$$ | Cutting rate at ACA sequences on mRNA by MazF dimer |
$$f$$ | The probability of distinction of ACA sequencess in each mRNA |
$$f_{mRNA_{RFP}}$$ | The number of ACA sequences in \(mRNA_{RFP}\) |
$$f_{mRNA_{GFP}}$$ | The number of ACA sequences in \(mRNA_{GFP}\) |
$$f_{mRNA_{RhlI}}$$ | The number of ACA sequences in \(mRNA_{RhlI}\) |
$$f_{mRNA_{LasI}}$$ | The number of ACA sequences in \(mRNA_{LasI}\) |
$$f_{mRNA_{MazF}}$$ | The number of ACA sequences in \(mRNA_{MazF}\) |
$$f_{mRNA_{MazE}}$$ | The number of ACA sequences in \(mRNA_{MazE}\) |
$$\alpha$$ | Translation rate of Protein |
$$k_{Di_{MazF}}$$ | Formation rate of MazF dimer |
$$k_{-Di_{MazF}}$$ | Dissociation rate of MazF dimer |
$$k_{Di_{MazE}}$$ | Formation rate of MazE dimer |
$$k_{-Di_{MazE}}$$ | Dissociation rate of MazE dimer |
$$k_{Hexa}$$ | Formation rate of Maz hexamer |
$$k_{-Hexa}$$ | Dissociation rate of Maz hexamer |
$$p_{C4}$$ | Production rate of C4HSL by RhlI |
$$p_{C12}$$ | Production rate of 3OC12HSL by LuxI |
$$D$$ | Decomposition rate of 3OC12HSL by AmiE |
$$d$$ | Degradation rate of mRNA |
$$d_{RFP}$$ | Degradation rate of RFP |
$$d_{GFP}$$ | Degradation rate of GFP |
$$d_{RhlI}$$ | Degradation rate of RhlI |
$$d_{LasI}$$ | Degradation rate of LasI |
$$d_{MazF}$$ | Degradation rate of MazF |
$$d_{DiMazF}$$ | Degradation rate of MazF dimer |
$$d_{MazE}$$ | Degradation rate of MazE |
$$d_{DiMazE}$$ | Degradation rate of MazE dimer |
$$d_{Hexa}$$ | Degradation rate of Maz Hexamer |
$$d_{C4}$$ | Degradation rate of C4HSL |
$$d_{C12}$$ | Degradation rate of 3OC12HSL |
$$d_{AmiE}$$ | Degradation rate of AmiE |
2-2. Results
As a result, we obtained and confirmed the desirable behavior of the whole system by modifying and improving parts. As described below, our simulation showed appropriate transition of fluorescence for the story.
![](https://static.igem.org/mediawiki/2016/7/7f/T--Tokyo_Tech--2-1-1.png)
Fig. 5-1.Time-dependent change of the concentrations of fluorescence proteins
![](https://static.igem.org/mediawiki/2016/7/79/T--Tokyo_Tech--4koma.png)
Fig. 5-2.The story of "Snow White"
In the blue area of Fig, the fluorescence intensity starts to increase. Snow White coli’s fluorescence intensity exceeds that of the Queen coli’s. It represents Snow White got fairer more and more.
In the pink area of Fig, C12 is being synthesized thanks to the C4’s appearance. And the concentration of C12 increases. As a result, the toxin MazF augments inside of Snow White coli and the Queen coli, suppressing the increment of fluorescent proteins. It is as if the Mirror’s answer transformed the Queen into a Witch, so she can give Snow White a poisoned apple.
In the green area of Fig, C12 increases and the MazF inside Snow White coli induced by it increases even more. So the GFP exceeds the RFP. It looks as if Snow White bit the apple, sinking into unconsciousness promptly.
In the yellow area of Fig, the AmiE synthesized by the introduced Prince coli decomposes the C12 molecules so the amount of MazF inside Snow White diminishes and the amount of C4 increases. It looks as if Prince lifted Snow White and she opened her eyes.
3. Fitting
TEST
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 reproduces 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$$0.0001$$ | |
$$p_{C4}$$ | $$0.0001$$0.001$$ | |
$$p_{C12}$$ | $$0.0001$$0.001$$ | |
$$α$$ | $$0.01<α<0.2$$ | $$0.01$$ |
$$d_{AmiE}$$ | $$0.001$$0.001$$ | |
As a result, we obtained the following parameter ranges.
Parameter | Value |
$$D$$ | $$0.0056 |
$$p_{C4}$$ | $$0.0029 |
$$p_{C12}$$ | $$0.001 |
$$α$$ | $$0.01<α<0.16$$ |
$$d_{AmiE}$$ | $$0.001 |
4-2. the Prince coli should be put in in the middle
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.
![](https://static.igem.org/mediawiki/2016/1/1f/T--Tokyo_Tech--Population-miss.png)
Fig. ??.Number of individuals when the Prince coli is introduced from the beginning
![](https://static.igem.org/mediawiki/2016/f/f1/T--Tokyo_Tech--AHL-miss.png)
Fig. ??.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) so the AmiE the Prince coli produces augments and the decomposition of C12 also overly occurs. So C12 is almost inexistent in the medium (Fig).
![](https://static.igem.org/mediawiki/2016/e/e7/T--Tokyo_Tech--Population.png)
Fig. ??.Number of individuals when the Prince coli is introduced at t = 700
![](https://static.igem.org/mediawiki/2016/4/4f/T--Tokyo_Tech--AHL.png)
Fig. ??.AHL concentrations when the Prince coli is introduced at t = 700
On the other hand, if we introduce the Prince coli at t=700, the number of Prince coli does not increment much (Fig ), so C12 can exist until t=700 and then decreases thanks to the augment of AmiE (Fig ).
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 whether by using a combination of the current promoters we would be able to reproduce the story with our gene circuit, we performed some simulations based on the results of our assays.
![](https://static.igem.org/mediawiki/2016/6/68/T--Tokyo_Tech--projct_model1.png)
Fig. 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 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 on Prhl so its expression level increases.
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.
![](https://static.igem.org/mediawiki/2016/9/9f/T--Tokyo_Tech--Production-AHL.png)
Fig. Concentrations of GFP and RFP dependencies of production rate of C4 by RhlI
Each line corresponds to the transition of the concentration of RFP/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 reproduce 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.
![](https://static.igem.org/mediawiki/2016/0/05/T--Tokyo_Tech--Production-AHL2.png)
Fig. 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 reproduce the story.
4-5. Translation rate of protein
Translation rate affects the production of protein.
![](https://static.igem.org/mediawiki/2016/0/03/T-Tokyo_Tech--Translation.png)
Fig. Concentrations of GFP and RFP dependencies of translation rate of proteins
Each line corresponds to the transition of the concentration of RFP/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 reproduce 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.
![](https://static.igem.org/mediawiki/2016/0/08/T--Tokyo_Tech--Decomposition.png)
Fig. Concentrations of GFP and RFP dependencies of decomposition rate of C12 by AmiE
Each line corresponds to the transition of the concentration of RFP/GFP with a certain value of D. 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 reproduce the story.
4-7. Degradation rate of AmiE
![](https://static.igem.org/mediawiki/2016/d/d4/T--Tokyo_Tech--DegradationAmiE.png)
Fig. Concentration of GFP and RFP dependencies of degradation rate of AmiE
Each line corresponds to the transition of the concentration of RFP/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 reproduced.
![](https://static.igem.org/mediawiki/2016/a/a0/T--Tokyo_Tech--DegradationGFP.png)
Fig. Relation of degradation rate of GFP and RFP
The story is reproduced only if the degradation rate of RFP and GFP are the same.
5. Software
5-1. Abstract
![](https://static.igem.org/mediawiki/2016/8/8b/T--Tokyo_Tech--ACA_Dwarfs.jpg)
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.
Therewith 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 suspend to express, at the desirable timing, by expression of MazF.
This software also evades using rare codons, so you don’t have to worry about them.
5-2. Key achievement
・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 mazE system