Team:OUC-China/Project Modeling

From 2014.igem.org

(Difference between revisions)
 
(43 intermediate revisions not shown)
Line 192: Line 192:
     <div id="mainContent" style="width:1024px;background-color:#d2e9fd;padding:10px;margin:0 auto">
     <div id="mainContent" style="width:1024px;background-color:#d2e9fd;padding:10px;margin:0 auto">
         <h1 id="Overview">Overview</h1>
         <h1 id="Overview">Overview</h1>
-
             <p>In this part,we analysed the chain suicide loop of cell lysis and establish the mathematical model,which is divided into the following three parts.First,mechanism analysis is to analyse the mechanism of the relationships among aTc, tet protein and cracked cells in suicide loop.Second, through Modeling and solving we identified the reasonable aTc concentration range, and established the ODE models according to the reaction mechanism,eventually got  the relationship of cell lysis , aTc concentration and the time.We also analysed the experimental results using pBAD instead of aTc.Last,we gave 2 promotions of the model,which respectively are the general mathematical description of the chain loop and the new circuit regulated by quorum principle.</p>
+
             <p>In this part, we analysed the chain suicide loop of cell lysis and establish the mathematical model, which is divided into the following three parts. First, mechanism analysis is to analyse the mechanism of the relationships among anhydrotetracycline(aTc), tet protein and cracked cells in suicide loop. Second, through Modeling and solving we identified the reasonable aTc concentration range, and established the ODE models according to the reaction mechanism, eventually got the relationship of cell lysis, aTc concentration and the time. We also analysed the experimental results using L-arabinose instead of aTc. Last, we gave 2 promotions of the model, which respectively are the general mathematical description of the chain loop and the new circuit regulated by quorum principle.</p>
       <h1 id="Mechanism_Analysis">Mechanism analysis</h1>
       <h1 id="Mechanism_Analysis">Mechanism analysis</h1>
-
           <h4>Background<h4>
+
           <h2 class="text-primary">Background</h2>
-
               <p>In order to make both TAT- H4 fusion protein produced by the engineering bacterium  and the transfected plasmid to be released in full,we designed corresponding suicide gene in the transfected plasmid,which has to express at the appropriate time. We used chain suicide loop in our project,which is simple and direct.However, the inductive substance aTc itself as a kind of antibiotics will affect the cell itself in the process of adding, and thus make factors complicated to affect the reaction result, which is not conducive to the analysis of the model.</p>
+
               <p>In order to make both TAT- H4 fusion protein produced by the engineering bacterium and the transfected plasmid to be released in full, we designed corresponding suicide gene in the transfected plasmid, which has to express at the appropriate time. We used chain suicide loop in our project, which is simple and direct. However, the inductive substance aTc itself as a kind of antibiotics will affect the cell itself in the process of adding, and thus make factors complicated to affect the reaction result, which is not conducive to the analysis of the model.</p>
-
           <h2>Introduction of the suicide loop</h2>  
+
           <h2 class="text-primary">Introduction of the suicide loop</h2>  
-
               <p></p>
+
               <p><img src="https://static.igem.org/mediawiki/2014/1/14/OUC-China_Project_Modeling_1.jpg"style="display:block;margin:0 auto" /></p>
-
               <p>On the one hand,tetR protein itself can repress the promoter R0040, that is to say,under normal conditions R0040 promoter is repressed and cell lysis does not happen.But the gene expression of tetR protein is restrained when people add the aTc , thus the number of tetR proteins in the cells will fall,and the repression of promoter R0040 respective will drop,leading cell lysis device to start and the cells to crack.On the other hand,aTc itself also can lead to cell lysis which brings difficulties to the analysis of the experimental result. As a result,we should first identify the aTc concentration range to ensure  aTc itself make less effect on cells, and then we can do quantitative analysis of the experimental process within this range.</p>
+
               <p>On the one hand, tetR protein itself can repress the promoter R0040, that is to say, under normal conditions R0040 promoter is repressed and cell lysis does not happen. But the gene expression of tetR protein is restrained when people add the aTc, thus the number of tetR proteins in the cells will fall, and the repression of promoter R0040 respective will drop, leading cell lysis device to start and the cells to crack. On the other hand, aTc itself also can lead to cell lysis which brings difficulties to the analysis of the experimental result. As a result, we should first identify the aTc concentration range to ensure  aTc itself make less effect on cells, and then we can do quantitative analysis of the experimental process within this range.</p>
-
          <h2>Two main problems</h2>
+
            <h2 class="text-primary">Two main problems</h2>
-
               <p>According to the analysis above, we put forward the following two major problems to solve:
+
               <p>According to the analysis above, we put forward the following two major problems to solve:</p>
-
Problem 1: We need to identify a reasonable aTc concentration range to make aTc itself have less effect on the engineering bacteria but also can effectively induce the expression of the lysis gene.
+
<h4 class="text-muted">Problem 1:</h4>
-
problem 2: We need to establish a reasonable ODE model to describe the relations of each material in the process of lysis and to guide the biological experiments.</p>
+
<p> We need to identify a reasonable aTc concentration range to make aTc itself have less effect on the engineering bacteria but also can effectively induce the expression of the lysis gene.</p>
 +
