Difference between revisions of "Team:HokkaidoU Japan/Model"

(Prototype team page)
 
 
(43 intermediate revisions by 3 users not shown)
Line 1: Line 1:
 
{{HokkaidoU_Japan}}
 
{{HokkaidoU_Japan}}
 +
{{Team:HokkaidoU_Japan/CSS}}
 +
{{HokkaidoU_Japan/footer}}
 +
 
<html>
 
<html>
  
 +
<br>
 +
<div id="Modeling"><img src="https://static.igem.org/mediawiki/2016/f/fd/T--HokkaidoU_Japan--modeling.png"
 +
width="300px" height="100px" alt="Modeling"></div>
  
<div class="column full_size judges-will-not-evaluate">
 
<h3>★  ALERT! </h3>
 
<p>This page is used by the judges to evaluate your team for the <a href="https://2016.igem.org/Judging/Awards#SpecialPrizes">Best Model award</a>. </p>
 
  
  
<p> Delete this box in order to be evaluated for this medal. See more information at <a href="https://2016.igem.org/Judging/Pages_for_Awards/Instructions"> Instructions for Pages for awards</a>.</p>
+
<span class="nomal2">
</div>
+
  
 +
<br>In our project, we tried to stabilize protein structure
 +
by circularization using SAR. About the detailed concept of
 +
circularization, please read <a href="https://2016.igem.org/Team:HokkaidoU_Japan/Circularization">circularization page</a>. To think about the power
 +
of stabilization by
 +
circularization using SAP, we used HP (Hydrophobic-Polar) model.
 +
HP model is a kind of simplified protein folding model and in this model,
 +
protein chain is given as zig-zag stick on 2D lattice. Each residue has
 +
the characteristic H or P (Hydrophobic or Polar). To calculate the stability
 +
of protein structures, in this model, if an H residue is next to another H
 +
residue without covalent bond, it decreases free energy because of hydrophobic
 +
interaction. In our model, the decreased energy by each hydrophobic interaction
 +
is defined as -E<span class="sitatuki">H</span>. We added another
 +
characteristic SAR into this model. Through thinking about this model,
 +
we can simply think about the effect of SAR reflected as the effect to
 +
probability to fold as native state. We thought SAR interaction is so strong,
 +
so in the case we add SAR at N terminus and C terminus, both ends are set
 +
next to each other in the model. So, let's think about the simplest model.
  
 +
<br>The simplest model is the model with the number of residue is
 +
4 and the sequence is HPPH. In this case, without SAR, the number
 +
of states is 4, excluding enantiomers and rotamers. The possible states
 +
and the energy are listed below. K<span class="sitatuki">B</span> is Boltzmann constant 1.38064852*10<sup>-23</sup> (J/K), T is temperature (K).
  
  
  
 +
<br>
 +
<center>
 +
<table style="border-style: none">
 +
<tr align="center">
 +
<td style="border-style: none;">
 +
      <tr>
 +
<td style="border-style:none; float:center">
 +
<img src="https://static.igem.org/mediawiki/2016/d/da/T--HokkaidoU_Japan--Basicpart_Model_fig1.png" width="600px" height="auto" alt="Fig_1"></td>
 +
      </tr>
 +
      <tr>
 +
<td style="border-style: none;" align="center"><span class="small">Fig. 1
  
 +
</span></td>
 +
      </tr>
 +
</table>
 +
</center>
 +
<br clear="all">
  
<div class="column full_size">
 
<h2> Modeling</h2>
 
<p>Mathematical models and computer simulations provide a great way to describe the function and operation of BioBrick Parts and Devices. Synthetic Biology is an engineering discipline, and part of engineering is simulation and modeling to determine the behavior of your design before you build it. Designing and simulating can be iterated many times in a computer before moving to the lab. This award is for teams who build a model of their system and use it to inform system design or simulate expected behavior in conjunction with experiments in the wetlab.</p>
 
  
  
