BACKGROUND
rna-binding proteins  are designed to efficiently recognize specific rna sequences after they are derived from the dna sequences.
protein-rna interactions are fundamental to cellular processes, including the assembly and function of ribonucleoprotein particles , such as ribosomes and spliceosomes and the post-transcriptional regulation of gene products.
for satisfying diverse functional requirements, rna binding proteins are composed of multiple blocks of rna-binding domains  presented in various structural arrangements to provide versatile functionality  <cit> .
although rna structure is hierarchical, that is, the primary sequence determines the secondary structure which, in turns, determines tertiary structure, the tertiary structure of rna is not as stable as secondary structure and is hard to predict  <cit> .
however, sequence conservations in rna-binding domains have been discovered in rna-binding proteins  <cit> .
with the recent growth of protein-rna complexes in the protein data bank   <cit>  and the nucleic acid database   <cit> , structural analysis on rna-binding pockets  <cit>  and the themes of rna-protein recognition  <cit>  have been investigated as well.
most recent works on predicting rna-binding residues used support vector machine  with protein evolutionary information from protein sequence.
wang and brown  developed the web service, bindn  <cit> , to predict dna and rna binding sites using sequence features to represent structural characteristics including relative solvent accessible surface area, side chain pka, hydrophobicity index and molecular mass of an amino acid.
tong et al.
<cit>  proposed the hybrid risp  method by adjusting cutoff value of svm discrimination function to improve prediction performance.
kumar et al.
developed pprint  <cit>  by using evolutionary profiles of the position-specific scoring matrices  and amino acid composition while they also adjusted cutoff value of svm discrimination function to improve prediction performance.
wang et al.
developed printr  <cit>  by using additional structural information from protein-rna complexes.
cheng et al.
developed rnaprob  <cit>  by smoothing pssm profiles with consideration of the correlation and dependency from the neighboring residues for each amino acid in a protein.
spriggs et al.
<cit>  developed the piranha by using support vector machine with a pssm profile and three amino acid properties, including interface propensity , predicted solvent accessibility  and hydrophobicity  for recognizing rna-binding residues  <cit> .
other machine learning approaches such as neural network and naïve bayes classifier have also been applied to predict rna-binding residues.
jeong et al.
<cit>  applied artificial neural network -based method with amino acid sequence and predicted secondary structure information and improved the performance by using post-processing procedures such as state-shifting and filtering isolated interacting residues from prediction.
improved version by jeong et al.
<cit>  used evolutionary information extracted from psi-blast profiles and clustalw alignment.
terribilini et al.
<cit>  applied a naïve bayes classifier with amino acid sequence information for predicting rna interacting residues and presented the results through the web service rnabindr  <cit> .
the ability to computationally predict rna-binding residues in a rna-binding protein can help biologists reveal site-directed mutagenesis in wet-lab experiments.
caragea et al.
<cit>  explored the problem of assessing the performance of classifiers trained on macromolecular sequence data, with the emphasis on cross-validation and data selection methods.
in comparison of window-based k-fold cross-validation and sequence-based k-fold cross-validation, window-based cross-validation can yield overly optimistic estimates of the performance of classifier relative to the estimates obtained using sequence-based cross-validation.
rnaprob, bindn, risp, printr and piranha are predictors that report performance window-based k-fold cross-validation while pprint and rnabindr report performance with sequence-based k-fold cross-validation.
the predictors evaluated with window-based k-fold cross-validation have superior performance than those with sequence-based k-fold cross-validation.
the reason is that data instances in the testing fold would be predicted by data instances with sub-sequence identity higher than 25% in the training fold in window-based k-fold cross-validation.
therefore, in data with class imbalance, the metrics that measure the classification performance must be chosen carefully.
matthew's correlation coefficient , f-score and f <dig> -score are widely applied to assess the prediction performance.
mcc is used to measure prediction quality with the consideration of both under- and over-predictions.
f-score and f <dig> -score are used to assess balanced prediction quality on both positive class and negative class.
in this article, we proposed the prediction framework “proterna” with the combination of svm-based classifier with evolutionary profiles and conserved residues discovery by sequence conservation for identifying rna-interacting residues in a rna-binding protein.
in the svm-based classifier, we use features including position-specific scoring matrix computed by psi-blast and secondary structure information predicted by psipred as feature vectors  <cit> .
to exploit the sequence conservation information, wildspan  <cit>  , which is developed to discover functional signatures and diagnostic patterns of proteins directly from a set of unaligned protein sequences, is incorporated.
the most distinguishing feature of wildspan is that it links short motifs  with large flexible gaps to deliver the most frequently observed discontinuous patterns present in related proteins.
wildspan has been embedded in many applications  <cit>  to discover functionally important residues; therefore, we apply wildspan to discover conserved residues as rna-binding residues in a protein sequence to improve prediction performance on detecting more rna-binding residues.
the independent testing dataset collected for performance evaluation contains  <dig> testing rna-binding proteins with less than 30% sequence identity against with training data.
in the independent testing dataset, proterna has been able to deliver overall accuracy of  <dig> %, mcc of  <dig> , f-score of  <dig> , and f <dig> -score of  <dig> .
we emphasize mcc, f-score and f <dig> -score because it provides the biochemist with a confidence level for designing an experiment to confirm whether a predicted binding residue is really involved in interaction with the rna.
RESULTS
in this section, we will report the experiments conducted to evaluate the performance of our proposed approach, proterna with the combination of svm-based classifier with evolutionary profiles and conserved residues discovery by sequence conservation.
in order to avoid bias, we repeated 5-fold cross-validation procedure  <dig> times to observe prediction performance on the training dataset rb <dig> .
for each run, we applied sequence-based 5-fold cross-validation; therefore, protein chains will be randomly divided into  <dig> folds: one fold for testing and remaining  <dig> folds for training.
for this study, libsvm  was used for data training and classification and wildspan was used for detecting conserved residues from homologous protein sequences.
we use independent testing dataset containing  <dig> protein chains for comparing proterna with other predictors such as piranha, pprint, bindn, and prip.
performance evaluation by five-fold cross-validation
in order to avoid bias, we repeated 5-fold cross-validation  <dig> times to observe prediction performance and experimental result was shown in table  <dig>  only using svm-based classifier, proternasvm delivers overall sensitivity of  <dig> %, specificity of  <dig> %, precision of  <dig> %, accuracy of  <dig> %, mcc of  <dig> , f-score of  <dig>  and f <dig> -score of  <dig> .
since the experiments are repeated  <dig> times for reducing prediction bias, standard deviation for each assessment is also listed.
the results have been obtained using the training parameters, c =  <dig>  γ = 2- <dig>  which give better results than other values for prediction of rna-binding residues.
then wildspan only outputs patterns for each input protein chain once so there is no information about standard deviation for each assessment.
proternawildspan delivers overall sensitivity of  <dig> %, specificity of  <dig> %, precision of  <dig> %, accuracy of  <dig> %, mcc of  <dig> , f-score of  <dig>  and f <dig> -score of  <dig> .
after combining prediction results by proternasvm and proternawildspan, proterna delivers overall sensitivity of  <dig> %, specificity of  <dig> %, precision of  <dig> %, accuracy of  <dig> %, mcc of  <dig> , f-score of  <dig>  and f <dig> -score of  <dig> .
as reported by towfic et al.
<cit> , over half  of the rnas are rrnas in the dataset of rb <dig>  according to their study, table  <dig> shows the distribution of different categories of rnas on rna-binding residues.
rrna is the major group that contains about 38% positive samples in rrna group.
remaining groups presents highly imbalanced class dataset, containing about 10% positive sample in average.
if the predictor tries to predict all samples as negative class exclusive of rrna group, the predictor may gain better performance in assessment but provide no clues for biologists.
table  <dig> describes the average prediction performance of  <dig> runs of 5-fold cross-validation; however, we only choose one of the repeated experiments that had a performance that is close to the average performance for detailed analysis in table  <dig>  as shown in table  <dig> proternawildspan predicts an equal amount of true positives and false positives on average.
previous research on studying rna-binding domains revealed that rna binding proteins are composed of multiple blocks of rna-binding domains to provide versatile functionality.
therefore, conserved residues in the same rna-binding domain from different rna-binding proteins would not always interact with a specific rna.
furthermore, while combining prediction results predicted by proternasvm and proternawildspan, proternawildspan detected additional rna-binding residues that proternasvm didn’t predict.
as we known, rrna is the major group among the training dataset.
comparing the amount of rna-binding proteins in terms of interacting target , we find that trna generally has the most interaction partners followed by mrna and rrna has the least partners.
proternasvm tends to predict negative for proteins in the mrna group and over-predict either positive class or negative class in trna group.
however, proternawildspan shows no different between categories of rnas because of discovered homologous proteins in swiss-prot.
in addition, proternawildspan detects conserved residues as binding residues that cover regions that proternasvm doesn’t predict; therefore, we apply wildspan to detect conserved residues because these conserved residues have higher probability to play roles in interacting rnas.
comparison with other predictors by independent testing
only predictors that predict rna-binding residues from protein primary sequence information were selected for performance comparison.
in addition, risp did not respond with any prediction results after submitting the jobs and printr is unavailable.
according to the designed framework of rnabindr, if there is an exact matched protein chain in protein data bank, rnabindr will return the actual rna-binding residues of protein-rna complex.
therefore, it is difficult to distinguish whether the returned result is actually binding or predicted binding so rnabindr will be excluded.
finally, table  <dig> shows the prediction performance of proterna in comparison with piranha, pprint, bindn, and prip.
while ordering prediction performance in terms of mcc, proterna delivers better performance than other predictors in accuracy, mcc, and f <dig> -score.
*order by mcc.
2qbe_t
2der_b
2g4b_a
2qbe_t
2g4b_a
2g4b_a
2der_b
2der_b
2g4b_a
2der_b
2qbe_t
2qbe_t
2qbe_t
2qbe_t
2qbe_t
2qbe_t
2der_b
2i82_c
2der_b
2g4b_a
2v47_c
2der_b
2v47_c
2i82_c
2v47_c
2i82_c
2i82_c
2qbe_t
2g4b_a
2g4b_a
2der_b
2g4b_a
2v47_c
values in bold indicate listing in the top  <dig> by at least  <dig> predictors.
values in bold and italics indicate listing in the top  <dig> by at least  <dig> predictors.
CONCLUSIONS
this article presents the design of a sequence based predictor aiming to identify the rna-binding residues in a rna-binding protein by machine learning and pattern mining approaches.
rna-binding proteins play different roles while interacting with different categories of rnas to represent diverse functions.
however, rna-binding proteins are accommodated by the presence of multiple copies of these rna-binding domains presented in various structural arrangements to expand the functional repertoire of rna-binding proteins.
therefore, it is still difficult to predict rna-binding residues in a rna-binding protein.
furthermore, predicting rna-binding residues in a rna-binding protein can help biologists reveal site-directed mutagenesis in wet-lab experiments.
in the experiments reported in this article, proterna used not only evolutionary profile with predicted secondary structure but also sequence conservation information.
although these conserved residues can be functional conserved residues or structural conserved residues, they also provide clues to indicate the important residues in a protein sequence.
in the independent testing dataset, proterna has been able to deliver overall accuracy of  <dig> %, mcc of  <dig> , f-score of  <dig> , and f <dig> -score of  <dig> .
it is anticipated that the prediction accuracy delivered by proterna will continue to improve as the number of protein-rna complexes deposited in the pdb continues to grow and the number of training samples that can be exploited continues to increase accordingly.
nevertheless, it is the computational biologists’ primary interest to develop more advanced prediction mechanisms.
in this respect, we believe that, as the number of protein-rna complexes deposited in the pdb increases, we can obtain more insights about the key physiochemical properties that play essential roles in protein-rna interactions and then we will be able to develop more advanced prediction mechanisms accordingly.
in addition, we will exploit the experiences learned in this study in order to design specific predictors for other families of proteins interacting with rna.
we believe that different families of proteins may have very different characteristics.
therefore, concerning a specific type of proteins, a specifically-designed predictor should be able to deliver superior performance in compared to a general-purpose predictor.
materials and methods
datasets
we used rb <dig> as the training dataset for predicting rna-binding residues in a protein collected by terribilini et al., containing  <dig> non-redundant protein chains with resolution better than  <dig>  Å in the pdb solved by x-ray crystallography  <cit> .
no two protein chains has a sequence identity greater than 30%.
based on the cut-off distance of  <dig> Å, a total of  <dig>  amino acids are in rb <dig>  which contains  <dig>  rna-binding residues and  <dig>  non-binding residues.
the list of pdb ids of the training dataset, rb <dig>  is shown in table  <dig>
in order to evaluate prediction performance among different prediction models, we collected a new independent testing dataset by extracting all structures of protein-rna complexes from the pdb that were added after january  <dig>  protein chains with a resolution better than  <dig>  Å and sequence length of protein chain longer than  <dig> amino acids will be reserved.
we then performed a redundancy reduction using blastclust  <cit>  to ensure that none of the chains showed a sequence similarity of more than 30% within the dataset and also in the training dataset; therefore,  <dig> protein-rna complexes were selected to create a dataset called rb <dig>  the list of pdb ids in rb <dig> are shown in table  <dig>  based on the cut-off distance of  <dig> Å, a total of  <dig>  amino acids are in rb <dig>  which contains  <dig> rna-binding residues and  <dig>  non-binding residues.
framework for prediction rna-interacting residues
where x is the entry value in 20xn matrix of pssm .
with the consideration of structure information, we also used secondary structure information predicted by psipred; the predicted secondary structure information consists of three probability values that represent helix, sheet and coil respectively  = ).
in addition, each residue was labelled based on whether it is involved in binding with the rna or not.
therefore, for each residue in a protein sequence, we construct a  <dig> *  <dig> =  <dig> dimensional feature factor ; the  <dig> dimensions include  <dig> features from pssm,  <dig> features from psipred and a boundary flag.
as shown in figure  <dig>  the detail data flow and feature vector preparation for svm-based classifier is addressed.
the best parameters selected for predicting rna-binding residues is decided by 5-fold cross-validation.
in the part of wildspan , for protein-based mining suggested by the authors, at most  <dig> unique homologous proteins with sequence identity ranged from 30% to 90% are required by searching against swiss-prot sequence database with psi-blast .
then we applied default parameter to obtain patterns by wildspan.
wildspan can’t generate any pattern if there are not enough homologous proteins selected from swiss-prot protein sequence database or too similar homologous proteins.
significance and performance evaluation
the predictions made for the testing instances are compared with the defined class labels  to evaluate the predictor.
the accuracy is defined as
where tp is the number of true positives ; tn is the number of true negatives ; fp is the number of false positives  and fn is the number of false negatives .
matthew's correlation coefficient  is defined as follows:
mcc is used to measure prediction performance with the consideration of both under- and over-predictions, where mcc =  <dig> denotes a perfect prediction, mcc =  <dig> indicates a completely random assignment, and mcc = - <dig> means a perfectly reverse correlation.
since the data for rna-binding residue prediction is skewed, so-called class imbalanced data, the accuracy alone may be misleading.
the predictor can achieve 85% accuracy by simply predicting all residues as negative for datasets where the positive to negative sample ratio is 1: <dig>  therefore, prediction performance on positive class and negative class should be assessed individually.
metrics of the specificity and sensitivity can help predictors to know their prediction performance on positive and negative samples respectively.
the sensitivity is used to measure the prediction capability of positive samples; the specificity is used to measure the prediction capability of negative samples.
specificity and sensitivity are defined as follows:
in addition, precision and fβ-score are also defined as follows:
precision is used to assess prediction power on positive class.
f-score  is the harmonic mean of precision and sensitivity if β =  <dig>  f <dig> -score weights precision twice as much as sensitivity if β =  <dig> .
list of abbreviations
rbp: rna-binding protein, rbd: rna-binding domain; rnp: ribonucleoprotein particle; pssm: position-specific scoring matrix; ndb: nucleic acid database; pdb: protein data bank; svm: support vector machine.
competing interests
the authors declare that they have no competing interests.
authors' contributions
yfh and lyc developed and implemented the overall framework and drafted the manuscript; cch provided valuable suggestion on experiments based on his previous works on predicting dna-binding residues.
ckh conceived of the study and participated in its design and coordination and helped to draft this manuscript.
all authors have read and approved the final manuscript.

 rna-binding proteins are designed to efficiently recognize specific rna sequences after they are derived from the dna sequences. protein-rna interactions are fundamental to cellular processes, including the assembly and function of ribonucleoprotein particles, such as ribosomes and spliceosomes and the post-transcriptional regulation of gene products. for satisfying diverse functional requirements, rna binding proteins are composed of multiple blocks of rna-binding domains presented in various structural arrangements to provide versatile functionality. although rna structure is hierarchical, that is, the primary sequence determines the secondary structure which, in turns, determines tertiary structure, the tertiary structure of rna is not as stable as secondary structure and is hard to predict. however, sequence conservations in rna-binding domains have been discovered in rna-binding proteins. with the recent growth of protein-rna complexes in the protein data bank and the nucleic acid database, structural analysis on rna-binding pockets and the themes of rna-protein recognition have been investigated as well. most recent works on predicting rna-binding residues used support vector machine with protein evolutionary information from protein sequence. wang and brown developed the web service, bindn, to predict dna and rna binding sites using sequence features to represent structural characteristics including relative solvent accessible surface area, side chain pka, hydrophobicity index and molecular mass of an amino acid. tong et al. proposed the hybrid risp method by adjusting cutoff value of svm discrimination function to improve prediction performance. kumar et al. developed pprint by using evolutionary profiles of the position-specific scoring matrices and amino acid composition while they also adjusted cutoff value of svm discrimination function to improve prediction performance. wang et al. developed printr by using additional structural information from protein-rna complexes. cheng et al. developed rnaprob by smoothing pssm profiles with consideration of the correlation and dependency from the neighboring residues for each amino acid in a protein. spriggs et al. developed the piranha by using support vector machine with a pssm profile and three amino acid properties, including interface propensity, predicted solvent accessibility and hydrophobicity for recognizing rna-binding residues. other machine learning approaches such as neural network and naïve bayes classifier have also been applied to predict rna-binding residues. jeong et al. applied artificial neural network -based method with amino acid sequence and predicted secondary structure information and improved the performance by using post-processing procedures such as state-shifting and filtering isolated interacting residues from prediction. improved version by jeong et al. used evolutionary information extracted from psi-blast profiles and clustalw alignment. terribilini et al. applied a naïve bayes classifier with amino acid sequence information for predicting rna interacting residues and presented the results through the web service rnabindr. the ability to computationally predict rna-binding residues in a rna-binding protein can help biologists reveal site-directed mutagenesis in wet-lab experiments. caragea et al. explored the problem of assessing the performance of classifiers trained on macromolecular sequence data, with the emphasis on cross-validation and data selection methods. in comparison of window-based k-fold cross-validation and sequence-based k-fold cross-validation, window-based cross-validation can yield overly optimistic estimates of the performance of classifier relative to the estimates obtained using sequence-based cross-validation. rnaprob, bindn, risp, printr and piranha are predictors that report performance window-based k-fold cross-validation while pprint and rnabindr report performance with sequence-based k-fold cross-validation. the predictors evaluated with window-based k-fold cross-validation have superior performance than those with sequence-based k-fold cross-validation. the reason is that data instances in the testing fold would be predicted by data instances with sub-sequence identity higher than 25% in the training fold in window-based k-fold cross-validation. therefore, in data with class imbalance, the metrics that measure the classification performance must be chosen carefully. matthew's correlation coefficient, f-score and f -score are widely applied to assess the prediction performance. mcc is used to measure prediction quality with the consideration of both under- and over-predictions. f-score and f -score are used to assess balanced prediction quality on both positive class and negative class. in this article, we proposed the prediction framework “proterna” with the combination of svm-based classifier with evolutionary profiles and conserved residues discovery by sequence conservation for identifying rna-interacting residues in a rna-binding protein. in the svm-based classifier, we use features including position-specific scoring matrix computed by psi-blast and secondary structure information predicted by psipred as feature vectors. to exploit the sequence conservation information, wildspan, which is developed to discover functional signatures and diagnostic patterns of proteins directly from a set of unaligned protein sequences, is incorporated. the most distinguishing feature of wildspan is that it links short motifs with large flexible gaps to deliver the most frequently observed discontinuous patterns present in related proteins. wildspan has been embedded in many applications to discover functionally important residues; therefore, we apply wildspan to discover conserved residues as rna-binding residues in a protein sequence to improve prediction performance on detecting more rna-binding residues. the independent testing dataset collected for performance evaluation contains testing rna-binding proteins with less than 30% sequence identity against with training data. in the independent testing dataset, proterna has been able to deliver overall accuracy of %, mcc of, f-score of, and f -score of. we emphasize mcc, f-score and f -score because it provides the biochemist with a confidence level for designing an experiment to confirm whether a predicted binding residue is really involved in interaction with the rna.  in this section, we will report the experiments conducted to evaluate the performance of our proposed approach, proterna with the combination of svm-based classifier with evolutionary profiles and conserved residues discovery by sequence conservation. in order to avoid bias, we repeated 5-fold cross-validation procedure times to observe prediction performance on the training dataset rb. for each run, we applied sequence-based 5-fold cross-validation; therefore, protein chains will be randomly divided into folds: one fold for testing and remaining folds for training. for this study, libsvm was used for data training and classification and wildspan was used for detecting conserved residues from homologous protein sequences. we use independent testing dataset containing protein chains for comparing proterna with other predictors such as piranha, pprint, bindn, and prip. performance evaluation by five-fold cross-validation in order to avoid bias, we repeated 5-fold cross-validation times to observe prediction performance and experimental result was shown in table only using svm-based classifier, proternasvm delivers overall sensitivity of %, specificity of %, precision of %, accuracy of %, mcc of, f-score of and f -score of. since the experiments are repeated times for reducing prediction bias, standard deviation for each assessment is also listed. the results have been obtained using the training parameters, c = γ = 2- which give better results than other values for prediction of rna-binding residues. then wildspan only outputs patterns for each input protein chain once so there is no information about standard deviation for each assessment. proternawildspan delivers overall sensitivity of %, specificity of %, precision of %, accuracy of %, mcc of, f-score of and f -score of. after combining prediction results by proternasvm and proternawildspan, proterna delivers overall sensitivity of %, specificity of %, precision of %, accuracy of %, mcc of, f-score of and f -score of. as reported by towfic et al., over half of the rnas are rrnas in the dataset of rb according to their study, table shows the distribution of different categories of rnas on rna-binding residues. rrna is the major group that contains about 38% positive samples in rrna group. remaining groups presents highly imbalanced class dataset, containing about 10% positive sample in average. if the predictor tries to predict all samples as negative class exclusive of rrna group, the predictor may gain better performance in assessment but provide no clues for biologists. table describes the average prediction performance of runs of 5-fold cross-validation; however, we only choose one of the repeated experiments that had a performance that is close to the average performance for detailed analysis in table as shown in table proternawildspan predicts an equal amount of true positives and false positives on average. previous research on studying rna-binding domains revealed that rna binding proteins are composed of multiple blocks of rna-binding domains to provide versatile functionality. therefore, conserved residues in the same rna-binding domain from different rna-binding proteins would not always interact with a specific rna. furthermore, while combining prediction results predicted by proternasvm and proternawildspan, proternawildspan detected additional rna-binding residues that proternasvm didn’t predict. as we known, rrna is the major group among the training dataset. comparing the amount of rna-binding proteins in terms of interacting target, we find that trna generally has the most interaction partners followed by mrna and rrna has the least partners. proternasvm tends to predict negative for proteins in the mrna group and over-predict either positive class or negative class in trna group. however, proternawildspan shows no different between categories of rnas because of discovered homologous proteins in swiss-prot. in addition, proternawildspan detects conserved residues as binding residues that cover regions that proternasvm doesn’t predict; therefore, we apply wildspan to detect conserved residues because these conserved residues have higher probability to play roles in interacting rnas. comparison with other predictors by independent testing only predictors that predict rna-binding residues from protein primary sequence information were selected for performance comparison. in addition, risp did not respond with any prediction results after submitting the jobs and printr is unavailable. according to the designed framework of rnabindr, if there is an exact matched protein chain in protein data bank, rnabindr will return the actual rna-binding residues of protein-rna complex. therefore, it is difficult to distinguish whether the returned result is actually binding or predicted binding so rnabindr will be excluded. finally, table shows the prediction performance of proterna in comparison with piranha, pprint, bindn, and prip. while ordering prediction performance in terms of mcc, proterna delivers better performance than other predictors in accuracy, mcc, and f -score. *order by mcc. 2qbe_t 2der_b 2g4b_a 2qbe_t 2g4b_a 2g4b_a 2der_b 2der_b 2g4b_a 2der_b 2qbe_t 2qbe_t 2qbe_t 2qbe_t 2qbe_t 2qbe_t 2der_b 2i82_c 2der_b 2g4b_a 2v47_c 2der_b 2v47_c 2i82_c 2v47_c 2i82_c 2i82_c 2qbe_t 2g4b_a 2g4b_a 2der_b 2g4b_a 2v47_c values in bold indicate listing in the top by at least predictors. values in bold and italics indicate listing in the top by at least predictors.  this article presents the design of a sequence based predictor aiming to identify the rna-binding residues in a rna-binding protein by machine learning and pattern mining approaches. rna-binding proteins play different roles while interacting with different categories of rnas to represent diverse functions. however, rna-binding proteins are accommodated by the presence of multiple copies of these rna-binding domains presented in various structural arrangements to expand the functional repertoire of rna-binding proteins. therefore, it is still difficult to predict rna-binding residues in a rna-binding protein. furthermore, predicting rna-binding residues in a rna-binding protein can help biologists reveal site-directed mutagenesis in wet-lab experiments. in the experiments reported in this article, proterna used not only evolutionary profile with predicted secondary structure but also sequence conservation information. although these conserved residues can be functional conserved residues or structural conserved residues, they also provide clues to indicate the important residues in a protein sequence. in the independent testing dataset, proterna has been able to deliver overall accuracy of %, mcc of, f-score of, and f -score of. it is anticipated that the prediction accuracy delivered by proterna will continue to improve as the number of protein-rna complexes deposited in the pdb continues to grow and the number of training samples that can be exploited continues to increase accordingly. nevertheless, it is the computational biologists’ primary interest to develop more advanced prediction mechanisms. in this respect, we believe that, as the number of protein-rna complexes deposited in the pdb increases, we can obtain more insights about the key physiochemical properties that play essential roles in protein-rna interactions and then we will be able to develop more advanced prediction mechanisms accordingly. in addition, we will exploit the experiences learned in this study in order to design specific predictors for other families of proteins interacting with rna. we believe that different families of proteins may have very different characteristics. therefore, concerning a specific type of proteins, a specifically-designed predictor should be able to deliver superior performance in compared to a general-purpose predictor. materials and methods datasets we used rb as the training dataset for predicting rna-binding residues in a protein collected by terribilini et al., containing non-redundant protein chains with resolution better than Å in the pdb solved by x-ray crystallography. no two protein chains has a sequence identity greater than 30%. based on the cut-off distance of Å, a total of amino acids are in rb which contains rna-binding residues and non-binding residues. the list of pdb ids of the training dataset, rb is shown in table in order to evaluate prediction performance among different prediction models, we collected a new independent testing dataset by extracting all structures of protein-rna complexes from the pdb that were added after january protein chains with a resolution better than Å and sequence length of protein chain longer than amino acids will be reserved. we then performed a redundancy reduction using blastclust to ensure that none of the chains showed a sequence similarity of more than 30% within the dataset and also in the training dataset; therefore, protein-rna complexes were selected to create a dataset called rb the list of pdb ids in rb are shown in table based on the cut-off distance of Å, a total of amino acids are in rb which contains rna-binding residues and non-binding residues. framework for prediction rna-interacting residues where x is the entry value in 20xn matrix of pssm. with the consideration of structure information, we also used secondary structure information predicted by psipred; the predicted secondary structure information consists of three probability values that represent helix, sheet and coil respectively = ). in addition, each residue was labelled based on whether it is involved in binding with the rna or not. therefore, for each residue in a protein sequence, we construct a * = dimensional feature factor; the dimensions include features from pssm, features from psipred and a boundary flag. as shown in figure the detail data flow and feature vector preparation for svm-based classifier is addressed. the best parameters selected for predicting rna-binding residues is decided by 5-fold cross-validation. in the part of wildspan, for protein-based mining suggested by the authors, at most unique homologous proteins with sequence identity ranged from 30% to 90% are required by searching against swiss-prot sequence database with psi-blast. then we applied default parameter to obtain patterns by wildspan. wildspan can’t generate any pattern if there are not enough homologous proteins selected from swiss-prot protein sequence database or too similar homologous proteins. significance and performance evaluation the predictions made for the testing instances are compared with the defined class labels to evaluate the predictor. the accuracy is defined as where tp is the number of true positives; tn is the number of true negatives; fp is the number of false positives and fn is the number of false negatives. matthew's correlation coefficient is defined as follows: mcc is used to measure prediction performance with the consideration of both under- and over-predictions, where mcc = denotes a perfect prediction, mcc = indicates a completely random assignment, and mcc = - means a perfectly reverse correlation. since the data for rna-binding residue prediction is skewed, so-called class imbalanced data, the accuracy alone may be misleading. the predictor can achieve 85% accuracy by simply predicting all residues as negative for datasets where the positive to negative sample ratio is 1: therefore, prediction performance on positive class and negative class should be assessed individually. metrics of the specificity and sensitivity can help predictors to know their prediction performance on positive and negative samples respectively. the sensitivity is used to measure the prediction capability of positive samples; the specificity is used to measure the prediction capability of negative samples. specificity and sensitivity are defined as follows: in addition, precision and fβ-score are also defined as follows: precision is used to assess prediction power on positive class. f-score is the harmonic mean of precision and sensitivity if β = f -score weights precision twice as much as sensitivity if β =. list of abbreviations rbp: rna-binding protein, rbd: rna-binding domain; rnp: ribonucleoprotein particle; pssm: position-specific scoring matrix; ndb: nucleic acid database; pdb: protein data bank; svm: support vector machine. competing interests the authors declare that they have no competing interests. authors' contributions yfh and lyc developed and implemented the overall framework and drafted the manuscript; cch provided valuable suggestion on experiments based on his previous works on predicting dna-binding residues. ckh conceived of the study and participated in its design and coordination and helped to draft this manuscript. all authors have read and approved the final manuscript. 
