Assignment and NOE analysis of 2'-hydroxyl protons in RNA: implications for stabilization of RNA A-form duplexes

Synthesis and properties of fully-modified 4'-selenoRNA, an endonuclease-resistant RNA analog

A large number of chemically modified oligonucleotides (ONs) have been developed for RNA-based technologies. In most modified RNAs, the characteristic 2'-hydroxyl (2'-OH) groups are removed to enhance both nuclease resistance and hybridization ability. However, the importance of the 2'-OH group in RNA structure and function is well known. Here, we report the synthesis and properties of 4'-selenoRNA in which all four nucleoside units retain the 2'-OH groups but contain a selenium atom instead of an oxygen atom at the 4'-position of the furanose ring. 4'-SelenoRNA has enhanced ability to form duplexes with RNA, and high endonuclease resistance despite the presence of the 2'-OH groups. X-ray crystallography analysis showed that the 4'-selenoRNA duplex adopts an A-conformation, similar to natural RNA, although one 4'-selenocytidine residue has unusual South-type sugar puckering. Furthermore, preliminary studies using 4'-seleno-modified siRNAs suggest that 4'-selenoRNA may be applicable to RNA interference technology. Collectively, our results raise the possibility of a new class of modified RNA in which 2'-OH groups do not need to be masked.

Sensitivity enhancement of homonuclear multidimensional NMR correlations for labile sites in proteins, polysaccharides, and nucleic acids

Multidimensional TOCSY and NOESY are central experiments in chemical and biophysical NMR. Limited efficiencies are an intrinsic downside of these methods, particularly when targeting labile sites. This study demonstrates that the decoherence imparted on these protons through solvent exchanges can, when suitably manipulated, lead to dramatic sensitivity gains per unit time in the acquisition of these experiments. To achieve this, a priori selected frequencies are encoded according to Hadamard recipes, while concurrently subject to looped selective inversion or selective saturation procedures. Suitable processing then leads to protein, oligosaccharide and nucleic acid cross-peak enhancements of ≈200–1000% per scan, in measurements that are ≈10-fold faster than conventional counterparts. The extent of these gains will depend on the solvent exchange and relaxation rates of the targeted sites; these gains also benefit considerably from the spectral resolution provided by ultrahigh fields, as corroborated by NMR experiments at 600 MHz and 1 GHz. The mechanisms underlying these experiments’ enhanced efficiencies are analyzed on the basis of three-way polarization transfer interplays between the water, labile and non-labile protons, and the experimental results are rationalized using both analytical and numerical derivations. Limitations as well as further extensions of the proposed methods, are also discussed.

Here, the authors present an approach that enhances the sensitivity of basic 2D biomolecular NMR experiments like NOESY and TOCSY, when carried out in polysaccharides, proteins and nucleic acids. This method combines principles associated to quantum Anti-Zeno Effects and advanced data acquisition methods based on Hadamard multiplexing.

Intrinsic Contributions of 2'-Hydroxyl to the Hydration of Nucleosides at the Monomeric Level

