Database proton NMR chemical shifts for RNA signal assignment and validation

The Biological Magnetic Resonance Data Bank contains NMR chemical shift depositions for 132 RNAs and RNA-containing complexes. We have analyzed the ¹H NMR chemical shifts reported for non-exchangeable protons of residues that reside within A-form helical regions of these RNAs. The analysis focused on the central base pair within a stretch of three adjacent base pairs (BP triplets), and included both Watson–Crick (WC; G:C, A:U) and G:U wobble pairs. Chemical shift values were included for all 4³ possible WC-BP triplets, as well as 137 additional triplets that contain one or more G:U wobbles. Sequence-dependent chemical shift correlations were identified, including correlations involving terminating base pairs within the triplets and canonical and non-canonical structures adjacent to the BP triplets (i.e. bulges, loops, WC and non-WC BPs), despite the fact that the NMR data were obtained under different conditions of pH, buffer, ionic strength, and temperature. A computer program (RNAShifts) was developed that enables convenient comparison of RNA ¹H NMR assignments with database predictions, which should facilitate future signal assignment/validation efforts and enable rapid identification of non-canonical RNA structures and RNA-ligand/protein interaction sites.     

Do the 16 mer, 5'-GUGGUCUGAUGAGGCC-3' and the 25 mer, 5'-GGCCGAAACUCGUAAGAGUCACCAC-3', form a hammerhead ribozyme structure in physiological conditions? An NMR and UV thermodynamic study

In a wide range of salt concentrations, 10-30 mM phosphate buffer containing up to 0.5 M Li2SO4 and 300 mM NaCl, 7.5 mM Mg2+, pH 5.5-7.5, a mixture of the 16 mer and the 25 mer RNA strands does not form a hammerhead in any amount detectable by NMR at 600 MHz. The imino-, amino-, aromatic- and anomeric protons in the NMR spectra of both the 16 mer and the 25 mer RNA have been assigned separately. Both the 16 mer and the 25 mer RNA both take up very stable hairpin structures, and when mixed together there is no major change of conformation in neither oligo-RNA.     

DO THE 16 MER, 5′-GUGGUCUGAUGAGGCC-3′ AND THE 25 MER, 5′-GGCCGAAACUCGUAAGAGUCACCAC-3′, FORM A HAMMERHEAD RIBOZYME STRUCTURE IN PHYSIOLOGICAL CONDITIONS? AN NMR AND UV THERMODYNAMIC STUDY

In a wide range of salt concentrations, 10–30 mM phosphate buffer containing up to 0.5 M Li₂SO₄ and 300 mM NaCl, 7.5 mM Mg²⁺, pH 5.5–7.5, a mixture of the 16 mer and the 25 mer RNA strands does not form a hammerhead in any amount detectable by NMR at 600 MHz. The imino-, amino-, aromatic- and anomeric protons in the NMR spectra of both the 16 mer and the 25 mer RNA have been assigned separately. Both the 16 mer and the 25 mer RNA both take up very stable hairpin structures, and when mixed together there is no major change of conformation in neither oligo-RNA.     

Structure of SL4 RNA from the HIV-1 packaging signal.

The NMR-based structure is described for an RNA model of stem-loop 4 (SL4) from the HIV-1 major packaging domain. The GAGA tetraloop adopts a conformation similar to the classic GNRA form, although there are differences in the details. The type II tandem G.U pairs have a combination of wobble and bifurcated hydrogen bonds where the uracil 2-carbonyl oxygen is hydrogen-bonded to both G,H1 and G,H2. There is the likelihood of a Na(+) ion coordinated to the four carbonyl oxygens in the major groove for these G.U pairs and perhaps to the N7 lone pairs of the G bases as well. A continuous stack of five bases extends over nearly the whole length of the stem to the base of the loop in the RNA 16mer: C15/U14/G13/G5/C6. There is no evidence for a terminal G.A pair; instead, G1 appears quite unrestrained, and A16 stacks on both C15 and G2. Residues G2 through G5 exhibit broadened resonances, especially G3 and U4, suggesting enhanced mobility for the 5'-side of the stem. The structure shows G2/G3/U4 stacking along the same strand, nearly isolated from interaction with the other bases. This is probably an important factor in the signal broadening and apparent mobility of these residues and the low stability of the 16mer hairpin against thermal denaturation.     

