Identification of base-triples in RNA using comparative sequence analysis

Abstract

Comparative sequence analysis has proven to be a very efficient tool for the determination of RNA secondary structure and certain tertiary interactions. However, base-triples, an important RNA structural element, cannot be predicted accurately from sequence data. We show here that the poor base correlations observed at base-triple positions are the result of two factors. (1) Base covariation is not as strictly required in triples as it is in Watson-Crick pairs. (2) Base-triple structures are less conserved among homologous molecules. A particularity of known triple-helical regions is the presence of multiple base correlations that do not reflect direct pairing. We suggest that natural mutations in base-triples create structural changes that require compensatory mutations in adjacent base-pairs and triples to maintain the triple-helix conformation. On the basis of these observations, we devised two new measures of association that significantly enhance the base-triple signal in correlation studies. We evaluated correlations between base-pairs and single stranded bases, and correlations between adjacent base-pairs. Positions that score well in both analyses are the best triple candidates. This procedure correctly identifies triples, or interactions very close to the proposed triples, in type I and type II tRNAs and in the group I intron.