Diversity of base-pair conformations and their occurrence in rRNA structure and RNA structural motifs

In addition to the canonical base-pairs comprising the standard Watson-Crick (C:G and U:A) and wobble U:G conformations, an analysis of the base-pair types and conformations in the rRNAs in the high-resolution crystal structures of the Thermus thermophilus 30S and Haloarcula marismortui 50S ribosomal subunits has identified a wide variety of non-canonical base-pair types and conformations. However, the existing nomenclatures do not describe all of the observed non-canonical conformations or describe them with some ambiguity. Thus, a standardized system is required to classify all of these non-canonical conformations appropriately. Here, we propose a new, simple and systematic nomenclature that unambiguously classifies base-pair conformations occurring in base-pairs, base-triples and base-quadruples that are associated with secondary and tertiary interactions. This system is based on the topological arrangement of the two bases and glycosidic bonds in a given base-pair. Base-pairs in the internal positions of regular secondary structure helices usually form with canonical base-pair groups (C:G, U:A, and U:G) and canonical conformations (C:G WC, U:A WC, and U:G Wb). In contrast, non-helical base-pairs outside of regular structure helices usually have non-canonical base-pair groups and conformations. In addition, many non-helical base-pairs are involved in RNA motifs that form a defined set of non-canonical conformations. Thus, each rare non-canonical conformation may be functionally and structurally important. Finally, the topology-based isostericity of base-pair conformations can rationalize base-pair exchanges in the evolution of RNA molecules.     

Dynalign: an algorithm for finding the secondary structure common to two RNA sequences

With the rapid increase in the size of the genome sequence database, computational analysis of RNA will become increasingly important in revealing structure-function relationships and potential drug targets. RNA secondary structure prediction for a single sequence is 73 % accurate on average for a large database of known secondary structures. This level of accuracy provides a good starting point for determining a secondary structure either by comparative sequence analysis or by the interpretation of experimental studies. Dynalign is a new computer algorithm that improves the accuracy of structure prediction by combining free energy minimization and comparative sequence analysis to find a low free energy structure common to two sequences without requiring any sequence identity. It uses a dynamic programming construct suggested by Sankoff. Dynalign, however, restricts the maximum distance, M, allowed between aligned nucleotides in the two sequences. This makes the calculation tractable because the complexity is simplified to O(M(3)N(3)), where N is the length of the shorter sequence.The accuracy of Dynalign was tested with sets of 13 tRNAs, seven 5 S rRNAs, and two R2 3' UTR sequences. On average, Dynalign predicted 86.1 % of known base-pairs in the tRNAs, as compared to 59.7 % for free energy minimization alone. For the 5 S rRNAs, the average accuracy improves from 47.8 % to 86.4 %. The secondary structure of the R2 3' UTR from Drosophila takahashii is poorly predicted by standard free energy minimization. With Dynalign, however, the structure predicted in tandem with the sequence from Drosophila melanogaster nearly matches the structure determined by comparative sequence analysis.     

Solution structure of a DNA double helix incorporating four consecutive non-Watson-Crick base-pairs

A series of DNA 21-mers containing a variety of the 4 x 4 internal loop sequence 5'-CAAG-3'/3'-ACGT-5' were studied using nuclear magnetic resonance (NMR) methodology and distance geometry (DG)/molecular dynamics (MD) approaches. Such oligomers exhibit excellent resolution in the NMR spectra and reveal many unusual NOEs (nuclear Overhauser effect) that allow for the detailed characterization of a DNA hairpin incorporating a track of four different non-Watson-Crick base-pairs in the stem. These include a wobble C.A base-pair, a sheared A.C base-pair, a sheared A.G base-pair, and a wobble G.T base-pair. Significantly different twisting angles were observed between the base-pairs in internal loop that results with excellent intra-strand and inter-strand base stacking within the four consecutive mismatches and the surrounding canonical base-pairs. This explains why it melts at 52 degrees C even though five out of ten base-pairs in the stem adopt non-Watson-Crick pairs. However, the 4 x 4 internal loop still fits into a B-DNA double helix very well without significant change in the backbone torsion angles; only zeta torsion angles between the tandem sheared base-pairs are changed to a great extent from the gauche(-) domain to the trans domain to accommodate the cross-strand base stacking in the internal loop. The observation that several consecutive non-canonical base-pairs can stably co-exist with Watson-Crick base-pairs greatly increases the limited repertoire of irregular DNA folds and reveals the possibility for unusual structural formation in the functionally important genomic regions that have potential to become single-stranded.     

