Synthesis of a thymidine phosphoramidite labelled with 13C at C6: relaxation studies of the loop region in a 13C labelled DNA hairpin

Isotope Labels Combined with Solution NMR Spectroscopy Make Visible the Invisible Conformations of Small-to-Large RNAs

RNA is central to the proper function of cellular processes important for life on earth and implicated in various medical dysfunctions. Yet, RNA structural biology lags significantly behind that of proteins, limiting mechanistic understanding of RNA chemical biology. Fortunately, solution NMR spectroscopy can probe the structural dynamics of RNA in solution at atomic resolution, opening the door to their functional understanding. However, NMR analysis of RNA, with only four unique ribonucleotide building blocks, suffers from spectral crowding and broad linewidths, especially as RNAs grow in size. One effective strategy to overcome these challenges is to introduce NMR-active stable isotopes into RNA. However, traditional uniform labeling methods introduce scalar and dipolar couplings that complicate the implementation and analysis of NMR measurements. This challenge can be circumvented with selective isotope labeling. In this review, we outline the development of labeling technologies and their application to study biologically relevant RNAs and their complexes ranging in size from 5 to 300 kDa by NMR spectroscopy.

Synthesis and incorporation of 13C-labeled DNA building blocks to probe structural dynamics of DNA by NMR

We report the synthesis of atom-specifically 13C-modified building blocks that can be incorporated into DNA via solid phase synthesis to facilitate investigations on structural and dynamic features via NMR spectroscopy. In detail, 6-13C-modified pyrimidine and 8-13C purine DNA phosphoramidites were synthesized and incorporated into a polypurine tract DNA/RNA hybrid duplex to showcase the facile resonance assignment using site-specific labeling. We also addressed micro- to millisecond dynamics in the mini-cTAR DNA. This DNA is involved in the HIV replication cycle and our data points toward an exchange process in the lower stem of the hairpin that is up-regulated in the presence of the HIV-1 nucleocapsid protein 7. As another example, we picked a G-quadruplex that was earlier shown to exist in two folds. Using site-specific 8-13C-2'deoxyguanosine labeling we were able to verify the slow exchange between the two forms on the chemical shift time scale. In a real-time NMR experiment the re-equilibration of the fold distribution after a T-jump could be monitored yielding a rate of 0.012 min-1. Finally, we used 13C-ZZ-exchange spectroscopy to characterize the kinetics between two stacked X-conformers of a Holliday junction mimic. At 25°C, the refolding process was found to occur at a forward rate constant of 3.1 s-1 and with a backward rate constant of 10.6 s-1.

Preparation of selective and segmentally labeled single-stranded DNA for NMR by self-primed PCR and asymmetrical endonuclease double digestion

We demonstrate a new, efficient and easy-to-use method for enzymatic synthesis of (stereo-)specific and segmental (13)C/(15)N/(2)H isotope-labeled single-stranded DNA in amounts sufficient for NMR, based on the highly efficient self-primed PCR. To achieve this, new approaches are introduced and combined. (i) Asymmetric endonuclease double digestion of tandem-repeated PCR product. (ii) T4 DNA ligase mediated ligation of two ssDNA segments. (iii) In vitro dNTP synthesis, consisting of in vitro rNTP synthesis followed by enzymatic stereo-selective reduction of the C2' of the rNTP, and a one-pot add-up synthesis of dTTP from dUTP. The method is demonstrated on two ssDNAs: (i) a 36-nt three-way junction, selectively (13)C(9)/(15)N(3)/(2)H((1',2'',3',4',5',5''))-dC labeled and (ii) a 39-nt triple-repeat three-way junction, selectively (13)C(9)/(15)N(3)/(2)H((1',2'',3',4',5',5''))-dC and (13)C(9)/(15)N(2)/(2)H((1',2'',3',4',5',5''))-dT labeled in segment C20-C39. Their NMR spectra show the spectral simplification, while the stereo-selective (2)H-labeling in the deoxyribose of the dC-residues, straightforwardly provided assignment of their C1'-H2' and C2'-H2' resonances. The labeling protocols can be extended to larger ssDNA molecules and to more than two segments.

Dynamics of large elongated RNA by NMR carbon relaxation