<h5> Inspiration </h5>
+
<br>Only the first one is stable and its energy is -E<span class="sitatuki">H</span>.
<p>
+
Because it's most stable, we thought it's native state.
Here are a few examples from previous teams:
+
The probability to fold as native state (P<sup>wild</sup><span class="sitatuki">HPPH</span>) is below.
</p>
+
<ul>
+
<li><a href="https://2014.igem.org/Team:ETH_Zurich/modeling/overview">ETH Zurich 2014</a></li>
+
<li><a href="https://2014.igem.org/Team:Waterloo/Math_Book">Waterloo 2014</a></li>
+
</ul>
+
  
 +
<br>
 +
<center>
 +
<table style="border-style: none">
 +
<tr align="center">
 +
<td style="border-style: none;">
 +
      <tr>
 +
<td style="border-style:none; float:center">
 +
<img src="https://static.igem.org/mediawiki/2016/6/6b/T--HokkaidoU_Japan--siki_1.png" width="400px" height="auto" alt="Formula_1"></td>
 +
      </tr>
 +
      <tr>
 +
<td style="border-style: none"; align="center"><h2>
  
 +
</h2></td>
 +
      </tr>
 +
</table>
 +
</center>
 +
<br clear="all">
 +
 +
 +
 +
<br>To calculate this, we used canonical ensemble from statistical mechanics. The probability of causing state <span class="italic">i</span> is calculated through the function below.
 +
E<span class="sitatuki">i</SARn> is the energy of state <span class="italic">i</Span>, E<SARn class="sitatuki">j</span> is the energy of the state <span class="italic">j</span>.
 +
 +
 +
 +
<br>
 +
<center>
 +
<table style="border-style: none">
 +
<tr align="center">
 +
<td style="border-style: none;">
 +
      <tr>
 +
<td style="border-style:none; float:center">
 +
<img src="https://static.igem.org/mediawiki/2016/1/10/T--HokkaidoU_Japan--siki1.jpg" width="300px" height="auto" alt="Formula_2"></td>
 +
      </tr>
 +
      <tr>
 +
<td style="border-style: none"; align="center"><h2>
 +
 +
</h2></td>
 +
      </tr>
 +
</table>
 +
</center>
 +
<br clear="all">
 +
 +
 +
 +
<br>But with SAR, the number of states is 1 and the state is the most stable one.
 +
 +
 +
<br>
 +
<center>
 +
<table style="border-style: none">
 +
<tr align="center">
 +
<td style="border-style: none;">
 +
      <tr>
 +
<td style="border-style:none; float:center">
 +
<img src="https://static.igem.org/mediawiki/2016/3/31/T--HokkaidoU_Japan--Basicpart_Model_fig2.png" width="250px" height="auto" alt="Fig_2"></td>
 +
      </tr>
 +
      <tr>
 +
<td style="border-style: none"; align="center"><span class="small">Fig. 2
 +
 +
</span></td>
 +
      </tr>
 +
</table>
 +
</center>
 +
<br clear="all">
 +
 +
 +
 +
 +
<br>The possibility to fold native conformation (P<span class="sitatuki">HPPH</span> <sup>SAR</sup>) is of course 1.
 +
 +
 +
<br>
 +
<center>
 +
<table style="border-style: none">
 +
<tr align="center">
 +
<td style="border-style: none;">
 +
      <tr>
 +
<td style="border-style:none; float:center">
 +
<img src="https://static.igem.org/mediawiki/2016/c/c9/T--HokkaidoU_Japan--siki3.jpg" width="300px" height="auto" alt="Formula_3"></td>
 +
      </tr>
 +
      <tr>
 +
<td style="border-style: none"; align="center"><h2>
 +
 +
</h2></td>
 +
      </tr>
 +
</table>
 +
</center>
 +
<br clear="all">
 +
 +
 +
 +
<br>Compared with both models, we can obviously think that thanks to the addition of SAR, we can increase the probability to fold correctly; the stability of native state is definitely increased.
 +
 +
<br>
 +
<center>
 +
<table style="border-style: none">
 +
