Team:Utah State/Modeling

From 2009.igem.org

(Difference between revisions)
 
(12 intermediate revisions not shown)
Line 8: Line 8:
<!DOCTYPE html PUBLIC "-//W3C//DTD XHTML 1.0 Transitional//EN" "http://www.w3.org/TR/xhtml1/DTD/xhtml1-transitional.dtd">
<!DOCTYPE html PUBLIC "-//W3C//DTD XHTML 1.0 Transitional//EN" "http://www.w3.org/TR/xhtml1/DTD/xhtml1-transitional.dtd">
<html xmlns="http://www.w3.org/1999/xhtml">
<html xmlns="http://www.w3.org/1999/xhtml">
 +
<head>
<head>
<meta http-equiv="Content-Type" content="text/html; charset=utf-8" />
<meta http-equiv="Content-Type" content="text/html; charset=utf-8" />
Line 84: Line 85:
         <td id="nav"><a href="https://2009.igem.org/Team:Utah_State/Project"><font size = 4>PROJECT</font></a></td>
         <td id="nav"><a href="https://2009.igem.org/Team:Utah_State/Project"><font size = 4>PROJECT</font></a></td>
       </tr>
       </tr>
-
      <tr>
+
        <tr>
         <td id="nav"><a href="https://2009.igem.org/Team:Utah_State/Parts"><font size = 4>BIOBRICKS</font></a></td>
         <td id="nav"><a href="https://2009.igem.org/Team:Utah_State/Parts"><font size = 4>BIOBRICKS</font></a></td>
       </tr>
       </tr>
       <tr>
       <tr>
       <td width="172" id="ana"><span class="currentPage"><font size = 4>MODELING</font></span>
       <td width="172" id="ana"><span class="currentPage"><font size = 4>MODELING</font></span>
 +
        <a href="#parameters">Parameters</a><br />
 +
        <a href="#simulations">Simulations</a><br />
 +
        <a href="#references">References</a>
       </tr>
       </tr>
       <tr>
       <tr>
Line 97: Line 101:
       </tr>
       </tr>
       <tr>
       <tr>
-
       <td id="nav"><a href="https://2009.igem.org/Team:Utah_State/Achievements"><font size = 4>SUCCESS</font></a></td>
+
       <td id="nav"><a href="https://2009.igem.org/Team:Utah_State/Achievements"><font size = 4>JUDGING</font></a></td>
 +
      </tr>
 +
    <tr>
 +
      <td id="nav"><a href="https://2009.igem.org/Team:Utah_State/Contact"><font size = 4>CONTACT</font></a></td>
       </tr>
       </tr>
     </table></td>
     </table></td>
Line 111: Line 118:
   </font>
   </font>
   <p class = "class">  
   <p class = "class">  
-
There is currently limited information on how to properly select which secretion pathway should be used for a given protein (Choi, 2004). The construction of models that evaluate protein secretion may provide a useful framework for studying and attempting to predict secretion efficiency. Models of the Sec and Tat translocation pathways were made using MATLAB’s Simbiology toolbox. They were evaluated using the embedded ordinary differential equation solver in the software. Assumptions were made for both models due to lack of time scale information for individual steps for these translocation mechanisms. However, both provide a flexible framework that can become more detailed as additional information is found in literature or discovered in the laboratory. Within both models, a protein is first generated and then carried out of the cytoplasm to the periplasm. Only protein species were tracked, as species involved with making the protein are currently not as useful to monitor.
+
There are currently few guidelines for selecting which signal peptide (which determine the secretion pathway taken) should be used for a given recombinant protein (Choi, 2004). The construction of models that evaluate protein secretion may provide a useful framework for studying and attempting to predict factors that affect secretion efficiency. Models of the Sec and Tat translocation pathways were made using MATLAB’s Simbiology toolbox. They were evaluated using the embedded ordinary differential equation solver in the software. Assumptions were made for both models due to lack of time scale information for individual steps for these translocation mechanisms. However, both provide a flexible framework that can become more detailed as additional information is found in literature or discovered in the laboratory. Within both models, a protein is first generated and then carried out of the cytoplasm to the periplasm. Only protein species were tracked, as species involved with making the protein are currently not as useful to monitor.
-
   </P>
+
   </P> </td></tr>
 +
<tr>
 +
    <td>
 +
      <a name="parameters">
 +
    </td>
 +
    </tr>
 +
    <tr>
 +
    <td>
 +
    <b><i><font size="3" face="Helvetica, Arial, San Serif" color =#000033>
 +
Model Parameters
 +
    </font></b></i> 
  <p class = "class">  
  <p class = "class">  