<h4 class="text-muted">problem 2: </h4>
 +
<p>We need to establish a reasonable ODE model to describe the relations of each material in the process of lysis and to guide the biological experiments.</p>
       <h1 id="Modeling_And_Solving">Modeling and solving</h1>
       <h1 id="Modeling_And_Solving">Modeling and solving</h1>
-
         <h2> Model assumes</h2>
+
         <h2 class="text-primary"> Model assumes</h2>
-
             <p> assume 1: There is only one factor that influences the suicide process in the experiment without any other affection.</p>
+
             <h4 class="text-muted"> assume 1: </h4>
-
             <p>assume 2: There is only one process of the chain reaction without other bifurcate structures.</p>
+
                  <p>There is only one factor that influences the suicide process in the experiment without any other affection.</p>
-
         <h2>Symbol Description</h2>
+
             <h4 class="text-muted"> assume 2: </h4>
-
             <p></p>
+
            <p>There is only one process of the chain reaction without other bifurcate structures.</p>
-
         <h2>Problem 1:To identify the aTc concentration range</h2>   
+
         <h2 class="text-primary">Symbol Description</h2>
 +
             <p><img src="https://static.igem.org/mediawiki/2014/1/16/OUC-China_Project_Modeling_2.png"style="display:block;margin:0 auto" /></p>
 +
         <h2 class="text-primary">Problem 1:To identify the aTc concentration range</h2>   
             <p>In order to find out reasonable aTc concentration range to make sure the cell death is mainly associated with the suicide mechanism.Therefore we designed a set of control experiment to To detect the reasonable aTc concentration.</p>
             <p>In order to find out reasonable aTc concentration range to make sure the cell death is mainly associated with the suicide mechanism.Therefore we designed a set of control experiment to To detect the reasonable aTc concentration.</p>
-
            <p>图 待画</p>
+
             <p>We suppose there are two groups of <em>E.coli</em> cells with similar cell activity (A and B). Engineered <em>E. coli</em> cells of group A contained artificial plasmid, and the aTc can trigger the cell lysis. While another <em>E. coli</em> cells of group B does not contain artificial plasmid, whose growth is directly affected by aTc itself. According to assume 1, there is only one factor that influences the suicide process in the experiment without any other affection. Therefore <em>E.coli</em> cells of group B will be directly affected by aTc and result in cell lysis, while engineered <em>E.coli</em> cells of group A  will be affected by both aTc and cell lysis.<p>
-
             <p>We suppose there are two groups of E.coli cells with similar cell activity (A and B).Engineered E. coli cells of group A contained artificial plasmid,and the aTc can trigger the cell lysis.While another E. coli cells of group B does not contain artificial plasmid,whose growth directly affected by aTc itself.According to assume 1, there is only one factor that influences the suicide process in the experiment without any other affection.Therefore E.coli cells of group B will be directly affected by aTc and result in cell lysis, while engineered E.coli cells of group A  will be affected by both aTc and cell lysis.<p>
+
             <p>In order to accurately determine the suitable aTc concentration,we introduced a quantity <img src="https://static.igem.org/mediawiki/2014/a/ab/OUC-China_Project_Modeling_3.png"/>,which is the ratio of numbers of cells of group B and group A per unit time under different aTc concentrations:</p>
-
             <p>In order to accurately determine the suitable aTc concentration,we introduced a quantity图,which is the ratio of numbers of cells of group B and group A per unit time under different aTc concentrations:</p>
+
             <p><img src="https://static.igem.org/mediawiki/2014/7/79/%E5%9B%BE%E7%89%8712.png"style="display:block;margin:0 auto" /></p>
-
             <p></p>
+
             <p>The formula shows that the main reason of cell death can be attributed to the influence of suicide system when <img src="https://static.igem.org/mediawiki/2014/a/ab/OUC-China_Project_Modeling_3.png"/>is larger .So far, we have presented a set of method to detect reasonable aTc concentration range.According to the results UCB has given in 2012 iGEM,it's better when aTc concentration range is between 0 and 140 nanogram per milliliter.The numerical simulation diagram is given.</p>
-
             <p>The formulas shows that the main reason of cell death can be attributed to the influence of suicide system when 图is larger .So far, we have presented a set of method to detect reasonable aTc concentration range.According to the results UCB has given in 2012 iGEM,it's better when aTc concentration range is between 0 and 140 nanogram per milliliter.The numerical simulation diagram is given.</p>
+
               <p><img src="https://static.igem.org/mediawiki/2014/9/96/OUC-China_Project_Modeling_126.png"style="display:block;margin:0 auto" /></p>
-
               <p></p>
+
               <p>We carried out on the numerical fitting by using this result and got the figure :</p>
               <p>We carried out on the numerical fitting by using this result and got the figure :</p>
-
               <p></p>
+
               <p> <img src="https://static.igem.org/mediawiki/2014/b/b5/OUC-China_Project_Modeling_euler%26yulei_we_are_a_team.png"style="display:block;margin:0 auto" /></p>
