Difference between revisions of "Task 6: Protein structure prediction from evolutionary sequence variation"

From Bioinformatikpedia
Line 1: Line 1:
  +
== Results ==
 
   
 
[[lab journal]]
 
[[lab journal]]
  +
   
 
== Contact Prediction ==
 
== Contact Prediction ==
Line 94: Line 95:
   
 
<figure id="cm_1a6z">
 
<figure id="cm_1a6z">
[[File:ContactMap_1a6z_new.png|right|thumb|300px| '''Figure 3:''' Predicted (red) and reference(grey) contacts of the MHC I domain (left box) and the Ig domain (right box) of the hemochromatosis protein.]]
+
[[File:ContactMap_1a6z_new.png|right|thumb|300px| '''Figure 4:''' Predicted (red) and reference(grey) contacts of the MHC I domain (left box) and the Ig domain (right box) of the hemochromatosis protein.]]
 
</figure>
 
</figure>
 
<figure id="cm_hras">
 
<figure id="cm_hras">
[[File:ContactMap_5p21.png|left|thumb|300px| '''Figure 4:''' Predicted (red) and reference(grey) contacts of the ras protein.]]
+
[[File:ContactMap_5p21.png|left|thumb|300px| '''Figure 3:''' Predicted (red) and reference(grey) contacts of the ras protein.]]
 
</figure>
 
</figure>
   
Line 116: Line 117:
   
 
<figure id="score_correlation">
 
