Difference between revisions of "Team:ZJUT-China/Model"

 
(10 intermediate revisions by 2 users not shown)
Line 127: Line 127:
 
#bodyContent h1, #bodyContent h2, #bodyContent h3, #bodyContent h4, #bodyContent h5 { margin-bottom: 0px; }
 
#bodyContent h1, #bodyContent h2, #bodyContent h3, #bodyContent h4, #bodyContent h5 { margin-bottom: 0px; }
  
.judges-will-not-evaluate { border: 4px solid #e4dede; padding: 2% !important; width: 92%!important;}
 
  
  
Line 569: Line 568:
  
 
}
 
}
 +
p{
 +
  text-indent:2em!important;
 +
}
 +
 
.homepage_img{
 
.homepage_img{
 
   text-align: center;
 
   text-align: center;
Line 839: Line 842:
 
}
 
}
 
/* 底部 */
 
/* 底部 */
.footer{
+
  .footer{
  background: #eee;
+
    background: #eee;
  margin-top: 30px;
+
    margin-top: 30px;
  padding: 30px 0;
+
    padding: 30px 0;
  width: 100%;
+
    width: 100%;
}
+
  }
.foottop{
+
  .foottop{
  height: 230px;
+
    height: 330px;
  width: 100%;
+
    width: 100%;
}
+
  }
.footbottom{
+
  .footbottom{
  height:50px;
+
    height:50px;
  width: 100%;
+
    width: 100%;
}
+
  }
.footbottom .note{
+
  .footbottom .note{
  padding-top: 40px;
+
    padding-top: 40px;
}
+
  }
.footmessage-left,
+
  .footmessage-left,
.footmessage-right,
+
  .footmessage-left2,
.footmessage-middle {
+
  .footmessage-right,
  width: 30%;
+
  .footmessage-middle {
  float: left;
+
    width: 24%;
  height: 100%;
+
    float: left;
}
+
    height: 100%;
.footmessage-right, img{
+
  }
   max-width: 100%;
+
  .footmessage-middle img{
  max-height: 100%;
+
    width: 60%;
  padding-top: 10px;
+
    padding-top: 60px;
   padding-left: 30px;
+
   }
}
+
  .footmessage-right img{
.footmessage-left img{
+
    width: 80%;
  padding-top: 22px;
+
    max-height: 100%;
  padding-left: 30px;
+
    padding-top: 10px;
  max-width: 100%;
+
    padding-left: 30px;
  max-height: 100%;
+
   }
}
+
  .footmessage-left2 img{
.footmessage-right h3{
+
    width: 80%;
  font-size: 30px;
+
    padding-left: 30px;
}
+
    padding-top: 15px;
.footmessage-right p{
+
  }
   font-size: 15px!important;
+
  .footmessage-left img{
}
+
    padding-top: 60px;
 +
    padding-left: 30px;
 +
    width: 80%;
 +
    max-height: 100%;
 +
  }
 +
  .footmessage-right h3{
 +
    font-size: 30px;
 +
  }
 +
  .footmessage-right{
 +
    padding-left: 30px;
 +
   }
 +
  .footmessage-right p{
 +
    font-size: 15px!important;
 +
  }
 
.col-10 table{
 
.col-10 table{
 
   align:center!important;
 
   align:center!important;
 
}
 
}
 
.modelingimg{
 
.modelingimg{
   width: 80%;
+
   width: 100%;
   padding-left: 10%;
+
  text-align: center;
 +
   padding-top: 20px;
 +
  padding-bottom: 20px;
 +
}
 +
.modelingimg img{
 +
  max-width: 60%;
 
}
 
}
 
.note2{
 
.note2{
 
   color:#708090;
 
   color:#708090;
 
   font-size: 18px;
 
   font-size: 18px;
 +
  width: 80%;
 +
  padding-left: 10%;
 +
}
 +
.note3{
 +
    color:#708090;
 +
    font-size: 18px;
 +
    width: 60%;
 +
    padding-left: 20%;
 +
    line-height: 30px;
 +
 
}
 
}
 
.modelingmath{
 
.modelingmath{
Line 898: Line 929:
 
   margin-bottom: 15px;
 
   margin-bottom: 15px;
 
}
 
}
 +
.footmessage-middle img{
 +
      padding-left: 60px;
 +
        }
 
</style>
 
</style>
  