-
               <p>Thus we can analyse and conclude that, 图 is larger when the aTc concentration is from 80 to 110 grams per milliliter, at the same time,it illustrates the suicide mechanism we designed works well when the aTc concentration is in this range. That is to say, it can be more powerful to kill cells and control numbers of cells.</P>
+
               <p>Thus we can analyse and conclude that, <img src="https://static.igem.org/mediawiki/2014/a/ab/OUC-China_Project_Modeling_3.png"/> is larger when the aTc concentration is from 80 to 110 grams per milliliter, at the same time, it illustrates the suicide mechanism we designed works well when the aTc concentration is in this range. That is to say, it can be more powerful to kill cells and control numbers of cells.</P>
-
           <h2>Problem 2:ODE Model</h2>
+
           <h2 class="text-primary">Problem 2:ODE Model</h2>
-
           <p>Anhydrotetracycline aTc inhibits the expression of tetR gene,so that in the part of the tetR protein producing, we deduced that the productive rate of the tetR protein is inversely proportional to aTc concentration,and got the equation. </p>
+
           <p>aTc inhibits the expression of tetR gene,so that in the part of the tetR protein producing, we deduced that the productive rate of the tetR protein is inversely proportional to aTc concentration,and got the equation. </p>
-
           <p></p>
+
           <p><img src="https://static.igem.org/mediawiki/2014/f/f2/OUC-China_Project_Modeling_113.png"style="display:block;margin:0 auto" /></p>
           <p>At the same time, tetR protein also inhibits its production,therefore we can identify:</p>
           <p>At the same time, tetR protein also inhibits its production,therefore we can identify:</p>
-
           <p></p>
+
           <p><img src="https://static.igem.org/mediawiki/2014/a/a4/OUC-China_Project_Modeling_114.png"style="display:block;margin:0 auto" /></p>
           <p>Combined with the two equations we got:</p>
           <p>Combined with the two equations we got:</p>
-
           <p><p>
+
           <p><img src="https://static.igem.org/mediawiki/2014/b/bf/OUC-China_Project_Modeling_115.png"style="display:block;margin:0 auto" /><p>
<p>TetR protein produced by the first step in the chain reaction can inhibit the promoter of lysis unit to express,so we described the relation between the rate of lysis and tetR protein concentration in terms of the equation:</p>
<p>TetR protein produced by the first step in the chain reaction can inhibit the promoter of lysis unit to express,so we described the relation between the rate of lysis and tetR protein concentration in terms of the equation:</p>
-
           <p></p>
+
           <p><img src="https://static.igem.org/mediawiki/2014/d/d6/OUC-China_Project_Modeling_17.png"style="display:block;margin:0 auto" /></p>
           <p>Sp is the amount of the cracked cells, is the reaction order of the combination between tetR protein and the lysis device which need obtaining by the experiment. as an unknown variable, needs to be discussed below.Because aTc itself can lead to cell lysis,we determined that C is a function of aTc concentration and can stimulate cell lysis, so we concluded that:</p>
           <p>Sp is the amount of the cracked cells, is the reaction order of the combination between tetR protein and the lysis device which need obtaining by the experiment. as an unknown variable, needs to be discussed below.Because aTc itself can lead to cell lysis,we determined that C is a function of aTc concentration and can stimulate cell lysis, so we concluded that:</p>
-
           <p></p>
+
           <p><img src="https://static.igem.org/mediawiki/2014/4/4a/OUC-China_Project_Modeling_118.png"style="display:block;margin:0 auto" /></p>
           <p>So in the part of cell lysis,we set up the following equation to describe the reaction of suicide loop:</p>
           <p>So in the part of cell lysis,we set up the following equation to describe the reaction of suicide loop:</p>
-
           <p></p>   
+
           <p><img src="https://static.igem.org/mediawiki/2014/4/45/OUC-China_Project_Modeling_119.png"style="display:block;margin:0 auto" /></p>   
           <p>Simultaneous equations above to establish system of equations:</p>
           <p>Simultaneous equations above to establish system of equations:</p>
-
           <p></p>
+
           <p><img src="https://static.igem.org/mediawiki/2014/4/4a/OUC-China_Project_Modeling_120.png"style="display:block;margin:0 auto" /></p>
           <p>System of equations (3) can describe the reaction mechanism of the suicide system. First we solved equation (1),and obtained the relation between the tetR protein concentration under the different aTc concentration and the time:</p>
           <p>System of equations (3) can describe the reaction mechanism of the suicide system. First we solved equation (1),and obtained the relation between the tetR protein concentration under the different aTc concentration and the time:</p>
-
           <p></p>
+
           <p><img src="https://static.igem.org/mediawiki/2014/5/5b/OUC-China_Project_Modeling_121.png"style="display:block;margin:0 auto" /></p>
           <p>Among this a2 controls the initial reaction rate,and c1 is an undetermined constant.We substituted equation (4) into equation (2) to solve the differential equation, and obtained the function of the amount of cracked cells :</p>
           <p>Among this a2 controls the initial reaction rate,and c1 is an undetermined constant.We substituted equation (4) into equation (2) to solve the differential equation, and obtained the function of the amount of cracked cells :</p>
-
           <p></p>
