Difference between revisions of "Homology Based Structure Predictions Hemochromatosis"

From Bioinformatikpedia
(Evaluation)
(Evaluation)
Line 169: Line 169:
   
 
==== Evaluation ====
 
==== Evaluation ====
  +
   
 
<figtable id="scores_native">
 
<figtable id="scores_native">
Line 174: Line 175:
 
|-
 
|-
 
! style="text-align:left; border-style: solid; border-width: 0 0 2px 0" |Model
 
! style="text-align:left; border-style: solid; border-width: 0 0 2px 0" |Model
! style="text-align:left; border-style: solid; border-width: 0 0 2px 0" |Common residues (TM)
+
! style="text-align:left; border-style: solid; border-width: 0 0 2px 0" |Common residues
 
! style="text-align:left; border-style: solid; border-width: 0 0 2px 0" |TM-Score
 
! style="text-align:left; border-style: solid; border-width: 0 0 2px 0" |TM-Score
 
! style="text-align:left; border-style: solid; border-width: 0 0 2px 0" |GDT-TS-Score
 
! style="text-align:left; border-style: solid; border-width: 0 0 2px 0" |GDT-TS-Score
 
! style="text-align:left; border-style: solid; border-width: 0 0 2px 0" |GDT-HA-Score
 
! style="text-align:left; border-style: solid; border-width: 0 0 2px 0" |GDT-HA-Score
! style="text-align:left; border-style: solid; border-width: 0 0 2px 0" |Common residues (SAP)
+
! style="text-align:left; border-style: solid; border-width: 0 0 2px 0" |TM-Align
  +
! style="text-align:left; border-style: solid; border-width: 0 0 2px 0" |Nr. of atoms
 
! style="text-align:left; border-style: solid; border-width: 0 0 2px 0" |Weighted RMSd
 
! style="text-align:left; border-style: solid; border-width: 0 0 2px 0" |Weighted RMSd
 
|-
 
|-
Line 186: Line 188:
 
| style="text-align:center; border-style: solid; border-width: 0 0 0 0" |0.0846
 
| style="text-align:center; border-style: solid; border-width: 0 0 0 0" |0.0846
 
| style="text-align:center; border-style: solid; border-width: 0 0 0 0" |0.0487
 
| style="text-align:center; border-style: solid; border-width: 0 0 0 0" |0.0487
  +
| style="text-align:center; border-style: solid; border-width: 0 0 0 0" |0.83298
 
| style="text-align:center; border-style: solid; border-width: 0 0 0 0" |272
 
| style="text-align:center; border-style: solid; border-width: 0 0 0 0" |272
 
| style="text-align:center; border-style: solid; border-width: 0 0 0 0" |2.200
 
| style="text-align:center; border-style: solid; border-width: 0 0 0 0" |2.200
Line 194: Line 197:
 
| style="text-align:center; border-style: solid; border-width: 0 0 0 0" |0.0846
 
| style="text-align:center; border-style: solid; border-width: 0 0 0 0" |0.0846
 
| style="text-align:center; border-style: solid; border-width: 0 0 0 0" |0.0506
 
| style="text-align:center; border-style: solid; border-width: 0 0 0 0" |0.0506
  +
| style="text-align:center; border-style: solid; border-width: 0 0 0 0" |0.83358
 
| style="text-align:center; border-style: solid; border-width: 0 0 0 0" |272
 
| style="text-align:center; border-style: solid; border-width: 0 0 0 0" |272
 
| style="text-align:center; border-style: solid; border-width: 0 0 0 0" |2.325
 
| style="text-align:center; border-style: solid; border-width: 0 0 0 0" |2.325
Line 202: Line 206:
 
| style="text-align:center; border-style: solid; border-width: 0 0 0 0" |0.0827
 
| style="text-align:center; border-style: solid; border-width: 0 0 0 0" |0.0827
 
| style="text-align:center; border-style: solid; border-width: 0 0 0 0" |0.0432
 
| style="text-align:center; border-style: solid; border-width: 0 0 0 0" |0.0432
  +
| style="text-align:center; border-style: solid; border-width: 0 0 0 0" |0.85494
 
| style="text-align:center; border-style: solid; border-width: 0 0 0 0" |272
 
| style="text-align:center; border-style: solid; border-width: 0 0 0 0" |272
 
| style="text-align:center; border-style: solid; border-width: 0 0 0 0" |1.642
 
