Skip to main content

VIRsiRNApred: a web server for predicting inhibition efficacy of siRNAs targeting human viruses

Abstract

Background

Selection of effective viral siRNA is an indispensable step in the development of siRNA based antiviral therapeutics. Despite immense potential, a viral siRNA efficacy prediction algorithm is still not available. Moreover, performances of the existing general mammalian siRNA efficacy predictors are not satisfactory for viral siRNAs. Therefore, we have developed “VIRsiRNApred” a support vector machine (SVM) based method for predicting the efficacy of viral siRNA.

Methods

In the present study, we have employed a new dataset of 1725 viral siRNAs with experimentally verified quantitative efficacies tested under heterogeneous experimental conditions and targeting as many as 37 important human viruses including HIV, Influenza, HCV, HBV, SARS etc. These siRNAs were divided into training (T1380) and validation (V345) datasets. Important siRNA sequence features including mono to penta nucleotide frequencies, binary pattern, thermodynamic properties and secondary structure were employed for model development.

Results

During 10-fold cross validation on T1380 using hybrid approach, we achieved a maximum Pearson Correlation Coefficient (PCC) of 0.55 between predicted and actual efficacy of viral siRNAs. On V345 independent dataset, our best model achieved a maximum correlation of 0.50 while existing general siRNA prediction methods showed PCC from 0.05 to 0.18. However, using leave one out cross validation PCC was improved to 0.58 and 0.55 on training and validation datasets respectively. SVM performed better than other machine learning techniques used like ANN, KNN and REP Tree.

Conclusion

VIRsiRNApred is the first algorithm for predicting inhibition efficacy of viral siRNAs which is developed using experimentally verified viral siRNAs. We hope this algorithm would be useful in predicting highly potent viral siRNA to aid siRNA based antiviral therapeutics development. The web server is freely available at http://crdd.osdd.net/servers/virsirnapred/.

Background

Viruses like Influenza, Hepatitis, Dengue, SARS coronavirus (SARS-CoV), Human Immunodeficiency Virus (HIV) etc. remain a public health concern worldwide due to their emerging and re-emerging nature [1]. Due to lack of therapeutics against majority of viruses, there is always a need to develop more effective antiviral agents [2]. Lately, RNA interference (RNAi) has emerged as a potential therapeutic tool for targeting human viruses. RNAi or gene silencing is a process by which sequence-specific degradation of mRNA takes place [3]. In this process, long dsRNA precursors are chopped into shorter (19–23 resides) units by a ribonuclease enzyme called dicer. These short interfering RNAs (siRNAs) possess two terminal nucleotide 3′ overhangs. After that a ribonucleoprotein machinery called RNA induced silencing complex (RISC) incorporates one of the siRNA strands and cleaves the complementary target mRNA using ATP [4].

Researchers have extensively used RNAi process to target a number of viral genes to suppress their expression [5, 6]. siRNAs targeting different regions of the HIV genome in infected cells showed promising results in inhibiting viral replication [7, 8]. Also, siRNAs targeting the influenza virus nucleocapsid and RNA transcriptase genes restricted its transcription and replication [9, 10]. Similarly, siRNAs directed against the Hepatitis B virus surface regions prevented the virus production [11]. siRNAs employed against SARS-CoV envelope and genes were able to effectively block their expression [12]. In another study, siRNAs targeting Dengue virus genes were able to impede the viral infection [13]. In addition, siRNAs have been shown to curb many other viruses like Human papillomavirus (HPV) [14], West Nile virus (WNV) [15] etc.

RNAi methodology has many desirable features to use as antiviral agents. It can target diverse types of viral genomes, whether it be double/single stranded DNA or RNA, which make it a suitable candidate for broad-spectrum antiviral therapy [16]. Also, siRNA aims at small length of the target mRNA instead of a functional domain of a protein, therefore, even a small viral genome can lend many targetable regions [5]. Further, many siRNAs may be expressed simultaneously to increase inhibition in a coordinated manner [17]. They are also harnessed to degrade mRNAs which generate disease causing proteins [6]. siRNA based drugs have also entered the clinical trials for various human diseases e.g. kidney disorders, LDL lowering, ocular/retinal disorders, cancer and viral diseases [18]. Specifically for Respiratory syncytial virus (RSV), ALN-RSV01 has completed Phase II trial. Yet another drug named SPC3649, developed by Santaris Pharma targets Hepatitis C virus (HCV) and was also under Phase II clinical trial [19, 20]. Many researchers have further reviewed the importance of RNAi in inhibiting viral infections [20, 21].

To predict effective siRNAs, a number of general siRNA design rules as well prediction algorithms been developed since every siRNA designed against a given mRNA are not equally effective [22, 23]. The earliest guidelines for siRNA design were based on presence and/or absence of specific nucleotide residues at different positions in siRNA as proposed by Elbashir [24], Reynolds [23], Ui-Tei [25], Amarzguioui [26] and Jagla [27]. Naito et al. have used the guidelines of Ui-Tei, Reynolds and Amarzguioui to select effective siRNAs against viruses [28].

Subsequently, machine learning techniques like SVM, ANN etc. have been utilized to predict effective mammalian siRNAs [29, 30] or their efficacy [31, 32]. Performance of these methods was better than the general siRNA design rules [33]. Saetrom used boosted genetic programming to predict mammalian siRNA efficacy on a modest dataset of 581 siRNAs [34]. Later Huesken et al. reported the screen method of functional siRNAs by using an artificial neural network [35]. Holen reported siRNA rules based on apparent overrepresentation or underrepresentation of certain nucleotides at different positions [36]. Researchers have utilised many siRNA features like nucleotide composition [37], thermodynamic parameters [38], nucleotide position [39] etc. to predict the efficiency of siRNAs. Many other workers have also used a combination of features using different machine learning techniques to predict mammalian siRNA efficacy [4042].

