|
DuplexFold is used to predict lowest free energy structures containing two strands. The structures will not contain intramolecular pairs.
Note that output is written to a CT file where the sequences are concatenated, with an intermolecular linker between them ("III").
USAGE: DuplexFold <seq file 1> <seq file 2> <ct file> [options]
| <seq file 1> |
The name of a sequence file containing input data for the first sequence.
Note that lowercase nucleotides are forced single-stranded in structure prediction.
|
| <seq file 2> |
The name of a sequence file containing input data for the second sequence.
Note that lowercase nucleotides are forced single-stranded in structure prediction.
|
| <ct file> |
The name of a CT file to which output will be written. |
| -d, -D, --DNA |
Specify that the sequence is DNA, and DNA parameters are to be used.
Note that the backbone type must be the same for both sequences.
Default is to use RNA parameters.
|
| -h, -H, --help |
Display the usage details message. |
| -a, -A, --alphabet |
Specify the name of a folding alphabet and associated
nearest neighbor parameters. The alphabet is the prefix
for the thermodynamic parameter files, e.g. "rna" for RNA
parameters or "dna" for DNA parameters or a custom
extended/modified alphabet. The thermodynamic parameters
need to reside in the at the location indicated by
environment variable DATAPATH.
The default is "rna" (i.e. use RNA parameters). This
option overrides the --DNA flag. |
| -l, -L, --loop |
Specify the maximum number of unpaired nucleotides in an internal or bulge loop.
Default is 6 unpaired nucleotides.
|
| -m, -M, --maximum |
Specify a maximum number of structures. Note that suboptimal structures are generated until either the maximum number of structures are reached or the maximum percent difference is reached (below).
Default is 10 structures.
|
| -p, -P, --percent |
Specify a maximum percent difference in folding free energy change in suboptimal structures. Note that suboptimal structures are generated until either the maximum percent free energy difference is reached or until the maximum number of structures is reached (above).
Default is 40 percent (specified as 40, not 0.4).
|
| -t, -T, --temperature |
Specify the temperature at which calculation takes place in Kelvin.
Default is 310.15 K, which is 37 degrees C.
|
| -w, -W, --window |
Specify a window size for generating suboptimal structures. Smaller windows (down to zero) will result in a larger set of structures that are similar. Larger windows result in fewer structures that are more different in predicted pairs.
Default is 0 nucleotides.
|
-
Reuter, J.S. and Mathews, D.H.
"RNAstructure: software for RNA secondary structure prediction and analysis."
BMC Bioinformatics, 11:129. (2010).
-
Piekna-Przybylska, D., DiChiacchio, L., Mathews, D.H. and Bambara, R.A.
"A sequence similar to tRNA3Lys gene is embedded in HIV-1 U3/R and promotes minus strand transfer."
Nature Structural & Molecular Biology, 17:83-89. (2010).
|
