Predicting RNA secondary structure using profile stochastic context-free grammars and phylogenic analysis
Stochastic context-free grammars (SCFGs) have been applied to predicting RNA secondary structure. The prediction of RNA secondary structure can be facilitated by incorporating with comparative sequence analysis. However, most of existing SCFG-based methods lack explicit phylogenic analysis of homologous RNA sequences, which is probably the reason why these methods are not ideal in practical application. Hence, we present a new SCFG-based method by integrating phylogenic analysis with the newly defined profile SCFG. The method can be summarized as: 1) we define a new profile SCFG, M, to depict consensus secondary structure of multiple RNA sequence alignment; 2) we introduce two distinct hidden Markov models, lambda and lambda', to perform phylogenic analysis of homologous RNA sequences. Here, lambda is for non-structural regions of the sequence and lambda' is for structural regions of the sequence, 3) we merge lambda and lambda' into M to devise a combined model for prediction of RNA secondary structure. We tested our method on data sets constructed from the Rfam database. The sensitivity and specificity of our method are more accurate than those of the predictions by Pfold.
rna secondary structure stochastic context-free grammar phylogenic analysis