<figure id="score_correlation">
[[File:Score_correlation_mhc1.png|right|thumb|300px|The distribution of scores for TP and FP predicted contact pairs of MHC I is shown. Note that only high scoring pairs are shown in this plot.]]
+
[[File:Score_correlation_mhc1.png|right|thumb|300px|'''Figure 5:'''' The distribution of scores for TP and FP predicted contact pairs of MHC I is shown. Note that only high scoring pairs are shown in this plot.]]
 
</figure>
 
</figure>
   
Line 124: Line 125:
   
 
*Can you determine evolutionary hot spots, i.e. functionally important residues? Compare to conserved sites in the MSA. Compare with your results from task 7 (when you are working on task 7, i.e. this is a task for the future).
 
*Can you determine evolutionary hot spots, i.e. functionally important residues? Compare to conserved sites in the MSA. Compare with your results from task 7 (when you are working on task 7, i.e. this is a task for the future).
{| class="basic2"
+
{| class="colBasic2"
 
! colspan="5" | MHC || colspan="5" | Ig C1-set
 
! colspan="5" | MHC || colspan="5" | Ig C1-set
 
|-
 
|-
Line 177: Line 178:
   
 
{|
 
{|
| [[File:Ras_di_76.png|center|thumb|300px|DI score and 76 EC constraints]] || [[File:Ras_di_123.png|center|thumb|300px|DI score and 123 EC constraints]] || [[File:Ras_di_189.png|center|thumb|300px|DI score and 189 EC constraints]]
+
| 76 EC constraints || 123 EC constraints || 189 EC constraints
  +
|-
  +
| DI scoring || [[File:Ras_di_76.png|center|thumb|300px]] || [[File:Ras_di_123.png|center|thumb|300px]] || [[File:Ras_di_189.png|center|thumb|300px]]
 
|-
 
|-
| [[File:Ras_plm_76.png|center|thumb|300px|PLM score and 76 EC constraints]] || [[File:Ras_plm_123.png|center|thumb|300px|PLM score and 123 EC constraints]] || [[File:Ras_plm_189.png|center|thumb|300px|PLM score and 189 EC constraints]]
+
| PLM scoring || [[File:Ras_plm_76.png|center|thumb|300px]] || [[File:Ras_plm_123.png|center|thumb|300px]] || [[File:Ras_plm_189.png|center|thumb|300px]]
 
|–
 
|–
|+ style="caption-side: bottom; text-align: left" |<font size=2>'''Table 3:''' ).
+
|+ style="caption-side: bottom; text-align: left" |<font size=2>'''Figure 6:''' Contact maps from the EVfold predictions for Ras are shown.
 
|}
 
|}
 
</figtable>
 
</figtable>

Revision as of 21:48, 27 August 2013


lab journal


Contact Prediction

<css> table.colBasic2 { margin-left: auto; margin-right: auto; border: 1px solid black; border-collapse:collapse; width: 40%; }

.colBasic2 th,td { padding: 3px; border: 1px solid black; }

.colBasic2 td { text-align:left; }

/* for orange try #ff7f00 and #ffaa56 for blue try #005fbf and #aad4ff

maria's style blue: #adceff grey: #efefef

  • /

.colBasic2 tr th { background-color:#efefef; color: black;} .colBasic2 tr:first-child th { background-color:#adceff; color:black;}

table.basic2 { margin-left: auto; margin-right: auto; border: 0px solid black; border-collapse:collapse; width: 40%; }

.basic2 th,td { padding: 3px; border: 1px solid black; }

.basic2 th { border-bottom: 2px solid black; background-color: #fff; }

.basic2 td { text-align:left; }

.basic2 tr:first-child th {

 border-top: 0;

} .basic2 tr:last-child td {

 border-bottom: 0;

} .basic2 tr td:first-child, .basic2 tr th:first-child {

 border-left: 0;

} .basic2 tr td:last-child, .basic2 tr th:last-child {

 border-right: 0;

} </css>


<figtable id="hfe_score_dist" >

Ras
MHC I domain
Ig C1-set domain
Figure 1: CN-score distributions for all three domains.

</figtable>

domain length seq. in alignment reference HS pairs TP FP TP -rate
Ras 160 21151 5P21 65 53 12 0.82
MHC I 174 25167 1A6Z_A 69 29 40 0.42
Ig C1-set 76 16509 1A6Z_A 15 9 6 0.6
Table 1: Overview of the domains of interest.


<figure id="cm_1a6z">

Figure 4: Predicted (red) and reference(grey) contacts of the MHC I domain (left box) and the Ig domain (right box) of the hemochromatosis protein.

</figure> <figure id="cm_hras">

Figure 3: Predicted (red) and reference(grey) contacts of the ras protein.

</figure>

  • Why are the scores of residues close in sequence amongst the highest? Why are the pairs distant in sequence (n>5) more interesting for structure prediction?

It lies in the nature of proteins, that residues that are close in sequence, are also close in structure. Consequently, they are evolutionary coupled and show covariation in the multiple sequence alignment. The pairs that are at least five residues apart in sequence, are more interesting for structure prediction, because they contain more information about the overall topology of the protein, i.e. they reduce the space of possible protein conformations more than pairs that are close in sequence.

  • Look at the values, range and distribution of scores.

For the MHC I domain of HFE_HUMAN, the score distribution is shown in table <xr id="hfe_score_dist" />. The values range from -0.94 to 2.57 with the mean at -0.07. The score distribution corresponds to a slightly right skewed normal distribution, where most values are in the range of -1 to 1. Only 0.5% of scores have a value above 1. Thus, scores with a value greater than one can be considered as high scoring.


  • How many of the high-scoring pairs are true or false positives? Does this correlate with the value of the score? Visualize the predicted contacts together with the crystal structure contacts in a contact map plot.

Table <xr id="hs_table"> shows, that TP-rate can range from 0.8 for Ras to 0.4 for the MHC I domain. Since the number of sequences in the multiple alignment is above 15 000 for all three domains, the TP-rate among the high scoring pairs depends not only on the number of sequences in the alignment, but also on the actual sequence at hand. As discussed in <ref name="EVfold_method"> Marks DS, Colwell LJ, Sheridan R, Hopf TA, Pagnani A, et al. (2011) Protein 3D Structure Computed from Evolutionary Sequence Variation. PLoS ONE 6(12): e28766. doi:10.1371/journal.pone.0028766 </ref> , possible confounding factors could be pyhlogenetic bias or functional constraints from interactions with other molecules.

<figure id="score_correlation">

Figure 5:' The distribution of scores for TP and FP predicted contact pairs of MHC I is shown. Note that only high scoring pairs are shown in this plot.

</figure>

The correlation between CN score and TP/FP contact is not very good as indicated by a Pearson correlation coefficient of 0.354 and the overlap between the boxplots in figure <xr id="score_correlation">.


  • Can you determine evolutionary hot spots, i.e. functionally important residues? Compare to conserved sites in the MSA. Compare with your results from task 7 (when you are working on task 7, i.e. this is a task for the future).
MHC Ig C1-set
AA Pos. Norm. Score SNP Conservation AA Pos Norm. Score SNP Conservation
V 59 9.97 Val->Met (DC) 2 C 225 9.09 - 10
C 124 8.20 - 2 W 239 6.53 - 10
L 33 7.08 - 2 C 282 6.04 Cys -> Tyr (DC) 2
Y 140 6.96 - 3 H 286 5.66 - 5
V 120 5.96 - 0 P 232 5.16 - 10
S 27 5.64 - 0 I 235 4.40 - 9
L 91 5.54 - 0 M 237 3.94 - 9
A 37 5.00 Ala->Val (nDC) 0 I 268 3.91 - 1
L 44 4.99 - 0 W 267 3.51 - 0
G 51 4.71 - 3 D 261 3.10 - 8
Table 2: Residues with the top ten normalised scores for the MHC_I and Ig C1-set domains. For each position, it is shown whether there is a known SNP and if the SNP is known to be disease causing(DC) or not (nDC). The conservation in the Pfam alignment is indicated by a value from 0(not conserved) to 10(highly conserved).


  • Here, the DI score is given. Compare the top 50 DI and CN (from freecontact) scores. How large is the overlap (>80%)?

For RAS_HUMAN, only 20(40%) of the top 50 scores overlaped.



imm EVFOLD

Structural Models

The structural models for Ras were calculated with three different numbers of evolutionary constraints: 76(40%), 123(65%) and 189 (100%).

<figure id="rmsd_evfold">

Comparison of Ca-RMSD values distributions between the different numbers of ECs used.

</figure>



<figtable id="ras_fold_cm">

76 EC constraints 123 EC constraints 189 EC constraints
DI scoring
Ras di 76.png
Ras di 123.png
Ras di 189.png
PLM scoring
Ras plm 76.png
Ras plm 123.png
Ras plm 189.png
Figure 6: Contact maps from the EVfold predictions for Ras are shown.

</figtable>

<figtable id="ras_fold_pymol">

DI scoring, Ca-RMSD=3.75 A
PLM scoring, Ca-RMSD=2.54 A
Figure x: The ras crystal structure(green,5p21) is shown in comparison to the EVfold models with best RMSD from DI scoring(cyan) and PLM scoring(purple). ).

|} </figtable>

Discussion

<references />