Skip to main content

RNA-binding residues prediction using structural features



RNA-protein complexes play an essential role in many biological processes. To explore potential functions of RNA-protein complexes, it’s important to identify RNA-binding residues in proteins.


In this work, we propose a set of new structural features for RNA-binding residue prediction. A set of template patches are first extracted from RNA-binding interfaces. To construct structural features for a residue, we compare its surrounding patches with each template patch and use the accumulated distances as its structural features. These new features provide sufficient structural information of surrounding surface of a residue and they can be used to measure the structural similarity between the surface surrounding two residues. The new structural features, together with other sequence features, are used to predict RNA-binding residues using ensemble learning technique.


The experimental results reveal the effectiveness of the proposed structural features. In addition, the clustering results on template patches exhibit distinct structural patterns of RNA-binding sites, although the sequences of template patches in the same cluster are not conserved. We speculate that RNAs may have structure preferences when binding with proteins.


RNA-protein complexes play an essential role in many biological processes, such as RNA splicing [1], translation [2, 3] and post-transcriptional gene regulation [4, 5], etc. Many large ribonucleoproteins (e.g. ribosome) are also RNA-protein complexes. In addition, certain proteins carry out specific functions such as repairing damaged RNAs [6] and editing transcribed RNAs [7]. Currently, there are 1,542 RNA-binding proteins in the human body which have been manually curated [8].

To understand the molecular mechanisms of the protein-RNA recognition, it is important to identify RNA-binding residues from target proteins. Obviously, the identification is straightforward if structures of RNA-protein complexes have been known in advance. However, it is expensive and time consuming to determine the structure of an RNA-protein complex through biological experiments. Considering this, people are resorting to computational methods which can quickly and accurately predict RNA-binding residues.

In previous studies, RNA-binding residues prediction was mostly based on sequence features of proteins. It is because that protein sequences are much easier to obtain than their structures in RNA-protein complexes. For sequence-based methods, the commonly used features include position-specific scoring matrix (PSSM) [914], hydrophobicity [15], electrostatics [16, 17], side chain environment [13, 18, 19], residue interface propensity [1921], and residue accessibility [12, 15].

In the last decades, abundant of 3D structures of RNA-protein complexes are emerging. According to the records in Protein Data Bank [22], there are only 491 RNA-protein complexes deposited before 2010. By August 2014, the number has been dramatically increased to 1560. As a consequence, researchers begin to explore new structural features which can improve the accuracy of RNA-binding residue prediction.

Compared with the efforts spent on exploring sequence features, studies focused on structural features are rather limited [23]. Kim et al. [24] proposed a new structural feature, namely residue doublet interface propensity. It describes the pairing preference of amino acids in protein-RNA interfaces. Chen and Lim [16] predicted RNA-binding residues based on irregular surface patches and clefts on the target proteins. The irregular surfaces and clefts were composed by the most conserved and electrostatically stabilized residues. Therefore, these surface patches and clefts could provide useful RNA-binding information and consequently improve the prediction accuracy. In [25], spatial patches and topological patches on protein surfaces, which were represented as contact graphs, were used to predict RNA-binding residues for target proteins. Instead of proposing new structural features, Towfic et al. [26] adopted several existing structural features, which included surface roughness of a residue [27] and the CX value [28]. Another method, DRNA [29], took a different strategy of using structural information to predict RNA-binding sites. It first aligned the structure of the target protein to the template proteins from an RNA-protein complex library and then predicted RNA-binding residues based on results of structure alignment. There were also some works focusing on utilizing secondary structure elements [30].

In the framework of RNA-binding residue prediction, a protein residue can be represented by a feature vector which is composed by a set of sequence and/or structural features. For protein residues whose binding status have been determined, their feature vectors can be constructed and used to train classification models. RNA binding propensity or status of a target residue can be predicted by the well trained classification models. Currently, a number of classification models have been adopted in RNA-binding residue prediction. One popular classification model is support vector machine (SVM) which has been adopted by several powerful servers such as BindN+, PiRaNhA, etc. [1214, 18, 3133]. Other popular classifiers include Naïve Bayes classifier [34, 35] and neural network [30, 36, 37]. In addition to using a single classifier, ensemble classification models which take the advantage of multiple classifiers are also exploited and they have been shown to greatly improve the prediction accuracy [15]. Besides using sequence/structural features, certain methods predict RNA-binding residues directly from amino acid sequences [34, 38].