Line 968: Line 1,002:
 
             <div class="dropdown-content">
 
             <div class="dropdown-content">
 
               <a href="https://2018.igem.org/Team:ZJUT-China/Human_Practices">Integrated hp</a>
 
               <a href="https://2018.igem.org/Team:ZJUT-China/Human_Practices">Integrated hp</a>
               <a href="https://2018.igem.org/Team:ZJUT-China/Public_Engagement">Public Engagement</a>
+
               <a href="https://2018.igem.org/Team:ZJUT-China/Public_Engagement">Public Engagement & Education</a>
              <a href="https://2018.igem.org/Team:ZJUT-China/Education">Education</a>
+
 
             </div>
 
             </div>
 
           </div>
 
           </div>
Line 995: Line 1,028:
 
     <h1>Abstract</h1>
 
     <h1>Abstract</h1>
 
     <p>
 
     <p>
       In order to make sure that our system could be common used by the researchers who want to clean the arg after experiments, we developed a model to Calculate the point when all arg are cleaned and the bacteria suicide. We got test data of our parts and then built rate equations. And next we used those rate equations to simulate how our system work and tried to find a better ratio of the vector’s copy number.
+
       In order to improve the feasibility of our system utilized by other researchers , we have tried our best to develop a model to predict the time point when the all antibiotic resistance genes(ARGs) are cut off and the cell auto-lysis is going to begin. This model is made up of some rate equations which are based on experimental data of our system. Then we are able to simulate the operation of our system and gain the information we need from the analog results.
 
     </p>
 
     </p>
     <h1>System modeling Var1.0</h1>
+
     <h1>System modeling Ver1.0</h1>
     <div class="modelingmath">
+
    For the purpose of calculating the time point when the all ARGs are cut off and the cell auto-lysis is going to begin,we have tried to construct a model used by the following equations.
            d[express-cas]/dt=[vector]*vm/(1+[light_signal])
+
     <div class="modelingimg">
       <br>  d[cut off]/dt=kf*[sg]*[cas9]*[arg]/(km+[arg])
+
       <img src="https://static.igem.org/mediawiki/2018/9/9a/T--ZJUT-China--modeling10181.png" alt="">
    </div>
+
    <div class="note2">
+
        Our ODE model of the light-controlled arg cutting system ver1.0 . We assumed that the expression rate is proportional to the copy number of the vectors, and vm means the mixture rate of one vector express the cas9, [light signal] we use is a dimensionless parameter here, and the cut-off rate of cas9 conform the Michaelis-Menten equation.
+
 
     </div>
 
     </div>
 +
    <div class="note3">Our ODE model of the light control ARG cutting system ver1.0.
 +
        <br>[Cas9]: The relative number of Cas9 protein
 +
        <br>[vector]: The relative copy number of the vector
 +
        <br>[light signal]: The relative light strength
 +
      <br>  [ARG]: The relative number of ARGs
 +
        <br>[sg]: The relative number of sgRNA
 +
      </div>
 
     <h1>Result</h1>
 
     <h1>Result</h1>
 
     <p>
 
     <p>
       This is the result we got by the ODE model of the light-controlled arg cutting system ver1.0
+
       The result that we gain from the system modeling ver1.0 is shown below.
 
     </p>
 
     </p>
 
     <div class="modelingimg">
 
     <div class="modelingimg">
       <img src="https://static.igem.org/mediawiki/2018/6/6c/T--ZJUT-China--modeling1.png" alt="">
