Difference between revisions of "Team:UNOTT/Modelling"

Line 4: Line 4:
 
<head>
 
<head>
 
   <meta charset="UTF-8">
 
   <meta charset="UTF-8">
 
+
 
 
   <script src="https://s.codepen.io/assets/libs/modernizr.js" type="text/javascript"></script>
 
   <script src="https://s.codepen.io/assets/libs/modernizr.js" type="text/javascript"></script>
  
Line 52: Line 52:
 
}
 
}
  
.image {  
+
.image {
   position: center;  
+
   position: center;
 
   width: 100%;
 
   width: 100%;
 
}
 
}
Line 74: Line 74:
 
}
 
}
  
h2 span {  
+
h2 span {
   color: white;  
+
   color: white;
   font: 100px/100px Helvetica, Sans-Serif;  
+
   font: 100px/100px Helvetica, Sans-Serif;
   letter-spacing: 10px;
+
   letter-spacing: 10px;
 
   background: rgb(0, 0, 0); /* fallback color */
 
   background: rgb(0, 0, 0); /* fallback color */
 
   background: rgba(0, 0, 0, 0.7);
 
   background: rgba(0, 0, 0, 0.7);
   padding: 10px;  
+
   padding: 10px;
 
}
 
}
h2 {  
+
h2 {
   position: absolute;  
+
   position: absolute;
   top: 350px;  
+
   top: 350px;
   left: 0;  
+
   left: 0;
 
   width: 100%;
 
   width: 100%;
 
   font-family: 'Karla', sans-
 
   font-family: 'Karla', sans-
 
}
 
}
h3 span {  
+
h3 span {
   color: white;  
+
   color: white;
   font: 20px/20px Helvetica, Sans-Serif;  
+
   font: 20px/20px Helvetica, Sans-Serif;
   letter-spacing: 0px;
+
   letter-spacing: 0px;
 
   background: rgb(0, 0, 0); /* fallback color */
 
   background: rgb(0, 0, 0); /* fallback color */
 
   background-color: #D74214;
 
   background-color: #D74214;
   padding: 10px;  
+
   padding: 10px;
 
   float: right;
 
   float: right;
 
   position: relative;
 
   position: relative;
Line 101: Line 101:
 
}
 
}
 
.container2 {
 
.container2 {
 
+
 
 
   width: 80%;
 
   width: 80%;
 
   margin: auto;
 
   margin: auto;
Line 108: Line 108:
 
   margin-right: 100px;
 
   margin-right: 100px;
 
   padding-top: 250px;
 
   padding-top: 250px;
+
 
  
 
}
 
}
Line 231: Line 231:
  
 
<p>A major problem the project faced is that the comparison process of the fluorescence proteins wouldn't be possible to be investigated with all combinations as it would take too long. </p>
 
<p>A major problem the project faced is that the comparison process of the fluorescence proteins wouldn't be possible to be investigated with all combinations as it would take too long. </p>
<p> &nbsp; </p>  
+
<p> &nbsp; </p>
 
<p> To answer this problem, the team will attempt to model the fluorescence spectra over time expressed by the proteins given different. First, the type of gene expression would need to be identified and then, would be modified to considered the effects of inhibition and finally, be applied over time to see how much expression would occur at a certain time period. The team will use Mathematical modeling such as Ordinary Differential Equations because they are easy to convert into programming in order to build components for the simulation.</p>
 
<p> To answer this problem, the team will attempt to model the fluorescence spectra over time expressed by the proteins given different. First, the type of gene expression would need to be identified and then, would be modified to considered the effects of inhibition and finally, be applied over time to see how much expression would occur at a certain time period. The team will use Mathematical modeling such as Ordinary Differential Equations because they are easy to convert into programming in order to build components for the simulation.</p>
 
<p> &nbsp; </p>
 
<p> &nbsp; </p>
Line 238: Line 238:
 
</a>
 
</a>
 
</div>
 
</div>
 
 
 
 
  
  
Line 258: Line 254:
  
 
<sup> Figure 1 </sup>
 
<sup> Figure 1 </sup>
$$  sfGFP \underset{Transcriptin}{\rightarrow} mRNA \underset{Translation}{\rightarrow} sfGFP $$  
+
$$  sfGFP \underset{Transcriptin}{\rightarrow} mRNA \underset{Translation}{\rightarrow} sfGFP $$
<p> The equation above describes the process of which the gene undergoes transcription to produce mRNA. The mRNA carries the genetic information copied from the DNA which codes for protein. The expression of protein, can therefore, be measured by the fluorescence which is the desired output of the system. </p>  
+
<p> The equation above describes the process of which the gene undergoes transcription to produce mRNA. The mRNA carries the genetic information copied from the DNA which codes for protein. The expression of protein, can therefore, be measured by the fluorescence which is the desired output of the system. </p>
  
 
<sup> Figure 2 </sup>
 
