Difference between revisions of "Task 3 - Sequence-based predictions"
Line 14: | Line 14: | ||
== Secondary structure == |
== Secondary structure == |
||
− | Use ReProf to predict secondary structure for your protein. Apply ReProf also to these proteins: |
+ | Use ReProf to predict secondary structure for your protein. Apply ReProf also to these proteins (given are UniProt IDs): |
* P10775 |
* P10775 |
||
* Q9X0E6 |
* Q9X0E6 |
||
Line 32: | Line 32: | ||
* P35462 |
* P35462 |
||
+ | PolyPhobius is installed in <code>/mnt/project/pracstrucfunc12/polyphobius</code>. |
||
− | PolyPhobius is installed in <code>/mnt/project/pracstrucfunc12/polyphobius</code>. In contrast to its precursor Phobius, PolyPhobius uses homology informtion for the prediction. First, you have to execute a blast search. PolyPhobius distributed its own perl script for this purpose: <code>blastget</code> (<code>/mnt/project/pracstrucfunc12/polyphobius/blastget</code>). Usage: <code>blastget -h</code>. Use only the -db and -ix parameters. Input is the fasta sequence of the above given proteins. Use SwissProt (<code>/mnt/project/pracstrucfunc12/data/swissprot/uniprot_sprot</code>) as database and <code>/mnt/project/pracstrucfunc12/data/index_pp/uniprot_sprot.idx</code> as index. |
||
+ | |||
+ | 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: <code>blastget</code> (<code>/mnt/project/pracstrucfunc12/polyphobius/blastget</code>). Usage: <code>blastget -h</code>. Use only the -db and -ix parameters. Input is the fasta sequence of the above given proteins. Use SwissProt (<code>/mnt/project/pracstrucfunc12/data/swissprot/uniprot_sprot</code>) as database and <code>/mnt/project/pracstrucfunc12/data/index_pp/uniprot_sprot.idx</code> as index. |
||
+ | |||
+ | Use the blastget output to create a MSA using Kalign (<code>/mnt/opt/T-Coffee/bin//kalign</code>). |
||
+ | |||
Revision as of 03:15, 8 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 preditiction methods use evolutionary information, i.e. sequence similiarity. Sequence alignments are therefore often a prerequisite for the predictions.
Contents
Theoretical background talks
The introductory talks will give an introduction to sequence-based protein predictions. In particular:
- secondary structure
- disorder
- transmembrane helices
- GO terms
Where to run the jobs
- You can log in to the student computer pool:
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 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 README here:
/opt/iupred/README
.
Compare the results to the DisProt database.
Transmembrane helices
Use PolyPhobius to predict transmembrane helices for your protein and for the follwoing proteins:
- P35462
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
).
Give a brief description of the theory and the algorithm, if possible
What is predicted? Describe the features in some detail
What information is required for the predictions?
Apply the prediction methods to your protein and explain how to do this
You may try out other protein sequences (this is required when specified)
Present, describe and discuss the results
Look for other methods, for example here: http://expasy.org/tools/
You may try out more methods