Protein-RNA interaction interfaces are reported to be composed by several clusters of positive charged residues scattered on protein surfaces [17, 39]. Besides positive charged residues on the protein surface, protein-RNA interfaces also comprise binding pockets/cavities [16, 19]. We assume that RNA-binding sites are assembled by certain patches with specific shapes, which are regarded as binding units. Based on this assumption, we compile a set of template patches which participate in protein-RNA interactions. The templates are extracted from surfaces around RNA-binding residues from the training set. They are further grouped into several clusters. Representative patches, which are centers of clusters, are identified. Based on the representative patches, a set of new structural features can be constructed for a residue. We first extract surface patches around the residue. Then we compare its surface patches with each representative patch and use the accumulated distances as structural features. The dimension of structural features is the same as the number of representative patches. These new features provide sufficient structural information of surrounding surface of a residue. The structural similarity between the environments of two residues can be transformed into the similarity between their structural features. If distances between a target residue and RNA-binding residues are small, it is more likely to bind RNA molecules. If not, it tends to be a non-RNA-binding residue. These new features, together with other features such as PSSM and residue propensity, are used to predict RNA-binding residues. Using machine learning techniques, specifically the ensemble learning technique, the combined features exhibit a good discrimination power for RNA-binding residue prediction.

Results and discussion

Experimental results

Proteins in the dataset are divided into four groups and a 4-fold cross validation is adopted to estimate the prediction performance using the new structural features as well as other sequence features (e.g. PSSM scores, residue interface propensity). In each fold, one group is used as the test set and the other three are combined as the training set. An ensemble classifier is trained using feature vectors and class labels (binding or non-binding) of residues from the training set. RNA-binding scores are computed for residues in the test set using their feature vectors by the well-trained classifier.

In addition to evaluating the prediction performance using all features (i.e. the new structural features, PSSM scores, and interface propensity), we evaluate the prediction performance using different combinations of features, i.e. using structural features, PSSM scores, and the combined features which include structural features and PSSM scores. The prediction performances using different combinations of features are shown in Table 1. The performance using PSSM scores is slightly better than using structural features in terms of area under the curve (AUC) and Matthews correlation coefficient (MCC) (AUCPSSM = 0.64 and MCCPSSM = 0.19; AUCStructuralFeatures = 0.62 and MCCStructuralFeatures = 0.18). However, based on some other performance metrics, structural features outperform PSSM scores. For example, the precision for structural features is 0.46 while for PSSM scores it is 0.41.

Table 1 Prediction performances using different features

When combining two types of features together, the performance can be greatly improved. Compared with using PSSM scores, the AUC value of using combined features is increased from 0.64 to 0.67. The MCC value is increased from 0.19 to 0.24. When the feature of residue interface propensity is introduced, the prediction performance is slightly improved. The AUC value is increased from 0.67 to 0.68. The MCC value is increased from 0.24 to 0.26. Precision, accuracy, and specificity are also increased from 0.46/0.66/0.71 to 0.48/0.68/0.76, respectively. However, the values of sensitivity and F-score are decreased.

We also evaluated our method on different categories of proteins. The corresponding scores of different performance metrics are shown in Table 2. When evaluated by AUC scores, our method achieves better performances on aptermer, ribosomal, and small classes, on which AUC values are larger than 0.68. When evaluated using MCC values, our method achieves better performances on aptermer and splicing classes, on which MCC values are larger than 0.26.

Table 2 Prediction performances on different protein categories

We select four types of proteins to show the prediction results of our method. Proteins 4J1D:D, 1ASY:A, 1FXL:A, and 3 V24:V belong to viral, tRNA, mRNA, and ribosomal classes, respectively. Figure 1 shows the prediction results of our method for four proteins. Residues of true positive, false positive, false negative and true negative are shown in red, green, yellow, and grey, respectively. The numbers of true positives, false positives, false negatives and true negatives in each protein is also given in the caption of Fig. 1. It can be seen that, most RNA-binding sites in the above four proteins are successfully identified by our method.

Fig. 1

Examples of RNA-binding residue prediction results. RNA-binding residue prediction for (a) 4J1G:D (viral), (b) 1ASY:A (tRNA), (c) 1FXL:A (mRNA), and (d) 3 V24:V (ribosomal RNA). True positive, false positive, false negative and true negative residues are shown in red, green, yellow and grey, respectively. (a) The numbers of tp, fp, fn, and tn in 4J1G:D are 38, 21, 73, and 95, respectively; (b) the numbers of tp, fp, fn, and tn in 1ASY:A are 61, 89, 89, and 251, respectively; (c) the numbers of tp, fp, fn, and tn in 1FXL:A are 85, 79, 0, and 3, respectively; (d) the numbers of tp, fp, fn, and tn 3 V24:V are 34, 17, 0, and 2, respectively

Comparison with other methods

To show the effectiveness of our method, we make a comparison with another three publicly available methods (BindN+ [13], PPRInt [14], and DRNA [29]). These methods showed better performances over other seven methods evaluated by Puton et al. [23]. We evaluate their performances on a new dataset, namely RB344 of PRIDB dataset [40]. RB344 is a non-redundant dataset which is much larger than the dataset used in [23]. The prediction performances of different methods evaluated using proteins from RB344 dataset are listed in Table 3. Receiver operating characteristic (ROC) curves are shown in Fig. 2.