When we used best of these existing algorithms for prediction of viral siRNA, their performance was far from satisfactory. It could be because none of these methods were developed using virus specific siRNAs. Therefore, this prompted us to make virus specific siRNA efficacy prediction algorithm using experimentally validated viral siRNA dataset. Here we employed machine learning techniques to predict the inhibition efficacy of siRNAs. This will help the researchers to select the best siRNAs for use as potential therapeutics against important human viruses.

Methods

Algorithm development

Data collection

We have chosen our dataset from Viral siRNA database (VIRsiRNAdb) which contains over 1358 siRNA sequences targeting different human viruses and HIV siRNA database (HIVsirDB) having over 750 entries. From these databases we selected 927 and 240 sequences of 19mer which were having numerical (quantitative) efficacy. In addition 1204 sequences from patents were also used. In addition 67 more siRNAs were collected from the literature. This combined dataset of 2294 siRNAs was reduced to 1725 sequences after removal of redundant sequences. From this dataset we randomly selected 345 sequences for validation (V345) and the rest 1380 sequences were used in training (T1380).

siRNA sequence features

In total 12 different parameters were employed including four siRNA sequence features as nucleotide frequencies, binary pattern, thermodynamic properties, secondary structure features and their combinations or hybrids.

1. Nucleotide Frequencies

In the past many scientists have utilised nucleotide composition to predict siRNA potency. [31, 40, 42]. Nucleotide frequency is the number of each nucleotide in a siRNA. Since the length of our siRNAs was constant, we used nucleotide frequency instead of composition. The objective of calculating nucleotide frequencies of siRNA sequences is to transform any length of nucleotide sequence to fixed length feature vectors. This is important while using machine-learning techniques because it requires fixed length pattern. The information of each siRNA can be encapsulated to a vector of 4, 16, 64, 256, 1024 etc. multi-dimensions using frequencies of its mononucleotide, dinucleotide, triinucleotide, tetranucleotide and pentanucleotide sub sequences respectively.

2. Binary Pattern

Previous reports have used nucleotide positions to predict the RNAi activity of siRNAs [31, 39]. We employed binary pattern to extract siRNA features based on the occupancy of nucleotides at each position of siRNA sequences. Four binary patterns used for each nucleotide are follows A = 1000, C = 0100, G = 0010, U = 0001 and this resulted in generation of 76 patterns for each 19-mer siRNA.

3. Thermodynamic Properties

Thermodynamic stability of the siRNA has been considered by previous workers as an important feature in siRNA design [37]. They were calculated by the method of Khvorova et al. [43]. The thermodynamic dimensions correspond to the Gibbs free energy stability of the nucleotide pairs of the siRNAs. In total 21 features were used to calculate the energies of the different sets of interacting nucleotides. These features include the binding free energies and stabilities of the folded structures.

4. Secondary Structure

Secondary structure represents the capability of a single molecule of nucleic acid sequence to form intra molecular contacts, thereby stabilizing certain sequence parts as double stranded. Secondary structure was calculated using RNAfold programme of Vienna RNA package as implemented by Peek [31]. It predicts the minimum free energy (MFE) secondary structure and equilibrium base-pairing probabilities of single sequences [44]. The secondary structures predicted by RNAfold are depicted as brackets and dots indicating ‘paired’ and ‘unpaired’ nucleotides respectively. These categorical attributes are made readable to SVM by converting them into numeric data representing the structural features as vectors.

5. Hybrid Approaches

In hybrid approach, besides the features being used individually, combinations of siRNA sequence parameters were used in order to increase the performance of the prediction method as done by earlier researchers [27, 42]. We have used four hybrid methods notably mono-di-tri-tetra-penta-binary which makes a vector of total 1440 and mono-di-tri-tetra-penta-binary-thermo which makes a vector of total 1489 features.

6. Leave one out cross validation (LOOCV)

In this method each siRNA is kept for testing iteratively while remaining viral siRNAs are used for training the respective predictive models. Besides LOOCV, we have also carried out Leave one virus out cross validation (LOVOCV) strategy. In this approach, siRNAs from each virus are iteratively excluded and SVM is trained on the remaining virus siRNAs followed by testing on the excluded siRNAs of that individual virus.

7. Viral siRNA target conservation

Target site conservation analysis of all experimental 1725 viral siRNA sequences was done by matching each siRNA sequence with its respective reference viral genome sequences taken from NCBI. For this purpose, we have used ALIGN0 algorithm [45], which computes the alignment of two DNA sequences without penalizing for end-gaps.

Algorithm and server implementation