Structures of Non-canonical Tandem Base Pairs in RNA Helices: Review

Abstract

The structures of tandem non-canonical base pairs, a frequently recurring motif in RNA molecules, are reviewed and analysed. The tandem non-canonical base pair motifs can be roughly divided in three groups, containing seven subgroups based on their base pairing patterns and local geometries. Structural details and helical parameters that can be used to numerically distinguish between the subgroups are tabulated. Remarkably, while the individual helical twists of the tandem and adjacent base pair steps can be substantially smaller or larger than the typical A-form value of 32.7°, the average value is close to A-form. This and other striking regularities resulting from compensating geometrical adjustments, important for understanding and predicting the configurations of non-canonical base pairs geometries are discussed.     

Structure and thermodynamics of metal binding in the P5 helix of a group I intron ribozyme.

The solution structure of an RNA hairpin modelling the P5 helix of a group I intron, complexed with Co(NH3)63+, has been determined by nuclear magnetic resonance. Co(NH3)63+, which possesses a geometry very close to Mg(H2O)62+, was used to identify and characterize a Mg2+binding site in the RNA. Strong and positive intermolecular nuclear Overhauser effect (NOE) cross-peaks define a specific complex in which the Co(NH3)63+molecule is in the major groove of tandem G.U base-pairs. The structure of the RNA is characterized by a very low twist angle between the two G.U base-pairs, providing a flat and narrowed major groove. The Co(NH3)63+, although highly localized, is free to rotate to hydrogen bond in several ways to the O4 atoms of the uracil bases and to N7 and O6 of the guanine bases. Negative and small NOE cross-peaks to other protons in the sequence reveal a non-specific or delocalized interaction, characterized by a high mobility of the cobalt ion. Mn2+titrations of P5 show specific broadening of protons of the G.U base-pairs that form the metal ion binding site, in agreement with the NOE data from Co(NH3)63+. Binding constants for the interaction of Co(NH3)63+and of Mg2+to P5 were determined by monitoring imino proton chemical shifts during titration of the RNA with the metal ions. Dissociation constants are on the order of 0.1 mM for Co(NH3)63+and 1 mM for Mg2+. Binding studies were done on mutants with sequences corresponding to the three orientations of tandem G.U base-pairs. The affinities of Co(NH3)63+and Mg2+for the tandem G.U base-pairs depend strongly on their sequences; the differences can be understood in terms of the different structures of the corresponding metal ion-RNA complexes. Substitution of G.C or A.U for G.U pairs also affected the binding, as expected. These structural and thermodynamic results provide systematic new information about major groove metal ion binding in RNA.     

The 2 A structure of helix 6 of the human signal recognition particle RNA

BACKGROUND: The mammalian signal recognition particle (SRP) is an essential cytoplasmic ribonucleoprotein complex involved in targeting signal-peptide-containing proteins to the endoplasmic reticulum. Assembly of the SRP requires protein SRP19 to bind first to helix 6 of the SRP RNA before the signal-peptide-recognizing protein, SRP54, can bind to helix 8 of the RNA. Helix 6 is closed by a GGAG tetraloop, which has been shown to form part of the SRP19-binding site. RESULTS: The high-resolution (2.0 A) structure of a fragment of human SRP RNA comprising 29 nucleotides of helix 6 has been determined using the multiple anomalous dispersion (MAD) method and bromine-labelled RNA. In the crystal the molecule forms 28-mer duplexes rather than the native monomeric hairpin structure, although two chemically equivalent 11 base pair stretches of the duplex represent the presumed native structure. The duplex has highly distorted A-RNA geometry caused by the occurrence of several non-Watson-Crick base pairs. These include a 5'-GGAG-3'/3'-GAGG-5' purine bulge (which replaces the tetraloop) and a 5'-AC-3'/3'-CA-5' tandem mismatch that, depending on the protonation state of the adenine bases, adopts a different conformation in the two native-like parts of the structure. The structure also shows the 2'3'-cyclic phosphate reaction product of the hammerhead ribozyme cleavage reaction. CONCLUSIONS: The 29-mer RNA is the first RNA structure of the human SRP and provides some insight into the binding mode of SRP19. The observed strong irregularities of the RNA helix make the major groove wide enough and flat enough to possibly accommodate an alpha helix of SRP19. The variety of non-canonical base pairs observed enlarges the limited repertoire of irregular RNA folds known to date and the observed conformation of the 2'3'-cyclic phosphate containing Ade29 is consistent with the current understanding of the hammerhead ribozyme reaction mechanism.     

