Incorporating chemical modification constraints into a dynamic programming algorithm for prediction of RNA secondary structure.
A dynamic programming algorithm for prediction of RNA secondary structure has been revised to accommodate folding constraints determined by chemical modification and to include free energy increments for coaxial stacking of helices when they are either adjacent or separated by a single mismatch. Furthermore, free energy parameters are revised to account for recent experimental results for terminal mismatches and hairpin, bulge, internal, and multibranch loops. To demonstrate the applicability of this method, in vivo modification was performed on 5S rRNA in both Escherichia coli and Candida albicans with 1-cyclohexyl-3-(2-morpholinoethyl) carbodiimide metho-p-toluene sulfonate, dimethyl sulfate, and kethoxal. The percentage of known base pairs in the predicted structure increased from 26.3% to 86.8% for the E. coli sequence by using modification constraints. For C. albicans, the accuracy remained 87.5% both with and without modification data. On average, for these sequences and a set of 14 sequences with known secondary structure and chemical modification data taken from the literature, accuracy improves from 67% to 76%. This enhancement primarily reflects improvement for three sequences that are predicted with <40% accuracy on the basis of energetics alone. For these sequences, inclusion of chemical modification constraints improves the average accuracy from 28% to 78%. For the 11 sequences with <6% pseudoknotted base pairs, structures predicted with constraints from chemical modification contain on average 84% of known canonical base pairs.
Mathews
David H
DH
Center for Human Genetics and Molecular Pediatric Disease, The Aab Institute of Biomedical Sciences, University of Rochester School of Medicine and Dentistry, 601 Elmwood Avenue, Box 703, Rochester, NY 14642, USA.
Disney
Matthew D
MD
Childs
Jessica L
JL
Schroeder
Susan J
SJ
Zuker
Michael
M
Turner
Douglas H
DH
D000465:Algorithms; D020137:Base Pair Mismatch; D001483:Base Sequence; D002176:Candida albicans / Q000235:genetics; D017931:DNA Primers; D004926:Escherichia coli / Q000235:genetics; D008969:Molecular Sequence Data; D009690:Nucleic Acid Conformation; D012329:RNA, Bacterial / Q000737:chemistry*; D012331:RNA, Fungal / Q000737:chemistry*
D017931:DNA Primers; D012329:RNA, Bacterial; D012331:RNA, Fungal