+
           <p><img src="https://static.igem.org/mediawiki/2014/9/9f/OUC-China_Project_Modeling_122.png"style="display:block;margin:0 auto" /></p>
-
           <p>We can carry on the numerical simulation By (4) (5) .In order to calculate conveniently we set图,and drew a diagram about the relation between the tetR protein concentration under the different aTc concentration and the time,which is as shown in the figure 1:</p>
+
           <p>We can carry on the numerical simulation By (4) (5) .In order to calculate conveniently we set<img src="https://static.igem.org/mediawiki/2014/f/f4/OUC-China_Project_Modeling_14.png"/>,and drew a diagram about the relation between the tetR protein concentration under the different aTc concentration and the time,which is as shown in the figure 1:</p>
-
           <p></p>
+
           <p><img src="https://static.igem.org/mediawiki/2014/2/2d/OUC-China_Project_Modeling_11.png"style="display:block;margin:0 auto" /></p>
-
           <p>From the diagram above we can see that the tetR concentration falls as aTc concentration rises in the range of 80 to 110 .We can further analyse the variation law of the number of cracked cells with time under different aTc concentration in the case that the tetR concentration is known.As well,in order to calculate conveniently, we set图and drew a diagram about it,which is as shown in figure 2:</p>
+
           <p>From the diagram above we can see that the tetR concentration falls as aTc concentration rises in the range of 80 to 110 .We can further analyse the variation law of the number of cracked cells with time under different aTc concentration in the case that the tetR concentration is known.As well,in order to calculate conveniently, we set<img src="https://static.igem.org/mediawiki/2014/e/ea/OUC-China_Project_Modeling_128.png"/>and drew a diagram about it,which is as shown in figure 2:</p>
-
           <p></p>
+
           <p><img src="https://static.igem.org/mediawiki/2014/e/e6/OUC-China_Project_Modeling_b2.png"style="display:block;margin:0 auto" /></p>
-
           <p>The diagram above shows that the amount of the cracked engineering bacterium cells increases with time grows and the trend becomes smooth.It also increases with the aTc concentration increases.  
+
           <p>The diagram above shows that the amount of the cracked engineering bacterium cells increases with time grows and the trend becomes smooth.It also increases with the anhydrotetracycline(aTc) concentration increases.</p>
-
To sum up, we have completed the quantitative analysis  in every process of suicide system, and given the  mathematical description of the product and changes in the biological reaction process. We also obtained the reasonable scope of aTc concentration and  the number of cracked cells with time under different aTc concentration,and thus completed the establishment and solution of the model.
+
<p>To sum up, we have completed the quantitative analysis  in every process of suicide system, and given the  mathematical description of the product and changes in the biological reaction process. We also obtained the reasonable scope of aTc concentration and the number of cracked cells with time under different aTc concentration,and thus completed the establishment and solution of the model.
-
Through the model we can theoretically give the mathematic expression of the aTc concentration and the amount of cracked cells, but the coefficients need to be determined in the subsequent experiments,after which we can estimate the amount of cracked cells via the aTc additive amount to predict the results of experiments.And we can also control the aTc additive amount via the amount of cracked cells as a threshold value.</p>
+
Through the model we can theoretically give the mathematic expression of the aTc concentration and the amount of cracked cells, but the coefficients need to be determined in the subsequent experiments, after which we can estimate the amount of cracked cells via the aTc additive amount to predict the results of experiments. And we can also control the aTc additive amount via the amount of cracked cells as a threshold value.</p>
-
         <h2>Analysis of the experimental results using pBAD instead of aTc</h2>
+
         <h2 class="text-primary">Analysis of the experimental results using L-arabinose instead of aTc</h2>
-
         <p>We used pBAD to replace aTc and then conducted the experiment.pBAD can induce lysis device but not kill cells directly.Experimental datas are as follows: </p>
+
         <p>We used L-arabinose to replace aTc and then conducted the experiment.L-arabinose can induce lysis device but not kill cells directly.Experimental datas are as follows: </p>
-
          <p></p>
+
<style>
 +
.centertable td { text-align: center; }
 +
</style>
 +
<h4 style="text-align:center">Absorbance at 600nm</h4>
 +
<table class="table table-striped table-hover centertable" style="border:2px solid #a0a0a0">
 +
<thead>
 +
<tr>
 +
  <td rowspan="2" style="width:350px"><strong>Concentration of L-arabinose</strong></td>
 +
  <td colspan="3"><strong>Time after induction</strong></td>
 +
</tr>
 +
<tr>
 +
  <td>3 hours</td>
 +
  <td>6 hours</td>
 +
  <td>9 hours</td>
 +
</tr>
 +
<tr>
 +
  <td>1μM</td>
 +
  <td>0.322</td>
 +
  <td>0.745</td>
 +
  <td>0.901</td>
 +
</tr>
 +
<tr>
 +
  <td>3μM</td>
 +
  <td>0.247</td>
 +
  <td>0.741</td>
 +
  <td>0.927</td>
 +
</tr>
 +
<tr>
 +
  <td>6μM</td>
 +
  <td>0.209</td>
 +
  <td>0.733</td>
 +
  <td>0.891</td>
 +