In both models, parameters for the manufacture of protein were the same and found using averaged values for <i>E. coli</i>. The average length of an mRNA for <i>E. coli</i> is 1100 nucleotides and an <i>E. coli</i> cell transcribes at an average rate of 70 nucleotides per second. A simple rate calculation determines that a strand of mRNA is made approximately every 15.7 seconds. This value was then used to yield a first order reaction rate of (1/15.7) 1/second.  The translation rate was found in a similar manner. Given an average protein size in <i>E. coli</i> as 360 amino acids and the average translation rate is 40 amino acids per second, the first order translation rate constant is (1/9) 1/s (Institute for Biomolecular Design, 2008). To initiate the model, the concentration of the gene was set at 1, which was locked at that value over the course of the simulation. The median half life for mRNA was found to be 3.7 minutes, after which it will degrade (Milo, n. d.). For the Tat-dependent mechanism, the protein is fully folded prior to translocation (Mergulhão, 2005). For this folding time, a place holder value of 10 minutes was used as an approximate order of magnitude, although the value can easily be changed depending on parameters of the protein in question. Degradation of protein was modeled as the average time it takes the cell to reproduce, which is given as half an hour (Institute for Biomolecular Design, 2008).
In both models, parameters for the manufacture of protein were the same and found using averaged values for <i>E. coli</i>. The average length of an mRNA for <i>E. coli</i> is 1100 nucleotides and an <i>E. coli</i> cell transcribes at an average rate of 70 nucleotides per second. A simple rate calculation determines that a strand of mRNA is made approximately every 15.7 seconds. This value was then used to yield a first order reaction rate of (1/15.7) 1/second.  The translation rate was found in a similar manner. Given an average protein size in <i>E. coli</i> as 360 amino acids and the average translation rate is 40 amino acids per second, the first order translation rate constant is (1/9) 1/s (Institute for Biomolecular Design, 2008). To initiate the model, the concentration of the gene was set at 1, which was locked at that value over the course of the simulation. The median half life for mRNA was found to be 3.7 minutes, after which it will degrade (Milo, n. d.). For the Tat-dependent mechanism, the protein is fully folded prior to translocation (Mergulhão, 2005). For this folding time, a place holder value of 10 minutes was used as an approximate order of magnitude, although the value can easily be changed depending on parameters of the protein in question. Degradation of protein was modeled as the average time it takes the cell to reproduce, which is given as half an hour (Institute for Biomolecular Design, 2008).
Line 119: Line 136:
In addition to the discussed values regarding the manufacture of protein, parameters for the Sec and Tat pathway are required. A literature review yielded no specific rate constants or time parameters for the individual steps involved with either process. However, one article stated that the Sec and Tat pathways take a few seconds and a few minutes, respectively, to translocate protein to the periplasmic space (Mergulhão, 2005). Accordingly, the Sec and Tat pathways were modeled to require 3 seconds and 3 minutes for protein translocation, respectively. First order rate constants were determined by taking the inverse of these values in seconds. Images of the model diagrams are seen below. The different pathways after the formation of protein can be activated and deactivated as more detailed information comes available about the translocation and secretion processes. The Sec secretion pathway also includes the SRP secretion pathway as a possible option. The constructed models for the Sec and Tat pathways are given as Figures 1 and 2, respectively.
In addition to the discussed values regarding the manufacture of protein, parameters for the Sec and Tat pathway are required. A literature review yielded no specific rate constants or time parameters for the individual steps involved with either process. However, one article stated that the Sec and Tat pathways take a few seconds and a few minutes, respectively, to translocate protein to the periplasmic space (Mergulhão, 2005). Accordingly, the Sec and Tat pathways were modeled to require 3 seconds and 3 minutes for protein translocation, respectively. First order rate constants were determined by taking the inverse of these values in seconds. Images of the model diagrams are seen below. The different pathways after the formation of protein can be activated and deactivated as more detailed information comes available about the translocation and secretion processes. The Sec secretion pathway also includes the SRP secretion pathway as a possible option. The constructed models for the Sec and Tat pathways are given as Figures 1 and 2, respectively.
</p>
</p>
-
<b><i> PUT IN FIGURES</b></i>
+
 
 +
<!---FIgure1--->
 +
<div align="center"><img src="https://static.igem.org/mediawiki/2009/2/2e/Model1.png"  align = "middle" height="230" style="padding:.5px; border-style:solid; border-color:#999" alt="Advisors"> </div>
 +
<br>
 +
<div align="center"><font size="2.5" face="Helvetica, Arial, San Serif" color =#231f20>
 +
