An essential non-Watson-Crick base pair motif in 3'UTR to mediate selenoprotein translation.

The SECIS element is an RNA hairpin in the 3'UTR of selenoprotein mRNAs required for decoding UGA selenocysteine codons. Our experimentally derived 2D structure model for the SECIS RNA revealed the conservation of four consecutive non-Watson-Crick base pairs, with a central G.A/A.G tandem. The present study was dedicated to gaining insight into the role of this quartet of base pairs. The effects of mutations introduced into the SECIS quartet of the glutathione peroxidase (GPx) cDNA, an enzyme with selenocysteine in its active center, were reported in vivo by the GPx activity. The detrimental consequence of an all-Watson-Crick mutant quartet disclosed the paramount importance of the non-Watson-Crick base pairs for GPx activity. Next, structure probing established that base pair changes in the central G.A/A.G tandem, predicted by the model to be structurally unfavorable, effectively led to local opening of the helix at the quartet. A concomitant abolition of GPx activity was observed, arising from translational impairment of full-length GPx. In contrast, an isosteric base pair replacement in the tandem did not affect base pairing in the quartet, leading to an almost wt GPx activity. Collectively, the data provided conclusive evidence for the functional relevance of these non-Watson-Crick base pairs in vivo, thus identifying a noncanonical RNA motif crucial to SECIS function in mediating selenoprotein translation. Within the quartet, the prominent requirement for the central G.A/A.G tandem is highlighted, our previous structural model and the mutagenesis data presented here strongly arguing in favor of a sheared arrangement for the G.A base pairs. The SECIS RNA is therefore another member to be added to the growing list of RNAs containing building blocks of non-Watson-Crick base pairs, required for structure and/or function.     

Analysis of RNA motifs

RNA motifs are directed and ordered stacked arrays of non-Watson-Crick base pairs forming distinctive foldings of the phosphodiester backbones of the interacting RNA strands. They correspond to the 'loops' - hairpin, internal and junction - that intersperse the Watson-Crick two-dimensional helices as seen in two-dimensional representations of RNA structure. RNA motifs mediate the specific interactions that induce the compact folding of complex RNAs. RNA motifs also constitute specific protein or ligand binding sites. A given motif is characterized by all the sequences that fold into essentially identical three-dimensional structures with the same ordered array of isosteric non-Watson-Crick base pairs. It is therefore crucial, when analyzing a three-dimensional RNA structure in order to identify and compare motifs, to first classify its non-Watson-Crick base pairs geometrically.     

Conformational features of the four successive non-Watson-Crick base pairs in RNA duplex.

A tridecaribonucleotide, r(UGAGCUUCGGCUC) doesn't form hairpin or interior loop and forms a double helix of 12 base pairs including the four successive nonstandard base pairs, U.G-U.C-C.U-G.U, in the crystal. Non-Watson-Crick base pairs, G.U and U.C are nicely incorporated in RNA duplex maintaining the regular A-form backbone. There exist the good overlapping between base pairings, U.G and U.C, so as to stabilize the nonstandard base pair track. Hydrogen bond networks involving water molecules in the major and minor grooves to stabilize this mismatch base pairing array, are observed and its conformational features are described.     

An alternating sheared AA pair and elements of stability for a single sheared purine-purine pair flanked by sheared GA pairs in RNA

A previous NMR structure of the duplex 5'GGU GGA GGCU/PCCG AAG CCG5' revealed an unusually stable RNA internal loop with three consecutive sheared GA pairs. Here, we report NMR studies of two duplexes, 5'GGU GGA GGCU/PCCA AAG CCG5' (replacing the UG pair with a UA closing pair) and 5'GGU GAA GGCU/PCCG AAG CCG5' (replacing the middle GA pair with an AA pair). An unusually stable loop with three consecutive sheared GA pairs forms in the duplex 5'GGU GGA GGCU/PCCA AAG CCG5'. The structure contrasts with that reported for this loop in the crystal structure of the large ribosomal subunit of Deinococcus radiodurans [Harms, J., Schluenzen, F., Zarivach, R., Bashan, A., Gat, S., Agmon, I., Bartels, H., Franceschi, F., and Yonath, A. (2001) Cell 107, 679-688]. The middle AA pair in the duplex 5'GGU GAA GGCU/PCCG AAG CCG5' rapidly exchanges orientations, resulting in alternative base stacking and pseudosymmetry with exclusively sheared pairs. The U GAA G/G AAG C internal loop is 2.1 kcal/mol less stable than the U GGA G/G AAG C internal loop at 37 degrees C. Structural, energetic, and dynamic consequences upon functional group substitutions within related 3 x 3 and 3 x 6 internal loops are also reported.     

RNA hairpin loop stability depends on closing base pair.

