- Open Access
SiRNA silencing efficacy prediction based on a deep architecture
© The Author(s). 2018
- Published: 24 September 2018
Small interfering RNA (siRNA) can be used to post-transcriptional gene regulation by knocking down targeted genes. In functional genomics, biomedical research and cancer therapeutics, siRNA design is a critical research topic. Various computational algorithms have been developed to select the most effective siRNA, whereas the efficacy prediction accuracy is not so satisfactory. Many existing computational methods are based on feature engineering, which may lead to biased and incomplete features. Deep learning utilizes non-linear mapping operations to detect potential feature pattern and has been considered perform better than existing machine learning method.
In this paper, to further improve the prediction accuracy and facilitate gene functional studies, we developed a new powerful siRNA efficacy predictor based on a deep architecture. First, we extracted hidden feature patterns from two modalities, including sequence context features and thermodynamic property. Then, we constructed a deep architecture to implement the prediction. On the available largest siRNA database, the performance of our proposed method was measured with 0.725 PCC and 0.903 AUC value. The comparative experiment showed that our proposed architecture outperformed several siRNA prediction methods.
The results demonstrate that our deep architecture is stable and efficient to predict siRNA silencing efficacy. The method could help select candidate siRNA for targeted mRNA, and further promote the development of RNA interference.
- Deep learning
In 1988, Fire first introduced RNA interference (RNAi) [1–3], and now it has been found that this mechanism can be detected in many eukaryotic systems, such as mammals, fungi, plants and invertebrates . Small interfering RNA (siRNA) is the production of RNAi, which can induce instant target gene knockdown . RNAi is a vital tool for researching gene function [5–7] and can be used as an effective therapeutic method in the treatment of virus and cancer [8–10].
The gene silencing efficacy of RNAi relies on siRNA design, and many efforts are being made in this area. In early days, several sets of empirical rules to select effective siRNA were proposed according to experimental data. These rules are mainly based on GC content , base preferences at specific positions [12, 13], thermodynamic stability , internal structure  and target mRNA secondary structure . However, these rules are summarized from small scale dataset and can hardly reach our acceptable level. With the accumulation of validated siRNAs, machine learning has been used in effective siRNA recognition. ‘Biopredsi’ is a classical siRNA efficacy prediction method, which are based upon artificial neural network algorithm . Besides, a major siRNA dataset was supplied by Huesken et al. The dataset includes 2431 siRNAs, which were built by high-throughput analysis technology. It has been truly admitted that this dataset is very helpful for the construction of other siRNA efficacy prediction methods . As another artificial neural network, ThermoComposition-21  includes both composition and thermodynamic features. The simple linear method was also used in this area. The method proposed by Jean-Philippe Vert  used two kinds of siRNA sequence features as feature set. One is the nucleotides present at each position in the siRNA sequence, the other is the global content of the siRNA in short motifs. It is an accurate and easily interpretable model, and according the experimental results the prediction accuracy of Biopredsi is as accurate as it. Another linear regression model was constructed by nucleotide preference scores .
The siRNA efficacy prediction accuracy cannot make us satisfied though the considerable efforts. The reason is the prediction results of most machine algorithms are highly dependent upon the siRNA features, including sequence feature, thermodynamic feature, secondary structure feature, etc. Most of these features are biased and incomplete feature vectors since they are produced by the traditional feature engineering way which is reliant on expert knowledge, and the prediction ability will be limited. Recently, a frontier machine learning algorithm, deep learning, has aroused the attention of researchers. It has been proved that deep learning performed better than the existing machine learning method. Different from the traditional machine learning methods, deep learning framework can conduct the prediction in a data-driven way.
In this paper, we constructed a new siRNA efficacy model based on deep learning algorithm. Firstly, we extracted hidden feature patterns from two modalities, including sequence context features and thermodynamic property. Then we merge them to implement the prediction. For the sequence context features, we utilized convolution layers to automatically learn motif encoding features. In the convolution layers, convolution kernels can be seen as motif detectors, and the potential feature pattern of siRNA multimode motif can be automatically learned by a data-driven method. This method is more abstract and more conductive to prediction and more closely to the essence. The experimental results showed that our deep architecture performed better than the current siRNA efficacy prediction methods in terms of prediction accuracy.
For siRNA efficacy prediction, we collected 4067 siRNA samples from the dataset of Huesken(2431) , Reynolds(248) , Vickers (80) , Haborth(44) , Takayuki(702) , Ui-Tei (62)  and siRNAdb(500) .
In this paper, we divided these siRNA sequences into two datasets by random partition, a training dataset (3660) and a testing dataset (407).
Encoding of siRNA sequences
There are two encoding method in our paper to transform the siRNA sequences into quantized biological descriptors.
Sequence context features
The research which has been reported showed that the sequence context outside the target region effected the efficacy of siRNA . In this paper, the 21 + 2n sequence of n upstream and downstream flanking nucleotides around binding region together with the targeted sites were intercepted.
When the length of the flanking region is less than n, the corresponding positions will be encoded to 0.05. The encoding method maps the sequence to a sparse coding and quantifies nucleotides according their relative position.
The thermodynamic properties for siRNA efficacy prediction
Number of features
Stability of hybridization formed between siRNA and mRNA
Differential thermodynamic stability of siRNA duplex ends
thermodynamic parameter of every two base pairs along the siRNA duplex antisense strand
The deep model construction
In our deep architecture, there are a convolutional layer and a pooling layer in view of the sample size and computational complexity. In the convolutional layer, there are multiple convolution kernels, which have different sizes. These convolution kernels can be seen as motif detectors, which can help us find the motifs playing important role in siRNA efficacy prediction. And through the convolutional operation we can get the corresponding motif encoding features. Most of the existing siRNA efficacy prediction methods encode the sequence according to the experience, but our deep architecture is different from them. The features of our method are trained by siRNA datasets. The feature extraction method has more information, guidance quality and usability. Then the pooling layer can select the most representative motif feature pattern as the feature representation.
The design of motif detector
To explore the potential feature pattern included in the siRNA sequence, we designed various convolution kernels. The large-scale training samples are used to correct the weights of convolution kernels by back-propagation algorithm, which guarantee that we can obtain the effective feature pattern.
In this formula, S represents the sequence of flanking nucleotides around binding region together with the targeted sites and M is the m × 4 convolution kernel. xk is the neuron of convolutional layer (1 ≤ k ≤ 22-m), and δk is the learning rate for correcting weights. The convolution result is a (22-m) × 1 matrix, which represents the feature pattern of every multimode motif.
Because there are various convolution kernels in the convolution layer, the output of pooling layer is 2d-dimensional vector, where d is the number of convolution kernels.
Assessment of the prediction system
To assess the model efficacy, we adopted two indices, including Pearson Correlation Coefficient (PCC) and the area under the ROC curve (AUC).
AUC is used extensively to measure the overall performance of prediction model. A higher PCC and AUC indicate the model performs well.
In this section, we will interpret our experimental results of different parameters. In every experiment, 10-fold cross-validation is conducted to obtain the best parameters.
The influence of the length of flanking nucleotides on prediction result
Figure 2 showed when n equals to 20, the prediction result achieves the best performance. The results indicated that our model needs more sequence information to detect more useful deep features.
The infuence of hyper parameters on prediction result
This part mainly discussed the influence of different hyper parameters on prediction result. In our deep architecture, there are three hyper parameters directly affecting the model’s robustness and deciding the structure of network, including the size of convolution kernel, activation function and learning rate. There comparative experiments were conducted to search the optimal hyper parameters for our deep architecture.
The size of convolution kernel
As shown in Fig. 3, different convolution kernels influence the prediction results. When m equals to 15, the prediction result achieves the best performance. The result indicates that the convolution kernels we designed could learn the effective feature pattern from the input modality. Next, we analyze the effect of m on prediction result. Figure 3 shows that when m is increasing, the prediction result of corresponding deep neural network becomes larger, but when m is larger than 15, the prediction result lower. We speculate that the reason could be the convolution kernels with smaller size only discover the information associated with low-mode motif and neglect the contribution of high-mode motif and the whole sequence feature. And when the size of convolution kernel is becoming increasing, the feature representation of high-mode motif will be detected and the contribution of low-mode motif will be neglected. The result in either case can give rise to the decrease of prediction result. Therefore, we should choose the reasonable size of convolution kernel to achieve effective motif feature learning. In this paper, we designed a convolution kernel set, and PCC of the convolution kernels contained in the set are higher than 0.6, which guarantee the learning feature pattern has adequate discriminating ability.
From Fig. 4 we can see that the prediction results with different combinations of activation function are diverse and the first combination performs best (ReLU + Sigmoid). The result shows that in the convolution layer the better choice is ReLU function. The reason may be that sparsity is added into the output feature of convolution layer by ReLU, and this way can enhance the nonzero neurons’ information. In the DNN layer, sigmoid can be used as activation function because it can summarize the contribution of all feature representation. Besides, output of sigmoid is from 0 to 1, which is consistent with the range of siRNA efficacy.
Figure 5 shows that when the learning rate equals to 0.1, we can get the best result. Besides, it can be found that when the learning rate equals to 0.5, the result is lowest. It shows that the neural network has lost the optimal weights and sank into relative extremum. Then when the learning rate equals to the other two values, PCC and AUC are relatively low. The reason may be that when the iterative is 1000 the network has slower convergence speed and cannot get the best weight. According to the prediction accuracy and training time, our deep architecture set the learning rate to 0.1.
Compared with other algorithms
From what has discussed above, our deep architecture has 15 convolution kernels with the size from 6 × 4 to 20 × 4. Through the convolution operation, we got 15 feature maps with size (22-m) × 1, each of which then was processed by max pooling and average pooling with (22-m) pooling size respectively in the following pooling stage. Thus, after such pooling operation, each input was transformed into 2 × 15 × 1 vector. The 15 kernels with different size are transformed into a 30-dimensional vector in the pooling layer. The activate function of convolutional layer is ReLU, the activation function of DNN layer is sigmoid, and the learning rate is 0.1. There are 25 nerouns in DNN layer.
From Fig. 6, it can be found out that our deep architecture performed best, reaching at 0.725 PCC and 0.903 AUC.
The most probable reason is that Biopredsi, DSIR and siRNApred are the traditional machine learning methods, which belong to the feature engineering way and rely on expert knowledge. And our deep learning methods can supply non-linear mapping operations and multiple layer networks to detect potential complex patterns and generate homogenous deep representations for prediction tasks. Therefore the performance of Biopredsi, DSIR and siRNApred are less than our deep architecture.
And we can find that the performance of our method is better than CNN. The method CNN used the feature of siRNA sequence and developed a convolutional neural network including a convolution layer and a pooling layer. Because the sequence cannot fully reflect the siRNA properties and the efficacy of siRNA strongly depends on the thermodynamic stability profile of the siRNA duplex, we used siRNA context feature and thermodynamic properties and added a DNN layer to combine the two types of feature. Because their components depict the feature of siRNA from different points of view, the fully connected DNN structure could interconnect all factors for their joint effect in its hidden states.
As a common molecular tool, siRNA can research gene function and be used as an effective therapeutic method in the treatment disease. Numerous methods have been developed to design active siRNA. However, the siRNA efficacy prediction accuracy cannot make us satisfied. In this study, we proposed a new siRNA efficacy prediction method based on a deep architecture. Comparing with the existing method Biopredsi, DSIR, siRNApred and CNN, our method performs best. The results show that our deep architecture could tap the contribution of siRNA context sequence and thermodynamic properties on efficacy prediction. Besides, our method can extract the valuable information contained in the feature pattern. Finally, the data-driven feature learning pattern outweighs the learning pattern which mainly depends on the expert knowledge.
This work was supported by the National Natural Science Foundation of China (11372155, 61802057), the open project program of Key Laboratory of Symbolic Computation and Knowledge Engineering of Ministry of Education, Jilin University under Grant No.93K172018K02, the National Key Research and Development Program of China(2017YFD0502001) and the Fundamental Research Funds for Industrial Innovation of Jilin Province, China (2018C039-3).
About this supplement
This article has been published as part of BMC Genomics Volume 19 Supplement 7, 2018: Selected articles from the IEEE BIBM International Conference on Bioinformatics & Biomedicine (BIBM) 2017: genomics. The full contents of the supplement are available online at https://bmcgenomics.biomedcentral.com/articles/supplements/volume-19-supplement-7.
YH initiated and designed the study. FH and HY conducted the data analysis. YH drafted the manuscript. YC and YL participated in experiment design, result interpretation and manuscript preparation. All authors read and approved the final manuscript.
Ethics approval and consent to participate
Consent for publication
The authors consent to publish this paper to BMC Genomics.
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Open AccessThis article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.
- Timmons L, Fire A. Specific interference by ingested dsRNA. Nature. 1998;395(6705):854.View ArticleGoogle Scholar
- Montgomery MK, Xu S, Fire A. RNA as a target of double-stranded RNA-mediated genetic interference in Caenorhabditis elegans. Proc Natl Acad Sci U S A. 1998;95(26):15502–7.View ArticleGoogle Scholar
- Elbashir SM. Duplexes of 21-nucleotide RNAs mediate RNA interference in cultured mammalian cells. Nature. 2001;411(6836):494–8.View ArticleGoogle Scholar
- Novina CD, Sharp PA. The RNAi revolution. Nature. 2004;430(6996):161–4.View ArticleGoogle Scholar
- Baulcombe DC. RNA as a target and an initiator of post-transcriptional gene silencing in trangenic plants. Plant Mol Biol. 1996;32(1–2):79.View ArticleGoogle Scholar
- Cogoni C, Irelan JT, Schumacher M, Schmidhauser TJ, Selker EU, Macino G. Transgene silencing of the al-1 gene in vegetative cells of Neurospora is mediated by a cytoplasmic effector and does not depend on DNA-DNA interactions or DNA methylation. EMBO J. 1996;15(12):3153–63.View ArticleGoogle Scholar
- 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(6669):806.View ArticleGoogle Scholar
- 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(24):3364–70.View ArticleGoogle Scholar
- Sharp PA. siRNA-directed inhibition of HIV-1 infection. Nat Med. 2002;8(7):681–6.View ArticleGoogle Scholar
- Resnier P, Montier T, Mathieu V, Benoit JP, Passirani C. A review of the current status of siRNA nanomedicines in the treatment of cancer. Biomaterials. 2013;34(27):6429.View ArticleGoogle Scholar
- Elbashir SM, Harborth J, Weber K, Tuschl T. Analysis of gene function in somatic mammalian cells using small interfering RNAs. Methods. 2002;26(2):199–213.View ArticleGoogle Scholar
- Reynolds A, Leake D, Boese Q, Scaringe S, Marshall WS, Khvorova A. Rational siRNA design for RNA interference. Nat Biotechnol. 2004;22(3):326–30.View ArticleGoogle Scholar
- Amarzguioui M, Prydz H. An algorithm for selection of functional siRNA sequences. Biochem Biophys Res Commun. 2004;316(4):1050–8.View ArticleGoogle Scholar
- Khvorova A, Reynolds A, Jayasena SD. Functional siRNAs and miRNAs exhibit strand bias. Cell. 2003;115(2):209–16.View ArticleGoogle Scholar
- Uitei K, Naito Y, Takahashi F, Haraguchi T, Ohkihamazaki 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(3):936.View ArticleGoogle Scholar
- Schubert S, Grünweller A, Erdmann VA, Kurreck J. Local RNA target structure influences siRNA efficacy: systematic analysis of intentionally designed binding regions. J Mol Biol. 2005;348(4):883.View ArticleGoogle Scholar
- 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(23):995–1001.View ArticleGoogle Scholar
- He F, Han Y, Gong J, Song J, Wang H, Li Y. Predicting siRNA efficacy based on multiple selective siRNA representations and their combination at score level. Sci Rep. 2017;7:44836.View ArticleGoogle Scholar
- Shabalina SA, Spiridonov AN, Ogurtsov AY. Computational models with thermodynamic and composition features improve siRNA design. BMC Bioinformatics. 2006;7View ArticleGoogle Scholar
- Vert JP, Foveau N, Lajaunie C, Vandenbrouck Y. An accurate and interpretable model for siRNA efficacy prediction. BMC Bioinformatics. 2006;7(1):520.View ArticleGoogle Scholar
- 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. 2006, 35(18):e123.View ArticleGoogle Scholar
- Vickers TA, Koo S, Bennett CF, Crooke ST, Dean NM, Baker BF. Efficient reduction of target RNAs by small interfering RNA and RNase H-dependent antisense agents. A comparative analysis. J Biol Chem. 2003;278(9):7108.View ArticleGoogle Scholar
- 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(2):83–105.View ArticleGoogle Scholar
- Katoh T, Suzuki T. Specific residues at every third position of siRNA shape its efficient RNAi activity. Nucleic Acids Res. 2007;35(4):e27.View ArticleGoogle Scholar
- Uitei K, Naito Y, Takahashi F, Haraguchi T, Ohkihamazaki 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(3):936–48.View ArticleGoogle Scholar
- Chalk AM, Warfinge RE, Georgiihemming P, Sonnhammer EL. siRNAdb: a database of siRNA sequences. Nucleic Acids Res. 2005;33(Database issue):D131–4.View ArticleGoogle Scholar
- Liu L, Li QZ, Lin H, Zuo YC. The effect of regions flanking target site on siRNA potency. Genomics. 2013;102(4):215–22.View ArticleGoogle Scholar
- Xia T, Jr SLJ, 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 base pairs. Biochemistry. 1998;37(42):14719.View ArticleGoogle Scholar
- Han Y, Liu Y, Zhang H, He F, Shu C, Dong L. Utilizing selected Di- and trinucleotides of siRNA to predict RNAi activity. Comput Math Methods Med. 2017;8:5043984.Google Scholar
- Han Y, He F, Tan X, Yu H. Effective small interfering RNA design based on convolutional neural network. IEEE Int Conf Bioinform Biomed. 2017;2017:16–21.Google Scholar