A Method for Predicting Common Structures of Homologous RNAs

Abstract

We have developed a procedure, composed of a set of computer programs, for predicting common RNA structures of homologous sequences. Given a set of homologous RNAs, these programs perform a multiple sequence alignment, generate a list of possible helical stems that are thermodynamically favored in RNA folding from a selected individual sequence, establish a conserved stem list by inspecting the equivalent base pairings and/or conserved helical stems from the derived alignment of homologous RNAs, and build common RNA secondary structures with the maximum scores (i.e., compensatory base changes and number of base pairs, etc.). The approach is a combination of phylogenetic and thermodynamic methods and has been applied to the prediction of common folding structures of the 5′ untranslated regions in a number of positive RNA viruses.