article information

Jundishapur Journal of Microbiology: Vol.9, issue 4; e60047
published online: April 23, 2016
article type: Research Article
received: November 3, 2015
revised: December 5, 2015
accepted: February 16, 2016
DOI: https://doi.org/10.5812/jjm.34304
Crossmark

Potential siRNA Molecules for Nucleoprotein and M2/L Overlapping Region of Respiratory Syncytial Virus: In Silico Design

authors:

avatar Somayeh Shatizadeh Malekshahi 1 , avatar Ehsan Arefian 2 , avatar Vahid Salimi 1 , avatar Talat Mokhtari Azad 1 , avatar Jila Yavarian 1 , *

Virology Department, School of Public Health, Tehran University of Medical Sciences, Tehran, IR Iran
Department of Microbiology, School of Biology, College of Science, University of Tehran, Tehran, IR Iran

how to cite: Shatizadeh Malekshahi S, Arefian E, Salimi V, Mokhtari Azad T, Yavarian J. Potential siRNA Molecules for Nucleoprotein and M2/L Overlapping Region of Respiratory Syncytial Virus: In Silico Design. Jundishapur J Microbiol. 2016;9(4):e60047. https://doi.org/10.5812/jjm.34304.

Abstract

Background:

Human respiratory syncytial virus (RSV) is a leading cause of severe lower respiratory tract disease in the pediatric population, elderly and in immunosuppressed individuals. Respiratory syncytial virus is also responsible for bronchiolitis, pneumonia, and chronic obstructive pulmonary infections in all age groups. With this high disease burden and the lack of an effective RSV treatment and vaccine, there is a clear need for discovery and development of novel, effective and safe drugs to prevent and treat RSV disease. The most innovative approach is the use of small interfering RNAs (siRNAs) which represent a revolutionary new concept in human therapeutics. The nucleoprotein gene of RSV which is known as the most conserved gene and the M2/L mRNA, which encompass sixty-eight overlapping nucleotides, were selected as suitable targets for siRNA design.

Objectives:

The present study is aimed to design potential siRNAs for silencing nucleoprotein and an overlapping region of M2-L coding mRNAs by computational analysis.

Materials and Methods:

Various computational methods (target alignment, similarity search, secondary structure prediction, and RNA interaction calculation) have been used for siRNA designing against different strains of RSV.

Results:

In this study, seven siRNA molecules were rationally designed against the nucleoprotein gene and validated using various computational methods for silencing different strains of RSV. Additionally, three effective siRNA molecules targeting the overlapping region of M2/L mRNA were designed.

Conclusions:

This approach provides insight and a validated strategy for chemical synthesis of an antiviral RNA molecule which meets many sequence features for efficient silencing and treatment at the genomic level.

1. Background

Respiratory syncytial virus (RSV) is an enveloped, single-stranded, negative-sense RNA virus with a genome of ∼15,200 bases and a member of the genus Pneumovirus in the Paramyxoviridae family. Based on antigenic and genetic variability, RSV is classified into two groups, A and B (1, 2). Respiratory Syncytial Virus is the major cause of lower respiratory tract infections (LRTIs) with considerable morbidity and mortality in young children, the elderly, immunocompromised individuals, and patients with chronic lung diseases (1, 3, 4). There is some evidence supporting the theory that early LRTI with RSV raises the chance of developing chronic asthma-like respiratory symptoms during childhood (5). With this high disease burden and the lack of an effective RSV treatment and vaccine, there is a clear need for discovery and development of novel, effective and safe drugs to prevent and treat RSV diseases (4, 6, 7).

