|
OligoWalk is used to predict the accessible regions in an RNA target to oligonucleotide hybridization. It calculates thermodynamic features of sense-antisense hybidization and predicts the free energy changes of oligonucleotides binding to a target RNA. The secondary structures of the oligomer and target mRNA are considered in the OligoWalk algorithm. It can be used to design efficient siRNA targeting a given mRNA sequence.
USAGE: OligoWalk <sequence file> <report file> [options]
Required parameters:
| <sequence file> |
The name of a sequence file (SEQ, FASTA) or structure file (CT, DBN) containing the input sequence. |
| <report file> |
The name of an output file to which report results will be written. |
| -d, --dna |
Specify that the oligomers are DNA. (The target sequence is still assumed to be RNA). This causes DNA:RNA hybridization energies to be used.
(In the absense of this flag, the default behavior is to assume both the oligomers and the target sequence are RNA) |
| -h, --help |
Display the usage details message. |
| -v, --version |
Display the RNAstructure package version information. |
| --html |
Write the report in HTML format instead of plain text. |
| --no-header |
The default behavior is to include a summary of the parameters and options used in the
calulation directly in the output report file. This option DISABLES that behavior. (I.e. the header information will not be included in the report file if this flag is present.) |
| -score |
Score the siRNA prefilter. |
| -w, --webserver |
Use special output for running OligoWalk as a webserver. This implies HTML=true and it sends the header (which summarizes the parameters and options used in the calulation) to STDOUT instead of including them in the report. It also turns off progress updates. |
| -write |
Write sav files to save time in test mode. |
| -c, -co, --conc |
Specify the molar concentration of oligos. e.g. "1.5E-6", "1.5uM", or "0.0000015" (1.5 micromolar).
This may be used in conjunction with the '--unit' flag. See '--unit for more examples.
Unit Abbreviations: mM=10-3, uM=10-6, nM=10-9, pM=10-12 |
| -u, --unit |
Specifies a power-of-ten factor to multiply the concentration by.
This parameter is not required, because it is entirely possible to specify the concentration directly in scientific notation (e.g. "--conc 4.2E-3" or "--conc 4.2mM". However there may be some circumstances when it is convenient to separate the concentration from the units, e.g. "--conc 4.2 --unit -3" or "--conc 4.2 --unit mM"
This might be useful, for example, when the program is run as a backend to a graphical or web interface and the user has selected the units separately from the concentration (from a drop-down list etc).
Examples:
- 2.5 millimolar can be specified in any of the following ways:
- --conc 2.5mM
- --conc 2.5E-3
- --conc 2.5 --unit -3
- --conc 2.5 --unit mM
- 5x10-5 M can be specified in any of the following ways:
- --conc 5E-5
- --conc 50uM
- --conc 5 --unit -5
- --conc 50 --unit uM
- --conc 0.05 --unit mM
Unit Abbreviations: mM=10-3, uM=10-6, nM=10-9, pM=10-12 |
| -l, --length |
Length of oligomers for hybridization. |
| -dist |
Limit the maximum distance between nucleotides that can pair. |
| -st, --start |
Truncate the target sequence to end at this nucleotide position. This can improve runtime since only a subset of the original sequence will be folded and scanned. Use this in conjunction with the --end parameter. |
| -en, --end |
Truncate the target sequence (see the --start parameter). |
| -fi, --filter |
Whether to use the siRNA prefilter to prefill functional siRNA.
0=No Prefilter (default); 1=Use Prefilter |
| -fold |
Only fold a fragment with the specified size (plus the length of the oligomer), which is centered on the binding region.
When FOLDSIZE > 1, MODE (-m) and SUBOPTIMAL (-s) must both be 2. |
| -ss, --from |
Start position of scanning on folded region of target. |
| -se, --to |
End position of scanning on folded region of target. |
| -m, --mode |
Determines how target structure is used:
1=Break Local Structure.
2=Refold target RNA after oligo binding.
3=No target structure considered.
|
| -s, --suboptimal |
Determines suboptimal structure options:
0=Do not consider suboptimal structures.
1=Use AllSub to generate suboptimal structures.
2=Use Partition Function to generate all suboptimal structures.
3=Use a heuristic method for both oligo-free and oligo-bound RNA.
4=Use stochastic sampling to generate 1000 structures.
|
| -sh, --shape |
Specify a SHAPE data
file to be used to generate restraints. These
restraints specifically use SHAPE pseudoenergy restraints.
Default is no SHAPE data file specified. |
| -test |
Perform self-tests. The parameter should be a list of space-separated test numbers, e.g.: -test '1 2 5' |
- J. Lu and D.H. Mathews. (2008) "Efficient siRNA selection using hybridization thermodynamics"
Nucleic Acids Research, 36:640-647.
- Reynolds, A., Leake, D., Boese, Q., Scaringe, S., Marshall, W.S. and Khvorova, A. (2004) Rational siRNA design for RNA interference. Nat Biotechnol., 22, 326-330.
- Lu, Z.J., Turner, D.H. and Mathews, D.H. (2006) A set of nearest neighbor parameters for predicting the enthalpy change of RNA secondary structure formation. Nucleic Acids Res., 34, 4912-4924.
- D.H. Mathews, M.D. Disney, J.L. Childs, S.J. Schroeder, M. Zuker, and D.H. Turner. "Incorporating Chemical Modification Constraints into a Dynamic Programming Algorithm for Prediction of RNA Secondary Structure"
Proceedings of the National Academy of Sciences USA, 101, 7287-7292, (2004).
- D.H. Mathews, J. Sabina, M. Zuker, and D. H. Turner "Expanded Sequence Dependence of Thermodynamic Parameters: Improves Prediction of RNA Secondary Structure" Journal of Molecular Biology, 288, 911-940, (1999).
- N. Sugimoto, S. Nakano, A. Katoh, H. Nakamura, T. Ohmichi, M. Yoneyama, and M. Sasaki
"Thermodynamic Parameters to Predict the Stability of RNA/DNA Hybrid Duplexes"
Biochemistry, 34, 11211-11216, (1995).
|