Nudged elastic band calculation of minimal energy paths for the conformational change of a GG non-canonical pair

The nudged elastic band (NEB) technique has been implemented in AMBER to calculate low-energy paths for conformational changes. A novel simulated annealing protocol that does not require an initial hypothesis for the path is used to sample low-energy paths. This was used to study the conformational change of an RNA cis Watson-Crick/Hoogsteen GG non-canonical pair, with one G syn around the glycosidic bond and the other anti. A previous solution structure, determined by NMR-constrained modeling, demonstrated that the GG pairs change from (syn)G-(anti)G to (anti)G-(syn)G in the context of duplex r(GCAGGCGUGC) on the millisecond timescale. The set of low-energy paths found by NEB show that each G flips independently around the glycosidic bond, with the anti G flipping to syn first. Guanine bases flip without opening adjacent base-pairs by protruding into the major groove, accommodated by a transient change by the ribose to C2'-exo sugar pucker. Hydrogen bonds between bases and the backbone, which lower the energetic barrier to flipping, are observed along the path. The results show the plasticity of RNA base-pairs in helices, which is important for biological processes, including mismatch repair, protein recognition, and translation. The modeling of the GG conformational change also demonstrates that NEB can be used to discover non-trivial paths for macromolecules and therefore NEB can be used as an exploratory method for predicting putative conformational change paths.     

AA.AG@helix.ends: A:A and A:G base-pairs at the ends of 16 S and 23 S rRNA helices.

This study reveals that AA and AG oppositions occur frequently at the ends of helices in RNA crystal and NMR structures in the PDB database and in the 16 S and 23 S rRNA comparative structure models, with the G usually 3' to the helix for the AG oppositions. In addition, these oppositions are frequently base-paired and usually in the sheared conformation, although other conformations are present in NMR and crystal structures. These A:A and A:G base-pairs are present in a variety of structural environments, including GNRA tetraloops, E and E-like loops, interfaced between two helices that are coaxially stacked, tandem G:A base-pairs, U-turns, and adenosine platforms. Finally, given structural studies that reveal conformational rearrangements occurring in regions of the RNA with AA and AG oppositions at the ends of helices, we suggest that these conformationally unique helix extensions might be associated with functionally important structural rearrangements.     

Compensatory mutations demonstrate that P8 and P6 are RNA secondary structure elements important for processing of a group I intron.

Compensatory mutations have been constructed which demonstrate that P8 and P6, two of nine proposed base-pairing interactions characteristic of group I introns, exist within the folded structure of the Tetrahymena thermophila rRNA intervening sequence, and that these secondary structure elements are important for splicing in E. coli and self-splicing in vitro. Two-base mutations in the 5' and 3' segments of P8 are predicted to disrupt P8 and a strong splicing-defective phenotype is observed in each case. A compensatory four-base mutation in P8 is predicted to restore pairing, and results in the restoration of splicing activity to nearly wild type levels. Thus, we conclude that P8 exists and is essential for splicing. In contrast to the strong phenotypes generally exhibited by mutations which disrupt RNA secondary structure, a two-base mutation in L8, the loop between P8[5'] and P8[3'], results in only a slight decrease in splicing activity. We also tested P6, a pairing which is proposed to consist of only two base-pairs in this intron. A two-base mutation in P6[3'] reduces splicing activity to a greater extent than does a two-base mutation in P6[5']. Comparison of the activities of these mutants and a compensatory P6 four-base mutant support the existence of P6, and suggest that the P6 pairing may be particularly important in the exon ligation step of splicing.     

The UAA/GAN internal loop motif: a new RNA structural element that forms a cross-strand AAA stack and long-range tertiary interactions

Analysis of aligned RNA sequences and high-resolution crystal structures has revealed a new RNA structural element, termed the UAA/GAN motif. Found in internal loops of the 23 S rRNA, as well as in RNase P RNA and group I and II introns, this six-nucleotide motif adopts a distinctive local structure that includes two base-pairs with non-canonical conformations and three conserved adenine bases, which form a cross-strand AAA stack in the minor groove. Most importantly, the motif invariably forms long-range tertiary contacts, as the AAA stack typically forms A-minor interactions and the flipped-out N nucleotide forms additional contacts that are specific to the structural context of each loop. The widespread presence of this motif and its propensity to form long-range contacts suggest that it plays a critical role in defining the architectures of structured RNAs.     

Cycling of the Sm-like protein Hfq on the DsrA small regulatory RNA

