Strand-specific 5'-O-methylation of siRNA duplexes controls guide strand selection and targeting specificity

Small interfering RNAs (siRNAs) and microRNAs (miRNAs) guide catalytic sequence-specific cleavage of fully or nearly fully complementary target mRNAs or control translation and/or stability of many mRNAs that share 6-8 nucleotides (nt) of complementarity to the siRNA and miRNA 5' end. siRNA- and miRNA-containing ribonucleoprotein silencing complexes are assembled from double-stranded 21- to 23-nt RNase III processing intermediates that carry 5' phosphates and 2-nt overhangs with free 3' hydroxyl groups. Despite the structural symmetry of a duplex siRNA, the nucleotide sequence asymmetry can generate a bias for preferred loading of one of the two duplex-forming strands into the RNA-induced silencing complex (RISC). Here we show that the 5'-phosphorylation status of the siRNA strands also acts as an important determinant for strand selection. 5'-O-methylated siRNA duplexes refractory to 5' phosphorylation were examined for their biases in siRNA strand selection. Asymmetric, single methylation of siRNA duplexes reduced the occupancy of the silencing complex by the methylated strand with concomitant elimination of its off-targeting signature and enhanced off-targeting signature of the phosphorylated strand. Methylation of both siRNA strands reduced but did not completely abolish RNA silencing, without affecting strand selection relative to that of the unmodified siRNA. We conclude that asymmetric 5' modification of siRNA duplexes can be useful for controlling targeting specificity.     

Development of a software tool and criteria evaluation for efficient design of small interfering RNA

RNA interference can be used as a tool for gene silencing mediated by small interfering RNAs (siRNA). The critical step in effective and specific RNAi processing is the selection of suitable constructs. Major design criteria, i.e., Reynolds’s design rules, thermodynamic stability, internal repeats, immunostimulatory motifs were emphasized and implemented in the siRNA design tool. The tool provides thermodynamic stability score, GC content and a total score based on other design criteria in the output. The viability of the tool was established with different datasets. In general, the siRNA constructs produced by the tool had better thermodynamic score and positional properties. Comparable thermodynamic scores and better total scores were observed with the existing tools. Moreover, the results generated had comparable off-target silencing effect. Criteria evaluations with additional criteria were achieved in WEKA.     

Rational siRNA design for RNA interference

Short-interfering RNAs suppress gene expression through a highly regulated enzyme-mediated process called RNA interference (RNAi). RNAi involves multiple RNA-protein interactions characterized by four major steps: assembly of siRNA with the RNA-induced silencing complex (RISC), activation of the RISC, target recognition and target cleavage. These interactions may bias strand selection during siRNA-RISC assembly and activation, and contribute to the overall efficiency of RNAi. To identify siRNA-specific features likely to contribute to efficient processing at each step, we performed a systematic analysis of 180 siRNAs targeting the mRNA of two genes. Eight characteristics associated with siRNA functionality were identified: low G/C content, a bias towards low internal stability at the sense strand 3'-terminus, lack of inverted repeats, and sense strand base preferences (positions 3, 10, 13 and 19). Further analyses revealed that application of an algorithm incorporating all eight criteria significantly improves potent siRNA selection. This highlights the utility of rational design for selecting potent siRNAs and facilitating functional gene knockdown studies.     

A comparison of siRNA efficacy predictors

Short interfering RNA (siRNA) efficacy prediction algorithms aim to increase the probability of selecting target sites that are applicable for gene silencing by RNA interference. Many algorithms have been published recently, and they base their predictions on such different features as duplex stability, sequence characteristics, mRNA secondary structure, and target site uniqueness. We compare the performance of the algorithms on a collection of publicly available siRNAs. First, we show that our regularized genetic programming algorithm GPboost appears to have a higher and more stable performance than other algorithms on the collected datasets. Second, several algorithms gave close to random classification on unseen data, and only GPboost and three other algorithms have a reasonably high and stable performance on all parts of the dataset. Third, the results indicate that the siRNAs' sequence is sufficient input to siRNA efficacy algorithms, and that other features that have been suggested to be important may be indirectly captured by the sequence.     

Focusing on RISC assembly in mammalian cells

RISC (RNA-induced silencing complex) is a central protein complex in RNAi, into which a siRNA strand is assembled to become effective in gene silencing. By using an in vitro RNAi reaction based on Drosophila embryo extract, an asymmetric model was recently proposed for RISC assembly of siRNA strands, suggesting that the strand that is more loosely paired at its 5′ end is selectively assembled into RISC and results in target gene silencing. However, in the present study, we were unable to establish such a correlation in cell-based RNAi assays, as well as in large-scale RNAi data analyses. This suggests that the thermodynamic stability of siRNA is not a major determinant of gene silencing in mammalian cells. Further studies on fork siRNAs showed that mismatch at the 5′ end of the siRNA sense strand decreased RISC assembly of the antisense strand, but surprisingly did not increase RISC assembly of the sense strand. More interestingly, measurements of melting temperature showed that the terminal stability of fork siRNAs correlated with the positions of the mismatches, but not gene silencing efficacy. In summary, our data demonstrate that there is no definite correlation between siRNA stability and gene silencing in mammalian cells, which suggests that instead of thermodynamic stability, other features of the siRNA duplex contribute to RISC assembly in RNAi.     