NMR studies of the structure and Mg2+ binding properties of a conserved RNA motif of EMCV picornavirus IRES element

The structure and Mg²⁺ binding properties of a conserved 75mer RNA motif of the internal ribosome entry site (IRES) element of encephalomyocarditis virus picornavirus have been investigated by ¹H-NMR and UV melting experiments. The assignment of the imino proton resonances with characteristic chemical shift dispersion for canonical and non-canonical base pairs confirmed the predicted secondary structure of the 75mer and its fragments. Addition of Mg²⁺ resulted in a dramatic increase in apparent melting temperature, with the 75mer RNA registering the biggest increase, from 63 to 80°C, thus providing evidence for enhanced stability arising from Mg²⁺ binding. Similarly, addition of Mg²⁺ induced selective changes to the chemical shifts of the imino protons of a GCGA tetraloop in the 75mer, that is essential for IRES activity, thereby highlighting a possible structural role for Mg²⁺ in the folding of the 75mer. Significantly, the same protons show retarded exchange to water solvent, even at elevated temperature, which suggest that Mg²⁺ induces a conformational rearrangement of the 75mer. Thus, we propose that Mg²⁺ serves two important roles: (i) enhancing thermodynamic stability of the 75mer RNA (and its submotifs) via non-specific interactions with the phosphate backbone and (ii) promoting the folding of the 75mer RNA by binding to the GCGA tetraloop.     

Crystal structure of an RNA duplex r(G GCGC CC)2 with non-adjacent G*U base pairs