Support vector machines (SVMs) were trained with the selected sequence features to predict siRNA potency. SVM allows choosing a number of parameters and kernels The SVMlight software package (available at http://svmlight.joachims.org/) was used to construct SVM classifiers. In this study, we used the radial basis function (RBF) kernel:

k x ¯ , y ¯ = exp γ ǁ x ¯ y ¯ ǁ 2

where x̅ and y̅ are two data vectors, and γ is a training parameter.

Artificial Neural Network (ANN) was implemented using Stuttgart Neural Network Simulator (SNNS) package (available at http://www.ra.cs.uni-tuebingen.de/SNNS/) while K-Nearest Neighbour (KNN) and Reduced Error Pruning (REP) Tree algorithms were implemented using Weka machine learning software suite (available at http://www.cs.waikato.ac.nz/ml/weka/).

The server is implemented on Red Hat Linux and Apache (2.2.17) in back-end and front-end of web interface is implemented with PHP (5.2.14).

Validation

In order to evaluate performance of our models, we used Pearson’s correlation coefficient (R). All models were evaluated using ten-fold cross validation technique.

R = n n = 1 n E i act E i pred n = 1 n E i act n = 1 n E i pred n n = 1 n E i act 2 ( n = 1 n E i act ) 2 n n = 1 n E i pred 2 ( n = 1 n E i pred ) 2

Where n is the size of test set, Eipred and Eiact is the predicted and actual efficacy respectively.

Results

Performance evaluation of predictive models during 10-fold cross validation

Viral siRNA prediction modelshave been developed using various siRNA sequence features including mono, di, tri, tetra, penta nucleotide frequencies, binary pattern of nucleotides, thermodynamic properties and secondary structure. During 10-fold cross validation using SVM we achieved a maximum correlation of 0.19, 0.32, 0.42, 0.43, 0.46, 0.19, 0.26, 0.07 respectively for the above-mentioned individual features, between predicted and actual efficacy of viral siRNAs. We employed different features like nucleotide frequency (mono to penta) and positional (binary) as well as structural and thermodynamic features on viral siRNA training dataset consisting of 1380 (T1380) sequences. During 10-fold cross validation using SVM, the Pearson Correlation Coefficient (PCC) increased from 0.19 to 0.46 while moving from individual mono to pentanucleotide frequency features respectively. The binary feature performed similar to mononucleotide frequency. The thermodynamic properties showed marginal correlations with siRNA inhibition while secondary structural features did not work. Also, we found that hybrid models (frequency-binary and frequency-binary-thermo) performed better compared to other features used individually. Amongst the hybrid models, PCC was found to increase by combining frequency, binary and thermodynamic feature vectors (0.55) as compared to hybrid frequency (0.48), hybrid frequency-binary (0.50) or hybrid frequency-binary-thermo-secondary structural features (0.53) as shown in Table 1.

Table 1 Ten-fold cross validation performance of predictive models on viral siRNA dataset of 1380 sequences (T 1380 ) using SVM, ANN, KNN and REP Tree machine learning techniques

We also used other machine learning algorithms like ANN, KNN and REP Tree to check their performance on the experimental viral siRNA data using the above mentioned features. During 10-fold cross validation, PCC increased from 0.10 to a maximum of 0.30 from mono to pentanucleotide frequency in all the three cases. Here also the binary and secondary structural features did not perform well, confirming our earlier results using SVM. However while using a combination of features (frequency-binary and frequency-binary-thermo) in hybrid models, the correlation performance increased from 0.30 to 0.46 for ANN, 0.31 to 0.48 for KNN and 0.31 to 0.35 for REP Tree (Table 1). However SVM performed better compared to other machine learning algorithms in all the cases and was thus chosen for model development.

Performance evaluation of predictive models on independent 345 viral siRNA dataset

Besides 10-fold cross validation, we also checked the performance of our algorithms on independent dataset of 345 siRNAs (V345) on the above models as shown in the Table 2. Here again, the binary and secondary structural features did not perform well while the model based on thermodynamic features showed a marginal correlation of 0.19. As expected, the hybrid frequency, hybrid frequency-binary and hybrid frequency-binary-thermo-secondary structural features gave a better correlation with PCC values 0.48, 0.48 and 0.45 respectively. The best performing hybrid model with frequency-binary-thermo feature vectors gave a correlation of 0.50 and thus performed better on the validation dataset. Other machine learning techniques performed in a similar trend but their correlation was less as compared to SVM. Their best correlations were achievedwhile using the hybrid frequency-binary-thermo model. The PCCs were 0.33 for ANN, 0.35 for KNN and 0.31 for REP Tree (Table 2).

Table 2 Evaluation of performance of predictive models on validation dataset of 345 viral siRNAs(V 345 )

Comparison with existing siRNA prediction algorithms

We also compared our models with existing siRNA efficacy prediction servers (some servers are not currently working or have become outdated) although they are not optimized to predict viral siRNA efficacy and are rather used for general mammalian siRNA prediction. These prediction servers were mainly developed using two mammalian siRNA datasets of Saetrom (581 siRNA) [34] and Huesken (2431 siRNA) [35]. Saetrom dataset was tested under heterogeneous experimental conditions while Huesken dataset was tested under homogenous experimental conditions. On both of these datasets (19 mer), we have also developed earlier a ‘siRNApred’ algorithm with PCC of 0.56 and 0.68 respectively during 10-fold cross validation [46]. We checked performance of ‘siRNApred’ on viral siRNA datasets of V345, T1380and (V + T)1725 and observed PCC of 0.10, 0.16 and 0.14 respectively. Detailed comparisons with other methods are shown in Tables 3 and 4.

Table 3 Comparison of VIRsiRNApred with existing siRNA efficacy prediction algorithms developed using heterogeneous siRNA dataset
Table 4 Comparison of VIRsiRNApred with existing siRNA efficacy prediction methods developed using mammalian homogeneous siRNA dataset

The siRNA efficacy prediction methods developed using heterogeneous siRNA datasets are listed in Table 3. Their correlation performance during cross validation on the training set varied between 0.46 and 0.56 and for the test set their PCC was between 0.40 and 0.52 which shows the performance of these methods is almost similar. However, no siRNA efficacy prediction web server based on heterogeneous siRNA datasets was presently available/working except our ‘siRNApred’ which showed maximum PCC of 0.10 on V345 viral siRNA dataset.

The siRNA efficacy prediction algorithms developed using homogeneous siRNA datasets are given in Table 4. Their PCCs ranged from 0.66 to 0.85 on training during 10-fold cross validation and 0.55 to 0.71 on validation datasets. However on viral siRNA dataset V345 their performance dropped drastically in the range of 0.05 to 0.18. These results showed that the mammalian siRNA efficacy prediction methods developed using the homogeneous datasets performed better than those developed using the heterogeneous datasets. However, all such general mammalian siRNA efficacy prediction methods did not perform appropriately in predicting viral siRNA inhibition.

Performance evaluation using Leave one out cross validation (LOOCV)

We also checked the performance of our method using LOOCV technique by employing earlier mentioned siRNA features (Table 5). Here also the performance was increased while using mono to penta nucleotide frequency and further using hybrid features. The performance is improved modestly in comparison to 10-fold cross validation. Using this technique we achieved a maximum correlation of 0.58 on the training dataset T1380 and 0.55 on the validation dataset V345 using the hybrid model combining hybrid nucleotide frequency, binary and thermodynamic features.

Table 5 Performance of the SVM models using leave one out cross validation (LOOCV) method

Performance evaluation using Leave one virus out cross validation (LOVOCV)

To further check our SVM based method on the best performing features combining hybrid composition (mono to penta nucleotide frequency), binary and thermo properties for each virus in the 1725 viral siRNA dataset, we used LOVOCV method. The results are shown in Table 6. In the LOVOCV method, 252 siRNAs from Influenza A Virus were excluded and the model trained with remaining 1473 siRNAs from 36 viruses. It showed a PCC of 0.48 and 0.46 on the training and validation dataset respectively. Similar performance PCC is observed for different viruses like Influenza A Virus (0.48 and 0.46), HCV (0.51 and 0.44), SARS (0.53 and 0.48), Measles Virus (0.56 and 0.51), JEV (0.58 and 0.51) etc. on the training and validation dataset respectively. Overall, the training dataset performance during 10-fold cross validation ranged from PCC value of a minimum 0.43 to a maximum 0.58 with an average 0.54. While the validation performance ranged from PCC 0.40 to 0.53 with an average 0.48.

Table 6 Performance of the SVM model for each virus in the 1725 viral siRNA dataset using leave one virus out cross validation (LOVOCV) method

Viral siRNA target conservation

Results of the viral siRNA target site conservation among reference viral genome sequences is provided in the Additional file 1: Figure S2. The analysis shows the number of nucleotide differences or mismatches (0, 1, 2, 3, 4 and >4) between each siRNA and the respective reference viral genome sequences in the alignment using Align0 algorithm. Overall percentages of the 0, 1, 2, 3, 4 and >4 mismatches were 36.75, 15.64, 9.01, 8.08, 7.76 and 23.68% respectively as shown in Additional file 1: Figure S2 (a). We have also checked the conservation of 322 highly effective siRNAs with inhibition above 80% and an equal number of least effective siRNAs with inhibition less than 4.0% from the main dataset. The effective and ineffective siRNAs had similar number of 1, 2, 3 and 4 mismatches with more differences for 0 and >4 mismatches as shown in Additional file 1: Figure S2 (b).

Webserver

The web server is freely available via the url http://crdd.osdd.net/servers/virsirnapred. The strategy used to develop VIRsiRNApred models is shown in Figure 1. Besides, overview of VIRsiRNApred web server functionality including input and output interfaces are shown in Figure 2. To predict the best performing siRNAs to silence a specific viral gene, the user needs to paste a fasta sequence of the corresponding mRNA/gene region and choose the desired model and click submit. The sequence should not contain non-nucleotide characters or symbols. The output shows the recursive siRNAs chopped from the mRNA & their inhibition. We have also displayed useful links pointing to BLAST, alignment and off-target score for each siRNA. The user can sort the results in increasing/decreasing order. The output also shows graph depicting the inhibition and off-target scores in a pictorial manner (Figure 2).

Figure 1
figure 1

Workflow of the VIRsiRNApred model development.

Figure 2
figure 2

Web server and its functionality (top) submit page (bottom) result output.

Analysis tools

The analysis tools help the user to select the best antiviral siRNAs for the desired gene/region. First, the user needs to sort the siRNAs with “inhibition” and choose the siRNAs with high inhibition values. Then the user needs to look for the off-target using the two provided options: a) the user should select the siRNAs with least number of seeds/off-targets which are also shown pictorially by different colours as green (minimum), yellow (medium) and red (maximum) and b) the user can BLAST each siRNA sequence against human genome to ascertain that it does not affect other human genes. Besides, VIRsiRNAdb-BLAST helps the user to find similar siRNAs reported in the VIRsiRNAdb database.

