An efficient algorithm for systematic analysis of nucleotide strings suitable for siRNA design

<p>Abstract</p> <p>Background</p> <p>The "off-target" silencing effect hinders the development of siRNA-based therapeutic and research applications. Existing solutions for finding possible locations of siRNA seats within a large database of genes are either too slow, miss a portion of the targets, or are simply not designed to handle a very large number of queries. We propose a new approach that reduces the computational time as compared to existing techniques.</p> <p>Findings</p> <p>The proposed method employs tree-based storage in a form of a modified truncated suffix tree to sort all possible short string substrings within given set of strings (i.e. transcriptome). Using the new algorithm, we pre-computed a list of the best siRNA locations within each human gene ("siRNA seats"). siRNAs designed to reside within siRNA seats are less likely to hybridize off-target. These siRNA seats could be used as an input for the traditional "set-of-rules" type of siRNA designing software. The list of siRNA seats is available through a publicly available database located at <url>http://web.cos.gmu.edu/~gmanyam/siRNA_db/search.php</url></p> <p>Conclusions</p> <p>In attempt to perform top-down prediction of the human siRNA with minimized off-target hybridization, we developed an efficient algorithm that employs suffix tree based storage of the substrings. Applications of this approach are not limited to optimal siRNA design, but can also be useful for other tasks involving selection of the characteristic strings specific to individual genes. These strings could then be used as siRNA seats, as specific probes for gene expression studies by oligonucleotide-based microarrays, for the design of molecular beacon probes for Real-Time PCR and, generally, any type of PCR primers.</p

An RNAi in silico approach to find an optimal shRNA cocktail against HIV-1

<p>Abstract</p> <p>Background</p> <p>HIV-1 can be inhibited by RNA interference <it>in vitro </it>through the expression of short hairpin RNAs (shRNAs) that target conserved genome sequences. <it>In silico </it>shRNA design for HIV has lacked a detailed study of virus variability constituting a possible breaking point in a clinical setting. We designed shRNAs against HIV-1 considering the variability observed in naïve and drug-resistant isolates available at public databases.</p> <p>Methods</p> <p>A Bioperl-based algorithm was developed to automatically scan multiple sequence alignments of HIV, while evaluating the possibility of identifying dominant and subdominant viral variants that could be used as efficient silencing molecules. Student t-test and Bonferroni Dunn correction test were used to assess statistical significance of our findings.</p> <p>Results</p> <p>Our <it>in silico </it>approach identified the most common viral variants within highly conserved genome regions, with a calculated free energy of ≥ -6.6 kcal/mol. This is crucial for strand loading to RISC complex and for a predicted silencing efficiency score, which could be used in combination for achieving over 90% silencing. Resistant and naïve isolate variability revealed that the most frequent shRNA per region targets a maximum of 85% of viral sequences. Adding more divergent sequences maintained this percentage. Specific sequence features that have been found to be related with higher silencing efficiency were hardly accomplished in conserved regions, even when lower entropy values correlated with better scores. We identified a conserved region among most HIV-1 genomes, which meets as many sequence features for efficient silencing.</p> <p>Conclusions</p> <p>HIV-1 variability is an obstacle to achieving absolute silencing using shRNAs designed against a consensus sequence, mainly because there are many functional viral variants. Our shRNA cocktail could be truly effective at silencing dominant and subdominant naïve viral variants. Additionally, resistant isolates might be targeted under specific antiretroviral selective pressure, but in both cases these should be tested exhaustively prior to clinical use.</p

Comparisons between Chemical Mapping and Binding to Isoenergetic Oligonucleotide Microarrays Reveal Unexpected Patterns of Binding to the Bacillus subtilis RNase P RNA Specificity Domain†

ABSTRACT: Microarrays with isoenergetic pentamer and hexamer 20-O-methyl oligonucleotide probes with LNA (locked nucleic acid) and 2,6-diaminopurine substitutions were used to probe the binding sites on theRNase P RNA specificity domain of Bacillus subtilis. Unexpected binding patterns were revealed. Because of their enhanced binding free energies, isoenergetic probes can break short duplexes, merge adjacent loops, and/or induce refolding. This suggests new approaches to the rational design of short oligonucleotide therapeutics but limits the utility of microarrays for providing constraints for RNA structure determination. The microarray results are compared to results from chemical mapping experiments, which do provide constraints. Results from both types of experiments indicate that the RNase P RNA folds similarly in 1MNaþ and 10 mMMg2þ. Binding of RNA to RNA is important for many natural func-tions, includingproteinsynthesis (1,2), translationregulation (3,4), gene silencing (5, 6), metabolic regulation (7), RNAmodification (8, 9), etc. (10-13). Binding of oligonucleotides toRNAs is impor-tant for therapeutic approaches, such as siRNA, ribozymes, and antisense therapy (14, 15).Much remains to bediscovered, however, of the rules for predicting binding sites andpotential therapeutics