Thermodynamic parameters are reported for hairpin formation in 1 M NaCl by RNA sequences of the type GGXAUAAUAYCC, where X and Y are CG, GC, AU, UA, GU, or UG. A nearest neighbor analysis of the data indicates the free energy change for loop formation at 37 degrees C, delta degrees Gl,37, averages 3.4 kcal/mol for hairpin loops closed with C.G, G.C, and G.U pairs. In contrast, delta G degree l,37 averages 4.6 kcal/mol for loops closed with A.U, U.A, or U.G pairs. Thus the stability of an RNA hairpin depends on the closing base pair. The hairpin with a GA mismatch that is formed by GGCGUAAUAGCC is more stable than the corresponding hairpin with an AA mismatch. Thus hairpin stability also depends on loop sequence. These effects are not included in current algorithms for prediction of RNA structure from sequence.     

Motif prediction in ribosomal RNAs Lessons and prospects for automated motif prediction in homologous RNA molecules

The traditional way to infer RNA secondary structure involves an iterative process of alignment and evaluation of covariation statistics between all positions possibly involved in basepairing. Watson-Crick basepairs typically show covariations that score well when examples of two or more possible basepairs occur. This is not necessarily the case for non-Watson-Crick basepairing geometries. For example, for sheared (trans Hoogsteen/Sugar edge) pairs, one base is highly conserved (always A or mostly A with some C or U), while the other can vary (G or A and sometimes C and U as well). RNA motifs consist of ordered, stacked arrays of non-Watson-Crick basepairs that in the secondary structure representation form hairpin or internal loops, multi-stem junctions, and even pseudoknots. Although RNA motifs occur recurrently and contribute in a modular fashion to RNA architecture, it is usually not apparent which bases interact and whether it is by edge-to-edge H-bonding or solely by stacking interactions. Using a modular sequence-analysis approach, recurrent motifs related to the sarcin-ricin loop of 23S RNA and to loop E from 5S RNA were predicted in universally conserved regions of the large ribosomal RNAs (16S- and 23S-like) before the publication of high-resolution, atomic-level structures of representative examples of 16S and 23S rRNA molecules in their native contexts. This provides the opportunity to evaluate the predictive power of motif-level sequence analysis, with the goal of automating the process for predicting RNA motifs in genomic sequences. The process of inferring structure from sequence by constructing accurate alignments is a circular one. The crucial link that allows a productive iteration of motif modeling and realignment is the comparison of the sequence variations for each putative pair with the corresponding isostericity matrix to determine which basepairs are consistent both with the sequence and the geometrical data.     

The multiple flavors of GoU pairs in RNA

Wobble GU pairs (or GoU) occur frequently within double‐stranded RNA helices interspersed within the standard G═C and A─U Watson‐Crick pairs. However, other types of GoU pairs interacting on their Watson‐Crick edges have been observed. The structural and functional roles of such alternative GoU pairs are surprisingly diverse and reflect the various pairings G and U can form by exploiting all the subtleties of their electronic configurations. Here, the structural characteristics of the GoU pairs are updated following the recent crystallographic structures of functional ribosomal complexes and the development in our understanding of ribosomal translation.     

A three-dimensional RNA motif in Potato spindle tuber viroid mediates trafficking from palisade mesophyll to spongy mesophyll in Nicotiana benthamiana

Cell-to-cell trafficking of RNA is an emerging biological principle that integrates systemic gene regulation, viral infection, antiviral response, and cell-to-cell communication. A key mechanistic question is how an RNA is specifically selected for trafficking from one type of cell into another type. Here, we report the identification of an RNA motif in Potato spindle tuber viroid (PSTVd) required for trafficking from palisade mesophyll to spongy mesophyll in Nicotiana benthamiana leaves. This motif, called loop 6, has the sequence 5'-CGA-3'...5'-GAC-3' flanked on both sides by cis Watson-Crick G/C and G/U wobble base pairs. We present a three-dimensional (3D) structural model of loop 6 that specifies all non-Watson-Crick base pair interactions, derived by isostericity-based sequence comparisons with 3D RNA motifs from the RNA x-ray crystal structure database. The model is supported by available chemical modification patterns, natural sequence conservation/variations in PSTVd isolates and related species, and functional characterization of all possible mutants for each of the loop 6 base pairs. Our findings and approaches have broad implications for studying the 3D RNA structural motifs mediating trafficking of diverse RNA species across specific cellular boundaries and for studying the structure-function relationships of RNA motifs in other biological processes.     

The G x U wobble base pair. A fundamental building block of RNA structure crucial to RNA function in diverse biological systems