Subsequently, to check the conservation of siRNA among different viral strains, our web server provides the “siTarConserve” tool. In this tool, user provided siRNA sequence is matched against selected reference viral genome sequences. There are options to search for conservation using individual virus as well as its family. The option to analyse conserved siRNA target regions using BLAST and Smith-Waterman algorithm has also been incoporated. This tool will be useful in selecting those siRNAs which are highly conserved in reference viral strains. Also while predicting siRNA efficacy, the user can check for siRNA conservation by clicking on the siTarConserve link provided against each siRNA on the prediction result page.

Discussion

Researchers working in the field of siRNA based therapeutics for viruses have generated a vast amount of data over the years. However, bioinformatics resources in the field were lacking and there was no viral siRNA efficacy prediction method available. In this direction, we have recently developed a comprehensive viral siRNA database “VIRsiRNAdb” [48] and another “HIVsirDB” [49] exclusively for HIV. Now we have developed VIRsiRNApred -a viral siRNA efficacy prediction algorithm.

Although many mammalian siRNA prediction algorithms have been developed in the past [33], these methods either classify a siRNA as effective/non-effective [29] or predict the inhibition efficacy of a siRNA [31, 32]. However, there is limited success in predicting siRNA efficacy due to limited size and diversity of available siRNA datasets [50].