| style="text-align:center; border-style: solid; border-width: 0 0 0 0" |1.642
Line 210: Line 215:
 
| style="text-align:center; border-style: solid; border-width: 0 0 0 0" |0.0790
 
| style="text-align:center; border-style: solid; border-width: 0 0 0 0" |0.0790
 
| style="text-align:center; border-style: solid; border-width: 0 0 0 0" |0.0423
 
| style="text-align:center; border-style: solid; border-width: 0 0 0 0" |0.0423
  +
| style="text-align:center; border-style: solid; border-width: 0 0 0 0" |0.79841
 
| style="text-align:center; border-style: solid; border-width: 0 0 0 0" |271
 
| style="text-align:center; border-style: solid; border-width: 0 0 0 0" |271
 
| style="text-align:center; border-style: solid; border-width: 0 0 0 0" |2.302
 
| style="text-align:center; border-style: solid; border-width: 0 0 0 0" |2.302
Line 218: Line 224:
 
| style="text-align:center; border-style: solid; border-width: 0 0 0 0" |0.0836
 
| style="text-align:center; border-style: solid; border-width: 0 0 0 0" |0.0836
 
| style="text-align:center; border-style: solid; border-width: 0 0 0 0" |0.0460
 
| style="text-align:center; border-style: solid; border-width: 0 0 0 0" |0.0460
  +
| style="text-align:center; border-style: solid; border-width: 0 0 0 0" |0.49103
 
| style="text-align:center; border-style: solid; border-width: 0 0 0 0" |269
 
| style="text-align:center; border-style: solid; border-width: 0 0 0 0" |269
 
| style="text-align:center; border-style: solid; border-width: 0 0 0 0" |2.166
 
| style="text-align:center; border-style: solid; border-width: 0 0 0 0" |2.166
Line 226: Line 233:
 
| style="text-align:center; border-style: solid; border-width: 0 0 0 0" |0.0607
 
| style="text-align:center; border-style: solid; border-width: 0 0 0 0" |0.0607
 
| style="text-align:center; border-style: solid; border-width: 0 0 0 0" |0.0349
 
| style="text-align:center; border-style: solid; border-width: 0 0 0 0" |0.0349
  +
| style="text-align:center; border-style: solid; border-width: 0 0 0 0" |0.40162
 
| style="text-align:center; border-style: solid; border-width: 0 0 0 0" |272
 
| style="text-align:center; border-style: solid; border-width: 0 0 0 0" |272
 
| style="text-align:center; border-style: solid; border-width: 0 0 0 0" |3.705
 
| style="text-align:center; border-style: solid; border-width: 0 0 0 0" |3.705
Line 234: Line 242:
 
| style="text-align:center; border-style: solid; border-width: 0 0 0 0" |0.0892
 
| style="text-align:center; border-style: solid; border-width: 0 0 0 0" |0.0892
 
| style="text-align:center; border-style: solid; border-width: 0 0 0 0" |0.0496
 
| style="text-align:center; border-style: solid; border-width: 0 0 0 0" |0.0496
  +
| style="text-align:center; border-style: solid; border-width: 0 0 0 0" |0.81698
 
| style="text-align:center; border-style: solid; border-width: 0 0 0 0" |267
 
| style="text-align:center; border-style: solid; border-width: 0 0 0 0" |267
 
| style="text-align:center; border-style: solid; border-width: 0 0 0 0" |2.374
 
| style="text-align:center; border-style: solid; border-width: 0 0 0 0" |2.374
Line 242: Line 251:
 
| style="text-align:center; border-style: solid; border-width: 0 0 0 0" |0.0882
 
| style="text-align:center; border-style: solid; border-width: 0 0 0 0" |0.0882
 
| style="text-align:center; border-style: solid; border-width: 0 0 0 0" |0.0496
 
| style="text-align:center; border-style: solid; border-width: 0 0 0 0" |0.0496
  +
| style="text-align:center; border-style: solid; border-width: 0 0 0 0" |0.86512
 
| style="text-align:center; border-style: solid; border-width: 0 0 0 0" |267
 
| style="text-align:center; border-style: solid; border-width: 0 0 0 0" |267
 
| style="text-align:center; border-style: solid; border-width: 0 0 0 0" |1.524
 
| style="text-align:center; border-style: solid; border-width: 0 0 0 0" |1.524
Line 253: Line 263:
 