<sup> Figure 2 </sup>
Line 270: Line 266:
 
<p>  We can apply Law of Mass Action combine both equations for the concentration of protein and mRNA over time. This model can be described as: </p>
 
<p>  We can apply Law of Mass Action combine both equations for the concentration of protein and mRNA over time. This model can be described as: </p>
  
<sup> Figure 3 </sup>  
+
<sup> Figure 3 </sup>
 
$$ mRNA = k_{1} -d _{1 } mRNA  $$
 
$$ mRNA = k_{1} -d _{1 } mRNA  $$
 
$$ Protein = k_{2} \cdot  mRNA - d_{2} \cdot Protein $$
 
$$ Protein = k_{2} \cdot  mRNA - d_{2} \cdot Protein $$
Line 279: Line 275:
 
   <li>Protein is the concentration of Protein</li>
 
   <li>Protein is the concentration of Protein</li>
 
   <li>k<sub> 1 </sub> is the constitutive transcription rate. This represents the number of mRNA molecules produced per gene, per unit of time.</li>
 
   <li>k<sub> 1 </sub> is the constitutive transcription rate. This represents the number of mRNA molecules produced per gene, per unit of time.</li>
<li> d <sub> 1 </sub>  is  the  mRNA  degradation  rate </li>  
+
<li> d <sub> 1 </sub>  is  the  mRNA  degradation  rate </li>
 
   <li>k<sub> 2 </sub> is the translation rate.  This represents the number of protein molecules produced per mRNA molecule, per unit of time.</li>
 
   <li>k<sub> 2 </sub> is the translation rate.  This represents the number of protein molecules produced per mRNA molecule, per unit of time.</li>
<li> d <sub> 2 </sub>  is  the  protein  degradation  rate. </li>  
+
<li> d <sub> 2 </sub>  is  the  protein  degradation  rate. </li>
 
</ul>
 
</ul>
<br> </br>  
+
<br> </br>
  
<p> This is important because we can use this model to calculate the concentration of proteins we can expect over time. This is useful as we can use this information to calculate the total emitted light spectra during the time period which is what we are looking for in our system. However, the constants and variables are individual for each protein and which means parameters for each protein would need to be found. These constants were found using literature <sup> 3 </sup> (for GFP) and lab results (the rest.) </p>  
+
<p> This is important because we can use this model to calculate the concentration of proteins we can expect over time. This is useful as we can use this information to calculate the total emitted light spectra during the time period which is what we are looking for in our system. However, the constants and variables are individual for each protein and which means parameters for each protein would need to be found. These constants were found using literature <sup> 3 </sup> (for GFP) and lab results (the rest.) </p>
  
 
<br> </br>
 
<br> </br>
Line 297: Line 293:
 
</a>
 
</a>
 
</div>
 
</div>
    </div>
 
 
</div>
 
</div>
  
 
     <div class="expandable-box">
 
     <div class="expandable-box">
       <h4 style="color: #ffffff; font-weight: bold; font-size: 30px;">STEP 3: Promoter Library</h4><center><img class="icons2" src="https://static.igem.org/mediawiki/2017/f/f1/T--UNOTT--Promoterpool.png" style="width:300px;height:auto;"></center>
+
 
 +
       <h4 style="color: #ffffff; font-weight: bold; font-size: 30px;">Gene Transcription Regulation by Repressors (CRISPRi) - Concentration over Time</h4><center><img class="icons2" src="https://static.igem.org/mediawiki/2017/f/f1/T--UNOTT--Promoterpool.png" style="width:300px;height:auto;"></center>
 
   <div id="clear4" style="display: none;">
 
   <div id="clear4" style="display: none;">
      <p></p>
+
    <p> The next step in developing our simulation was to calculate our protein concentration at any given time when using CRISPRi. Discussion with wet-lab revealed our method would be using CRISPRi as a repressor, which works by inhibiting the expression of one or more genes by binding to the promoter region <sup> 1 </sup>. The expanded mRNA and Protein concentration models from the Constitutive Gene Expression Model <sup> 2 </sup> were modified to include the element of repression from the CRISPRi inhibition. </p>
      <link href="https://fonts.googleapis.com/css?family=Roboto" rel="stylesheet">
+
  
 +
    $$ Gene \overset{Repressor}{\rightarrow} mRNA \rightarrow Protein  $$
  
<div class="clickbutton">
+
    $$  mRNA \underset{Degradation}{\rightarrow} \oslash  $$
<a href="https://2017.igem.org/Team:UNOTT/Experiments">
+
            <p>BACK</p>
+
</a>
+
</div>
+
  
 +
    $$  sfGFP \underset{Degradation}{\rightarrow} \oslash  $$
  
