ShapeKnots is used to predict a set of structures that can contain pseudoknots. It can take SHAPE mapping data as a restraint to increase structure prediction accuracy.
USAGE: ShapeKnots <seq file> <ct file> [options]
<seq file> |
The name of a sequence file containing input data.
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. |
-h, -H, --help |
Display the usage details message. |
-v, -V, --version |
Display version and copyright information for this interface. |
-dsh, -DSH, --DSHAPE |
Specify a differential SHAPE data file to be used to generate restraints in addition to SHAPE restraints specified by --SHAPE. These restraints specifically use SHAPE pseudoenergy restraints where the offset is zero. Also, the function for calculating the pseudo free energy for nucleotide i is = (differential slope) * (differential SHAPE for nucleotide i). These pseudoenergies are added to those generated with the --SHAPE option.
Default is no differential SHAPE data file specified. |
-dsm, -DSM, --DSHAPEslope |
Specify a slope used with differential SHAPE restraints.
Default is 2.11 kcal/mol. |
-dso, -DSO, --doubleOffset |
Specify a double-stranded offset file, which adds energy bonuses to particular double-stranded nucleotides.
Default is to have no double-stranded offset specified.
|
-im, -IM, --IMaximum |
Specify a maximum number of internally generated structures for each call of the dynamic programming algorithm. Note that suboptimal structures are generated until either the maximum number of structures is reached or the maximum percent difference is reached (below). This is not the maximum number of structures provided to the user, which is controlled by –m, -M, --maximum.
Default is 100 structures. |
-ip, -IP, --IPercent |
Specify a maximum percent difference in folding free energy change for internally generated suboptimal structures for each call of the dynamic programming algorithm. For example, 20 would indicate 20%. This is not the maximum percent difference in energy for structures provided to the user, which is controlled by –p, -P, --percent.
Default is 20%. |
-iw, -IW, --IWindow |
Specify a window size for the internally generated suboptimal structures for each call of the dynamic programming algorithm. This is not the window for structures provided to the user, which is controlled by –w, -W, --window.
Default is determined by the length of the sequence. |
-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 20 structures.
|
-p, -P, --percent |
Specify a maximum percent difference in folding free energy change for generating suboptimal structures. For example, 20 would indicate 20%.
Default is determined by the length of the sequence.
|
-p1, -P1, --Penalty1 |
Specify a pseudoknot penalty 1.
Default is 0.35 kcal/mol. |
-p2, -P2, --Penalty2 |
Specify a pseudoknot penalty 2.
Default is 0.65 kcal/mol. |
-ph, -PH, --PseudoknottedHelices |
Specify maximum number of helices to be processed.
Default is 100 helices. |
-sh, -SH, --SHAPE |
Specify a SHAPE data file to be used to generate restraints. These restraints specificially use SHAPE pseudoenergy restraints.
Default is no SHAPE data file specified.
|
-si, -SI, --SHAPEintercept |
Specify an intercept used with SHAPE restraints.
Default is -0.6 kcal/mol.
|
-sm, -SM, --SHAPEslope |
Specify a slope used with SHAPE restraints.
Default is 1.8 kcal/mol.
|
-sso, -SSO, --singleOffset |
Specify a single-stranded offset file, which adds energy bonuses to particular single-stranded nucleotides.
Default is to have no single-stranded offset specified.
|
-w, -W, --window |
Specify a window size.
Default is determined by the length of the sequence.
|
- Hajdin, C.E., Bellaousov, S., Higgins, W., Leonard, C.W., Mathews, D.H. and Weeks, K.M.
“Accurate SHAPE-directed RNA secondary structure modeling, including pseudoknots.” Proc. Natl. Acad. Sci. U.S.A., In Press. (2013).
-
Reuter, J.S. and Mathews, D.H.
"RNAstructure: software for RNA secondary structure prediction and analysis."
BMC Bioinformatics, 11:129. (2010).
-
Deigan, K.E., Li, T.W., Mathews, D.H. and Weeks, K.M.
"Accurate SHAPE-directed RNA structure determination."
Proc. Natl. Acad. Sci. U.S.A., 106:97-102. (2009).
|