</tr>
 +
<tr>
 +
  <td>10μM</td>
 +
  <td>0.191</td>
 +
  <td>0.750</td>
 +
  <td>0.924</td>
 +
</tr>
 +
<tr>
 +
  <td>100μM</td>
 +
  <td>0.183</td>
 +
  <td>0.721</td>
 +
  <td>0.848</td>
 +
</tr>
 +
<tr>
 +
  <td>1M</td>
 +
  <td>0.222</td>
 +
  <td>0.548</td>
 +
  <td>0.687</td>
 +
</tr>
 +
<tr>
 +
  <td>10mM</td>
 +
  <td>0.132</td>
 +
  <td>0.631</td>
 +
  <td>0.685</td>
 +
</tr>
 +
</table>
           <p>According to the data we drew the bar chart and the line chart: </p>
           <p>According to the data we drew the bar chart and the line chart: </p>
-
           <p></p>
+
           <p><img src="https://static.igem.org/mediawiki/2014/8/80/OUC-China_Project_Modeling_19.png"style="display:block;margin:0 auto" /></p>
-
           <p></p>
+
           <p><img src="https://static.igem.org/mediawiki/2014/3/3d/OUC-China_Project_Modeling_20.jpg"/style="display:block;margin:0 auto" /><P>In order to demonstrate the reliability of the experiment, we used one-way analysis of variance to analyse the experiment above and got the diagram below: </p>
-
          <p>In order to demonstrate the reliability of the experiment, we used one-way analysis of variance to analyse the experiment above and got the diagram below: </p>
+
           <p><img src="https://static.igem.org/mediawiki/2014/b/bd/OUC-China_Project_Modeling_yu29.png"style="display:block;margin:0 auto" /></p>
-
           <p></p>
+
           <p>Through one-way analysis of variance about the OD600 after 3 hours, 6 hours and 9 hours, we can see that the addition of L-arabinose significantly induced the cell lysis and the effect is remarkable.Thus we got the conclusion that L-arabinose can lead to cell lysis and have a good experimental effect.</p>
-
           <p>Through one-way analysis of variance about the OD600 after 3 hours, 6 hours and 9 hours, we can see that the addition of pBAD significantly induced the cell lysis and the effect is remarkable.Thus we got the conclusion that pBAD can lead to cell lysis and have a good experimental effect.</p>
+
       <h1 id="Promotion">Promotion</h1>
       <h1 id="Promotion">Promotion</h1>
-
           <h2>The analysis of the model</h2>
+
           <h2 class="text-primary">The analysis of the model</h2>
-
             <p>The model can well reflect the suicide mechanism and The effects of aTc content on the cell cracking,which has a certain practical value.At the same time,it can be concluded as a method suitable for more experiments,and complete the quantitative description of the experiment.However,the model is relatively simple,only considering the effects of the single variable on the cells, but in reality, the activities of the cell and the changes of its number not only controlled by one variable, and it is not given a better strategy of the random interference factors in the experiment. So it need improving in the future research.</p>
+
             <p>The model can well reflect the suicide mechanism and the effects of aTc content on the cell cracking, which has a certain practical value.At the same time,it can be concluded as a method suitable for more experiments,and complete the quantitative description of the experiment. However, the model is relatively simple, only considering the effects of the single variable on the cells, but in reality, the activities of the cell and the changes of its number not only controlled by one variable, and it is not given a better strategy of the random interference factors in the experiment. So it needs improving in the future research.</p>
-
          <h2>Promotions of the model</h2>
+
        <h2 class="text-primary">Promotions of the model</h2>
-
           <h3>Promotion 1:</h3>
+
           <h4 class="text-muted">Promotion 1:</h4>
<p>On the basis of the experiment, we found that experiments which meet the topological structure  shown in figure 3 can take advantage of this model to describe.</p>
<p>On the basis of the experiment, we found that experiments which meet the topological structure  shown in figure 3 can take advantage of this model to describe.</p>
-
             <p></p>
+
             <p><img src="https://static.igem.org/mediawiki/2014/2/25/OUC-China_Project_Modeling_b3.png"style="display:block;margin:0 auto" /></p>
             <p>Input is a single variable which has influence on cells,and it can promote the expression of a specific purpose gene to produce the protein which can make cell lysis.At the same time,many kinds of inputs are antibiotic substances, therefore the input itself can cause cracking of cells directly.For such experiments it can be described by the equation:</p>
             <p>Input is a single variable which has influence on cells,and it can promote the expression of a specific purpose gene to produce the protein which can make cell lysis.At the same time,many kinds of inputs are antibiotic substances, therefore the input itself can cause cracking of cells directly.For such experiments it can be described by the equation:</p>
-
             <p></p>
+
             <p><img src="https://static.igem.org/mediawiki/2014/e/eb/OUC-China_Project_Modeling_b23.png"style="display:block;margin:0 auto" /></p>
             <p>Thus we can accomplish to solve a class of problems and promote the models.</p>
             <p>Thus we can accomplish to solve a class of problems and promote the models.</p>
-
            <h3>Promotion 2:</h3>
+
          <h4 class="text-muted">Promotion 2:</h4>
