Difference between revisions of "Gaucher Disease: Task 03 - Sequence-based predictions"
(→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''' |
|'''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
Contents
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 | - | |||
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 | ||||
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 | ||||
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 | ||||
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/