Conformational analysis of the B and Z forms of the d(m5C-G)3 and d(br5C-G)3 hexamers in solution. A 300-MHz and 500-MHz NMR study

A general strategy exploiting m5C duplex-remodelling effect for selective detection of RNA and DNA m5C methyltransferase activity in cells

RNA:5-methylcytosine (m⁵C) methyltransferases are currently the focus of intense research following a series of high-profile reports documenting their physiological links to several diseases. However, no methods exist which permit the specific analysis of RNA:m⁵C methyltransferases in cells. Herein, we described how a combination of biophysical studies led us to identify distinct duplex-remodelling effects of m⁵C on RNA and DNA duplexes. Specifically, m⁵C induces a C3′-endo to C2′-endo sugar-pucker switch in CpG RNA duplex but triggers a B-to-Z transformation in CpG DNA duplex. Inspired by these different ‘structural signatures’, we developed a m⁵C-sensitive probe which fluoresces spontaneously in response to m⁵C-induced sugar-pucker switch, hence useful for sensing RNA:m⁵C methyltransferase activity. Through the use of this probe, we achieved real-time imaging and flow cytometry analysis of NOP2/Sun RNA methyltransferase 2 (NSUN2) activity in HeLa cells. We further applied the probe to the cell-based screening of NSUN2 inhibitors. The developed strategy could also be adapted for the detection of DNA:m⁵C methyltransferases. This was demonstrated by the development of DNA m⁵C-probe which permits the screening of DNA methyltransferase 3A inhibitors. To our knowledge, this study represents not only the first examples of m⁵C-responsive probes, but also a new strategy for discriminating RNA and DNA m⁵C methyltransferase activity in cells.

5-Methylcytosine containing CG decamer as Z-DNA embedded sequence for a potential Z-DNA binding protein probe

Attempting to elucidate biological significance of the left-handed Z-DNA is a research challenge due to Z-DNA potential role in many diseases. Discovery of Z-DNA binding proteins has ignited the interest in search for Z-DNA functions. Biosensor with Z-DNA forming probe can be useful to study the interaction between Z-DNA conformation and Z-DNA binding proteins. In this study, 5-methylcytosine (^mC) containing CG decamers were characterized for their suitability to form Z-DNA and to be used in Z-DNA forming probe. The 5'-thiol oligonucleotide embedded with 5'-^mCG^mCG^mCG^mCG^m CG-3' was designed and developed as a potential Z-DNA forming probe for Z-DNA binding protein screening.

Conformational preferences of modified uridines: comparison of AMBER derived force fields

The widespread occurrence of modified residues in RNA sequences necessitates development of accurate parameters for these modifications for reliable modeling of RNA structure and dynamics. A comprehensive set of parameters for the 107 naturally occurring RNA modifications was proposed by Aduri et al. (J. Chem. Theory Comput. 2007, 3, 1464-1475) for the AMBER FF99 force field. In this work, we tested these parameters on a set of modified uridine residues, namely, dihydrouridine, 2-thiouridine, 4-thiouridine, pseudouridine, and uridine-5-oxyacetic acid, by performing molecular dynamics and replica exchange molecular dynamics simulations of these nucleosides. Although our simulations using the FF99 force field did not, in general, reproduce the experimentally observed conformational characteristics well, combination of the parameter set with recent revisions of the FF99 force field for RNA showed noticeable improvement for some of the nucleosides.

Solution structure of the TnAn DNA duplex GCCGTTAACGCG containing the HpaI restriction site

