Difference between revisions of "Protein structure prediction from evolutionary sequence variation (Phenylketonuria)"
(→PAH) |
(→H-Ras) |
||
Line 117: | Line 117: | ||
<figure id="cn_distr_ras"> |
<figure id="cn_distr_ras"> |
||
− | [[File:2pah_rash_CN.png|thumb|'''<caption>''' CN-score distribution of H-Ras calculated with <code>freecontact</code> for all (green) and extracted residue pairs with a distance of at least six amino acids (purple).</caption>]] |
+ | [[File:2pah_rash_CN.png|thumb|300px|'''<caption>''' CN-score distribution of H-Ras calculated with <code>freecontact</code> for all (green) and extracted residue pairs with a distance of at least six amino acids (purple).</caption>]] |
</figure> |
</figure> |
||
Revision as of 15:04, 15 June 2013
Contents
Summary
...
Multiple Sequence Alignment
Lab journal
The multiple alignment of Ras (PF00071) was downloaded from Pfam.
For our protein two domains are included (ACT and Biopterin domain).
Results
Freecontact
is a program that takes a multiple alignment and can calculate the mutual information score and the corrected norm contact score.
- MI:...
- CN:...
H-Ras
<figtable id="cn-ras">
Ten results with highest CN for all and for extracted pairs | ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|
|
|
|
The ten results of freecontact
calculated for H-RAS with highest CN-value first for all pairs and second only for pairs with a distance of at least five residues. The third table represents the ten best residue pairs with DI score calculated with EVcouplings
. All three tables show the position of the first residue in column one with its corresponding amino acid (column 2) and of the second residue in column three with its amino acid in column four. The next one represents the mutual information score (MI) and the last one the corrected norm contact score(CN).
</figtable>
<figure id="cn_distr_ras">
</figure>
Comparing the best ten CN-scores for all and for residue pairs with at least a distant of more than five amino acids it is remarkable that for all pairs the ten best all are directly neighboring amino acids and have almost twice of the CN-score. The range of the CN-score goes from -0.65 to 3.40 for the extracted residue pairs. The distribution can be viewed in <xr id="cn_distr_ras"/>. Additionally the distribution of all residue pairs is shown. They are very similar, nevertheless the range of all pairs goes up to 6.01. For the DI-score calculation using EVcouplings
a range between 0.00 and 0.22 was achieved. Five of the ten best results are in common for extracted residue pairs found with freecontact
and EVcouplings
, which are 10(G)-16(K), 11(A)-92(D), 81(V)-116(N), 82(F)-141(Y) and 130(A)-141(Y).
PAH
<figtable id="cn-pah">
Ten results with highest CN for all and for extracted pairs | ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|
|
|
|
The ten results of freecontact
calculated for the biopterin domain of PAH with highest CN-value first for all pairs and second only for pairs with a distance of at least five residues. The third table represents the best ten DI-scores also for extracted pairs calculated by EVcouplings. All three tables show the position of the first residue in column one with its corresponding amino acid (column 2) and of the second residue in column three with its amino acid in column four. The next one represents the mutual information score (MI) and the last one the corrected norm contact score(CN).
</figtable>
<figure id="cn_distr_pah">
</figure>
For the biopterin-domain the CN-scores of freecontact
for all residue pairs are ranged between -1.04 and 3.44, for extracted pairs they are between -1.04 and 2.46 and for the DI scores calculated with EVcouplings
can be found in a range between 0.00 and 0.25. In <xr id="cn_distr_pah"/> the distribution of the CN-scores for all and extracted residue pairs can be viewed.
Comparing the ten best results for the extracted residue pairs calculated with freecontact
and with EVcouplings
you can see that four results can be found in both, which are the residue pairs 342(A)-354(L), 257(G)-264(H), 282(D)-290(H) and 192(K)-221(E).
Discussion
Calculate structural model
H-Ras
<figure id="ras_evfold">
</figure>
PAH
<figure id="biopterin_evfold">
</figure>
References
<references/>