<b>Figure 1.</b> Sec Pathway Model Diagram
 +
</div>
 +
<br>
 +
<br>
 +
<div align="center"><img src="https://static.igem.org/mediawiki/2009/b/b8/Model2.png"  align = "middle" height="199" style="padding:.5px; border-style:solid; border-color:#999" alt="Advisors"> </div>
 +
<br>
 +
<div align="center"><font size="2.5" face="Helvetica, Arial, San Serif" color =#231f20>
 +
<b>Figure 2.</b>  Tat Pathway Model Diagram
 +
</div>
 +
<br>
 +
<br>
 +
 
<p class = "class">  
<p class = "class">  
-
Simulations for the Sec pathways were run for 5400 seconds (1.5 hours), which allowed the number of periplasmic to reach approximately steady state. The simulation for the Sec pathway secretion model is shown below as Figure 3.
+
</td></tr>
-
</p>
+
<tr>
-
<b><i> PUT IN FIGURES </b></i>
+
    <td>
-
<p class = "class">
+
      <a name="simulations">
 +
    </td>
 +
    </tr>
 +
    <tr>
 +
    <td>
 +
    <b><i><font size="3" face="Helvetica, Arial, San Serif" color =#000033>
 +
Simulations
 +
    </font></b></i>
 +
<p class = "class">Simulations for the Sec pathways were run for 5400 seconds (1.5 hours), which allowed the number of periplasmic to reach approximately steady state. The simulation for the Sec pathway secretion model is shown below as Figure 3.  
Protein generated by the <i>E. coli</i> cell reaches steady state quickly, as it is both produced and translocated in to the periplasm rapidly. The periplasmic protein reaches approximately 90 proteins upon achieving steady state. The long settling time and high number of periplasmic protein are a result of the long protein degradation time.</p>
Protein generated by the <i>E. coli</i> cell reaches steady state quickly, as it is both produced and translocated in to the periplasm rapidly. The periplasmic protein reaches approximately 90 proteins upon achieving steady state. The long settling time and high number of periplasmic protein are a result of the long protein degradation time.</p>
 +
 +
<br>
 +
<div align="center"><img src="https://static.igem.org/mediawiki/2009/a/a3/ModifiedGraph1.jpg"  align = "middle" height="400" style="padding:.5px; border-style:solid; border-color:#999" alt="Figure 3"> </div>
 +
<br>
 +
<div align="center"><font size="2.5" face="Helvetica, Arial, San Serif" color =#231f20>
 +
<b>Figure 3.</b>  Simulation of protein production and translocation by the Sec pathway
 +
</div>
 +
<br>
 +
<p class= "class">Simulations for the Tat pathway were run for 3600 seconds (1 hour).  The shorter length of simulation relative to the Sec model is used despite having significantly longer protein translocation time because the time-consuming folding step limits the reaction. The results of the ordinary differential equation simulation are shown below in Figure 4. The amount of periplasmic protein is lower than cytoplasmic and folded protein values upon reaching steady state. This is due to the rate-limiting folding step, which also leads to an expected high concentration of folded protein as well.
<p class= "class">Simulations for the Tat pathway were run for 3600 seconds (1 hour).  The shorter length of simulation relative to the Sec model is used despite having significantly longer protein translocation time because the time-consuming folding step limits the reaction. The results of the ordinary differential equation simulation are shown below in Figure 4. The amount of periplasmic protein is lower than cytoplasmic and folded protein values upon reaching steady state. This is due to the rate-limiting folding step, which also leads to an expected high concentration of folded protein as well.
</p>
</p>
 +
 +
<br>
 +
<div align="center"><img src="https://static.igem.org/mediawiki/2009/d/d2/ModifiedGraph2.jpg"  align = "middle" height="400" style="padding:.5px; border-style:solid; border-color:#999" alt="Figure 3"> </div>
 +
<br>
 +
<div align="center"><font size="2.5" face="Helvetica, Arial, San Serif" color =#231f20>
 +
<b>Figure 3.</b>  Simulation of protein production and translocation by the Sec pathway
 +
</div>
 +
<br>
 +
<br>
 +
 +