The solution structure of the self-complementary DNA duplex [d(GCCGTTAACGGC)]2, which contains the HpaI restriction site GTTAAC, has been elucidated by two-dimensional NMR, distance geometry (DG), and NOE back-calculation methods. Initial distance constraints were determined by polynomial fitting the two-spin initial NOE rates; backbone constraints from NOE and J-coupling observations (Kim et al., 1992) were included. RMSDs between initial-distance-refined structures derived from random-embedded DG, A-DNA, and B-DNA starting structures were all in the range 0.5-1.0 A, indicating good convergence properties of the algorithm, regardless of the starting structure. A semiautomatic back-calculation refinement procedure was developed and used to generate more refined structures for which the BKCALC-simulated NOE volumes matched the experimental data. The six final structures refined from various starting structures exhibit very good agreement with the experimental data (R values = 0.18) and converge well to within 0.8-A RMSD differences for the central 8 base pairs. The torsion and pseudorotation phase angles were found to be well determined by the data, and the local helical parameters for each base step converged quite well. The final structures show that the central T6-A7 step is somewhat underwound (twist angle ca. 29 degrees), with a large negative cup and a normal (wide) minor groove width, while the T5-T6 and A7-A8 steps have a partially narrowed minor groove.

The importance of the specific Z-DNA structure and polyamines in carcinogenesis: fact or fiction

In this work some aspects of carcinogenesis are given. The importance of the emergence of Z or H DNA structure in the gene, or in the flanking gene sequences for the gene deletion and unusual gene recombination, is discussed. Some considerations on the role of selective pressure (of polyamines, of Mg2+, of the various levels of topoisomerase II, and of ATP) in the process of oncogene amplification, are given too.

Deoxyribose conformation in [d(GTATATAC)]2: evaluation of sugar pucker by simulation of double-quantum-filtered COSY cross-peaks

Exchangeable and nonexchangeable proton and phosphorus resonances (11.75 T) of [d(GTATATAC)]2 in aqueous solution were assigned by using proton two-dimensional nuclear Overhauser effect (2D NOE) spectra, homonuclear proton double-quantum-filtered COSY (2QF-COSY) spectra, proton spin-lattice relaxation time measurements, and 31P1H heteronuclear shift correlation spectra. Due to the large line widths, it was not possible to directly extract vicinal proton coupling constant values from any spectrum including ECOSY or 2QF-COSY. However, comparison of quantitative 2QF-COSY spectral simulations with experimental spectra enabled elucidation of coupling constants. The scope and limitations of this approach were explored by computation and by use of experimental data. It was found that proton line widths exhibit some variability from one residue to the next as well as from one proton to the next within a residue and the exact line width is critical to accurate evaluation of coupling constants. Experimental 2QF-COSY spectra were not consistent with a rigid deoxyribose conformation for any of the nucleotide residues. A classical two-state model, with rapid jumps between C2'-endo (pseudorotation angle P = 162 degrees) and C3'-endo (P = 9 degrees) conformations, was able to account for the spectral characteristics of terminal residue sugars: 60% C2'-endo and 40% C3'-endo. However, the 2QF-COSY cross-peaks from the -TATATA- core could be simulated only if the classical two-state model was altered such that the dominant conformer had a pseudorotation angle at 144 degrees instead of 162 degrees. In this case, the major conformer amounted to 80-85%. Alternatively, the spectral data were consistent with a three-state model in which C2'-endo and C3'-endo conformations had the largest and smallest populations, respectively, but a third conformer corresponding to C1'-exo (P = 126 degrees) was present, consistent with recent molecular dynamics calculations. This alternative yielded populations of 50% (P = 162 degrees), 35% (P = 126 degrees), and 15% (P = 9 degrees) for the -TATATA- sugars. The spectral results indicate little variation of sugar pucker between T and A. Small differences in cross-peak component intensities and characteristic spectral distortions, however, do suggest some unquantified variation. 31P1H heteronuclear chemical shift correlation spectra manifested alternating chemical shifts and coupling constants suggestive of phosphodiester backbone conformational differences between TA and AT junctions.

1H NMR assignments and conformational analysis of the oligoribonucleotides CA, CAU, CAUG, ACAUG, and UCAUG: observation of pyrimidine H5-H1' long-range scalar couplings

