ABSTRACT
The structure of the oligonucleotide d(A-G)10 in 0.6 M Na+, pH 7.0 has been investigated with UV resonance Raman (UVRR) spectroscopy. Variable wavelength excitation was used to distinguish the spectral contributions of dG and dA residues. Both classes of residues show UVRR hyperchromism with increasing temperature, reflecting unstacking of the bases. The dG residues melt relatively cooperatively with a Tm of ~42oC. Unstacking is non-cooperative for the dA residues, increasing linearly between 4 and 80oC. G-tetrads at low temperature are indicated by UVRR frequency shifts of modes associated with C6=O and C2-NH2 of the dG residues, and of vibrations involving N7, all sites of H-bonding. However, there are no indications of interbase H-bonds for the dA residues, showing they do not form H-bonded tetrads. Most of the bases are oriented anti about the glycosyl bond, but at 4oC a fraction of the residues are syn. These results, together with the findings by Shiber et al. [Shiber,M.C., Braswell,E.H., Klump,H. and Fresco,J.R. (1996) Nucleic Acids Res. 24, 5004-5012] that d(A-G)10 under comparable conditions has the molecular weight of a dimer, support a model in which two hairpins interact to form a helical structure with G-tetrads and intercalated dA residues.
Approximately 0.4% of the mammalian genome consists of homopurine.homopyrimidine sequences (2 ) which have been implicated in genomic processes such as transcriptional regulation and recombination (3 ,4 ). Concentration and pH strongly influence the structure of these sequences; e.g. at acidic pH homopurine. homopyrimidine tracts can form intramolecular triplexes (3 ,5 ). Similarly, dG-rich sequences, as in Oxytricha telomere DNA, are thought to form tetraplex structures stabilized by H-bonds that involve the same H-bonding sites on dG residues as in Hoogsteen and Watson-Crick base pairs (Fig. 1 , top) (6 ,7 ).
Telomere DNA (8 -10 ) is typically characterized by multiple repeats of two to five consecutive dG residues that occur at the ends of eukaryotic chromosomes. It has long been known that homopolymers of rG as well as 5'-GMP and other guanosine derivatives form helices with G-tetrads (11 -15 ). It has been shown that dG-rich oligomers, modeled after telomeric repeat sequences, adopt tetraplex structures which are thermodynamically and kinetically stable (8 ). This ability of telomere DNA to form tetraplex-type structures in vitro and the associated stability of these structures has elicited considerable interest because of their possible involvement in recombinational events (16 ,17 ) and in telomere function (18 ,19 ).
G-tetrad structures are characterized by a cyclic network of eight H-bonds (Fig. 1 , top), which confer a high degree of stability, the half-life for unfolding being on the order of hours (20 ). Additionally, NMR characterization has revealed very slow rates for imino proton exchange, and in certain cases imino protons have still been observed at 90oC (21 ). Interestingly, chemical modification and native gel electrophoresis experiments (22 ,23 ) indicate that telomeric sequences containing adenine form tetraplexes that are less stable than those formed by non-adenine-containing sequences. However, the presence of adenine in the telomeric repeat of humans and other vertebrates, d(TnAGGG) (9 ), raises questions regarding the effect of dA residues on tetraplex structure and stability and their possible role in telomere function.
2'-deoxyadenosine-5'-monophosphate-free acid (dAMP) and 2'-deoxyguanosine-5'-monophosphate-disodium salt (dGMP) (Peninsula Laboratories) were used without further purification. d(A-G)10 was synthesized, purified to homogeneity (24 ) and desalted (26 ). All concentrations are reported in total nucleotide residues.
d(A-G)10 (2 mM) was in 0.3 M Na2SO4/0.01 M cacodylate or phosphate (pH 7.0). dAMP (1 mM) and dGMP (1 mM) were in 0.3 M Na2SO4/0.01 M cacodylate (pH 7.0). All samples (2 ml) were directly titrated at 25oC using a combination pH electrode (MI-410, Microelectrodes, Inc.) and filtered using a 0.2 mm Nylon-66 membrane (Rainin). Samples in D2O solvent were obtained by freezing samples in H2O, lyophilizing and dissolving at 4oC in D2O. Oligonucleotide integrity was verified after irradiation (27 ).
Laser radiation was generated by frequency doubling the output of a 300 Hz excimer-pumped dye laser (Lambda Physik LPX130/FL3002). To generate 240 and 250 nm excitation, coumarin 480 (Exciton) was used as the lasing dye. For 210 nm excitation, stilbene 420 (Exciton) was the lasing dye. UV radiation at an average power of 0.3 mW was incident on the sample using a 135o backscattering geometry. The UVRR spectrometer has been described previously (28 ). Samples in a stirred quartz cuvette were maintained at specified temperatures with an aluminum block cell holder and circulating bath, with sample temperature monitored by a copper-constantan thermocouple in the block.
Spectra resulting from 0.5 to 1 h of data collection were corrected for the response of the monochromator and the detector by measuring the throughput of a deuterium lamp. Curve-fitting of the spectra was accomplished using Labcalc software (Galactic Instruments Corp.). The curves were generated assuming a 50% Lorentzian and 50% Gaussian line shape with a fixed width of 20 cm-1. The quality of the fit was assessed by both visual inspection and examination of the [chi]2 values resulting from the deviations between the fit simulations and the actual data.
The 250 nm-excited UVRR spectrum of the homopurine oligomer d(A-G)10 (Fig. 2 ) closely resembles the co-added monomer spectrum, except that the intensity is diminished by ~65% relative to the monomer intensities at 4oC (Table 1 and 2 ). The bands arise from in-plane ring stretching vibrations, and the suppression of Raman intensity results from base stacking interactions (29 ).
Involvement of the C6 carbonyl. The nature of base-pairing was investigated with 210 nm excitation, in order to resonantly enhance scattering from the carbonyl moieties of the dG residues (31 ,33 ). The C6=O stretching coordinate is the main contributor to the d(A-G)10 band seen at 1680 cm-1 in H2O (Fig. 7 ) or at 1654 cm-1 in D2O (Fig. 8 ); the D2O downshift reflects a substantial contribution from the NH bend of the adjacent imine group (32 ). Both bands shift up substantially when the temperature is raised from 4 to 80oC, 7 cm-1 in H2O (Fig. 7 ) and 13 cm-1 in D2O (Fig. 8 ). These upshifts imply a significant strengthening of the C6=O bond at high temperature, consistent with a loss of H-bonding (34 ). dGMP by itself shows an upshift in this band (data not shown), but the extent of the shift between 4 and 80oC is only 5 cm-1 in either H2O or D2O. Weakening of the H-bonding to H2O molecules can account for the temperature effect on dGMP. The larger shifts seen for d(A-G)10 reflect disruption of stronger H-bonds, consistent with G.G base pairing at low temperature. The same magnitude of shift was also seen with dGMP, which is known to form helically wound G-tetrads (11 ) when the concentration is brought to 10 mg/ml (data not shown). We note that in helical structures of poly(G) and poly(dG-dC), the C6=O band splits and the shift pattern is more complex because of strong interbase coupling between C=O modes of adjacent residues (35 -37 ). The absence of these effects in d(A-G)10 indicates that adjacent bases along the strands are not coupled vibrationally, consistent with the uncorrelated fluctuations of dG and dA residues implied by the melting behavior, and that the bases themselves deviate significantly from coplanarity (32 ).
Figure
Figure
Involvement of the exocyclic NH2. The exocyclic -NH2 scissors mode gives rise to weakly enhanced UVRR bands near 1600 cm-1 for both dA and dG residues in DNA (31 ). In D2O, exchange of NH2 for ND2 shifts the mode to ~1195 cm-1 for dA residues and ~1210 cm-1 for dG residues. It is therefore revealing that the UVRR difference spectrum of d(A-G)10 between 80 and 4oC in D2O (Fig. 8 ) shows a bisignate band with a trough at 1235 cm-1 and a peak at 1198 cm-1. We interpret this signal as arising from an upshift of the -ND2 scissors mode of the dG residues at low temperature, due to stronger H-bonding, an interaction which is expected to increase the N-D bending force constant. Comparable upshifts have been reported for 9-ethyladenine in proton-accepting solvents and support our assignment of frequency upshifts to H-bond formation (38 ). The H-bonding effect is less obvious in the H2O spectra, because of overlap between the G and A -NH2 scissors modes. Nevertheless, a close examination of the 1600 cm-1 region (Fig. 9 ) reveals a significant broadening of the composite band at low temperature. Curve resolution yields two components, 1598 and 1609 cm-1, at 4oC, but only one component at 1601 cm-1, at 80oC. We interpret the broadening at low temperature as an upshift in the -NH2 scissors frequency of the dG residues due to H-bonding, just as in the D2O spectrum.
Figure
Involvement of guanine N7. The spectra also provide evidence for H-bonding to N7 of the dG residues. When excited at 240 nm, the strong band at 1484 cm-1 shifts to 1487 cm-1 as the temperature is raised from 4 to 80oC (Fig. 5 , inset). This band is due mainly to a G ring mode involving N7=C8 stretching and C8-H bending coordinates (Fig. 2 , inset) (31 ,34 ). The frequency of this mode is sensitive to the C8 substituent and to any molecular interactions at N7 (34 ,37 ). Such a downshift was also observed for 3',5' GpG (37 ) at low temperature (1490 to 1481 cm-1 between 80 and 22oC) and attributed to N7[middot][middot][middot]H2N intermolecular H-bonding in helical aggregates of stacked GpG tetramers. Likewise, downshifts of this band were observed in helices of dGMP gels and of poly(dG) (29 ). It is important to note that with 250 nm excitation, where an A ring mode dominates the 1485 cm-1 band, there is no temperature effect on the frequency (Fig. 3 ). The temperature effect observed with 240 nm excitation is therefore attributed to H-bond involvement by the N7 atoms of G at low temperature.
Further evidence for H-bonding at N7 of G is detected in the course of melting in D2O (Fig. 8 ); a band at 1338 cm-1 shifts to 1345 cm-1 upon thermal denaturation. A new band at 1348 cm-1 in the 80-4oC difference spectrum additionally confirms this frequency upshift (Fig. 8 ). Fodor et al. (31 ) assigned the G vibrational mode observed at 1364 cm-1 in H2O to in-phase stretching motions of the triene linkage: C8=N7-C5=C4-N3=C2. The intensity of this mode is amplified with 210 nm excitation and its frequency is lowered in D2O (31 ,39 ); thus, the band observed at 1338 cm-1 in d(A-G)10 is assigned to this stretching vibration. 15N substitution in model compounds has shown that N7 motion significantly affects the frequency of this band (39 ). Hence, we suggest that the upshift of this mode in the oligomer at higher temperature reflects a loss of H-bonding at N7.
UVRR melting studies were also examined for evidence of H-bonding of the dA residues. Previous studies (27 ,38 ) have shown that H-bonding or protonation at the acceptor sites of dA residues (N1, N3 and N7) induces an upshift of ring frequencies above 1400 cm-1 with concomitant downshifts of the modes at 1335 and 1310 cm-1. For d(A-G)10, when dA residues are preferentially selected by 250 nm excitation, the frequencies of these bands remain constant during thermal denaturation, implying that H-bonding at these sites is extremely weak or non-existent. As described above, the downshift of the -NH2 scissors mode upon deuteration and subsequent thermally-induced wavenumber shift implicates -NH2 groups of dG but not dA residues in H-bonding at low temperature (Fig. 8 ). Figure 1 (bottom) depicts a model for A-tetrads that is isostructural with the proposed G-tetrads, but it requires H-bonding between N1 and -NH2 of dA residues of adjacent strands. The UVRR studies show no evidence for H-bonding at those sites, nor at any other H-bond acceptor sites on the dA residues. Therefore, we conclude that H-bonded A-tetrads are not present in d(A-G)10.
As a consequence of vibrational coupling between base and ribose rings, UVRR spectra in the 600-900 cm-1 region (Fig. 10 ) contain bands arising from ring vibrations, which are characteristic of glycosyl bond orientation and sugar pucker. The bands observed at 682 and 863 cm-1 in the 4oC spectrum reflect dG residues which are in an anti orientation with respect to deoxyribose rings in the C2'-endo conformation (Fig. 10 ) (32 ,40 ). When the temperature is raised, these bands become broader and weaker, reflecting increased flexibility of the backbone and a broadening of the range of accessible conformers. Figure 11 shows that the diminution of the 683 cm-1 band peak height exhibits the same cooperative melting behavior as do the high frequency G mode enhancements (Fig. 6 ). Thus, the conformational flexibility of the dG residues correlates with unstacking. The 4oC spectrum (Fig. 10 ) also shows weak bands at 628 and 847 cm-1 positions characteristic of dG residues in a C3'-endo/syn configuration. Thus, a fraction of the dG residues in d(A-G)10 may adopt syn orientations at low temperature, but by 20oC, these bands have disappeared. As the UV Raman cross sections for these bands are not well known, the population of dG residues that are C3'-endo/syn cannot be accurately estimated; however, the observation of both bands confirms that some of the dG residues do adopt the syn configuration.
Figure
The 730 cm-1 band in the d(A-G)10 spectrum is characteristic of dA residues in an anti orientation with C2'-endo sugar conformation (41 ). This band becomes stronger with increasing temperature, just as do the high-frequency bands of the dA residues (Fig. 4 ). However, the 730 cm-1 band augmentation is not linear over the whole temperature range (Fig. 11 ), but only begins when the temperature is raised above 20oC. One possible interpretation is that some of the dA residues are in the syn orientation at 4oC. Such residues are expected to contribute to the 628 cm-1 UVRR band (41 ), along with any dG syn residues. Since the 628 cm-1 band disappears when the temperature is raised to 20oC, the syn residues (being in a higher energy conformation) melt before the main population of dA residues, which remain in the C2'-endo/anti configuration. Theseresidues, which had melted out of the structure, would sample a mixture of configurations and therefore, not contribute to any particular Raman band. Both melting stages contribute to the high frequency A mode hyperchromism, which has an overall linear temperature dependence. This Raman band is composed primarily of pyrimidine modes, and with 250 nm excitation exhibits little conformational sensitivity (39 ).
Figure
The UVRR observations presented here are fully consistent with the presence of G-tetrads in d(A-G)10 at low temperature. Complex formation from two strands (1 ) implicates base-base association between strands, but the uncorrelated melting of the dA and dG residues would seem to rule out H-bonding interactions between them. Moreover, the involvement of the -NH2 substituent of dG but not of dA residues in interbase H-bonding is implicated by low-temperature upshifts in only the G-associated NH2 or ND2 scissors mode frequencies; and H-bond acceptance by the C6=O carbonyl groups and N7 atoms of G is indicated by frequency downshifts of appropriate modes. These H-bond interactions are entirely consistent with the G-tetrad model in Figure 1 (top), in which the C6=O and N1H groups form H-bonds, as do the exocyclic -NH2 groups and the N7 atoms.
Some uncertainty exists regarding the assignment of the -NH2 scissors vibration. In a recent report (42 ), the spectral effects of D2O exchange of the amino and imino protons in the oligomer (dG)12 were examined and the 1603 cm-1 band was attributed to N1H in-plane bending, while a band at 1644 cm-1 (not observed in our data) was attributed to scissoring of the NH2 amino group. We note that the ~1604 and 1640 cm-1 bands both arise primarily from ring vibrations which have variable contributions from N1H bending and -NH2 scissoring motions (43 ). In the case of adenine and its N9-derivatives, the relative contribution of -NH2 scissoring to the modes at ~1604 and 1635 cm-1 strongly depends upon the degree of H-bonding. -NH2 scissoring motion primarily contributes to the lower frequency band in weak to moderate H-bonding solutions, but to the higher frequency band in stronger H-bonding solutions or in the crystalline state (38 ). In the case of (dG)12, as for linear homo G tetraplexes in general, the H-bonding is considerably stronger and the melting more cooperative (e.g. ref. 12 ) than in d(A-G)10 (1 ), which argues for a less stable structure for d(A-G)10 and consequently, weaker H-bonds. It is possible that in guanine as well as adenine the relative -NH2 contribution to the 1640 and 1604 cm-1 modes is strongly dependent on the degree of H-bonding and could be reversed for the two sequences, d(A-G)10 and (dG)12. However, we note that even if the band observed arises mainly from N1H bending rather than -NH2 scissors vibrations, the detected frequency shifts would reflect H-bond interactions of the C6=O and N1H groups and would still be consistent with G-tetrad formation.
Figure
In contrast, the UVRR spectra exhibit no evidence for interbase H-bonding of the dA residues in d(A-G)10, even at low temperature, as shown by the relatively temperature-insensitive frequencies of dA vibrational modes. These findings are consistent with the conclusion (23 ) that in telomeric sequences containing dA residues, there are G-tetrads but not A-tetrads.
Melting of the G-tetrads in d(A-G)10 occurs more or less cooperatively, with a 42oC transition temperature, and leads to unstacking of the dG residues, as seen in the UVRR hyperchromism of the G ring modes. Melting also disorders the glycosyl bond orientations and the sugar pucker, as revealed by the conformationally-sensitive G vibrations. Because of the absence of anchoring H-bond interactions, the dA residues melt gradually. Nevertheless, it is apparent that the dA residues, like the dG residues, contribute to the UVRR hypochromism at low temperature. This must mean that at low temperature the dA residues incur some stabilization by intercalating and stacking with their nearest neighbor G-tetrads.
The fact that the dA residues do not form H-bonded tetrads is probably because they cannot form a tetrad of sufficient stability that is isostructural with the G-tetrads (Fig. 1 ). Hence, as indicated by the UV melting profiles of Shiber et al. (1 ), the dA residues melt gradually over the whole temperature range; at low temperatures by becoming extrahelical, then by stacking with nearest neighbor dG residues as those melt, and then at still higher temperatures by unstacking of the dA-dG stacked neighbors.
Although UVRR observations shed light on the nature of base stacking and the interbase interactions, they do not directly identify the relative orientation of the strands. From the results of Shiber et al. (1 ) it is possible to distinguish between the four-stranded linear and two-stranded hairpin duplex structures, and it is the hairpin duplex structure that is consistent with experimental observations (Fig. 12 ). Our UVRR observations are fully compatible with this model. Although Rippe et al. (44 ) proposed parallel-stranded duplexes for alternating d(G-A) sequences, the indication that some of the residues are in the syn orientation (bands at 628 and 847 cm-1 in Fig. 10 ) strongly supports an antiparallel strand orientation. Moreover, the model of Rippe et al. (44 ) does not anticipate our findings regarding the presence of G-tetrads. From a Raman study of d(G4T4G4), (45 ) it was concluded that syn orientations are associated only with the antiparallel hairpin conformation and not with a linear parallel conformation. Other spectroscopic investigations (8 ) have also shown that bases involved in the same tetrad have opposite glycosyl torsions if the bases are on antiparallel strands. In those studies, however, G-residues adopt a C2'-endo/syn configuration and not the C3'-endo/syn configuration observed in the present study. Definitive evidence for the orientation of the strands must await further spectroscopic investigation.
We thank Horst Klump for helpful discussions. This work was supported by grants from the National Institutes of Health (GM 42936) to J.R.F. and (GM 25158) to T.G.S. M.C.S. was supported in part by a Sterling postdoctoral fellowship, and I.M. by a NRSA fellowship (GM 14324).
This is paper No. 25 in the series Polynucleotides, of which the last is ref. 42 .
*To whom correspondence should be addressed. Tel: +1 609 258 3927; Fax: +1 609 258 6730; Email: jrfresco@princeton.edu or spiro@princeton.edu
+Present Address: Department of Molecular Biology and Biochemistry, Wesleyan University, Middletown, CT 06459-0175, USA
REFERENCES
Return

