comparison rnaup.xml @ 0:d14bc79748a3 draft

planemo upload for repository https://github.com/bgruening/galaxytools/tree/master/tools/rna_tools/vienna_rna commit 0065dafe7bbd382bb995b28cc4089c9e4f4eeeb9
author rnateam
date Tue, 06 Dec 2016 12:34:24 -0500
parents
children
comparison
equal deleted inserted replaced
-1:000000000000 0:d14bc79748a3
1 <tool id="viennarna_rnaup" name="@EXECUTABLE@" version="@VERSION@.0">
2 <description>Calculate the thermodynamics of RNA-RNA interactions</description>
3 <macros>
4 <token name="@EXECUTABLE@">RNAup</token>
5 <import>macros.xml</import>
6 </macros>
7 <expand macro="requirements" />
8 <expand macro="stdio" />
9 <expand macro="version_command" />
10 <command>
11 <![CDATA[
12 RNAup < '$input' > '$output'
13 -T$temperature
14 --dangles=$dangling
15 $constraint
16 --ulength=$ulength
17 --contributions=#echo ''.join(str($contributions).split(','))#
18 --window=$window
19 $includeboth
20 $interactionpairwise
21 $interactionfirst
22 --extend5=$extendfive
23 --extend3=$extendthree
24 #if $varExists('$advancedOptions.nolp')
25 --pfScale=$advancedOptions.pfscale
26 $advancedOptions.noconversion
27 $advancedOptions.nolp
28 $advancedOptions.nogu
29 $advancedOptions.noclosinggu
30 $advancedOptions.notetra
31 #end if
32 && tar -cf '$accesibilitiesFile' *.out
33 ]]>
34 </command>
35
36 <inputs>
37 <param format="fasta" name="input" type="data" label="Fasta file"/>
38 <param name="temperature" type="float" value="37.0" label="temperature [°C]" help="-T"/>
39 <param name="dangling" type="select" label="how to treat dangling end energies" help="-d">
40 <option value="2" selected="true">unpaired bases participate in all dangling ends (2)</option>
41 <option value="0">ignore dangling ends (0)</option>
42 <option value="1">unpaired bases participate in one dangling end only (1)</option>
43 <option value="3">allow coaxial stacking (3)</option>
44 </param>
45 <param name="constraint" type="boolean" truevalue="--constraint" falsevalue="" checked="false" label="Constraints for the secondary structure" help="--constraint"/>
46 <param name="ulength" type="integer" min="0" value="4" label="Length of the unstructured region." help="--ulength"/>
47 <param name="contributions" type="select" multiple="true" display="checkboxes" label="Specify which contributions to the probabilty of being unpaired are listed in the output" argument="--contributions">
48 <option value="S" selected="true">sum of all</option>
49 <option value="H">unpaired within hairpin loop</option>
50 <option value="I">unpaired within interior loop</option>
51 <option value="M">unpaired within multiloop</option>
52 <option value="E">unpaired within exterior loop</option>
53 <validator type="no_options" message="Please select at least one contribution."/>
54 </param>
55 <param name="window" type="integer" min="0" value="25" label="Maximal length of the region of interaction." help="--window"/>
56 <param name="includeboth" type="boolean" truevalue="--include_both" falsevalue="" checked="false" label="Include the probability of unpaired regions in both RNAs" help="--include_both"/>
57 <param name="interactionpairwise" type="boolean" truevalue="--interaction_pairwise" falsevalue="" checked="false" label="Activate pairwise interaction mode" help="--interaction_pairwise"/>
58 <param name="interactionfirst" type="boolean" truevalue="--interaction_first" falsevalue="" checked="false" label="Activate pairwise interaction mode, compare always with the first sequence" help="--interaction_first"/>
59 <param name="extendfive" type="integer" min="0" value="0" label="Extend the region of interaction in the target to some residues on the 5' side." help="--extend5"/>
60 <param name="extendthree" type="integer" min="0" value="0" label="Extend the region of interaction in the target to some residues on the 3' side." help="--extend3"/>
61 <conditional name="advancedOptions">
62 <param name="advancedSelector" type="select" label="advanced options">
63 <option value="basic">basic Options</option>
64 <option value="advanced">advanced Options</option>
65 </param>
66 <when value="advanced">
67 <param name="pfscale" type="float" value="1.07" label="Use scale*mfe as an estimate for the ensemble free energy " help="--pfScale"/>
68 <param name="nolp" type="boolean" truevalue="" falsevalue="--noLP" checked="true" label="Allow lonely base-pairs" help="(--noLP)"/>
69 <param name="nogu" type="boolean" truevalue="" falsevalue="--noGU" checked="true" label="Allow GU pairing" help="--noGU"/>
70 <param name="noclosinggu" type="boolean" truevalue="" falsevalue="--noClosingGU" checked="true" label="Allow GU pairing at the ends" help="Allow pairing of G and U at the ends of helices. --noClosingGU"/>
71 <param name="notetra" type="boolean" truevalue="" falsevalue="--noTetra" checked="true" label="Allow stabilization for loops, hairpins etc." help=" Include special tabulated stabilizing energies for tri-, tetra- and hexaloop hairpins. Mostly for testing. (--noTetra)"/>
72 <param name="noconversion" type="boolean" truevalue="" falsevalue="--noconv" checked="true" label="Convert Thymine to Uracil (T -> U)" help="Avoids confusion with DNA sequences (--noconv)"/>
73 </when>
74 <when value="basic">
75 </when>
76 </conditional>
77 </inputs>
78 <outputs>
79 <data format="txt" name="output"/>
80 <data format="tar" name="accesibilitiesFile"/>
81 </outputs>
82 <tests>
83 <test>
84 <param name="input" value="rnaup_input1.fa"/>
85 <output name="output" file="rnaup_result1.txt"/>
86 </test>
87 </tests>
88 <help>
89 <![CDATA[
90
91 **RNAup**
92
93 RNAup calculates the thermodynamics of RNA−RNA interactions, by decomposing the binding into two stages. (1) First the probability that a potential binding sites remains unpaired (equivalent to the free energy needed to open the site) is computed. (2) Then this accessibility is combined with the interaction energy to obtain the total binding energy. All calculations are done by computing partition functions over all possible conformations.
94
95 RNAup provides two different modes: By default RNAup computes accessibilities, in terms of the free energies needed to open a region (default length 4). It prints the region of highest accessibility and its opening energy to stdout, opening energies for all other regions are written to a file.
96
97 In interaction mode the interaction between two RNAs is calculated. It is invoked if the input consists of two sequences concatenated with an "&". RNAup assumes that the longer RNA is a structured target sequence while the shorter one is an unstructured small RNA. Additionally, for every position along the target sequence the best free energy of binding for an interaction that includes this position is written to the the output file. Output to stdout consists of the location and free energy, dG, for the optimal region of interaction. The binding energy dG is also split into its components the interaction energy dGint and the opening energy dGu_l (and possibly dGu_s for the shorter sequence).
98 In addition we print the optimal interaction structure as computed by RNAduplex for this region. Note that it can happen that the RNAduplex computed optimal interaction does not coincide with the optimal RNAup region. If the two predictions don’t match the structure string is replaced by a run of ".".
99
100 -----
101
102 **Input format**
103
104 RNAup requires one input file
105
106 - Fasta file
107
108 For the interaction mode, the sequences of the query and target have to be concatenated by an '&', or the --interaction-pairwise or --interaction-first flags have to be set.
109
110 ------
111
112 **Outputs**
113
114 - most accesible region of each sequence
115 - accesibilties of each sequence is written in a file and these files are bundled together in a tar file
116
117
118 ]]>
119 </help>
120 <expand macro="citations" />
121 </tool>