<div class="clickbutton1">
+
    <h5> This change can be applied to the Law of Mass Action <sup> 3 </sup> : </h5>
<a href="https://2017.igem.org/Team:UNOTT/Results">
+
    $$ m = k_{1} \cdot \frac{k^{n}}{k^{n} + R^{n}}- d_{1}m $$
            <h3><span>RESULTS</span></h3>
+
    $$ p = k_{2} m - d_{2}p $$
</a>
+
 
</div>
+
    <p>Where...</p>
 +
 
 +
    <p>m is mRNA concentration, p is Protein concentration, R is Repressor, k1 is Max Transcription Rate, k is the Repression Coefficient, n is number of repressors that need to cooperatively bind the promoter to trigger the inhibition of gene expression (Hill Coefficient), R is Repressor, d1 is mRNA degradation rate, d2 is Protein degradation rate </p>
 +
 
 +
    <p> The value for these constants and variables were taken from literature and calculating them <sup> 4 </sup> but later, adjusted to the lab results.</p>
 +
 
 +
    <sup> Figure 6 </sup>
 +
    <br> </br>
 +
    <img src="https://static.igem.org/mediawiki/2017/7/73/T--UNOTT--InhibitedAndNon.png" class="border" height="350" width="550" style= position: fixed; align=center; >
 +
 
 +
    <br> </br>
 +
    <p> Figure 6 shows the structure which underwent CRISPRi inhibition are expected to produce lower concentration of the protein whose expression were are inhibiting. This is important as it means the team can calculate concentration of proteins which are inhibited and compare them to the control conditions as well as giving the correct concentration for the simulation. </p>
 +
 
 +
    <p> Furthermore, by having a model which can calculate the protein concentration at any given time, we can deduce how much fluorescence is being emitted at that time period by the bacteria </p>
 +
    <p> <sup> 4 </sup> See Relationship between Max Fluorescence and Protein Concentration </p>
 +
   
 
     </div>
 
     </div>
 
   </div>
 
   </div>
Line 344: Line 352:
 
       <h4 style="color: #ffffff; font-weight: bold; font-size: 30px;">STEP 5: Freeze Drying & Revival</h4><center><img class="icons4" src="https://static.igem.org/mediawiki/2017/7/73/T--UNOTT--FreezeDrying.png" style="width:150px;height:auto;"></center>
 
       <h4 style="color: #ffffff; font-weight: bold; font-size: 30px;">STEP 5: Freeze Drying & Revival</h4><center><img class="icons4" src="https://static.igem.org/mediawiki/2017/7/73/T--UNOTT--FreezeDrying.png" style="width:150px;height:auto;"></center>
 
       <div id="clear6" style="display: none;">
 
       <div id="clear6" style="display: none;">
       <p>For <i>Key. coli</i> to work as intended and not deteriorate we need need to things to occur:  
+
       <p>For <i>Key. coli</i> to work as intended and not deteriorate we need need to things to occur:
 
<br>
 
<br>
 
<ol>
 
<ol>
 
<li><span>The <i> E. coli </i> cells must be kept inactive so that nutrients is not depleted causing the transformed cells to die </span></li>
 
<li><span>The <i> E. coli </i> cells must be kept inactive so that nutrients is not depleted causing the transformed cells to die </span></li>
<li><span> The <i> E. coli </i> cells must be able to be activated after inactivation to allow the fluorescent genes to be expressed to give the key its unique fluorescent code which will allow access to the appliance.</span></li>  
+
<li><span> The <i> E. coli </i> cells must be able to be activated after inactivation to allow the fluorescent genes to be expressed to give the key its unique fluorescent code which will allow access to the appliance.</span></li>
 
</ol>
 
</ol>
 
</p>
 
</p>
Line 380: Line 388:
 
</a>
 
</a>
 
</div>
 
</div>
   
+
 
 
     </div>
 
     </div>
 
</div>
 
</div>
Line 486: Line 494:
 
   font-weight: 400;
 
   font-weight: 400;
 
   text-align: center;
 
   text-align: center;
 
+
 
 
}
 
}
 
.blocks {
 
.blocks {
Line 550: Line 558:
 
}
 
}
  
