Search

RNA thermodynamics play an important role in determining the two- and three-dimensional structures of RNA. Internal loops of the sequence 5'-GMNU/3'-UNMG are relatively unstable thermodynamically. Here, five duplexes with GU-flanked 2 × 2 nucleotide internal loops were structurally investigated to reveal determinants of their instability. The following internal loops were investigated: 5'-GCAU/3'-UACG, 5'-UUCG/3'-GCUU, 5'-GCUU/3'-UUCG, 5'-GUCU/3'-UCUG, and 5'-GCCU/3'-UCCG. Two-dimensional nuclear magnetic resonance spectra indicate the absence of GU wobble base pairing in 5'-GCUU/3'-UUCG, 5'-GUCU/3'-UCUG, and 5'-GCCU/3'-UCCG. The 5'-GCUU/3'-UUCG loop has an unusual conformation of the GU base pairs, in which U's O2 carbonyl forms a bifurcated hydrogen bond with G's amino and imino protons. The internal loop of 5'-GUCU/3'-UCUG displays a shifted configuration in which GC pairs flank a U-U pair and several U's are in fast exchange between positions inside and outside the helix. In contrast, 5'-GCAU/3'-UACG and 5'-UUCG/3'-GCUU both have the expected GU wobble base pairs flanking the internal loop. Evidently, GU base pairs flanking internal loops are more likely to display atypical structures relative to Watson-Crick base pairs flanking internal loops. This appears to be more likely when the G of the GU pair is 5' to the loop. Such unusual structures could serve as recognition elements for biological function and as benchmarks for structure prediction methods.

RNA internal loops are often important sites for folding and function. Residues in internal loops can have pKa values shifted close to neutral pH because of the local structural environment. A series of RNA internal loops were studied at different pH by UV absorbance versus temperature melting experiments and imino proton nuclear magnetic resonance (NMR). A stabilizing CA pair forms at pH 7 in the CG/AA and CA/AA nearest neighbors when the CA pair is the first noncanonical pair (loop-terminal pair) in 3 x 3 nucleotide and larger size-symmetric internal loops. These CG/AA and CA/AA nearest neighbors, with CA adjacent to a closing Watson-Crick pair, are further stabilized when the pH is lowered from 7 to 5.5. The results are consistent with a significantly larger fraction (from approximately 20% at pH 7 to approximately 90% at pH 5.5) of adenines being protonated at the N1 position to form stabilizing wobble CA+ pairs adjacent to a sheared GA or AA pair. The noncanonical pair adjacent to the GA pair in CG/AA can either stabilize or destabilize the loop, consistent with the sequence-dependent thermodynamics of GA pairs. No significant pH-dependent stabilization is found for most of the other nearest neighbor combinations involving CA pairs (e.g., CA/AG and AG/CA), which is consistent with the formation of various nonwobble pairs observed in different local sequence contexts in crystal and NMR structures. A revised free-energy model, including stabilization by wobble CA+ pairs, is derived for predicting stabilities of medium-size RNA internal loops.

The Watson-Crick-like isoG-isoC (iGiC) pair, with the amino and carbonyl groups transposed relative to the Watson-Crick GC pair, provides an expanded alphabet for understanding interactions that shape nucleic acid structure. Here, thermodynamic stabilities of tandem GA pairs flanked by iGiC pairs are reported along with the NMR structures of the RNA self-complementary duplexes (GCiGGAiCGCA)2 and (GGiCGAiGCCA)2. A sheared GA pairing forms in (GCiGGAiCGCA)2, and an imino GA pairing forms in (GGiCGAiGCCA)2. The structures contrast with the formation of tandem imino and sheared GA pairs flanked by GC pairs in the RNA self-complementary duplexes (GCGGACGC)2 and (GGCGAGCC)2, respectively. In both iGiC duplexes, Watson-Crick-like hydrogen bonds are formed between iG and iC, and iGiC substitutions result in less favorable loop stability. The results provide benchmarks for testing computations of molecular interactions that shape RNA three-dimensional structure.

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.

Internal loops in RNA are important for folding and function. Consecutive noncanonical pairs can form in internal loops having at least two nucleotides on each side. Thermodynamic and structural insights into such internal loops should improve approximations for their stabilities and predictions of secondary and three-dimensional structures. Most natural internal loops are purine rich. A series of oligoribonucleotides that form purine-rich internal loops of 5-10 nucleotides, including kink-turn loops, were studied by UV melting, exchangeable proton and phosphorus NMR. Three consecutive GA pairs with the motif 5' Y GGA/3' R AAG or GGA R 3'/AAG Y 5' (i.e., 5' GGA 3'/3' AAG 5' closed on at least one side with a CG, UA, or UG pair with Y representing C or U and R representing A or G) stabilize internal loops having 6-10 nucleotides. Certain motifs with two consecutive GA pairs are also stabilizing. In internal loops with three or more nucleotides on each side, the motif 5' U G/3' G A has stability similar to 5' C G/3' G A. A revised model for predicting stabilities of internal loops with 6-10 nucleotides is derived by multiple linear regression. Loops with 2 x 3 nucleotides are predicted well by a previous thermodynamic model.

