Analysis by siRNA_profile program displays novel thermodynamic characteristics of highly functional siRNA molecules
© Muhonen et al; licensee BioMed Central Ltd. 2008
Received: 13 August 2007
Accepted: 21 May 2008
Published: 21 May 2008
Here we report the improved results of a new siRNA design program and analysis tool called siRNA_profile that reveals an additional criterion for bioinformatic search of highly functional siRNA sequences.
We retrospectively analysed over 2400 siRNA sequences from 34 genes and with known efficacies to categorize factors that differentiate highly, moderately and non-functional siRNA sequences in more detail. We tested the biological relevance of siRNA_profile in CHO cells stably expressing human TRACP.
The highly functional siRNA molecules exhibited lower overall stabilities than non-functional siRNAs after taking into consideration all the nucleotides from 5'-terminus to the 3'-terminus along the siRNA molecule, in addition to the 5'-section of the antisense strand and the region between 9–14 nucleotides as previously has been acknowledged. Comparison of the siRNA_profile program to five other programs resulted in a wide range of selected siRNA sequences with diverse gene silencing capacities, even when the target was only 197 nucleotides long. Six siRNA design programs selected 24 different siRNA sequences, and only 6 of them were selected by two or more programs. The other 18 sequences were individually selected by these six programs.
Low general stability of dsRNA plays a significant role in the RNAi pathway and is a recommended criterion to consider, in addition to 5'-instability, internal instability, nucleotide preferences and target mRNA position, when designing highly efficient siRNAs.
RNA interference (RNAi) is a gene silencing mechanism where short interfering RNA (siRNAs) and microRNA (miRNAs) molecules inhibit the transcription and translation of target genes in a sequence-specific manner [1–3]. siRNAs are exogenously produced ~21 nucleotides long, double stranded RNA molecules with complete complementarity to the target sequence. miRNAs are a family of endogenously encoded small non-coding RNAs, derived by processing of short RNA hairpins, that can inhibit the translation of mRNAs bearing partial complementarity to the target sequences. RNAi has been acknowledged as a practical tool for new drug target discovery and RNAi drug development in mammalian cells . Therefore designing highly functional siRNA molecules has become an essential part of RNAi methodology. We have developed a novel and user-friendly siRNA design algorithm siRNA_profile with multiple options for minimizing the identification of non-functional, unspecific and immunostimulatory siRNA molecules. The analysis of functional and non-functional siRNA molecules were done by the siRNA_profile program to demonstrate in more detail the characteristics of highly functional siRNA molecules to help scientists in their search for theoretically and biologically efficient siRNAs.
Conventional siRNA design criteria
The rapid development of RNAi applications has revealed the need for efficient and specific siRNA design and analysis tools to maximize the efficiency while minimizing possible side-effects . Computational methods and neural networks are tools approaching ideal siRNA design; however, so far none of them are perfect. Currently, thermodynamic characteristics of functional siRNA molecules guide siRNA design strategy. First, it has been shown that thermodynamic differences in the base-pairing stabilities of the 5'-ends of both siRNA and miRNA molecules play a critical role in determining which strand initiates RNA induced silencing complex (RISC) activation [6, 7]. To achieve an efficient RNAi effect, activated RISC should be able to silence multiple copies of the target mRNA. A second criterion for effective siRNAs is low internal stability in the cleavage region between 9 to 14 nucleotides (calculated from the 5'-terminus) of the antisense strand. This is believed to have a critical role in mRNA cleavage and it may also help to release RISC from the cleaved target . Third, there are findings of the nucleotide preferences over the length of siRNA sequence [8, 9].
siRNA_profile program design
The siRNA_profile program and a full text including additional data and a printable help page are available in the siRNA_profile program web page . Briefly, the siRNA_profile program is based on findings of the asymmetric differences between functional and non-functional siRNAs [6, 7, 11] and on our studies of positional nucleotide differences and average dsRNA stability along the siRNA antisense strand. We have also incorporated some recommendations of Elbashir, S.M., et al.  and developed a novel, interactive and user-friendly siRNA design algorithm with multiple options for minimizing unspecific siRNA design. Previously, it has been recognized that 5'-UGUGU-3' motifs have immunostimulatory potential in synthetic siRNA molecules, in addition to CpG motifs in single stranded RNAs and DNA oligomers [12–14]. We have incorporated a sensor able to recognize immunostimulatory motifs intending to avoid unnecessary false phenotypes by siRNA molecules and, additionally, by anti-miRNA oligomers. In addition, a scoring system was integrated into the siRNA_profile program. It was adjusted based on our findings on nucleotide, purine and pyrimidine distribution along functional and non-functional siRNA sequences [the scoring system described in siRNA_profile web page].
The siRNA_profile program uses free energy values for calculation of average internal stability profiles. The average internal stability profiles were calculated as the sum of stability in five nucleotide windows by using the nearest-neighbour method as previously described [7, 11]. Our program browses and calculates the target sequence from the antisense point of view, from the 5' - to the 3'- direction. However, it does not utilize nucleotides on the mRNA beyond the 3'- end of the siRNA as previously has been described , because this calculation method may easily lead to false bending of the profile. This unique and beneficial feature of the siRNA_profile program allows the calculation of the factual, not the false dsRNA thermodynamic characters.
The siRNA_profile program was originally developed on Red Hat Linux 9.0 and compiled to an executable file. The algorithm was written in C-programming language and the CGI application was written in Perl programming language. The application uses the user-given input values from the siRNA_profile website to analyze and display the results.
To demonstrate the usefulness of siRNA_profile program as a siRNA candidate search service, the biological efficacies of siRNA molecules targeting tartrate-resistant acid phosphatase (TRACP) experimentally designed by siRNA_ profile were validated in a Chinese Hamster Ovary (CHO) cell line stably over-expressing human TRACP under a CMV promoter  [data shown in siRNA_profile web page].
Comparison of the siRNA_profile program and five other siRNA design programs
Number of hits
In conclusion, our results suggest that low general stability of dsRNA plays a significant role in the RNAi pathway and is recommended to be used as an additional guideline for highly efficient siRNA design. Additionally, all known immunostimulatory motifs to date are highlighted in colour to prevent analysis of false phenotypes and the output of the siRNA_profile program displays the average internal energy profiles for each siRNA candidate. According to our results, the siRNA sequences exhibiting both favourable asymmetric properties and low overall stabilities according to Figure 1 are recommended for gene knock down experiments.
Availability and requirements
siRNA_profile program availability: http://bonebiology.utu.fi/pimaki/main.html
The authors thank Dr. Pentti Riikonen (Department of Information Technology, University of Turku), Anne Seppänen and Janne Harjunpää for providing valuable help during development of the siRNA_profile algorithm and the web page application. The authors also thank Dr. Tiina Laitala-Leinonen (Bone Biology Research Consortium, University of Turku) for her assistance in formatting the text for publication.
- Lagos-Quintana M, Rauhut R, Lendeckel W, Tuschl T: Identification of novel genes coding for small expressed RNAs. Science. 2001, 294 (5543): 853-858. 10.1126/science.1064921.View ArticlePubMedGoogle Scholar
- Elbashir SM, Lendeckel W, Tuschl T: RNA interference is mediated by 21- and 22-nucleotide RNAs. Genes Dev. 2001, 15 (2): 188-200. 10.1101/gad.862301.PubMed CentralView ArticlePubMedGoogle Scholar
- Fire A, Xu S, Montgomery MK, Kostas SA, Driver SE, Mello CC: Potent and specific genetic interference by double-stranded RNA in Caenorhabditis elegans. Nature. 1998, 391 (6669): 806-811. 10.1038/35888.View ArticlePubMedGoogle Scholar
- Soutschek J, Akinc A, Bramlage B, Charisse K, Constien R, Donoghue M, Elbashir S, Geick A, Hadwiger P, Harborth J, John M, Kesavan V, Lavine G, Pandey RK, Racie T, Rajeev KG, Rohl I, Toudjarska I, Wang G, Wuschko S, Bumcrot D, Koteliansky V, Limmer S, Manoharan M, Vornlocher HP: Therapeutic silencing of an endogenous gene by systemic administration of modified siRNAs. Nature. 2004, 432 (7014): 173-178. 10.1038/nature03121.View ArticlePubMedGoogle Scholar
- Pei Y, Tuschl T: On the art of identifying effective and specific siRNAs. Nat Methods. 2006, 3 (9): 670-676. 10.1038/nmeth911.View ArticlePubMedGoogle Scholar
- Schwarz DS, Hutvagner G, Du T, Xu Z, Aronin N, Zamore PD: Asymmetry in the assembly of the RNAi enzyme complex. Cell. 2003, 115 (2): 199-208. 10.1016/S0092-8674(03)00759-1.View ArticlePubMedGoogle Scholar
- Khvorova A, Reynolds A, Jayasena SD: Functional siRNAs and miRNAs exhibit strand bias. Cell. 2003, 115 (2): 209-216. 10.1016/S0092-8674(03)00801-8.View ArticlePubMedGoogle Scholar
- Reynolds A, Leake D, Boese Q, Scaringe S, Marshall WS, Khvorova A: Rational siRNA design for RNA interference. Nat Biotechnol. 2004, 22 (3): 326-330. 10.1038/nbt936.View ArticlePubMedGoogle Scholar
- Jagla B, Aulner N, Kelly PD, Song D, Volchuk A, Zatorski A, Shum D, Mayer T, De Angelis DA, Ouerfelli O, Rutishauser U, Rothman JE: Sequence characteristics of functional siRNAs. Rna. 2005, 11 (6): 864-872. 10.1261/rna.7275905.PubMed CentralView ArticlePubMedGoogle Scholar
- siRNA_profile program [http://bonebiology.utu.fi/pimaki/main.html].Google Scholar
- Freier SM, Kierzek R, Jaeger JA, Sugimoto N, Caruthers MH, Neilson T, Turner DH: Improved free-energy parameters for predictions of RNA duplex stability. Proc Natl Acad Sci U S A. 1986, 83 (24): 9373-9377. 10.1073/pnas.83.24.9373.PubMed CentralView ArticlePubMedGoogle Scholar
- Hemmi H, Takeuchi O, Kawai T, Kaisho T, Sato S, Sanjo H, Matsumoto M, Hoshino K, Wagner H, Takeda K, Akira S: A Toll-like receptor recognizes bacterial DNA. Nature. 2000, 408 (6813): 740-745. 10.1038/35047123.View ArticlePubMedGoogle Scholar
- Judge AD, Sood V, Shaw JR, Fang D, McClintock K, MacLachlan I: Sequence-dependent stimulation of the mammalian innate immune response by synthetic siRNA. Nat Biotechnol. 2005, 23 (4): 457-462. 10.1038/nbt1081.View ArticlePubMedGoogle Scholar
- Sugiyama T, Gursel M, Takeshita F, Coban C, Conover J, Kaisho T, Akira S, Klinman DM, Ishii KJ: CpG RNA: identification of novel single-stranded RNA that stimulates human CD14+CD11c+ monocytes. J Immunol. 2005, 174 (4): 2273-2279.View ArticlePubMedGoogle Scholar
- Huesken D, Lange J, Mickanin C, Weiler J, Asselbergs F, Warner J, Meloon B, Engel S, Rosenberg A, Cohen D, Labow M, Reinhardt M, Natt F, Hall J: Design of a genome-wide siRNA library using an artificial neural network. Nat Biotechnol. 2005, 23 (8): 995-1001. 10.1038/nbt1118.View ArticlePubMedGoogle Scholar
- Janckila AJ, Parthasarathy RN, Parthasarathy LK, Seelan RS, Yam LT: Stable expression of human tartrate-resistant acid phosphatase isoforms by CHO cells. Clin Chim Acta. 2002, 326 (1-2): 113-122. 10.1016/S0009-8981(02)00280-2.View ArticlePubMedGoogle Scholar
- Henschel A, Buchholz F, Habermann B: DEQOR: a web-based tool for the design and quality control of siRNAs. Nucleic Acids Res. 2004, 32 (Web Server issue): W113-20. 10.1093/nar/gkh408.PubMed CentralView ArticlePubMedGoogle Scholar
- Naito Y, Yamada T, Ui-Tei K, Morishita S, Saigo K: siDirect: highly effective, target-specific siRNA design software for mammalian RNA interference. Nucleic Acids Res. 2004, 32 (Web Server issue): W124-9. 10.1093/nar/gkh442.PubMed CentralView ArticlePubMedGoogle Scholar
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.