Team:Tianjin/Modeling/Calculation

From 2012.igem.org

(Difference between revisions)
Line 1: Line 1:
{{:Team:Tianjin/frame/modeling}}
{{:Team:Tianjin/frame/modeling}}
 +
<html>
 +
<script type="text/javascript">
 +
$(document).ready(function()
 +
{
 +
//slides the element with class "menu_body" when mouse is over the paragraph
 +
$("#secondpane p.menu_head").mouseover(function()
 +
    {
 +
    $(this).next("div.menu_body").slideDown(500).siblings("div.menu_body").slideUp("slow");
 +
        $(this).siblings();
 +
});
 +
});
 +
</script>
 +
 +
<style type="text/css">
 +
#secondpane {color:#fff;
 +
width: 180px;
 +
margin-left:5px;
 +
}
 +
.menu_list {
 +
position:fixed;top:0;}
 +
.menu_head {
 +
padding: 5px 10px;
 +
cursor: pointer;
 +
position: relative;
 +
border:1px solid #fff;
 +
margin:1px;
 +
    font-weight:bold;
 +
    background: #ffbe10;
 +
color:#fff
 +
}
 +
.menu_head:hover
 +
{
 +
cursor: pointer;
 +
color: #fff;
 +
background: #d39f10;
 +
text-decoration: none;
 +
background-position: 0 0;
 +
}
 +
.menu_body {
 +
display:none;
 +
}
 +
.menu_children {
 +
background: #aaa;
 +
cursor: pointer;
 +
padding: 5px 10px;
 +
border-top: 1px solid #fff;
 +
border-bottom: 1px solid #fff;
 +
border-left: 1px solid #fff;
 +
border-right: 1px solid #fff;
 +
}
 +
.menu_body a{
 +
  display:block;
 +
  color:#fff;
 +
  padding-left:10px;
 +
  font-weight:bold;
 +
  text-decoration:none;
 +
}
 +
.menu_children:hover{
 +
cursor: pointer;
 +
  color: #fff;
 +
  text-decoration:none;
 +
  background:#777777;
 +
  }
 +
 +
</style>
 +
<div style="float:left;">
 +
  <div id="secondpane"> <!--Code for menu starts here-->
 +
<p class="menu_head">Header-1</p>
 +
<div class="menu_body">
 +
        <div class="menu_children">
 +
<a href="#">Link-1</a>
 +
        </div>
 +
        <div class="menu_children">
 +
<a href="#">Link-2</a>
 +
        </div>
 +
        <div class="menu_children">
 +
<a href="#">Link-3</a>
 +
        </div>
 +
</div>
 +
<p class="menu_head">Header-2</p>
 +
<div class="menu_body">
 +
        <div class="menu_children">
 +
<a href="#">Link-1</a>
 +
        </div>
 +
        <div class="menu_children">
 +
<a href="#">Link-2</a>
 +
        </div>
 +
        <div class="menu_children">
 +
<a href="#">Link-3</a>
 +
        </div>
 +
</div>
 +
<p class="menu_head">Header-3</p>
 +
<div class="menu_body">
 +
        <div class="menu_children">
 +
<a href="#">Link-1</a>
 +
        </div>
 +
        <div class="menu_children">
 +
<a href="#">Link-2</a>
 +
        </div>
 +
        <div class="menu_children">
 +
<a href="#">Link-3</a>
 +
        </div>
 +
      </div>
 +
  </div>      <!--Code for menu ends here-->
 +
</div>
 +
</html>
 +
<div id="text-content">
 +
<br><center><span style="font-size:46px;font-family:Cambria;margin-top:10px;line-height:80%">Calculation and Derivation of the protein expression amount model in three states</span></center>
<br><center><span style="font-size:46px;font-family:Cambria;margin-top:10px;line-height:80%">Calculation and Derivation of the protein expression amount model in three states</span></center>

Revision as of 02:32, 23 September 2012



Calculation and Derivation of the protein expression amount model in three states


Contents

Overview

Three problems came up while we start calculating: what sequences to input, which method to use and how cogent the result will be. As for sequence, we input the both the SD and the protein coding sequence.

Our goal is to calculate the total ΔG of each reaction and then predict the amount of protein expressed. There are several softwares dealing with DNA or RNA base-pairing progress, such as NUPACK, RBS-Calculator, and Vienna RNA etc. After comparison, we decide to use RBS-Calculator.

Due to the complication of translation progress and our lack of insight in this issue, the results of our modeling can’t be very precise. But at least it should have the precision of order of magnitude.

Input of Our Calculation

RBS sequence

  • RFP, normal RBS

ATTTCACACATACTAGAGAAAGAGGAGAAATACTAGATGGCTTCCTCCGAAGACGTTATCAAAGAGTT

CATGCGTT

  • RFP, orthogonal RBS

ATTTCACACATGTTCCGTACTAGATGGCTTCCTCCGAAGACGTTATCAAAGAGTTCATGCGTT

  • GFP, normal RBS

TACTAGAGAAAGAGGAGAAATACTAGATGCGTAAAGGAGAAGAACTTTTCACTGGAGTTGTCCCAAT

TCTTGTT

16S rRNA sequence

  • normal 16S: ACCTCCTTA
  • orthogonal 16S: ACGGAACTA

Formula Derivation

Basic Idea

For model design, please refer to Design part. Data of the curve Er-time is obtained from Ec=K∙Er experiments. The function of our model is to work out the proportion factor 'K'.

Basic Assumption

  1. The expression of RFP and GFP are independent.
  2. The expression of the two proteins are determined by the percentage of normal and orthogonal ribosomes rather than the number of the two ribosomes.
  3. The growth curve of bacteria do not change significantly after the transferred into orthogonal protein expression system.

Formula Derivation

equation1

We start from Formula 1 and 2. In Formula 1, m stands for the number of mRNA transcript, Rtot is the total number of ribosomes, β is the apparent Boltzmann constant, ∆Gtot is the total change of Gibbs free energy, k is proportion factor.

Because the GFP and RFP coding sequence are on the same mRNA transcript, the values of m of GFP and RFP are always same. In different state and time, the total number of ribosomes varies. We assume that at the same time Rtot of different state remain same. Calculation of DGtot is the main job of this model. The proportion factor 'k' represents all unknown factors. Here we assume that as for the same protein in deferent state 'k' varies little.