|-
 
|-
 
! style="text-align:left; border-style: solid; border-width: 0 0 2px 0" |Model
 
! style="text-align:left; border-style: solid; border-width: 0 0 2px 0" |Model
! style="text-align:left; border-style: solid; border-width: 0 0 2px 0" |Common residues (TM)
+
! style="text-align:left; border-style: solid; border-width: 0 0 2px 0" |Common residues
 
! style="text-align:left; border-style: solid; border-width: 0 0 2px 0" |TM-Score
 
! style="text-align:left; border-style: solid; border-width: 0 0 2px 0" |TM-Score
 
! style="text-align:left; border-style: solid; border-width: 0 0 2px 0" |GDT-TS-Score
 
! style="text-align:left; border-style: solid; border-width: 0 0 2px 0" |GDT-TS-Score
 
! style="text-align:left; border-style: solid; border-width: 0 0 2px 0" |GDT-HA-Score
 
! style="text-align:left; border-style: solid; border-width: 0 0 2px 0" |GDT-HA-Score
! style="text-align:left; border-style: solid; border-width: 0 0 2px 0" |Common residues (SAP)
+
! style="text-align:left; border-style: solid; border-width: 0 0 2px 0" |TM-Align
  +
! style="text-align:left; border-style: solid; border-width: 0 0 2px 0" |Nr. of atoms
 
! style="text-align:left; border-style: solid; border-width: 0 0 2px 0" |Weighted RMSd
 
! style="text-align:left; border-style: solid; border-width: 0 0 2px 0" |Weighted RMSd
 
|-
 
|-
Line 265: Line 276:
 
| style="text-align:center; border-style: solid; border-width: 0 0 0 0" |0.0800
 
| style="text-align:center; border-style: solid; border-width: 0 0 0 0" |0.0800
 
| style="text-align:center; border-style: solid; border-width: 0 0 0 0" |0.0469
 
| style="text-align:center; border-style: solid; border-width: 0 0 0 0" |0.0469
  +
| style="text-align:center; border-style: solid; border-width: 0 0 0 0" |0.85238
 
| style="text-align:center; border-style: solid; border-width: 0 0 0 0" |272
 
| style="text-align:center; border-style: solid; border-width: 0 0 0 0" |272
 
| style="text-align:center; border-style: solid; border-width: 0 0 0 0" |1.986
 
| style="text-align:center; border-style: solid; border-width: 0 0 0 0" |1.986
Line 273: Line 285:
 
| style="text-align:center; border-style: solid; border-width: 0 0 0 0" |0.0800
 
| style="text-align:center; border-style: solid; border-width: 0 0 0 0" |0.0800
 
| style="text-align:center; border-style: solid; border-width: 0 0 0 0" |0.0478
 
| style="text-align:center; border-style: solid; border-width: 0 0 0 0" |0.0478
  +
| style="text-align:center; border-style: solid; border-width: 0 0 0 0" |0.83836
 
| style="text-align:center; border-style: solid; border-width: 0 0 0 0" |272
 
| style="text-align:center; border-style: solid; border-width: 0 0 0 0" |272
 
| style="text-align:center; border-style: solid; border-width: 0 0 0 0" |2.142
 
| style="text-align:center; border-style: solid; border-width: 0 0 0 0" |2.142
Line 281: Line 294:
 
| style="text-align:center; border-style: solid; border-width: 0 0 0 0" |0.0827
 
| style="text-align:center; border-style: solid; border-width: 0 0 0 0" |0.0827
 
| style="text-align:center; border-style: solid; border-width: 0 0 0 0" |0.0450
 
| style="text-align:center; border-style: solid; border-width: 0 0 0 0" |0.0450
  +
| style="text-align:center; border-style: solid; border-width: 0 0 0 0" |0.86026
 
| style="text-align:center; border-style: solid; border-width: 0 0 0 0" |272
 
| style="text-align:center; border-style: solid; border-width: 0 0 0 0" |272
 
| style="text-align:center; border-style: solid; border-width: 0 0 0 0" |1.720
 
| style="text-align:center; border-style: solid; border-width: 0 0 0 0" |1.720
Line 289: Line 303:
 
| style="text-align:center; border-style: solid; border-width: 0 0 0 0" |0.0790
 
| style="text-align:center; border-style: solid; border-width: 0 0 0 0" |0.0790
 
| style="text-align:center; border-style: solid; border-width: 0 0 0 0" |0.0414
 