Mammalian siRNA efficacy prediction methods were initially developed using siRNA tested under heterogeneous experimental conditions like Saetrom (581 siRNA), Shabalina (653 siRNA) or Holen (526 siRNA) etc. and achieved a PCC of 0.46-0.56 and subsequently methods were developed using siRNAs tested under homogenous experimental conditions like Heusken (2431 siRNA) and Katoch (702 siRNA) to achieve a PCC of 0.66-0.80. Though the mammalian siRNA prediction methods developed on homogenous siRNA datasets performed better than those developed on heterogeneous siRNA datasets, however all such methods did not perform well in predicting viral siRNA efficacy with PCC as low as 0.05-0.18.

We have developed our algorithm using a heterogeneous dataset of 1725 experimentally validated siRNA collected from 161 studies and targeting more than 200 mRNAs of over 40 diverse human viruses. Earlier Reynolds [23] has reported that 16.3% of randomly selected siRNAs have greater than 80% efficacy in targeting two mammalian genes (human cyclophilin B and firefly luciferase). Similarly, in 1725 viral siRNA dataset, only 17.1% were able to inhibit the virus above the 80% level despite majority of these siRNAs were designed using existing siRNA prediction methods before their testing in the wet-lab. This further verifies that the available mammalian siRNA efficacy prediction methods do not perform satisfactorily in predicting siRNAs targeting viruses.

Mammalian siRNA design rules are suitable for getting effective sequences for specific genes but they might sometimes prove unsuitable for selecting sequences for many other genes [51]. Therefore if these rules were used to design siRNAs for say, viral genes, many sequences might be wrongly selected as potential candidates. Even within mammalian siRNA prediction methods, the nucleotide positional preferences are not consistent (Additional file 1: Figure S1) with increase in heterogeneity of the data. Hence, designing functional siRNAs that target viral sequences is problematic because of their extraordinarily high genetic diversity. Therefore we need an algorithm to be developed on diverse viral siRNA dataset to tackle this issue.

However, in the past there have been attempts to predict viral siRNA. Naito et al. in 2006 published siVirus to classify a given viral sequence as siRNA based on earlier published guidelines of Ui-Tei et al., Reynolds et al. and Amarzguioui et al. for a few viruses [28]. Later, ElHefnawi in 2011 also provided a guideline to design siRNA against Influenza A virus using earlier published guidelines of Tuschl et al., Reynolds et al., Amarzguioui et al., Ui-Tei et al., and Hsieh et al., to select the siRNA [52]. However these methods just predict the viral siRNA as effective or non-effective based on earlier siRNA design rules. But they did not utilize any machine learning based methods which perform better than general siRNA design guidelines [33]. Moreover these methods were not developed using any experimental viral siRNA dataset and also do not predict quantitative efficacy of siRNA.

In the present study we made an effort to enhance the in silico prediction of highly effective human viral siRNA using manually curated datasets. Our model employed various siRNA sequence features particularly nucleotide frequencies, binary patterns, thermodynamic features, secondary structure and their hybrids as these parameters play a vital role in predicting efficient viral siRNAs. Mononucleotide composition can, for example, account for GC contents separately while higher nucleotide compositions can account for presence of any long GC stretch. Similarly trinucleotide composition can be useful in codon preference considerations. Studies have suggested that thermodynamic features of nucleotides at the beginning and at the end of the siRNA strands are important to their potency [53]. An increase in Pearson Correlation Coefficient (PCC) from 0.19 to 0.55 was observed using different approaches. The increase in performance could be attributed to the fact that as we increase the length of siRNA fragment (e.g., from 4 patterns in mono to 1024 patterns in pentanucleotide frequencies) more and more input information is available for processing not only in terms of composition but also other features like base stacking and positional effect to a significant extent. However secondary structural features of nucleotides showed minimal performance during training and were thus not included in the final model. Hence, the sequence frequencies combined with binary and thermodynamic features as hybrid approach we could achieve a maximum PCC of 0.58 during LOOCV. Thus, individual siRNA features provide marginal but cumulative increase in the probability of selecting a potent siRNA which is consistent with previous findings [32, 33].

The LOVOCV method showed similar performance for most of the viruses; however Performance PCC improved with increasing number of siRNAs in the training dataset. Since viruses are known to be genetically diverse [54] some variation in performance is expected, however, the hybrid features allow different siRNA features to complement each other during the training process, which eventually makes the efficacy prediction better than when the features are used individually. LOVOCV performances are similar to that of PCC (0.55 and 0.50) of 10-fold cross validation of T1380 training and V345 validation datasets using best hybrid siRNA features reported in Table 1. Despite exclusion of individual virus siRNAs from the training dataset, SVM model predicts siRNA appropriately for that particular virus. Therefore, VIRsiRNApred would function as a general viral siRNA efficacy predictor even for viruses currently not included in the training and testing. Besides, we also checked the performance of general mammalian siRNA predictors (siRNApred and siExplorer) on the above datasets and their performance PCC ranged from 0.10 to 0.16 only. This again accentuate that VIRsiRNApred perform competently for prediction viral siRNA efficacy.

Since viruses are known to be genetically diverse [54], some researchers have used conserved target sites to design siRNAs [28]. ElHefnawi has used this approach to design siRNAs against the influenza A virus [52]. Similarly Rosales et al. predicted siRNAs against NS4B and NS5 of Dengue virus [55] and Raza et al. predicted siRNAs against HA and NA genes of Influenza virus A based on sequence conservation and alignment [56]. Importance of selection of conserved regions targeted by siRNA in HIV-1 has been discussed by Naito et al. [28, 57].

