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

 
(43 intermediate revisions by 4 users not shown)
Line 2: Line 2:
  
 
<html>
 
<html>
<style>
+
<body style="width:1263px; margin-left:auto;margin-right:auto;background-color:white" onscroll="scrollEvent()">
    body {
+
        width: 100%;
+
        margin-left: auto;
+
        margin-right: auto;
+
    }
+
 
+
    #body {
+
        width: 1263px;
+
        margin-left: auto;
+
        margin-right: auto;
+
        background-color: aliceblue;
+
    }
+
 
+
    #comic p {
+
        background-color: aliceblue;
+
        text-align: center;
+
    }
+
 
+
    #comic img {
+
        margin-top: 50px;
+
        width: 1263px;
+
        height: 620px;
+
        margin-left: auto;
+
        margin-right: auto;
+
    }
+
 
+
    /*just to eliminate any margin */
+
    #basic_description {
+
        margin: 0 0 0 0;
+
    }
+
 
+
        /*also eliminate any margin*/
+
        #basic_description p {
+
            margin: 0px;
+
        }
+
 
+
    /*put footer in the center!*/
+
    #_footer p {
+
        text-align: center;
+
    }
+
 
+
    .description {
+
        width: 1000px;
+
        margin-left: auto;
+
        margin-right: auto;
+
        font-family: Arial, Helvetica, sans-serif;
+
        font-size: large;
+
        padding-bottom: 30px;
+
    }
+
 
+
        .description li, .description p {
+
            color: #404040;
+
        }
+
 
+
        .description h3 {
+
            font-family: 'Arial, Helvetica, sans-serif';
+
            font-size: 30px;
+
            font-weight: 100px;
+
            line-height: 30px;
+
        }
+
 
+
        .description p {
+
            text-indent: 25px !important;
+
            font-family: Arial, Helvetica, sans-serif !important;
+
            font-size: 25px !important;
+
        }
+
 
+
        .description #ref li {
+
            font-style: italic;
+
            font-size: x-small;
+
            color: black;
+
        }
+
 
+
        .description table {
+
            margin-top: 20px;
+
            margin-left: 20px;
+
            margin-right: 20px;
+
            border: thin solid black;
+
            caption-side: top;
+
            border-spacing: 0px;
+
        }
+
 
+
        .description td, .description th {
+
            border: thin dotted gray;
+
            padding: 10px;
+
            background-color: aliceblue;
+
            font-weight: 200;
+
            line-height: 30px;
+
        }
+
 
+
 
+
        .description #pricetable table {
+
            text-align: center;
+
            margin-left: 200px;
+
            width: 600px;
+
        }
+
 
+
    .basic_description .description #price th {
+
        padding: 15px;
+
    }
+
 
+
    .description #price caption {
+
        padding-top: 8px;
+
    }
+
 
+
    .insert {
+
        text-align: center;
+
        background-color: aliceblue;
+
    }
+
 
+
        .insert p img {
+
            width: 1200px;
+
            margin-left: -100px;
+
        }
+
 
+
    /* 左侧列表 */
+
    .paper-sidenav {
+
        position: fixed;
+
        left: 0.1em;
+
        top: 15em;
+
        text-align: left;
+
        width: 20em;
+
        font-family: Arial, Helvetica, sans-serif;
+
        background-color: rgba(8,8,8,0);
+
    }
+
 
+
        .paper-sidenav ul {
+
            list-style-type: none;
+
            list-style: none;
+
            padding-inline-start: 0.5em;
+
            margin-left: 0px;
+
        }
+
 
+
            .paper-sidenav ul p {
+
                color: black;
+
                font-size: 1.6em !important;
+
                font-weight: 500;
+
            }
+
 
+
            .paper-sidenav ul a:link {
+
                text-decoration: none;
+
                color: black;
+
            }
+
 
+
            .paper-sidenav ul a:visited {
+
                text-decoration: none;
+
                color: black;
+
            }
+
 
+
            .paper-sidenav ul a:hover {
+
                text-decoration: none;
+
                color: blue;
+
                background-color: #f9f9f9;
+
                font-size: larger;
+
            }
+
 
+
            .paper-sidenav ul a:active {
+
                text-decoration: none;
+
                color: black;
+
            }