| style="text-align:center; border-style: solid; border-width: 0 0 0 0" |0.0414
  +
| style="text-align:center; border-style: solid; border-width: 0 0 0 0" |0.82124
 
| style="text-align:center; border-style: solid; border-width: 0 0 0 0" |271
 
| style="text-align:center; border-style: solid; border-width: 0 0 0 0" |271
 
| style="text-align:center; border-style: solid; border-width: 0 0 0 0" |2.059
 
| style="text-align:center; border-style: solid; border-width: 0 0 0 0" |2.059
Line 297: Line 312:
 
| style="text-align:center; border-style: solid; border-width: 0 0 0 0" |0.0836
 
| style="text-align:center; border-style: solid; border-width: 0 0 0 0" |0.0836
 
| style="text-align:center; border-style: solid; border-width: 0 0 0 0" |0.0441
 
| style="text-align:center; border-style: solid; border-width: 0 0 0 0" |0.0441
  +
| style="text-align:center; border-style: solid; border-width: 0 0 0 0" |0.48549
 
| style="text-align:center; border-style: solid; border-width: 0 0 0 0" |269
 
| style="text-align:center; border-style: solid; border-width: 0 0 0 0" |269
 
| style="text-align:center; border-style: solid; border-width: 0 0 0 0" |2.279
 
| style="text-align:center; border-style: solid; border-width: 0 0 0 0" |2.279
Line 305: Line 321:
 
| style="text-align:center; border-style: solid; border-width: 0 0 0 0" |0.0653
 
| style="text-align:center; border-style: solid; border-width: 0 0 0 0" |0.0653
 
| style="text-align:center; border-style: solid; border-width: 0 0 0 0" |0.0358
 
| style="text-align:center; border-style: solid; border-width: 0 0 0 0" |0.0358
  +
| style="text-align:center; border-style: solid; border-width: 0 0 0 0" |0.40088
 
| style="text-align:center; border-style: solid; border-width: 0 0 0 0" |272
 
| style="text-align:center; border-style: solid; border-width: 0 0 0 0" |272
 
| style="text-align:center; border-style: solid; border-width: 0 0 0 0" |3.655
 
| style="text-align:center; border-style: solid; border-width: 0 0 0 0" |3.655
Line 313: Line 330:
 
| style="text-align:center; border-style: solid; border-width: 0 0 0 0" |0.0873
 
| style="text-align:center; border-style: solid; border-width: 0 0 0 0" |0.0873
 
| style="text-align:center; border-style: solid; border-width: 0 0 0 0" |0.0487
 
| style="text-align:center; border-style: solid; border-width: 0 0 0 0" |0.0487
  +
| style="text-align:center; border-style: solid; border-width: 0 0 0 0" |0.83662
 
| style="text-align:center; border-style: solid; border-width: 0 0 0 0" |269
 
| style="text-align:center; border-style: solid; border-width: 0 0 0 0" |269
 
| style="text-align:center; border-style: solid; border-width: 0 0 0 0" |2.088
 
| style="text-align:center; border-style: solid; border-width: 0 0 0 0" |2.088
Line 321: Line 339:
 
| style="text-align:center; border-style: solid; border-width: 0 0 0 0" |0.0901
 
| style="text-align:center; border-style: solid; border-width: 0 0 0 0" |0.0901
 
| style="text-align:center; border-style: solid; border-width: 0 0 0 0" |0.0515
 
| style="text-align:center; border-style: solid; border-width: 0 0 0 0" |0.0515
  +
| style="text-align:center; border-style: solid; border-width: 0 0 0 0" |0.86877
 
| style="text-align:center; border-style: solid; border-width: 0 0 0 0" |269
 
| style="text-align:center; border-style: solid; border-width: 0 0 0 0" |269
 
| style="text-align:center; border-style: solid; border-width: 0 0 0 0" |1.492
 
| style="text-align:center; border-style: solid; border-width: 0 0 0 0" |1.492

Revision as of 11:16, 1 June 2012

Hemochromatosis>>Task 4: Homology based structure predictions

Riddle of the task

After endless battles and deadly traps you have finally reached the tomb's final chamber. As you enter it you notice that there is no sign of the treasures that were promised by the old map you found months ago. Suddenly you hear a loud noise behind you and a solid wall of stone blocks the only entrance into the room. You are trapped! You look around and notice something on the walls. On the left wall are four runes in an ancient language. Luckily its the same language as the notes on the map you deciphered and they translate into four single letters:

  • C N O I