to satisfy diverse functional requirements, rna binding proteins are composed of multiple blocks of rna-binding domains presented in various structural arrangements to provide versatile functions. rna-binding proteins have diverse functions while interacting with different categories of rnas because these proteins are composed of multiple copies of rna-binding domains presented in various structural arrangements to expand the functional repertoire of rna-binding proteins. rbps are designed to efficiently recognize specific rna sequences after it is derived from the dna sequence. in the independent testing dataset, proterna has been able to deliver overall accuracy of %, mcc of, f-score of, and f -score of. the ability to computationally predict rna-binding residues in a rna-binding protein can help biologists reveal important site-directed mutagenesis in wet-lab experiments. furthermore, predicting rna-binding residues in a rna-binding protein can help biologists reveal important site-directed mutagenesis in wet-lab experiments. although these conserved residues can be either functionally conserved residues or structurally conserved residues, they provide clues on the important residues in a protein sequence. this article presents the design of a sequence-based predictor aiming to identify the rna-binding residues in a rna-binding protein by combining machine learning and pattern mining approaches.
 rna-binding proteins play crucial roles in post-transcriptional control of rna. rbps are designed to efficiently recognize specific rna sequences after it is derived from the dna sequence. to satisfy diverse functional requirements, rna binding proteins are composed of multiple blocks of rna-binding domains presented in various structural arrangements to provide versatile functions. the ability to computationally predict rna-binding residues in a rna-binding protein can help biologists reveal important site-directed mutagenesis in wet-lab experiments.  the proposed prediction framework named “proterna” combines a svm-based classifier with conserved residue discovery by wildspan to identify the residues that interact with rna in a rna-binding protein. although these conserved residues can be either functionally conserved residues or structurally conserved residues, they provide clues on the important residues in a protein sequence. in the independent testing dataset, proterna has been able to deliver overall accuracy of %, mcc of, f-score of, and f -score of.  this article presents the design of a sequence-based predictor aiming to identify the rna-binding residues in a rna-binding protein by combining machine learning and pattern mining approaches. rna-binding proteins have diverse functions while interacting with different categories of rnas because these proteins are composed of multiple copies of rna-binding domains presented in various structural arrangements to expand the functional repertoire of rna-binding proteins. furthermore, predicting rna-binding residues in a rna-binding protein can help biologists reveal important site-directed mutagenesis in wet-lab experiments. 26– september asia pacific bioinformatics network ninth international conference on bioinformatics tokyo, japan 