In the overall conservation analysis of 1725 viral siRNAs with their respective reference genomes only around 37% were fully conserved (0 mismatch). It could be because of high viral genome heterogeneity. Besides that the researchers often target a particular viral strain and the target site may not be conserved among all other strains. However conservation analysis of highly and least effective siRNAs showed that number highly effective siRNAs are more in the fully conserved regions (0 mismatch) whereas the number of ineffective siRNAs is more in the least conserved regions (>4 mismatches). This suggests that siRNAs selected in the conserved regions tend to be more effective. However siRNAs designed from highly conserved virus genome regions as available in our experimental dataset display all ranges of siRNA inhibition efficacies ranging from 0-100%. This implies that not all siRNAs chosen from conserved regions are highly effective as many other factors like nucleotide frequencies, binary and thermodynamics features etc. contribute to the efficacy of a siRNA besides conservation as used in the VIRsiRNApred development. Nevertheless, selecting potential siRNAs from conserved genome regions is advantageous as it will target multiple viral strains.

This regression model can be useful in selecting siRNA molecules targeting viral genes for therapeutic purpose. Using our web-server researchers can target any given viral mRNA and get a list of highly performing siRNAs which have a greater chance to fully degrade the viral gene. In addition the ‘VIRsiRNApred’ server also offers an updated list of the experimentally validated siRNAs along with their numerical inhibition values reported in scientific literature that have been used against diverse types of viral genes. Also our web-server offers predicted siRNAs against important genes of HIV, Influenza, HCV, HBV and SARS. As siRNAs have been used against many viruses targeting a variety of genes and a few are already in different stages of clinical trials so they can serve as useful tools to develop potential antiviral drugs.

Currently the data availability for viral siRNAs comes from a wide range of studies due to which the overall performance of the model is affected. As more homogeneous data is generated from high-throughput studies, we will be interested in updating our prediction models in accordance with the new information to further improve the predictive power of our algorithm.

Conclusions

In this study we have reported the first viral siRNA efficacy prediction algorithm developed on experimentally verified viral siRNAs targeting as many as 37 diverse human viruses since existing general mammalian siRNA prediction methods are not able to effectively predict viral siRNA activity. VIRsiRNApred web-server will be helpful to select potent virus inhibiting siRNAs that can increase the knock down success rate and thus shorten the validation time in the development of antiviral siRNA therapeutics.

Abbreviations

RNAi:

RNA interference

siRNA:

Small interfering RNA

VIRsiRNApred:

Viral siRNA predictor

LOOCV:

Leave on out cross validation

LOVOCV:

Leave on virus out cross validation

SVM:

Support vector machine

RBF:

Radial basis function

ANN:

Artificial neural network

SNNS:

Stuttgart neural network simulator

KNN:

K-Nearest Neighbour

REP:

Reduced error pruning

V345:

Validation dataset

T1380:

Training dataset

siTarConserve:

siRNA Target conservation tool.