Internal loops in RNA are important for folding and function. Many folding motifs are internal loops containing GA base pairs, which are usually thermodynamically stabilizing, i.e., contribute favorable free energy to folding. Understanding the sequence dependence of folding stability and structure in terms of molecular interactions, such as hydrogen bonding and base stacking, will provide a foundation for predicting stability and structure. Here, we report the NMR structure of the oligonucleotide duplex, 5'GGUGGAGGCU3'/3'PCCGAAGCCG5' (P = purine), containing an unusually stable and relatively abundant internal loop, 5'GGA3'/3'AAG5'. This loop contains three consecutive sheared GA pairs (trans Hoogsteen/Sugar edge AG) with separate stacks of three G's and three A's in a row. The thermodynamic consequences of various nucleotide substitutions are also reported. Significant destabilization of approximately 2 kcal/mol at 37 degrees C is found for substitution of the middle GA with AA to form 5'GAA3'/3'AAG5'. This destabilization correlates with a unique base stacking and hydrogen-bonding network within the 5'GGA3'/3'AAG5' loop. Interestingly, the motifs, 5'UG3'/3'GA5' and 5'UG3'/3'AA5', have stability similar to 5'CG3'/3'GA5' even though UG and UA pairs are usually less stable than CG pairs. Consecutive sheared GA pairs in the 5'GGA3'/3'AAG5' loop are preorganized for potential tertiary interactions and ligand binding.

Many internal loops that form tertiary contacts in natural RNAs have GU closing pairs; examples include the tetraloop receptor and P1 helix docking site in group I introns. Thus, thermodynamic parameters of internal loops with GU closing pairs can contribute to the prediction of both secondary and tertiary structure. Oligoribonucleotide duplexes containing small internal loops with GU closing pairs were studied by optical melting, one-dimensional imino proton NMR, and one-dimensional phosphorus NMR. The thermodynamic stabilities of asymmetric internal loops with GU closing pairs relative to those of loops with GC closing pairs may be explained by hydrogen bonds. In contrast, the free energy increments for symmetric internal loops of two noncanonical pairs with GU closing pairs relative to loops with GC closing pairs show much more sequence dependence. Imino proton and phosphorus NMR spectra suggest that some GA pairs adjacent to GU closing pairs may form an overall thermodynamically stable but non-A-form conformation.

The symmetric, tandem GU mismatch motifs, and , which only differ in the mismatch order, have an average difference in thermodynamic stability of 2 kcal/mol at 37 degrees C. Thermodynamic studies of duplexes containing these motifs indicate the effect is largely localized to the mismatches and adjacent base pairs. The three-dimensional structures of two representative duplexes, (rGGAGUUCC)2 and (rGGAUGUCC)2, were determined by two-dimensional NMR and a simulated annealing protocol. Local deviations are similar to other intrahelical GU mismatches with little effect on backbone torsion angles and a slight overtwisting between the base pair 5' of the G of the mismatch and the mismatch itself. Comparisons of the resulting stacking patterns along with electrostatic potential maps suggest that interactions between highly negative electrostatic regions between base pairs may play a role in the observed thermodynamic differences.

There is a phylogenetically conserved G.U pair at the 5'-splice site of group I introns. When this is mutagenized to a G-C pair, splicing of these introns is greatly reduced. We have used a ribozyme derived from the Tetrahymena group I intron to compare the binding and reactivity of oligonucleotides that form either a G.U or a G-C pair at this position. Ribozyme binding of oligonucleotides at 42 degrees C was measured by native gel electrophoresis and equilibrium dialysis. Binding of GGCCCUCC (C(-1)P), which base-pairs with the ribozyme guide sequence to form a G-C at the cleavage site, was 10-fold weaker than the binding of GGCCCUCU (U(-1)P), which maintains the conserved G.U pair at the cleavage site. This is surprising since a terminal G-C enhances the binding between oligonucleotides by 20-fold relative to a terminal G.U. Thermal denaturation studies indicate that C(-1)P and several analogs with deoxy substitutions bind the guide-sequence oligonucleotide, GGAGGGAAA, as strongly as they bind the ribozyme. In contrast, U(-1)P binds 240-fold more strongly to the ribozyme than to GGAGGGAAA, a difference that is decreased by deoxy substitutions. Thus, while U(-1)P binds the ribozyme through a combination of base-pairing and specific 2-OH and other tertiary interactions, C(-1)P may bind by base-pairing alone. The substrate GGCCCUCCAAAAA (C(-1)S) is cleaved 100-fold more slowly than GGCCCUCUAAAAA (U(-1)S) and also has a higher propensity to be cleaved at the wrong nucleotide position.(ABSTRACT TRUNCATED AT 250 WORDS)