+
       <img src="https://static.igem.org/mediawiki/2018/3/31/T--ZJUT-China--modeling10182.png" alt="">
 +
    </div>
 +
    <h1>System modeling Ver1.1</h1>
 +
      <p>
 +
        Our system modeling ver1.0 is limited by problems that we are unable to verify the validity of the model through the experiments.Because the quantization of cas9 protein and ARGs is difficult.Therefore,we have tried to build up another model which is named as System modeling Ver1.1. This model can be used to predict the time point which is same as the System modeling Ver1.0 by growth curve. We are able to verify our model by experiments because the growth curve is easily gained by experiments.In addition,we can also prove the validity of system model ver1.0 for the reason that the equations of this model is based on the system model ver1.0.
 +
      </p>
 +
      <div class="modelingimg">
 +
        <img src="https://static.igem.org/mediawiki/2018/0/0a/T--ZJUT-China--modeling10183.png" alt="">
 +
      </div>
 +
      <div class="note3">
 +
        Our ODE model of the light control ARG cutting system ver1.1.(the experimental group with the system and ARG)
 +
        <br>[Cas9]: The relative number of Cas9 protein
 +
        <br>[vector]: The relative copy number of the vector
 +
        <br>[ara]: The concentration of arabinose used to induced expression
 +
        <br>[ARG]: The relative number of ARGs
 +
        <br>[sg]: The relative number of sgRNA
 +
        <br>[cell]: The OD600 of viable cell
 +
      </div>
 +
      <div class="modelingimg">
 +
        <img src="https://static.igem.org/mediawiki/2018/8/8c/T--ZJUT-China--modeling10184.png" alt="">
 +
      </div>
 +
        <div class="note3">
 +
          Our ODE model of the light control ARG cutting system ver1.1.(the control group with ARG but without the system)
 +
          <br>[cell-c]: The OD600 of viable cell
 +
        </div>
 +
      <h1>Result</h1>
 +
    <div class="modelingimg">
 +
      <img src="https://static.igem.org/mediawiki/2018/b/bd/T--ZJUT-China--modeli20.png" alt="">
 +
    </div>
 +
    <div class="note">
 +
      When time=5h we put in different concentrations of Arabinose in the medium.
 +
      <br>The points are the data get tested by wet lab and the lines are the data simulated by computer.
 +
    </div>
 +
    <p>Our ODE model after fit data:</p>
 +
    <div class="modelingimg">
 +
      <img src="https://static.igem.org/mediawiki/2018/c/c8/T--ZJUT-China--modeling10186.png" alt="">
 +
    </div>
 +
    <div class="note3">
 +
      Our ODE model of the light control ARG cutting system ver1.1.(the experimental group)
 +
    <br>  [Cas9]: The relative number of Cas9 protein
 +
    <br>  [vector]: The relative copy number of the vector
 +
    <br>  [ara]: The concentration of arabinose used to induced expression
 +
    <br>  [ARG]: The relative number of ARGs
 +
    <br>  [sg]: The relative number of sgRNA
 +
    <br>  [cell]: The OD600 of viable cell
 +
    </div>
 +
    <div class="modelingimg">
 +
      <img src="https://static.igem.org/mediawiki/2018/1/18/T--ZJUT-China--modeling10187.png" alt="">
 +
    </div>
 +
    <div class="note3">
 +
      Our ODE model of the light control ARG cutting system ver1.1.(the control group)
 +
      <br>[cell]: The OD600 of viable cell
 
     </div>
 
     </div>
    <h1>System modeling Ver3.0</h1>
 
 
     <p>
 
     <p>
       This presents an issue for researchers and factors who wish to make use of our system of the arg killing parts in a more complex environment. In order to address this issue, we decided to develop a mathematical model of how light signal effects influence the time of cutting off all args and bacteria's suicide induced by our arg killing system.
+
       By fitting the experimental data, we find that our system is feasible to predict the completion time of all the ARGs cleavage.Furthermore, we determine some constants based on experimental data which further improves our system and meets the requirements of the construction of the next model.
      In the light-controlled arg cutting system ver3.0 we want to add 3 repressors, one more sgRNA, and a lysin gene to make our cells can cut off all the args and then suicide in  chronological order which is designed.
+
    <br>
+
 
     </p>
 
     </p>
    <div class="modelingmath">
+
        <h1>System modeling Ver2.0</h1>
      d[express_r1]/t=kf*[vector_1]/(vm+[light signal])
+
    <p>
      <br>d[degrade_r1]/t=vm*[repressor_1]/([repressor_1]+km)
+
       However, if researchers would like to use our system in the complex environment,System modeling Ver1.1 can not meet their requirements. In order to address this issue, we have managed to develop a mathematical model considering the influence of light signal and some other factors.In other words,we regard the factors which include the influence of three repressors, sgRNA and lysin as parameters of system modeling ver2.0 .
      <br>d[express_r2]/t=kf*[vector_2]
+
    </p>
      <br>d[degrade_r2]/t=vm*[repressor_2]/([repressor_2]+km)
+
    <div class="modelingimg">
       <br>d[express_r3]/t=kf*[vector_1]
+
       <img src="https://static.igem.org/mediawiki/2018/c/c5/T--ZJUT-China--modeling10188.png" alt="">
      <br>d[degrade_r3]/t=vm*[repressor_3]/([repressor_3]+km)