Table 3 Performance comparison
Fig. 2

Performance comparison using ROC curves. ROC curves of prediction results of our method, BindN+, and PPRInt

The AUC values of BindN+, PPRInt, and our method are the same, which is 0.68. When evaluated using MCC metric, BindN+ shows the best performance with its value 0.28. Our method and PPRInt have the same MCC value 0.26, which is a bit worse than BindN+. DRNA achieves the worst performance with its value of 0.22. BindN+ has the highest precision and F-score and DRNA has the highest accuracy and specificity. Our method has the highest sensitivity.


RNA-binding residue prediction is essential for understanding mechanisms of RNA recognition for proteins. In this work, we propose a set of new structural features for RNA-binding residues prediction. Specifically, we construct a set of 3-aa and 2-aa template patches around binding residues and group them into 60 clusters (40 clusters of 3-aa patches and 20 clusters of 2-aa patches). Representative patches which are centers of 60 clusters are identified and used for construct structural features. For each residue, its surrounding patches are extracted. Accumulated distances of surrounding patches to representative patches are computed and comprise a set of structural features. These new structural features, together with other features such as PSSM score and residue interface propensity are used to describe the RNA-binding specificity of the target residue. To accurately predict RNA-binding residues, our method adopts ensemble learning technique whose effectiveness in solving classification problems has been shown. The experimental results reveal that our method achieves a high performance on a benchmark test.

The clustering on 3-aa and 2-aa template patches show distinct structural patterns of RNA-binding sites. The clustering results indicate that RNAs have structure preference when binding with proteins. Currently these template patches are confined to 3-aa and 2-aa patches. In the future, we will extract larger template patches from binding sites and explore their structural patterns.



The dataset used in our experiments is RB344 of PRIDB dataset [40]. RB344 is a non-redundant dataset which contains 344 proteins belonging to 13 categories: RNAse, SRP, aptamer, dsRNA, exosome, mRNA, ribosomal, small, snRNP, splicing, tRNA, viral, and other. Global sequence alignment was performed on the dataset using the needle function provided by the emboss suite [41]. The sequence identities in RB344 are smaller than 30 %. RNA-binding residues were determined using two definitions: (i) a residue whose any atom is within a 5 Å distance of any atom in a nucleotide; and (ii) residues involved in van der Waals, hydrogen-bonding, hydrophobic or electrostatic interactions with nucleotides [40]. Any amino acid residue satisfying the above definitions are regarded as RNA-binding residue. RNA-protein complexes in the dataset are shown in Table 4.

Table 4 Protein-RNA complexes in RB344 dataset

Identification of protein surface residues

To determine protein surface residues, accessible areas will be computed first. If its accessible area is larger than zero, the residue is considered as a surface residue. Otherwise it is a non-surface residue. The accessible areas can be calculated using VMD software [42] with the probe radius of 1.4 Å.

Shape descriptor for protein residues

The backbone of an amino acid is defined by four atoms: N, CA, C, and O. The center of the side-chain is defined as the mean of coordinates of all heavy atoms on the side-chain. However, the side-chain of glycine only has a hydrogen substituent. Therefore, the hydrogen is used as the side-chain center of glycine. An amino acid is represented by atoms from its backbone and the point of its side-chain center.

Template patches construction

The neighbors of a protein residue are defined as amino acids which are located within a certain distance (d N ) on the protein surface. The distance between two amino acids is the smallest distance between their atoms. A patch is called 2-aa/3-aa patch if it is composed by two/three residues, respectively. Suppose a residue has a set of neighbors {n 1, n 2, …, n k }. When k ≥ 2, two neighbors can be selected from the neighbor set and construct a 3-aa patch with the residue. The total number of 3-aa patches constructed from the neighbor set is C2 k. When k = 1, a 2-aa patch can be constructed which is composed by the residue and its neighbor. The situation of k = 0 is not considered currently because we assume that interaction interfaces are areas consisting of two or more residues. If a surface residue in the training set is known to interact with RNAs, the patches constructed from its neighbors are regarded as positive patches. Only positive patches are used as template patches. We obtained 175,989 3-aa template patches and 122 2-aa template patches from RB344 dataset when d N equals 3 Å.

Structural similarity between patches

Because each amino acid can be represented as a set of atoms from the backbone and the center of its side chain, a 3-aa/2-aa patch can be represented as the assembly of representative points from all its member residues. When comparing shapes of two surface patches of the same size, i.e. both of them are 3-aa or 2-aa patches, they are treated as rigid objects. The structural similarity between two patches can be measured by the sum of Euclidean distances between the corresponding points after rotation and translation (i.e. the least-squares distance between two sets of points). Suppose a patch contains m (m{2,3}) residues, each of which is composed by n points, the least-squares distance between patch X and Y can be computed using Eq. (1).

