Skip Navigation


Nucleic Acids Research Advance Access originally published online on May 15, 2008
Nucleic Acids Research 2008 36(11):3738-3745; doi:10.1093/nar/gkn266
This Article
Right arrow Abstract Freely available
Right arrow Print PDF (230K) Freely available
Right arrow Screen PDF (248K) Freely available
Right arrowOA All Versions of this Article:
36/11/3738    most recent
gkn266v1
Right arrow Alert me when this article is cited
Right arrow Alert me if a correction is posted
Services
Right arrow Email this article to a friend
Right arrow Similar articles in this journal
Right arrow Similar articles in PubMed
Right arrow Alert me to new issues of the journal
Right arrow Add to My Personal Archive
Right arrow Download to citation manager
Right arrow Commercial Re-use Guidelines
for Open Access NAR Content
Google Scholar
Right arrow Articles by Lu, Z. J.
Right arrow Articles by Mathews, D. H.
Right arrow Search for Related Content
PubMed
Right arrow PubMed Citation
Right arrow Articles by Lu, Z. J.
Right arrow Articles by Mathews, D. H.
Social Bookmarking
Add to CiteULike   Add to Connotea   Add to Del.icio.us  
What's this?

Nucleic Acids Research, 2008, Vol. 36, No. 11 3738-3745
© 2008 The Author(s)
This is an Open Access article distributed under the terms of the Creative Commons Attribution Non-Commercial License (http://creativecommons.org/licenses/by-nc/2.0/uk/) which permits unrestricted non-commercial use, distribution, and reproduction in any medium, provided the original work is properly cited.

Computational Biology

Fundamental differences in the equilibrium considerations for siRNA and antisense oligodeoxynucleotide design

Zhi John Lu1 and David H. Mathews1,2,*

1Department of Biochemistry and Biophysics and 2Department of Biostatistics & Computational Biology, University of Rochester Medical Center, Box 712, 601 Elmwood Avenue, Rochester, NY 14642, USA

*To whom correspondence should be addressed. Tel: +01 585 275 1734; Fax: +01 585 275 6007; Email: david_mathews{at}urmc.rochester.edu

Received February 22, 2008. Revised April 11, 2008. Accepted April 21, 2008.


   ABSTRACT
TOP
 ABSTRACT
INTRODUCTION
 MATERIALS AND METHODS
 RESULTS
 DISCUSSION
 REFERENCES
 
Both siRNA and antisense oligodeoxynucleotides (ODNs) inhibit the expression of a complementary gene. In this study, fundamental differences in the considerations for RNA interference and antisense ODNs are reported. In siRNA and antisense ODN databases, positive correlations are observed between the cost to open the mRNA target self-structure and the stability of the duplex to be formed, meaning the sites along the mRNA target with highest potential to form strong duplexes with antisense strands also have the greatest tendency to be involved in pre-existing structure. Efficient siRNA have less stable siRNA–target duplex stability than inefficient siRNA, but the opposite is true for antisense ODNs. It is, therefore, more difficult to avoid target self-structure in antisense ODN design. Self-structure stabilities of oligonucleotide and target correlate to the silencing efficacy of siRNA. Oligonucleotide self-structure correlations to efficacy of antisense ODNs, conversely, are insignificant. Furthermore, self-structure in the target appears to correlate with antisense ODN efficacy, but such that more effective antisense ODNs appear to target mRNA regions with greater self-structure. Therefore, different criteria are suggested for the design of efficient siRNA and antisense ODNs and the design of antisense ODNs is more challenging.


   INTRODUCTION
TOP
 ABSTRACT
 INTRODUCTION
MATERIALS AND METHODS
 RESULTS
 DISCUSSION
 REFERENCES
 
Antisense oligonucleotides, such as siRNA or antisense oligodeoxynucleotides (ODNs), can silence gene expression (1). siRNA associate with the protein–RNA complex called the RNA-induced silencing complex (RISC) to cleave the target mRNA or attenuate the gene expression with the RNAi pathway (2–4). Antisense ODNs also bind to a complementary region of the target mRNA and generally inhibit expression by stimulating degradation of the mRNA via RNase H (5–7).

The silencing efficacies of RNAi and antisense ODNs are found to correlate with their sequence features. Efficient siRNA have preference for low G/C content, A at position 3, U at position 10, absence of G at position 13, absence of G or C at position 19, etc. (8–14). Antisense ODN silencing efficacy also correlates highly with some specific motifs of oligonucleotide sequence, such as CCAC and ACTG (15,16). Additionally, the local secondary structure of the target mRNA also influences the binding affinity of siRNA (17–20) and antisense ODNs (21–24).

In this study, predicted free energy changes of hybridization of both antisense ODNs and siRNA are compared to inhibition efficacy databases to demonstrate contrasts in the hybridization terms that influence efficacy. Free energy changes of hybridization of the antisense oligonucleotide to the mRNA target are calculated using the OligoWalk algorithm (25,26), which uses the equilibrium shown in Figure 1. The equilibrium includes self-structure terms, Formula , Formula and Formula , which correspond to the free energy change of opening intramolecular pairs in the oligonucleotide, intermolecular pairs in the oligonucleotide and base pairs in the hybridization region of the target, respectively.


Figure 1
View larger version (12K):
[in this window]
[in a new window]
[Download PowerPoint slide]
  		Figure 1. Equilibrium considered in the OligoWalk algorithm (25,26) for siRNA and antisense ODNs. The equilibrium constants, Kduplex, Ktarget structure, Kintraoligonucleotide, and Kinteroligonucleotide are related to Figure 1 and Figure 1 by {Delta}G° = –RT ln K, respectively. Self-folding in the target and self-structure in the oligonucleotide both compete with the formation of the oligonucleotide–target complex. Only RNA secondary structure interactions are considered in the calculations. The longer arrow for each equilibrium shows the generally favored direction of the equilibrium, i.e. a negative folding free energy change is predicted for an equilibrium favoring the direction of the longer arrow.

 
The stability of duplex hybridization between antisense sequence and target is found, for the first time, to be significantly correlated with the stability of the target mRNA's self-structure at the hybridization region for both siRNA and antisense ODNs. Duplex stability is also shown to be correlated with the oligonucleotide self-structure stability for both siRNA and antisense ODNs. Different preferences of duplex stability, however, are observed for siRNA and antisense ODNs. Because RNAi is attenuated by the unwinding cost of opening the siRNA duplex, efficient siRNA (or miRNA) usually have less stable sense–antisense duplexes (27). This is just the opposite for efficient antisense ODNs, where tight hybridization to the target is apparently required. Furthermore, in addition to duplex stability, siRNA silencing efficacy also significantly correlates with other terms such as the self-structure stabilities of siRNA and target mRNA. These correlations are not as strong for antisense ODN efficacy.


   MATERIALS AND METHODS
TOP
 ABSTRACT
 INTRODUCTION
 MATERIALS AND METHODS
RESULTS
 DISCUSSION
 REFERENCES
 
Prediction of self-structure of oligonucleotide and target
To quantify the accessibility of oligonucleotide and target mRNA for hybridization, a free energy change of self-structure is predicted for opening base pairs in the region of complementarity to the target. A partition function (Q) calculation (28,29) is used to predict the ensemble free energy change (26). For example, the free energy cost of opening the self-structure of a target binding site is calculated using:

Formula
where Q is a partition function that sums the equilibrium constants for all possible structures, s.

Formula
R is gas constant, T is absolute temperature, which was set to 310.15 K in this study, Qunconstrained is the partition function of the native target structures, Qconstrained is the partition function of the target structures in the state where the oligonucleotide is able to bind. To predict the constrained partition function, the calculation is performed with a constraint that nucleotides in the binding region are forced single stranded. In order to reduce calculation time, Formula is calculated with a partition function of the local structure on mRNA binding site, i.e. only a certain number of nucleotides centered at the binding region (800 nt) is folded (26). It was previously demonstrated that local folding of 800 nt does not significantly affect the accuracy of the accessibility prediction (26). If the binding site is located close to the 5' or 3' end of the target, the same size of region is folded beginning from the end of the sequence, which means the binding site is not centered on the folding region.

For the oligonucleotide, all self-structure must be broken during duplex formation with the target, so the self-structure free energy change is predicted with:

Formula
Both unimolecular and bimolecular self structure are considered for the oligonucleotide using appropriate partition functions (26).

Thermodynamic parameters
Folding free energy changes for individual structures are predicted using nearest-neighbor models. For RNA structures, the nearest neighbor parameters from Turner and co-workers are used (30). For DNA structures, the nearest neighbor parameters for DNA from the RNAstructure program (30) are used. In the case of ODN hybridization to RNA targets, DNA–RNA duplex parameters are used for helix formation (31).

Databases
The experimental data for gene silencing efficacy of oligonucleotides is derived from two databases. One is derived from an antisense ODN database, AOBase (32). 418 ODNs targeting 28 mRNA are used for this study. Thirty ODNs were removed from the original database because these sequences are not consistent in sequence with the Genbank database (33). The silencing efficacy of each oligonucleotide is represented as ln(A), the natural logarithm of Activity, which is defined as the ratio of gene expression after antisense silencing over the untreated control. For the correlation calculations, any value of activity that is <0 is reset to 0.1% and any value that is >100% is reset to 99.9%. Two hundred and fifteen antisense ODNs induced more than 50% gene silencing (silencing efficacy = 1 – Activity), 103 induced more than 70% and 30 induced more than 90%. The second database is an siRNA database of experiments from Huesken et al. (34) at Novartis, which contains efficacy data for 2431 siRNAs targeting 31 mRNA sequences on random positions. Two thousand siRNAs have silencing efficacy >50%, 1222 of them have efficacy >70%, 369 have efficacy >90%. The silencing efficacies reported in the siRNA database are transformed to Activity (Activity = 1 – silencing efficacy) in order to calculate ln(A).

Statistical analysis
Linear correlation coefficients (r) are calculated between the free energy changes of duplex formation and free energy changes for self-structure formation in both oligonucleotide and target mRNA. Correlations are also explored between ln(Activity) and thermodynamic features involved in the equilibrium of binding for both siRNA and antisense ODNs. The significance of each linear correlation (Table 1) is tested with a two-tailed t-test. The t-test is performed with the Statistics-Basic-0.42 Perl module downloaded from: http://www.cpan.org and the data analysis tool in Excel 2004 (Microsoft Inc). For this study, a P-value of the test <0.05 is considered to be a significant correlation, i.e. rejection of the null hypothesis that the correlation is by chance.


View this table:
[in this window]
[in a new window]

  		Table 1. Correlations between ln(A)a and free energy change terms for both siRNA and Antisense ODNs

 

   RESULTS
TOP
 ABSTRACT
 INTRODUCTION
 MATERIALS AND METHODS
 RESULTS
DISCUSSION
 REFERENCES
 
The OligoWalk algorithm (25,26) was developed to predict the affinity of a structured oligonucleotide to a structured RNA target using the equilibrium shown in Figure 1. The prediction explicitly considers self-structure of the oligonucleotide and target, quantified by free energy changes calculated with the nearest neighbor model (30,31,35). The formation of self-structure (Formula and Formula ) competes with the hybridization of the antisense oligonucleotide to the target, which is driven by the favorable free energy change of duplex formation (Formula ). In addition, the stability difference between the duplex's two ends (Formula ) was also calculated because it is well known that efficient siRNA prefer a less stable duplex at the 5' end of the antisense strand (8).

Duplex stability requirements are different for siRNA and antisense ODNs
In RNAi, the siRNA duplex needs to unwind for loading the antisense strand on RISC and the antisense–target duplex needs to unwind for multiple turnover. Therefore, a general rule of siRNA design is a requirement for a low G/C content in the oligonucleotide (12). It was also reported that sense–antisense duplexes of efficient siRNA (or miRNA) are less stable than inefficient siRNA in previous studies (27,36). In this study, the same trend was observed in the Novartis siRNA database (34) (Figure 2A). The average Formula (–33.0 ± 4.6 kcal/mol) of efficient siRNA (silencing efficacy is not <70%) is 2.8 kcal/mol more than the average Formula (–35.8 ± 5.7 kcal/mol) of inefficient siRNA (silencing efficacy is <50%). Antisense ODNs, however, do not have to destabilize the duplex formation to be efficient and, in contrast to siRNA, require stable binding to the target (Figure 2B). The difference between the average Formula of efficient ODNs (–26.1 ± 4.2 kcal/mol) and inefficient ODNs (–24.9 ± 5.9 kcal/mol) is –1.2 kcal/mol.


Figure 2
View larger version (34K):
[in this window]
[in a new window]
[Download PowerPoint slide]
  		Figure 2. Oligonucleotide–target duplex stabilities in siRNA and antisense ODNs databases. The histograms of free energy changes of oligonucleotide–target duplexes (Figure 2) for efficient oligonucleotides (silencing efficacy is not <70%) and inefficient oligonucleotides (silencing efficacy is <50%) are shown in (A). the siRNA data set (34) and (B) the antisense ODNs data set (32). The duplex free energy change (Figure 2) is plotted against ln(A) for the siRNA database in (C) and the antisense ODNs database in (D). In (E), ln(A) is plotted as a function of the per base pair duplex free energy change for the ODNs database. ln(A) is the natural logarithm of Activity, which is the fraction of the targeted mRNA expression after antisense silencing compared to the control.

 
The ln(A), natural logarithm of message activity, is plotted versus duplex free energy changes (Formula ) of all binding sites for siRNA and antisense ODNs, in Figure 2C and D, respectively. The correlation coefficient of Formula and ln(A) is negative (r = –0.250) for siRNA, yet positive (r = 0.160) for antisense ODNs (Table 1). This shows again that less stable duplex formation is preferred by efficient siRNA but more stable duplexes are preferred by efficient antisense ODNs. Because the ODNs range in length from 9 to 21 nt, the correlation was also tested for ln(A) as a function of ODN duplex free energy change per base pair (Figure 2E). The correlation coefficient is 0.181, with a P-value of 0.000207. This is an even stronger correlation than that between ln(A) and Formula , which suggests that it is more important for antisense ODN activity to have stronger pairing per base pair than to simply favor longer helices.

Effect of self-structure appears different for siRNA and antisense ODNs
The silencing efficacy by siRNA has been previously demonstrated to be influenced by the secondary structures of both the antisense oligonucleotide and target mRNA (19,37). Each of the thermodynamic features calculated by OligoWalk, Formula Formula and Formula , were previously shown to correlate with the gene-silencing efficacy by siRNA (26) (Table 1).

In this study, the same terms were calculated for 418 antisense ODNs with reported inhibition activities (32,38). Significant correlations were also found between ln(A) and both Formula and Formula (Table 1). The correlation between ln(A) and Formula is 0.141, which means that the more efficient antisense ODNs apparently anneal to regions of mRNA with more stable self-structure to be disrupted. This correlation is exactly opposite that for siRNAs and is counter-intuitive.

Furthermore, in contrast to siRNA, no significant correlations were observed between the free energy changes of oligonucleotide self-structure and the silencing efficacy of antisense ODNs. This is probably simply because of the wide range of ODN lengths. The correlations between oligonucleotide self-structure and ln(A) can be improved using a Formula cutoff, where only sequences having Formula ≤ –30 kcal/mol are considered. After the cutoff, the lengths of the remaining antisense ODNs vary less, as most of them have 20 or 21 nt. For this subset of antisense ODNs, the antisense efficacy is statistically significantly influenced by the self-structure of oligonucleotide (Table 1). This is consistent with previous findings for antisense ODNs (39).

In a previous study of antisense ODNs (39), the self-structure of target was poorly predicted by either optimal or suboptimal structure prediction, which are not as rigorous as the partition function calculation used here. The ambiguous correlation between antisense efficacy and Formula in previous studies also comes from the relationship of hybridized duplex stability and self-structure accessibility of target (below).

Correlation between hybridized duplex stability and self-structure accessibility
To understand the basis of the different influence of target self-structure on siRNA and antisense ODNs, the relationship between Formula and the self-structure folding free energy changes was explored (Figure 3A and B). It was found that the duplex free energy change correlates significantly with each of the self-structure folding free energy changes (Formula and Formula ) for oligonucleotides in both the siRNA and antisense ODN database (Table 1). This correlation indicates that sequences that form stronger duplexes also tend to have stronger self-structures, both for the antisense sequences and for the target mRNA.


Figure 3
View larger version (29K):
[in this window]
[in a new window]
[Download PowerPoint slide]
  		Figure 3. Correlations between free energy change of hybridized duplex (Figure 3) and free energy cost of opening target base pairs for hybridization (Figure 3). The Figure 3 values were calculated with a partition function with a folding size of 800 nt centered on the binding site. (A) For the siRNA data set (34), the correlation coefficient is 0.5946 and the t-test P-value is 2.22 x 10–16. (B) For the antisense ODNs data set (32), the correlation coefficient is 0.5097 and the t-test P-value is 4.44 x 10–16. (C) For a full scan of an mRNA sequence (Genbank ID: X61940 [GenBank] , length: 1933 bases) from the 5' end to 3' end, the correlation coefficient is 0.6187 and the t-test P-value is 3.95 x 10–30.

 
To control for whether this correlation is a result of a selection bias in the design of antisense sequences in the databases, it was tested for all 19mer antisense sequences in a complete scan of an mRNA (Genbank ID: X61940 [GenBank] , length: 1933 bases) (Figure 3C). Again, the duplex-binding stability significantly correlates with the cost of opening the local self-structure of the target.

These correlations explain the apparent correlation that efficient antisense ODNs preferentially hybridize to targets with stronger self-structure. The strong correlation between target self-structure and duplex stability suggests the true preference for reduced target self-structure is obscured for ODNs because of the strong requirement for greater duplex stability. For siRNA, the correlation is readily observed because the requirement for reduced stability in the duplex also leads to a tendency for less target self-structure.

Differing equilibria for RNAi and antisense ODNs
In the initial step of RNA interference, the siRNA duplex needs to unwind (Figure 1), so the equilibrium constant in the direction of the necessary product is 1/Kduplex. Subsequently, the antisense strand hybridizes to mRNA, with the equilibrium constant for product of Kduplex. The cost of opening the mRNA self-tructure is 1/Ktarget structure. The overall equilibrium, including these three effects, relates to the log of activity:

Formula
 Because of the positive correlation between the hybridized duplex's stability (Formula , Kduplex) and the target structure's accessibility (Formula , Ktarget structure), the proportionality is then:

Formula
This suggests that siRNA design is simple because less stable duplexes target less stable target mRNA self-structures and efficient siRNA require both of these considerations at the same time.

For antisense ODNs, however, the opposite trend emerges because there is no duplex unwinding step involved in the inhibition mechanism. The cost of opening self-structure of target and oligonucleotide competes with the formation of hybridized duplex for antisense ODNs. When the self-structure thermodynamics are compared with ln(A) for antisense ODNs, the self-structure stability correlates with the ln(A) (Table 1), but in an unintuitive manner. The hybridized duplex stability apparently accounts the most for the efficacy of antisense ODNs. Therefore, in contrast to siRNA design, the requirement of stable duplex hybridization and unstable self-structure of target simultaneously makes design difficult.


   DISCUSSION
TOP
 ABSTRACT
 INTRODUCTION
 MATERIALS AND METHODS
 RESULTS
 DISCUSSION
REFERENCES
 
This study explores the underlying differences between the binding thermodynamics of RNAi and antisense ODNs. The preference of functional siRNA for low G/C content has been noted previously (27,40) and this leads to a lower stability for Formula (Figure 1). It is possible that the free energy cost to unwind the siRNA is more important than the stability of oligonucleotide–target duplex. This does not apply to antisense ODNs because antisense ODNs are delivered as a single-stranded agent. Another explanation is that turnover of RISC may be facilitated by having a lower duplex affinity between the siRNA and target. The cleavage mechanism of RISC has been well studied (41–44). RISC is an endonuclease that makes a single cleavage with preference to the middle of the mRNA binding site (10 nt from the 5' end of the siRNA) (41,45). The cleaved mRNA are released from RISC (41) and, presumably, the cleavage products are degraded in a common RNA degradation pathway because they do not have either the poly(A) tail or the 5' cap (45,46). The antisense siRNA in RISC is then intact for another round of cleavage (46). We speculate that it is possible that RISC needs to open the base pairs between the siRNA and target mRNA strand in order to release the siRNA and RISC before degradation of the mRNA. This would lead to a preference for reduced binding strength by siRNA.

In contrast, functional antisense ODNs are known to prefer a stronger duplex affinity. In the antisense mechanism, RNase H binds to an RNA–DNA duplex and degrades the RNA. Although RNase H belongs to a nucleotidyl-transferase super family of enzymes that includes RISC (47), it may have a different process of cleavage. Experimental evidence suggests that RNase H degrades the RNA of a hybrid DNA–RNA duplex in a processive manner (48). The entire portion of the RNA strand in complex with the antisense ODN is probably degraded by RNase H and release of antisense ODN is facilitated regardless of the strength of antisense–target duplexes. Therefore, a propensity for strong duplex formation is important because it would favor target binding.

A number of studies have addressed the rational design of siRNA (12,34,49) and antisense ODNs (16,39,50), but these studies did not consider the structure features involved in the antisense binding using our rigorous partition function method. It has been demonstrated that including self-structure terms of siRNA and target mRNA helps the selection of efficient siRNA (26). The correlations found in this study show that different thermodynamic features could also be considered to improve the design of antisense ODNs. Contributions from multiple features of antisense ODNs need to be considered in order to find an optimized combination for an efficient candidate.

Another important factor in design of effective oligonucleotides is the accessibility of the target self-structure, which competes with the hybridization of the oligonucleotide to the target. The paradox demonstrated here is that the sequence features conducive to a stronger formation duplex also contribute to less binding accessibility because of self-structure of the target. This is observed as the positive correlation between the free energy changes of duplex formation and self-structure. Because siRNA favors less stable duplexes, it is easy to simultaneously avoid target structure in siRNA design. For antisense ODNs design, however, it is difficult to design strong duplexes that will bind to regions with little self-structure. This means it is fundamentally more difficult to design antisense ODNs than siRNA.


   ACKNOWLEDGEMENTS
 
This work was supported by National Institutes of Health grant R01GM076485 to D.H.M. D.H.M. is an Alfred P. Sloan Foundation Research Fellow. The authors thank Professor Douglas H. Turner from the Department of Chemistry, University of Rochester, for discussions. The authors also thank two anonymous reviewers for constructive comments. Funding to pay the Open Access publication charges for this article was provided by the National Institutes of Health.

Conflict of interest statement. None declared.


   REFERENCES
TOP
 ABSTRACT
 INTRODUCTION
 MATERIALS AND METHODS
 RESULTS
 DISCUSSION
 REFERENCES
 

  1. Scherer LJ, Rossi JJ. Approaches for the sequence-specific knockdown of mRNA. Nat. Biotechnol. (2003) 21:1457–1465.[CrossRef][Web of Science][Medline]

  2. Mello CC, Conte D. Jr. Revealing the world of RNA interference. Nature (2004) 431:338–342.[CrossRef][Medline]

  3. Filipowicz W. RNAi: the nuts and bolts of the RISC machine. Cell (2005) 122:17–20.[CrossRef][Web of Science][Medline]

  4. Ameres SL, Martinez J, Schroeder R. Molecular basis for target RNA recognition and cleavage by human RISC. Cell (2007) 130:101–112.[Medline]

  5. Dias N, Stein CA. Antisense oligonucleotides: basic concepts and mechanisms. Mol. Cancer Ther. (2002) 1:347–355.[Free Full Text]

  6. Crooke ST. Molecular mechanisms of action of antisense drugs. Biochim. Biophys. Acta (1999) 1489:31–44.[Medline]

  7. Taylor MF, Wiederholt K, Sverdrup F. Antisense oligonucleotides: a systematic high-throughput approach to target validation and gene function determination. Drug Discov. Today (1999) 4:562–567.[CrossRef][Web of Science][Medline]

  8. Khvorova A, Reynolds A, Jayasena SD. Functional siRNAs and miRNAs exhibit strand bias. Cell (2003) 115:209–216.[CrossRef][Web of Science][Medline]

  9. Schwarz DS, Hutvagner G, Du T, Xu Z, Aronin N, Zamore PD. Asymmetry in the assembly of the RNAi enzyme complex. Cell (2003) 115:199–208.[CrossRef][Web of Science][Medline]

  10. Amarzguioui M, Prydz H. An algorithm for selection of functional siRNA sequences. Biochem. Biophys. Res. Commun. (2004) 316:1050–1058.[CrossRef][Web of Science][Medline]

  11. Harborth J, Elbashir SM, Vandenburgh K, Manninga H, Scaringe SA, Weber K, Tuschl T. Sequence, chemical, and structural variation of small interfering RNAs and short hairpin RNAs and the effect on mammalian gene silencing. Antisense Nucleic Acid Drug Dev. (2003) 13:83–105.[CrossRef][Web of Science][Medline]

  12. Reynolds A, Leake D, Boese Q, Scaringe S, Marshall WS, Khvorova A. Rational siRNA design for RNA interference. Nat. Biotechnol. (2004) 22:326–330.[CrossRef][Web of Science][Medline]

  13. Ui-Tei K, Naito Y, Takahashi F, Haraguchi T, Ohki-Hamazaki H, Juni A, Ueda R, Saigo K. Guidelines for the selection of highly effective siRNA sequences for mammalian and chick RNA interference. Nucleic Acids Res. (2004) 32:936–948.[Abstract/Free Full Text]

  14. Yuan B, Latek R, Hossbach M, Tuschl T, Lewitter F. siRNA Selection Server: an automated siRNA oligonucleotide prediction server. Nucleic Acids Res. (2004) 32:W130–W134.[Abstract/Free Full Text]

  15. Matveeva OV, Tsodikov AD, Giddings M, Freier SM, Wyatt JR, Spiridonov AN, Shabalina SA, Gesteland RF, Atkins JF. Identification of sequence motifs in oligonucleotides whose presence is correlated with antisense activity. Nucleic Acids Res. (2000) 28:2862–2865.[Abstract/Free Full Text]

  16. Camps-Valls G, Chalk AM, Serrano-Lopez AJ, Martin-Guerrero JD, Sonnhammer EL. Profiled support vector machines for antisense oligonucleotide efficacy prediction. BMC Bioinformatics (2004) 5:135.[CrossRef][Medline]

  17. Bohula EA, Salisbury AJ, Sohail M, Playford MP, Riedemann J, Southern EM, Macaulay VM. The efficacy of small interfering RNAs targeted to the type 1 insulin-like growth factor receptor (IGF1R) is influenced by secondary structure in the IGF1R transcript. J. Biol. Chem. (2003) 278:15991–15997.[Abstract/Free Full Text]

  18. Far RK, Sczakiel G. The activity of siRNA in mammalian cells is related to structural target accessibility: a comparison with antisense oligonucleotides. Nucleic Acids Res. (2003) 31:4417–4424.[Abstract/Free Full Text]

  19. Schubert S, Grunweller A, Erdmann VA, Kurreck J. Local RNA target structure influences siRNA efficacy: systematic analysis of intentionally designed binding regions. J. Mol. Biol. (2005) 348:883–893.[CrossRef][Web of Science][Medline]

  20. Westerhout EM, Berkhout B. A systematic analysis of the effect of target RNA structure on RNA interference. Nucleic Acids Res. (2007) 35:4322–4330.[Abstract/Free Full Text]

  21. Vickers TA, Wyatt JR, Freier SM. Effects of RNA secondary structure on cellular antisense activity. Nucleic Acids Res. (2000) 28:1340–1347.[Abstract/Free Full Text]

  22. Sohail M, Southern EM. Selecting optimal antisense reagents. Adv. Drug Deliv. Rev. (2000) 44:23–34.[CrossRef][Web of Science][Medline]

  23. Wagner RW, Matteucci MD, Grant D, Huang T, Froehler BC. Potent and selective inhibition of gene expression by an antisense heptanucleotide. Nat. Biotechnol. (1996) 14:840–844.[CrossRef][Web of Science][Medline]

  24. Milner N, Mir KU, Southern EM. Selecting effective antisense reagents on combinatorial oligonucleotide arrays. Nat. Biotech. (1997) 15:537–541.[CrossRef][Web of Science][Medline]

  25. Mathews DH, Burkard ME, Freier SM, Wyatt JR, Turner DH. Predicting oligonucleotide affinity to nucleic acid targets. RNA (1999) 5:1458–1469.[Abstract]

  26. Lu ZJ, Mathews DH. Efficient siRNA selection using hybridization thermodynamics. Nucleic Acids Res. (2008) 36:640–647.[Abstract/Free Full Text]

  27. Shabalina SA, Spiridonov AN, Ogurtsov AY. Computational models with thermodynamic and composition features improve siRNA design. BMC Bioinformatics (2006) 7:65.[CrossRef][Medline]

  28. McCaskill JS. The equilibrium partition function and base pair probabilities for RNA secondary structure. Biopolymers (1990) 29:1105–1119.[CrossRef][Web of Science][Medline]

  29. Mathews DH. Using an RNA secondary structure partition function to determine confidence in base pairs predicted by free energy minimization. RNA (2004) 10:1178–1190.[Abstract/Free Full Text]

  30. Mathews DH, Disney MD, Childs JL, Schroeder SJ, Zuker M, Turner DH. Incorporating chemical modification constraints into a dynamic programming algorithm for prediction of RNA secondary structure. Proc. Natl Acad. Sci. USA (2004) 101:7287–7292.[Abstract/Free Full Text]

  31. Sugimoto N, Nakano S, Katoh M, Matsumura A, Nakamuta H, Ohmichi T, Yoneyama M, Sasaki M. Thermodynamic parameters to predict stability of RNA/DNA hybrid duplexes. Biochemistry (1995) 34:11211–11216.[CrossRef][Web of Science][Medline]

  32. Bo X, Lou S, Sun D, Shu W, Yang J, Wang S. Selection of antisense oligonucleotides based on multiple predicted target mRNA structures. BMC Bioinformatics (2006) 7:122.[CrossRef][Medline]

  33. Benson DA, Karsch-Mizrachi I, Lipman DJ, Ostell J, Wheeler DL. GenBank. Nucleic Acids Res. (2007) 35:D21–D25.[Abstract/Free Full Text]

  34. Huesken D, Lange J, Mickanin C, Weiler J, Asselbergs F, Warner J, Meloon B, Engel S, Rosenberg A, Cohen D, et al. Design of a genome-wide siRNA library using an artificial neural network. Nat. Biotechnol. (2005) 23:995–1001.[CrossRef][Web of Science][Medline]

  35. Xia T, SantaLucia J. Jr, Burkard ME, Kierzek R, Schroeder SJ, Jiao X, Cox C, Turner DH. Thermodynamic parameters for an expanded nearest-neighbor model for formation of RNA duplexes with Watson–Crick pairs. Biochemistry (1998) 37:14719–14735.[CrossRef][Web of Science][Medline]

  36. Ichihara M, Murakumo Y, Masuda A, Matsuura T, Asai N, Jijiwa M, Ishida M, Shinmi J, Yatsuya H, Qiao S, et al. Thermodynamic instability of siRNA duplex is a prerequisite for dependable prediction of siRNA activities. Nucleic Acids Res. (2007) 35:e123.[Abstract/Free Full Text]

  37. Heale BS, Soifer HS, Bowers C, Rossi JJ. siRNA target site secondary structure predictions using local stable substructures. Nucleic Acids Res. (2005) 33:e30.[Abstract/Free Full Text]

  38. Bo X, Lou S, Sun D, Yang J, Wang S. AOBase: a database for antisense oligonucleotides selection and design. Nucleic Acids Res. (2006) 34:D664–D667.[Abstract/Free Full Text]

  39. Matveeva OV, Mathews DH, Tsodikov AD, Shabalina SA, Gesteland RF, Atkins JF, Freier SM. Thermodynamic criteria for high hit rate antisense oligonucleotide design. Nucleic Acids Res. (2003) 31:4989–4994.[Abstract/Free Full Text]

  40. Holen T, Amarzguioui M, Wiiger MT, Babaie E, Prydz H. Positional effects of short interfering RNAs targeting the human coagulation trigger tissue factor. Nucleic Acids Res. (2002) 30:1757–1766.[Abstract/Free Full Text]

  41. Martinez J, Tuschl T. RISC is a 5' phosphomonoester-producing RNA endonuclease. Genes Dev. (2004) 18:975–980.[Abstract/Free Full Text]

  42. Rivas FV, Tolia NH, Song JJ, Aragon JP, Liu J, Hannon GJ, Joshua-Tor L. Purified Argonaute2 and an siRNA form recombinant human RISC. Nat. Struct. Mol. Biol. (2005) 12:340–349.[CrossRef][Web of Science][Medline]

  43. Song JJ, Smith SK, Hannon GJ, Joshua-Tor L. Crystal structure of Argonaute and its implications for RISC slicer activity. Science (2004) 305:1434–1437.[Abstract/Free Full Text]

  44. Tolia NH, Joshua-Tor L. Slicer and the argonautes. Nat. Chem. Biol. (2007) 3:36–43.[CrossRef][Web of Science][Medline]

  45. Elbashir SM, Martinez J, Patkaniowska A, Lendeckel W, Tuschl T. Functional anatomy of siRNAs for mediating efficient RNAi in Drosophila melanogaster embryo lysate. EMBO J. (2001) 20:6877–6888.[CrossRef][Web of Science][Medline]

  46. Haley B, Zamore PD. Kinetic analysis of the RNAi enzyme complex. Nat. Struct. Mol. Biol. (2004) 11:599–606.[CrossRef][Web of Science][Medline]

  47. Nowotny M, Gaidamakov SA, Crouch RJ, Yang W. Crystal structures of RNase H bound to an RNA/DNA hybrid: substrate specificity and metal-dependent catalysis. Cell (2005) 121:1005–1016.[CrossRef][Web of Science][Medline]

  48. Gaidamakov SA, Gorshkova II, Schuck P, Steinbach PJ, Yamada H, Crouch RJ, Cerritelli SM. Eukaryotic RNases H1 act processively by interactions through the duplex RNA-binding domain. Nucleic Acids Res. (2005) 33:2166–2175.[Abstract/Free Full Text]

  49. Ladunga I. More complete gene silencing by fewer siRNAs: transparent optimized design and biophysical signature. Nucleic Acids Res. (2007) 35:433–440.[Abstract/Free Full Text]

  50. Far RK, Leppert J, Frank K, Sczakiel G. Technical improvements in the computational target search for antisense oligonucleotides. Oligonucleotides (2005) 15:223–233.[CrossRef][Web of Science][Medline]


Add to CiteULike CiteULike   Add to Connotea Connotea   Add to Del.icio.us Del.icio.us    What's this?



This Article
Right arrow Abstract Freely available
Right arrow Print PDF (230K) Freely available
Right arrow Screen PDF (248K) Freely available
Right arrowOA All Versions of this Article:
36/11/3738    most recent
gkn266v1
Right arrow Alert me when this article is cited
Right arrow Alert me if a correction is posted
Services
Right arrow Email this article to a friend
Right arrow Similar articles in this journal
Right arrow Similar articles in PubMed
Right arrow Alert me to new issues of the journal
Right arrow Add to My Personal Archive
Right arrow Download to citation manager
Right arrow Commercial Re-use Guidelines
for Open Access NAR Content
Google Scholar
Right arrow Articles by Lu, Z. J.
Right arrow Articles by Mathews, D. H.
Right arrow Search for Related Content
PubMed
Right arrow PubMed Citation
Right arrow Articles by Lu, Z. J.
Right arrow Articles by Mathews, D. H.
Social Bookmarking
Add to CiteULike   Add to Connotea   Add to Del.icio.us  
What's this?