Although many reports have revealed structural differences between DNA and RNA at the polymeric level, there are no comparative studies with 2'-deoxyribonucleoside and ribonucleoside to explore the role of the 2'-OH group at the monomeric level under the same conditions. Inspired by this, herein, the intrinsic contributions of the 2'-OH group in the nucleoside have been systematically investigated by directly solving the single-crystal structures of 2'-deoxy-2-aminoadenosine (1), 2-aminoadenosine (2), and 2-aminoarabinofuranosyladenine (3) in water. The 2'-OH group not only influenced the conformation and base-pair pattern of the single-nucleoside molecule, but also played a fundamental role in the entire supramolecular structure. Interestingly, compound 1, which did not contain the 2'-OH group, displayed strong hydration, whereas 2 and 3 (with the 2'-OH group in the opposite direction) exhibited no hydration, which was completely different from that observed in nucleic acids. Meanwhile, compound 1 trapped water molecules to form unique trihydrol moieties, which further served as the backbone to construct the simplest double-chain DNA-like structures. To this end, to investigate the effect of the biological environment on these unique structures, the solvent was changed from water to phosphate-buffered saline (PBS). Surprisingly, such a subtle adjustment led to entirely different superstructures, consisting of 2D lamellar structures in water and 3D porous structures in PBS. These large morphological differences could be attributed to delicate ion hydration, which was also confirmed through variable-temperature X-ray analysis, SEM, and intermolecular interaction energy calculations. In summary, this study comprehensively investigated the intrinsic contributions of 2'-hydroxyl to the hydration of nucleosides at the monomeric level; this is helpful to further understand the differences in DNA/RNA and the impact of their surrounding environment.

NMR Observation of Intermolecular Hydrogen Bonds between Protein Tyrosine Side-Chain OH and DNA Phosphate Groups

Hydrogen bonds between protein side-chain hydroxyl (OH) and phosphate groups are one of the most common types of intermolecular hydrogen bonds in protein-DNA/RNA complexes. Using NMR spectroscopy, we identified and characterized the hydrogen bonds between tyrosine side-chain OH and DNA phosphate groups in a protein-DNA complex. These OH groups exhibited relatively slow hydrogen-exchange rates and sizable scalar couplings between hydroxyl 1H and DNA phosphate 31P nuclei across the hydrogen bonds. Information about intermolecular hydrogen bonds facilitates investigations of the DNA/RNA recognition by the protein.

Effect of the sulfur atom on S2 and S4 positions of the uracil ring in different DNA:RNA hybrid microhelixes with three nucleotide base pairs

The effect of the sulphur atom on the uracil ring was analyzed in different DNA:RNA microhelixes with three nucleotide base-pairs, including uridine, 2-thiouridine, 4-thiouridine, 2,4-dithiouridine, cytidine, adenosine and guanosine. Distinct backbone and helical parameters were optimized at different density functional (DFT) levels. The Watson-Crick pair with 2-thiouridine appears weaker than with uridine, but its interaction with water molecules appears easier. Two types of microhelixes were found, depending on the H-bond of H2' hydroxyl atom: A-type appears with the ribose ring in ³ E-envelope C_3' -endo, and B-type in ² E-envelope C_2' -endo. B-type is less common but it is more stable and with higher dipole-moment. The sulphur atoms significantly increase the dipole-moment of the microhelix, as well as the rise and propeller twist parameters. Simulations with four Na atoms H-bonded to the phosphate groups, and further hydration with explicit water molecules were carried out. A re-definition of the numerical value calculation of several base-pair and base-stacking parameters is suggested.

Direct Observation of Carbohydrate Hydroxyl Protons in Hydrogen Bonds with a Protein

Hydroxyl proton resonances of uniformly ¹³C-labeled Manα(1-2)Manα(1-2)ManαOMe (Man₃) bound to cyanovirin-N (CV-N) were detected at ambient temperature in aqueous solution by NMR spectroscopy. The directions of the hydroxyl groups were determined on the basis of NOEs, and a previously unknown hydrogen-bonding network between Man₃ and CV-N was discovered. This is the first report on detecting hydroxyl protons of a protein-bound carbohydrate in aqueous solution by NMR. Approaches such as those presented here may open the door for accurately determining intermolecular hydrogen bonds in carbohydrate-protein complexes.

Chemical shifts in nucleic acids studied by density functional theory calculations and comparison with experiment

NMR chemical shifts are highly sensitive probes of local molecular conformation and environment and form an important source of structural information. In this study, the relationship between the NMR chemical shifts of nucleic acids and the glycosidic torsion angle, χ, has been investigated for the two commonly occurring sugar conformations. We have calculated by means of DFT the chemical shifts of all atoms in the eight DNA and RNA mono-nucleosides as a function of these two variables. From the DFT calculations, structures and potential energy surfaces were determined by using constrained geometry optimizations at the BP86/TZ2P level of theory. The NMR parameters were subsequently calculated by single-point calculations at the SAOP/TZ2P level of theory. Comparison of the (1)H and (13)C NMR shifts calculated for the mono-nucleosides with the shifts determined by NMR spectroscopy for nucleic acids demonstrates that the theoretical shifts are valuable for the characterization of nucleic acid conformation. For example, a clear distinction can be made between χ angles in the anti and syn domains. Furthermore, a quantitative determination of the χ angle in the syn domain is possible, in particular when (13)C and (1)H chemical shift data are combined. The approximate linear dependence of the C1' shift on the χ angle in the anti domain provides a good estimate of the angle in this region. It is also possible to derive the sugar conformation from the chemical shift information. The DFT calculations reported herein were performed on mono-nucleosides, but examples are also provided to estimate intramolecularly induced shifts as a result of hydrogen bonding, polarization effects, or ring-current effects.

Probing RNA dynamics via longitudinal exchange and CPMG relaxation dispersion NMR spectroscopy using a sensitive 13C-methyl label

The refolding kinetics of bistable RNA sequences were studied in unperturbed equilibrium via (13)C exchange NMR spectroscopy. For this purpose a straightforward labeling technique was elaborated using a 2'-(13)C-methoxy uridine modification, which was prepared by a two-step synthesis and introduced into RNA using standard protocols. Using (13)C longitudinal exchange NMR spectroscopy the refolding kinetics of a 20 nt bistable RNA were characterized at temperatures between 298 and 310K, yielding the enthalpy and entropy differences between the conformers at equilibrium and the activation energy of the refolding process. The kinetics of a more stable 32 nt bistable RNA could be analyzed by the same approach at elevated temperatures, i.e. at 314 and 316 K. Finally, the dynamics of a multi-stable RNA able to fold into two hairpin- and a pseudo-knotted conformation was studied by (13)C relaxation dispersion NMR spectroscopy.

Effect of pseudouridylation on the structure and activity of the catalytically essential P6.1 hairpin in human telomerase RNA

Telomerase extends the 3'-ends of linear chromosomes by adding conserved telomeric DNA repeats and is essential for cell proliferation and genomic stability. Telomerases from all organisms contain a telomerase reverse transcriptase and a telomerase RNA (TER), which together provide the minimal functional elements for catalytic activity in vitro. The RNA component of many functional ribonucleoproteins contains modified nucleotides, including conserved pseudouridines (Ψs) that can have subtle effects on structure and activity. We have identified potential Ψ modification sites in human TER. Two of the predicted Ψs are located in the loop of the essential P6.1 hairpin from the CR4-CR5 domain that is critical for telomerase catalytic activity. We investigated the effect of P6.1 pseudouridylation on its solution NMR structure, thermodynamic stability of folding and telomerase activation in vitro. The pseudouridylated P6.1 has a significantly different loop structure and increase in stability compared to the unmodified P6.1. The extent of loop nucleotide interaction with adjacent residues more closely parallels the extent of loop nucleotide evolutionary sequence conservation in the Ψ-modified P6.1 structure. Pseudouridine-modification of P6.1 slightly attenuates telomerase activity but slightly increases processivity in vitro. Our results suggest that Ψs could have a subtle influence on human telomerase activity via impact on TER-TERT or TER-TER interactions.

High-resolution NMR structure of an RNA model system: the 14-mer cUUCGg tetraloop hairpin RNA

We present a high-resolution nuclear magnetic resonance (NMR) solution structure of a 14-mer RNA hairpin capped by cUUCGg tetraloop. This short and very stable RNA presents an important model system for the study of RNA structure and dynamics using NMR spectroscopy, molecular dynamics (MD) simulations and RNA force-field development. The extraordinary high precision of the structure (root mean square deviation of 0.3 A) could be achieved by measuring and incorporating all currently accessible NMR parameters, including distances derived from nuclear Overhauser effect (NOE) intensities, torsion-angle dependent homonuclear and heteronuclear scalar coupling constants, projection-angle-dependent cross-correlated relaxation rates and residual dipolar couplings. The structure calculations were performed with the program CNS using the ARIA setup and protocols. The structure quality was further improved by a final refinement in explicit water using OPLS force field parameters for non-bonded interactions and charges. In addition, the 2'-hydroxyl groups have been assigned and their conformation has been analyzed based on NOE contacts. The structure currently defines a benchmark for the precision and accuracy amenable to RNA structure determination by NMR spectroscopy. Here, we discuss the impact of various NMR restraints on structure quality and discuss in detail the dynamics of this system as previously determined.

NMR assignments of HIV-2 TAR RNA

We report nearly complete assignment for all (1)H, (13)C, (31)P, and (15)N resonances in the 30-nucleotide stem-loop HIV-2 TAR RNA located at the 5' end of all viral mRNAs.

Liquid-crystal NMR structure of HIV TAR RNA bound to its SELEX RNA aptamer reveals the origins of the high stability of the complex

Transactivation-response element (TAR) is a stable stem-loop structure of HIV RNA, which plays a crucial role during the life cycle of the virus. The apical loop of TAR acts as a binding site for essential cellular cofactors required for the replication of HIV. High-affinity aptamers directed against the apical loop of TAR have been identified by the SELEX approach. The RNA aptamers with the highest affinity for TAR fold as hairpins and form kissing complexes with the targeted RNA through loop-loop interactions. The aptamers with the strongest binding properties all possess a GA base pair combination at the loop-closing position. Using liquid-crystal NMR methodology, we have obtained a structural model in solution of a TAR-aptamer kissing complex with an unprecedented accuracy. This high-resolution structure reveals that the GA base pair is unilaterally shifted toward the 5' strand and is stabilized by a network of intersugar hydrogen bonds. This specific conformation of the GA base pair allows for the formation of two supplementary stable base-pair interactions. By systematic permutations of the loop-closing base pair, we establish that the identified atomic interactions, which form the basis for the high stability of the complex, are maintained in several other kissing complexes. This study rationalizes the stabilizing role of the loop-closing GA base pairs in kissing complexes and may help the development or improvement of drugs against RNA loops of viruses or pathogens as well as the conception of biochemical tools targeting RNA hairpins involved in important biological functions.

RNA structure determination by NMR

This chapter reviews the methodologies for RNA structure determination by liquid-state nuclear magnetic resonance (NMR). The routine production of milligram quantities of isotopically labeled RNA remains critical to the success of NMR-based structure studies. The standard method for the preparation of isotopically labeled RNA for structural studies in solution is in vitro transcription from DNA oligonucleotide templates using T7 RNA polymerase and unlabeled or isotopically labeled nucleotide triphosphates (NTPs). The purification of the desired RNA can be performed by either denaturing polyacrylamide gel electrophoresis (PAGE) or anion-exchange chromatography. Our basic strategy for studying RNA in solution by NMR is outlined. The topics covered include RNA resonance assignment, restraint collection, and the structure calculation process. Selected examples of NMR spectra are given for a correctly folded 30 nucleotide-containing RNA.

2'-hydroxyl proton positions in helical RNA from simultaneously measured heteronuclear scalar couplings and NOEs

The 2'-hydroxyl group in RNA plays an important structural role; it defines hydration in the minor groove, impacts thermodynamic stability of RNA, and often participates in RNA catalysis. To better study this important functional group in RNA, we describe a constant-time HMQC-IPAP-NOESY 3D NMR experiment. It simultaneously yields highly resolved 13C-separated NOEs from ribose protons to the 2'OH proton, as well as E.COSY-type measurement of JC-OH couplings, thereby permitting a quantitative study of the orientation of the 2'OH proton in RNA. The observed NOE patterns indicate that the 2'OH bonds in A-form helical RNA are primarily oriented toward the base domain, as further supported by small (1.3 +/- 0.7 Hz) 3JC1'-2'OH and relatively large (4.2 +/- 0.7 Hz) 3JC3'-2'OH couplings. The constant-time HMQC-IPAP-NOESY is suitable for measurement of interactions of rapidly exchanging protons in proteins and nucleic acids.

Stacking effects on local structure in RNA: changes in the structure of tandem GA pairs when flanking GC pairs are replaced by isoG-isoC pairs

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.

NMR studies of the conformational effect of single and double 3'-S-phosphorothiolate substitutions within deoxythymidine trinucleotides

NMR spectroscopy has been used to investigate the conformational effects of single and two consecutive 3'-S-phosphorothiolate modifications within a deoxythymidine trinucleotide. The presence of a single 3'-phosphorothioate modification shifts the conformation of the sugar ring it is attached to, from a mainly south to north pucker; this effect is also transmitted to the 3'-neighbour deoxyribose. This transmission is thought to be caused by favourable stacking of the heterocyclic bases. Similar observations have been made previously by this group. When two adjacent modifications are present, the conformations of the attached deoxyribose rings are again shifted almost completely to the north, however, there is no transmission to the 3' deoxyribose ring. Base proton chemical shift analysis and molecular modelling have been used to aid elucidation of the origin of this feature. The observation for the dimodified sequence is consistent with our previously reported results for a related system in which spaced modifications are more thermodynamically stable than consecutive ones.

Dissecting non-canonical interactions in frameshift-stimulating mRNA pseudoknots

A variety of powerful NMR experiments have been introduced over the last few years that allow for the direct identification of different combinations of donor and acceptor atoms involved in hydrogen bonds in biomolecules. This ability to directly observe tertiary structural hydrogen bonds in solution tremendously facilitates structural studies of nucleic acids. We show here that an adiabatic HNN-COSY pulse scheme permits observation and measurement of J(N,N) couplings for nitrogen sites that are separated by up to 140 ppm in a single experiment at a proton resonance frequency of 500 MHz. Crucial hydrogen bond acceptor sites in nucleic acids, such as cytidine N3 nitrogens, can be unambiguously identified even in the absence of detectable H41 and H42 amino protons using a novel triple-resonance two-dimensional experiment, denoted H5(C5C4)N3. The unambiguous identification of amino nitrogen donor and aromatic nitrogen acceptor sites associated with both major groove as well as minor groove triple base pairs reveal the details of hydrogen bonding networks that stabilize the complex architecture of frameshift-stimulating mRNA pseudoknots. Another key tertiary interaction involving a 2'-OH hydroxyl proton that donates a hydrogen bond to an aromatic nitrogen acceptor in a cis Watson-Crick/sugar edge interaction can also be directly detected using a quantitative J(H,N) 1H,15N-HSQC experiment.

Influence of the 2'-hydroxyl group conformation on the stability of A-form helices in RNA

The 2'-hydroxyl group plays fundamental roles in both the structure and the function of RNA, and is the major determinant of the conformational and thermodynamic differences between RNA and DNA. Here, we report a conformational analysis of 2'-OH groups of the HIV-2 TAR RNA by means of NMR scalar coupling measurements in solution. Our analysis supports the existence of a network of water molecules spanning the minor groove of an RNA A-form helix, as has been suggested on the basis of a high-resolution X-ray study of an RNA duplex. The 2'-OH protons of the lower stem nucleotides of the TAR RNA project either towards the O3' or towards the base, where the 2'-OH group can favorably participate in H-bonding interactions with a water molecule situated in the nucleotide base plane. We observe that the k(ex) rate of the 2'-OH proton with the bulk solvent anti-correlates with the base-pair stability, confirming the involvement of the 2'-OH group in a collective network of H-bonds, which requires the presence of canonical helical secondary structure. The methodology and conformational analysis presented here are broadly applicable and facilitate future studies aimed to correlate the conformation of the 2'-OH group with both the structure and the function of RNA and RNA-ligand complexes.