1
# $Id: rnaheat.pm,v 1.6 2006/07/04 22:23:35 mauricio Exp $
2
# BioPerl module for Bio::Tools::Run::PiseApplication::rnaheat
4
# Cared for by Catherine Letondal <letondal@pasteur.fr>
6
# For copyright and disclaimer see below.
8
# POD documentation - main docs before the code
12
Bio::Tools::Run::PiseApplication::rnaheat
20
Bio::Tools::Run::PiseApplication::rnaheat
24
VIENNARNA RNAheat - calculate specific heat of RNAs (Hofacker, Stadler)
28
I.L. Hofacker, W. Fontana, P.F. Stadler, S. Bonhoeffer, M. Tacker, P. Schuster (1994) Fast Folding and Comparison of RNA Secondary Structures. Monatshefte f. Chemie 125: 167-188
30
J.S. McCaskill (1990) The equilibrium partition function and base pair binding probabilities for RNA secondary structures, Biopolymers 29: 11051119 D. Adams (1979) The hitchhiker's guide to the galaxy, Pan Books, London
37
http://bioweb.pasteur.fr/seqanal/interfaces/rnaheat.html
38
for available values):
47
Lowest temperature, default is 0C (-Tmin)
50
Highest temperature, default is 100C (-Tmax)
53
Calculate partition function every stepsize degrees C. Default is 1C (-h)
56
Produces a smoother curve by increasing ipoints (-m)
59
Do not include special stabilizing energies for certain tetraloops (-4)
62
How to treat dangling end energies for bases adjacent to helices in free ends and multiloops (-d)
65
Do not allow GU pairs (-noGU)
68
Do not allow GU pairs at the end of helices (-noCloseGU)
71
Non standard pairs (comma seperated list) (-nsp)
77
Energy parameters for the artificial ABCD... alphabet (-e)
85
User feedback is an integral part of the evolution of this and other
86
Bioperl modules. Send your comments and suggestions preferably to
87
the Bioperl mailing list. Your participation is much appreciated.
89
bioperl-l@bioperl.org - General discussion
90
http://bioperl.org/wiki/Mailing_lists - About the mailing lists
94
Report bugs to the Bioperl bug tracking system to help us keep track
95
of the bugs and their resolution. Bug reports can be submitted via the
98
http://bugzilla.open-bio.org/
102
Catherine Letondal (letondal@pasteur.fr)
106
Copyright (C) 2003 Institut Pasteur & Catherine Letondal.
109
This module is free software; you can redistribute it and/or modify
110
it under the same terms as Perl itself.
114
This software is provided "as is" without warranty of any kind.
122
http://bioweb.pasteur.fr/seqanal/interfaces/rnaheat.html
126
Bio::Tools::Run::PiseApplication
130
Bio::Tools::Run::AnalysisFactory::Pise
134
Bio::Tools::Run::PiseJob
141
package Bio::Tools::Run::PiseApplication::rnaheat;
145
use Bio::Tools::Run::PiseApplication;
147
@ISA = qw(Bio::Tools::Run::PiseApplication);
152
Usage : my $rnaheat = Bio::Tools::Run::PiseApplication::rnaheat->new($location, $email, @params);
153
Function: Creates a Bio::Tools::Run::PiseApplication::rnaheat object.
154
This method should not be used directly, but rather by
155
a Bio::Tools::Run::AnalysisFactory::Pise instance.
156
my $factory = Bio::Tools::Run::AnalysisFactory::Pise->new();
157
my $rnaheat = $factory->program('rnaheat');
159
Returns : An instance of Bio::Tools::Run::PiseApplication::rnaheat.
164
my ($class, $location, $email, @params) = @_;
165
my $self = $class->SUPER::new($location, $email);
167
# -- begin of definitions extracted from /local/gensoft/lib/Pise/5.a/PerlDef/rnaheat.pm
169
$self->{COMMAND} = "rnaheat";
170
$self->{VERSION} = "5.a";
171
$self->{TITLE} = "VIENNARNA";
173
$self->{DESCRIPTION} = "RNAheat - calculate specific heat of RNAs";
175
$self->{OPT_EMAIL} = 0;
177
$self->{AUTHORS} = "Hofacker, Stadler";
179
$self->{DOCLINK} = "http://bioweb.pasteur.fr/docs/gensoft-na.html#VIENNARNA";
181
$self->{REFERENCE} = [
183
"I.L. Hofacker, W. Fontana, P.F. Stadler, S. Bonhoeffer, M. Tacker, P. Schuster (1994) Fast Folding and Comparison of RNA Secondary Structures. Monatshefte f. Chemie 125: 167-188",
185
"J.S. McCaskill (1990) The equilibrium partition function and base pair binding probabilities for RNA secondary structures, Biopolymers 29: 11051119 D. Adams (1979) The hitchhiker's guide to the galaxy, Pan Books, London",
188
$self->{_INTERFACE_STANDOUT} = undef;
189
$self->{_STANDOUT_FILE} = undef;
191
$self->{TOP_PARAMETERS} = [
200
$self->{PARAMETERS_ORDER} = [
202
"seq", # RNA Sequences File
203
"control", # Control options
204
"temp_min", # Lowest temperature, default is 0C (-Tmin)
205
"temp_max", # Highest temperature, default is 100C (-Tmax)
206
"stepsize", # Calculate partition function every stepsize degrees C. Default is 1C (-h)
207
"ipoints", # Produces a smoother curve by increasing ipoints (-m)
208
"tetraloops", # Do not include special stabilizing energies for certain tetraloops (-4)
209
"dangling", # How to treat dangling end energies for bases adjacent to helices in free ends and multiloops (-d)
210
"input", # Input parameters
211
"noGU", # Do not allow GU pairs (-noGU)
212
"noCloseGU", # Do not allow GU pairs at the end of helices (-noCloseGU)
213
"nsp", # Non standard pairs (comma seperated list) (-nsp)
214
"parameter", # Parameter file (-P)
215
"energy", # Energy parameters for the artificial ABCD... alphabet (-e)
221
"rnaheat" => 'String',
223
"control" => 'Paragraph',
224
"temp_min" => 'Integer',
225
"temp_max" => 'Integer',
226
"stepsize" => 'Integer',
227
"ipoints" => 'Integer',
228
"tetraloops" => 'Switch',
229
"dangling" => 'Excl',
230
"input" => 'Paragraph',
232
"noCloseGU" => 'Switch',
234
"parameter" => 'InFile',
236
"readseq" => 'String',
242
"perl" => '"RNAheat"',
245
"perl" => '" < $value" ',
250
"perl" => '(defined $value && $value ne $vdef)? " -Tmin $value" : ""',
253
"perl" => '(defined $value && $value ne $vdef)? " -Tmax $value" : ""',
256
"perl" => '(defined $value && $value ne $vdef)? " -h $value" : ""',
259
"perl" => '(defined $value && $value ne $vdef)? " -m $value" : ""',
262
"perl" => '($value)? " -4" : ""',
265
"perl" => '($value)? " $value" : ""',
270
"perl" => '($value)? " -noGU" : ""',
273
"perl" => '($value)? " -noCloseGU" : ""',
276
"perl" => '($value)? " -nsp $value" : "" ',
279
"perl" => '($value)? " -P $value" : ""',
282
"perl" => '($value)? " -e $value" : ""',
285
"perl" => '"/local/gensoft/lib/ViennaRNA/readseq -f=19 -a $seq > $seq.tmp && (cp $seq $seq.orig && mv $seq.tmp $seq) ; "',
290
$self->{FILENAMES} = {
319
$self->{BY_GROUP_PARAMETERS} = [
343
$self->{ISHIDDEN} = {
363
$self->{ISCOMMAND} = {
383
$self->{ISMANDATORY} = {
405
"seq" => "RNA Sequences File",
406
"control" => "Control options",
407
"temp_min" => "Lowest temperature, default is 0C (-Tmin)",
408
"temp_max" => "Highest temperature, default is 100C (-Tmax)",
409
"stepsize" => "Calculate partition function every stepsize degrees C. Default is 1C (-h)",
410
"ipoints" => "Produces a smoother curve by increasing ipoints (-m)",
411
"tetraloops" => "Do not include special stabilizing energies for certain tetraloops (-4)",
412
"dangling" => "How to treat dangling end energies for bases adjacent to helices in free ends and multiloops (-d)",
413
"input" => "Input parameters",
414
"noGU" => "Do not allow GU pairs (-noGU)",
415
"noCloseGU" => "Do not allow GU pairs at the end of helices (-noCloseGU)",
416
"nsp" => "Non standard pairs (comma seperated list) (-nsp)",
417
"parameter" => "Parameter file (-P)",
418
"energy" => "Energy parameters for the artificial ABCD... alphabet (-e)",
423
$self->{ISSTANDOUT} = {
445
"control" => ['temp_min','temp_max','stepsize','ipoints','tetraloops','dangling',],
446
"dangling" => ['','only unpaired bases can participate in at most one dangling end','-d','-d: ignores dangling ends altogether','-d2','-d2: the check is ignored, this is the default for partition function folding.',],
447
"input" => ['noGU','noCloseGU','nsp','parameter','energy',],
448
"energy" => ['1','1: use energy parameters for GC pairs','2','2: use energy parameters for AU pairs',],
455
$self->{SEPARATOR} = {
472
"rnaheat" => { "perl" => '1' },
473
"seq" => { "perl" => '1' },
474
"control" => { "perl" => '1' },
475
"temp_min" => { "perl" => '1' },
476
"temp_max" => { "perl" => '1' },
477
"stepsize" => { "perl" => '1' },
478
"ipoints" => { "perl" => '1' },
479
"tetraloops" => { "perl" => '1' },
480
"dangling" => { "perl" => '1' },
481
"input" => { "perl" => '1' },
482
"noGU" => { "perl" => '1' },
483
"noCloseGU" => { "perl" => '1' },
484
"nsp" => { "perl" => '1' },
485
"parameter" => { "perl" => '1' },
486
"energy" => { "perl" => '1' },
488
"perl" => 'defined $rnafold || defined $rnasubopt',
496
'(! (defined $rnafold || defined $rnaeval || defined $rnainverse) && ($dangling eq "-d2") && ($dangling = "") && 0)' => "no message",
506
$self->{WITHPIPEOUT} = {
514
$self->{WITHPIPEIN} = {
538
$self->{ISSIMPLE} = {
558
$self->{PARAMFILE} = {
564
"The program fits a parabola to 2*ipoints+1 data points to calculate 2nd derivatives. Increasing this parameter produces a smoother curve. Default is 2.",
567
"How to treat \'dangling end\' energies for bases adjacent to helices in free ends and multiloops: Normally only unpaired bases can participate in at most one dangling end. With -d2 this check is ignored, this is the default for partition function folding (-p). -d ignores dangling ends altogether. Note that by default pf and mfe folding treat dangling ends differently, use -d2 (or -d) in addition to -p to ensure that both algorithms use the same energy model. The -d2 options is available for RNAfold, RNAeval, and RNAinverse only.",
570
"Allow other pairs in addition to the usual AU,GC,and GU pairs. pairs is a comma seperated list of additionally allowed pairs. If a the first character is a \'-\' then AB will imply that AB and BA are allowed pairs. e.g. RNAfold -nsp -GA will allow GA and AG pairs. Nonstandard pairs are given 0 stacking energy.",
573
"Read energy parameters from paramfile, instead of using the default parameter set. A sample parameterfile should accompany your distribution. See the RNAlib documentation for details on the file format.",
578
$self->{SCALEMIN} = {
582
$self->{SCALEMAX} = {
586
$self->{SCALEINC} = {
594
# -- end of definitions extracted from /local/gensoft/lib/Pise/5.a/PerlDef/rnaheat.pm
598
$self->_init_params(@params);
605
1; # Needed to keep compiler happy