+
 
+
    .description .topic-title {
+
        padding-top: 50px;
+
    }
+
 
+
    .reference {
+
        width: 1000px;
+
        margin-left: auto;
+
        margin-right: auto;
+
        font-family: Arial, Helvetica, sans-serif;
+
        padding-bottom: 30px;
+
    }
+
 
+
    .hidden-link {
+
        text-decoration: none;
+
    }
+
</style>
+
<body style="width:1263px; margin-left:auto;margin-right:auto;background-color:white">
+
 
     <div id="body">
 
     <div id="body">
 
         <!--default for floating navigation-->
 
         <!--default for floating navigation-->
         <div class="paper-sidenav">
+
         <div id="sidebar" class="paper-sidenav">
 
             <ul>
 
             <ul>
 
                 <li>
 
                 <li>
 
                     <ul>
 
                     <ul>
                         <li><a href="#section1">Bound Fraction &Unbound Fraction</a></li>
+
                         <li><a href="#section1">Bound Fraction & Unbound Fraction</a></li>
 
                         <li><a href="#section2">Modelling Equations</a></li>
 
                         <li><a href="#section2">Modelling Equations</a></li>
                         <li><a href="#section3">Methods</a></li>
+
                         <li><a href="#section3">Methods&Results</a></li>
                         <li><a href="#section4">Results</a></li>
+
                         <li><a href="#section4">Discussion</a></li>
                        <li><a href="#section5">Discussion</a></li>
+
 
                     </ul>
 
                     </ul>
 
                 </li>
 
                 </li>
Line 198: Line 19:
 
         <div id="comic">
 
         <div id="comic">
 
             <p>
 
             <p>
                 <img src= alt="model" />
+
                 <!-- <img src= alt="model" /> -->
 
             </p>
 
             </p>
 
         </div>
 
         </div>
 
         <div class="description">
 
         <div class="description">
             <h1 style="text-align: center;">Modelling</h1>
+
             <h1>Modelling</h1>
 
             <div class="topic-title" id="section1">
 
             <div class="topic-title" id="section1">
                 <h3>Bound Fraction &Unbound Fraction</h3>
+
                 <h3>Bound Fraction & Unbound Fraction</h3>
 
                 <p> By defining a simple equation where P is the promoter region and X is the ligand that
 
                 <p> By defining a simple equation where P is the promoter region and X is the ligand that
can bind with P, namely a repressor in this case: </p>
+
can bind with P, which is an activator or a repressor: </p>
 
               <div style="text-align:center;"><img src="https://static.igem.org/mediawiki/2018/c/c4/T--RDFZ-China--reaction.jpeg"/></div>
 
               <div style="text-align:center;"><img src="https://static.igem.org/mediawiki/2018/c/c4/T--RDFZ-China--reaction.jpeg"/></div>
               <p> Intuitively, the factional saturation of the promoter is denoted by Ps , which is the fraction of
+
               <p> The factional saturation of the promoter is denoted by Ps , which is the fraction of
 
occupied binding sites: </p>
 
occupied binding sites: </p>
 
               <div style="text-align:center;"><img src="https://static.igem.org/mediawiki/2018/3/3d/T--RDFZ-China--equation11.jpeg"/style="text-align:center;"></div>
 
               <div style="text-align:center;"><img src="https://static.igem.org/mediawiki/2018/3/3d/T--RDFZ-China--equation11.jpeg"/style="text-align:center;"></div>