On the opposite wall you can see four simple symbols:

  • a triangle
  • a square
  • a circle
  • and a diamond (dt. Raute)

After further investigation you notice that the four symbols can be pushed into the wall, but you don't know what would happen and which one(s) to push.

What do you do?


Some hints:

  • You only have to push one button and only once.

Short Task Description

Detailed description: Homology based structure predictions

Protocol

A protocol with a description of the data acquisition and other scripts used for this task is available here.

PDB templates

In order to find templates for our models we performed several searches for homologs with COMA and HHPred. We also reused the sequences from Task 2. However none of these methods yielded homologs with a sequence identity above 40% (except 1a6z which is HFE itself) that could be mapped to a PDB structure. The best results for COMA and HHPred are listed in <xr id="coma_t"/> and <xr id="hhpred_t"/> respectively. Therefore we could not generate models with a sequence identity above 80% and the 40%-80% range was limited to its lower bound.

<figtable id="coma_t">

PDB ID e-Value Identities Positives
1a6z_A 1.00E-63 100% 100%
1t7v_A 1.20E-63 34% 62%
3nwm_A 1.60E-57 30% 57%
1frt_A 4.90E-65 28% 66%
2wy3_A 2.80E-63 26% 66%
3ov6_A 2.70E-55 20% 67%
1u58_A 1.00E-54 19% 59%
3d2u_A 7.50E-60 16% 70%
3dbx_A 9.70E-59 16% 71%
3it8_D 3.40E-52 15% 65%
TODO: COMA: Top10 (e-Value) sorted by Identities.

</figtable>

<figtable id="hhpred_t">

PDB ID e-Value Identities Similarity
1a6z_A 1.80E-69 100% 1.623
1k5n_A 2.80E-68 40% 0.725
1s7q_A 5.80E-78 37% 0.655
1t7v_A 7.40E-68 36% 0.702
3p73_A 1.10E-69 35% 0.638
3bev_A 1.00E-69 34% 0.684
2yf1_A 2.40E-74 32% 0.617
1zs8_A 7.30E-68 30% 0.553
2wy3_A 3.30E-68 29% 0.496
1cd1_A 1.60E-67 21% 0.394
TODO: HHpred: Top10 (e-Value) sorted by Identities.

</figtable>


Model generation


Modeller


Evaluation

<figtable id="scores_native">

Model Common residues TM-Score GDT-TS-Score GDT-HA-Score TM-Align Nr. of atoms Weighted RMSd
1K5N_2d 272 0.1686 0.0846 0.0487 0.83298 272 2.200
1K5N_simple 272 0.1649 0.0846 0.0506 0.83358 272 2.325
1ZS8_2d 272 0.1698 0.0827 0.0432 0.85494 272 1.642
1ZS8_simple 272 0.1550 0.0790 0.0423 0.79841 271 2.302
2IAD_2d 272 0.1725 0.0836 0.0460 0.49103 269 2.166
2IAD_simple 272 0.1337 0.0607 0.0349 0.40162 272 3.705
3DBX_2d 272 0.1742 0.0892 0.0496 0.81698 267 2.374
3DBX_simple 272 0.1684 0.0882 0.0496 0.86512 267 1.524
TODO: NATIVE.

</figtable>

<figtable id="scores_complex">

Model Common residues TM-Score GDT-TS-Score GDT-HA-Score TM-Align Nr. of atoms Weighted RMSd
1K5N_2d 272 0.1674 0.0800 0.0469 0.85238 272 1.986
1K5N_simple 272 0.1638 0.0800 0.0478 0.83836 272 2.142
1ZS8_2d 272 0.1715 0.0827 0.0450 0.86026 272 1.720
1ZS8_simple 272 0.1550 0.0790 0.0414 0.82124 271 2.059
2IAD_2d 272 0.1706 0.0836 0.0441 0.48549 269 2.279
2IAD_simple 272 0.1318 0.0653 0.0358 0.40088 272 3.655
3DBX_2d 272 0.1737 0.0873 0.0487 0.83662 269 2.088
3DBX_simple 272 0.1687 0.0901 0.0515 0.86877 269 1.492
TODO: COMPLEX.

</figtable>


SwissModel


Evaluation


I-Tasser


Evaluation