Difference between revisions of "Gaucher Disease: Task 03 - Sequence-based predictions"

From Bioinformatikpedia
(Transmembrane Helices)
(Transmembrane Helices)
Line 198: Line 198:
 
|-
 
|-
 
|align="center" |TMH Topology
 
|align="center" |TMH Topology
|''' ''' <span style="color:#045FB4"> 35-56</span><br>''' <span style="color:#045FB4"> 71-89</span><br><br>113-136<br>''' <span style="color:#045FB4"> 157-178</span><br>''' <span style="color:#045FB4"> 190-205</span><br><br>232-252'''
+
|''' <span style="color:#045FB4"> 35-56</span><br> <span style="color:#045FB4"> 71-89</span><br><br>113-136<br> <span style="color:#045FB4"> 157-178</span><br> <span style="color:#045FB4"> 190-205</span><br><br>232-252'''
 
|'''34-58<br> 70-91<br><br> 115-136<br>156-177<br>188-208<br><br> 231-252'''
 
|'''34-58<br> 70-91<br><br> 115-136<br>156-177<br>188-208<br><br> 231-252'''
 
|'''34-56<br> 70-88<br><span style="color:#a30909"> 98-107</span><br> 112-136<br> 156-178<br> 189-203<br> <span style="color:#a30909"> 214-223</span><br> 231-252'''
 
|'''34-56<br> 70-88<br><span style="color:#a30909"> 98-107</span><br> 112-136<br> 156-178<br> 189-203<br> <span style="color:#a30909"> 214-223</span><br> 231-252'''
Line 240: Line 240:
 
|-
 
|-
 
|more information
 
|more information
  +
