"Well-determined" regions in RNA secondary structure prediction: analysis of small subunit ribosomal RNA.

Recent structural analyses of genomic RNAs from RNA coliphages suggest that both well-determined base paired helices and well-determined structural domains that are identified by "energy dot plot" analysis using the RNA folding package mfold, are likely to be predicted correctly. To test these observations with another group of large RNAs, we have analyzed 15 ribosomal RNAs. Published secondary structure models that were derived by comparative sequence analysis were used to evaluate the predicted structures. Both the optimal predicted fold and the predicted "energy dot plot" of each sequence were examined. Each prediction was obtained from a single computer run on an entire ribosomal RNA sequence. All predicted base pairs in optimal foldings were examined for agreement with proven base pairs in the comparative models. Our analyses show that the overall correspondence between the predicted and comparative models varied for different RNAs and ranges from a low of 27% to high of 70%, with a mean value of 49%. The correspondence improves to a mean value of 81% when the analysis is limited to well-determined helices. In addition to well-determined helices, large well-determined structural domains can be observed in "energy dot plots" of some 16S ribosomal RNAs. The predicted domains correspond closely with structural domains that are found by the comparative method in the same RNAs. Our analyses also show that measuring the agreement between predicted and comparative secondary structure models underestimates the reliability of structural prediction by mfold.

[1]  Michael Zuker,et al.  Optimal computer folding of large RNA sequences using thermodynamics and auxiliary information , 1981, Nucleic Acids Res..

[2]  R. Gutell,et al.  Detailed analysis of the higher-order structure of 16S-like ribosomal ribonucleic acids. , 1983, Microbiological reviews.

[3]  Ignacio Tinoco,et al.  A dynamic programming algorithm for finding alternative RNA secondary structures , 1986, Nucleic Acids Res..

[4]  D. Turner,et al.  Improved free-energy parameters for predictions of RNA duplex stability. , 1986, Proceedings of the National Academy of Sciences of the United States of America.

[5]  C. E. Longfellow,et al.  Improved parameters for prediction of RNA structure. , 1987, Cold Spring Harbor symposia on quantitative biology.

[6]  D. Turner,et al.  RNA structure prediction. , 1988, Annual review of biophysics and biophysical chemistry.

[7]  Jih-Hsiang Chen,et al.  A program for predicting significant RNA secondary structures , 1988, Comput. Appl. Biosci..

[8]  D. Turner,et al.  Improved predictions of secondary structures for RNA. , 1989, Proceedings of the National Academy of Sciences of the United States of America.

[9]  S. Le,et al.  A method for assessing the statistical significance of RNA folding. , 1989, Journal of theoretical biology.

[10]  M. Zuker.  On finding all suboptimal foldings of an RNA molecule. , 1989, Science.

[11]  Ross A. Overbeek,et al.  Structure detection through automated covariance search , 1990, Comput. Appl. Biosci..

[12]  D. Turner,et al.  Predicting optimal and suboptimal secondary structure for RNA. , 1990, Methods in enzymology.

[13]  J. McCaskill.  The equilibrium partition function and base pair binding probabilities for RNA secondary structure , 1990, Biopolymers.

[14]  D. Turner,et al.  Stabilities of consecutive A.C, C.C, G.G, U.C, and U.U mismatches in RNA internal loops: Evidence for stable hydrogen-bonded U.U and C.C.+ pairs. , 1991, Biochemistry.

[15]  D. Turner,et al.  A comparison of optimal and suboptimal RNA secondary structures predicted by free energy minimization with structures determined by phylogenetic comparison. , 1991, Nucleic acids research.

[16]  Ross A. Overbeek,et al.  The ribosomal database project , 1992, Nucleic Acids Res..

[17]  A Two-dimensional Model at the Nucleotide Level for the Central Hairpin of Coliphage Qβ RNA , 1993 .

[18]  M. Zuker,et al.  Structural analysis by energy dot plot of a large mRNA. , 1993, Journal of molecular biology.

[19]  R. Gutell,et al.  A compilation of large subunit (23S and 23S-like) ribosomal RNA structures: 1993. , 1992, Nucleic acids research.

[20]  R. Gutell,et al.  Lessons from an evolving rRNA: 16S and 23S rRNA structures from a comparative perspective. , 1994, Microbiological reviews.

[21]  R. Gutell,et al.  Collection of small subunit (16S- and 16S-like) ribosomal RNA structures: 1994. , 1993, Nucleic acids research.

[22]  A. E. Walter,et al.  Coaxial stacking of helixes enhances binding of oligoribonucleotides and improves predictions of RNA folding. , 1994, Proceedings of the National Academy of Sciences of the United States of America.

[23]  M. Zuker.  Prediction of RNA secondary structure by energy minimization. , 1994, Methods in molecular biology.