<tr align="center">
 +
<td style="border-style: none;">
 +
      <tr>
 +
<td style="border-style:none; float:center">
 +
<img src="https://static.igem.org/mediawiki/2016/9/93/T--HokkaidoU_Japan--siki4.jpg" width="190px" height="auto" alt="Formula_4"></td>
 +
      </tr>
 +
      <tr>
 +
<td style="border-style: none"; align="center"><h2>
 +
 +
</h2></td>
 +
      </tr>
 +
</table>
 +
</center>
 +
<br clear="all">
 +
 +
 +
 +
 +
<br>Let's think about more complicated case. The number of residue is 6 and the sequence is HPPHPH. The possible states are shown below. Also, we excluded enantiomers and rotamers.
 +
 +
 +
<br>
 +
<center>
 +
<table style="border-style: none">
 +
<tr align="center">
 +
<td style="border-style: none;">
 +
      <tr>
 +
<td style="border-style:none; float:center">
 +
<img src="https://static.igem.org/mediawiki/2016/9/9e/T--HokkaidoU_Japan--Model_fig3.png" width="900px" height="auto" alt="Fig_3"></td>
 +
      </tr>
 +
      <tr>
 +
<td style="border-style: none"; align="center"><span class="small">Fig. 3
 +
 +
</span></td>
 +
      </tr>
 +
</table>
 +
</center>
 +
<br clear="all">
 +
 +
 +
 +
 +
<br>As we did in the simplest model, we
 +
thought the most stable state is the native state;
 +
the native state has -2E<span class="sitatuki">H</span>
 +
as its energy. In this case, the possibility to fold as
 +
native structure (P<span class="sitatuki">HPPHPH</span><sup>wild</sup>) is below.
 +
 +
 +
<br>
 +
<center>
 +
<table style="border-style: none">
 +
<tr align="center">
 +
<td style="border-style: none;">
 +
      <tr>
 +
<td style="border-style:none; float:center">
 +
<img src="https://static.igem.org/mediawiki/2016/7/75/T--HokkaidoU_Japan--siki5.jpg" width="450px" height="auto" alt="Formula_5"></td>
 +
      </tr>
 +
      <tr>
 +
<td style="border-style: none"; align="center"><h2>
 +
 +
</h2></td>
 +
      </tr>
 +
</table>
 +
</center>
 +
<br clear="all">
 +
 +
 +
 +
 +
 +
<br>With SAR, the possible states are shown below.
 +
 +
 +
 +
<br>
 +
<center>
 +
<table style="border-style: none">
 +
<tr align="center">
 +
<td style="border-style: none;">
 +
      <tr>
 +
<td style="border-style:none; float:center">
 +
<img src="https://static.igem.org/mediawiki/2016/1/10/T--HokkaidoU_Japan--Model_fig4.png" width="550px" height="auto" alt="Fig_4"></td>
 +
      </tr>
 +
      <tr>
 +
<td style="border-style: none"; align="center"><span class="small">Fig. 4
 +
 +
</span></td>
 +
      </tr>
 +
</table>
 +
</center>
 +
<br clear="all">
 +
 +
 +
 +
 +
<br>The probability to fold as native structure (P<span class="sitatuki">HPPHPH</span><sup>SAR</sup>) is below.
 +
 +
 +
 +
<br>
 +
<center>
 +
<table style="border-style: none">
 +
<tr align="center">
 +
<td style="border-style: none;">
 +
      <tr>
 +
<td style="border-style:none; float:center">
 +
<img src="https://static.igem.org/mediawiki/2016/7/7b/T--HokkaidoU_Japan--siki6.jpg" width="400px" height="auto" alt="Formula_6"></td>
 +
      </tr>
 +
      <tr>
 +
<td style="border-style: none"; align="center"><h2>
 +
 +
</h2></td>
 +
      </tr>
 +
</table>
 +
</center>
 +
<br clear="all">
 +
 +
 +
 +
 +
<br>As we have shown in the simplest case, by the addition of SAR, the probability to fold correctly is definitely increased.
 +
 +
 +
