Difference between revisions of "Homology based structure prediction (Phenylketonuria)"
(→Single template modeling) |
(→Model calculation) |
||
Line 4: | Line 4: | ||
== Model calculation == |
== Model calculation == |
||
[[Lab Journal - Task 5 (PAH) #xxx|Lab journal]] <br> |
[[Lab Journal - Task 5 (PAH) #xxx|Lab journal]] <br> |
||
+ | In our last task we already found some protein structures we now can discover a bit more. Therefore we first build a dataset with four structures with a sequence identity bigger than 60% and another also with four structures and a sequence similarity between 19 and 36% relative to our reference structure 2pah. |
||
+ | |||
=== Modeller === |
=== Modeller === |
||
+ | ... |
||
==== Single template modeling ==== |
==== Single template modeling ==== |
||
− | |||
{| |
{| |
||
|<figure id="mod_1j8u">[[File:Modeller 2pah 1j8u.png|thumb|left|'''<caption>'''Structure of target 2pah (green) and template 1j8u (purple) superimposed. For the creation of the superimposed structure <code>Modeller</code> was used.</caption>]] </figure> |
|<figure id="mod_1j8u">[[File:Modeller 2pah 1j8u.png|thumb|left|'''<caption>'''Structure of target 2pah (green) and template 1j8u (purple) superimposed. For the creation of the superimposed structure <code>Modeller</code> was used.</caption>]] </figure> |
Revision as of 14:03, 6 June 2013
Contents
Summary
...
Model calculation
Lab journal
In our last task we already found some protein structures we now can discover a bit more. Therefore we first build a dataset with four structures with a sequence identity bigger than 60% and another also with four structures and a sequence similarity between 19 and 36% relative to our reference structure 2pah.
Modeller
...
Single template modeling
</figure> </figure> </figure><figure id="mod_1j8u"> |
<figure id="mod_2phm"> |
<figure id="mod_3luy"> |
1j8u: quite good:
The LGA-RMSD of 1j8u is 0.73 and has a sequence identity of 100% to 2pah, for that reason we also modelled 2phm (<xr id="mod_2phm"/>) with a LGA-RMSD of 0.88 and sequence identity of 89.7%.
3luy has a LGA-RMSD of 3.30 and a sequence identity of 22%.
...
Multiple alignments
</figure> </figure> </figure><figure id="mod_1j8u"> |
<figure id="mod_2phm"> |
<figure id="mod_3luy"> |
Swissmodel
SwissModel ... <ref name="swissmodel"> Konstantin Arnold, Lorenza Bordoli, Jürgen Kopp and Torsten Schwede (2006): "The SWISS-MODEL workspace: a web-based environment for protein structure homology modelling". Bioinformatics Vol.22(2): 195-201. doi:10.1093/bioinformatics/bti770 </ref>
</figure> </figure> </figure><figure id="swiss_1j8u"> |
<figure id="swiss_1j8u"> |
<figure id="swiss_3luy"> |
...
iTasser
iTasser... <ref name="itasser"> Yang Zhang (2008): "[zhanglab.ccmb.med.umich.edu/papers/2008_2.pdf I-TASSER server for protein 3D structure prediction]". BMC Bioinformatics Vol.9: 40. doi:10.1186/1471-2105-9-40 </ref>
We used different cutoff values for the modelling:
- for 1J8U we did not set any cutoff since the sequence has a identity of 100%
- for 2PHM we used a cutoff of 90% (all templates are excluded, which have a sequence identity >90% to the query protein)
- for 3LUY we used a cutoff of 35% (exclude all templates which have a sequence identity >35% to the query protein)
Model evaluation
Scores
//TODO How: Calculate the C_alpha RMSD of the models (use /mnt/project/pracstrucfunc13/bin/sap).??? <figtable id="eval_scores">
... | |||||||
---|---|---|---|---|---|---|---|
Template | Seq.Identity | GDT | Cα RMSD | ... | |||
Swiss-Model | |||||||
1j8u | 100 | ? | ? | ||||
2phm | 91.93? | ? | ? | ||||
3luy | 16.79? | ... | ? | ||||
iTasser | |||||||
1j8u | 100 | ? | |||||
2phm | 89? | ? | |||||
3luy | 22 | ? | |||||
Modeller - Single template modeling | |||||||
1j8u | 100 | ? | ? | ||||
2phm | 89? | ? | ? | ||||
3luy | 22 | ? | ? | ||||
Modeller - Multiple alignments | |||||||
1j8u,2phm,1qey,1phz | >60 | ? | ? | ||||
1j8u,2phm,3luy,1wyp | ... | ? | ? | ||||
3luy,1wyp,2v27,2qmx | <30 | ? | ? |
</figtable> Swiss-Model:
- 1j8u:
- Ligands in the template: FE2: 2, H4B: 2.
- Ligands in the model: FE2: 2
- 2phm:
- Ligands in the template: FE: 2.
- Ligands in the model: FE: 2
- 3luy:
- Ligands in the template: PPY: 1.
- Ligands in the model: none.
References
<references/>