.border {  border-width: 2px;  border-color: #02263E;  border-style: solid; }  
+
.border {  border-width: 2px;  border-color: #02263E;  border-style: solid; }
  
 
.blocks-content__content {
 
.blocks-content__content {
Line 717: Line 725:
  
  
     <main>  
+
     <main>
 
<h1>Modelling</h1>
 
<h1>Modelling</h1>
  
Line 725: Line 733:
  
 
<p>A major problem the project faced is that the comparison process of the fluorescence proteins wouldn't be possible to be investigated with all combinations as it would take too long. </p>
 
<p>A major problem the project faced is that the comparison process of the fluorescence proteins wouldn't be possible to be investigated with all combinations as it would take too long. </p>
<p> &nbsp; </p>  
+
<p> &nbsp; </p>
 
<p> To answer this problem, the team will attempt to model the fluorescence spectra over time expressed by the proteins given different. First, the type of gene expression would need to be identified and then, would be modified to considered the effects of inhibition and finally, be applied over time to see how much expression would occur at a certain time period. The team will use Mathematical modeling such as Ordinary Differential Equations because they are easy to convert into programming in order to build components for the simulation.</p>
 
<p> To answer this problem, the team will attempt to model the fluorescence spectra over time expressed by the proteins given different. First, the type of gene expression would need to be identified and then, would be modified to considered the effects of inhibition and finally, be applied over time to see how much expression would occur at a certain time period. The team will use Mathematical modeling such as Ordinary Differential Equations because they are easy to convert into programming in order to build components for the simulation.</p>
 
<p> &nbsp; </p>
 
<p> &nbsp; </p>
 
<p> As a side project, the team will also investigate into whether our method is random and unique by investigating how many combinations we can make and whether we can accurately predict which combination will occur. </p>
 
<p> As a side project, the team will also investigate into whether our method is random and unique by investigating how many combinations we can make and whether we can accurately predict which combination will occur. </p>
<br> </br>  
+
<br> </br>
  
 
     </main>
 
     </main>
    <!-- Component starts here-->
+
 
 
     <ul class="blocks-names">
 
     <ul class="blocks-names">
 
       <li class="blocks__name">Constitutive Gene Expression </li>
 
       <li class="blocks__name">Constitutive Gene Expression </li>
Line 752: Line 760:
 
<center>
 
<center>
 
<br> </br> <br> </br> <br> </br> <br> </br> <br> </br> <br> </br> <br> </br> <br> </br><br> </br> <br> </br><br> </br> <br> </br> <br> </br> <br> </br> <br> </br> <br> </br> <br> </br> <br> </br> <br> </br> <br> </br><br> </br> <br> </br><br> </br> <br> </br> <br> </br> <br> </br> <br> </br> <br> </br> <br> </br> <br> </br> <br> </br> <br> </br><br> </br> <br> </br><br> </br><br> </br> <br> </br><br> </br> <br> </br><br> </br><br> </br> <br> </br><br> </br> <br> </br><br> </br>
 
<br> </br> <br> </br> <br> </br> <br> </br> <br> </br> <br> </br> <br> </br> <br> </br><br> </br> <br> </br><br> </br> <br> </br> <br> </br> <br> </br> <br> </br> <br> </br> <br> </br> <br> </br> <br> </br> <br> </br><br> </br> <br> </br><br> </br> <br> </br> <br> </br> <br> </br> <br> </br> <br> </br> <br> </br> <br> </br> <br> </br> <br> </br><br> </br> <br> </br><br> </br><br> </br> <br> </br><br> </br> <br> </br><br> </br><br> </br> <br> </br><br> </br> <br> </br><br> </br>
<br> </br> <br> </br> <br> </br> <br> </br> <br> </br> <br> </br> <br> </br>  
+
<br> </br> <br> </br> <br> </br> <br> </br> <br> </br> <br> </br> <br> </br>
  
 
<h2> Constitutive Gene Expression For Protein and mRNA Expression over Time </h2> <i class="blocks__content-close fa fa-times"></i>
 
<h2> Constitutive Gene Expression For Protein and mRNA Expression over Time </h2> <i class="blocks__content-close fa fa-times"></i>
Line 762: Line 770:
  
 
<sup> Figure 1 </sup>
 
<sup> Figure 1 </sup>
$$  sfGFP \underset{Transcriptin}{\rightarrow} mRNA \underset{Translation}{\rightarrow} sfGFP $$  
+
$$  sfGFP \underset{Transcriptin}{\rightarrow} mRNA \underset{Translation}{\rightarrow} sfGFP $$
<p> The equation above describes the process of which the gene undergoes transcription to produce mRNA. The mRNA carries the genetic information copied from the DNA which codes for protein. The expression of protein, can therefore, be measured by the fluorescence which is the desired output of the system. </p>  
+
<p> The equation above describes the process of which the gene undergoes transcription to produce mRNA. The mRNA carries the genetic information copied from the DNA which codes for protein. The expression of protein, can therefore, be measured by the fluorescence which is the desired output of the system. </p>
  
 
<sup> Figure 2 </sup>
 
<sup> Figure 2 </sup>
Line 774: Line 782:
 
<p>  We can apply Law of Mass Action combine both equations for the concentration of protein and mRNA over time. This model can be described as: </p>
 
<p>  We can apply Law of Mass Action combine both equations for the concentration of protein and mRNA over time. This model can be described as: </p>
  
<sup> Figure 3 </sup>  
+
<sup> Figure 3 </sup>
 
$$ mRNA = k_{1} -d _{1 } mRNA  $$
 
$$ mRNA = k_{1} -d _{1 } mRNA  $$
 
$$ Protein = k_{2} \cdot  mRNA - d_{2} \cdot Protein $$
 
$$ Protein = k_{2} \cdot  mRNA - d_{2} \cdot Protein $$
Line 783: Line 791:
 
   <li>Protein is the concentration of Protein</li>
 
   <li>Protein is the concentration of Protein</li>
 
   <li>k<sub> 1 </sub> is the constitutive transcription rate. This represents the number of mRNA molecules produced per gene, per unit of time.</li>
 
   <li>k<sub> 1 </sub> is the constitutive transcription rate. This represents the number of mRNA molecules produced per gene, per unit of time.</li>
<li> d <sub> 1 </sub>  is  the  mRNA  degradation  rate </li>  
+
<li> d <sub> 1 </sub>  is  the  mRNA  degradation  rate </li>
 
   <li>k<sub> 2 </sub> is the translation rate.  This represents the number of protein molecules produced per mRNA molecule, per unit of time.</li>
 
   <li>k<sub> 2 </sub> is the translation rate.  This represents the number of protein molecules produced per mRNA molecule, per unit of time.</li>
<li> d <sub> 2 </sub>  is  the  protein  degradation  rate. </li>  
+
<li> d <sub> 2 </sub>  is  the  protein  degradation  rate. </li>
 
</ul>
 
</ul>
<br> </br>  
+
<br> </br>
  
<p> This is important because we can use this model to calculate the concentration of proteins we can expect over time. This is useful as we can use this information to calculate the total emitted light spectra during the time period which is what we are looking for in our system. However, the constants and variables are individual for each protein and which means parameters for each protein would need to be found. These constants were found using literature <sup> 3 </sup> (for GFP) and lab results (the rest.) </p>  
+
<p> This is important because we can use this model to calculate the concentration of proteins we can expect over time. This is useful as we can use this information to calculate the total emitted light spectra during the time period which is what we are looking for in our system. However, the constants and variables are individual for each protein and which means parameters for each protein would need to be found. These constants were found using literature <sup> 3 </sup> (for GFP) and lab results (the rest.) </p>
  
 
<br> </br>
 
<br> </br>
Line 809: Line 817:
 
<br> </br>
 
<br> </br>
 
<br> </br>
 
<br> </br>
<br> </br> <br> </br> <br> </br> <br> </br> <br> </br> <br> </br> <br> </br> <br> </br><br> </br> <br> </br><br> </br> <br> </br> <br> </br> <br> </br> <br> </br> <br> </br> <br> </br> <br> </br> <br> </br> <br> </br><br> </br> <br> </br><br> </br> <br> </br> <br> </br><br> </br> <br> </br> <br> </br><br> </br> <br> </br> <br> </br><br> </br> <br> </br>  
+
<br> </br> <br> </br> <br> </br> <br> </br> <br> </br> <br> </br> <br> </br> <br> </br><br> </br> <br> </br><br> </br> <br> </br> <br> </br> <br> </br> <br> </br> <br> </br> <br> </br> <br> </br> <br> </br> <br> </br><br> </br> <br> </br><br> </br> <br> </br> <br> </br><br> </br> <br> </br> <br> </br><br> </br> <br> </br> <br> </br><br> </br> <br> </br>
 
<h2> Absorption and Emission Wavelengths From Given Concentrations of sfGFP, mRFP & ECFP </h2><i class="blocks__content-close fa fa-times"></i>
 
<h2> Absorption and Emission Wavelengths From Given Concentrations of sfGFP, mRFP & ECFP </h2><i class="blocks__content-close fa fa-times"></i>
 
<p> After concluding the general scheme we would be using, the team evaluated the selection of proteins. The proteins selected for the system use fluorescence, indicating they take in a light at a certain wavelength, and re-emit it at a different wavelength. This has to be considered because it informs the wet-lab in knowing which wavelengths are required to produce a spectra as well as highlighting the importance of considering any side effects from producing the spectra such as light being reabsorbed and re-emitted at a different wavelength / color, which would result in the spectra being similar to each other rather than unique. </p>
 
<p> After concluding the general scheme we would be using, the team evaluated the selection of proteins. The proteins selected for the system use fluorescence, indicating they take in a light at a certain wavelength, and re-emit it at a different wavelength. This has to be considered because it informs the wet-lab in knowing which wavelengths are required to produce a spectra as well as highlighting the importance of considering any side effects from producing the spectra such as light being reabsorbed and re-emitted at a different wavelength / color, which would result in the spectra being similar to each other rather than unique. </p>
 
<p> In order to save time and program a model, the team used Semrock's Online Fluorescence graph maker <sup> 1 </sup> which operated by taking in the expected Absorption wavelengths and emitting the Emission wavelengths expected by sfGFP (green), mRFP (red) and ECFP (blue) proteins. This was done through the Web App on the website. Furthermore, they provided the raw data in a text file format which was useful as it allows the team to read the data into a stand alone program. </p>
 
<p> In order to save time and program a model, the team used Semrock's Online Fluorescence graph maker <sup> 1 </sup> which operated by taking in the expected Absorption wavelengths and emitting the Emission wavelengths expected by sfGFP (green), mRFP (red) and ECFP (blue) proteins. This was done through the Web App on the website. Furthermore, they provided the raw data in a text file format which was useful as it allows the team to read the data into a stand alone program. </p>
  
<sup> Figure 4 </sup>  
+
<sup> Figure 4 </sup>
  
 
<img src="https://static.igem.org/mediawiki/2017/8/8f/T--UNOTT--SpectrumAbsoprtionEM.png" class="border" width="550" height="300" style= position: fixed; align=center;>
 
<img src="https://static.igem.org/mediawiki/2017/8/8f/T--UNOTT--SpectrumAbsoprtionEM.png" class="border" width="550" height="300" style= position: fixed; align=center;>
Line 833: Line 841:
 
         <div class="content__close"></div>
 
         <div class="content__close"></div>
 
<center>
 
<center>
<br> </br> <br> </br> <br> </br> <br> </br> <br> </br> <br> </br> <br> </br> <br> </br><br> </br> <br> </br><br> </br> <br> </br> <br> </br> <br> </br> <br> </br> <br> </br> <br> </br> <br> </br> <br> </br> <br> </br><br> </br> <br> </br><br> </br> <br> </br> <br> </br> <br> </br> <br> </br> <br> </br> <br> </br> <br> </br> <br> </br> <br> </br><br> </br> <br> </br><br> </br><br> </br>  
+
<br> </br> <br> </br> <br> </br> <br> </br> <br> </br> <br> </br> <br> </br> <br> </br><br> </br> <br> </br><br> </br> <br> </br> <br> </br> <br> </br> <br> </br> <br> </br> <br> </br> <br> </br> <br> </br> <br> </br><br> </br> <br> </br><br> </br> <br> </br> <br> </br> <br> </br> <br> </br> <br> </br> <br> </br> <br> </br> <br> </br> <br> </br><br> </br> <br> </br><br> </br><br> </br>
 
<h2>Gene Transcription Regulation by Repressors (CRISPRi) - Concentration over Time</h2> <i class="blocks__content-close fa fa-times"></i>
 
<h2>Gene Transcription Regulation by Repressors (CRISPRi) - Concentration over Time</h2> <i class="blocks__content-close fa fa-times"></i>
 
<p> The next step in developing our simulation was to calculate our protein concentration at any given time when using CRISPRi. Discussion with wet-lab revealed our method would be using CRISPRi as a repressor, which works by inhibiting the expression of one or more genes by binding to the promoter region <sup> 1 </sup>. The expanded mRNA and Protein concentration models from the Constitutive Gene Expression Model <sup> 2 </sup> were modified to include the element of repression from the CRISPRi inhibition. </p>
 
<p> The next step in developing our simulation was to calculate our protein concentration at any given time when using CRISPRi. Discussion with wet-lab revealed our method would be using CRISPRi as a repressor, which works by inhibiting the expression of one or more genes by binding to the promoter region <sup> 1 </sup>. The expanded mRNA and Protein concentration models from the Constitutive Gene Expression Model <sup> 2 </sup> were modified to include the element of repression from the CRISPRi inhibition. </p>
Line 843: Line 851:
 
$$  sfGFP \underset{Degradation}{\rightarrow} \oslash  $$
 
$$  sfGFP \underset{Degradation}{\rightarrow} \oslash  $$
  
<h5> This change can be applied to the Law of Mass Action <sup> 3 </sup> : </h5>  
+
<h5> This change can be applied to the Law of Mass Action <sup> 3 </sup> : </h5>
 
$$ m = k_{1} \cdot \frac{k^{n}}{k^{n} + R^{n}}- d_{1}m $$
 
$$ m = k_{1} \cdot \frac{k^{n}}{k^{n} + R^{n}}- d_{1}m $$
$$ p = k_{2} m - d_{2}p $$  
+
$$ p = k_{2} m - d_{2}p $$
  
 
<p>Where...</p>
 
<p>Where...</p>
Line 851: Line 859:
 
<p>m is mRNA concentration, p is Protein concentration, R is Repressor, k1 is Max Transcription Rate, k is the Repression Coefficient, n is number of repressors that need to cooperatively bind the promoter to trigger the inhibition of gene expression (Hill Coefficient), R is Repressor, d1 is mRNA degradation rate, d2 is Protein degradation rate </p>
 
<p>m is mRNA concentration, p is Protein concentration, R is Repressor, k1 is Max Transcription Rate, k is the Repression Coefficient, n is number of repressors that need to cooperatively bind the promoter to trigger the inhibition of gene expression (Hill Coefficient), R is Repressor, d1 is mRNA degradation rate, d2 is Protein degradation rate </p>
  
<p> The value for these constants and variables were taken from literature and calculating them <sup> 4 </sup> but later, adjusted to the lab results.</p>  
+
<p> The value for these constants and variables were taken from literature and calculating them <sup> 4 </sup> but later, adjusted to the lab results.</p>
  
 
<sup> Figure 6 </sup>
 
<sup> Figure 6 </sup>
<br> </br>
+
<br> </br>
 
<img src="https://static.igem.org/mediawiki/2017/7/73/T--UNOTT--InhibitedAndNon.png" class="border" height="350" width="550" style= position: fixed; align=center; >
 
<img src="https://static.igem.org/mediawiki/2017/7/73/T--UNOTT--InhibitedAndNon.png" class="border" height="350" width="550" style= position: fixed; align=center; >
  
<br> </br>  
+
<br> </br>
 
<p> Figure 6 shows the structure which underwent CRISPRi inhibition are expected to produce lower concentration of the protein whose expression were are inhibiting. This is important as it means the team can calculate concentration of proteins which are inhibited and compare them to the control conditions as well as giving the correct concentration for the simulation. </p>
 
<p> Figure 6 shows the structure which underwent CRISPRi inhibition are expected to produce lower concentration of the protein whose expression were are inhibiting. This is important as it means the team can calculate concentration of proteins which are inhibited and compare them to the control conditions as well as giving the correct concentration for the simulation. </p>
  
Line 869: Line 877:
  
 
       <li class="blocks-content__content">
 
       <li class="blocks-content__content">
     
+
 
 
         <ul class="blocks-content">
 
         <ul class="blocks-content">
 
         <div class="content__close"></div>
 
         <div class="content__close"></div>
 
<center>
 
<center>
<br> </br> <br> </br> <br>  
+
<br> </br> <br> </br> <br>
  
 
<h2>Relationship between Max Fluorescence and Protein Concentration</h2> <i class="blocks__content-close fa fa-times"></i>
 
<h2>Relationship between Max Fluorescence and Protein Concentration</h2> <i class="blocks__content-close fa fa-times"></i>
Line 880: Line 888:
 
<p> This means that an equation must be developed to find out what the intensity of fluorescence would be at that certain time. This consisted of of calculating the protein concentration at the time and using real life lab data of the fluorescence at that time period, the team could map that intensity to the protein concentration at that time. </p>
 
<p> This means that an equation must be developed to find out what the intensity of fluorescence would be at that certain time. This consisted of of calculating the protein concentration at the time and using real life lab data of the fluorescence at that time period, the team could map that intensity to the protein concentration at that time. </p>
  
<p> When the fluorescence data received from the wet lab were graphed, a model was constructed, refined and optimised to demonstrate the trends shown from the real data gained from the labs. Originally, the data from the lab was the Fluorescence against Time but by using the Gene Transcription Regulation by Repressors model developed earlier <sup> 1 </sup>, the team was able to estimate the protein concentration at a certain time periods.  </p>  
+
<p> When the fluorescence data received from the wet lab were graphed, a model was constructed, refined and optimised to demonstrate the trends shown from the real data gained from the labs. Originally, the data from the lab was the Fluorescence against Time but by using the Gene Transcription Regulation by Repressors model developed earlier <sup> 1 </sup>, the team was able to estimate the protein concentration at a certain time periods.  </p>
<sup> Figure 7 </sup>  
+
<sup> Figure 7 </sup>
<br> </br>  
+
<br> </br>
 
<img src="https://static.igem.org/mediawiki/2017/1/19/T--UNOTT--ProteinConcVsFluorescence.png" class="border" height="550" width="800" style= position: fixed; align=center; >
 
<img src="https://static.igem.org/mediawiki/2017/1/19/T--UNOTT--ProteinConcVsFluorescence.png" class="border" height="550" width="800" style= position: fixed; align=center; >
<br> </br>  
+
<br> </br>
<p> Due to time constraints, rather than implementing the relationship directly from lab data, the data was fitted using a Polynomial Fit of Order 3 using Excel and an equation was calculated from these. These equations were directly plugged into the simulation. However, this is inaccurate as the R squared value was ... , suggesting that it doesn't fully capture the data trend. </p>  
+
<p> Due to time constraints, rather than implementing the relationship directly from lab data, the data was fitted using a Polynomial Fit of Order 3 using Excel and an equation was calculated from these. These equations were directly plugged into the simulation. However, this is inaccurate as the R squared value was ... , suggesting that it doesn't fully capture the data trend. </p>
  
 
<p> These relationships were implemented into the simulation to give the expected spectra produced by each protein. This highlights another use: by adding or subtracting values from our fit, we can create a threshold for our Keys. This was essential when developing the Raw Data Simulator. <sup> 2 </sup></p>
 
<p> These relationships were implemented into the simulation to give the expected spectra produced by each protein. This highlights another use: by adding or subtracting values from our fit, we can create a threshold for our Keys. This was essential when developing the Raw Data Simulator. <sup> 2 </sup></p>
Line 896: Line 904:
 
       </li>
 
       </li>
  
     
+
 
 
       <li class="blocks-content__content">
 
       <li class="blocks-content__content">
     
+
 
 
         <ul class="blocks-content">
 
         <ul class="blocks-content">
 
         <div class="content__close"></div>
 
         <div class="content__close"></div>
Line 909: Line 917:
  
 
<p> However, in this proof of concept, order is irrelevant as the gene is either inhibited (1) or not (0). Using </p>
 
<p> However, in this proof of concept, order is irrelevant as the gene is either inhibited (1) or not (0). Using </p>
$$ n ^ r $$  
+
$$ n ^ r $$
 
<p> Where n = 2 and r = 3, this gives us a total combination of 2<sup> 3 </sup> {1,1,1} {1,1,0} {1,0,1} {1,0,0} {0,1,1} {0,1,0} {0,0,1} {0,0,0} </p>
 
<p> Where n = 2 and r = 3, this gives us a total combination of 2<sup> 3 </sup> {1,1,1} {1,1,0} {1,0,1} {1,0,0} {0,1,1} {0,1,0} {0,0,1} {0,0,0} </p>
  
<p> Randomness comes from the fact the system relies on Brownian Motion <sup> 1 </sup>, a random process to create these combinations.</p>  
+
<p> Randomness comes from the fact the system relies on Brownian Motion <sup> 1 </sup>, a random process to create these combinations.</p>
  
<p> However, in order for a movement to fall under Brownian Motion, it must fulfill a condition where the process must have continuous paths. This is not true as once the structures begin to form, the paths stop  (they do not collide off each other elastically, but rather, combine.) Furthermore, the bacterium might become biased towards options that put less metabolic stress on the bacterium, which results in selection. Alternatively, using metabolites to undergo transposition can improve randomness. <sup> 2 </sup> </p>  
+
<p> However, in order for a movement to fall under Brownian Motion, it must fulfill a condition where the process must have continuous paths. This is not true as once the structures begin to form, the paths stop  (they do not collide off each other elastically, but rather, combine.) Furthermore, the bacterium might become biased towards options that put less metabolic stress on the bacterium, which results in selection. Alternatively, using metabolites to undergo transposition can improve randomness. <sup> 2 </sup> </p>
  
<p> In order to aid, with the wet lab in what combinations they can expect, the team developed an Excel Spreadsheet where a user can simply input details of the construction and it would show what construction it would look like </p>  
+
<p> In order to aid, with the wet lab in what combinations they can expect, the team developed an Excel Spreadsheet where a user can simply input details of the construction and it would show what construction it would look like </p>
  
 
<p> Members of the public are encouraged to try it out and can use it to help with identifying how their spectra would look if they used the same proteins the project used </p>
 
<p> Members of the public are encouraged to try it out and can use it to help with identifying how their spectra would look if they used the same proteins the project used </p>
Line 922: Line 930:
 
<a href="https://github.com/BurgundyIsAPublicEnemy/iGEMNotts2017/tree/master/Models">Excel Spreadsheet</a>
 
<a href="https://github.com/BurgundyIsAPublicEnemy/iGEMNotts2017/tree/master/Models">Excel Spreadsheet</a>
  
<br> </br>  
+
<br> </br>
 
<p> <sup> 1 </sup> Refer to https://statistics.stanford.edu/sites/default/files/EFS%20NSF%20149.pdf </p>
 
<p> <sup> 1 </sup> Refer to https://statistics.stanford.edu/sites/default/files/EFS%20NSF%20149.pdf </p>
 
<p> <sup> 2 </sup> Refer https://link.springer.com/book/10.1007%2F978-1-4612-0459-6 for more information about Brownian Motion and Random Walk. </p>
 
<p> <sup> 2 </sup> Refer https://link.springer.com/book/10.1007%2F978-1-4612-0459-6 for more information about Brownian Motion and Random Walk. </p>

Revision as of 19:31, 31 October 2017





MODELLING

Overview

Constitutive Gene Expression For Protein and mRNA Expression over Time

Gene Transcription Regulation by Repressors (CRISPRi) - Concentration over Time

STEP 4: Random Ligations

STEP 5: Freeze Drying & Revival

STEP 6: CRISPRi & gRNA Efficiency








//----------------------------------------------------------------------------