<br>
 +
<center>
 +
<table style="border-style: none">
 +
<tr align="center">
 +
<td style="border-style: none;">
 +
      <tr>
 +
<td style="border-style:none; float:center">
 +
<img src="https://static.igem.org/mediawiki/2016/e/e6/T--HokkaidoU_Japan--siki7.jpg" width="200px" height="auto" alt="Formula_7"></td>
 +
      </tr>
 +
      <tr>
 +
<td style="border-style: none"; align="center"><h2>
 +
 +
</h2></td>
 +
      </tr>
 +
</table>
 +
</center>
 +
<br clear="all">
 +
 +
 +
 +
 +
<br>As we have shown, circularization using SAR
 +
can stabilize protein native structure. However,
 +
we should be careful about SAR' characteristic;
 +
SAR can limit the structure by circularization,
 +
but of course, if the stabilized structure is different
 +
from native structure, the addition of SAR means that
 +
it increase the stability of the denatured structure.
 +
This can be shown in the model. If we add SAR to the ends
 +
of HPHPHHPPPHHH model, the most stable structure is changed.
 +
 +
 +
<br>
 +
<center>
 +
<table style="border-style: none">
 +
<tr align="center">
 +
<td style="border-style: none;">
 +
      <tr>
 +
<td style="border-style:none; float:center">
 +
<img src="https://static.igem.org/mediawiki/2016/f/fb/T--HokkaidoU_Japan--Model_fig5.png" width="550px" height="auto" alt="Fig_5"></td>
 +
      </tr>
 +
      <tr>
 +
<td style="border-style: none"; align="center"><span class="small">Fig. 5
 +
 +
</span></td>
 +
      </tr>
 +
</table>
 +
</center>
 +
<br clear="all">
 +
 +
 +
 +
<br>As shown above, native state's free energy is -5E<span class="sitatuki">H</span>, but stabilized structure's lowest energy is only -3E<span class="sitatuki">H</span>. This means that if we want to stabilize a protein with circularization using SAR, we have to be careful about the difference between its native structure and stabilized structure. If they have huge difference, we have to add linkers not to break its native structure by circularization using SAR.
 +
 +
</span>
 +
 +
<br>
 +
<br>
 +
<div id="Reference"><img src="https://static.igem.org/mediawiki/2016/8/86/T--HokkaidoU_Japan--reference.png"
 +
width="270px" height="90px" alt="reference"></div>
 +
 +
<span class="nomal2">
 +
<br>[1] Dill K.A. (1985). "Theory for the folding and stability of globular proteins". Biochemistry. 24(6): 1501-9. doi:10.1021/bi00327a032. PMID 3986190.
 +
 +
<br>[2] Rob Philips. (2008) "Physical Biology Of the Cell". Garland Science
 +
</span>
 +
 +
<br>
 +
<br>
 +
 +
 +
</div>
 +
<!-- 2016contents 閉じる -->
 +
 +
<!--begin footer-->
 +
<br>
 +
<div id="footer" style="background-color: #97D3E3; position:relative;width:100%";>
 +
<div id="footer-logo">
 +
<br>
 +
<br>
 +
 +
<a href="https://www.facebook.com/igem.hokkaido.u.bio">
 +
<img style="height:50px;position:relative;" src="https://static.igem.org/mediawiki/2016/9/96/T--HokkaidoU_Japan--facebook.png" alt="Facebook"></a>
 +
 +
<a href="https://twitter.com/igem_hokkaidou">
 +
<img style="height:50px;position:relative;" src="https://static.igem.org/mediawiki/2016/1/13/T--HokkaidoU_Japan--twitter.png" alt="Twitter"></a>
 +
 +
<a href="mailto:igemhokkaidou@gmail.com">
 +
<img style="height:50px;position:relative;" src="https://static.igem.org/mediawiki/2016/4/4e/T--HokkaidoU_Japan--mail.png" alt="E-mail"></a>
 +
 +