Conformational analysis and 1H NMR spectral assignments have been carried out using COSY and RELAY methods for a series of related oligoribonucleotides including two pentamers with 5'-dangling bases. Intraresidue long-range five bond scalar coupling was observed between pyrimidine H5 and H1' protons in the COSY-45 spectra and this feature was useful for both assignment purposes and conformational analysis. The ribose ring conformations were predominantly C3'-endo with the C2'-endo population increasing at the 3'-terminus. The 5'-dangling bases were not stacked efficiently, exhibiting lower % C3'-endo values than their 3'-nearest neighbors. Backbone torsion angle population. beta t, gamma +, epsilon t, were determined using 1H-1H, 1H-31P, and 13C-31P coupling constants. From beta t and gamma + populations the U3-G4 step in CAUG was found to be less efficiently stacked than the C1-A2 and A2-U3 steps. This observation in solution is consistent with the fiber diffraction A-RNA model (S. Arnott, D.W.L. Hukins, S.D. Dover, W. Fuller and A.R. Hodgson, J. Mol. Biol. 81, 107-122, 1973) which also predicts poor stacking in a U-G dinucleotide. The epsilon t populations were greater than 65% for all C3'-O3' bonds and consistent with a right-handed A-RNA helix.

Z-RNA: the solution NMR structure of r(CGCGCG)

Left-handed double-helical Z-RNA has been studied using the ribohexanucleotide pentaphosphate r(CpGpCpGpCpG). One-dimensional and two-dimensional proton nmr experiments were used to probe the structural details of the left-handed helix in concentrated sodium perchlorate solution. In 1M NaClO4 the RNA adopts the normal A-form double helix, and in 6M NaClO4 it is nearly all in the Z form. In 4M NaClO4 it exists as nearly equal parts of A form and Z form. Resonances corresponding to both A and Z form appear in the nmr spectrum, indicating that the duplex exchanges slowly between forms. Spin-spin coupling constants between protons in the ribose rings were used to determine the sugar-pucker conformations of the individual nucleotides. Quantitative nuclear Overhauser experiments were used to determine proton-proton distances within the nucleoside, and from these distances values for the glycosidic torsion angle were determined. The results show that the cytidines adopt C2'-endo sugar puckers (S type) with pseudo-rotation phase values (P) of approximately 165 degrees. The bases are in the anti conformation, with chi values of approximately -140 degrees. The internal guanosines adopt C3'-endo sugar puckers (N type) with P approximately 18 degrees, while the 3'-terminal guanosine ribose exists in an equilibrium between S- and N-type conformations. All three guanosine bases adopt the syn conformation, with chi approximately 70 degrees. The results indicate that the solution structure of Z-RNA is very similar to that of Z-DNA.

A Monte Carlo simulation study of the influence of internal motions on the molecular conformation deduced from two-dimensional NMR experiments

Monte Carlo methods have been used to simulate internal motions of aromatic protons of an oligonucleotide at the nanosecond time scale. Each proton is allowed to fluctuate about its equilibrium position. The longitudinal cross-relaxation rates of such a system of spins have been determined by computing the appropriate correlation functions. Then the interproton distances have been deduced according to the procedure generally used in two-dimensional nmr techniques (nuclear Overhauser effect spectroscopy--NOESY) and compared to the true values. The influence of the amplitude A and of the internal rotational diffusion constant Dint characterizing the dynamics of the system has been checked for in-phase and for uncorrelated motions. It is shown that for the investigated models the distances deduced from NOESY experiments may be under- or overestimated, depending strongly on the values of A and Dint. Furthermore, the cross-relaxation rate of a couple of protons is very sensitive to the correlation level of the motions of both protons.

Non-contiguous regions of Z-DNA in a DNA dodecamer

The conformation of the self-complimentary DNA dodecamer d(br5CGbr5CGAATTbr5CGbr5CG) has been investigated in a variety of salt and solvent conditions by one and two-dimensional 1H NMR. In low salt aqueous solutions, the molecule forms a regular B-DNA structure similar to the unmodified dodecamer. However, in aqueous solution containing high salt concentration and methanol, the dodecamer adopts a structure in which the br5CGbr5CG ends of the molecule are in a Z-DNA like conformation and the AATT region is neither standard B-DNA nor Z-DNA. The implications of these results for the structure of junctions between B and Z-DNA and the sequence specificity of Z-DNA are discussed.