Line 216: Line 37:
 
               <p>Then we rearrange this equation and get: </p>
 
               <p>Then we rearrange this equation and get: </p>
 
               <div style="text-align:center;"><img src="https://static.igem.org/mediawiki/2018/6/63/T--rdfZ-China--equation13.jpeg"/></div>
 
               <div style="text-align:center;"><img src="https://static.igem.org/mediawiki/2018/6/63/T--rdfZ-China--equation13.jpeg"/></div>
 +
              <p> Where kd is the dissociation constant of this binding event. Now we substitute equation 1.3 into equation 1.1 and hence obtain: </p>
 +
              <div style="text-align:center"><img src="https://static.igem.org/mediawiki/2018/7/77/T--rdfZ-China--equation14.jpeg"/></div>
 +
              <p>Assuming that the sum of the bound and the unbound equals one, the unbound fraction Pu is thus: </p>
 +
              <div style="text-align:center;"><img src="https://static.igem.org/mediawiki/2018/f/f2/T--rdfZ-China--equation15.jpeg"/></div>
 +
              <p>Bound fraction is considered for activation in gene regulation and unbound fraction for repressive regulation. </p>
 
             </div>
 
             </div>
 
             <div class="topic-title" id="section2">
 
             <div class="topic-title" id="section2">
 
                 <h3>Modelling Equations</h3>
 
                 <h3>Modelling Equations</h3>
                 <p> We provided team Tsinghua with IPTG assay with Ptac plasmid backbone in Tsinghua’s IPTG induction devices in their NEON system.Also we provided AHL for their Lux system.Tsinghua was generous to provide us some lab space and equipments to work on our projects. We worked close, both physically and intellectually, to create maximum impacts on our team projects.</p>
+
                <div style="text-align:center;"><img src="https://static.igem.org/mediawiki/2018/e/ec/T--RDFZ-China--Equation_1.png"/></div>
 +
                 <p>
 +
                The chemical AHL (A) is required by LuxR to activate transcription. The active form for LuxR it consists of a dimer of the
 +
                complex LuxR-AHL ((Y:A)2). Thus we can consider the production of this dimer as an elementary reaction: 2A+2Y --> (Y:A)2,
 +
                whose rate depends upon A, the concentration of AHL, and the concentration of unbound LuxR monomer, Runbound.
 +
                </p>
 +
                <div style="text-align:center;"><img src="https://static.igem.org/mediawiki/2018/a/a4/T--RDFZ-China--Equation_2.png"/></div>
 +
                <p>
 +
                Where synthesis is the rate of synthesis of the dimer, k1 is the binding coefficient RT is the total concentration of LuxR
 +
                monomer, and R*(t) is the concentration of AHL-LuxR dimer.
 +
                And k2R*(t) is simply the rate of degradation of the dimer, where k2 is the degradation coefficient of AHL-LuxR.
 +
                </p>
 +
                <div style="text-align:center;"><img src="https://static.igem.org/mediawiki/2018/b/bc/T--RDFZ-China--Equation_3.png"/></div>
 +
                <p>
 +
                For any given mRNA:
 +
                </p>
 +
                <div style="text-align:center;"><img src="https://static.igem.org/mediawiki/2018/a/ad/T--RDFZ-China--Equation_4.png"></div>
 +
                <p>
 +
                Where m(t) is the concentration of mRNA, k0 is the amount of mRNA produced, δm is the degradation rate of mRNA.
 +
                For PhIF protein translated from the given mRNA:
 +
                </p>
 +
                <div style="text-align:center;"><img src="https://static.igem.org/mediawiki/2018/2/20/T--RDFZ-China--Equation_5.png"/></div>
 +
                <p>
 +
                Where P(t) is the concentration of PhIF protein, k3 is the translation rate, and δp is the degradation rate of PhIF
 +
                protein.
 +
                Then substitute equation 1.4 into the calculation of k0 (for the activation of Plux promoter by AHL-LuxR correspond to  
 +
                bound fraction):
 +
                </p>
 +
                <div style="text-align:center;"><img src="https://static.igem.org/mediawiki/2018/d/db/T--RDFZ-China--Equation_6.png"/></div>
 +
                <p>
 +
                Where a is the transcription rate, k2 is the equivalent of kd in equation 1.4 and 1.5.
 +
                </p>
 +
                <div style="text-align:center;"><img src="https://static.igem.org/mediawiki/2018/3/32/T--RDFZ-China--Equation_7.png"/></div>
 +
                <p>
 +
                For sfGFP protein translated from the given mRNA:
 +
                </p>
 +
                <div style="text-align:center;"><img src="https://static.igem.org/mediawiki/2018/5/50/T--RDFZ-China--Equation_8.png"/></div>
 +
                <p>
 +
                Where G(t) is the concentration of sfGFP protein, k3 is the translation rate, and δg is the degradation rate of PhIF
 +
                protein.
 +
                Then substitute equation 1.5 into the calculation of k0 in this case (for the repression of PPhIF promoter by PhIF protein
 +
                correspond to the unbound fraction):
 +
                </p>
 +
                <div style="text-align:center;"><img src="https://static.igem.org/mediawiki/2018/2/24/T--RDFZ-China--Equation_9.png"/></div>
 +
                <p>
 +
                Where a is the transcription rate, k2 is the equivalent of kd in equation 1.4 and 1.5.
 +
                </p>
 +
 
 +
 
 
             </div>
 
             </div>
 
             <div class="topic-title" id="section3">
 
             <div class="topic-title" id="section3">
                 <h3>Methods</h3>
