|
1
|
1 <tool id="antarna" name="antaRNA" version="1.1">
|
|
|
2 <description>
|
|
|
3 Ant Colony Optimized RNA Sequence Design
|
|
|
4 </description>
|
|
|
5 <requirements>
|
|
|
6 <requirement type="package" version"2.2.12">pkiss</requirement>
|
|
|
7 <requirement type="package" version"3.2.5">rnashapes</requirement>
|
|
|
8 <requirement type="package" version"1.9">numpy</requirement>
|
|
2
|
9 <requirement type="package" version"1.1">antarna</requirement>
|
|
1
|
10 </requirements>
|
|
|
11 <stdio>
|
|
|
12 <exit_code level="fatal" range="1:"/>
|
|
|
13 </stdio>
|
|
|
14 <version_command>python antaRNA.py --version</version_command>
|
|
|
15 <command interpreter="python"><![CDATA[antaRNA.py
|
|
|
16 #if $Cstr and $Cstr is not None:
|
|
|
17 -Cstr "$Cstr"
|
|
|
18 #end if
|
|
|
19
|
|
|
20 #if $tGC and $tGC is not None:
|
|
|
21 -tGC $tGC
|
|
|
22 #end if
|
|
|
23
|
|
|
24 #if $n and $n is not None:
|
|
|
25 -n $n
|
|
|
26 #end if
|
|
|
27 $GU
|
|
|
28
|
|
|
29 #if $s and $s is not None:
|
|
|
30 -s $s
|
|
|
31 #end if
|
|
|
32
|
|
|
33 #if $ip and $ip is not None:
|
|
|
34 -ip $ip
|
|
|
35 #end if
|
|
|
36
|
|
|
37 #if $r and $r is not None:
|
|
|
38 -r $r
|
|
|
39 #end if
|
|
|
40
|
|
|
41 #if $CC and $CC is not None:
|
|
|
42 -CC $CC
|
|
|
43 #end if
|
|
|
44
|
|
|
45 #if $aTC and $aTC is not None:
|
|
|
46 -aTC $aTC
|
|
|
47 #end if
|
|
|
48 $p
|
|
|
49 $pkPar
|
|
|
50
|
|
|
51 #if $strategy and $strategy is not None:
|
|
|
52 --strategy $strategy
|
|
|
53 #end if
|
|
|
54
|
|
|
55 #if $tGCmax and $tGCmax is not None:
|
|
|
56 -tGCmax $tGCmax
|
|
|
57 #end if
|
|
|
58
|
|
|
59 #if $tGCvar and $tGCvar is not None:
|
|
|
60 -tGCvar $tGCvar
|
|
|
61 #end if
|
|
|
62
|
|
|
63 #if $t and $t is not None:
|
|
|
64 -t $t
|
|
|
65 #end if
|
|
|
66
|
|
|
67 #if $P and $P is not None:
|
|
|
68 -P $P
|
|
|
69 #end if
|
|
|
70
|
|
|
71 #if $of and $of is not None:
|
|
|
72 -of $of
|
|
|
73 #end if
|
|
|
74
|
|
|
75 #if $Cseq and $Cseq is not None:
|
|
|
76 -Cseq $Cseq
|
|
|
77 #end if
|
|
|
78
|
|
|
79 #if $l and $l is not None:
|
|
|
80 -l $l
|
|
|
81 #end if
|
|
|
82
|
|
|
83 #if $name and $name is not None:
|
|
|
84 --name $name
|
|
|
85 #end if
|
|
|
86
|
|
|
87 #if $a and $a is not None:
|
|
|
88 -a $a
|
|
|
89 #end if
|
|
|
90
|
|
|
91 #if $b and $b is not None:
|
|
|
92 -b $b
|
|
|
93 #end if
|
|
|
94
|
|
|
95 #if $er and $er is not None:
|
|
|
96 -er $er
|
|
|
97 #end if
|
|
|
98
|
|
|
99 #if $Cstrw and $Cstrw is not None:
|
|
|
100 -Cstrw $Cstrw
|
|
|
101 #end if
|
|
|
102
|
|
|
103 #if $Cgcw and $Cgcw is not None:
|
|
|
104 -Cgcw $Cgcw
|
|
|
105 #end if
|
|
|
106
|
|
|
107 #if $Cseqw and $Cseqw is not None:
|
|
|
108 -Cseqw $Cseqw
|
|
|
109 #end if
|
|
|
110 $v
|
|
|
111 $ov
|
|
|
112 > $default]]></command>
|
|
|
113 <inputs>
|
|
|
114 <param label="Structure constraint using RNA dotbracket notation with fuzzy block constraint. (-Cstr)" name="Cstr" type="text"/>
|
|
|
115 <param label="Objective target GC content in [0,1]. (-tGC)" name="tGC" type="float" value="0"/>
|
|
|
116 <param label="Number of sequences which shall be produced. (-n)" name="n" type="integer" value="0"/>
|
|
|
117 <param checked="false" label="Allowing GU base pairs. (-GU)" name="GU" type="boolean" truevalue="-GU" falsevalue=""/>
|
|
|
118 <param default="none" label="Provides a seed value for the used pseudo random number generator. (-s)" name="s" type="text"/>
|
|
|
119 <param default="s" label="Select the improving method. h=hierarchical, s=score_based. (-ip)" name="ip" type="text"/>
|
|
|
120 <param label="Amount of maximal terrain resets, until the best solution is retuned as solution. (-r)" name="r" type="integer" value="5"/>
|
|
|
121 <param label="Delimits the convergence count criterion for a reset. (-CC)" name="CC" type="integer" value="130"/>
|
|
|
122 <param label="Delimits the amount of internal ants for termination convergence criterion for a reset. (-aTC)" name="aTC" type="integer" value="50"/>
|
|
|
123 <param checked="false" label="Switch to pseudoknot based prediction using pKiss. Check the pseudoknot parameter usage!!! (-p)" name="p" type="boolean" truevalue="-p" falsevalue=""/>
|
|
|
124 <param checked="false" label="Enable optimized parameters for the usage of pseudo knots (Further parameter input ignored). (-pkPar)" name="pkPar" type="boolean" truevalue="-pkPar" falsevalue=""/>
|
|
|
125 <param default="A" label="Defining the pKiss folding strategy. (--strategy)" name="strategy" type="text"/>
|
|
|
126 <param label="Provides a maximum tGC value [0,1] for the case of uniform distribution sampling. The regular tGC value serves as minimum value. (-tGCmax)" name="tGCmax" type="float" value="-1.0"/>
|
|
|
127 <param label="Provides a tGC variance (sigma square) for the case of normal distribution sampling. The regular tGC value serves as expectation value (mu). (-tGCvar)" name="tGCvar" type="float" value="-1.0"/>
|
|
|
128 <param label="Provides a temperature for the folding algorithms. (-t)" name="t" type="float" value="37.0"/>
|
|
|
129 <param default="" label="Changes the energy parameterfile of RNAfold. If using this explicitly, please provide a suitable energy file delivered by RNAfold. (-P)" name="P" type="text"/>
|
|
|
130 <param default="STDOUT" label="Provide a path and an output file, e.g. "/path/to/the/target_file". (-of)" name="of" type="text"/>
|
|
|
131 <param default="" label="Sequence constraint using RNA nucleotide alphabet {A,C,G,U} and wild-card "N". (-Cseq)" name="Cseq" type="text"/>
|
|
|
132 <param label="Sets the level of allowed influence of sequence constraint on the structure constraint [0:no influence; 3:extensive influence]. (-l)" name="l" type="integer" value="1"/>
|
|
|
133 <param default="antaRNA_" label="Defines a name which is used in the sequence output. (--name)" name="name" type="text"/>
|
|
|
134 <param label="Sets alpha, probability weight for terrain path influence. [0,1] (-a)" name="a" type="float" value="1.0"/>
|
|
|
135 <param label="Sets beta, probability weight for terrain pheromone influence. [0,1] (-b)" name="b" type="float" value="1.0"/>
|
|
|
136 <param label="Pheromone evaporation rate. (-er)" name="er" type="float" value="0.2"/>
|
|
|
137 <param label="Structure constraint quality weighting factor. [0,1] (-Cstrw)" name="Cstrw" type="float" value="0.5"/>
|
|
|
138 <param label="GC content constraint quality weighting factor. [0,1] (-Cgcw)" name="Cgcw" type="float" value="5.0"/>
|
|
|
139 <param label="Sequence constraint quality weighting factor. [0,1] (-Cseqw)" name="Cseqw" type="float" value="1.0"/>
|
|
|
140 <param checked="false" label="Displayes intermediate output. (-v)" name="v" type="boolean" truevalue="-v" falsevalue=""/>
|
|
|
141 <param checked="false" label="Prints additional output to the headers of the produced sequences. (-ov)" name="ov" type="boolean" truevalue="-ov" falsevalue=""/>
|
|
|
142 </inputs>
|
|
|
143 <outputs>
|
|
|
144 <data format="txt" hidden="false" name="default"/>
|
|
|
145 </outputs>
|
|
|
146 <help><![CDATA[
|
|
|
147 .
|
|
|
148
|
|
|
149 ===========================
|
|
|
150 antaRNA - ant assembled RNA
|
|
|
151 ===========================
|
|
|
152
|
|
|
153 - antaRNA uses the VIENNNA RNA Package
|
|
|
154 - specifically it uses RNAfold and RNAdistance to calculate energies of and distances between secondary structures (version 2.1.x)
|
|
|
155 - for the parametrization of antaRNA the version 2.1.3 of the ViennaRNA package was used
|
|
|
156
|
|
|
157
|
|
|
158 - For questions and remarks please feel free to contact us at http://www.bioinf.uni-freiburg.de/
|
|
|
159
|
|
|
160 Example calls:
|
|
|
161
|
|
|
162 - python antaRNA.py --Cstr "...(((...)))..." --tGC 0.5 -n 2
|
|
|
163 - python antaRNA.py --Cstr ".........AAA(((...)))AAA........." --tGC 0.5 -n 10 --output_file /path/to/antaRNA_TESTRUN -ov
|
|
|
164 - python antaRNA.py --Cstr "BBBBB....AAA(((...)))AAA....BBBBB" --Cseq "NNNNANNNNNCNNNNNNNNNNNGNNNNNNUNNN" --tGC 0.5 -n 10
|
|
|
165 ]]></help>
|
|
|
166 </tool>
|
|
|
167
|