$$ {d}_{LS}\left(X,Y\right)=\sqrt{\underset{s,R,t}{ \min }{\displaystyle \sum_{i=1}^m{\displaystyle \sum_{j=1}^n}{\left\Vert {x}_{ij}-\left(sR\left({y}_{ij}\right)+t\right)\right\Vert}^2}} $$

where s is a scale factor, R is a rotation matrix, and t is a translation vector. x ij and y ij are the j-th point from i-th residue of X and Y respectively. The optimal solution of s, R, and t for Eq. (1) is:

$$ R=V{U}^T,\;s=\sqrt{\frac{{\displaystyle \sum_{i=1}^n}{\left\Vert x{\hbox{'}}_i\right\Vert}^2}{{\displaystyle \sum_{i=1}^n}{\left\Vert y{\hbox{'}}_i\right\Vert}^2}},\;t=\overline{x}-sR\left(\overline{y}\right) $$

In Eq. (2), \( \overline{x} \) and \( \overline{y} \) are the centroids of X and Y. Matrices U and V are obtained by singular value decomposition: Y'X'T = UΣV T, where X' = {x' ij } i=1,…,m; j=1,…,n and Y' = {y' ij } i=1,…,m; j=1,…,n are obtained by subtracting \( \overline{x} \) and \( \overline{y} \) from the points, i.e. \( x{\hbox{'}}_{ij}={x}_{ij}-\overline{x} \) and \( y{\hbox{'}}_{ij}={y}_{ij}-\overline{y} \); i = 1,…,m and j = 1,…,n. Details of this optimal solution can be found in [43]. In our problem, we assume that there is no scale change between two similar patches. Therefore, the scale factor s is set to 1 instead of using the value in Eq. (2).

To compute the least-squares distance, the correspondence between points from two objects should be known in advance. However, the correspondence between the residues from two patches has not been determined yet. Therefore, orders of residues in patch X are permuted to create different correspondences to residues in patch Y. Once the correspondence between residues from two patches has been determined, the correspondence between their representative points will be automatically determined. To compare two 3-aa patches, there are 6 ways (P3 3) of correspondences and to compare two 2-aa patches, there are 2 ways (P2 2) of correspondences. We compute the least-squares distances between two patches using different correspondences and the minimum one is defined as the structural similarity (d SS ) between two patches. d SS can be computed using Eq. (3) according to its definition.

$$ {d}_{SS}\left(X,Y\right)=\underset{i\in \left\{1,2,\cdots, {P}_m^m\right\}}{ \min }{d}_{LS}\left({X}^{(i)},Y\right) $$

In Eq. (3), X (i) is the i-th way of reordering residues in patch X and there are Pm m ways (m is the number of residues in a patch) of residue reordering in all.

Clustering template patches

It’s difficult to use all template patches to construct structural features. Therefore, we select some representative ones from them so that the dimension of structural features can be acceptable. We group the extracted 3-aa and 2-aa template patches using complete-linkage hierarchical clustering. The distance metric used in clustering algorithm is the least-squares distance shown in Eq. (3). The cluster dendrograms of 3-aa patches and 2-aa patches extracted from the training set in a fold of cross validation are shown in Fig. 3. Hierarchical clustering using single- and average-linkage is also performed. However, the resulting dendrograms have ladder shapes. It indicates that these two methods are not suitable for clustering template patches. The final clusters represent distinct structural patterns of template patches. They are more or less similar to protein structural motifs but are much smaller. They can be regarded as binding units of interaction interfaces of proteins and are used to describe RNA-binding surfaces.

Fig. 3

Hierarchical clustering on 3-aa and 2-aa template patches. Hierarchical clustering with complete-linkage on (a) 3-aa and (b) 2-aa template patches

In each fold of cross validation, there are ~130,000 3-aa template patches constructed from the training set. We randomly selected 10,000 3-aa template patches and performed hierarchical clustering. The selected 3-aa patches and all 2-aa patches are further grouped into 40 and 20 clusters. In each cluster, the centroid patch, which has the smallest sum-of-square distance to other members, is also determined. The centroid patches are regarded as the representative patches.

Patches in each cluster reveal distinct structural patterns. For example, in one cluster, three amino acids are arranged in a linear way (see Fig. 4(a)). While in another, they are placed like the head of a fork (see Fig. 4(b)). The sequences of patches in each cluster are not conserved. However, their structures are quite similar. It indicates that template patches have specific structural patterns and RNAs may have structure preference when binding with proteins.

Fig. 4

Structures from two clusters of 3-aa template patches. (a) Four template patches in Cluster 1; (b) four template patches in Cluster 4

The new structural features construction