The crystal structure of a self-complementary RNA duplex r(GGGCGCUCC)2with non-adjacent G*U and U*G wobble pairs separated by four Watson-Crick base pairs has been determined to 2.5 A resolution. Crystals belong to the space group R3; a = 33.09 A,alpha = 87.30 degrees with a pseudodyad related duplex in the asymmetric unit. The structure was refined to a final Rworkof 17.5% and Rfreeof 24.0%. The duplexes stack head-to-tail forming infinite columns with virtually no twist at the junction steps. The 3'-terminal cytosine nucleosides are disordered and there are no electron densities, but the 3' penultimate phosphates are observed. As expected, the wobble pairs are displaced with guanine towards the minor groove and uracil towards the major groove. The largest twist angles (37.70 and 40.57 degrees ) are at steps G1*C17/G2*U16 and U7*G11/C8*G10, while the smallest twist angles (28.24 and 27.27 degrees ) are at G2*U16/G3*C15 and C6*G12/U7*G11 and conform to the pseudo-dyad symmetry of the duplex. The molecule has two unequal kinks (17 and 11 degrees ) at the wobble sites and a third kink at the central G5 site which may be attributed to trans alpha (O5'-P), trans gamma (C4'-C5') backbone conformations. The 2'-hydroxyl groups in the minor groove form inter-column hydrogen bonding, either directly or through water molecules.     

Structural dynamics of double-helical RNAs composed of CUG/CUG- and CUG/CGG-repeats

Human genetic trinucleotide repeat expansion diseases (TREDs) are characterized by triplet repeat expansions, most frequently found as CNG-tracts in genome. At RNA level, such expansions suggestively result in formation of double-helical hairpins that become a potential source for small RNAs involved in RNA interference (RNAi). Here, we present three crystal structures of RNA fragments composed of triplet repeats CUG and CGG/CUG, as well as two crystal structures of same triplets in a protein-bound state. We show that both 20mer pG(CUG)(6)C and 19mer pGG(CGG)(3)(CUG)(2)CC form A-RNA duplexes, in which U·U or G·U mismatches are flanked/stabilized by two consecutive Watson-Crick G·C base pairs resulting in high-stacking GpC steps in every third position of the duplex. Despite interruption of this regularity in another 19mer, p(CGG)(3)C(CUG)(3), the oligonucleotide still forms regular double-helical structure, characterized, however, by 12?bp (rather than 11?bp) per turn. Analysis of newly determined molecular structures reveals the dynamic aspects of U·U and G·U mismatching within CNG-repetitive A-RNA and in a protein-bound state, as well as identifies an additional mode of U·U pairing essential for its dynamics and sheds the light on possible role of regularity of trinucleotide repeats for double-helical RNA structure. Findings are important for understanding the structural behavior of CNG-repetitive RNA double helices implicated in TREDs.     

The identity determinants required for the discrimination between tRNAGlu and tRNAAsp by glutamyl-tRNA synthetase from Escherichia coli.

We previously elucidated the major determinant set for Escherichia coli tRNAGlu identity (U34, U35, C36, A37, G1*C72, U2*A71, U11*A24, U13*G22**Alpha46, and Delta47) and showed that the set is sufficient to switch the identity of tRNAGln to Glu [Sekine, S., Nureki, O., Sakamoto, K., Niimi, T., Tateno, M., Go, M., Kohno, T., Brisson, A., Lapointe, J. & Yokoyama, S. (1996) J. Mol. Biol. 256, 685-700]. In the present study, we attempted to switch the identity of tRNAAsp, which has a sequence similar to that of tRNAGlu, and consequently possesses many nucleotide residues corresponding to the Glu identity determinants (U35, C36, A37, G1*C72, and U11*A24). A simple transplantation of the rest of the major determinants (U34, U2*A71, U13*G22**Alpha46, and Delta47) to the framework of tRNAAsp did not result in a sufficient switch of the tRNAAsp identity to Glu. To confer an optimal glutamate accepting activity to tRNAAsp, two other elements, C4*G69 in the middle of the acceptor stem and C12*G23**C9 in the augmented D helix, were required. Consistently, the two base pairs, C4*G69 and C12*G23, in tRNAGlu had been shown to exist in the interface with glutamyl-tRNA synthetase (GluRS) by phosphate-group footprinting. We also found the two elements in the framework of tRNAGln, and determined that their contributions successfully changed the identity of tRNAGln to Glu in the previous study. By the identity-determinant set (C4*G69 and C12*G23**C9 in addition to U34, U35, C36, A37, G1*C72, U2*A71, U11*A24, U13*G22**Alpha46, and Delta47) the activity of GluRS was optimized and efficient discrimination from the noncognate tRNAs was achieved.     

Thermodynamics of triple helix formation: spectrophotometric studies on the d(A)10.2d(T)10 and d(C+3T4C+3).d(G3A4G3).d(C3T4C3) triple helices.

We have stabilized the d(A)10.2d(T)10 and d(C+LT4C+3).d(G3A4G3).d(C3T4C3) triple helices with either NaCl or MgCl2 at pH 5.5. UV mixing curves demonstrate a 1:2 stoichiometry of purine to pyrimidine strands under the appropriate conditions of pH and ionic strength. Circular dichroic titrations suggest a possible sequence-independent spectral signature for triplex formation. Thermal denaturation profiles indicate the initial loss of the third strand followed by dissociation of the underlying duplex with increasing temperature. Depending on the base sequence and ionic conditions, the binding affinity of the third strand for the duplex at 25 degrees C is two to five orders of magnitude lower than that of the two strands forming the duplex. Thermodynamic parameters for triplex formation were determined for both sequences in the presence of 50 mM MgCl2 and/or 2.0 M NaCl. Hoogsteen base pairs are 0.22-0.64 kcal/mole less stable than Watson-Crick base pairs, depending on ionic conditions and base composition. C+.G and T.A Hoogsteen base pairs appear to have similar stability in the presence of Mg2+ ions at low pH.     

Protection of amide protons in folding intermediates of ribonuclease A measured by pH-pulse exchange curves

pH-pulse exchange curves have been measured for samples taken during the folding of ribonuclease A. The curve gives the number of protected amide protons remaining after a 10-s pulse of exchange at pHs from 6.0 to 9.5, at 10 degrees C. Amide proton exchange is base catalyzed, and the rate of exchange increases 3000-fold between pH 6.0 and pH 9.5. The pH at which exchange occurs depends on the degree of protection against exchange provided by structure. Pulse exchange curves have been measured for samples taken at three times during folding, and these are compared to the pulse exchange curves of N, the native protein, of U, the unfolded protein in 4 M guanidinium chloride, and of IN, the native-like intermediate obtained by the prefolding method of Schmid. The results are used to determine whether folding intermediates are present that can be distinguished from N and U and to measure the average degree of protection of the protected protons in folding intermediates. The amide (peptide NH) protons of unfolded ribonuclease A were prelabeled with 3H by a previous procedure that labels only the slow-folding species. Folding was initiated at pH 4.0, 10 degrees C, where amide proton exchange is slower than the folding of the slow-folding species. Samples were taken at 0-, 10-, and 20-s folding, and their pH-pulse exchange curves were measured.(ABSTRACT TRUNCATED AT 250 WORDS)     

Thermodynamics and kinetics for base-pair opening in the P1 duplex of the Tetrahymena group I ribozyme

The thermodynamics and kinetics for base-pair opening of the P1 duplex of the Tetrahymena group I ribozyme were studied by NMR hydrogen exchange experiments. The apparent equilibrium constants for base pair opening were measured for most of the imino protons in the P1 duplex using the base catalysts NH3, HPO4(2-) or TRIS. These equilibrium constants were also measured for several modified P1 duplexes, and the C-2.G23 base pair was the most stable base pair in all the duplexes. The conserved U-1*G22 base pair is required for activity of the ribozyme and the data here show that this wobble base pair destabilizes neighboring base pairs on only one side of the wobble. A 2'-OMe modification on the U-3 residue stabilized its own base pair but had little effect on the neighboring base pairs. Three base pairs, U-1*G22, C-2*G23 and A2*U21 showed unusual equilibrium constants for opening and possible implications of the opening thermodynamics of these base pairs on the undocking rates of the P1 helix with catalytic core are discussed.     

An estimate of the nearest neighbor base-pair content of 5S RNA using CD and absorption spectroscopy.

We analyzed the CD and uv absorption spectra of 5S RNA from Escherichia coli using the method developed in the preceding paper. The analysis of spectra of 5S RNA at 20 degrees C in 0.1M NaClO4, 2.5 mM Na+ (phosphate), pH 7.0, and 0.5 mM MgSO4 gave 7 +/- 3.6 A.U base pairs, 25 +/- 3.6 G.C base pairs, and 7.5 +/- 3.6 G.U base pairs. Estimates of nearest neighbor base pairs were more consistent with the Pieler-Erdmann and the Gewirth-Moore secondary structure models than with the Fox-Woese or the Burns-Luoma-Marshall models. We also examined the structure of 5S RNA as a function of temperature. The melting profile exhibited two transitions--one at about 35 degrees C and one above 50 degrees C. Our spectral data showed that helices I and II were stable during the first transition, and agreed with other data that helix III was the most likely helix to have melted. The results from this in-depth study of 5S RNA indicate that our method of analysis should be useful for studying the secondary structures of other small, unmodified RNAs.     

The identification of the A-type RNA helices in a 55mer RNA by selective incorporation of deuterium-labelled nucleotide residues (Uppsala NMR-window concept)

The 55-nt long RNA, modelling a three-way junction, with non-uniformly incorporated deuterated nucleotides has been synthesised in a pure form. The NMR-window part in this partially deuterated 55mer RNA consists of natural non-enriched nucleotide blocks at the three-way junction (shown in a square box in Fig. 2), whereas all other nucleotides of the rest of the molecule are partially deuterated (> 97 atom% 2H at C2', C3', C5', C5, and approximately 50 atom% 2H at C4'). The secondary structure of this 55mer RNA was determined by 2D 1H NOESY spectroscopy in D2O or in 10% D2O-H2O mixture. The use of deuterated building blocks in the specific region of the 55mer RNA allowed us to identify two distinct A-type RNA helices in a straightforward manner by observing connectivities of H1' with the basepaired imino and the aromatic H2 of all adenosine nucleotides as the first step for the determination of its tertiary structure in a cost- and time-effective manner without employing any 13C/15N labelling. These two decameric helices involve 40 nucleotides, for which all non-exchangeable H1', H6, H2, H8 and H5 protons (all 40 H1', all 40 H6 or H8 aromatics, all seven H2 of adenine nucleotide and all four non-deuterated H5 of cytosines) as well as all 16 exchangeable imino protons (with the exception of four terminal basepairs) and 16 amino protons of cytosines have been assigned. Since all aromatic-H2', H3' as well as H5'/5' crosspeaks from partially deuterated residues have been eliminated from the NMR spectra, the observation of natural nucleotide residues in the NMR window part has essentially been simplified. It has been found that the crosspeaks from the natural nucleotides located at the three-way junction in the NMR-window part show different degrees of line-broadening, thereby indicating that the various nucleotide residues have very different mobilities with respect to themselves as well as compared to other nucleotides in the helices. The assignment of H2' and H3' in the NMR-window part has been made based on NOESY and DQF-COSY crosspeaks. It is noteworthy that, even in this preliminary study, it has been possible to identify 10 H2' out of total 14 and 9 H3' out of 14. The data show that expanded AU containing a tract of 55mer RNA does not self-organise into a tight third helix, as the two decameric A-type helices, across the three-way junction which is evident from the absence of any additional imino protons, except those that already have been assigned for the two decameric helices.     

Dynamic behavior of the imino protons of the gamma OR3 17mer in H2O solution studied by high-resolution NMR

The imino proton resonances of gamma OR3 17mer in water were observed at 500 MHz with the time-shared Redfield pulse train. All of the 17 imino proton resonances could be assigned specifically to individual base pairs by utilizing the trace of NOE connectivities between the imino and adenine C2H protons and between imino protons themselves. AT1 and 17 showed abnormally high chemical shifts in comparison with the other AT pairs. On raising the temperature, broadening of the signal occurred in a sequential manner from the terminals except for AT10 and AT11, which were broadened at a lower temperature than GC12. The relaxation rates of the imino protons were measured by the inversion recovery method. The rates at higher temperatures represent the exchange rates of the imino protons. From the temperature dependences, activation energies of about 15 kcal/mol for the AT imino protons and 23-26 kcal/mol for the GC imino protons were obtained.     

Calorimetric unfolding of the bimolecular and i-motif complexes of the human telomere complementary strand, d(C(3)TA(2))(4)

A combination of spectroscopic and calorimetric techniques is used to determine the unfolding thermodynamics of the complexes formed by the complementary sequence of the human telomere, d(C(3)TA(2))(4), in the pH range of 4.2 to 6. Calorimetric melting curves show biphasic transitions; both transitions are shifted to higher temperatures as the pH is decreased, indicative of cytosine protonation, which favors the formation of C*C(+) base pairs. Furthermore, the transition temperature, T(M), of the lower transition depends on strand concentration, while the T(M) of the higher transition is independent of strand concentration, indicating the following sequential melting: bimolecular complex(s)-->intramolecular complex-->random coil. The thermodynamic profiles for the formation of each complex, bimolecular and i-motif reveals small favorable free energy terms resulting from favorable enthalpy-unfavorable entropy compensations, uptake of protons, marginal uptake of counterions (i-motif) and marginal release of water molecules (i-motif). Furthermore, an enthalpy of 3.2 kcal/mol (bimolecular complex) and 5.0 kcal/mol (i-motif) is estimated for a single C*C(+)/C*C(+) base-pair stack.     

Optical properties and base pairing of E. coli 5S RNA

The ultraviolet absorption, optical rotatory dispersion, circular dichroism, and infrared absorption spectra of renatured 5S RNA have been measured at pH 7.0 in 0.1M NaCl at 25° and used to obtain four independent estimates of the number of base pairs. These four estimations are in reaonable agreement and average values of 28 ± 4 G.C and 13 ± 4 A.U. base pairs.     

Multiple folding pathways for the P4-P6 RNA domain

We recently described site-specific pyrene labeling of RNA to monitor Mg(2+)-dependent equilibrium formation of tertiary structure. Here we extend these studies to follow the folding kinetics of the 160-nucleotide P4-P6 domain of the Tetrahymena group I intron RNA, using stopped-flow fluorescence with approximately 1 ms time resolution. Pyrene-labeled P4-P6 was prepared using a new phosphoramidite that allows high-yield automated synthesis of oligoribonucleotides with pyrene incorporated at a specific 2'-amino-2'-deoxyuridine residue. P4-P6 forms its higher-order tertiary structure rapidly, with k(obs) = 15-31 s(-1) (t(1/2) approximately 20-50 ms) at 35 degrees C and [Mg(2+)] approximately 10 mM in Tris-borate (TB) buffer. The folding rate increases strongly with temperature from 4 to 45 degrees C, demonstrating a large activation enthalpy DeltaH(double dagger) approximately 26 kcal/mol; the activation entropy DeltaS(double dagger) is large and positive. In low ionic strength 10 mM sodium cacodylate buffer at 35 degrees C, a slow (t(1/2) approximately 1 s) folding component is also observed. The folding kinetics are both ionic strength- and temperature-dependent; the slow phase vanishes upon increasing [Na(+)] in the cacodylate buffer, and the kinetics switch completely from fast at 30 degrees C to slow at 40 degrees C. Using synchrotron hydroxyl radical footprinting, we confirm that fluorescence monitors the same kinetic events as hydroxyl radical cleavage, and we show that the previously reported slow P4-P6 folding kinetics apply only to low ionic strength conditions. One model to explain the fast and slow folding kinetics postulates that some tertiary interactions are present even without Mg(2+) in the initial state. The fast kinetic phase reflects folding that is facilitated by these interactions, whereas the slow kinetics are observed when these interactions are disrupted at lower ionic strength and higher temperature.     

Nuclear magnetic resonance study of hydrogen-bonded ring protons in oligonucleotide helices involving classical and nonclassical base pairs.

A study of the exchangeable ring nitrogen protons in aqueous solutions of oligonucleotide complexes involving Watson-Crick base pairs as well as Hoogsteen pairs and other nonclassical hydrogen bonding schemes shows that resolvable resonances in the low-field (-10 to -16 ppm from sodium 4,4-dimethyl-4-silapentanesulfonate) region can be detected in a variety of structures other than double stranded helices. Ring nitrogen proton resonances arising from the following hydrogen-bonding situations are reported: (1) AT and GC Watson-Crick base pairs in a self-complementary octanucleotide, dApApApGpCpTpTpT; (2) U-A-U base triples in complexes between oligo-U15 and AMP; (3) C-G-C+ base triples in complexes between oligo-C17 and GMP at acid pH; (4) s4U-A-s4U base triples in complexes between oligo-s4U15 and AMP, all of which involve both Watson-Crick and Hoogsteen base pairing to form triplexes; (5) C-C+ base pairing between protonated and unprotonated C residues in oligo-C17 at acid pH; and (6) I4 base quadruples in the four strand association among oligo-I at high salt. The behavior of the dA3G-CT3 helix is consistent with both fraying of the terminal base pairs and presence of intermediate states as the helix opens. In the monomer-oligomer complexes, under the conditions used here, the exchange appears to be governed by the dissociation rate of monomer from the complex. These findings suggest that those tertiary structure hydrogen bonds in tRNA involving ring nitrogen protons should have representative resonances in the low-field (11-16 ppm) proton NMR region in H2O.     

Affinities of packaging domain loops in HIV-1 RNA for the nucleocapsid protein

To design anti-nucleocapsid drugs, it is useful to know the affinities the protein has for its natural substrates under physiological conditions. Dissociation equilibrium constants are reported for seven RNA stem-loops bound to the mature HIV-1 nucleocapsid protein, NCp7. The loops include SL1, SL2, SL3, and SL4 from the major packaging domain of genomic RNA. The binding assay is based on quenching the fluorescence of tryptophan-37 in the protein by G residues in the single-stranded loops. Tightly bound RNA molecules quench nearly all the fluorescence of freshly purified NCp7 in 0.2 M NaCl. In contrast, when the GGAG-tetraloop of tight-binding SL3 is replaced with UUCG or GAUA, quenching is almost nil, indicating very low affinity. Interpreting fluorescence titrations in terms of a rapidly equilibrating 1:1 complex explains nearly all of the experimental variance for the loops. Analyzed in this way, the highest affinities are for 20mer SL3 and 19mer SL2 hairpin constructs (K(d) = 28 +/- 3 and 23 +/- 2 nM, respectively). The 20mer stem-UUCG-loop and GAUA-loop constructs have <0.5% of the affinity for NCp7 relative to SL3. Affinities relative to SL3 for the other stem-loops are the following: 10% for a 16mer construct to model SL4, 30% for a 27mer model of the 9-residue apical loop of SL1, and 20% for a 23mer model of a 1 x 3 asymmetric internal loop in SL1. A 154mer construct that includes all four stem-loops binds tightly to NCp7, with the equivalent of three NCp7 molecules bound with high affinity per RNA; it is also possible that two strong sites and several weaker ones combine to give the appearance of three strong sites.