Conformational studies of d-(AAAAATTTTT)2 using constraints from nuclear overhauser effects and from quantitative analysis of the cross-peak fine structures in two-dimensional 1H nuclear magnetic resonance spectra

The conformation at the dA-dT junction in d-(AAAAATTTTT)2 was investigated by using a variety of phase-sensitive two-dimensional nuclear magnetic resonance experiments at 500 MHz for detailed studies of the deoxyribose ring puckers. Conformational constraints were collected from two-dimensional nuclear Overhauser enhancement spectra recorded with short mixing times and from quantitative simulations of the cross-peaks in two-dimensional correlated spectra. Overall, the decamer duplex adopts a conformation of the B-DNA type, and for dA4 and dA5 the pseudorotation phase angle P is in the standard range 150-180 degrees. The deoxyribose puckers for the other nucleotides deviate significantly from the standard B-DNA structure. Spectrum simulations assuming either static deviations from standard B-DNA or a simple two-state dynamic equilibrium between the C2'-endo and C3'-endo forms of the deoxyribose were used to analyze the experimental data. It was thus found that the ring pucker for dT6 deviates from the regular C2'-endo form of B-DNA by a static distortion, with the pseudorotation phase angle P in the range 100-130 degrees, and a similar value of P is indicated for dT7. For the peripheral base pairs dynamic distortions of the C2'-endo form of the deoxyribose were found. In agreement with recent papers on related duplexes containing (dA)n tracts, we observed prominent nuclear Overhauser effects between adenine-2H and deoxyribose-1'H, which could be largely due to pronounced propeller twisting as observed in the crystal structures of (dA)n-containing compounds.

An NMR study of polymorphous behaviour of the mismatched DNA octamer d(m5C-G-m5C-G-A-G-m5C-G) in solution. The B-duplex and hairpin forms

By means of one- and two-dimensional NMR spectroscopy the solution structures of the partly self-complementary octamer d(m5C-G-m5C-G-A-G-m5C-G) were investigated. It is shown that this DNA fragment, under conditions of high DNA concentration (8 mM DNA) and/or high ionic strength prefers to adopt a duplex structure. At low DNA concentration (0.4 mM DNA), the duplex exists in a 1:1 slow equilibrium with a monomeric hairpin form. Addition of salt destabilizes the hairpin structure in favour of the dimer. At high temperatures the hairpin form, as well as the dimer structure, exist in a fast equilibrium with the random-coil form. For the hairpin/random-coil equilibrium a Tm of 329 K and a delta H degree of -121 kJ.mol-1 were deduced. These thermodynamic parameters are independent of the DNA concentration, as is expected for a monomeric structure. For the dimer to coil transition a Tm of 359 K (1 M DNA) and a delta H degree of -285 kJ.mol duplex-1 were derived. The thermodynamic data of the hairpin-coil transition mutually agree with those recently reported for the hairpin to random coil equilibrium of the DNA octamer d(m5C-G-m5C-G-T-G-m5C-G) [Orbons, L. P. M., van der Marel, G. A., van Boom, J. H. & Altona, C. (1987) J. Biomol. Struct. Dyns. 4, 939-963]. It is demonstrated that the dimer structure exhibits B-DNA characteristics, as is witnessed by the NOESY experiments and the analysis of the proton-proton coupling data. It is shown that the base-pair formation of the G x A mismatches is anti-anti. A comparison of 1H and 31P chemical-shift data of the title compound with those of a well-characterized B-DNA structure reveals large differences in the dm5C(3)-dG(4)-dA(5) part of the mismatched dimer structure. These differences apparently indicate some major local structural changes due to the incorporation of the G x A mismatches. Under the most extreme conditions used (i.e. up to 3 M NaCl or 75% CH3OH in the presence of 10 mM MgCl2) no Z-DNA structure was observed. It is shown that the structural features of the hairpin form of the title compound mimic those of the hairpin structure of d(m5C-G-m5C-G-T-G-m5C-G). An energy-minimized model of the hairpin form is given.

Synthesis and proton-NMR studies of oligonucleotides containing an apurinic (AP) site