The RSV gene order is 3’-NS1-NS2-N-P-M-SH-G-F-M2-L-5’ (8). There are several conserved motifs at the boundaries of each gene, including the gene start and gene end signals which direct transcription initiation and termination, respectively (9). Viral proteins are produced by RSV-specific, RNA-dependent RNA polymerase (RdRp) that is packaged into the nucleocapsid (10). The nucleocapsid is the necessary unit for RNA replication and production of progeny virus particles (5). Respiratory syncytial virus protein sequence diversity studies showed that the surface protein G is highly divergent between strains and the internal proteins, such as neucleoprotein (N), are highly conserved, with more than 90% homology (5, 11). One distinctive feature of RSV is the 68 nucleotide overlap in the M2 and L genes, so that the M2 gene-end (GE) signal lies downstream of the L gene-start (GS) signal (9). The invention of RNA interference (RNAi), an endogenous and ubiquitous pathway, leads to revolutionary new concepts in human therapeutics (12). Since the discovery of RNAi, there has been rapid progress toward its use as a therapeutic approach against human diseases (13, 14).

There are two separate pathways of RNAi mechanisms, the siRNA (small interfering RNA) and the microRNA which can be used to effectively suppress the expression of genes in a sequence specific manner (15). The siRNA pathway, which exploits only cytoplasmic processes, is now the primary pathway utilized in RNAi-based drug development (12). siRNA-based drugs have distinct advantages over conventional small molecule or protein-based drugs, including high specificity, higher potency, and reduced toxicity (16, 17). Several key criteria are considered for selection of potent and functional siRNA molecules, such as the GC content of the sequence, specific nucleotides in distinct positions, sequence conservation, thermodynamic parameters of siRNA duplexes, and circumventing possible off-target effects (18).

Several studies have shown that siRNA can function as miRNA and regulate unintended transcripts. Such off-target effects via seed complementarities in their 3’UTRs can produce false positives in siRNA screens. In this regard, accurate prediction of miRNA-like off-target effects is important to mitigate undesired results (19, 20). siRNAs are loaded into the RNA-induced silencing complex (RISC) to cleave target mRNAs that share sequence identity with the siRNA. Within argonaute2-directed RISC, the siRNA unwind and antisense (guide) strand directs RISC to target the mRNA sequence after the sense (passenger) strand is released and cleaved by cellular nucleases (21, 22).

2. Objectives

The present study aimed to design potential siRNAs for silencing nucleoprotein and an overlapping region of M2-L coding mRNAs of RSV by computational analysis.

3. Materials and Methods

3.1. Sequence Retrieval and Multiple Sequence Alignment

