A New Free Energy Rule Toward Psuedo-knotted RNA Structure Prediction
Jian Zhang,Joe Dundas,Ming Lin,Rong Chen,Wei Wang,Jie Liang
DOI: https://doi.org/10.1016/j.bpj.2008.12.3025
IF: 3.4
2009-01-01
Biophysical Journal
Abstract:Accurate free energy estimation is essential for the prediction of RNA secondary structure. The widely Turner model works well for nested structures. For pseudoknotted RNAs, there is no effective free energy rules. We present a new free energy estimation method that goes beyond the Turner model based on the physical consideration that the free energy of loops (including pseudoknotted loops) are mostly determined by the loop length and the end-end distance. Given a specific secondary structure (refer to as topology hereafter), we grow the RNA in three dimensional space and search among all possible arrangement of helices for stable structure. The free energy contribution of helices is estimated from the Turner model, and the free energy of loops is calculated based on our new approach. We assess our free energy model by testing its ability in selecting native structure from decoy topologies. For a testing set of 50-pseudoknotted RNAs, we create thousands of competitive decoy topologies and calculate their free energy. On average, our energy model provides an average of 7-8 percent improvement in term of sensitivity and specificity. In a blind pseudoknot structure prediction without inserting the native topology into the candidate list, our prediction has an average of 0.79 sensitives and 0.75 specificity. This performance is comparable to that of TdFOLD (0.76, 0.79) and better than Hotknots (0.69, 0.72), ILM(0.73, 0.69) and PKNOTS (0.78,0.73). Our free energy estimation method is fundamentally different from the others: it is based on the growing of conformations in three dimensional space and is well suited for predicting pseudoknotted structures. In addition, our method also gives coarse 3D conformation of a given secondary structure, which can lead to accurate 3D ab initio predicted structure upon further refinement.