+
                 <h3>Methods&Results</h3>
                 <p> This year, we provided SJTU-BioX-Shanghai with plasmids expressing acoustic reporter genes arg1 and arg2. They were applied to report the location of cancer foci combining with ultrasound in their project “E.Coli for Colon Health Observation”. We also shared protocols and equipments regarding non-invasive bacteria monitoring in our section of gas vesicles.</p>
+
                 <p>
 +
              To simulate our model, we use the Matlab ODE suite to solve three differential equations simultaneously. A breif diagram of the regulatory network can be found below:                         
 +
              </p>
 +
              <div style="text-align:center;"><img src="https://static.igem.org/mediawiki/2018/1/1e/T--RDFZ-China--quorum_sensing_model.jpeg"/></div>
 +
              <p>
 +
              By adjusting the initial dose of AHL, ranging from 10^-5 M to 10^-12 M , we can get a variance of dose response below:
 +
              </p>
 +
              <div style="text-align:center;"><img src="https://static.igem.org/mediawiki/2018/a/a2/T--RDFZ-China--parameterScanGif.gif"/></div>
 +
              <p>From the simulation results, we can see that between [AHL]=10^-10 and 10^-11 M, the concentration of the sfGFP(which represents ccdb) rises above 5*10^-8M which gives the half maximal killing rate, meaning a large proportion of bacteria in the fermentor would die. Therefore, the worker in the factory should always be aware not having [AHL] below that range.</p>
 +
              <div style="text-align:center;"><img src="https://static.igem.org/mediawiki/2018/1/17/T--RDFZ-China--AHL_parameter_scan.jpeg"/width="80%;"></div>
 +
              <p>A local sensitivity analysis is also carried out to help us increase the robustness of the model. Two influential factors are selected to perform a further parameter scan.</p>
 +
              <div style="text-align:center;"><img src="https://static.igem.org/mediawiki/2018/9/9a/T--RDFZ-China--sensitivity_matrix.jpeg"/width="80%;"></div>
 +
              <div style="text-align:center;"><img src="https://static.igem.org/mediawiki/2018/1/14/T--RDFZ-China--sensitivity_figure.jpeg"/width="80%;"></div>
 +
              <div style="text-align:center;padding:20px,0,10px,0;"><img src="https://static.igem.org/mediawiki/2018/6/6d/T--RDFZ-China--scanning_figure_K1.jpeg"/width="50%;"><img src="https://static.igem.org/mediawiki/2018/8/83/T--RDFZ-China--scanning_figure_K2.jpeg"/width="50%;"></div>
 +
              <p>The scanning results suggest that the variation of the parameter, K1 and K2,  by 10% respectively will not affect the trend.</p> 
 +
              <p>Furthermore, our results from the wetlab indicates a similar trend, corroborating our results from modelling</p>
 +
                <div style="text-align:center;"><img src="https://static.igem.org/mediawiki/2018/2/28/T--RDFZ-China--Lux-characterization.png"/width="80%;"></div>
 
             </div>
 
             </div>
 
             <div class="topic-title" id="section4">
 
             <div class="topic-title" id="section4">
                <h3>Results</h3>
 
                <p> Zeyu Tang from the team, RDFZ-China, kindly provided help on modeling for team GBC. To be specific, Zeyu gave them one document about previous work he has done on the fermentor model, which contains deterministic ODEs for predicting the productivity of engineered bacteria in the fermentor. In addition, Zeyu wrote another document which illustrates the deduction of Hill equation, helping them create their model simulation. </p>
 
                <p> Knowing GreatBay was constructing a constitutive expressing GFP into the bacterial genome, we enquired for it and used it as a capacity monitor for characterizing the expression burden. Greatbay kindly gave it to us as a gift.
 
            </div>
 
            <div class="topic-title" id="section5">
 
 
                 <h3>Discussion</h3>
 
                 <h3>Discussion</h3>
                 <p> Our team leader in 2017, Bowen Xiao, did a report on strategies preparing the iGEM competition; then the two teams shared resources and had a conversation on the biggest problems encounter during iGEM. Team BFSUICC-China in September met with us to discuss project updates and Jamboree preparation. </p>