More complete gene silencing by fewer siRNAs: transparent optimized design and biophysical signature

Highly accurate knockdown functional analyses based on RNA interference (RNAi) require the possible most complete hydrolysis of the targeted mRNA while avoiding the degradation of untargeted genes (off-target effects). This in turn requires significant improvements to target selection for two reasons. First, the average silencing activity of randomly selected siRNAs is as low as 62%. Second, applying more than five different siRNAs may lead to saturation of the RNA-induced silencing complex (RISC) and to the degradation of untargeted genes. Therefore, selecting a small number of highly active siRNAs is critical for maximizing knockdown and minimizing off-target effects. To satisfy these needs, a publicly available and transparent machine learning tool is presented that ranks all possible siRNAs for each targeted gene. Support vector machines (SVMs) with polynomial kernels and constrained optimization models select and utilize the most predictive effective combinations from 572 sequence, thermodynamic, accessibility and self-hairpin features over 2200 published siRNAs. This tool reaches an accuracy of 92.3% in cross-validation experiments. We fully present the underlying biophysical signature that involves free energy, accessibility and dinucleotide characteristics. We show that while complete silencing is possible at certain structured target sites, accessibility information improves the prediction of the 90% active siRNA target sites. Fast siRNA activity predictions can be performed on our web server at http://optirna.unl.edu/.     

Chemically Modified Oligonucleotide–Increased Stability Negatively Correlates with Its Efficacy Despite Efficient Electrotransfer

Despite great potential for disease treatment, small interfering RNA (siRNA) development has been hampered due to its poor stability and the lack of efficient delivery method. To overcome the sensitivity, new generations of chemically modified oligonucleotides have been developed such as the locked nucleic acid (LNA). LNA substitution in an siRNA sequence (siLNA) is supposed to increase its stability and its affinity for its complementary sequence. The purpose of this study was to evaluate the potential benefit of an anti-GFP siLNA using the biophysical delivery method electropermeabilization. We used two types of electrical conditions: electrochemotherapy (ECT), a condition for efficient transfer of small molecules in clinics, and electrogenotherapy (EGT), a condition for efficient transfer of macromolecules. We first confirmed that siLNA was indeed more stable in mouse serum than unmodified siRNA. After determining the ECT and EGT optimal electrical parameters for a human colorectal carcinoma cell line (HCT-116) expressing eGFP, we showed that modifications of siRNA do not interfere with electrotransfer efficiency. However, despite its higher stability and its high electrotransfer efficacy, siLNA was less efficient for eGFP silencing compared to the electrotransferred, unmodified siRNA regardless of the electrical conditions used. Our study highlighted the care that is needed when designing chemically modified oligonucleotides.     

Gene silencing activity of siRNA molecules containing phosphorodithioate substitutions

Chemically synthesized small interfering RNAs (siRNAs) have been widely used to identify gene function and hold great potential in providing a new class of therapeutics. Chemical modifications are desired for therapeutic applications to improve siRNA efficacy. Appropriately protected ribonucleoside-3'-yl S-[β-(benzoylmercapto)ethyl]pyrrolidino-thiophosphoramidite monomers were prepared for the synthesis of siRNA containing phosphorodithioate (PS2) substitutions in which the two non-bridging oxygen atoms are replaced by sulfur atoms. A series of siRNAs containing PS2 substitutions have been strategically designed, synthesized, and evaluated for their gene silencing activities. These PS2-siRNA duplexes exhibit an A-form helical structure similar to unmodified siRNA. The effect of PS2 substitutions on gene silencing activity is position-dependent, with certain PS2-siRNAs showing activity significantly higher than that of unmodified siRNA. The relative gene silencing activities of siRNAs containing either PS2 or phosphoromonothioate (PS) linkages at identical positions are variable and depend on the sites of modification. 5'-Phosphorylation of PS2-siRNAs has little or no effect on gene silencing activity. Incorporation of PS2 substitutions into siRNA duplexes increases their serum stability. These results offer preliminary evidence of the potential value of PS2-modified siRNAs.     

Silencing activity of 2'-O-methyl modified anti-MDR1 siRNAs with mismatches in the central part of the duplexes