References

  1. Nichol ST, Arikawa J, Kawaoka Y: Emerging viral diseases. Proc Natl Acad Sci USA. 2000, 97: 12411-12412.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  2. Qureshi A, Thakur N, Kumar M: HIPdb: a database of experimentally validated HIV inhibiting peptides. PLoS One. 2013, 8: e54908-

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  3. 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: 806-811.

    Article  CAS  PubMed  Google Scholar 

  4. Filipowicz W: RNAi: the nuts and bolts of the RISC machine. Cell. 2005, 122: 17-20.

    Article  CAS  PubMed  Google Scholar 

  5. Leonard JN, Schaffer DV: Antiviral RNAi therapy: emerging approaches for hitting a moving target. Gene Ther. 2006, 13: 532-540.

    Article  CAS  PubMed  Google Scholar 

  6. Haasnoot J, Berkhout B: Nucleic acids-based therapeutics in the battle against pathogenic viruses. Handb Exp Pharmacol. 2009, 189: 243-263.

    Article  CAS  PubMed  Google Scholar 

  7. Singh SK: RNA interference and its therapeutic potential against HIV infection. Expert Opin Biol Ther. 2008, 8: 449-461.

    Article  CAS  PubMed  Google Scholar 

  8. Rossi JJ: RNAi as a treatment for HIV-1 infection. Biotechniques. 2006, 40: 25-29.

    Article  Google Scholar 

  9. Ge Q, McManus MT, Nguyen T, Shen CH, Sharp PA, Eisen HN, Chen J: RNA interference of influenza virus production by directly targeting mRNA for degradation and indirectly inhibiting all viral RNA transcription. Proc Natl Acad Sci USA. 2003, 100: 2718-2723.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  10. Nakazawa M, Kadowaki SE, Watanabe I, Kadowaki Y, Takei M, Fukuda H: PA subunit of RNA polymerase as a promising target for anti-influenza virus agents. Antiviral Res. 2008, 78: 194-201.

    Article  CAS  PubMed  Google Scholar 

  11. Wu KL, Zhang X, Zhang J, Yang Y, Mu YX, Liu M, Lu L, Li Y, Zhu Y, Wu J: Inhibition of Hepatitis B virus gene expression by single and dual small interfering RNA treatment. Virus Res. 2005, 112: 100-107.

    Article  CAS  PubMed  Google Scholar 

  12. Meng B, Lui YW, Meng S, Cao C, Hu Y: Identification of effective siRNA blocking the expression of SARS viral envelope E and RDRP genes. Mol Biotechnol. 2006, 33: 141-148.

    Article  CAS  PubMed  Google Scholar 

  13. Stein DA, Perry ST, Buck MD, Oehmen CS, Fischer MA, Poore E, Smith JL, Lancaster AM, Hirsch AJ, Slifka MK: Inhibition of dengue virus infections in cell cultures and in AG129 mice by a small interfering RNA targeting a highly conserved sequence. J Virol. 2011, 85: 10154-10166.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  14. Bousarghin L, Touze A, Gaud G, Iochmann S, Alvarez E, Reverdiau P, Gaitan J, Jourdan ML, Sizaret PY, Coursaget PL: Inhibition of cervical cancer cell growth by human papillomavirus virus-like particles packaged with human papillomavirus oncoprotein short hairpin RNAs. Mol Cancer Ther. 2009, 8: 357-365.

    Article  CAS  PubMed  Google Scholar 

  15. Kumar P, Lee SK, Shankar P, Manjunath N: A single siRNA suppresses fatal encephalitis induced by two different flaviviruses. PLoS Med. 2006, 3: e96-

    Article  PubMed Central  PubMed  Google Scholar 

  16. Haasnoot J, Westerhout EM, Berkhout B: RNA interference against viruses: strike and counterstrike. Nat Biotechnol. 2007, 25: 1435-1443.

    Article  CAS  PubMed  Google Scholar 

  17. Chen Y, Mahato RI: siRNA pool targeting different sites of human hepatitis B surface antigen efficiently inhibits HBV infection. J Drug Target. 2008, 16: 140-148.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  18. Liu G, Wong-Staal F, Li QX: Development of new RNAi therapeutics. Histol Histopathol. 2007, 22: 211-217.

    CAS  PubMed  Google Scholar 

  19. Haussecker D: The Business of RNAi Therapeutics in 2012. Mol Ther Nucleic Acids. 2012, 1: e8-

    Article  PubMed Central  PubMed  Google Scholar 

  20. Shah PS, Schaffer DV: Antiviral RNAi: translating science towards therapeutic success. Pharm Res. 2011, 28: 2966-2982.

    Article  CAS  PubMed  Google Scholar 

  21. Umbach JL, Cullen BR: The role of RNAi and microRNAs in animal virus replication and antiviral immunity. Genes Dev. 2009, 23: 1151-1164.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  22. Takasaki S, Kawamura Y, Konagaya A: Selecting effective siRNA sequences by using radial basis function network and decision tree learning. BMC Bioinformatics. 2006, 7 (5): S22-

    Article  PubMed Central  PubMed  Google Scholar 

  23. Reynolds A, Leake D, Boese Q, Scaringe S, Marshall WS, Khvorova A: Rational siRNA design for RNA interference. Nat Biotechnol. 2004, 22: 326-330.

    Article  CAS  PubMed  Google Scholar 

  24. 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: 494-498.

    Article  CAS  PubMed  Google Scholar 

  25. 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.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  26. Amarzguioui M, Prydz H: An algorithm for selection of functional siRNA sequences. Biochem Biophys Res Commun. 2004, 316: 1050-1058.

    Article  CAS  PubMed  Google Scholar 

  27. Jagla B, Aulner N, Kelly PD, Song D, Volchuk A, Zatorski A, Shum D, Mayer T, De Angelis DA, Ouerfelli O: Sequence characteristics of functional siRNAs. RNA. 2005, 11: 864-872.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  28. Naito Y, Ui-Tei K, Nishikawa T, Takebe Y, Saigo K: siVirus: web-based antiviral siRNA design software for highly divergent viral sequences. Nucleic Acids Res. 2006, 34: W448-W450.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  29. Teramoto R, Aoki M, Kimura T, Kanaoka M: Prediction of siRNA functionality using generalized string kernel and support vector machine. FEBS Lett. 2005, 579: 2878-2882.

    Article  CAS  PubMed  Google Scholar 

  30. Mysara M, Garibaldi JM, Elhefnawi M: MysiRNA-designer: a workflow for efficient siRNA design. PLoS One. 2011, 6: e25642-

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  31. Peek AS: Improving model predictions for RNA interference activities that use support vector machine regression by combining and filtering features. BMC Bioinformatics. 2007, 8: 182-

    Article  PubMed Central  PubMed  Google Scholar 

  32. Vert JP, Foveau N, Lajaunie C, Vandenbrouck Y: An accurate and interpretable model for siRNA efficacy prediction. BMC Bioinforma. 2006, 7: 520-

    Article  Google Scholar 

  33. McQuisten KA, Peek AS: Comparing artificial neural networks, general linear models and support vector machines in building predictive models for small interfering RNAs. PLoS One. 2009, 4: e7522-

    Article  PubMed Central  PubMed  Google Scholar 

  34. Saetrom P: Predicting the efficacy of short oligonucleotides in antisense and RNAi experiments with boosted genetic programming. Bioinformatics. 2004, 20: 3055-3063.

    Article  CAS  PubMed  Google Scholar 

  35. Huesken D, Lange J, Mickanin C, Weiler J, Asselbergs F, Warner J, Meloon B, Engel S, Rosenberg A, Cohen D: Design of a genome-wide siRNA library using an artificial neural network. Nat Biotechnol. 2005, 23: 995-1001.

    Article  CAS  PubMed  Google Scholar 

  36. Holen T: Efficient prediction of siRNAs with siRNArules 1.0: an open-source JAVA approach to siRNA algorithms. RNA. 2006, 12: 1620-1625.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  37. Shabalina SA, Spiridonov AN, Ogurtsov AY: Computational models with thermodynamic and composition features improve siRNA design. BMC Bioinforma. 2006, 7: 65-

    Article  Google Scholar 

  38. Ichihara M, Murakumo Y, Masuda A, Matsuura T, Asai N, Jijiwa M, Ishida M, Shinmi J, Yatsuya H, Qiao S: Thermodynamic instability of siRNA duplex is a prerequisite for dependable prediction of siRNA activities. Nucleic Acids Res. 2007, 35: e123-

    Article  PubMed Central  PubMed  Google Scholar 

  39. Katoh T, Suzuki T: Specific residues at every third position of siRNA shape its efficient RNAi activity. Nucleic Acids Res. 2007, 35: e27-

    Article  PubMed Central  PubMed  Google Scholar 

  40. Pan WJ, Chen CW, Chu YW: siPRED: predicting siRNA efficacy using various characteristic methods. PLoS One. 2011, 6: e27602-

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  41. Filhol O, Ciais D, Lajaunie C, Charbonnier P, Foveau N, Vert JP, Vandenbrouck Y: DSIR: assessing the design of highly potent siRNA by testing a set of cancer-relevant target genes. PLoS One. 2012, 7: e48057-

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  42. Sciabola S, Cao Q, Orozco M, Faustino I, Stanton RV: Improved nucleic acid descriptors for siRNA efficacy prediction. Nucleic Acids Res. 2013, 41: 1383-1394.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  43. Khvorova A, Reynolds A, Jayasena SD: Functional siRNAs and miRNAs exhibit strand bias. Cell. 2003, 115: 209-216.

    Article  CAS  PubMed  Google Scholar 

  44. Hofacker IL: Vienna RNA secondary structure server. Nucleic Acids Res. 2003, 31: 3429-3431.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  45. Myers EW, Miller W: Optimal alignments in linear space. Comput Appl Biosci. 1988, 4: 11-17.

    CAS  PubMed  Google Scholar 

  46. Kumar M, Lata S, Raghava GPS: siRNApred: SVM based method for predicting efficacy value of siRNA. Proceedings of the first international conference on Open Source for Computer Aided Drug Discovery (OSCADD). 2009, Chandigarh: CSIR-IMTECH

    Google Scholar 

  47. Jiang P, Wu H, Da Y, Sang F, Wei J, Sun X, Lu Z: RFRCDB-siRNA: improved design of siRNAs by random forest regression model coupled with database searching. Comput Methods Programs Biomed. 2007, 87: 230-238.

    Article  PubMed  Google Scholar 

  48. Thakur N, Qureshi A, Kumar M: VIRsiRNAdb: a curated database of experimentally validated viral siRNA/shRNA. Nucleic Acids Res. 2012, 40: D230-D236.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  49. Tyagi A, Ahmed F, Thakur N, Sharma A, Raghava GP, Kumar M: HIVsirDB: a database of HIV inhibiting siRNAs. PLoS One. 2011, 6: e25917-

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  50. Liu Q, Xu Q, Zheng VW, Xue H, Cao Z, Yang Q: Multi-task learning for cross-platform siRNA efficacy prediction: an in-silico study. BMC Bioinformatics. 2010, 11: 181-

    Article  PubMed Central  PubMed  Google Scholar 

  51. Takasaki S: Efficient prediction methods for selecting effective siRNA sequences. Comput Biol Med. 2010, 40: 149-158.

    Article  CAS  PubMed  Google Scholar 

  52. ElHefnawi M, Hassan N, Kamar M, Siam R, Remoli AL, El-Azab I, AlAidy O, Marsili G, Sgarbanti M: The design of optimal therapeutic small interfering RNA molecules targeting diverse strains of influenza A virus. Bioinformatics. 2011, 27: 3364-3370.

    Article  CAS  PubMed  Google Scholar 

  53. Liu Q, Zhou H, Zhu R, Xu Y, Cao Z: Reconsideration of in silico siRNA design from a perspective of heterogeneous data integration: problems and solutions. Brief Bioinform. 2012, Epub ahead of print

    Google Scholar 

  54. Shackelton LA, Parrish CR, Holmes EC: Evolutionary basis of codon usage and nucleotide composition bias in vertebrate DNA viruses. J Mol Evol. 2006, 62: 551-563.

    Article  CAS  PubMed  Google Scholar 

  55. 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: 519-522.

    Article  PubMed Central  PubMed  Google Scholar 

  56. 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: 340-343.

    Article  PubMed Central  PubMed  Google Scholar 

  57. Naito Y, Nohtomi K, Onogi T, Uenishi R, Ui-Tei K, Saigo K, Takebe Y: Optimal design and validation of antiviral siRNA for targeting HIV-1. Retrovirology. 2007, 4: 80-

    Article  PubMed Central  PubMed  Google Scholar 