<p>On this basis, we further promoted the system regulated by quorum sensing principle,which makes the new suicide circuit work better.
<p>On this basis, we further promoted the system regulated by quorum sensing principle,which makes the new suicide circuit work better.
Symbol Description</p>
Symbol Description</p>
-
               <p></p>
+
               <p><img src="https://static.igem.org/mediawiki/2014/f/f6/OUC-China_Project_Modeling_24.png"style="display:block;margin:0 auto" /></p>
               <p>The following response structure is given: </p>
               <p>The following response structure is given: </p>
-
               <p></p>
+
               <p><img src="https://static.igem.org/mediawiki/2014/b/b6/OUC-China_Project_Modeling_b5.png"style="display:block;margin:0 auto" /></p>
               <p>This system can better control lysis compared with the chain reaction before,and we gave the ODE model to describe the reactions:</p>    
               <p>This system can better control lysis compared with the chain reaction before,and we gave the ODE model to describe the reactions:</p>    
-
               <p></p>
+
               <p><img src="https://static.igem.org/mediawiki/2014/8/89/OUC-China_Project_Modeling_yu24.png"style="display:block;margin:0 auto" /></p>
-
<p>We used numerical method to solve the equations. According to the experience,we set </p>
+
<p>We used numerical method to solve the equations. According to the experience,we set<img src="https://static.igem.org/mediawiki/2014/d/df/OUC-China_Project_Modeling_yu25.png"/></p>
<p>Diagram is as follows when PH=10:</p>   
<p>Diagram is as follows when PH=10:</p>   
-
<p></p>
+
<p><img src="https://static.igem.org/mediawiki/2014/d/de/OUC-China_Project_Modeling_yu6.png"style="display:block;margin:0 auto" /></p>
<p>We set PH = 1 under the premise of other coefficients invariant staying constant and got the diagram below: </p>
<p>We set PH = 1 under the premise of other coefficients invariant staying constant and got the diagram below: </p>
-
<p></p>
+
<p><img src="https://static.igem.org/mediawiki/2014/a/a2/OUC-China_Project_Modeling_yu7.png"style="display:block;margin:0 auto" /></p>
<p>Thus we can see that PH can effectively affect the content of three kinds of variables , and then control cell lysis.  
<p>Thus we can see that PH can effectively affect the content of three kinds of variables , and then control cell lysis.  
-
In oder to reflect  the interaction of three kinds of variables intuitivly,we can draw the diagram among the three,which is as follows (PH= 1) </p>
+
In oder to reflect  the interaction of three kinds of variables intuitivly,we can draw the diagram among the three,which is as follows (PH=1) </p>
-
<p>And we can also drew the diagram of any two variables:</p>
+
<p><img src="https://static.igem.org/mediawiki/2014/c/cd/OUC-China_Project_Modeling_yuyu8.png"style="display:block;margin:0 auto" /></p>
-
<p></p>
+
<p>And we can also draw the diagram of any two variables:</p>
-
<p></p>
+
<p><img src="https://static.igem.org/mediawiki/2014/4/4a/OUC-China_Project_Modeling_yuyuyu9.png"style="display:block;margin:0 auto" /></p>
-
<p></p>
+
<p><img src="https://static.igem.org/mediawiki/2014/5/59/OUC-China_Project_Modeling_bibibibibibi10.png"style="display:block;margin:0 auto" /></p>
 +
<p><img src="https://static.igem.org/mediawiki/2014/7/72/%E5%9B%BE%E7%89%8711.png"style="display:block;margin:0 auto" /></p>
 +
<h3>References</h3> 
 +
<p>[1] Thomas Lederer,Martin Kintrup,Masayuki Takahashi, Phaik-Eng Sum,George A. Ellestad,and Wolfgang Hillen, Tetracycline Analogs Affecting Binding to Tn10-Encoded Tet Repressor Trigger the Same Mechanism of Induction, <em>Biochemistry</em> 1996, 35, 7439-7446.</p>
 +
<p>[2] Casonya M. Johnson and Robert F. Schleif,In Vivo Induction Kinetics of the Arabinose Promoters in Escherichia coli, <em>Journal Of Bacteriology</em>, June 1995, p. 3438–3442.</p>
     </div>
     </div>
Line 316: Line 379:
         $('#childNavigator').affix({ offset: { top: 285 } });
         $('#childNavigator').affix({ offset: { top: 285 } });
     </script>
     </script>
 +