In order to elucidate the conformational properties of base-deleted oligodeoxyribonucleotides, the molecules d-CpS(pCpG)n (n = 1,2; S = sugar) were synthesized by the phosphotriester method and characterized by 1H-NMR spectroscopy. Complete assignment of all non-exchangeable proton resonances of both compounds was obtained by 1D- and 2D-NMR techniques. In combination with computer simulation, these spectra yielded proton-proton and proton-phosphorus coupling constants of high accuracy. These data provide valuable information about the sugar and the backbone conformation. It appears that d-Cp1Sp2Cp3G4 does not form a duplex under any of the conditions studied. On the contrary, the base-deleted hexamer d-Cp1Sp2Cp3Gp4Cp5G6 occurs as a right-handed' staggered' DNA duplex at 280 K: the core of this duplex is formed by the residues C(3)-G(6); two 'dangling' residues C(1) and S(2) are located at the two 5'-ends of the duplex. The assignment of the corresponding imino proton resonances for [d-CpS(pCpG)2]2 was based on their thermal behavior: the line broadening of these resonances was studied as a function of temperature. The chemical shift and the number of imino proton resonances accord well with the number and type of Watson-Crick base pairs which can be formed in the staggered duplex described above. Thermodynamic parameters of duplex formation were obtained from an analysis of the chemical shift versus temperature profiles of aromatic base and H-1' protons. It is suggested that the cytosine ring of C(1) stacks, at least part of the time, with the guanine ring on the nucleotide residue, G(6), situated in the complementary strand. The binding of Lys-Trp-Lys to [d-CpS(pCpG)2]2 as well as to [d-CpGpCpG]1 was investigated. It is concluded that the indole ring of the tryptophan residue probably stacks on top of the 3'-terminal guanine base of both duplexes, but not on the nucleic acid bases next to the apurinic (AP) site.

Two-dimensional 1H NMR study of d(br5C-G)3 in the Z-form. Self association and flexibility of the left-handed double helix

The Z conformation of the auto complementary hexanucleoside pentaphosphate d(br5C-G)3 in 1 M NaClO4 solution has been investigated by using 2D NMR techniques. NOESY experiments performed at different temperatures show that the oligonucleotide exhibits end-to-end associations at room temperature. The conformation of the hexanucleotide molecules is very similar to that found in the crystal which was described by Chevrier et al. (J. Mol. Biol., 1986, 188, 707-719) as a Z-I form. When the temperature is increased the aggregates are dissociated and a conformational change is observed which is interpreted as a Z-I in equilibrium Z-II transition.

Influence of the base sequence on the conformational behaviour of DNA polynucleotides in solution