|[http://bioinf.cs.ucl.ac.uk/psipred/result/18e493a6-c560-11e2-9aac-00163e110593 MemsatSVM]
|
 
 
|
 
|
 
|[http://opm.phar.umich.edu/protein.php?search=2d57 2D57]
 
|[http://opm.phar.umich.edu/protein.php?search=2d57 2D57]

Revision as of 15:32, 27 May 2013

Secondary Structure

In this task secondary structure is predicted using ReProf and PsiPred and compared to DSSP structure assignment. ReProf uses a fasta sequence or a PSI-BLAST PSSM for prediction, PsiPred a fasta sequence and DSSP server needs a PDB file in order to use the 3D coordinates of atoms. The predictions were made for the proteins below. If several PDB structures are available, the one covering the longest sequence and with the best resolution was chosen. For glucosylceramidase the structure 1OGS was used (as in the previous week task).

Uniprot PDB
Entry Protein name Origin Length Entry Method Resolution (Å) Chain Positions
P10775 Ribonuclease inhibitor pig 456 AA 2BNH X-ray 2.30 A 1-456
Q9X0E6 Divalent-cation tolerance protein CutA bacterium Thermotoga maritima 101 AA 1KR4 X-ray 1.40 A 1-101
Q08209 Serine/threonine-protein phosphatase 2B catalytic subunit alpha isoform, EC=3.1.3.16 human 521 AA 1AUI X-ray 2.10 A 1-521
P04062 Glucosylceramidase/acid-beta-glucosidase, EC=3.2.1.45 human 536 AA 1OGS X-ray 2.00 A/B 40-536

TODO: What features are predicted? Discuss the results for your protein and the example proteins. Using the predictions, what could you learn about your protein and the example proteins? Compare to the available knowledge in UniProt, PDB, DisProt, OPM, PDBTM, Pfam...

Disorder

Transmembrane Helices

Four Proteins, including the Gaucher's disease causing Protein, where analysed under reference by transmembrane helices. The used prediction tools differ in their analysing features. While Polyphobius only differs between residues being part of a transmembrane helix or being inside/outside of the cytoplama, Memsat-SVM also predicts re-entrant helices and pore-linig helices. Due to the fact that pore-lining helices are also transmembrane helices, this kind of helices is detected of both prediction tools. In case of re-entrant helices both programms differ. Polyphobius takes this helix as a membrane helix, but Memsat-SVM predicts a re-entrant helix outside the membrane. Therefore the number of the transmembrane hlices, predicted by the tools, differ from each other. In this case also the C-terminal or the N-terminal may be predicted in a different area because of an additional helix.

Comparison of TMH for P04062
Prediction Assignment
Memsat SVM Polyphobius OPM PDMTM
# of TMH 1 - - -
TMH Topology 456-471 - - -
N-terminal extracellular extracellular extracellular -
C-terminal cytoplasmic extracellular extracellular -
Signal peptide 1-34 1-40 - -
Re-entrant Helix - - - -
Pore-lining Helix 1 - - -
Graphical position
Cartoon P04062.png
Graphik P04062.png
1ogs.png
-
more information MemsatSVM 1OGS 1OGS is not in the PDBTM


Comparison of TMH for Q9YDF8
Prediction Assignment
Memsat SVM Polyphobius OPM PDMTM
# of TMH 6 7 5 4
TMH Topology 43-59
72-90
101-118
128-143
163-184

221-245
42-60
68-88
108-129
137-157
163-184
196-213
224-244
25-46
55-78
86-97
100-107
117-148
27-50
55-75
88-107
118-142
N-terminal cytoplasmic extracellular cytoplasmic extracellular
C-terminal cytoplasmic cytoplasmic cytoplasmic extracellular
Signal peptide - -
Re-entrant Helix 188-217 -
Pore-lining Helix 1 -
Graphical position
Cartoon Q9YDF8.png
Graphik Q9YDF8.png
1ors-Q9YDF.png
1ors lm.png
more information MemsatSVM 1ORS 1ORS
Comparison of TMH for P47863
Prediction Assignment
Memsat SVM Polyphobius OPM PDMTM
# of TMH 6 6 8 (per chain)
TMH Topology 35-56
71-89

113-136
157-178
190-205

232-252
34-58
70-91

115-136
156-177
188-208

231-252
34-56
70-88
98-107
112-136
156-178
189-203
214-223
231-252
N-terminal cytoplasmic cytoplasmic cytoplasmic
C-terminal cytoplasmic cytoplasmic extracellular
Signal peptide 1-20
Re-entrant Helix 93-109
209-225
Pore-lining Helix 3
Graphical position
Graphik P47863.png
2d57.png
more information MemsatSVM 2D57
Comparison of TMH for P35462
Prediction Assignment
Memsat SVM Polyphobius OPM PDMTM
# of TMH 7 7
TMH Topology 30-55
66-88
105-126
150-170
188-212
329-352
367-386
34-52
67-91
101-126
150-170
187-209
330-351
363-386
N-terminal extracellular extracellular
C-terminal cytoplasmic
Signal peptide
Re-entrant Helix
Pore-lining Helix
Graphical position
Graphik P35462.png
3pbl.png
more information 3PBL

Signal Peptides

GO Terms

Discussion

Other available methods

Prediction of Tool Information
secondary structure GOR http://gor.bb.iastate.edu/
disorder DISOPRED2 http://bioinf.cs.ucl.ac.uk/psipred/
transmembrane helices MEMSAT3 http://bioinf.cs.ucl.ac.uk/psipred/
TMHMM http://www.cbs.dtu.dk/services/TMHMM/
signal peptides
GO terms

What else can/is be predicted from protein sequence alone

  • Fold recognition (profile based pGenTHREADER and rapid GenTHREADER)
  • Fold domain recognition (pDomTHREADER)
  • Protein domain prediction (DomPred)
  • Homology modelling (BioSerf v2.0)
  • Function prediction (eukaryotic function: FFPred v2.0)
  • Prediction of TM topology and helix packing (SVM-based MEMPACK)

http://bioinf.cs.ucl.ac.uk/psipred/

Which predictions can be improved considerably by structure-based approaches