The G x U wobble base pair is a fundamental unit of RNA secondary structure that is present in nearly every class of RNA from organisms of all three phylogenetic domains. It has comparable thermodynamic stability to Watson-Crick base pairs and is nearly isomorphic to them. Therefore, it often substitutes for G x C or A x U base pairs. The G x U wobble base pair also has unique chemical, structural, dynamic and ligand-binding properties, which can only be partially mimicked by Watson-Crick base pairs or other mispairs. These features mark sites containing G x U pairs for recognition by proteins and other RNAs and allow the wobble pair to play essential functional roles in a remarkably wide range of biological processes.     

39 Simulation of RNA tandem GA base pairs provides insights about the forces behind conformational preference

Conformational changes are important for RNA function. We used molecular mechanics with all-atom models to understand conformational preference in RNA tandem guanine–adenine (GA) base pairs. These tandem GA base pairs play important roles in determining the stability and structural dynamics of RNA tertiary structures. Previous solution structures showed that these tandem GA base pairs adopt either imino (cis-Watson-Crick/cis-Watson-Crick interaction) or sheared (trans-Hoogsteen/trans-Hoogsteen interaction) pairing depending upon the sequence and orientation of the adjacent base pairs. In our simulations, we modeled (GCGGACGC)₂ (Wu and Turner 1996) and (GCGGAUGC)₂ (Tolbert et al., 2007), experimentally preferred as imino and sheared, respectively. Besides the experimentally preferred conformation, we constructed models of the nonnative conformations by changing cytosine to uracil or uracil to cytosine. We used explicit solvent molecular dynamics and free energy calculation with umbrella sampling to measure the free energy deference of the experimentally preferred conformation and the nonnative conformations. A modification to ff10 was required, which allowed the guanine bases’ amino group to leave the base plane (Yildirim et al., 2009). With this modification, the RMSD of unrestrained simulations and the free energy surfaces are improved, suggesting the importance of electrostatic interactions by G amino groups in stabilizing the native structures.     

NMR structure of a 4 x 4 nucleotide RNA internal loop from an R2 retrotransposon: identification of a three purine-purine sheared pair motif and comparison to MC-SYM predictions

The NMR solution structure is reported of a duplex, 5'GUGAAGCCCGU/3'UCACAGGAGGC, containing a 4 × 4 nucleotide internal loop from an R2 retrotransposon RNA. The loop contains three sheared purine-purine pairs and reveals a structural element found in other RNAs, which we refer to as the 3RRs motif. Optical melting measurements of the thermodynamics of the duplex indicate that the internal loop is 1.6 kcal/mol more stable at 37°C than predicted. The results identify the 3RRs motif as a common structural element that can facilitate prediction of 3D structure. Known examples include internal loops having the pairings: 5'GAA/3'AGG, 5'GAG/3'AGG, 5'GAA/3'AAG, and 5'AAG/3'AGG. The structural information is compared with predictions made with the MC-Sym program.     

Crystal structure of an RNA octamer duplex r(CCCIUGGG)2 incorporating tandem I.U wobbles.

The crystal structure of the RNA octamer duplex r(CCCIUGGG)2has been elucidated at 2.5 A resolution. The crystals belong to the space group P21and have unit cell constants a = 33.44 A, b = 43.41 A, c = 49.39 A and beta = 104.7 degrees with three independent duplexes (duplexes 1-3) in the asymmetric unit. The structure was solved by the molecular replacement method and refined to an Rwork/Rfree of 0.185/0.243 using 3765 reflections between 8.0 and 2.5 A. This is the first report of an RNA crystal structure incorporating I.U wobbles and three molecules in the asymmetric unit. Duplex 1 displays a kink of 24 degrees between the mismatch sites, while duplexes 2 and 3 have two kinks each of 19 degrees and 27 degrees, and 24 degrees and 29 degrees, respectively, on either side of the tandem mismatches. At the I.U/U.I mismatch steps, duplex 1 has a twist angle of 33.9 degrees, close to the average for all base pair steps, but duplexes 2 and 3 are underwound, with twist angles of 24.4 degrees and 26.5 degrees, respectively. The tandem I.U wobbles show intrastrand purine-pyrimidine stacking but exhibit interstrand purine-purine stacking with the flanking C.G pairs. The three independent duplexes are stacked non-coaxially in a head-to-tail fashion to form infinite pseudo-continuous helical columns which form intercolumn hydrogen bonding interactions through the 2'-hydroxyl groups where the minor grooves come together.     

Structure of the ribozyme substrate hairpin of Neurospora VS RNA: a close look at the cleavage site

