0
+ − 1 <tool id="twistdna" name="TwistDNA" version="0.1.0">
+ − 2 <description>Twist-DNA allows visualization of results in standard genome browsers to compare DNA opening properties to other available datasets</description>
+ − 3 <requirements>
+ − 4 <requirement type="package" version="1.1">twistdna</requirement>
+ − 5 </requirements>
+ − 6 <command>
+ − 7 cp $configfile "input.dat";
+ − 8 mkdir param;
+ − 9 mkdir outputs;
+ − 10 cp $enthalpy "param/param_enthalpy.dat";
+ − 11 cp $entropy "param/param_entropy.dat";
+ − 12 TwistDNA < $input;
+ − 13 ##remove the first line, because there they are useless
+ − 14 tail -n+2 outputs/bubblethresh.bed | tr -s ' ' '\t' > $outfile1;
+ − 15 tail -n+2 outputs/openproba.bed | tr -s ' ' '\t' > $outfile2;
+ − 16
+ − 17 </command>
+ − 18 <stdio>
+ − 19 <exit_code range="1:" level="fatal" description="Error occurred. Please check Tool Standard Error" />
+ − 20 <exit_code range=":-1" level="fatal" description="Error occurred. Please check Tool Standard Error" />
+ − 21 </stdio>
+ − 22 <inputs>
+ − 23 <param format="fasta" name="input" type="data" label="Input Alignment File" />
+ − 24 <conditional name="check">
+ − 25 <param name="advanced" type="select" label="Show advanced Options" help="" checked="false">
+ − 26 <option value="yes">Yes</option>
+ − 27 <option value="no">No</option>
+ − 28 </param>
+ − 29 <when value="yes">
+ − 30 <param name="temperature" label="Temperature in Celsius" type="float" value="37" />
+ − 31 <param name="sigma" label="Superhelical density Sigma" type="float" value="-0.06" />
+ − 32 <param name="salt" label="Salt concentration CNa" type="float" value="0.1" />
+ − 33 <param name="small_bubble" label="Smallest bubble size bsizemin" type="integer" value="1" />
+ − 34 <param name="large_bubble" label="Largest bubble size bsizemax" type="integer" value="0" />
+ − 35 <param name="bubble_size_step" label="Bubble size step bstep" type="integer" value="1" />
+ − 36 <param name="threshold" label="Bubble storing threshold" type="float" value="-3.0" />
+ − 37 </when>
+ − 38 </conditional>
+ − 39 </inputs>
+ − 40 <configfiles>
+ − 41 <configfile name="configfile">
+ − 42 $getVar('check.temperature', '37') ::Temperature T in Celsius (default 37 C)
+ − 43 $getVar('check.sigma', '-0.06') ::Superhelical density Sigma (default -0.06)
+ − 44 $getVar('check.salt', '0.1') ::Salt concentration CNa (default 0.1 M)
+ − 45 $getVar('check.small_bubble', '1') ::Smallest bubble size bsizemin
+ − 46 $getVar('check.large_bubble', '0') ::Largest bubble size bsizemax
+ − 47 $getVar('check.bubble_size_step', '1') ::Bubble size step bstep
+ − 48 $getVar('check.threshold', '-3.0') ::Bubble storing threshold (default -3)
+ − 49 </configfile>
+ − 50 <configfile name="enthalpy">
+ − 51 -12.84 -12.52 -11.58 -13.72
+ − 52 -13.43 -13.00 -13.72 -15.89
+ − 53 -11.72 -13.67 -12.84 -13.43
+ − 54 -13.67 -17.07 -12.52 -13.00
+ − 55 </configfile>
+ − 56 <configfile name="entropy">
+ − 57 -11.29 -10.53 -10.16 -11.56
+ − 58 -11.38 -9.89 -11.56 -12.42
+ − 59 -10.89 -11.55 -11.29 -11.38
+ − 60 -11.55 -13.68 -10.53 -9.89
+ − 61 </configfile>
+ − 62 </configfiles>
+ − 63 <outputs>
+ − 64 <data name="outfile1" format="bed" />
+ − 65 <data name="outfile2" format="bed" />
+ − 66 </outputs>
+ − 67 <help>**What it does**
+ − 68 http://www.cbp.ens-lyon.fr/doku.php?id=developpement:productions:logiciels:twistdna
+ − 69
+ − 70 Local opening of the DNA double-helix is required in many fundamental biological processes and is in part controlled by the degree of superhelicity imposed in vivo by the protein machinery. In particular, positions of superhelically destabilized regions correlate with regulatory sites along the genome. Based on a self-consistent linearization of a thermodynamic model of superhelical DNA introduced by Benham, we have developed Twist-DNA, a program that predicts the locations of these regions by efficiently computing base-pair and bubble opening probabilities in genomic DNA. The program allows visualization of results in standard genome browsers to compare DNA opening properties to other available datasets.
+ − 71 </help>
+ − 72 </tool>