+
      <br>d[express_cas]/t=kf*[vector_1]*(vm-[repressor_1])
+
      <br>d[express_sg1]/t=0
+
      <br>d[express_sg2]/t=kf*[vector_2]*(vm-[repressor_2])
+
       <br>d[express_ly]/t=kf*[chromosome]*(vm-[repressor_3])
+
      <br>d[cut off_1]/t=kf*[arg_1]*[sg_1]*[cas9]/(km+[arg_1])
+
      <br>d[cut off_2]/t=kf*[arg_2]*[sg_2]*[cas9]/(km+[arg_2])
+
    </div>
+
    <div class="note2">
+
      Our ODE model of the light-controlled arg cutting system ver3.0. In this model we considered the expression and degrading rate of the repressor and the effect of the copy numbers of two vectors.
+
 
     </div>
 
     </div>
 +
    <div class="note3">
 +
      Our ODE model of the light control ARG cutting system ver2.0. In this model we considered the express and degrade rate of the Repressors and the effect of the copy number of 2vectors.
 +
    <br>  [Repressor-1]: The relative number of the first Repressor protein we used in the system
 +
    <br>  [Repressor-2]: The relative number of the second Repressor protein we used in the system
 +
    <br>  [Repressor-3]: The relative number of the third Repressor protein we used in the system
 +
    <br>  [vector-1]: The relative copy number of the vector with the light-controlled promoter, repressor gene-2, cas9 gene, sgRNA gene-1, ARG-2, repressor gene-3 on it
 +
    <br>  [vector-2]: The relative copy number of the vector with ARG-1, repressor gene-2 on it
 +
    <br>  [chromosome]: The relative copy number of the chromosome
 +
    <br>  [Cas9]: The relative number of Cas9 protein
 +
    <br>  [sg-1]: The relative number of the first sgRNA
 +
    <br>  [sg-2]: The relative number of the second sgRNA
 +
    <br>  [ARG-1]: The relative number of the first ARG
 +
    <br>  [ARG-2]: The relative number of the second ARG
 +
    <br>  [lysin]: The relative number of the lysin protein
 +
    <br>  [cell]: The OD600 of viable cell
 +
        </div>
 
     <p>
 
     <p>
       Through changing the copy numbers of two vectors in our system model, we could  get a ratio of the copy numbers letting  our system do nothing before we give the light signal, while having a higher reacting rate when we want it to work.
+
       Through change the copy number of two vectors in our system model, we can get a ratio of the copy number let our system do nothing before we give the light signal, while has a faster react rate when we want it work.
 
     </p>
 
     </p>
 
     <h1>Result</h1>
 
     <h1>Result</h1>
 
     <div class="modelingimg">
 
     <div class="modelingimg">
       <img src="https://static.igem.org/mediawiki/2018/4/4c/T--ZJUT-China--modeling3.png" alt="">
+
       <img src="https://static.igem.org/mediawiki/2018/1/1f/T--ZJUT-China--modeling10189.png" alt="">
 
     </div>
 
     </div>
     <div class="note">
+
     <div class="note">Our system spend 9t on cutting off all the ARGs and 23t for lysis, when time no less than 1t and the strength of light signal is 0.1.
      Figure 1:when Time>=1t the strength of light signal=0.1,it use 12.5t to cut off all antibiotic resistance gene.
+
 
     </div>
 
     </div>
 
     <div class="modelingimg">
 
     <div class="modelingimg">
       <img src="https://static.igem.org/mediawiki/2018/4/4a/T--ZJUT-China--modeling4.png" alt="">
+
       <img src="https://static.igem.org/mediawiki/2018/0/07/T--ZJUT-China--modeling101810.png" alt="">
 
     </div>
 
     </div>
     <div class="note">
+
     <div class="note">Our system spend 8t on cutting off all the ARGs and 21t for lysis, when time no less than 1t and the strength of light signal is 1.</div>
       Figure 2:when Time>=1t the strength of light signal=1,it use 12t to cut off all antibiotic resistance gene.
+
    <div class="modelingimg">
 +
       <img src="https://static.igem.org/mediawiki/2018/e/e8/T--ZJUT-China--modeling101811.png" alt="">
 +
    </div>
 +
    <div class="note">Our system spend 11t on cutting off all the ARGs and 23t for lysis, when time no less than 1t , the strength of light signal is 1 and the time of stopping lighting is 2t.
 
     </div>
 
     </div>
 
     <div class="modelingimg">
 
     <div class="modelingimg">
       <img src="https://static.igem.org/mediawiki/2018/f/fa/T--ZJUT-China--modeling5.png" alt="">