We present an NMR strategy for characterizing picosecond-to-nanosecond internal motions in uniformly 13C/15N-labeled RNAs that combines measurements of R1, R1rho, and heteronuclear 13C{1H} NOEs for protonated base (C2, C5, C6, and C8) and sugar (C1') carbons with a domain elongation strategy for decoupling internal from overall motions and residual dipolar coupling (RDC) measurements for determining the average RNA global conformation and orientation of the principal axis of the axially symmetric rotational diffusion. TROSY-detected pulse sequences are presented for the accurate measurement of nucleobase carbon R1 and R1rho rates in large RNAs. The relaxation data is analyzed using a model free formalism which takes into account the very high anisotropy of overall rotational diffusion (Dratio approximately 4.7), asymmetry of the nucleobase CSAs and noncollinearity of C-C, C-H dipolar and CSA interactions under the assumption that all interaction tensors for a given carbon experience identical isotropic internal motions. The approach is demonstrated and validated on an elongated HIV-1 TAR RNA (taum approximately 18 ns) both in free form and bound to the ligand argininamide (ARG). Results show that, while ARG binding reduces the amplitude of collective helix motions and local mobility at the binding pocket, it leads to a drastic increase in the local mobility of "spacer" bulge residues linking the two helices which undergo virtually unrestricted internal motions (S2 approximately 0.2) in the ARG bound state. Our results establish the ability to quantitatively study the dynamics of RNAs which are significantly larger and more anisotropic than customarily studied by NMR carbon relaxation.

Structural differences in the NOE-derived structure of G-T mismatched DNA relative to normal DNA are correlated with differences in (13)C relaxation-based internal dynamics

Detailed description of the characteristics of mismatched DNA that are distinct from normal DNA is vital to the understanding of how mismatch repair proteins are able to recognize and repair these DNA lesions. To this end, we have used nuclear Overhauser effect spectroscopy (NOESY)-based distance restraints and (13)C relaxation measurements to solve the solution structures and measure some of the internal dynamics of the G-T mismatched DNA oligomer d(CCATGCGTGG)(2) (GT) and its parent DNA sequence d(CCACGCGTGG)(2) (GC). In GT, the mismatched G7 is structurally perturbed much more than the mismatched T4 relative to their corresponding bases in GC. The degree of G7 displacement differs from previous high-resolution structures of G-T mismatch-containing B-DNA, suggesting a dependence of G-T mismatch-induced structural perturbation on sequence context. The internal dynamics of GC and GT differ on multiple timescales. The mismatched G7 of GT contains spins that decrease significantly in order in GT compared to GC, while spins in C6, T8, and A3 have significantly higher order in GT compared to GC. Linear correlations between helical parameters of GC and GT and the order of C-1' and aromatic methine carbon atoms relate differences in internal dynamics to the structures quantitatively. The dynamic differences between the normal and mismatched DNA signify changes in local flexibility that may be exploited by the mismatch repair system to bind mismatched DNA preferentially while ignoring normal DNA.

Solid-phase synthesis of selectively labeled DNA: applications for multidimensional nuclear magnetic resonance spectroscopy

The solid-phase chemical synthesis method has a strong advantage over the enzymatic method for preparing selectively labeled DNA oligomers. Atom-specific and fully labeled 2'-deoxynucleosides are economically prepared with routinely available isotope precursors using this synthetic route. Special DNA oligomers prepared by advanced labeling techniques are needed for advanced NMR applications, and chemical synthesis is the method of choice to respond to such demands. As a summary of this chapter, two tables are given. Table I lists the labeled nucleosides reported to be available by chemical syntheses. Table II lists the NMR studies using labeled DNA oligomers that were prepared by chemical syntheses.

DNA duplex dynamics: NMR relaxation studies of a decamer with uniformly 13C-labeled purine nucleotides

Dynamics in a DNA decamer duplex, d(CATTTGCATC). d(GATGCAAATG), were investigated via a detailed 13C NMR relaxation study. Every 2'-deoxyadenosine and 2'-deoxyguanidine was chemically enriched with 15% 13C and 98% 15N isotopes. Six nuclear relaxation parameters [R(13Cz), R(1Hz), R(2(1)Hz13Cz), R(13Cx), R(2(1)Hz13Cx) and steady-state 13C¿1H¿ NOE] were measured at 600 MHz and three were measured at 500 MHz (1H frequency) for the CH spin systems of sugar 1', 3', and 4' as well as base 8 and 2 positions. A dependence of relaxation parameter values on chemical position was clearly observed; however, no sequence-specific variation was readily evident within our experimental error of approximately 5-10%, except for 3' and 5' termini. It was demonstrated that the random 15% 13C enrichment effectively suppressed both scalar and dipolar contributions of the neighboring carbons and protons on the relaxation parameters. To analyze dynamics via all observed relaxation parameters, full spectral density mapping (1992, J. W. Peng and G. Wagner, J. Magn. Reson. 98, 308) and the "model-free" approach (1982, Lipari and Szabo, J. Am. Chem. Soc. 104, 4546) were applied complementarily. A linear correlation between three spectral density values, J(omegaC), J(omegaH - omegaC), and J(omegaH + omegaC) was observed in plots containing all measured values, but not for the other spectral density terms including J(0). These linear correlations reflect the effect of overall motion and similar internal motions for each CH vector in the decamer. The correlations yielded two correlation times, 3-4 ns and 10-200 ps. One value, 3-4 ns, corresponds to the value of 3.3 ns obtained for the overall isotropic tumbling correlation time determined from analysis of 13C T1/T2 ratios. The possibility of overall anisotropic tumbling was examined, but statistical analysis showed no advantage over the assumption of simple isotropic tumbling. Lack of correlations entailing J(0) implies that a relatively slow chemical exchange contributes to yielding of effective Jeff(0) values. Based on spectral density mapping and the T1/T2 ratio analysis, three basic assumptions were initially employed (and subsequently justified) for the model-free calculation: isotropic overall tumbling, one internal motion, and the presence of chemical exchange terms. Except for terminal residues, the order parameter S2 and the corresponding fast internal motion correlation time were determined to be about 0.8 +/- 0.1 and 20 +/- 20 ps, respectively, for the various CH vectors. Only a few differences were observed between or within sugars and bases. The internal motion is very fast (ps-ns time scale) and its amplitude restricted; e.g., assuming a simple wobble-in-a-cone model, the internal motion is restricted to an angular amplitude of +/-22. 5 degrees for each of the 1', 3', 4', 2, and 8 positions in the purine nucleotides in the entire duplex.

Distance measurements in nucleic acids using windowless dipolar recoupling solid state NMR

A windowless, homonuclear dipolar recoupling pulse sequence (DRAWS) is described and a theoretical basis for describing its recoupling performance is developed using numerical techniques. It is demonstrated that DRAWS recouples weak dipolar interactions over a broad range of experimental and molecular conditions. We discuss two spectroscopic control experiments, which help to take into account effects due to insufficient proton decoupling, relaxation, and static dipolar couplings to nearby 13C spins at natural abundance. Finally DRAWS is used in combination with selective 13C labeling to measure 13C-13C distances in five doubly labeled DNA dodecamers, [d(CGCGAAT*T*CGCG)]2, which contain the binding site for the restriction enzyme EcoRI. The longest distance reported is 4.8 A. In most cases the distances agree well with those derived from X-ray crystallographic data, although small changes in hydration level can result in relatively large changes in internuclear distances.

Selectively 13C-enriched DNA: evidence from 13C1' relaxation rate measurements of an internal dynamics sequence effect in the lac operator

In order to study some internal dynamic processes of the lac operator sequence, the 13C-labeled duplex 5'd(C0G1C2T3C4A5C6A7A8T9T10).d(A10A9T8T7G6T5G4A3G2C1G0)3' was used. The spreading of both the 1H1' and 13C1' resonances brought about an excellent dispersion of the 1H1'-13C1' correlations. The spin-lattice relaxation parameters R(Cz), R(Cx,y) and R(Hz --> Cz) were measured for each residue of the two complementary strands, except for the 3'-terminal residues which were not labeled. Variation of the relaxation rates was found along the sequence. These data were analyzed in the context of the model-free formalism proposed by Lipari and Szabo [(1982) J. Am. Chem. Soc., 104, 4546-4570] and extended to three parameters by Clore et al. [(1990) Biochemistry, 29, 7387-7401; and (1990) J. Am. Chem. Soc., 112, 4989-4991]. A careful analysis using a least-squares program showed that our data must be interpreted in terms of a three-parameter spectral density function. With this approach, the global correlation time was found to be the same for each residue. All the C1'-H1' fragments exhibited both slow (tau s = 1.5 ns) and fast (tau f = 20 ps) restricted libration motions (Ss2 = 0.74 to 1.0 and Sf2 = 0.52 to 0.96). Relaxation processes were described as governed by the motion of the sugar relative to the base and in terms of bending of the whole duplex. The possible role played by the special structure of the AATT sequence is discussed. No evident correlation was found between the amplitude motions of the complementary residues. The 5'-terminal residues showed large internal motions (S2 = 0.5), which describe the fraying of the double helix. Global examination of the microdynamical parameters Sf2 and Ss2 along the nucleotide sequence showed that the adenine residues exhibit more restricted fast internal motions (Sf2 = 0.88 to 0.96) than the others, whereas the measured relaxation rates of the four nucleosides in solution were mainly of dipolar origin. Moreover, the fit of both R(Cz) and R(Hz --> Cz) experimental relaxation rates using an only global correlation time for all the residues, gave evidence of a supplementary relaxation pathway affecting R(Cx,y) for the purine residues in the (5' --> 3') G4A3 and A10A9T8T7 sequences. This relaxation process was analyzed in terms of exchange stemming from motions of the sugar around the glycosidic bond on the millisecond time scale. It should be pointed out that these residues gave evidence of close contacts with the protein in the complex with the lac operator [Boelens et al. (1987) J. Mol. Biol., 193, 213-216] and that these motions could be implied in the lac-operator-lac-repressor recognition process.

13C-labeled oligodeoxyribonucleotides: a solution study of a CCAAT-containing sequence at the nuclear factor I recognition site of human adenovirus

The solution behavior of the single-stranded CCAAT-containing octamer 1, d(AGCCAATA), that comprises part of the nuclear factor I (NF-I) recognition site at the origin of replication of human adenovirus has been studied by nmr spectroscopy at 500 and 600 MHz. Proton resonance assignments for 1 were aided by selective 13C enrichment at C1' of A1 or A5. High-resolution 13C-1H heteronuclear multiple-bond coherence spectra of the 13C-labeled oligomers permitted the selective detection of furanosyl ring protons within each labeled residue due to short- and long-range 13C-1H couplings to the enriched C1'. The resulting assignments provided firm starting points in the interpretation of double quantum filtered correlated spectra, yielding information supplemented by total correlated spectroscopy (TOCSY) and rotating frame nuclear Overhauser effect spectroscopic data to completely assign the 1H-nmr spectrum of 1 and extract 3JHH values for furanose conformational analysis. Several 13C-1H spin-coupling constants within the 13C-enriched A1 or A5 residues were measured from cross-peak shifts in TOCSY spectra, and their signs determined by inspection of the relative orientations of these shifts. 1H-1H and 13C-1H spin-couplings both indicate a preference (> 75%) for south (C2'-endo) conformations by the furanosyl rings of 1.

Preparation and heteronuclear 2D NMR spectroscopy of a DNA dodecamer containing a thymidine residue with a uniformly 13C-labeled deoxyribose ring

[13C5]-2-Deoxy-D-ribose, synthesized from [13C6]-D-glucose (98% 13C), was coupled with thymine to give [1',2',3',4',5'-13C5]-thymidine (T) in an 18% overall yield. The thymidine was converted to the 3'-phosphoramidite derivative and was then incorporated into a dodecamer 5'-d(CGCGAATTCGCG)-3' by solid-phase DNA synthesis. Preparation of 0.24 mumole of the labeled dodecamer, which is sufficient for a single NMR sample, consumed only 25 mg of glucose. By virtue of the 13C labels, all of the 1H-1H vicinal coupling constants in the sugar moieties were accurately determined by HCCH-E.COSY.

Structure of a small RNA hairpin

The hairpin stem-loop form of the RNA oligonucleotide rCGC(UUU)GCG has been studied by NMR spectroscopy. In 10 mM phosphate buffer this RNA molecule forms a unimolecular hairpin with a stem of three base pairs and a loop of three uridines, as judged by both NMR and UV absorbance melting behavior. Distance and torsion angle restraints were determined using homonuclear proton-proton and heteronuclear proton-phosphorus 2-D NMR. These values were used in restrained molecular dynamics to determine the structure of the hairpin. The stem has characteristics of A-form geometry, although distortion from A-form occurs in the 3'-side of the stem, presumably to aid in accommodating the small loop. The loop nucleotides adopt C2'-endo conformations. NOE's strongly suggest stacking of the uracils with the stem, especially the first uracil on the 5'-side of the loop. The reversal of the chain direction in the loop seems to occur between U5 and U6. Loop structures produced by molecular dynamics simulations had a wide range of conformations and did not show stacking of the uracils. A flexible loop with significant dynamics is consistent with all the data.

Analysis of melting transitions of the DNA hairpins formed from the oligomer sequences d[GGATAC(X)4GTATCC] (X = A, T, G, C)

Optical melting transitions of the short DNA hairpins formed from the self-complementary DNA oligomers d[GGATACX4GTATCC] where X = A, T, G, or C measured in 100 mM NaCl are presented. A significant dependence of the melting transitions on loop sequence is observed and transition temperatures, tm, of the hairpins vary from 58.3 degrees C for the T4 loop hairpin to 55.3 degrees C for the A4 loop. A nearest-neighbor sequence-dependent theoretical algorithm for calculating melting curves of DNA hairpins is presented and employed to analyze the experimental melting transitions. Experimental melting curves were fit by adjustment of a single theoretical parameter, Fend(n), the weighting function for a hairpin loop comprised of n single-strand bases. Empirically determined values of Fend(n) provide an evaluation of the free-energy of hairpin loop formation and stability. Effects of heterogeneous nearest-neighbor sequence interactions in the duplex stem on hairpin loop formation were investigated by evaluating Fend(n) in individual fitting procedures using two of the published sets of nearest-neighbor stacking interactions in DNA evaluated in 100 mM NaCl and given by Wartell and Benight, 1985. In all cases, evaluated values of Fend(n) were obtained that provided exact theoretical predictions of the experimental transitions. Results of the evaluations indicate: (1) Evaluated free-energies of hairpin loop formation are only slightly dependent on loop sequences examined. At the transition temperature, Tm, the free-energy of forming a loop of four bases is approximately equal for T4, G4, or C4 loops and varies from 3.9 to 4.8 kcal/mole depending on the set of nearest-neighbor interactions employed in the evaluations. This result suggests, in light of the observed differences in stability between the T4, G4, and C4 loop hairpins, that sequence-dependent interactions between base residues of the loop are most likely not the source of the enhanced stability of a T4 loop.(ABSTRACT TRUNCATED AT 400 WORDS)

Thermodynamic stability of the 5' dangling-ended DNA hairpins formed from sequences 5'-(XY)2GGATAC(T)4GTATCC-3', where X, Y = A, T, G, C

Expressions for the partition function Q (T) of DNA hairpins are presented. Calculations of Q (T), in conjunction with our previously reported numerically exact algorithm [T. M. Paner, M. Amaratunga, M. J. Doktycz, and A. S. Benight (1990) Biopolymers, 29, 1715-1734], yield a numerical method to evaluate the temperature dependence of the transition enthalpy, entropy, and free energy of a DNA hairpin directly from its optical melting curve. No prior assumptions that the short hairpins melt in a two-state manner are required. This method is then applied in a systematic manner to investigate the stability of the six basepair duplex stem 5'-GGATAC-3' having four-base dangling single-strand ends with the sequences (XY)2, where X, Y = A, T, G, C, on the 5' end and a T4 loop on the 3' end. Results show that all dangling ends of the sample set stabilize the hairpin against melting. Increases in transition temperatures as great as 4.0 degrees C above the blunt-ended control hairpin were observed. The hierarchy of the hairpin transition temperatures is dictated by the identity of the first base of the dangling end adjoining the duplex in the order: purine greater than T greater than C. Calculated melting curves of every hairpin were fit to experimental curves by adjustment of a single parameter in the numerically exact theoretical algorithm. Exact fits were obtained in all cases. Experimental melting curves were also calculated assuming a two-state melting process. Equally accurate fits of all dangling-ended hairpin melting curves were obtained with the two-state model calculation. This was not the case for the melting curve of the blunt-ended hairpin, indicating the presence of a four-base dangling-end drives hairpin melting to a two-state process. Q (T) was calculated as a function of temperature for each hairpin using the theoretical parameters that provided calculated curves in exact agreement with the experimentally obtained optical melting curves. From Q (T), the temperature dependence of the transition enthalpy delta H, entropy delta S, and free energy delta G were calculated for every hairpin providing a quantitative assessment of the effects of dangling ends on hairpin thermodynamics. Comparisons of our results are made with those of the Breslauer group [M. Senior, R. A. Jones, and K. J. Breslauer (1988) Biochemistry 27, 3879-3885] on the T2 5' dangling-ended d(GC)3 duplexes.(ABSTRACT TRUNCATED AT 400 WORDS)