Given a residue r, structural features can be constructed in the following way. Firstly, all neighbours of r located within the distance of d N on the protein surface are identified. Suppose there are k neighbouring residues and they are denoted by {n 1, n 2,…, n k }. A set of patches {X 1,X 2,…,X K } can be constructed using {n 1, n 2,…, n k } and r: if k = 1, i.e. r has only one neighbour n 1, a 2-aa patch {X 1} which is simply composed of r and n 1 can be constructed. If k ≥ 2, several 3-aa patches which are composed of r and two of its neighbours can be constructed (the total number of 3-aa patches, K, equals to \( {C}_k^2 \)). To construct structural features for residue r, {X 1,X 2,…,X K } will be compared with each representative patch and accumulated distances to each representative patch will be computed.

Denote the set of representative patches by {Y 1,Y 2,…,Y L }. The length of the new feature vector L is the total number of 3-aa and 2-aa representative patches. In our work, L is 60 (because there are 40 3-aa representative patches and 20 2-aa representative patches). Rearrange representative patches and let {Y 1,Y 2,…,Y 40} be 3-aa representative patches and {Y 41,Y 42,…,Y 60} be 2-aa representative patches.

If k = 1, there will be only one 2-aa patch {X 1} surrounding r. X 1 will be compared with each 2-aa representative patch. Suppose f j is the distance of X 1 to Y j (40 < j ≤ L). Then

