|
1
|
1 <tool id="interproscan" name="Interproscan functional predictions of ORFs" version="1.2">
|
|
|
2 <description>Interproscan functional predictions of ORFs</description>
|
|
0
|
3 <command>
|
|
|
4 ## The command is a Cheetah template which allows some Python based syntax.
|
|
|
5 ## Lines starting hash hash are comments. Galaxy will turn newlines into spaces
|
|
|
6
|
|
|
7 ## create temp directory
|
|
|
8 #import tempfile, os
|
|
|
9 #set $tfile = tempfile.mkstemp()[1]
|
|
|
10
|
|
|
11 sed 's/ /_/g' $input > $tfile;
|
|
|
12
|
|
|
13 ## Hack, because interproscan does not seem to produce gff output even if it is configured
|
|
1
|
14 #if str($oformat) == "gff":
|
|
0
|
15 #set $tfile2 = tempfile.mkstemp()[1]
|
|
1
|
16 iprscan -cli -nocrc -i $tfile -o $tfile2 -goterms -seqtype p -altjobs -format raw -appl $appl 2>&1;
|
|
|
17 converter.pl -format gff3 -input $tfile2 -output $output;
|
|
|
18 rm $tfile2;
|
|
0
|
19 #else
|
|
1
|
20 iprscan -cli -nocrc -i $tfile -o $output -goterms -seqtype p -altjobs -format $oformat -appl $appl 2>&1;
|
|
0
|
21 #end if
|
|
|
22
|
|
1
|
23 rm $tfile
|
|
0
|
24
|
|
1
|
25 </command>
|
|
0
|
26 <inputs>
|
|
|
27 <param name="input" type="data" format="fasta" label="Protein Fasta File"/>
|
|
|
28
|
|
1
|
29 <param name="appl" type="select" format="text" label="Applications to run" help="Select your programm.">
|
|
0
|
30 <option value="blastprodom+coils+gene3d+hamap+hmmpanther+hmmpir+hmmpfam+hmmsmart+hmmtigr+fprintscan+patternscan+profilescan+superfamily+seg+signalp+tmhmm" selected="true">all</option>
|
|
|
31 <option value="seg">seg</option>
|
|
|
32 <option value="signalp">signalp</option>
|
|
|
33 <option value="tmhmm">tmhmm</option>
|
|
|
34 <option value="fprintscan">fprintscan</option>
|
|
|
35 <option value="patternscan">patternscan</option>
|
|
|
36 <option value="profilescan">profilescan</option>
|
|
|
37 <option value="superfamily">superfamily</option>
|
|
|
38 <option value="hmmpir">hmmpir</option>
|
|
|
39 <option value="hmmpfam">hmmpfam</option>
|
|
|
40 <option value="hmmsmart">hmmsmart</option>
|
|
|
41 <option value="hmmtigr">hmmtigr</option>
|
|
|
42 <option value="hmmpanther">hmmpanther</option>
|
|
|
43 <option value="hamap">hamap</option>
|
|
|
44 <option value="gene3d">gene3d</option>
|
|
|
45 <option value="coils">coils</option>
|
|
|
46 <option value="blastprodom">blastprodom</option>
|
|
1
|
47 </param>
|
|
0
|
48
|
|
|
49 <param name="oformat" type="select" label="Output format" help="Please select a output format.">
|
|
1
|
50 <option value="gff">gff</option>
|
|
|
51 <option value="raw" selected="true">raw</option>
|
|
0
|
52 <option value="txt">txt</option>
|
|
|
53 <option value="html">html</option>
|
|
|
54 <option value="xml">xml</option>
|
|
|
55 <option value="ebixml">EBI header on top of xml</option>
|
|
|
56 </param>
|
|
|
57
|
|
|
58 </inputs>
|
|
|
59 <outputs>
|
|
|
60
|
|
|
61 <data format="txt" name="output" label="Interproscan calculation on ${on_string}">
|
|
|
62 <change_format>
|
|
|
63 <when input="oformat" value="html" format="html"/>
|
|
|
64 <when input="oformat" value="xml" format="xml"/>
|
|
|
65 <when input="oformat" value="ebixml" format="xml"/>
|
|
|
66 <when input="oformat" value="gff" format="gff"/>
|
|
|
67 </change_format>
|
|
|
68 </data>
|
|
|
69
|
|
|
70 </outputs>
|
|
|
71 <requirements>
|
|
|
72 </requirements>
|
|
|
73 <help>
|
|
|
74 **What it does**
|
|
|
75
|
|
|
76 Interproscan is a batch tool to query the Interpro database. It provides annotations based on multiple searches of profile and other functional databases.
|
|
|
77 These include SCOP, CATH, PFAM and SUPERFAMILY.
|
|
|
78
|
|
|
79 **Input**
|
|
1
|
80
|
|
|
81 Required is a FASTA file containing ORF predictions. This file must NOT contain any spaces in the FASTA headers - any spaces will be convereted to underscores ``_`` by this tool before running with Interproscan.
|
|
0
|
82
|
|
|
83 **Output**
|
|
|
84
|
|
|
85 Example for the raw format.
|
|
1
|
86 The output will consist of a tabular file in galaxy with 14 columns and can be easily concatenated or filtered.
|
|
0
|
87
|
|
|
88 ====== ================================================================ ======================================================================
|
|
|
89 column example description
|
|
|
90 ====== ================================================================ ======================================================================
|
|
1
|
91 c1 NF00181542 the id of the input sequence
|
|
|
92 c2 27A9BBAC0587AB84 the crc64 (checksum) of the protein sequence (supposed to be unique)
|
|
|
93 c3 272 the length of the sequence (in AA)
|
|
0
|
94 c4 HMMPIR the anaysis method launched.
|
|
1
|
95 c5 PIRSF001424 the database members entry for this match
|
|
|
96 c6 Prephenate dehydratase the database member description for the entry
|
|
|
97 c7 1 the start of the domain match
|
|
|
98 c8 270 the end of the domain match
|
|
|
99 c9 6.5e-141 the evalue of the match (reported by member database method)
|
|
|
100 c10 T the status of the match (T: true, ?: unknown)
|
|
0
|
101 c11 06-Aug-2005 the date of the run.
|
|
1
|
102 c12 IPR008237 the corresponding InterPro entry (if iprlookup requested by the user)
|
|
|
103 c13 Prephenate dehydratase with ACT region the description of the InterPro entry
|
|
|
104 c14 Molecular Function:prephenate dehydratase activity (GO:0004664) the GO (gene ontology) description for the InterPro entry
|
|
0
|
105 ====== ================================================================ ======================================================================
|
|
|
106
|
|
|
107 **Database updates**
|
|
|
108
|
|
1
|
109 Typically these take place 2-3 times a year. Please contact your Galaxy administrator to update the databases.
|
|
0
|
110
|
|
1
|
111 -----
|
|
|
112 Tools
|
|
|
113 -----
|
|
0
|
114
|
|
1
|
115 **PROSITE patterns**::
|
|
0
|
116
|
|
|
117 Some biologically significant amino acid patterns can be summarised in
|
|
|
118 the form of regular expressions.
|
|
1
|
119 ScanRegExp (by Wolfgang.Fleischmann@ebi.ac.uk).
|
|
0
|
120
|
|
1
|
121 **PROSITE profiles**::
|
|
0
|
122
|
|
|
123 There are a number of protein families as well as functional or
|
|
|
124 structural domains that cannot be detected using patterns due to their extreme
|
|
|
125 sequence divergence, so the use of techniques based on weight matrices
|
|
|
126 (also known as profiles) allows the detection of such proteins or domains.
|
|
|
127 A profile is a table of position-specific amino acid weights and gap costs.
|
|
|
128 The profile structure used in PROSITE is similar to but slightly more general
|
|
1
|
129 (Bucher P. et al., 1996) than the one introduced by M. Gribskov and
|
|
0
|
130 co-workers.
|
|
|
131 pfscan from the Pftools package (by Philipp.Bucher@isrec.unil.ch).
|
|
|
132
|
|
1
|
133 **PRINTS**::
|
|
0
|
134
|
|
|
135 The PRINTS database houses a collection of protein family fingerprints.
|
|
|
136 These are groups of motifs that together are diagnostically more
|
|
|
137 powerful than single motifs by making use of the biological context inherent in a
|
|
|
138 multiple-motif method. The fingerprinting method arose from the need for
|
|
|
139 a reliable technique for detecting members of large, highly divergent
|
|
|
140 protein super-families.
|
|
1
|
141 FingerPRINTScan (Scordis P. et al., 1999).
|
|
0
|
142
|
|
1
|
143 **PFAM**::
|
|
0
|
144
|
|
|
145 Pfam is a database of protein domain families. Pfam contains curated
|
|
|
146 multiple sequence alignments for each family and corresponding hidden
|
|
1
|
147 Markov models (HMMs) (Eddy S.R., 1998).
|
|
0
|
148 Profile hidden Markov models are statistical models of the primary
|
|
|
149 structure consensus of a sequence family. The construction and use
|
|
|
150 of Pfam is tightly tied to the HMMER software package.
|
|
|
151 hmmpfam from the HMMER2.3.2 package (by Sean Eddy,
|
|
|
152 eddy@genetics.wustl.edu, http://hmmer.wustl.edu).
|
|
|
153
|
|
1
|
154 **PRODOM**::
|
|
0
|
155
|
|
|
156 ProDom is a database of protein domain families obtained by automated
|
|
|
157 analysis of the SWISS-PROT and TrEMBL protein sequences. It is useful
|
|
|
158 for analysing the domain arrangements of complex protein families and the
|
|
|
159 homology relationships in modular proteins. ProDom families are built by
|
|
|
160 an automated process based on a recursive use of PSI-BLAST homology
|
|
|
161 searches.
|
|
|
162 ProDomBlast3i.pl (by Emmanuel Courcelle emmanuel.courcelle@toulouse.inra.fr
|
|
|
163 and Yoann Beausse beausse@toulouse.inra.fr)
|
|
1
|
164 a wrapper on top of the Blast package (Altschul S.F. et al., 1997).
|
|
0
|
165
|
|
1
|
166 **SMART**::
|
|
0
|
167
|
|
|
168 SMART (a Simple Modular Architecture Research Tool) allows the
|
|
|
169 identification and annotation of genetically mobile domains and the
|
|
|
170 analysis of domain architectures. These domains are extensively
|
|
|
171 annotated with respect to phyletic distributions, functional class, tertiary
|
|
|
172 structures and functionally important residues. SMART alignments are
|
|
|
173 optimised manually and following construction of corresponding hidden Markov models (HMMs).
|
|
|
174 hmmpfam from the HMMER2.3.2 package (by Sean Eddy,
|
|
|
175 eddy@genetics.wustl.edu, http://hmmer.wustl.edu).
|
|
|
176
|
|
1
|
177 **TIGRFAMs**::
|
|
0
|
178
|
|
|
179 TIGRFAMs are a collection of protein families featuring curated multiple
|
|
|
180 sequence alignments, Hidden Markov Models (HMMs) and associated
|
|
|
181 information designed to support the automated functional identification
|
|
|
182 of proteins by sequence homology. Classification by equivalog family
|
|
|
183 (see below), where achievable, complements classification by orthologs,
|
|
|
184 superfamily, domain or motif. It provides the information best suited
|
|
|
185 for automatic assignment of specific functions to proteins from large
|
|
|
186 scale genome sequencing projects.
|
|
|
187 hmmpfam from the HMMER2.3.2 package (by Sean Eddy,
|
|
|
188 eddy@genetics.wustl.edu, http://hmmer.wustl.edu).
|
|
|
189
|
|
1
|
190 **PIR SuperFamily**::
|
|
0
|
191
|
|
|
192 PIR SuperFamily (PIRSF) is a classification system based on evolutionary
|
|
|
193 relationship of whole proteins.
|
|
|
194 hmmpfam from the HMMER2.3.2 package (by Sean Eddy,
|
|
|
195 eddy@genetics.wustl.edu, http://hmmer.wustl.edu).
|
|
|
196
|
|
1
|
197 **SUPERFAMILY**::
|
|
0
|
198
|
|
|
199 SUPERFAMILY is a library of profile hidden Markov models that represent
|
|
|
200 all proteins of known structure, based on SCOP.
|
|
|
201 hmmpfam/hmmsearch from the HMMER2.3.2 package (by Sean Eddy,
|
|
|
202 eddy@genetics.wustl.edu, http://hmmer.wustl.edu).
|
|
|
203 Optionally, predictions for coiled-coil, signal peptide cleavage sites
|
|
|
204 (SignalP v3) and TM helices (TMHMM v2) are supported (See the FAQs file
|
|
|
205 for details).
|
|
|
206
|
|
1
|
207 **GENE3D**::
|
|
0
|
208
|
|
1
|
209 Gene3D is supplementary to the CATH database. This protein sequence database
|
|
|
210 contains proteins from complete genomes which have been clustered into protein
|
|
|
211 families and annotated with CATH domains, Pfam domains and functional
|
|
|
212 information from KEGG, GO, COG, Affymetrix and STRINGS.
|
|
|
213 hmmpfam from the HMM2.3.2 package (by Sean Eddy,
|
|
|
214 eddy@genetics.wustl.edu, http://hmmer.wustl.edu).
|
|
0
|
215
|
|
1
|
216 **PANTHER**::
|
|
|
217
|
|
|
218 The PANTHER (Protein ANalysis THrough Evolutionary Relationships)
|
|
|
219 Classification System was designed to classify proteins (and their genes)
|
|
|
220 in order to facilitate high-throughput analysis.
|
|
|
221 hmmsearch from the HMM2.3.2 package (by Sean Eddy,
|
|
|
222 eddy@genetics.wustl.edu, http://hmmer.wustl.edu).
|
|
|
223 and blastall from the Blast package (Altschul S.F. et al., 1997).
|
|
0
|
224
|
|
1
|
225 ----------
|
|
|
226 References
|
|
|
227 ----------
|
|
|
228
|
|
|
229 Zdobnov EM, Apweiler R (2001)
|
|
|
230 InterProScan an integration platform for the signature-recognition methods in InterPro.
|
|
|
231 Bioinformatics 17, 847-848.
|
|
|
232 http://dx.doi.org/10.1093/bioinformatics/17.9.847
|
|
|
233
|
|
|
234 Quevillon E, Silventoinen V, Pillai S, Harte N, Mulder N, Apweiler R, Lopez R (2005)
|
|
|
235 InterProScan: protein domains identifier.
|
|
|
236 Nucleic Acids Research 33 (Web Server issue), W116-W120.
|
|
|
237 http://dx.doi.org/10.1093/nar/gki442
|
|
|
238
|
|
|
239 Hunter S, Apweiler R, Attwood TK, Bairoch A, Bateman A, Binns D, Bork P, Das U, Daugherty L, Duquenne L, Finn RD, Gough J, Haft D, Hulo N, Kahn D, Kelly E, Laugraud A, Letunic I, Lonsdale D, Lopez R, Madera M, Maslen J, McAnulla C, McDowall J, Mistry J, Mitchell A, Mulder N, Natale D, Orengo C, Quinn AF, Selengut JD, Sigrist CJ, Thimma M, Thomas PD, Valentin F, Wilson D, Wu CH, Yeats C. (2009)
|
|
|
240 InterPro: the integrative protein signature database.
|
|
|
241 Nucleic Acids Research 37 (Database Issue), D224-228.
|
|
|
242 http://dx.doi.org/10.1093/nar/gkn785
|
|
0
|
243
|
|
|
244
|
|
1
|
245 **Galaxy Wrapper Author**::
|
|
0
|
246
|
|
1
|
247 * Bjoern Gruening, Pharmaceutical Bioinformatics, University of Freiburg
|
|
|
248 * Konrad Paszkiewicz, Exeter Sequencing Service, University of Exeter
|
|
0
|
249
|
|
1
|
250 </help>
|
|
0
|
251 </tool>
|