NMR studies were carried out on samples of the non-self-complementary tetramers d(C-A-C-A), d(T-G-T-G), d(G-A-G-A) and d(T-C-T-C) and of 1:1 mixtures of the complementary tetramers d(C-A-C-A) X d(T-G-T-G) and d(G-A-G-A) X d(T-C-T-C) at two DNA concentrations and of the self-complementary octamers d(C-A-C-A-T-G-T-G) and d(G-A-G-A-T-C-T-C). Assignments, based upon one-dimensional NOE and homonuclear-decoupling and two-dimensional correlated and NOE spectroscopies are given of the resonances of most of the base and sugar protons. Chemical shift vs temperature profiles, constructed for all samples, yielded insight into the temperature- and concentration-dependent conformational behaviour of the compounds and were used to obtain thermodynamic parameters pertaining to the stacked-single-strand----random-coil and duplex----random-coil equilibria. Vicinal proton-proton couplings were analyzed in terms of the conformation of the deoxyribose rings in the single-stranded tetramers and duplexed octamers. The NOE patterns, chemical shift profiles, imino-proton resonances and coupling data revealed that the compounds adopt B-DNA-like structures. The ratio duplexed/stacked-single-strand/random coil depends upon external conditions as well as upon base sequence. The thermodynamic data indicate that: in terms of single-helical stacking, the R-R steps (Tm 321-328 K) appear more stable than the Y-R or R-Y steps (Tm 308-316 K) and the Y-Y steps score least (Tm 290-300 K), and the duplexes consisting of alternating, d(Y-R)n, strands are more stable, in terms of delta H degrees, compared to the d(R-R)n X d(Y-Y)n duplexes. The analyses of the couplings demonstrated that the sugars of the single-stranded tetramers and duplexed octamers occur as a blend of N- and S-type conformers, with a preference for the S-type (C2'-endo) sugar conformation: upon duplex formation, no significant shift in the N-type/S-type ratio was observed. The fraction S-type sugar conformation of a given residue, %S, in the stacked-single strands was found to depend upon the nature of its own base and that of the adjacent residues: sugars in an R-R stretch display high values of %S (90-100), whereas those in Y-Y stretches show relatively low values (approximately equal to 65).(ABSTRACT TRUNCATED AT 400 WORDS)

Conformational and model-building studies of the hairpin form of the mismatched DNA octamer d(m5C-G-m5C-G-T-G-m5C-G)

The hairpin form of the mismatched octamer d(m5C-G-m5C-G-T-G-m5C-G) was studied by means of NMR spectroscopy. In a companion study it is shown that the hairpin form of this DNA fragment consists of a structure with a stem of three Watson-Crick-type base pairs and a loop consisting of only two nucleotides. The non-exchangeable proton resonances were assigned by means of two-dimensional correlation spectroscopy and two-dimensional nuclear Overhauser effect spectroscopy. Proton-proton coupling constants were used for the conformational analysis of the deoxyribose ring and for some of the backbone torsion angles. From the two-dimensional NMR spectra and the coupling-constant analysis it is concluded that: (i) the stem of the hairpin exhibits B-DNA characteristics; (ii) the sugar rings are not conformationally pure, but display a certain amount of conformational flexibility; (iii) the stacking interaction in the stem of the hairpin is elongated from the 3'-side in a more or less regular fashion with the two loop nucleotides; (iv) at the 5'-side of the stem a stacking discontinuity occurs between the stem and the loop; (v) at the 5'-side of the stem the loop is closed by means of a sharp backbone turn which involves unusual gamma' and beta+ torsion angles in residue dG(6). The NMR results led to the construction of a hairpin-loop model which was energy-minimized by means of a molecular-mechanics program. The results clearly show that a DNA hairpin-loop structure in which the loop consists of only two nucleotides bridging the minor groove in a straightforward fashion, (i) causes no undue steric strain, and (ii) involves well-known conformational principles throughout the course of the backbone. The hairpin form of the title compound is compared with the hairpin form of d(A-T-C-C-T-A-T4-T-A-G-G-A-T), in which the central -T4- part forms a loop of four nucleotides. Both models display similarities as far as stacking interactions are concerned.

An NMR study of the polymorphous behavior of the mismatched DNA octamer d(m5C-G-m5C-G-T-G-m5C-G) in solution. The B, Z, and hairpin forms

The polymorphism exhibited by the mismatched octamer d(m5C-G-m5C-G-T-G-m5C-G), as a function of the temperature, DNA concentration and ionic strength, was investigated by means of NMR spectroscopy. It is shown that this partly self-complementary DNA fragment, under conditions of low DNA concentration (0.4 mM) and low ionic strength, exclusively prefers to adopt a monomeric hairpin form, which consists of a stem of three Watson-Crick-type base pairs and a loop of only two residues. This in striking contrast with earlier intimations in literature, which postulated that in oligonucleotides loop formations containing only two residues are sterically impossible. Moreover, the hairpin form displays an unusual stability in comparison with previously reported hairpins. A Tm of 332 K and a delta H degree of -130 kJ.mol-1 were calculated for the hairpin to random coil transition. At high DNA concentration (8 mM) and/or upon the addition of sodium chloride the hairpin form occurs in slow exchange with a B-DNA dimer structure (approximately 20% at 270 K, no added salt), which comprises two central GxT-mismatched base pairs with the bases as major tautomers. At higher ionic strength (greater than 100 mM NaCl), or upon the addition of methanol, a third species appears, which is in slow exchange with both the B dimer and the hairpin form. This third species could be identified with a Z DNA form, comprising two GxT mismatches with the bases as major tautomers, with the guanine bases syn and the cytosine and thymine bases anti.

NMR and CD studies on an oligonucleotide containing N4-methylcytosine

The hexamer d(CGm4CGCG) exists predominantly as a right handed B form helix at 20 degrees C in 150 mM NaCl, as shown by 2D NOE spectra. Under these conditions a minor species is also observed which corresponds to the single strand in slow exchange on a proton NMR time scale with the double strand. This exchange is unusually slow and separate resonances for the two species are seen up to 65 degrees C. At 50 degrees C the lifetime of the single strand species is 0.85 s. Under high salt conditions the hexamer is partly converted into the Z form, but the complete transition is only observed at 5M NaCl at -6 degrees C.

Influence of N6-methylation of residue A(5) on the conformational behaviour of d(C-C-G-A-A-T-T-C-G-G) in solution studied by 1H-NMR spectroscopy. 1. The duplex form

One- and two-dimensional NMR studies at 300 MHz and 500 MHz were carried out on the two oligonucleotides d(C-C-G-A-A-T-T-C-G-G) and d(C-C-G-A-m6A-T-T-C-G-G) in aqueous solution. NMR spectra were observed at 10 mM sample concentration over the temperature range 273-368 K. Assignments are given of the base, H1', H2', H2", H3' and of some H4' resonances, based upon a combination of two-dimensional correlation spectra (COSY) and two-dimensional nuclear Overhauser effect spectra (NOESY); imino-proton resonances were assigned with the aid of a two-dimensional NOE experiment. Chemical shift vs temperature profiles were constructed in order to gain insight into the influence of N6-methylation of residue A(5) on the temperature-dependent conformational behaviour of the decamer and to determine thermodynamic parameters for the duplex-to-coil transition. The NOESY spectra, the imino-proton spectra and the shift profiles of the two compounds, under conditions where each forms a B-DNA-type duplex, are very similar. This is taken to indicate that the influence of N6-methylation of residue A(5) on the local structure of the duplex must be small. However, the temperature dependence of the (non-)exchangeable proton resonances of the two compounds reveals that methylation slows down the duplex-single-strand exchange. Furthermore, a thermodynamic analysis of the two compounds indicates that N6-methylation slightly decreases the stability of the duplex relative to the monomeric forms (Tm is reduced from 332 K down to 325 K at 10 mM sample concentration). Proton-proton couplings were obtained by means of one-dimensional and two-dimensional NMR experiments and were used in a conformational analysis of the sugar ring of each residue of the two compounds in the duplex form. The analysis indicated that all sugar rings display conformational flexibility in the intact duplex: population S-type sugar conformation ranges from 70% to 100%. A more refined analysis of the sugar rings of the parent compound revealed a sequence-dependent variation of the sugar geometry. This variation does not follow well the trend predicted by the Calladine/Dickerson sigma 3-sum rule [Dickerson, R. E. (1983) J. Mol. Biol. 166, 419-441; Calladine, C. R. (1982) J. Mol. Biol. 161, 343-352]; moreover the actual variations appear to be smaller in solution than those expected on the basis of known X-ray structures.

Influence of N6-methylation of residue A(5) on the conformational behaviour of d(C-C-G-A-A-T-T-C-G-G) in solution studied by 1H-NMR spectroscopy. 2. The hairpin form

The solution conformations of the oligonucleotides d(C-C-G-A-A-T-T-C-G-G) and d(C-C-G-A-m6A-T-T-C-G-G) as a function of temperature and sample concentration were investigated by means of 1H-NMR spectroscopy. The NMR spectra revealed that, at certain combinations of temperature and low sample and salt concentration, both compounds exist as a B-DNA-type duplex slowly (on the 1H-NMR time scale) interconverting with a monomeric species. From chemical shift data and imino-proton spectra, it is concluded that the monomeric species consists of a mixture of a hairpin form in rapid equilibrium with the random-coil form. The double-helical stem of the hairpin is formed by the six terminal cytidine and guanine residues, whereas the four core residues, -A-(m6)A-T-T-, partake in the loop. Thermodynamic analysis of the chemical shift of the resonances of the monomeric species vs temperature profiles of the two decamers and mutual comparison of these profiles indicate the following: the influence of N6-methylation of residue A(5) upon the local structure of the hairpin must be small; methylation decreases the stability of the duplex relative to the monomeric species: the temperature at which the fraction duplex equals 0.5 was found to be 312 K for the parent compound and 305 K for the methylated decamer at 2 mM sample concentration; methylation does not significantly alter the stability of the hairpin form relative to the random coil form: the Tm of the hairp----n equilibrium random-coil equilibrium is 308 K for the parent compound and 306 K for the methylated decamer. A higher fraction hairpin-like structure for the N6-methylated compound is observed under identical conditions of temperature and sample concentration: at 300 K, 2 mM sample concentration, the fraction hairpin form is 0.12 for d(C-C-G-A-A-T-T-C-G-G) and 0.20 for d(C-C-G-A-m6A-T-T-C-G-G). This finding appears to be a consequence of the reduced stability of the methylated dimeric species relative to the monomeric species, and to depend upon the sodium-ion concentration: it becomes more pronounced under low-salt conditions.

Conformational analysis of the deoxyribofuranose ring in DNA by means of sums of proton-proton coupling constants: a graphical method

A graphical method is presented for the conformational analysis of the sugar ring in DNA fragments by means of proton-proton couplings. The coupling data required for this analysis consist of sums of couplings, which are referred to as sigma 1' (= J1'2' + J1'2''), sigma 2' (= J1'2' + J2'3' + J2'2''), sigma 2'' (= J1'2'' + J2''3' + J2'2'') and sigma 3' (= J2'3' + J2''3' + J3'4'). These sums of couplings correspond to the distance between the outer peaks of the H1', H2', H2'' and H3' [31P] resonances, respectively, (except for sigma 2' and sigma 2'' in the case of a small chemical shift difference between the H2' and H2'' resonances) and can often be obtained from 1H-NMR spectra via first-order measurement, obviating the necessity of a computer-assisted simulation of the fine structure of these resonances. Two different types of graphs for the interpretation of the coupling data are discussed: the first type of graph serves to probe as to whether or not the sugar ring occurs as a single conformer, and if so to analyze the coupling data in terms of the geometry of this sugar ring. In cases where the sugar ring does not occur as a single conformer, but as a blend of N- and S-type sugar puckers, the second type of graph is used to analyze the coupling data in terms of the geometry and population of the most abundant form. It is shown that the latter type of analysis can be carried out on the basis of experimental values for merely sigma 1',sigma 2' and sigma 2'', without any assumptions or restrictions concerning a relation between the geometry of the N- and S-type conformer. In addition, the question is discussed as to how insight can be gained into the conformational purity of the sugar ring from the observed fine structure of the H1' resonance. Finally, a comparison is made between experimental coupling data reported for single-stranded and duplex DNA fragments and covalent RNA-DNA hybrids on the one hand and the predicted couplings and sums of couplings presented in this paper on the other hand.

The B and Z forms of the d(m5C-G)3 and d(br5C-G)3 hexamers in solution. A 300-MHz and 500-MHz two-dimensional NMR study

The non-exchangeable proton resonances of the hexadeoxynucleoside pentakisphosphates d(m5C-G)3 and d(br5C-G)3 in the B form as well as in the Z form were assigned by means of two-dimensional correlated spectroscopy and two-dimensional nuclear Overhauser enhancement spectroscopy. The complete proton NMR spectrum of the B form of the methylated compound was assigned in a pure 2H2O solution as well as in a 2H2O/C2H3O2H mixed solvent, containing 5 mM MgCl2. In the latter solvent the B form occurs in slow equilibrium (on the NMR time scale) with the Z form, the resonances of which also were fully assigned. The proton resonances of the B and Z forms of the brominated fragment were assigned in a 2H2O/C2H3O2H solution containing 5 mM MgCl2. A new and general method is described for the sequential assignment of the non-exchangeable proton resonances of oligonucleotides in the Z form.