+
                 <p> In a real factory, there will be two fermentors, one for the production of AHL which helps keeping the bacteria in the second fermentor alive, one for the production of the main product. Theoretically, by measuring the OD in the first fermentor, the concentration of AHL can be inferred by using an algorithm shown below.</p>
 +
              <div style="text-align:center;"><img src="https://static.igem.org/mediawiki/2018/4/4a/T--RDFZ-China--concept_proof.png"/width="80%;"></div>
 +
              <div style="text-align:center;"><img src="https://static.igem.org/mediawiki/2018/8/8c/T--RDFZ-China--proof_concept_label.png"/width="80%;"></div>
 +
              <p> In conclusion, our model sheds light on the amount of concentration of AHL in the fermentor that is necessary for keeping the bacteria alive. The model simulation results correspond well to the data obtained from the wet-lab experiment. In addition, we also attempt to derive an algorithm to help measure the AHL concentration in the other fermentor. </p>
 
             </div>
 
             </div>
 
         </div>
 
         </div>
 +
      <div class="reference">
 +
                <div class="topic-title" id="">
 +
                    <h3>References</h3>
 +
                    <p>[1] Basu, S., Mehreja, R., Thiberge, S., Chen, M., & Weiss, R. (2004). Spatiotemporal control of gene expression with pulse-generating networks. Proceedings Of The National Academy Of Sciences, 101(17), 6355-6360. doi:10.1073/pnas.0307571101</p>
 +
                    <p>[2] Ingalls, Brian P. Mathematical Modeling in Systems Biology: an Introduction. MIT Press, 2013. </p>
 +
                    <p>[3]"Luxi - Acyl-Homoserine-Lactone Synthase - Aliivibrio Fischeri - Luxi Gene & Protein". 2018. Uniprot.Org. Accessed October 18 2018. https://www.uniprot.org/uniprot/P12747.</p>
 +
                    <p>[4]Rodrigo, G., Carrera, J., Elena, S. F., & Jaramillo, A. (2010). Robust dynamical pattern formation from a multifunctional minimal genetic circuit. BMC systems biology, 4(1), 48.</p>
 +
                    <p>[5] Smith, R., Tan, C., Srimani, J., Pai, A., Riccione, K., Song, H., & You, L. (2014). Programmed Allee effect in bacteria causes a tradeoff between population spread and survival. Proceedings Of The National Academy Of Sciences, 111(5), 1969-1974. doi:10.1073/pnas.1315954111 </p>
 +
                </div>
 +
            </div>
 
</body>
 
</body>
 +
<script>
 +
function scrollEvent(){
 +
if(document.documentElement.scrollTop > 300) {
 +
  document.getElementById("sidebar").style.display="block";
 +
} else {
 +
  document.getElementById("sidebar").style.display="none";
 +
}
 +
}
 +
</script>
 
</html>
 
</html>
 
{{RDFZ-China/footer}}
 
{{RDFZ-China/footer}}

Latest revision as of 10:07, 13 November 2018

Modelling

Bound Fraction & Unbound Fraction

By defining a simple equation where P is the promoter region and X is the ligand that can bind with P, which is an activator or a repressor:

The factional saturation of the promoter is denoted by Ps , which is the fraction of occupied binding sites:

At equilibrium of this reaction, the rate of forward reaction and that of backward reaction are the same:

Then we rearrange this equation and get:

Where kd is the dissociation constant of this binding event. Now we substitute equation 1.3 into equation 1.1 and hence obtain:

Assuming that the sum of the bound and the unbound equals one, the unbound fraction Pu is thus:

Bound fraction is considered for activation in gene regulation and unbound fraction for repressive regulation.

Modelling Equations