Two hundred and nineteen complete CDS of N genes of different human RSV strains were collected from the viral gene bank database available at National Centre for Biotechnological Information (http://www.ncbi.nlm.nih.gov/). The accession numbers of these complete CDS are available upon request. The accession number of the sequence used for retrieval of the sixty-eight nt overlapping region of the M2-L gene was AF254574.1. In comparison to the N gene, the M2-L overlapping region is highly conserved in different RSV strains, then just one RSV reference strain was used for siRNA designing. Figure 1 depicts the structure of the M2/L overlap. Using MEGA5, all the complete CDS (219) were aligned for the selection of conserved regions. As one mismatch between the target mRNA and siRNA could dramatically affect siRNA efficiency, this method ensures the design of siRNAs that target all the consensus regions of transcripts. The regions in the target site with the highest homology by the mRNAs are considered in the following analysis.

A, Part of the RSV Genome Containing the M2/L Overlap; B, the Sequence of the M2/L Overlap is Shown as Negative-Sense RNA. GS: Gene Start; GE: Gene End.
A, Part of the RSV Genome Containing the M2/L Overlap; B, the Sequence of the M2/L Overlap is Shown as Negative-Sense RNA. GS: Gene Start; GE: Gene End.
Figure 1.

A, Part of the RSV Genome Containing the M2/L Overlap; B, the Sequence of the M2/L Overlap is Shown as Negative-Sense RNA. GS: Gene Start; GE: Gene End.

3.2. Potential siRNA Designing

In order to design siRNAs, two forms of software were used: the Whitehead siRNA selection program (http://sirna.wi.mit.edu/home.php) and Invitrogen BLOCK-iTTMRNAi designer (https://rnaidesigner.lifetechnologies.com/rnaiexpress/). The Whitehead siRNA selection website provides the flexibility of using several common siRNA sequence patterns (e.g. A: AAN19TT, B: NAN19NN, C: N2(CG)N8(AU)(N8(AU)N2) or a custom pattern (23). All of the siRNAs were aligned to find the most common siRNAs. The siRNA design rules can be categorized into both sequence and structural rules. Briefly, sequence rules consider that the target region should be preferably 50 - 100 nt downstream of the start codon because those regions are usually occupied by diverse protein factors (24). The G/C content of the binding site and siRNA should be in the range of 45% - 65% which provides the necessary stability for the siRNA (18). Numerous sequence rules regarding the preference or avoidance of the specific positions of the sense or antisense strand of the duplex have been suggested (24). Structural rules specify the structural accessibility of the target site and the thermodynamic features of the siRNA/target duplex. Apart from using the siRNA design tool, the resulting siRNAs were manually analyzed according to the parameters recommended by Ui-Tei et al. (25), Reynold et al. (26) and Jagla et al. (27) for optimal design.

3.3. Similarity Search

A BLAST search was used to compare a query sequence of our designed siRNAs with human genomic and transcripts database (http://blast.ncbi.nlm.nih.gov/) to identify any off-target sequence similarity in other non-targeted genomes by applying expected threshold value 10 and BLOSUM 62 matrix as parameters. siRNA with 15 or more consecutive bases of homology with any other mRNAs were excluded from consideration.

3.4. GC Calculation and Secondary Structure Prediction

For GC content calculation of predicted siRNA, the Oligo Analyzer Calculator (https://eu.idtdna.com/calc/analyzer) was used. The GC content of the siRNA duplex is a candidate which might correlate with siRNA functionality. Too high GC content may slow down duplex unwinding and too low GC content may decrease the efficiency of target mRNA recognition and hybridization (28). The Mfold web server (http://www.mfold.rna.albany.edu/) was used for secondary structure prediction.

3.5. Calculation of RNA-RNA Interaction Through Thermodynamics

For computing the free energy of siRNA folding and studying the thermodynamics of interaction between the predicted siRNA and the target gene, the RNAcofold program was used. RNAcofold calculates the hybridization energy and base pairing form of two RNA sequences (29).

4. Results and Discussion

RNAi methodology can be selected to target viral genes with high levels of sequence conservation. Nucleoprotein (∼1175 bp) is known as the most conserved gene of RSV that encapsidates the genome for protection from nuclease. The encapsidated genomic RNA serves as a template for transcription and replication (30). In transcription, subgenomic mRNAs were produced by a stop-start mechanism and, in replication, genome-length positive strand antigenomes were produced by read-through synthesis (30, 31). A key feature of both the genome and antigenome templates is that they remain coated with N protein all the time which, aside from preventing the formation of secondary structures in the RNA also protects them from nuclease attack (32). It can be considered that silencing of this gene will abrogate its essential functions.

In the present study, about 219 gene sequences from different strains of RSV were obtained from the GenBank, NCBI. The sequences were subjected to the MEGA5 computational tool in order to find the consensus regions of selected sequences. Nucleoprotein conserved regions were taken into account for optimal siRNA selection. Whitehead siRNA selection software and the Invitrogen BLOCK-iTTMRNAi Designer were used to provide functional siRNA design. Whitehead software also predicts off-target effects based on the seed complementarity of the target site among related species.

Fifty-six siRNA suggested by the Whitehead server and ten siRNA suggested by Invitrogen were manually analyzed for optimal design (mentioned in part 3.2 in methods). According to this system, a score of six or more significantly increases the probability for gene silencing. Finally, seven siRNAs were found to meet the rules and these are presented in Table 1. The siRNAs were sorted on the basis of targeting the largest number of nucleoprotein genes. All siRNAs were subjected to BLASTN to ensure that the siRNAs were specific for RSV. All the siRNAs retain GC content within 36% - 48%. Too high and too low GC content can impede the loading of siRNAs into RISC and weakening hybridization between the siRNA and target mRNA, respectively (33).

Table 1.

Characteristics of Effective siRNA Molecules Targeting Nucleoprotein

Target No.Location of Target Within mRNAsiRNA Target Within mRNAPredicted siRNA Duplex CandidateLength and PatternGC, %Tm, °CFree Energy of Folding With TargetΔG Binding, k/calmiRNA TargetNo. of Targeted Sequences
1290 - 311UAACAACACAUCGUCAAGACAU22 (B, C)36.441.40.4-34.8204 (78)150
UAACAACACAUCGUCAAGACAU
CCAGAATACAGGCATGACT
2439 - 457CCAGAATACAGGCATGACT19 (C)47.451.80.6-33.9142 (63)141
CCAGAAUACAGGCAUGACU
AGUCAUGCCUGUAUUCUGG
3378 - 396CATTGAGATAGAATCTAGA19 (B,C)31.642.40.6-2895 (45)126
CAUUGAGAUAGAAUCUAGA
CAUUGAGAUAGAAUCUAGA
4825 - 845TAGTGTGCAAGCAGAAATGGAT22 (C)40.946.30.4-37.7116 (53)60
UAGUGUGCAAGCAGAAAUGGAU
AUCCAUUUCUGCUUGCACACUA
5245 - 264TACTCAGAGATGCGGGATAT20 (B, C)4547.60.4-36272 (87)36
UACUCAGAGAUGCGGGAUAU
AUAUCCCGCAUCUCUGAGUA
6434 - 453TAGCTCCAGAATACAGGCAT20 (C)45530.4-36.2142 (63)15
UAGCUCCUGAAUACAGGCAU
AUGCCTGUAUUCUGGAGCUA
7610 - 628GAC ATAGCC AAC AGC TTCT19 (C)47.452.20.6-33.8198 (77)15
GACAUAGCCAACAGCUUCU
AGAAGCUGUUGGCUAUGUC

Since viruses are known to be genetically diverse (34), some researchers have focused more on conserved target sites to design siRNAs. Rosales et al. have used this approach to design siRNAs against NS4B and NS5 of the Dengue virus (35). Likewise Raza et al. predicted siRNA against the conserved region of HA and NA genes of the Influenza A virus (36). Naito et al. addressed the importance of conserved regions of HIV-1 for siRNA targeting (37).

Ideally, designed siRNAs must be complementary to their target sequences. However, several studies have shown that siRNA can bind to the mRNAs in a miRNA-like way through partial complementarity, leading to undesirable side-effects (38, 39). Here, for each siRNA candidate, the number of genes that the 7mer miRNA seed could potentially bind to is displayed (Table 1, miRNA target). For example, 204 (78): 204 means the miRNA seed (7mer) could potentially bind to 204 genes; 78 means that 78% of all 7mers would potentially bind to those 204 genes. In other words, only 22% of all 7mers would bind to more than 204 genes. This can be very useful for explaining the off-target effects in the siRNA experimental results.

siRNA sequences that contain internal repeats or palindromes may reduce the effective silencing by forming fold-back structures (40). Relative stability and the propensity to form internal hairpins can be predicted by melting temperature (Tm). Sequences with high Tm are prone to form hairpin structures. So, duplexes lacking stable internal repeats are better silencers (Tm < 60°C) (26).

RNA structure prediction of target mRNA is an important tool for designing siRNA. The accessibility of the target site can be important for estimating the effectiveness of the small RNA for regulating gene expression (41). A linear correlation has been reported between increasing stability of secondary structures and decreasing cleavage efficiency. Indeed, secondary structures which limit the accessibility of the target site by direct blockage reduce RISC-mediated cleavage efficiency (42). Using Mfold software, the secondary structure of nucleoprotein was predicted and the interaction sites with the proposed siRNAs are shown in Figure 2. It can be concluded that siRNA numbers 1, 3, 5 and 6 have more accessibility for target recognition.

siRNA Target Sites (Yellow) on the Nucleoprotein mRNA Secondary Structure Predicted Using Mfold
The arrows and digits indicate the initial site of interaction and number of siRNA, respectively.

The arrows and digits indicate the initial site of interaction and number of siRNA, respectively.

Figure 2.

siRNA Target Sites (Yellow) on the Nucleoprotein mRNA Secondary Structure Predicted Using Mfold

RNAcofold computes the hybridization energy and base pairing pattern of interacting RNA molecules by a combination of thermodynamic and kinetic considerations. It provides an extension of McCaskill’s partition function algorithm to compute base pairing probabilities (29). So, for the interaction of a target mRNA and its predicted siRNA, free energy of binding is calculated. The target-siRNA complex product prefers negative overall ΔG. A stable duplex has more negative ΔG than an unstable one (Tables 1 and 2).

Table 2.

Characteristics of Effective siRNA Molecules Targeting Overlapping Region of M2/L mRNA

Target No.Location of Target Within mRNAsiRNA Target Within mRNAPredicted siRNA Duplex CandidateLength and PatternGC, %Tm, °CFree Energy of Folding With TargetΔG binding (k/cal)miRNA Target
110 - 32GATCCCATTATTAATGGAA2032 (C)530.4-28.7263 (86)
GAUCCCAUUAUUAAUGGAA
UUCCAUUAAUAAUGGGAUC
29 - 31GGATCCCATTATTAATGGA2237 (C)46.30.6-30182 (74)
GGAUCCCAUUAUUAAUGGA
UCCAUUAAUAAUGGGAUCC
335 - 57CTAATGTTTATCTAACCGA2032 (C)47.60.5-26.61 (1)
CUAAUGUUUAUCUAACCGA
UCGGUUAGAUAAACAUUAG

It is worth noting that targeting the overlapping regions of the transcripts simultaneously reduced the level of both transcripts (43). The overlapping region of M2/L mRNA encompasses 68 nt between M2 GE signal and L GS signal. We hypothesized that targeting these two proteins at the same time would likely abolish their function. The properties of the three designed siRNAs targeting the overlapping region of the M2/L gene are summarized in Table 2. A computer-assisted analysis of the secondary structure of the M2/L RNA sequence was performed by Mfold (Figure 3). The results showed that siRNA 1 and 2 have better accessibility to target mRNA.

siRNA Target Sites (Yellow) on the M2/L mRNA Secondary Structure Predicted Using Mfold
The arrows and digits indicate the initial site of interaction and number of siRNA, respectively.

The arrows and digits indicate the initial site of interaction and number of siRNA, respectively.

Figure 3.

siRNA Target Sites (Yellow) on the M2/L mRNA Secondary Structure Predicted Using Mfold

In conclusion, there is a major unmet medical need for an effective therapy of RSV infection and respiratory viruses may be advantageous to siRNA therapeutics. This paper focused on antivirals that utilize the siRNA pathway as a forefront of drug development. The possibility of turning off the pathogenic genes appears to be an appealing approach to slow down or stop the disease progress of a wide variety of clinical pathogens.

Acknowledgements

References

  • 1.

    Martinelli M, Frati ER, Zappa A, Ebranati E, Bianchi S, Pariani E, et al. Phylogeny and population dynamics of respiratory syncytial virus (Rsv) A and B. Virus Res. 2014;189:293-302. [PubMed ID: 24954788]. https://doi.org/10.1016/j.virusres.2014.06.006.

  • 2.

    Tapia LI, Shaw CA, Aideyan LO, Jewell AM, Dawson BC, Haq TR, et al. Gene sequence variability of the three surface proteins of human respiratory syncytial virus (HRSV) in Texas. PLoS One. 2014;9(3). ee90786. [PubMed ID: 24625544]. https://doi.org/10.1371/journal.pone.0090786.

  • 3.

    Malekshahi SS, Azad TM, Yavarian J, Shahmahmoodi S, Naseri M, Rezaei F. Molecular detection of respiratory viruses in clinical specimens from children with acute respiratory disease in Iran. Pediatr Infect Dis J. 2010;29(10):931-3. [PubMed ID: 20879092].

  • 4.

    Laganas VA, Dunn EF, McLaughlin RE, Tiong-Yip CL, Yuzhakov O, Isabella VM, et al. Characterization of novel respiratory syncytial virus inhibitors identified by high throughput screen. Antiviral Res. 2015;115:71-4. [PubMed ID: 25542974]. https://doi.org/10.1016/j.antiviral.2014.12.012.

  • 5.

    O’Donnell R. Features of respiratory syncytial virus. Paediatr Child Health. 2009;19(1):43-7. https://doi.org/10.1016/j.paed.2008.09.010.

  • 6.

    Ramilo O. Evolution of prophylaxis: MoAb, siRNA, vaccine, and small molecules. Paediatr Respir Rev. 2009;10 Suppl 1:23-5. [PubMed ID: 19651397]. https://doi.org/10.1016/S1526-0542(09)70011-9.

  • 7.

    Salimi V, Hennus MP, Mokhtari-Azad T, Shokri F, Janssen R, Hodemaekers HM, et al. Opioid receptors control viral replication in the airways. Crit Care Med. 2013;41(1):205-14. [PubMed ID: 23222260]. https://doi.org/10.1097/CCM.0b013e31826767a8.

  • 8.

    Borchers AT, Chang C, Gershwin ME, Gershwin LJ. Respiratory syncytial virus--a comprehensive review. Clin Rev Allergy Immunol. 2013;45(3):331-79. [PubMed ID: 23575961]. https://doi.org/10.1007/s12016-013-8368-9.

  • 9.

    Hardy RW, Harmon SB, Wertz GW. Diverse gene junctions of respiratory syncytial virus modulate the efficiency of transcription termination and respond differently to M2-mediated antitermination. J Virol. 1999;73(1):170-6. [PubMed ID: 9847319].

  • 10.

    Empey KM, Peebles RJ, Kolls JK. Pharmacologic advances in the treatment and prevention of respiratory syncytial virus. Clin Infect Dis. 2010;50(9):1258-67. [PubMed ID: 20235830]. https://doi.org/10.1086/651603.

  • 11.

    Faghihloo E, Yavarian J, Jandaghi NZ, Shadab A, Azad TM. Genotype circulation pattern of human respiratory syncytial virus in Iran. Infect Genet Evol. 2014;22:130-3. [PubMed ID: 24462624]. https://doi.org/10.1016/j.meegid.2014.01.009.

  • 12.

    DeVincenzo JP. Harnessing RNA interference to develop neonatal therapies: from Nobel Prize winning discovery to proof of concept clinical trials. Early Hum Dev. 2009;85(10 Suppl):S31-5. [PubMed ID: 19833462]. https://doi.org/10.1016/j.earlhumdev.2009.08.013.

  • 13.

    Elbashir SM, Harborth J, Lendeckel W, Yalcin A, Weber K, Tuschl T. Duplexes of 21-nucleotide RNAs mediate RNA interference in cultured mammalian cells. Nature. 2001;411(6836):494-8. [PubMed ID: 11373684]. https://doi.org/10.1038/35078107.

  • 14.

    Kim DH, Rossi JJ. Strategies for silencing human disease using RNA interference. Nat Rev Genet. 2007;8(3):173-84. [PubMed ID: 17304245]. https://doi.org/10.1038/nrg2006.

  • 15.

    DeVincenzo JP. The promise, pitfalls and progress of RNA-interference-based antiviral therapy for respiratory viruses. Antivir Ther. 2012;17(1 Pt B):213-25. [PubMed ID: 22311654]. https://doi.org/10.3851/IMP2064.

  • 16.

    Bora RS, Gupta D, Mukkur TK, Saini KS. RNA interference therapeutics for cancer: challenges and opportunities (review). Mol Med Rep. 2012;6(1):9-15. [PubMed ID: 22576734]. https://doi.org/10.3892/mmr.2012.871.

  • 17.

    Xu C, Wang J. Delivery systems for siRNA drug development in cancer therapy. Asian J Pharm Sci. 2015;10(1):1-12. https://doi.org/10.1016/j.ajps.2014.08.011.

  • 18.

    Hajeri PB, Singh SK. siRNAs: their potential as therapeutic agents--Part I. Designing of siRNAs. Drug Discov Today. 2009;14(17-18):851-8. [PubMed ID: 19540928]. https://doi.org/10.1016/j.drudis.2009.06.001.

  • 19.

    Jackson AL, Bartz SR, Schelter J, Kobayashi SV, Burchard J, Mao M, et al. Expression profiling reveals off-target gene regulation by RNAi. Nat Biotechnol. 2003;21(6):635-7. [PubMed ID: 12754523]. https://doi.org/10.1038/nbt831.

  • 20.

    Scacheri PC, Rozenblatt-Rosen O, Caplen NJ, Wolfsberg TG, Umayam L, Lee JC, et al. Short interfering RNAs can induce unexpected and divergent changes in the levels of untargeted proteins in mammalian cells. Proc Natl Acad Sci U S A. 2004;101(7):1892-7. [PubMed ID: 14769924]. https://doi.org/10.1073/pnas.0308698100.

  • 21.

    Cejka D, Losert D, Wacheck V. Short interfering RNA (siRNA): tool or therapeutic? Clin Sci (Lond). 2006;110(1):47-58. [PubMed ID: 16336204]. https://doi.org/10.1042/CS20050162.

  • 22.

    Nolte A, Ott K, Rohayem J, Walker T, Schlensak C, Wendel HP. Modification of small interfering RNAs to prevent off-target effects by the sense strand. N Biotechnol. 2013;30(2):159-65. [PubMed ID: 23124048]. https://doi.org/10.1016/j.nbt.2012.10.001.

  • 23.

    Yuan B, Latek R, Hossbach M, Tuschl T, Lewitter F. siRNA Selection Server: an automated siRNA oligonucleotide prediction server. Nucleic Acids Res. 2004;32(Web Server issue):W130-4. [PubMed ID: 15215365]. https://doi.org/10.1093/nar/gkh366.

  • 24.

    Lagana A, Shasha D, Croce CM. Synthetic RNAs for Gene Regulation: Design Principles and Computational Tools. Front Bioeng Biotechnol. 2014;2:65. [PubMed ID: 25566532]. https://doi.org/10.3389/fbioe.2014.00065.

  • 25.

    Ui-Tei K, Naito Y, Takahashi F, Haraguchi T, Ohki-Hamazaki H, Juni A, et al. Guidelines for the selection of highly effective siRNA sequences for mammalian and chick RNA interference. Nucleic Acids Res. 2004;32(3):936-48. [PubMed ID: 14769950]. https://doi.org/10.1093/nar/gkh247.

  • 26.

    Reynolds A, Leake D, Boese Q, Scaringe S, Marshall WS, Khvorova A. Rational siRNA design for RNA interference. Nat Biotechnol. 2004;22(3):326-30. [PubMed ID: 14758366]. https://doi.org/10.1038/nbt936.

  • 27.

    Jagla B, Aulner N, Kelly PD, Song D, Volchuk A, Zatorski A, et al. Sequence characteristics of functional siRNAs. RNA. 2005;11(6):864-72. [PubMed ID: 15923373]. https://doi.org/10.1261/rna.7275905.

  • 28.

    Amarzguioui M, Prydz H. An algorithm for selection of functional siRNA sequences. Biochem Biophys Res Commun. 2004;316(4):1050-8. [PubMed ID: 15044091]. https://doi.org/10.1016/j.bbrc.2004.02.157.

  • 29.

    Bernhart SH, Tafer H, Muckstein U, Flamm C, Stadler PF, Hofacker IL. Partition function and base pairing probabilities of RNA heterodimers. Algorithms Mol Biol. 2006;1(1):3. [PubMed ID: 16722605]. https://doi.org/10.1186/1748-7188-1-3.

  • 30.

    Fearns R, Peeples ME, Collins PL. Increased expression of the N protein of respiratory syncytial virus stimulates minigenome replication but does not alter the balance between the synthesis of mRNA and antigenome. Virology. 1997;236(1):188-201. [PubMed ID: 9299631]. https://doi.org/10.1006/viro.1997.8734.

  • 31.

    Collins PL, Hill MG, Cristina J, Grosfeld H. Transcription elongation factor of respiratory syncytial virus, a nonsegmented negative-strand RNA virus. Proc Natl Acad Sci U S A. 1996;93(1):81-5. [PubMed ID: 8552680].

  • 32.

    Cowton VM, McGivern DR, Fearns R. Unravelling the complexities of respiratory syncytial virus RNA synthesis. J Gen Virol. 2006;87(Pt 7):1805-21. [PubMed ID: 16760383]. https://doi.org/10.1099/vir.0.81786-0.

  • 33.

    Kanasty RL, Whitehead KA, Vegas AJ, Anderson DG. Action and reaction: the biological response to siRNA and its delivery vehicles. Mol Ther. 2012;20(3):513-24. [PubMed ID: 22252451]. https://doi.org/10.1038/mt.2011.294.

  • 34.

    Schneider WL, Roossinck MJ. Genetic diversity in RNA virus quasispecies is controlled by host-virus interactions. J Virol. 2001;75(14):6566-71. [PubMed ID: 11413324]. https://doi.org/10.1128/JVI.75.14.6566-6571.2001.

  • 35.

    Villegas-Rosales PM, Mendez-Tenorio A, Ortega-Soto E, Barron BL. Bioinformatics prediction of siRNAs as potential antiviral agents against dengue viruses. Bioinformation. 2012;8(11):519-22. [PubMed ID: 22829722]. https://doi.org/10.6026/97320630008519.

  • 36.

    Raza A, Shareef H, Salim H, Khushal R, Bokhari H. Selection of predicted siRNA as potential antiviral therapeutic agent against influenza virus. Bioinformation. 2011;6(9):340-3. [PubMed ID: 21814391].

  • 37.

    Naito Y, Nohtomi K, Onogi T, Uenishi R, Ui-Tei K, Saigo K, et al. Optimal design and validation of antiviral siRNA for targeting HIV-1. Retrovirology. 2007;4:80. [PubMed ID: 17996047]. https://doi.org/10.1186/1742-4690-4-80.

  • 38.

    Birmingham A, Anderson EM, Reynolds A, Ilsley-Tyree D, Leake D, Fedorov Y, et al. 3' UTR seed matches, but not overall identity, are associated with RNAi off-targets. Nat Methods. 2006;3(3):199-204. [PubMed ID: 16489337]. https://doi.org/10.1038/nmeth854.

  • 39.

    Jackson AL, Burchard J, Schelter J, Chau BN, Cleary M, Lim L, et al. Widespread siRNA "off-target" transcript silencing mediated by seed region sequence complementarity. RNA. 2006;12(7):1179-87. [PubMed ID: 16682560]. https://doi.org/10.1261/rna.25706.

  • 40.

    Kirchner R, Vogtherr M, Limmer S, Sprinzl M. Secondary structure dimorphism and interconversion between hairpin and duplex form of oligoribonucleotides. Antisense Nucleic Acid Drug Dev. 1998;8(6):507-16. [PubMed ID: 9918115].

  • 41.

    Schroeder SJ. Advances in RNA structure prediction from sequence: new tools for generating hypotheses about viral RNA structure-function relationships. J Virol. 2009;83(13):6326-34. [PubMed ID: 19369331]. https://doi.org/10.1128/JVI.00251-09.

  • 42.

    Ameres SL, Martinez J, Schroeder R. Molecular basis for target RNA recognition and cleavage by human RISC. Cell. 2007;130(1):101-12. [PubMed ID: 17632058]. https://doi.org/10.1016/j.cell.2007.04.037.

  • 43.

    Hu X, Hipolito S, Lynn R, Abraham V, Ramos S, Wong-Staal F. Relative gene-silencing efficiencies of small interfering RNAs targeting sense and antisense transcripts from the same genetic locus. Nucleic Acids Res. 2004;32(15):4609-17. [PubMed ID: 15333693]. https://doi.org/10.1093/nar/gkh790.