Pattern analysis of RNA secondary structure similarity and consensus of minimal-energy folding

We describe an automated procedure to search for consensus structures or substructures in a set of homologous or related RNA molecules. The procedure is based on the calculation of optimal and sub-optimal secondary structures using thermodynamic rules for base-pairing by energy-minimization. A linear representation of the secondary structures of the related RNAs is used so that they can be compared and classified using standard alignment and clusterings programs. We illustrate the method by means of two sets of homologous small RNAs, U2 and U3, and a set of alpha-globin mRNAs and show that biologically interesting consensus structures are obtained.