+
       <img src="https://static.igem.org/mediawiki/2018/4/47/T--ZJUT-China--modeling101812.png" alt="">
 
     </div>
 
     </div>
     <div class="note">
+
     <div class="note">Our system spend 11t on cutting off all the ARGs and 28t for lysis, when time no less than 1t , the strength of light signal is 0.01 and the time of stopping lighting is 2t.
      Figure 3:when Time>=1t the strength of light signal=1, and shutdown when Time=2, it use 12t to cut off.
+
 
     </div>
 
     </div>
    <h1>
 
      Discussion
 
    </h1>
 
 
     <p>
 
     <p>
       Our modeling and analysis focused on the aim of achieving a better theoretical ground to forecast how our system work after giving the light signal.
+
       From the above results,we can clearly see that our system react faster under a very strength light signal.What’ more,we find that just a quite weak pulsed light signal is enough to start the operation of our system. According to this conclusion,our system is safe because the accidental leakage of bacteria will cut off all the ARGs and autolysis for the natural light.
       From the test of ODE model of the light-controlled arg cutting system ver3.0, we found that giving a short pulse light signal can make almost same effects as we constant light signal, so we find a plan to just give a short pulse light signal but make our system still response as fast as it was before, as well as saving energy used to give the light signal.
+
    </p>
 +
    <h1>Discussion</h1>
 +
    <p>
 +
       From the test of system modeling ver2.0, we find that the effect of giving a short pulsed light signal is almost same as a constant light signal.Therefore, we decide to give a short pulsed light signal for energy conservation.All in all, our modeling and analysis was focused to achieve a better theoretical grounding of forecasting how our system works after giving the light signal.
 
     </p>
 
     </p>
 
     <div class="note1">
 
     <div class="note1">
 
       [1]Quantitative approaches to the study of bistability in the lac operon of Escherichia coli  J R Soc Interface. 2008 Aug 6; 5(Suppl 1): S29–S39.
 
       [1]Quantitative approaches to the study of bistability in the lac operon of Escherichia coli  J R Soc Interface. 2008 Aug 6; 5(Suppl 1): S29–S39.
 
       Published online 2008 Apr 15. doi:  10.1098/rsif.2008.0086.focus
 
       Published online 2008 Apr 15. doi:  10.1098/rsif.2008.0086.focus
       [2]Combinatorial transcriptional control of the lactose operon of Escherichia coli
+
       <br>[2]Combinatorial transcriptional control of the lactose operon of <i>Escherichia coli</i>
 
       Thomas Kuhlman, Zhongge Zhang, Milton H. Saier, Jr., Terence Hwa
 
       Thomas Kuhlman, Zhongge Zhang, Milton H. Saier, Jr., Terence Hwa
 
       Proc Natl Acad Sci U S A. 2007 Apr 3; 104(14): 6043–6048. Published online 2007 Mar 21. doi: 10.1073/pnas.0606717104
 
       Proc Natl Acad Sci U S A. 2007 Apr 3; 104(14): 6043–6048. Published online 2007 Mar 21. doi: 10.1073/pnas.0606717104
Line 1,074: Line 1,166:
 
</div>
 
</div>
 
   <!-- 底部信息 -->
 
   <!-- 底部信息 -->
<div class="footer">
+
  <div class="footer">
  <div class="foottop">
+
    <div class="foottop">
    <div class="footmessage-left">
+
      <div class="footmessage-left">
      <a href="http://www.zjut.edu.cn/"><img src="https://static.igem.org/mediawiki/2018/4/4c/T--ZJUT-China--%E6%B5%99%E5%B7%A5%E5%A4%A7.jpg"  width="348px" alt="工大首页"></a>
+
        <img src="https://static.igem.org/mediawiki/2018/4/4c/T--ZJUT-China--%E6%B5%99%E5%B7%A5%E5%A4%A7.jpg"  width="348px" alt="工大首页">
      <a href="http://www.cbb.zjut.edu.cn/"><img src="https://static.igem.org/mediawiki/2018/e/e8/T--ZJUT-China--%E6%B5%99%E5%B7%A5%E5%A4%A7%E7%94%9F%E5%B7%A5%E5%AD%A6%E9%99%A2.jpg" alt="生工首页"></a>