<div class = "footer">
 +
  <style>
 +
    .footer { padding: 10px; background-color: rgb(45,111,165); }
 +
    .footer > div { width: 418px; margin: 0 auto; padding-top: 50px; }
 +
    .footer > div > img { float: left; height: 100px; display: inline-block; position: relative; margin-left: 6px; margin-right: 6px; }
 +
    .footer > p { text-align: center; color: #ffffff }
 +
    .footer > p > a { color: #ffffff; }
 +
  </style>
 +
  <div>
 +
    <img src="https://static.igem.org/mediawiki/2014/a/ad/OUC-China_Index_OrgLogo.png" />
 +
    <img src="https://static.igem.org/mediawiki/2014/e/ec/OUC-China_Index_OrgTitle.jpg" />
 +
    <img src="https://static.igem.org/mediawiki/2014/9/96/OUC-China_Foot_ouclogo.png" />
 +
    <div style="clear:both"></div>
 +
  </div>
 +
  <p>contact us: <a href="mailto:oucigem@163.com">oucigem@163.com</a></p>
 +
</div>
<img class="anchor" onclick="document.documentElement.scrollTop = document.body.scrollTop = 0;" src="https://static.igem.org/mediawiki/2014/b/b6/ToTop.png">
<img class="anchor" onclick="document.documentElement.scrollTop = document.body.scrollTop = 0;" src="https://static.igem.org/mediawiki/2014/b/b6/ToTop.png">
</html>
</html>

Latest revision as of 03:12, 18 October 2014

Overview

In this part, we analysed the chain suicide loop of cell lysis and establish the mathematical model, which is divided into the following three parts. First, mechanism analysis is to analyse the mechanism of the relationships among anhydrotetracycline(aTc), tet protein and cracked cells in suicide loop. Second, through Modeling and solving we identified the reasonable aTc concentration range, and established the ODE models according to the reaction mechanism, eventually got the relationship of cell lysis, aTc concentration and the time. We also analysed the experimental results using L-arabinose instead of aTc. Last, we gave 2 promotions of the model, which respectively are the general mathematical description of the chain loop and the new circuit regulated by quorum principle.

Mechanism analysis

Background

In order to make both TAT- H4 fusion protein produced by the engineering bacterium and the transfected plasmid to be released in full, we designed corresponding suicide gene in the transfected plasmid, which has to express at the appropriate time. We used chain suicide loop in our project, which is simple and direct. However, the inductive substance aTc itself as a kind of antibiotics will affect the cell itself in the process of adding, and thus make factors complicated to affect the reaction result, which is not conducive to the analysis of the model.

Introduction of the suicide loop

On the one hand, tetR protein itself can repress the promoter R0040, that is to say, under normal conditions R0040 promoter is repressed and cell lysis does not happen. But the gene expression of tetR protein is restrained when people add the aTc, thus the number of tetR proteins in the cells will fall, and the repression of promoter R0040 respective will drop, leading cell lysis device to start and the cells to crack. On the other hand, aTc itself also can lead to cell lysis which brings difficulties to the analysis of the experimental result. As a result, we should first identify the aTc concentration range to ensure  aTc itself make less effect on cells, and then we can do quantitative analysis of the experimental process within this range.

Two main problems

According to the analysis above, we put forward the following two major problems to solve:

Problem 1:

We need to identify a reasonable aTc concentration range to make aTc itself have less effect on the engineering bacteria but also can effectively induce the expression of the lysis gene.

problem 2:

We need to establish a reasonable ODE model to describe the relations of each material in the process of lysis and to guide the biological experiments.

Modeling and solving

Model assumes

assume 1:

There is only one factor that influences the suicide process in the experiment without any other affection.

assume 2:

There is only one process of the chain reaction without other bifurcate structures.

Symbol Description

Problem 1:To identify the aTc concentration range

In order to find out reasonable aTc concentration range to make sure the cell death is mainly associated with the suicide mechanism.Therefore we designed a set of control experiment to To detect the reasonable aTc concentration.

We suppose there are two groups of E.coli cells with similar cell activity (A and B). Engineered E. coli cells of group A contained artificial plasmid, and the aTc can trigger the cell lysis. While another E. coli cells of group B does not contain artificial plasmid, whose growth is directly affected by aTc itself. According to assume 1, there is only one factor that influences the suicide process in the experiment without any other affection. Therefore E.coli cells of group B will be directly affected by aTc and result in cell lysis, while engineered E.coli cells of group A  will be affected by both aTc and cell lysis.

In order to accurately determine the suitable aTc concentration,we introduced a quantity ,which is the ratio of numbers of cells of group B and group A per unit time under different aTc concentrations:

The formula shows that the main reason of cell death can be attributed to the influence of suicide system when is larger .So far, we have presented a set of method to detect reasonable aTc concentration range.According to the results UCB has given in 2012 iGEM,it's better when aTc concentration range is between 0 and 140 nanogram per milliliter.The numerical simulation diagram is given.

We carried out on the numerical fitting by using this result and got the figure :

Thus we can analyse and conclude that,  is larger when the aTc concentration is from 80 to 110 grams per milliliter, at the same time, it illustrates the suicide mechanism we designed works well when the aTc concentration is in this range. That is to say, it can be more powerful to kill cells and control numbers of cells.

Problem 2:ODE Model

aTc inhibits the expression of tetR gene,so that in the part of the tetR protein producing, we deduced that the productive rate of the tetR protein is inversely proportional to aTc concentration,and got the equation.

At the same time, tetR protein also inhibits its production,therefore we can identify:

Combined with the two equations we got:

TetR protein produced by the first step in the chain reaction can inhibit the promoter of lysis unit to express,so we described the relation between the rate of lysis and tetR protein concentration in terms of the equation:

Sp is the amount of the cracked cells, is the reaction order of the combination between tetR protein and the lysis device which need obtaining by the experiment. as an unknown variable, needs to be discussed below.Because aTc itself can lead to cell lysis,we determined that C is a function of aTc concentration and can stimulate cell lysis, so we concluded that:

So in the part of cell lysis,we set up the following equation to describe the reaction of suicide loop:

Simultaneous equations above to establish system of equations:

System of equations (3) can describe the reaction mechanism of the suicide system. First we solved equation (1),and obtained the relation between the tetR protein concentration under the different aTc concentration and the time:

Among this a2 controls the initial reaction rate,and c1 is an undetermined constant.We substituted equation (4) into equation (2) to solve the differential equation, and obtained the function of the amount of cracked cells :

We can carry on the numerical simulation By (4) (5) .In order to calculate conveniently we set,and drew a diagram about the relation between the tetR protein concentration under the different aTc concentration and the time,which is as shown in the figure 1:

From the diagram above we can see that the tetR concentration falls as aTc concentration rises in the range of 80 to 110 .We can further analyse the variation law of the number of cracked cells with time under different aTc concentration in the case that the tetR concentration is known.As well,in order to calculate conveniently, we setand drew a diagram about it,which is as shown in figure 2:

The diagram above shows that the amount of the cracked engineering bacterium cells increases with time grows and the trend becomes smooth.It also increases with the anhydrotetracycline(aTc) concentration increases.

To sum up, we have completed the quantitative analysis  in every process of suicide system, and given the  mathematical description of the product and changes in the biological reaction process. We also obtained the reasonable scope of aTc concentration and the number of cracked cells with time under different aTc concentration,and thus completed the establishment and solution of the model. Through the model we can theoretically give the mathematic expression of the aTc concentration and the amount of cracked cells, but the coefficients need to be determined in the subsequent experiments, after which we can estimate the amount of cracked cells via the aTc additive amount to predict the results of experiments. And we can also control the aTc additive amount via the amount of cracked cells as a threshold value.

Analysis of the experimental results using L-arabinose instead of aTc

We used L-arabinose to replace aTc and then conducted the experiment.L-arabinose can induce lysis device but not kill cells directly.Experimental datas are as follows:

Absorbance at 600nm

Concentration of L-arabinose Time after induction
3 hours 6 hours 9 hours
1μM 0.322 0.745 0.901
3μM 0.247 0.741 0.927
6μM 0.209 0.733 0.891
10μM 0.191 0.750 0.924
100μM 0.183 0.721 0.848
1M 0.222 0.548 0.687
10mM 0.132 0.631 0.685

According to the data we drew the bar chart and the line chart:

In order to demonstrate the reliability of the experiment, we used one-way analysis of variance to analyse the experiment above and got the diagram below:

Through one-way analysis of variance about the OD600 after 3 hours, 6 hours and 9 hours, we can see that the addition of L-arabinose significantly induced the cell lysis and the effect is remarkable.Thus we got the conclusion that L-arabinose can lead to cell lysis and have a good experimental effect.

Promotion

The analysis of the model

The model can well reflect the suicide mechanism and the effects of aTc content on the cell cracking, which has a certain practical value.At the same time,it can be concluded as a method suitable for more experiments,and complete the quantitative description of the experiment. However, the model is relatively simple, only considering the effects of the single variable on the cells, but in reality, the activities of the cell and the changes of its number not only controlled by one variable, and it is not given a better strategy of the random interference factors in the experiment. So it needs improving in the future research.

Promotions of the model

Promotion 1:

On the basis of the experiment, we found that experiments which meet the topological structure  shown in figure 3 can take advantage of this model to describe.

Input is a single variable which has influence on cells,and it can promote the expression of a specific purpose gene to produce the protein which can make cell lysis.At the same time,many kinds of inputs are antibiotic substances, therefore the input itself can cause cracking of cells directly.For such experiments it can be described by the equation:

Thus we can accomplish to solve a class of problems and promote the models.

Promotion 2:

On this basis, we further promoted the system regulated by quorum sensing principle,which makes the new suicide circuit work better. Symbol Description

The following response structure is given:

This system can better control lysis compared with the chain reaction before,and we gave the ODE model to describe the reactions:

 

We used numerical method to solve the equations. According to the experience,we set

Diagram is as follows when PH=10:

We set PH = 1 under the premise of other coefficients invariant staying constant and got the diagram below:

Thus we can see that PH can effectively affect the content of three kinds of variables , and then control cell lysis. In oder to reflect  the interaction of three kinds of variables intuitivly,we can draw the diagram among the three,which is as follows (PH=1)

And we can also draw the diagram of any two variables:

References

 

[1] Thomas Lederer,Martin Kintrup,Masayuki Takahashi, Phaik-Eng Sum,George A. Ellestad,and Wolfgang Hillen, Tetracycline Analogs Affecting Binding to Tn10-Encoded Tet Repressor Trigger the Same Mechanism of Induction, Biochemistry 1996, 35, 7439-7446.

[2] Casonya M. Johnson and Robert F. Schleif,In Vivo Induction Kinetics of the Arabinose Promoters in Escherichia coli, Journal Of Bacteriology, June 1995, p. 3438–3442.