UV melting and imino proton NMR studies show that the stabilities and structures of tandem GA mismatches in RNA are dependent upon the closing base pairs around these mismatches. Internal loops of sequence 5'XGAY3'3'YAGXS' and 5'XAGY3'3'YGAX5' in the middle of octanucleotides have a range of stabilities over 5 kcal/mol when XY is a Watson-Crick or GU pair. The order of stabilities for these internal loops is 5'-GGAC-3' > UGAG, CGAG > AGAU > UGAA > GGAU. The motifs GGAC, UGAG, and CGAG are stabilizing, while the other GA motifs are destabilizing. The GAGC motif is more stable than CAGG and CGAG, but less stable than GGAC. Chemical shifts for imino protons suggest that the G imino proton of each GA mismatch in 5'-GGAC-3', 5'-GAGC-3', and 5'-CAGG-3' [SantaLucia, J., Jr., Kierzek, R., & Turner, D. H. (1990) Biochemistry 29, 8813-8819] is involved in a hydrogen bond to the base A, whereas in other 5'-XGAY-3' sequences, it is not involved in a hydrogen bond to the base A.

Thermodynamic parameters derived from optical melting studies are reported for duplex formation by a series of oligoribonucleotides containing G.U mismatches. The results are used to determine nearest-neighbor parameters for helix propagation by G.U mismatches. Surprisingly, the [formula; see text] nearest-neighbor free energy increment in unfavorable in the contexts [formula; see text], and [formula; see text] but favorable in the context [formula; see text]. This is a non-nearest-neighbor effect. In contrast, the [formula; see text] free energy increment is favorable and independent of context. Circular dichroism and imino proton NMR spectra of several sequences do not reveal an obvious structural basis for this dichotomy. For example, all the G.U mismatches have two slowly exchanging imino protons. The imino resonances for the G.U mismatches in GGAGUUCC, GUCGUGAC, and CCUGUAGG, however, broaden at lower temperature than the imino resonances for the interior Watson-Crick base pairs. In contrast, the imino resonances for the G.U mismatches in GGAUGUCC remain sharp at high temperature. The improved parameters for G.U mismatches should improve predictions of RNA structure from sequence.

UV melting, CD and NMR studies indicate rGCGAGCG and rGCAGGCG from unusually stable duplexes of type a and b. The observed delta G degree 37's in 1 M NaCl are -6.7 and -6.3 kcal/mol, respectively. For the related duplex, c, delta G degree 37 is -4.2 kcal/mol. The predicted delta G degree 37 from nearest-neighbor parameters (formula; see text) for all three duplexes is -4.7 kcal/mol (Freier, S.M., Kierzek, R., Jaeger, J.A., Sugimoto, N., Caruthers, M.H., Neilson, T., & Turner, D.H. (1986) Proc. Natl. Acad. Sci. U.S.A. 83, 9373-9377). The results suggest a special interaction in the duplexes containing GA mismatches. Presumably, this is hydrogen bonding between G and A. While the thermodynamics for (rGCGAGCG)2 and (rGCAGGCG)2 are similar, CD and the imino region of the proton NMR spectra indicate their structures are different. In particular, (rGCAGGCG)2 exhibits a CD spectrum typical of A-form geometry with a weak negative band at 280 nm. In contrast, the CD spectrum for (rGCGAGCG)2 has an intense positive band at 285 nm. The NMR spectrum of (rGCAGGCG)2 has a resonance corresponding to a hydrogen-bonded GA mismatch, while for (rGCGAGCG)2 no hydrogen-bonded imino proton is observed for the mismatch. The glycosidic torsion angles of the bases in the GA mismatches of (rGCAGGCG)2 and (rCGCAGGCG)2 are anti. Duplexes of type d, where X is A, G, or U, are more stable than e, and the stability differences are similar to those (formula; see text) observed for f versus g. Thus, 3'-dangling ends in this system make contributions to duplex stability that are similar to contributions observed with fully paired duplexes.

Thermodynamic parameters of helix formation were measured spectroscopically for 16 oligoribonucleotides containing either internal GU mismatches or the corresponding AU pairs. Internal GU mismatches stabilize each helix, but not as much as the corresponding AU pairs. The differences in the enthalpy and entropy changes of helix formation associated with replacing AU pairs with GU mismatches are less than previously realized. At both 25 and 37 degrees C, the decrease in helix stability associated with replacing an AU with a GU is also less than thought previously. Approximations are suggested for predicting the effects of GU mismatches on helix stability.

Thermodynamic parameters of helix formation were measured spectroscopically for seven hexaribonucleotides containing a GC tetramer core and G.U or other terminal mismatches. The free energies of helix formation are compared with those for the tetramer core alone and with those for the hexamer with six Watson-Crick base pairs. In 1 M NaCl, at 37 degrees C, the free energy of a terminal G.U mismatch is about equal to that of the corresponding A.U pair. Although other terminal mismatches studied add between -1.0 and -1.6 kcal/mol to delta G0 37 for helix formation, all are less stable than the corresponding Watson-Crick pairs. Comparisons of the stability increments for terminal G.U mismatches and G.C pairs suggest when stacking is weak the additional hydrogen bond in the G.C pair adds roughly -1 kcal/mol to the favorable free energy of duplex formation.