<p class = "class"> The accuracy of models, such as those constructed, will increase as more knowledge about the time step parameters for each of the mechanisms is determined. These improvements will allow for a better comparison of the two secretion methods, which can lead to better signal peptide selection (and corresponding pathway selection) for the protein in question. Future work could also include incorporating promoter strength into the model to affect transcription rates to provide a more accurate depiction of what is happening inside the cell. Better correlation between these parameters and the model will result in a more useful aid in when studying protein secretion.  
<p class = "class"> The accuracy of models, such as those constructed, will increase as more knowledge about the time step parameters for each of the mechanisms is determined. These improvements will allow for a better comparison of the two secretion methods, which can lead to better signal peptide selection (and corresponding pathway selection) for the protein in question. Future work could also include incorporating promoter strength into the model to affect transcription rates to provide a more accurate depiction of what is happening inside the cell. Better correlation between these parameters and the model will result in a more useful aid in when studying protein secretion.  
</p>
</p>
Line 135: Line 196:
     </tr>
     </tr>
   </tr>
   </tr>
 +
 +
<!---Go To Top of Page--->
 +
 +
  <tr>
 +
    <td>
 +
      <font size="3" face="Century Gothic, Arial, San Serif" color =#000033 <div align="center">
 +
  <b><a href="#top">Go to Top of Page</a></b>
 +
</div> </font>
 +
    </td>
 +
    </tr>
 +
</table>
</table>
Line 142: Line 214:
   <tr>
   <tr>
     <td>
     <td>
-
     <p><b><i><font size="3" face="Helvetica, Arial, San Serif" color =#000033>
+
</td></tr>
 +
<tr>
 +
     <td>
 +
      <a name="references">
 +
    </td>
 +
    </tr>
 +
    <tr>
 +
    <td>
 +
    <b><i><font size="3" face="Helvetica, Arial, San Serif" color =#000033>
References
References
-
     </font></p><hr></b></i>
+
     </font></b></i> <hr>
 +
   
     <p><font size="2.5" face="Tahoma, Helvetica, Sans-Serif" color =#000000>
     <p><font size="2.5" face="Tahoma, Helvetica, Sans-Serif" color =#000000>
<ul class = "circle">
<ul class = "circle">
-
<li>Institute for Biomolecular Design (2008, August 1) CyberCell Database. Retrieved October 2009.
+
<li>Institute for Biomolecular Design (2008, August 1) CyberCell Database. Retrieved October 2009. http://redpoll.pharmacy.ualberta.ca/CCDB/cgi-bin/STAT_NEW.cgi
-
<http://redpoll.pharmacy.ualberta.ca/CCDB/cgi-bin/STAT_NEW.cgi>
+
<li>Choi JH, Lee SY (2004) Secretory and extracellular production of recombinant proteins using <i>Escherichia coli</i>. Appl Microbiol Biotechnol 64:625-635</li>
<li>Choi JH, Lee SY (2004) Secretory and extracellular production of recombinant proteins using <i>Escherichia coli</i>. Appl Microbiol Biotechnol 64:625-635</li>
<li>Mergulhão FJM, Summers DK, Montier GA (2005) Recombinant protein secretion in <i>Escherichia coli</i>. Biotechnol Adv 23:177-202</li>
<li>Mergulhão FJM, Summers DK, Montier GA (2005) Recombinant protein secretion in <i>Escherichia coli</i>. Biotechnol Adv 23:177-202</li>
-
<li>Milo R (n.d.) BioNumbers: The Database of Useful Biological Numbers. Retrieved October 2009. <http://bionumbers.hms.harvard.edu/default.aspx></li>
+
<li>Milo R (n.d.) BioNumbers: The Database of Useful Biological Numbers. Retrieved October 2009. http://bionumbers.hms.harvard.edu/default.aspx</li>
</ul>
</ul>
Line 243: Line 323:
font size="2.5";
font size="2.5";
font-style: italic;
font-style: italic;
 +
}
 +