+
        <img src="https://static.igem.org/mediawiki/2018/e/e8/T--ZJUT-China--%E6%B5%99%E5%B7%A5%E5%A4%A7%E7%94%9F%E5%B7%A5%E5%AD%A6%E9%99%A2.jpg" alt="生工首页">
 +
      </div>
 +
      <div class="footmessage-left2">
 +
        <img src="https://static.igem.org/mediawiki/2018/f/f6/T--ZJUT-China--logo5.png" alt="">
 +
        <img src="https://static.igem.org/mediawiki/2018/d/db/HQ_IDT_logo.jpg" alt="" >
 +
        <img src="https://static.igem.org/mediawiki/2018/1/1e/T--ZJUT-China--logo6.png" alt="">
 +
        <img src="https://static.igem.org/mediawiki/2018/f/fe/New-mathworks_logo.jpg" alt="">
 +
      </div>
 +
      <div class="footmessage-middle">
 +
          <img src="https://static.igem.org/mediawiki/2018/b/be/T--ZJUT-China--%E5%85%AC%E4%BC%97%E5%8F%B7%E4%BA%8C%E7%BB%B4%E7%A0%81.jpg" alt="公众号二维码">
 +
      </div>
 +
      <div class="footmessage-right">
 +
        <h3>Contact</h3>
 +
        <p>
 +
        <br><br>
 +
        College of Biotechnology and Bioengineering
 +
        <br><br>Zhejiang University of Technology
 +
        <br><br>Chaowang Rd. 18, 310014, Hangzhou, China
 +
        <br><br>E-mail: cbb@zjut.edu.cn Tel: +86-571-88320391
 +
        </p>
 +
      </div>
 
     </div>
 
     </div>
     <div class="footmessage-middle">
+
     <div class="footbottom">
        <img src="https://static.igem.org/mediawiki/2018/b/be/T--ZJUT-China--%E5%85%AC%E4%BC%97%E5%8F%B7%E4%BA%8C%E7%BB%B4%E7%A0%81.jpg" alt="公众号二维码">
+
      <div class="note">
    </div>
+
       Copyright © 2018 ZJUT-China
    <div class="footmessage-right">
+
       </div>
       <h3>Contact</h3>
+
      <br><br>
+
      <p>College of Bioengineering
+
 
+
      <br><br>Zhejiang University of Technology
+
 
+
      <br><br>310000, Hangzhou
+
 
+
       </p>
+
 
     </div>
 
     </div>
 
   </div>
 
   </div>
  <div class="footbottom">
 
    <div class="note">
 
    Copyright © 2018 ZJUT-China
 
    </div>
 
  </div>
 
</div>
 

Latest revision as of 17:40, 17 October 2018

Team:ZJUT-China - 2018.igem.org

Loading menubar.....

Team:ZJUT-China

Abstract

In order to improve the feasibility of our system utilized by other researchers , we have tried our best to develop a model to predict the time point when the all antibiotic resistance genes(ARGs) are cut off and the cell auto-lysis is going to begin. This model is made up of some rate equations which are based on experimental data of our system. Then we are able to simulate the operation of our system and gain the information we need from the analog results.

System modeling Ver1.0

For the purpose of calculating the time point when the all ARGs are cut off and the cell auto-lysis is going to begin,we have tried to construct a model used by the following equations.
Our ODE model of the light control ARG cutting system ver1.0.
[Cas9]: The relative number of Cas9 protein
[vector]: The relative copy number of the vector
[light signal]: The relative light strength
[ARG]: The relative number of ARGs
[sg]: The relative number of sgRNA

Result

The result that we gain from the system modeling ver1.0 is shown below.

System modeling Ver1.1

Our system modeling ver1.0 is limited by problems that we are unable to verify the validity of the model through the experiments.Because the quantization of cas9 protein and ARGs is difficult.Therefore,we have tried to build up another model which is named as System modeling Ver1.1. This model can be used to predict the time point which is same as the System modeling Ver1.0 by growth curve. We are able to verify our model by experiments because the growth curve is easily gained by experiments.In addition,we can also prove the validity of system model ver1.0 for the reason that the equations of this model is based on the system model ver1.0.

