Difference between revisions of "Task 3 - Sequence-based predictions"
(→Transmembrane helices) |
(→Signal peptides) |
||
Line 58: | Line 58: | ||
You can run <code>signalp</code> in the student computer pool. This is right now version 3.0. However, Tim will update to version 4.0 soon. Please note which version you are using. You can also use the [http://www.cbs.dtu.dk/services/SignalP/ SignalP] server. This is version 4.0. For one of the example proteins I got a different prediction from SignalP 3.0 and 4.0. |
You can run <code>signalp</code> in the student computer pool. This is right now version 3.0. However, Tim will update to version 4.0 soon. Please note which version you are using. You can also use the [http://www.cbs.dtu.dk/services/SignalP/ SignalP] server. This is version 4.0. For one of the example proteins I got a different prediction from SignalP 3.0 and 4.0. |
||
− | You can use the [http://www.signalpeptide.de/index.php Signal Peptide Website] to look up the proteins. |
+ | You can use the [http://www.signalpeptide.de/index.php Signal Peptide Website] to look up the proteins. Check also for predicted transmembrane helices using PolyPhobius. |
== GO terms == |
== GO terms == |
Revision as of 14:23, 11 May 2012
In contrast to the vast amount of known protein sequences, information about structure and function is available for only very few proteins. Sequence-based predictions of protein features aim to decrease this gap. Many sequence-based prediction methods use evolutionary information. Sequence alignments are therefore often a prerequisite for the predictions.
Contents
Theoretical background talks
The talks will give an introduction to sequence-based protein predictions. In particular:
- secondary structure
- disorder
- transmembrane helices
- GO terms
The slides are available here [1] (sorry too large to upload to the wiki system).
Where to run the jobs
- You can log in to the student computer pool from outside:
i12k-biolab??.informatik.tu-muenchen.de
, where??
goes from 01 to 10. - Work in the student computer pool.
- You can also install the programs on your own computer.
Secondary structure
Use ReProf (available as Debian package from rostlab) to predict secondary structure for your protein. Apply ReProf also to these proteins (given are UniProt IDs):
- P10775
- Q9X0E6
- Q08209
Use fasta sequences for the prediction. You can find out about Reprof usage by running reprof
or reading the man page (man reprof
). Peter Hoenigschmig (hoenigschmid@rostlab.org
) would like to hear about anything that would improve the description or if anything seems unclear. For help, you can always ask us first.
Compare the ReProf results to PsiPred and [ http://mrs.cmbi.ru.nl/hsspsoap/ DSSP_server] (DSSP). Before you use DSSP, find out more about the example proteins (and yours) using UniProt and the PDB.
Disorder
Use IUPred to predict disorder for your protein. Apply IUPred to the example proteins given above, too (run iupred
). You can find a more information (README and example) here: /opt/iupred/
.
Compare the results to the information in the DisProt database.
Transmembrane helices
Use PolyPhobius to predict transmembrane helices for your protein and for the following proteins (UniProt IDs given):
- P35462
- Q9YDF8
- P47863
PolyPhobius is installed in /mnt/project/pracstrucfunc12/polyphobius/
.
In contrast to its precursor Phobius, PolyPhobius uses homology information for the prediction. First, you have to execute a blast search. PolyPhobius distributed its own perl script for this purpose: blastget
(/mnt/project/pracstrucfunc12/polyphobius/blastget
). Usage: blastget -h
. Use only the -db and -ix parameters. Input is the fasta sequence of the above given proteins. Use SwissProt (/mnt/project/pracstrucfunc12/data/swissprot/uniprot_sprot
) as database and /mnt/project/pracstrucfunc12/data/index_pp/uniprot_sprot.idx
as index.
Use the blastget output to create a MSA using Kalign (/mnt/opt/T-Coffee/bin/kalign
).
Use the MSA as input for PolyPhobius (/mnt/project/pracstrucfunc12/polyphobius/jphobius
). Usage: jphobius -h
. Do not forget the -poly parameter.
Compare the results to the membrane assignment of the structures for these proteins in OPM and/or PDBTM.
Signal peptides
Use SignalP to predict signal peptides for the following proteins:
- P02768
- P47863
- P11279
You can run signalp
in the student computer pool. This is right now version 3.0. However, Tim will update to version 4.0 soon. Please note which version you are using. You can also use the SignalP server. This is version 4.0. For one of the example proteins I got a different prediction from SignalP 3.0 and 4.0.
You can use the Signal Peptide Website to look up the proteins. Check also for predicted transmembrane helices using PolyPhobius.
GO terms
Use GOPET and ProtFun2.0 to predit GO terms for your protein. What can you learn about its function form sequence alone?
Use Pfam -> SEQUENCE SEARCH to find out more about the Pfam family of your protein.
Questions to answer
- 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, Pfam...
- Look for other methods to get an idea how many different are available. (You can also try out more methods.)
- What else could be predicted from protein sequence alone?