The thermodynamic properties of siRNA duplexes are important for their silencing activity. siRNAs with high thermodynamic stability of both the central part of the duplex and in the whole, usually display low silencing activity. Destabilization of the central part of the siRNA duplex could increase its silencing activity. However, mismatches located in the central part of the duplex could substantially decrease the amount of RNAi efficacy, hindering active RISC formation and function. In this study, we examined the impact of duplex destabilization by nucleotide substitutions in the central part (7-10 nt counting from the 5'-end of the antisense strand) of the nuclease-resistant siRNA on its silencing activity.     

Viral resistance mediated by shRNA depends on the sequence similarity and mismatched sites between the target sequence and siRNA

Viral resistance can be effectively induced in transgenic plants through their silencing machinery. Thus, we designed nine short hairpin RNAs (shRNA) constructs to target nuclear inclusion protein b (NIb), helper component proteinase (HC-Pro), cylindrical inclusion protein (CI) and viral protein genome linked (VPg) genes of Potato virus Y (PVY^N) and Tobacco etch virus (TEV-SD1). The shRNAs were completely complementary to the genes of PVY^N, and contained 1–3 nt mismatches to the genes of TEV-SD1. To study the specificity of gene silencing in shRNA-mediated viral resistance, the constructs were introduced into tobacco plants. The results of viral resistance assay revealed that these nine kinds of transgenic tobacco plants can effectively induce viral resistance against both PVY^N and TEV-SD1, and the shRNA construct targeting the NIb gene showed higher silencing efficiency. Northern blot and short interfering RNA (siRNA) analyses demonstrated that the viral resistance can be attributed to the degradation of the target RNA through the RNA silencing system. Correlation analysis of siRNA sequence characteristics with its activity suggested that the secondary structure stability of the antisense strand did not influence siRNA activity; 1 to 3 nt 5’ end of the sense strand caused a significant effect on siRNA activity where the first base such as U was favourable for silencing; the base mismatch between the siRNA and the target gene may be more tolerated in the 5’ end.     

Physicochemical characterization of anionic lipid-based ternary siRNA complexes

Physicochemical characterization is a useful tool in understanding lipoplex assemblies and their correlation to biological activity. Anionic lipid-based ternary siRNA complexes composed of anionic liposomes (DOPG/DOPE), calcium ions and siRNA, have recently been shown to be safe and efficient in a breast cancer cell culture model. In the present work, the effects of various formulation parameters such as liposome composition (DOPG/DOPE ratio) and anionic lipid/Ca(2+)/siRNA molar charge ratio, on the physicochemical attributes (particle size, surface charge, siRNA loading efficiency and serum stability) of these ternary anionic lipoplexes were evaluated. Particle size, siRNA loading efficiency and serum stability correlated with the in vitro silencing efficiency of these lipoplexes. For example, large lipoplex particles (5/2.5/1 anionic lipid/Ca(2+)/siRNA molar charge ratio) showed less efficient silencing while absolute serum stability and high siRNA loading (1.3/2.5/1 anionic lipid/Ca(2+)/siRNA molar charge ratio), exhibited maximum silencing in breast cancer cells. The physicochemical properties also indicated that the siRNA exists in the complexed and/or encapsulated form within the lipoplexes, depending on the anionic lipid/siRNA charge ratio. Based on these studies a model representing lipid-siRNA association within the anionic lipoplexes prepared under various formulation conditions is proposed. Physicochemical attributes can be utilized to estimate in vitro activity of lipid-siRNA complexes and understand their morphology.     

Designing highly active siRNAs for therapeutic applications

The discovery of RNA interference (RNAi) generated considerable interest in developing short interfering RNAs (siRNAs) for understanding basic biology and as the active agents in a new variety of therapeutics. Early studies showed that selecting an active siRNA was not as straightforward as simply picking a sequence on the target mRNA and synthesizing the siRNA complementary to that sequence. As interest in applying RNAi has increased, the methods for identifying active siRNA sequences have evolved from focusing on the simplicity of synthesis and purification, to identifying preferred target sequences and secondary structures, to predicting the thermodynamic stability of the siRNA. As more specific details of the RNAi mechanism have been defined, these have been incorporated into more complex siRNA selection algorithms, increasing the reliability of selecting active siRNAs against a single target. Ultimately, design of the best siRNA therapeutics will require design of the siRNA itself, in addition to design of the vehicle and other components necessary for it to function in vivo. In this minireview, we summarize the evolution of siRNA selection techniques with a particular focus on one issue of current importance to the field, how best to identify those siRNA sequences likely to have high activity. Approaches to designing active siRNAs through chemical and structural modifications will also be highlighted. As the understanding of how to control the activity and specificity of siRNAs improves, the potential utility of siRNAs as human therapeutics will concomitantly grow.     

Improved serum stability and biophysical properties of siRNAs following chemical modifications

Small interfering RNAs (siRNAs) which mediate sequence-specific gene silencing through RNA interference have been harnessed for therapeutic applications. Here, we found combinations of modest 2′-ribose siRNA modifications that permit effective gene silencing and improve siRNA stability in human serum. Introduction of as few as two different or identical 2′-ribose modifications at a specific position of either strand, which has been previously suggested to suppress unintended off-target effects, enhances the thermal and thermodynamic stability of siRNA duplexes, and does not alter the overall A-form helical geometry that is required for silencing activity. Our results provide a useful template for the design for therapeutic siRNAs. Ik Sang Cho and Jung Kim contributed equally to this work.