Small RNAs (sRNAs) regulate bacterial genes involved in environmental adaptation. This RNA regulation requires Hfq, a bacterial Sm-like protein that stabilizes sRNAs and enhances RNA-RNA interactions. To understand the mechanism of target recognition by sRNAs, we investigated the interactions between Hfq, the sRNA DsrA, and its regulatory target rpoS mRNA, which encodes the stress response sigma factor. Nuclease footprinting revealed that Hfq recognized multiple sites in rpoS mRNA without significantly perturbing secondary structure in the 5' leader that inhibits translation initiation. Base-pairing with DsrA, however, made the rpoS ribosome binding site fully accessible, as predicted by genetic data. Hfq bound DsrA four times more tightly than the DsrA.rpoS RNA complex in gel mobility-shift assays. Consequently, Hfq is displaced rapidly from its high-affinity binding site on DsrA by conformational changes in DsrA, when DsrA base-pairs with rpoS mRNA. Hfq accelerated DsrA.rpoS RNA association and stabilized the RNA complex up to twofold. Hybridization of DsrA and rpoS mRNA was optimal when Hfq occupied its primary binding site on free DsrA, but was inhibited when Hfq associated with the DsrA.rpoS RNA complex. We conclude that recognition of rpoS mRNA is stimulated by binding of Hfq to free DsrA sRNA, followed by release of Hfq from the sRNA.mRNA complex.     

Recognition and cleavage of the bacteriophage P1 packaging site (pac). II. Functional limits of pac and location of pac cleavage termini

Bacteriophage P1 initiates the processive packaging of its DNA at a unique site called pac. We show that a functional pac site is contained within a 161 base-pair segment of P1 EcoRI fragment 20. It extends from a position 71 base-pairs to a position 232 base-pairs from the EcoRI-22 proximal side of that fragment. The 3' and 5' pac termini are located centrally within that 161 base-pair region and are distributed over about a turn of the DNA helix. The DNA sequence of the terminus region is shown below, with the large arrows indicating the positions of termini that are frequently represented in the PI population and the small arrows indicating the positions of termini that are rarely represented in the P1 population. (Sequence: in text). Digestion of P1 virus DNA with EcoRI generates two major EcoRI-pac fragments, which differ in size by about five or six base-pairs. While the structure and position of the double-stranded pac ends of these fragments have not been determined precisely, the 5' termini at those ends probably correspond to the two major pac cleavage sites in the upper strand of the sequences shown above. The 161 base-pair pac site contains the hexanucleotide sequence 5'-TGATCAG-3' repeated four times at one end and three times at the other. Removal of just one of those elements from either the right or left ends of pac reduces pac cleavage by about tenfold. Moreover, the elements appear to be additive in their effect on pac cleavage, as removal of one and a half elements or all three elements from the right side of pac reduces pac cleavage 100-fold, and greater than 1000-fold, respectively.     

Drosophila telomeres: A variation on the telomerase theme

In Drosophila, the role of telomerase is carried out by three specialized retrotransposable elements, HeT-A, TART and TAHRE. Telomeres contain long tandem head-to-tail arrays of these elements. Within each array, the three elements occur in random, but polarized, order. Some are truncated at the 5' end, giving the telomere an enriched content of the large 3' untranslated regions, which distinguish these telomeric elements from other retrotransposons. Thus, Drosophila telomeres resemble other telomeres because they are long arrays of repeated sequences, albeit more irregular arrays than those produced by telomerase. The telomeric retrotransposons are reverse-transcribed directly onto the end of the chromosome, extending the end by successive transpositions. Their transposition uses exactly the same method by which telomerase extends chromosome ends--copying an RNA template. In addition to these similarities in structure and maintenance, Drosophila telomeres have strong functional similarities to other telomeres and, as variants, provide an important model for understanding general principles of telomere function and evolution.     

Efficient use of a 'dead-end' GA 5' splice site in the human fibroblast growth factor receptor genes

We have investigated use of a conserved non-canonical GA 5' splice site present in vertebrate fibroblast growth factor receptor (FGFR) genes. Despite previous studies suggesting that GA at the beginning of an intron is incompatible with splicing, we observe efficient utilization of this splice site for human FGFR1 gene constructs. We show that use of the GA splice site is dependent on both a conventional splice site six nucleotides upstream and sequence elements within the downstream intron. Furthermore, our results are consistent with competition between the tandem 5' splice sites being mediated by U6 snRNP, rather than U1 snRNP. Thus the GA 5' splice site represents an extension of the adjacent conventional 5' splice site, the first natural example of such a composite 5' splice site.