$$ {f}_j=\left\{\begin{array}{cc}\hfill {d}_{SS}\left({X}_1,{Y}_j\right)\hfill & \hfill 40<j\le L\hfill \\ {}\hfill 0\hfill & \hfill 1\le j\le 40\hfill \end{array}\right. $$

Because X 1 only contains two residues, it cannot be compared with 3-aa representative patches. The distances between X 1 and 3-aa representative patches are directly set zeros.

If k ≥ 2, each 3-aa patch X i (i = 1,…,K) is compared with Y j (1 ≤ j ≤ 40). f j can be computed using Eq. (5).

$$ {f}_j=\left\{\begin{array}{cc}\hfill {\displaystyle \sum_{i=1}^K{d}_{SS}\left({X}_i,{Y}_j\right)}\hfill & \hfill 1\le j\le 40\hfill \\ {}\hfill 0\hfill & \hfill 40<j\le L\hfill \end{array}\right. $$

f j (j = 1,…,40) is the accumulated distance of surrounding patches {X 1, X 2,…, X K } to the representative patch Y j . When 40 < j ≤ 60, f j is set zeroes because X i is a 3-aa patch which cannot be compared with 2-aa representative patches. In the end, a 60-dimension feature vector [f 1,…,f L ] can be constructed for the residue r.

The rationale of comparing {X 1,X 2,…, X K } with representative patches is as follows. The protein surface around a binding residue can be characterized by template patches. After clustering, template patches can be approximated by representative patches. Therefore, we can describe the protein surface surrounding a binding residue using the combination of representative patches. The problem is how to quantitatively measure the structural similarity of surfaces surrounding two residues. Considering that surfaces can be approximated by combinations of representative patches, we compute the accumulated distance of surrounding patches to each representative patch and denote it as a structural feature. If there are L representative patches, L features will be obtained. These structural features contain potential structural information. It can be seen that, for all residues, no matter RNA-binding or non-RNA-binding, their structural features can be constructed by computing the accumulated distance of surrounding patches to representative patches. Given a target residue r, if its surrounding surface is similar to the surfaces surrounding RNA-binding residues, its structural features will be more close to features of RNA-binding residues. Based on its structural features, r can be classified as an RNA-binding residue or a non-RNA-binding residue.

Other features used for RNA-binding residue prediction

In addition to the proposed structural features, other sequence features of amino acids are also introduced to describe RNA-binding property. Each residue in RNA-binding proteins is characterized by another two descriptors including: (i) PSSM which gives values of sequence conservation for residues using PSI-BLAST [44]; (ii) the residue interface propensity which describes the frequency of different types of amino acids occurring in the interaction interface than on the protein surface [23].

For residues from the training set and the test set, we can construct feature vectors which combine the new structural features and two additional sequence features. The dimension of all features is 81.

RNA-binding residue prediction using ensemble method

RNA-binding residue prediction can be regarded as a classification problem when feature vectors have been presented. In the learning process, a classification model can be learned using feature vectors and class labels of residues from the training set. Then, the classification model can be applied to predict binding propensities for residues in the test set. Compared with individual classifiers, ensemble classifiers have already been shown to produce better classification results [45, 46]. Specifically, in the problem of RNA-binding residue prediction, random forest, an ensemble classifier, has already been adopted and showed a high performance.

In our method, ensemble learning technique is also used. ENTOOL [47] is a package which integrates a series of classification algorithms, which include SVM, decision tree, ridge regression, Gaussian mixture models, multilayer perceptron, etc. In our work, models of ridge regression, perceptron, and multilayer perceptron are selected as constituent classifiers because they can achieve better performances than other classifiers in ENTOOL.

Methods for prediction performance evaluation

ENTOOL first performs five-fold cross-validation on the training residues to adjust parameters of the ensemble classifier and then predicts binding scores for target residues using the trained models. The predicted binding scores vary from −1 to 1. The larger the binding score, the higher binding propensity of the target residue.

By comparing the predicted scores with the true labels of those residues in the test set, four metrics can be computed: true positives (TP), true negatives (TN), false positives (FP), and false negatives (FN). Based on the four metrics, false positive rate (FPR) and true positive rate (TPR, which is also called sensitivity) can be computed (see Eq. (6)). ROC curve can be created by plotting FPR values against TPR values. Other performance metrics, such as AUC, accuracy, precision, specificity, F-score, and MCC can also be computed (see Eq. (6)).

$$ \begin{array}{l}\mathrm{T}\mathrm{P}\mathrm{R}=\mathrm{sensitivity}=\frac{\mathrm{TP}}{\mathrm{TP}+\mathrm{F}\mathrm{N}}\;\\ {}\mathrm{F}\mathrm{P}\mathrm{R}=\frac{\mathrm{FP}}{\mathrm{TN}+\mathrm{F}\mathrm{P}}\\ {}\mathrm{accuracy}=\frac{\mathrm{TP}+\mathrm{T}\mathrm{N}}{\mathrm{TP}+\mathrm{F}\mathrm{P}+\mathrm{F}\mathrm{N}+\mathrm{T}\mathrm{N}}\\ {}\mathrm{precision}=\frac{\mathrm{TP}}{\mathrm{TP}+\mathrm{F}\mathrm{P}}\\ {}\mathrm{specificity}=\frac{\mathrm{TN}}{\mathrm{TN}+\mathrm{F}\mathrm{P}}\\ {}\mathrm{F}\hbox{-} \mathrm{score}=\frac{2\mathrm{T}\mathrm{P}}{2\mathrm{T}\mathrm{P}+\mathrm{F}\mathrm{P}+\mathrm{F}\mathrm{N}}\\ {}\mathrm{M}\mathrm{C}\mathrm{C}=\frac{\mathrm{TP}\times \mathrm{T}\mathrm{N}-\mathrm{F}\mathrm{P}\times \mathrm{F}\mathrm{N}}{\sqrt{\left(\mathrm{T}\mathrm{P}+\mathrm{F}\mathrm{N}\right)\left(\mathrm{T}\mathrm{P}+\mathrm{F}\mathrm{P}\right)\left(\mathrm{T}\mathrm{N}+\mathrm{F}\mathrm{P}\right)\left(\mathrm{T}\mathrm{N}+\mathrm{F}\mathrm{N}\right)}}\end{array} $$



Area under the curve


False negatives


False positives


False positive rate


Matthews correlation coefficient


Position-specific scoring matrix


Receiver operating characteristic


Support vector machine


True negatives


True positives


True positive rate


  1. 1.

    Jurica MS, Moore MJ. Pre-mRNA splicing: awash in a sea of proteins. Mol Cell. 2003;12:5–14.

    CAS  Article  PubMed  Google Scholar 

  2. 2.

    Noller HF. RNA structure: reading the ribosome. Science. 2005;309:1508–14.

    CAS  Article  PubMed  Google Scholar 

  3. 3.

    Galicia-Vázquez G, Lindqvist L, Wang X, Harvey I, Liu J, Pelletier J. High-throughput assays probing protein–RNA interactions of eukaryotic translation initiation factors. Analytical Biochemistry. 2009;384:180–8.

    Article  CAS  PubMed  Google Scholar 

  4. 4.

    Glisovic T, Bachorik JL, Yong J, Dreyfuss G. RNA-binding proteins and post-transcriptional gene regulation. FEBS Lett. 2008;582:1977–86.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  5. 5.

    Keene JD. RNA regulons: coordination of post-transcriptional events. Nat Rev Genet. 2007;8:533–43.

    CAS  Article  PubMed  Google Scholar 

  6. 6.

    Aas PA, Otterlei M, Falnes PO, Vagbo CB, Skorpen F, Akbari M, et al. Human and bacterial oxidative demethylases repair alkylation damage in both RNA and DNA. Nature. 2003;421:859–63.

  7. 7.

    Bock R. Sense from nonsense: how the genetic information of chloroplasts is altered by RNA editing. Biochimie. 2000;82:549–57.

    CAS  Article  PubMed  Google Scholar 

  8. 8.

    Gerstberger S, Hafner M, Tuschl T. A census of human RNA-binding proteins. Nature Reviews Genetics. 2014;15:829–45.

    CAS  Article  PubMed  Google Scholar 

  9. 9.

    Ahmad S, Sarai A. PSSM-based prediction of DNA binding sites in proteins. BMC Bioinformatics. 2005;6:33.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  10. 10.

    Kuznetsov IB, Gou Z, Li R, Hwang S. Using evolutionary and structural information to predict DNA-binding residues in DNA-binding proteins. Proteins. 2006;64:19–27.

    CAS  Article  PubMed  Google Scholar 

  11. 11.

    Hwang S, Gou Z, Kuznetsov IB. DP-Bind: a web server for sequence-based prediction of DNA-binding residues in DNA-binding proteins. Bioinformatics. 2007;23:634–6.

    CAS  Article  PubMed  Google Scholar 

  12. 12.

    Murkami Y, Spriggs RV, Nakamura H, Jones S. PiRaNhA: a server for the computational prediction of RNA-binding residues in protein sequences. Nucleic Acids Res. 2010;38:W412–6.

    Article  CAS  Google Scholar 

  13. 13.

    Wang L, Huang C, Yang MQ, Yang JY. BindN+ for accurate prediction of DNA and RNA-binding residues from protein sequence features. BMC System Biology. 2010;4 Suppl 1:S3.

    Article  CAS  Google Scholar 

  14. 14.

    Kumar M, Gromiha MM, Raghava GPS. Prediction of RNA binding sites in a protein using SVM and PSSM profile. Proteins: Structure, Function and Bioinformatics. 2008;71:189–94.

    CAS  Article  Google Scholar 

  15. 15.

    Liu Z-P, Wu L-Y, Wang Y, Zhang X-S, Chen L. Prediction of protein-RNA binding sites by a random forest method with combined features. Bioinformatics. 2010;26:1616–22.

    CAS  Article  PubMed  Google Scholar 

  16. 16.

    Chen YC, Lim C. Predicting RNA-binding sites from the protein structure based on electrostatics, evolution and geometry. Nucleic Acids Res. 2008;36, e29.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  17. 17.

    Shazman S, Mandel-Gutfreund Y. Classifying RNA-binding proteins based on electrostatic properties. PLoS Comput Biol. 2008;4, e1000146.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  18. 18.

    Wang L, Brown SJ. BindN: a web-based tool for efficient prediction of DNA and RNA binding sites in amino acid sequences. Nucleic Acids Res. 2006;34:W243–8.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  19. 19.

    Ellis JJ, Broom M, Jones S. Protein-RNA interactions: structural analysis and functional classes. Proteins: Structure, Function, and Bioinformatics. 2007;66:903–11.

    CAS  Article  Google Scholar 

  20. 20.

    Jones S, Daley DT, Luscombe NM, Berman HM, Thornton JM. Protein-RNA interactions: a structural analysis. Nucleic Acids Res. 2001;29:943–54.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  21. 21.

    Lejeune D, Delsaux N, Charloteaux B, Thomas A, Brasseur R. Protein-nucleic acid recognition: statistical analysis of atomic interactions and influence of DNA structure. Proteins: Structure, Function, and Bioinformatics. 2005;61:258–71.

    CAS  Article  Google Scholar 

  22. 22.

    Berman HM, Westbrook J, Feng Z, Gilliland G, Bhat TN, Weissig H, et al. The Protein Data Bank. Nucleic Acids Res. 2000;28:235–42.

  23. 23.

    Puton T, Kozlowski L, Tuszynska I, Rother K, Bujnicki JM. Computational methods for prediction of protein–RNA interactions. Journal of Structural Biology. 2012;178:261–8.

    Article  CAS  Google Scholar 

  24. 24.

    Kim OT, Yura K, Go N. Amino acid residue doublet propensity in the protein–RNA interface and its application to RNA interface prediction. Nucleic Acids Res. 2006;34:6450–60.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  25. 25.

    Maetschke SR, Yuan Z. Exploiting structural and topological information to improve prediction of RNA-protein binding sites. BMC Bioinformatics. 2009;10:341.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  26. 26.

    Towfic F, Caragea C, Gemperline DC, Dobbs D, Honavar V. Struct-NB: predicting protein-RNA binding sites using structural features. Int J Data Min Bioinform. 2010;4:21–43.

    Article  PubMed  PubMed Central  Google Scholar 

  27. 27.

    Lewis M, Rees D. Fractal surfaces of proteins. Science. 1985;230:1163–5.

    CAS  Article  PubMed  Google Scholar 

  28. 28.

    Pintar A, Carugo O, Pongor S. Cx, an algorithm that identifies protruding atoms in proteins. Bioinformatics. 2002;18:980–4.

    CAS  Article  PubMed  Google Scholar 

  29. 29.

    Zhao H, Yang Y, Zhou Y. Structure-based prediction of RNA-binding domains and RNA-binding sites and application to structural genomics target. Nucleic Acids Res. 2010;39:3017–25.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  30. 30.

    Jeong E, Chung I, Miyano S. A neural network method for identification of RNA-interacting residues in protein. Genome Inform. 2004;15:105–16.

    CAS  PubMed  Google Scholar 

  31. 31.

    Cai YD, Lin SL. Support vector machines for predicting rRNA-, RNA-, and DNA-binding proteins from amino acid sequence. Biochim Biophis Acta-Proteins Proteomics. 2003;1648:127–33.

    CAS  Article  Google Scholar 

  32. 32.

    Han LY, Cai CZ, Lo SL, Chung MCM, Chen YZ. Prediction of RNA-binding proteins from primary sequence by a support vector machine approach. RNA. 2004;10:355–68.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  33. 33.

    Shao X, Tian Y, Wu L, Wang Y, Jing L, Deng N. Predicting DNA-, and RNA-binding proteins from sequences with kernel methods. J Thor Biol. 2009;258:289–93.

    CAS  Article  Google Scholar 

  34. 34.

    Terribilini M, Lee J, Yan C, Jernigan R, Honavar V, Dobbs D. Prediction of RNA-binding sites in proteins from amino acid sequence. RNA. 2006;16:1450–62.

    Article  CAS  Google Scholar 

  35. 35.

    Terribilini M, Sander J, Lee J-H, Zaback P, Jernigan RL, Honavar V, et al. RNABindR: a server for analysing and predicting RNA-binding sites in proteins. Nucleic Acids Res. 2007;35:W578–84.

  36. 36.

    Jeong E. A neural network method for identification of RNA-interacting residues in protein. Genome Inform. 2004;15:105–16.

    CAS  PubMed  Google Scholar 

  37. 37.

    Zhou H-X, Shan Y. Prediction of protein interaction sites from sequence profile and residue neighbour list. Proteins: Structures, Function, and Genetics. 2001;44:336–43.

    CAS  Article  Google Scholar 

  38. 38.

    Walia RR, Xue LC, Wilkins K, El-Manzalawy Y, Dobbs D, Honavar V. RNABindRPlus: A predictor that combines machine learning and sequence homology-based methods to improve the reliability of predicted RNA-binding residues in proteins. PLoS ONE. 2014;9, e97725.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  39. 39.

    Shazman S, Celniker G, Haber O, Glaser F, Mandel-Gutfreund Y. Patch Finder Plus (PFplus): a web server for extracting and displaying positive electrostatic patches on protein surfaces. Nucleic Acids Res. 2007;35:W526–30.

    Article  PubMed  PubMed Central  Google Scholar 

  40. 40.

    Lewis BA, Walia RR, Terribilini M, Ferguson J, Zheng C, Honavar V, et al. PRIDB: a protein-RNA interface database. Nucl Acids Res. 2011;39 Suppl 1:D277–82.

  41. 41.

    Rice P, Longden I, Bleasby A. EMBOSS: The European Molecular Biology Open Software Suite. Trends in Genetics. 2000;16:276–7.

    CAS  Article  PubMed  Google Scholar 

  42. 42.

    Humphrey W, Dalke A, Schulten K. VMD - Visual Molecular Dynamics. J Molec Graphics. 1996;14:33–8.

    CAS  Article  Google Scholar 

  43. 43.

    Arun KS, Huang TS, Blostein SD. Least-squares fitting of two 3-D point sets. IEEE Trans Pattern Anal Machine Intell. 1987;9:698–700.

    CAS  Article  Google Scholar 

  44. 44.

    Altschul S, Madden T, Schaffer A, Zhang J, Zhang Z, Miller W, et al. Gapped BLAST and PSI-BLAST: a new generation of protein database search programs. Nucleic Acids Res. 1997;25:3389–402.

  45. 45.

    Polikar R. Ensemble based systems in decision making. IEEE Circuits and Systems Magazine. 2006;6:21–45.

    Article  Google Scholar 

  46. 46.

    Rokach L. Ensemble-based classifiers. Artificial Intelligence Review. 2010;33:1–39.

    Article  Google Scholar 

  47. 47.

    ENTOOL []

Download references


This work was supported in part by the Natural Science Foundation of China under Grant 61303112 and 81373864, in part by the Key Laboratory of Intelligent Perception and Systems for High-Dimensional Information, Ministry of Education, Nanjing University of Science and Technology, Nanjing, China, under Grant 30920140122006.

Author information



Corresponding author

Correspondence to Ying Shen.

Additional information

Competing interests

The authors declare that they have no competing interests.

Authors’ contributions

YS conceived this project and designed the methodology, HR implemented the algorithm using Matlab and performed the statistical analysis. HR and YS drafted the manuscript. All authors have read and approved the final manuscript.

Rights and permissions

Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (, 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 ( applies to the data made available in this article, unless otherwise stated.

Reprints and Permissions

About this article

Verify currency and authenticity via CrossMark

Cite this article

Ren, H., Shen, Y. RNA-binding residues prediction using structural features. BMC Bioinformatics 16, 249 (2015).

Download citation


  • Protein-RNA interaction prediction
  • Structural information
  • Least-squares distance