The cleavage site of the Neurospora VS RNA ribozyme is located in a separate hairpin domain containing a hexanucleotide internal loop with an A-C mismatch and two adjacent G-A mismatches. The solution structure of the internal loop and helix la of the ribozyme substrate hairpin has been determined by nuclear magnetic resonance (NMR) spectroscopy. The 2 nt in the internal loop, flanking the cleavage site, a guanine and adenine, are involved in two sheared G.A base pairs similar to the magnesium ion-binding site of the hammerhead ribozyme. Adjacent to the tandem G.A base pairs, the adenine and cytidine, which are important for cleavage, form a noncanonical wobble A+-C base pair. The dynamic properties of the internal loop and details of the high-resolution structure support the view that the hairpin structure represents a ground state, which has to undergo a conformational change prior to cleavage. Results of chemical modification and mutagenesis data of the Neurospora VS RNA ribozyme can be explained in context with the present three-dimensional structure.     

Formation of sheared G:A base pairs in an RNA duplex modelled after ribozymes, as revealed by NMR.

The thermal stability and structure of an RNA duplex, r(GGACGAGUCC)2, the base sequence of which was modelled after both a hammerhead ribozyme and a lead ribozyme, were studied by CD and NMR. We previously demonstrated that the corresponding DNA duplex, d(GGACGAGTCC)2, formed unique 'sheared' G:A base pairs, where an amino proton, instead of an imino proton, of G is involved in the hydrogen bonding, and G and A bases are arranged 'side by side' instead of 'head to head' (Nucleic Acids Res. (1993) 21, 5418-5424). CD melting profiles showed that the RNA duplex is thermally more stable than the corresponding DNA duplex. NMR studies revealed that sheared G:A base pairs are formed in the RNA duplex, too, although the overall structure of the RNA is the A form, which differs from the B form taken on by the corresponding DNA. A model building study confirmed that sheared G:A base pairs can be accommodated in the double helical structure of the A form. A difference between the RNA and DNA duplexes in the stacking interaction involving G:A mismatch bases is also suggested. The demonstration that sheared G:A base pairs can be formed not only in DNA but also in RNA suggests that this base pairing plays an important role regarding the RNA structure.     

Isoalloxazine derivatives promote photocleavage of natural RNAs at G.U base pairs embedded within helices.

We have recently shown that isoalloxazine derivatives are able to photocleave RNA specifically at G.U base pairs embedded within a helical stack. The reaction involves the selective molecular recognition of G.U base pairs by the isoalloxazine ring and the removal of one nucleoside downstream of the uracil residue. Divalent metal ions are absolutely required for cleavage. Here we extend our studies to complex natural RNA molecules with known secondary and tertiary structures, such as tRNAs and a group I intron (td). G.U pairs were cleaved in accordance with the phylogenetically and experimentally derived secondary and tertiary structures. Tandem G.U pairs or certain G.U pairs located at a helix extremity were not affected. These new cleavage data, together with the RNA crystal structure, allowed us to perform molecular dynamics simulations to provide a structural basis for the observed specificity. We present a stable structural model for the ternary complex of the G. U-containing helical stack, the isoalloxazine molecule and a metal ion. This model provides significant new insight into several aspects of the cleavage phenomenon, mechanism and specificity for G. U pairs. Our study shows that in large natural RNAs a secondary structure motif made of an unusual base pair can be recognized and cleaved with high specificity by a low molecular weight molecule. This photocleavage reaction thus opens up the possibility of probing the accessibility of G.U base pairs, which are endowed with specific structural and functional roles in numerous structured and catalytic RNAs and interactions of RNA with proteins, in folded RNAs.     

Conformation of guanine-8-oxoadenine base pairs in the crystal structure of d(CGCGAATT(O8A)GCG).

The structure of the synthetic deoxydodecamer d(CGCGAATT(O8A)GCG)2 (O8A = 8-oxoadenine) has been determined by single-crystal X-ray diffraction techniques. The oligonucleotide crystallizes in the orthorhombic space group P2(1)2(1)2(1) with cell dimensions of a = 25.48 A, b = 41.84 A, and c = 64.91 A. The refinement has converged with an R-factor of 0.151 for 1119 reflections in the resolution range 8.0-2.25 A. Sixty-seven solvent molecules were located during the course of the refinement. The B-DNA helix consists of ten Watson-Crick base pairs and two guanine-8-oxoadenine (G.O8A) base pairs. In order to achieve hydrogen-bonding complementarity between the two bases, an unusual G(anti).O8A-(syn) wobble conformation is adopted. It is proposed that the G.O8A mispairs are held together by a network of four interbase hydrogen bonds which are the result of the formation of two reverse three-center hydrogen-bonding systems. These involve one carbonyl oxygen lone pair interacting with two hydrogen atoms. In a departure from previous observations of the characteristics of purine-purine anti-syn base pairs, lambda 1 and lambda 2, the angles between the glycosidic bonds and the C1'-C1' vector, are symmetric. A reassessment of the other purine-purine mispairs suggests that similar three-center hydrogen bonds may occur and make a contribution to stabilizing other base pairings.