Download references

Acknowledgments

Authors are thankful to Council of Scientific and Industrial Research (CSIR), (OLP0083), Department of Biotechnology (GAP001) and Institute of Microbial Technology for financially supporting the work.

Author information

Authors and Affiliations

Authors

Corresponding author

Correspondence to Manoj Kumar.

Additional information

Competing interests

The authors declare that they have no competing interests.

Authors’ contributions

MK conceived the approach, helped in analysis and interpretation of data, gave overall supervision to the project. MK, AQ wrote the manuscript. AQ, NT collected the data, implemented machine learning software and developed the web server. All of the authors read and approved the final manuscript.

Abid Qureshi, Nishant Thakur contributed equally to this work.

Electronic supplementary material

12967_2013_1691_MOESM1_ESM.doc

Additional file 1: Figure S1: Comparison of highly and least effective siRNAs using two sample logo. Figure S2. Viral siRNA target conservation. (DOC 126 KB)

Authors’ original submitted files for images

Below are the links to the authors’ original submitted files for images.

Authors’ original file for figure 1

Authors’ original file for figure 2

Rights and permissions

Open Access This article is published under license to BioMed Central Ltd. This is an Open Access article is distributed under the terms of the Creative Commons Attribution License ( https://creativecommons.org/licenses/by/2.0 ), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly credited. The Creative Commons Public Domain Dedication waiver ( https://creativecommons.org/publicdomain/zero/1.0/ ) applies to the data made available in this article, unless otherwise stated.

Reprints and permissions

About this article

Cite this article

Qureshi, A., Thakur, N. & Kumar, M. VIRsiRNApred: a web server for predicting inhibition efficacy of siRNAs targeting human viruses. J Transl Med 11, 305 (2013). https://doi.org/10.1186/1479-5876-11-305

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/1479-5876-11-305

Keywords