The chemical AHL (A) is required by LuxR to activate transcription. The active form for LuxR it consists of a dimer of the complex LuxR-AHL ((Y:A)2). Thus we can consider the production of this dimer as an elementary reaction: 2A+2Y --> (Y:A)2, whose rate depends upon A, the concentration of AHL, and the concentration of unbound LuxR monomer, Runbound.

Where synthesis is the rate of synthesis of the dimer, k1 is the binding coefficient RT is the total concentration of LuxR monomer, and R*(t) is the concentration of AHL-LuxR dimer. And k2R*(t) is simply the rate of degradation of the dimer, where k2 is the degradation coefficient of AHL-LuxR.

For any given mRNA:

Where m(t) is the concentration of mRNA, k0 is the amount of mRNA produced, δm is the degradation rate of mRNA. For PhIF protein translated from the given mRNA:

Where P(t) is the concentration of PhIF protein, k3 is the translation rate, and δp is the degradation rate of PhIF protein. Then substitute equation 1.4 into the calculation of k0 (for the activation of Plux promoter by AHL-LuxR correspond to bound fraction):

Where a is the transcription rate, k2 is the equivalent of kd in equation 1.4 and 1.5.

For sfGFP protein translated from the given mRNA:

Where G(t) is the concentration of sfGFP protein, k3 is the translation rate, and δg is the degradation rate of PhIF protein. Then substitute equation 1.5 into the calculation of k0 in this case (for the repression of PPhIF promoter by PhIF protein correspond to the unbound fraction):

Where a is the transcription rate, k2 is the equivalent of kd in equation 1.4 and 1.5.

Methods&Results

To simulate our model, we use the Matlab ODE suite to solve three differential equations simultaneously. A breif diagram of the regulatory network can be found below:

By adjusting the initial dose of AHL, ranging from 10^-5 M to 10^-12 M , we can get a variance of dose response below:

From the simulation results, we can see that between [AHL]=10^-10 and 10^-11 M, the concentration of the sfGFP(which represents ccdb) rises above 5*10^-8M which gives the half maximal killing rate, meaning a large proportion of bacteria in the fermentor would die. Therefore, the worker in the factory should always be aware not having [AHL] below that range.

A local sensitivity analysis is also carried out to help us increase the robustness of the model. Two influential factors are selected to perform a further parameter scan.

The scanning results suggest that the variation of the parameter, K1 and K2, by 10% respectively will not affect the trend.

Furthermore, our results from the wetlab indicates a similar trend, corroborating our results from modelling

Discussion

In a real factory, there will be two fermentors, one for the production of AHL which helps keeping the bacteria in the second fermentor alive, one for the production of the main product. Theoretically, by measuring the OD in the first fermentor, the concentration of AHL can be inferred by using an algorithm shown below.

In conclusion, our model sheds light on the amount of concentration of AHL in the fermentor that is necessary for keeping the bacteria alive. The model simulation results correspond well to the data obtained from the wet-lab experiment. In addition, we also attempt to derive an algorithm to help measure the AHL concentration in the other fermentor.

References

[1] Basu, S., Mehreja, R., Thiberge, S., Chen, M., & Weiss, R. (2004). Spatiotemporal control of gene expression with pulse-generating networks. Proceedings Of The National Academy Of Sciences, 101(17), 6355-6360. doi:10.1073/pnas.0307571101

[2] Ingalls, Brian P. Mathematical Modeling in Systems Biology: an Introduction. MIT Press, 2013.

[3]"Luxi - Acyl-Homoserine-Lactone Synthase - Aliivibrio Fischeri - Luxi Gene & Protein". 2018. Uniprot.Org. Accessed October 18 2018. https://www.uniprot.org/uniprot/P12747.

[4]Rodrigo, G., Carrera, J., Elena, S. F., & Jaramillo, A. (2010). Robust dynamical pattern formation from a multifunctional minimal genetic circuit. BMC systems biology, 4(1), 48.

[5] Smith, R., Tan, C., Srimani, J., Pai, A., Riccione, K., Song, H., & You, L. (2014). Programmed Allee effect in bacteria causes a tradeoff between population spread and survival. Proceedings Of The National Academy Of Sciences, 111(5), 1969-1974. doi:10.1073/pnas.1315954111

footer