Our ODE model of the light control ARG cutting system ver1.1.(the experimental group with the system and ARG)
[Cas9]: The relative number of Cas9 protein
[vector]: The relative copy number of the vector
[ara]: The concentration of arabinose used to induced expression
[ARG]: The relative number of ARGs
[sg]: The relative number of sgRNA
[cell]: The OD600 of viable cell
Our ODE model of the light control ARG cutting system ver1.1.(the control group with ARG but without the system)
[cell-c]: The OD600 of viable cell

Result

When time=5h we put in different concentrations of Arabinose in the medium.
The points are the data get tested by wet lab and the lines are the data simulated by computer.

Our ODE model after fit data:

Our ODE model of the light control ARG cutting system ver1.1.(the experimental group)
[Cas9]: The relative number of Cas9 protein
[vector]: The relative copy number of the vector
[ara]: The concentration of arabinose used to induced expression
[ARG]: The relative number of ARGs
[sg]: The relative number of sgRNA
[cell]: The OD600 of viable cell
Our ODE model of the light control ARG cutting system ver1.1.(the control group)
[cell]: The OD600 of viable cell

By fitting the experimental data, we find that our system is feasible to predict the completion time of all the ARGs cleavage.Furthermore, we determine some constants based on experimental data which further improves our system and meets the requirements of the construction of the next model.

System modeling Ver2.0

However, if researchers would like to use our system in the complex environment,System modeling Ver1.1 can not meet their requirements. In order to address this issue, we have managed to develop a mathematical model considering the influence of light signal and some other factors.In other words,we regard the factors which include the influence of three repressors, sgRNA and lysin as parameters of system modeling ver2.0 .

Our ODE model of the light control ARG cutting system ver2.0. In this model we considered the express and degrade rate of the Repressors and the effect of the copy number of 2vectors.
[Repressor-1]: The relative number of the first Repressor protein we used in the system
[Repressor-2]: The relative number of the second Repressor protein we used in the system
[Repressor-3]: The relative number of the third Repressor protein we used in the system
[vector-1]: The relative copy number of the vector with the light-controlled promoter, repressor gene-2, cas9 gene, sgRNA gene-1, ARG-2, repressor gene-3 on it
[vector-2]: The relative copy number of the vector with ARG-1, repressor gene-2 on it
[chromosome]: The relative copy number of the chromosome
[Cas9]: The relative number of Cas9 protein
[sg-1]: The relative number of the first sgRNA
[sg-2]: The relative number of the second sgRNA
[ARG-1]: The relative number of the first ARG
[ARG-2]: The relative number of the second ARG
[lysin]: The relative number of the lysin protein
[cell]: The OD600 of viable cell

Through change the copy number of two vectors in our system model, we can get a ratio of the copy number let our system do nothing before we give the light signal, while has a faster react rate when we want it work.

Result

Our system spend 9t on cutting off all the ARGs and 23t for lysis, when time no less than 1t and the strength of light signal is 0.1.
Our system spend 8t on cutting off all the ARGs and 21t for lysis, when time no less than 1t and the strength of light signal is 1.
Our system spend 11t on cutting off all the ARGs and 23t for lysis, when time no less than 1t , the strength of light signal is 1 and the time of stopping lighting is 2t.
Our system spend 11t on cutting off all the ARGs and 28t for lysis, when time no less than 1t , the strength of light signal is 0.01 and the time of stopping lighting is 2t.

From the above results,we can clearly see that our system react faster under a very strength light signal.What’ more,we find that just a quite weak pulsed light signal is enough to start the operation of our system. According to this conclusion,our system is safe because the accidental leakage of bacteria will cut off all the ARGs and autolysis for the natural light.

Discussion

From the test of system modeling ver2.0, we find that the effect of giving a short pulsed light signal is almost same as a constant light signal.Therefore, we decide to give a short pulsed light signal for energy conservation.All in all, our modeling and analysis was focused to achieve a better theoretical grounding of forecasting how our system works after giving the light signal.

[1]Quantitative approaches to the study of bistability in the lac operon of Escherichia coli J R Soc Interface. 2008 Aug 6; 5(Suppl 1): S29–S39. Published online 2008 Apr 15. doi:  10.1098/rsif.2008.0086.focus
[2]Combinatorial transcriptional control of the lactose operon of Escherichia coli Thomas Kuhlman, Zhongge Zhang, Milton H. Saier, Jr., Terence Hwa Proc Natl Acad Sci U S A. 2007 Apr 3; 104(14): 6043–6048. Published online 2007 Mar 21. doi: 10.1073/pnas.0606717104