<a href="http://igemhokkaidou.wordpress.com">
 +
<img style="height:50px;position:relative;" src="https://static.igem.org/mediawiki/2016/3/3b/T--HokkaidoU_Japan--icon1.png" alt="Web">
 +
</a>
 +
 +
</div>
 +
</div>
 
</div>
 
</div>
 +
<!--background-->
  
 
</html>
 
</html>

Latest revision as of 23:44, 19 October 2016

Team:HokkaidoU Japan - 2016.igem.org

 

Team:HokkaidoU Japan

\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\


Modeling

In our project, we tried to stabilize protein structure by circularization using SAR. About the detailed concept of circularization, please read circularization page. To think about the power of stabilization by circularization using SAP, we used HP (Hydrophobic-Polar) model. HP model is a kind of simplified protein folding model and in this model, protein chain is given as zig-zag stick on 2D lattice. Each residue has the characteristic H or P (Hydrophobic or Polar). To calculate the stability of protein structures, in this model, if an H residue is next to another H residue without covalent bond, it decreases free energy because of hydrophobic interaction. In our model, the decreased energy by each hydrophobic interaction is defined as -EH. We added another characteristic SAR into this model. Through thinking about this model, we can simply think about the effect of SAR reflected as the effect to probability to fold as native state. We thought SAR interaction is so strong, so in the case we add SAR at N terminus and C terminus, both ends are set next to each other in the model. So, let's think about the simplest model.
The simplest model is the model with the number of residue is 4 and the sequence is HPPH. In this case, without SAR, the number of states is 4, excluding enantiomers and rotamers. The possible states and the energy are listed below. KB is Boltzmann constant 1.38064852*10-23 (J/K), T is temperature (K).
Fig_1
Fig. 1


Only the first one is stable and its energy is -EH. Because it's most stable, we thought it's native state. The probability to fold as native state (PwildHPPH) is below.
Formula_1



To calculate this, we used canonical ensemble from statistical mechanics. The probability of causing state i is calculated through the function below. Ei is the energy of state i, Ej is the energy of the state j.
Formula_2



But with SAR, the number of states is 1 and the state is the most stable one.
Fig_2
Fig. 2


The possibility to fold native conformation (PHPPH SAR) is of course 1.
Formula_3



Compared with both models, we can obviously think that thanks to the addition of SAR, we can increase the probability to fold correctly; the stability of native state is definitely increased.
Formula_4



Let's think about more complicated case. The number of residue is 6 and the sequence is HPPHPH. The possible states are shown below. Also, we excluded enantiomers and rotamers.
Fig_3
Fig. 3


As we did in the simplest model, we thought the most stable state is the native state; the native state has -2EH as its energy. In this case, the possibility to fold as native structure (PHPPHPHwild) is below.
Formula_5



With SAR, the possible states are shown below.
Fig_4
Fig. 4


The probability to fold as native structure (PHPPHPHSAR) is below.
Formula_6



As we have shown in the simplest case, by the addition of SAR, the probability to fold correctly is definitely increased.
Formula_7



As we have shown, circularization using SAR can stabilize protein native structure. However, we should be careful about SAR' characteristic; SAR can limit the structure by circularization, but of course, if the stabilized structure is different from native structure, the addition of SAR means that it increase the stability of the denatured structure. This can be shown in the model. If we add SAR to the ends of HPHPHHPPPHHH model, the most stable structure is changed.
Fig_5
Fig. 5


As shown above, native state's free energy is -5EH, but stabilized structure's lowest energy is only -3EH. This means that if we want to stabilize a protein with circularization using SAR, we have to be careful about the difference between its native structure and stabilized structure. If they have huge difference, we have to add linkers not to break its native structure by circularization using SAR.


reference

[1] Dill K.A. (1985). "Theory for the folding and stability of globular proteins". Biochemistry. 24(6): 1501-9. doi:10.1021/bi00327a032. PMID 3986190.
[2] Rob Philips. (2008) "Physical Biology Of the Cell". Garland Science