p.margin2
 +
{
 +
margin-top:0px;
 +
margin-bottom:0px;
 +
margin-right:50px;
 +
margin-left:250px;
 +
color: #231f20;
 +
font-weight : none;
 +
font-family:Tahoma, Arial, Sans-serif;
 +
font size="2.5";
}
}
ul.circle {list-style-type:circle; color: #231f20;}
ul.circle {list-style-type:circle; color: #231f20;}
</style>
</style>
</head>
</head>

Latest revision as of 03:06, 22 October 2009

USU iGem Untitled Document

MODELING Parameters
Simulations
References
Modeling Secretion Mechanisms

There are currently few guidelines for selecting which signal peptide (which determine the secretion pathway taken) should be used for a given recombinant protein (Choi, 2004). The construction of models that evaluate protein secretion may provide a useful framework for studying and attempting to predict factors that affect secretion efficiency. Models of the Sec and Tat translocation pathways were made using MATLAB’s Simbiology toolbox. They were evaluated using the embedded ordinary differential equation solver in the software. Assumptions were made for both models due to lack of time scale information for individual steps for these translocation mechanisms. However, both provide a flexible framework that can become more detailed as additional information is found in literature or discovered in the laboratory. Within both models, a protein is first generated and then carried out of the cytoplasm to the periplasm. Only protein species were tracked, as species involved with making the protein are currently not as useful to monitor.

Model Parameters

In both models, parameters for the manufacture of protein were the same and found using averaged values for E. coli. The average length of an mRNA for E. coli is 1100 nucleotides and an E. coli cell transcribes at an average rate of 70 nucleotides per second. A simple rate calculation determines that a strand of mRNA is made approximately every 15.7 seconds. This value was then used to yield a first order reaction rate of (1/15.7) 1/second. The translation rate was found in a similar manner. Given an average protein size in E. coli as 360 amino acids and the average translation rate is 40 amino acids per second, the first order translation rate constant is (1/9) 1/s (Institute for Biomolecular Design, 2008). To initiate the model, the concentration of the gene was set at 1, which was locked at that value over the course of the simulation. The median half life for mRNA was found to be 3.7 minutes, after which it will degrade (Milo, n. d.). For the Tat-dependent mechanism, the protein is fully folded prior to translocation (Mergulhão, 2005). For this folding time, a place holder value of 10 minutes was used as an approximate order of magnitude, although the value can easily be changed depending on parameters of the protein in question. Degradation of protein was modeled as the average time it takes the cell to reproduce, which is given as half an hour (Institute for Biomolecular Design, 2008).

In addition to the discussed values regarding the manufacture of protein, parameters for the Sec and Tat pathway are required. A literature review yielded no specific rate constants or time parameters for the individual steps involved with either process. However, one article stated that the Sec and Tat pathways take a few seconds and a few minutes, respectively, to translocate protein to the periplasmic space (Mergulhão, 2005). Accordingly, the Sec and Tat pathways were modeled to require 3 seconds and 3 minutes for protein translocation, respectively. First order rate constants were determined by taking the inverse of these values in seconds. Images of the model diagrams are seen below. The different pathways after the formation of protein can be activated and deactivated as more detailed information comes available about the translocation and secretion processes. The Sec secretion pathway also includes the SRP secretion pathway as a possible option. The constructed models for the Sec and Tat pathways are given as Figures 1 and 2, respectively.

Advisors

Figure 1. Sec Pathway Model Diagram


Advisors

Figure 2. Tat Pathway Model Diagram


Simulations

Simulations for the Sec pathways were run for 5400 seconds (1.5 hours), which allowed the number of periplasmic to reach approximately steady state. The simulation for the Sec pathway secretion model is shown below as Figure 3. Protein generated by the E. coli cell reaches steady state quickly, as it is both produced and translocated in to the periplasm rapidly. The periplasmic protein reaches approximately 90 proteins upon achieving steady state. The long settling time and high number of periplasmic protein are a result of the long protein degradation time.


Figure 3

Figure 3. Simulation of protein production and translocation by the Sec pathway

Simulations for the Tat pathway were run for 3600 seconds (1 hour). The shorter length of simulation relative to the Sec model is used despite having significantly longer protein translocation time because the time-consuming folding step limits the reaction. The results of the ordinary differential equation simulation are shown below in Figure 4. The amount of periplasmic protein is lower than cytoplasmic and folded protein values upon reaching steady state. This is due to the rate-limiting folding step, which also leads to an expected high concentration of folded protein as well.


Figure 3

Figure 3. Simulation of protein production and translocation by the Sec pathway


The accuracy of models, such as those constructed, will increase as more knowledge about the time step parameters for each of the mechanisms is determined. These improvements will allow for a better comparison of the two secretion methods, which can lead to better signal peptide selection (and corresponding pathway selection) for the protein in question. Future work could also include incorporating promoter strength into the model to affect transcription rates to provide a more accurate depiction of what is happening inside the cell. Better correlation between these parameters and the model will result in a more useful aid in when studying protein secretion.

Go to Top of Page

References

  • Institute for Biomolecular Design (2008, August 1) CyberCell Database. Retrieved October 2009. http://redpoll.pharmacy.ualberta.ca/CCDB/cgi-bin/STAT_NEW.cgi
  • Choi JH, Lee SY (2004) Secretory and extracellular production of recombinant proteins using Escherichia coli. Appl Microbiol Biotechnol 64:625-635
  • Mergulhão FJM, Summers DK, Montier GA (2005) Recombinant protein secretion in Escherichia coli. Biotechnol Adv 23:177-202
  • Milo R (n.d.) BioNumbers: The Database of Useful Biological Numbers. Retrieved October 2009. http://bionumbers.hms.harvard.edu/default.aspx