Efficient parameter estimation for RNA secondary structure prediction
1Department of Computer Science, University of British Columbia, Vancouver BC V6T 1Z4, Canada and 2Department of Biochemistry & Biophysics and Department of Biostatistics & Computational Biology, University of Rochester Medical Center, Rochester NY 14642, USA
*To whom correspondence should be addressed.
 |
Abstract |
|---|
Motivation: Accurate prediction of RNA secondary structure from the base sequence is an unsolved computational challenge. The accuracy of predictions made by free energy minimization is limited by the quality of the energy parameters in the underlying free energy model. The most widely used model, the Turner99 model, has hundreds of parameters, and so a robust parameter estimation scheme should efficiently handle large data sets with thousands of structures. Moreover, the estimation scheme should also be trained using available experimental free energy data in addition to structural data.
Results: In this work, we present constraint generation (CG), the first computational approach to RNA free energy parameter estimation that can be efficiently trained on large sets of structural as well as thermodynamic data. Our CG approach employs a novel iterative scheme, whereby the energy values are first computed as the solution to a constrained optimization problem. Then the newly computed energy parameters are used to update the constraints on the optimization function, so as to better optimize the energy parameters in the next iteration. Using our method on biologically sound data, we obtain revised parameters for the Turner99 energy model. We show that by using our new parameters, we obtain significant improvements in prediction accuracy over current state of-the-art methods.
Availability: Our CG implementation is available at http://www.rnasoft.ca/CG/
Contact: andrones{at}cs.ubc.ca
CiteULike 
Connotea 
Del.icio.us What's this?
This article has been cited by other articles:
|
|
 |
|
  D. H. Turner and D. H. Mathews NNDB: the nearest neighbor parameter database for predicting stability of nucleic acid secondary structure Nucleic Acids Res., January 1, 2010; 38(suppl_1): D280 - D282. [Abstract] [Full Text] [PDF]
|
|
|
|
 |
|
 J. Zhang, J. Dundas, M. Lin, R. Chen, W. Wang, and J. Liang Prediction of geometrically feasible three-dimensional structures of pseudoknotted RNA through free energy estimation RNA, December 1, 2009; 15(12): 2248 - 2263. [Abstract] [Full Text] [PDF]
|
|
|
|
|
|
Z. J. Lu, J. W. Gloor, and D. H. Mathews Improved RNA secondary structure prediction by maximizing expected pair accuracy RNA, October 1, 2009; 15(10): 1805 - 1813. [Abstract] [Full Text] [PDF]
|
|
|
|
 |
|
 M. Hamada, K. Sato, H. Kiryu, T. Mituyama, and K. Asai Predictions of RNA secondary structure by combining homologous sequence information Bioinformatics, June 15, 2009; 25(12): i330 - i338. [Abstract] [Full Text] [PDF]
|
|
|
|
|
|
S. Smit, K. Rother, J. Heringa, and R. Knight From knotted to nested RNA structures: A variety of computational methods for pseudoknot removal RNA, March 1, 2008; 14(3): 410 - 416. [Abstract] [Full Text] [PDF]
|
|