Thermal difference spectra: a specific signature for nucleic acid stru
http://www.100md.com
《核酸研究医学期刊》
Laboratoire de Biophysique, Muséum National d'Histoire Naturelle USM503, INSERM U 565, CNRS UMR 5153 43 rue Cuvier, 75231 Paris cedex 05, France 1Department of Biochemistry, University of Mississippi Medical Center 2500 N. State St., Jackson, MS 39216-4505, USA 2James Graham Brown Cancer Center, Department of Medicine, Health Sciences Center, University of Louisville 529 South Jackson Street, Louisville KY 40202, USA
*To whom correspondence should be addressed. Tel: +33 1 40 79 36 89; Fax: +33 1 40 79 37 05; Email: faucon@mnhn.fr
ABSTRACT
We show that nucleic acid structures may be conveniently and inexpensively characterized by their UV thermal difference spectra. A thermal difference spectrum (TDS) is obtained for a nucleic acid by simply recording the ultraviolet absorbance spectra of the unfolded and folded states at temperatures above and below its melting temperature (Tm). The difference between these two spectra is the TDS. The TDS has a specific shape that is unique for each type of nucleic acid structure, a conclusion that is based on a comparison of >900 spectra from 200 different sequences. The shape of the TDS reflects the subtleties of base stacking interactions that occur uniquely within each type of nucleic acid structure. TDS provides a simple, inexpensive and rapid method to obtain structural insight into nucleic acid structures, which is applicable to both DNA and RNA from short oligomers to polynucleotides. TDS complements circular dichroism as a tool for the structural characterization of nucleic acids in solution.
INTRODUCTION
DNA and RNA may adopt a number of unusual conformations. Several laboratories have shown the existence and functional importance of nucleic acid structures that are fundamentally different from the classical, canonical B-form duplex. A few years ago, Alexander Rich wrote: ‘DNA comes in many forms’ (1). This is true at all levels of structure. Different base pairs can form (G.G, A.A, etc.) and the helicity and/or orientation of the strands can be altered such as in Z-DNA (2), parallel-stranded (3) and Hoogsteen duplexes. One can also observe nucleic acid structures composed of more than two strands such as triple- (4,5) and quadruple- (6–12) helices. All these structures rely upon the formation of specific base-association schemes (base pairs, triplets or quartets) (1). The stability of these various conformational forms can be measured by a thermal denaturation experiment. Heating a structured nucleic acid leads to changes in its ultraviolet absorbance, which reflect the conformational change of the molecule in solution, specifically a disruption of base-stacking interactions. Such melting experiments can be analyzed to determine the thermodynamic basis of DNA or RNA structure stability. In most cases, one may perform thermal denaturation experiments by recording absorbance at 260 nm as a function of temperature. However, it may also be preferable to record absorbance at other wavelengths. For example, quadruplex denaturation does not lead to large variation in absorbance at 260 nm. It is, however, possible to follow this denaturation at other wavelengths such as 295 nm (13). This observation prompted us to determine the optimal wavelengths for thermal denaturation experiments for each type structure. In the course of our investigations, we realized that thermal difference spectra (TDS) provided a distinctive spectroscopic signature for each nucleic acid structure. We describe the method and its utility here.
MATERIALS AND METHODS
Nucleic acids samples
Most oligonucleotides presented here were synthesized by Eurogentec (Belgium) or Genosys on the 40–1000 nmol scale and resuspended in 200–500 μl of bi-distilled water. We checked the purity of most samples by mass spectrometry, HPLC or gel electrophoresis. Polynucleotide samples were purchased from Amersham-Pharmacia (Piscataway, NJ) and were used without further purification.
Determination of extinction coefficients
Concentration of all oligodeoxynucleotides were estimated by UV absorption using published sequence-dependent extinction coefficients (14). In selected cases, molar extinction coefficients were determined by a colorimetric phosphate analysis (15). One should not underestimate the difficulty in obtaining reliable extinction coefficients for all oligonucleotides. For example, we measured for the human telomeric repeat sequence 5'-AGGG(TTAG3)3, a molar extinction coefficient of 196 479 ± 5675 M–1 cm–1 by a colorimetric phosphate assay. Using two commonly used algorithms, values for the molar extinction coefficient of 251 800 and 228 500 M–1 cm–1 were calculated. The calculated values differ by 16–28% from the more accurate analytically determined value.
Absorbance spectra
All experiments presented here were performed on either Kontron Uvikons 940/943 or a Cary 3E UV/Vis spectrophotometer using quartz cuvettes with an optical pathlength of 0.2–1 cm (16). Absorbance spectra were recorded in the 220–335 nm range, with a scan speed of 200–1000 nm/min, with a data interval of 1 nm. In most cases, using cuvettes of 0.05–1 cm pathlength, maximum absorbance of the sample was 0.5–2.0.
Choice of the upper and lower temperatures
A melting curve is usually measured first to learn the width of the transition. The difference in spectra weakly depends on the difference in temperature. Using the narrowest possible temperature range is, in principle, recommended. However, one can still have an excellent superposition for most structures if one compares various temperature windows for an i-DNA sample, using a number of different temperatures. Clearly, once the curves are normalized, the shapes are almost superimposable (Supplementary Figure S1B), demonstrating that the choice of the exact temperatures does not play a great role.
Experimental constraints
Reproducibility of the TDS is excellent even if spectra are recorded in a relatively crude way. There is no need for very slow scans and 1 pt/nm is enough for most applications. Any high-quality single- or double-beam spectrophotometer may be used, provided that its performance and stability over time in the far-UV region (and at high temperature) is checked. Either a Peltier thermoelectric temperature controller or a circulating water bath may be used, as long as a wide temperature range can be used. Appropriate precautions should be taken below room temperature to avoid condensation, and at high temperature to avoid evaporation and air bubbles. Thorough degassing of solutions helps to minimize bubble formation at high temperatures. Layering of solutions with mineral oil and tightly stoppering cuvettes help to minimize evaporation. Condensation at low temperatures can be minimized by blowing a stream of dry air on the outer cuvette surfaces. It is, of course, essential to choose a buffer that does not absorb light in the far-UV region. Cacodylate (pKa: 6.14 at 25°C) and acetate (pKa: 4.62) were preferred over phosphate or Tris buffer for experiments performed at neutral or slightly acidic pH. Their pKa is almost temperature independent.
RESULTS AND DISCUSSION
An example of the determination of a TDS is shown in Figure 1. The structure in this determination is the antiparallel ‘basket’ G-quadruplex form observed in Na+ solution for the human telomere repeat sequence 5'-AGGG(TTAGGG)3 (17). Figure 1A shows the circular dichroic spectra of the quadruplex. Figure 1B shows UV absorbance spectra at temperatures well above (90°C) and below (20°C) the melting temperature (56°C). Figure 1C shows the TDS obtained by subtracting the low temperature spectrum from the high temperature spectrum. In this case, the TDS is displayed as the difference between molar extinction coefficients between the unfolded and folded forms. Other methods for displaying the TDS that are convenient and informative are described below.
Figure 1 Example of a thermal difference spectrum. The structure in this determination is the antiparallel ‘basket’ G-quadruplex form observed in Na+ solution for the human telomere repeat sequence 5'-AGGG(TTAGGG)3. (A) Circular dichroic spectrum, expressed as molar ellipticity. (B) UV absorbance spectra (expressed as molar extinction, M–1 cm–1) at 20°C (solid line) and 90°C (dashed line). (C). Thermal difference spectrum resulting from the subtraction of the 20°C spectrum in (B) from the 90°C spectrum.
Examples of TDS for nine different structure types are shown in Figure 2. The curves are simply the result of the arithmetic difference between high and low temperature absorbance spectra (example shown in Figure 1B). For most experiments, high temperature was chosen >80°C, low temperature 20°C or below. These difference spectra were then normalized in order to facilitate the comparison of the spectral shapes, simply by dividing the raw data by its maximum value, so that the highest positive peak gets a Y-value of +1. This normalization, which unfortunately leads to a loss of information, gives in our hands the most reproducible results. In principle, molar extinction differences should also be interesting. However, our opinion is that, unless analytically determined extinction coefficients are available for all oligonucleotides (see Materials and Methods section), it is preferable and more reliable to use the normalization procedure. Analysis of the normalized TDS offers another advantage, as it allows to deal with partial transitions when (i) the Tm is too low or too high (no spectra of the pure folded/unfolded species) (ii) two transitions are overlapping. One may also analyze the shape of the ‘beginning’ of a transition (i.e. the part corresponding to partial dissociation at low temperature, below Tm) and compare it with the global shape or the ‘end’ (see Supplementary Figure S1 for an example).
Figure 2 Normalized differential absorbance signatures: (A) DNA self-complementary duplexes, 100% AT; (B) DNA self-complementary duplexes 100% GC; (C) Z-DNA; (D) Parallel-stranded DNA; (E) GA DNA duplexes; (F) Hoogsteen DNA duplexes; (G) i-DNA; (H) Pyrimidine triplexes; (I) DNA G-quadruplexes in Na+. The complete legend as well as expanded spectra and relevant references may be found as Supplementary Figure S2.
The resulting normalized TDS are concentration independent (data not shown). Except for pyrimidine triplexes (Figure 2H) and to a lesser extent G-quadruplexes (Figure 2I), one cannot but notice that all spectra within a structural family are strikingly similar (e.g. Figure 2A, B and G). Sometimes spectra are essentially identical (Figure D and E). Comparison of these different structures reveals reproducible differences in TDS shapes, allowing an investigator to determine the structural type by the spectrum alone.
The robustness of this method was analyzed by several control experiments (see Figure 3). A simple DNA self-complementary duplex (5'-d-AAATTTAATTAAATTT) was chosen; its differential absorbance spectra was determined at various pH (5.5–7.5), salt concentration, in the presence or absence of magnesium and with different monocations (sodium, lithium or potassium). All spectra were nearly superimposable, despite the fact that the Tm of the duplex may depend on these conditions: as long as the melting transition occurs in the temperature interval defined by the low and high-temperature spectra, very reproducible differences were found (Figure 3A). Similar observations were found for most structures, with a few notable exceptions. For G-quadruplexes, this profile was independent on pH, magnesium, ionic strength, but not on the nature of the monocation . A systematic difference was found between the spectra recorded in sodium and in potassium (Figure 3B). For i-DNA, the spectra appear to be different when working below pH 6: at this pH, a significant proportion of cytosines are protonated in the unfolded form, and C+ has a different absorbance from unprotonated C .
Figure 3 Experimental condition dependency of the differential absorbance shapes. (A) For a 16-base long self-complementary DNA duplex 5'-AAATTTAATTAAATTT; pH 5.5–7.2, with or without 10 mM MgCl2, in a sodium cacodylate buffer with 0.1–1 M NaCl, KCl or LiCl. All curves are superimposable. (B) For the 5'-AGGGTTAGGGTTAGGGTTAGGG intramolecular DNA quadruplex pH 6.0–7.2, in a sodium cacodylate buffer with 0.1–1 M NaCl (in blue), or 0.03–0.6 M KCl (in red). Note that CD spectra of the same sequence are shown as Supplementary Figure S3, and error bars at each wavelength (average of 11 independent measurements) as Supplementary Figure S4.
The first obvious application of this analysis is to help in the choice of appropriate wavelength(s) for a thermal denaturation experiment for more detailed thermodynamic studies. These results have been collected in Table 1. It should be noted that the relation between a differential peak and a ‘good’ wavelength is not always direct. In some cases, the observed hypochromicity does not result from the melting of a structure, but is the consequence of temperature-dependent extinction coefficient of the folded and/or the unfolded forms. For example, the positive peak at 273 nm observed in most quadruplexes (Figure 2I) results from a significant temperature-dependent increase of absorbance of the single strands. This is related to the classical problem of the sloping baseline analysis of melting profiles: all species exhibit temperature-dependent absorbance properties. One would then argue that the most convenient differential spectra should be isothermal; however, this requires the ability to record the absorbance of the folded and the unfolded states at the same temperature. This is possible in only a limited number of cases (21), when association and dissociation of the structure are very slow. Otherwise, T-jump or fast mixing experiments are required. Despite this limitation, this simple dual-temperature measurement is sufficient for most determinations.
Table 1 Useful wavelengths for each nucleic acid structure
These spectra provide information complementary to circular dichroism. The UV absorbance of nucleic acids from 200 to 300 nm is exclusively due to transitions of the planar purine and pyrimidine bases. The phosphate backbone begins to contribute at 190 nm. The CD spectrum is probing the asymmetry of bases. At wavelengths >200 nm, the asymmetry arises from CD induced into the bases as a result of stacking interactions and coupling with backbone transitions. Hypochromicity is mainly the result of base stacking (22). Therefore, one can expect that the relative positioning of the bases will play a crucial role in this differential spectra. The theory of hypochromicity is complex, the changes in light absorption occurring upon denaturation is attributed to the interaction between the dipoles induced in the chromophore by the light (22). It is striking that a number of very different structures involve a negative differential absorbance at 295 nm (Table 1). This might be the result of the significant contribution of n * transition moments at this wavelength: as these transition moments are parallel to the helix axis (23), increased stacking resulting from the formation of a folded structure may result in hyperchromism. However, the absorbance spectra of nucleic acid bases are deceptively simple: they actually correspond to a number of different * and n * transitions. Hypochromism at a given wavelength will not only depend on the intrinsic transition moments of each base, but also on the relative moment of the interacting bases (base-paired or stacked). Such complexity will make an accurate calculation of predicted differential absorbances a difficult task.
Although the dissociation of all G-quadruplexes analyzed so far show an inverted transition at 295 nm (13), this property is not unambiguously indicative of G-quadruplexes. Z-DNA (Figure 2C) i-motif (Figure 2G) Hoogsteen duplexes (Figure 2F) and pyrimidine triplex (Figure 2H) formation also leads to an inverted transition at 295 nm. In Figure 2F–H, this is the result of the differences in absorbance between C and C+ (20). A simple comparison of these spectra indicates that the global shape rather than the wavelength of the peak is specific for each structure (Figure 2).
These observations suggest that the simple recording of two absorbance spectra might be sufficient to obtain qualitative structural information on the identity of the folded form. This was a common accepted practice in the sixties, when the melting of polynucleotides was analyzed . Unfortunately, with a few notable exceptions (27–32), complete spectra are rarely presented in contemporary reports. We demonstrate here that in many cases this difference is very informative. Of course, this requires that the high and low temperatures correspond to the unfolded and folded forms, respectively. In other words, the Tm of the structure must lie between these two values. Although the spectra of different nucleic acids that fold into the same type of structure are not identical, significant homology was found between the spectra. The uniqueness of these spectra is sufficient, in most cases, to suggest a type of folding. However, we could not confidently distinguish the differential spectra of a purine-rich triplex (the so-called GA or GT motif) from the melting of a duplex (data not shown). This is in agreement with the conclusions that the formation of these triplexes occurs without significant hypochromicity (33). Similarly, no significant difference was found between a mismatched and a fully complementary DNA duplex (data not shown).
We have now collected >900 spectra, corresponding to >200 different nucleic acid sequences (mostly DNA), each tested 2–26 times under different experimental conditions. This is, to our knowledge, the first attempt to systematically collect thermal difference spectra as initially proposed 40 years ago by Fresco et al (34). Our goal is to first build an open access database of the spectra of nucleic acids and to propose a predictive algorithm that provides a probability function of a given structure. This should be helpful in the discovery of new nucleic acid folding schemes, and in the crude determination of the global folding of an oligonucleotide (which, for example, might have been found by SELEX). We are now applying this method to the analysis of trinucleotide repeats (35) and to samples which may adopt complex conformations, involving more than one simple structure. Efforts are also being made on RNA structures (and corresponding DNA strands carrying deoxyuracyl bases): RNA duplex and quadruplex structures exhibit related but distinct signatures from the corresponding DNA duplexes (data not shown) and quadruplexes (18,21), respectively. We also wish to analyze the signature of modified bases (such as 5-methyl cytosines) and the impact of modifying the sugar/phosphate backbone (such as 2'-O-methyl phosphorothioates, LNA or PNA).
SUPPLEMENTARY DATA
Supplementary Data are available at NAR Online.
ACKNOWLEDGEMENTS
We thank M. Rougée, T. Garestier, L. Guittat, P. Arimondo, P. Alberti, B. Saccà (MNHN, Paris, France) M. Mills (UCT, Cape Town, Republic of South Africa) for helpful discussions. This work was supported by an ARC grant (#3365 to J.L.M), a EU FP6 MOLCANCERMED grant (to J.L.M.), a fondation Jérome Lejeune research grant (to S.A.) and National Cancer Institute grant CA35635 (to J.B.C). Funding to pay the Open Access publication charges for this article was provided by INSERM.
Conflict of interest statement: None declared.
REFERENCES
Rich, A. (1993) DNA comes in many forms Gene, 135, 99–109 .
Rich, A. and Zhang, S.G. (2003) Z-DNA: the long road to biological function Nature Rev. Genet., 4, 566–572 .
Jovin, T.M., Rippe, K., Ramsing, N.B., Klement, R., Elhorst, W., Vojtisova, M. (1990) Parallel stranded DNA In Sarma, R.H. and Sarma, M.H. (Eds.). Structure & Methods, DNA & RNA, Schenectady, New York Adenine Press Vol 3, pp. 155–174 .
Moser, H.E. and Dervan, P.B. (1987) Sequence specific cleavage of double helical DNA by triple helix formation Science, 238, 645–650 .
Le Doan, T., Perrouault, L., Praseuth, D., Habhoub, N., Decout, J.-L., Thuong, N.T., Lhomme, J., Hélène, C. (1987) Sequence specific recognition, photocrosslinking and cleavage of the DNA double helix by an oligo thymidylate covalently linked to an azidoproflavine derivative Nucleic Acids Res., 15, 7749–7760 .
Gellert, M., Lipsett, M.N., Davies, D.R. (1962) Helix formation by guanylic acid Proc. Natl Acad. Sci. USA, 48, 2013–2018 .
Gehring, K., Leroy, J.L., Guéron, M. (1993) A tetrameric structure with protonated cytosine.cytosine base pairs Nature, 363, 561–565 .
Henderson, E., Hardin, C.C., Walk, S.K., Tinoco, I., Jr, Blackburn, E.H. (1987) Telomeric DNA oligonucleotides form novel intramolecular structures containing guanine.guanine base pairs Cell, 51, 899–908 .
Sundquist, W.I. and Klug, A. (1989) Telomeric DNA dimerizes by formation of guanine tetrads between hairpin loops Nature, 342, 825–829 .
Williamson, J.R., Raghuraman, M.K., Cech, T.R. (1989) Monovalent cation induced structure of telomeric DNA: the G-quartet model Cell, 59, 871–880 .
Kang, C.H., Berger, I., Lockshin, C., Radliff, R., Moyzis, R., Rich, A. (1994) Crystal structure of intercalated four stranded d(C3T) at 1.4 ? resolution Proc. Natl Acad. Sci. USA, 91, 11636–11640 .
Parkinson, G.N., Lee, M.P.H., Neidle, S. (2002) Crystal structure of parallel quadruplexes from human telomeric DNA Nature, 417, 876–880 .
Mergny, J.L., Phan, A.T., Lacroix, L. (1998) Following G-quartet formation by UV-spectroscopy FEBS Lett., 435, 74–78 .
Cantor, C.R., Warshaw, M.M., Shapiro, H. (1970) Oligonucleotide interactions. 3. Circular dichroism studies of the conformation of deoxyoligonucleotides Biopolymers, 9, 1059–1077 .
Plum, G.E. (2000) Determination of oligonucleotide molar extinction coefficients Current Protocols in Nucleic Acid Chemistry, New York Wiley pp. 7.3.1–7.3.17 .
Mergny, J.L. and Lacroix, L. (2003) Analysis of thermal melting curves Oligonucleotides, 13, 515–537 .
Wang, Y. and Patel, D.J. (1993) Solution structure of the human telomeric repeat d G-tetraplex Structure, 1, 263–282 .
Saccà, B., Lacroix, L., Mergny, J.L. (2005) The effect of chemical modifications on the thermal stability of different G-quadruplexes-forming oligonucleotides Nucleic Acids Res., 33, 1182–1192 .
Li, J., Correia, J.J., Wang, L., Trent, J.O., Chaires, J.B. (2005) Not so crystal clear: the structure of the human telomeric G-quadruplex in solution differs from that present in a crystal Nucleic Acids Res., 33, 4649–4659 .
Mergny, J.L., Lacroix, L., Han, X., Leroy, J.L., Hélène, C. (1995) Intramolecular folding of pyrimidine oligodeoxynucleotides into an i-DNA motif J. Am. Chem. Soc., 117, 8887–8898 .
Mergny, J.L., de Cian, A., Ghelab, A., Saccà, B., Lacroix, L. (2005) Kinetics of tetramolecular quadruplexes Nucleic Acids Res., 33, 81–94 .
Tinoco, I., Jr. (1960) Hypochromism in polynucleotides J. Am. Chem. Soc., 82, 4785–4790 .
Rich, A. and Kasha, M. (1960) The n->pi* transition in nucleic acids and polynucleotides J. Am. Chem. Soc., 82, 6197–6199 .
Riley, M., Maling, B., Chamberlin, M.J. (1966) Physical and chemical characterization of two and three stranded adenine thymine and adenine uracil homopolymer complexes J. Mol. Biol., 20, 359–389 .
Blake, R.D., Massoulié, J., Fresco, J.R. (1967) A spectral approach to the equilibria between polyriboadenylate and polyribouridylate and their complexes J. Mol. Biol., 30, 291–308 .
Martin, F.H., Uhlenbeck, O.C., Doty, P. (1971) Self complementary oligoribonucleotides: adenylic acid–uridylic acid block copolylers J. Mol. Biol., 57, 201–215 .
Ramsing, N.B. and Jovin, T.M. (1988) Parallel stranded duplex DNA Nucleic Acids Res., 16, 6659–6676 .
van de Sande, J.H. and Jovin, T.M. (1982) Z* DNA, the left handed helical form of poly d(G–C) in MgCl2-ethanol is biologically active EMBO J., 1, 115–120 .
Lavelle, L. and Fresco, J.R. (2003) Enhanced stabilization of the triplexes d(C+-T)(6):d(A-G)(6).d(C-T)(6), d(T)(21):d(A)21.d(T)(21) and poly r(U:A.U) by water structure-making solutes Biophys. Chem., 105, 701–720 .
Davis, T.M., McFail-Isom, L., Keane, E., Williams, L.D. (1998) Melting of a DNA hairpin without hyperchromism Biochemistry, 37, 6975–6978 .
Kankia, B.I. (2003) Mg2+-induced triplex formation of an equimolar mixture of poly(rA) and poly(rU) Nucleic Acids Res., 31, 5101–5107 .
Haq, I., Chowdhry, B.Z., Chaires, J.B. (1997) Singular value decomposition of 3-D DNA melting curves reveals complexity in the melting process Eur. Biophys. J., 26, 419–426 .
Mills, M., Arimondo, P., Lacroix, L., Garestier, T., Hélène, C., Klump, H.H., Mergny, J.L. (1999) Energetics of strand displacement reactions in triple helices: a spectroscopic study J. Mol. Biol., 291, 1035–1054 .
Fresco, J.R., Klotz, L.C., Richards, E.G. (1963) A new spectroscopic approach to the determination of helical secondary structure in ribonucleic acids Cold Spring Harbor Symp. Quant. Biol., 28, 83–90 .
Amrane, S., Saccà, B., Mills, M., Chauhan, M., Klump, H.H., Mergny, J.L. (2005) Length-dependent energetics of (CTG)n and (CAG)n trinucleotide repeats Nucleic Acids Res., 33, 4065–4077 .
Jose, D. and P?rschke, D. (2004) Dynamics of the B-A transition of DNA double-helices Nucleic Acids Res., 32, 2251–2258 .
Pohl, F.M. and Jovin, T.M. (1972) Salt-induced co-operative conformational change of a synthetic DNA: equilibrium and kinetic studies with poly(dG-dC) J. Mol. Biol., 67, 375–396 .(Jean-Louis Mergny*, Jing Li1, Laurent La)
*To whom correspondence should be addressed. Tel: +33 1 40 79 36 89; Fax: +33 1 40 79 37 05; Email: faucon@mnhn.fr
ABSTRACT
We show that nucleic acid structures may be conveniently and inexpensively characterized by their UV thermal difference spectra. A thermal difference spectrum (TDS) is obtained for a nucleic acid by simply recording the ultraviolet absorbance spectra of the unfolded and folded states at temperatures above and below its melting temperature (Tm). The difference between these two spectra is the TDS. The TDS has a specific shape that is unique for each type of nucleic acid structure, a conclusion that is based on a comparison of >900 spectra from 200 different sequences. The shape of the TDS reflects the subtleties of base stacking interactions that occur uniquely within each type of nucleic acid structure. TDS provides a simple, inexpensive and rapid method to obtain structural insight into nucleic acid structures, which is applicable to both DNA and RNA from short oligomers to polynucleotides. TDS complements circular dichroism as a tool for the structural characterization of nucleic acids in solution.
INTRODUCTION
DNA and RNA may adopt a number of unusual conformations. Several laboratories have shown the existence and functional importance of nucleic acid structures that are fundamentally different from the classical, canonical B-form duplex. A few years ago, Alexander Rich wrote: ‘DNA comes in many forms’ (1). This is true at all levels of structure. Different base pairs can form (G.G, A.A, etc.) and the helicity and/or orientation of the strands can be altered such as in Z-DNA (2), parallel-stranded (3) and Hoogsteen duplexes. One can also observe nucleic acid structures composed of more than two strands such as triple- (4,5) and quadruple- (6–12) helices. All these structures rely upon the formation of specific base-association schemes (base pairs, triplets or quartets) (1). The stability of these various conformational forms can be measured by a thermal denaturation experiment. Heating a structured nucleic acid leads to changes in its ultraviolet absorbance, which reflect the conformational change of the molecule in solution, specifically a disruption of base-stacking interactions. Such melting experiments can be analyzed to determine the thermodynamic basis of DNA or RNA structure stability. In most cases, one may perform thermal denaturation experiments by recording absorbance at 260 nm as a function of temperature. However, it may also be preferable to record absorbance at other wavelengths. For example, quadruplex denaturation does not lead to large variation in absorbance at 260 nm. It is, however, possible to follow this denaturation at other wavelengths such as 295 nm (13). This observation prompted us to determine the optimal wavelengths for thermal denaturation experiments for each type structure. In the course of our investigations, we realized that thermal difference spectra (TDS) provided a distinctive spectroscopic signature for each nucleic acid structure. We describe the method and its utility here.
MATERIALS AND METHODS
Nucleic acids samples
Most oligonucleotides presented here were synthesized by Eurogentec (Belgium) or Genosys on the 40–1000 nmol scale and resuspended in 200–500 μl of bi-distilled water. We checked the purity of most samples by mass spectrometry, HPLC or gel electrophoresis. Polynucleotide samples were purchased from Amersham-Pharmacia (Piscataway, NJ) and were used without further purification.
Determination of extinction coefficients
Concentration of all oligodeoxynucleotides were estimated by UV absorption using published sequence-dependent extinction coefficients (14). In selected cases, molar extinction coefficients were determined by a colorimetric phosphate analysis (15). One should not underestimate the difficulty in obtaining reliable extinction coefficients for all oligonucleotides. For example, we measured for the human telomeric repeat sequence 5'-AGGG(TTAG3)3, a molar extinction coefficient of 196 479 ± 5675 M–1 cm–1 by a colorimetric phosphate assay. Using two commonly used algorithms, values for the molar extinction coefficient of 251 800 and 228 500 M–1 cm–1 were calculated. The calculated values differ by 16–28% from the more accurate analytically determined value.
Absorbance spectra
All experiments presented here were performed on either Kontron Uvikons 940/943 or a Cary 3E UV/Vis spectrophotometer using quartz cuvettes with an optical pathlength of 0.2–1 cm (16). Absorbance spectra were recorded in the 220–335 nm range, with a scan speed of 200–1000 nm/min, with a data interval of 1 nm. In most cases, using cuvettes of 0.05–1 cm pathlength, maximum absorbance of the sample was 0.5–2.0.
Choice of the upper and lower temperatures
A melting curve is usually measured first to learn the width of the transition. The difference in spectra weakly depends on the difference in temperature. Using the narrowest possible temperature range is, in principle, recommended. However, one can still have an excellent superposition for most structures if one compares various temperature windows for an i-DNA sample, using a number of different temperatures. Clearly, once the curves are normalized, the shapes are almost superimposable (Supplementary Figure S1B), demonstrating that the choice of the exact temperatures does not play a great role.
Experimental constraints
Reproducibility of the TDS is excellent even if spectra are recorded in a relatively crude way. There is no need for very slow scans and 1 pt/nm is enough for most applications. Any high-quality single- or double-beam spectrophotometer may be used, provided that its performance and stability over time in the far-UV region (and at high temperature) is checked. Either a Peltier thermoelectric temperature controller or a circulating water bath may be used, as long as a wide temperature range can be used. Appropriate precautions should be taken below room temperature to avoid condensation, and at high temperature to avoid evaporation and air bubbles. Thorough degassing of solutions helps to minimize bubble formation at high temperatures. Layering of solutions with mineral oil and tightly stoppering cuvettes help to minimize evaporation. Condensation at low temperatures can be minimized by blowing a stream of dry air on the outer cuvette surfaces. It is, of course, essential to choose a buffer that does not absorb light in the far-UV region. Cacodylate (pKa: 6.14 at 25°C) and acetate (pKa: 4.62) were preferred over phosphate or Tris buffer for experiments performed at neutral or slightly acidic pH. Their pKa is almost temperature independent.
RESULTS AND DISCUSSION
An example of the determination of a TDS is shown in Figure 1. The structure in this determination is the antiparallel ‘basket’ G-quadruplex form observed in Na+ solution for the human telomere repeat sequence 5'-AGGG(TTAGGG)3 (17). Figure 1A shows the circular dichroic spectra of the quadruplex. Figure 1B shows UV absorbance spectra at temperatures well above (90°C) and below (20°C) the melting temperature (56°C). Figure 1C shows the TDS obtained by subtracting the low temperature spectrum from the high temperature spectrum. In this case, the TDS is displayed as the difference between molar extinction coefficients between the unfolded and folded forms. Other methods for displaying the TDS that are convenient and informative are described below.
Figure 1 Example of a thermal difference spectrum. The structure in this determination is the antiparallel ‘basket’ G-quadruplex form observed in Na+ solution for the human telomere repeat sequence 5'-AGGG(TTAGGG)3. (A) Circular dichroic spectrum, expressed as molar ellipticity. (B) UV absorbance spectra (expressed as molar extinction, M–1 cm–1) at 20°C (solid line) and 90°C (dashed line). (C). Thermal difference spectrum resulting from the subtraction of the 20°C spectrum in (B) from the 90°C spectrum.
Examples of TDS for nine different structure types are shown in Figure 2. The curves are simply the result of the arithmetic difference between high and low temperature absorbance spectra (example shown in Figure 1B). For most experiments, high temperature was chosen >80°C, low temperature 20°C or below. These difference spectra were then normalized in order to facilitate the comparison of the spectral shapes, simply by dividing the raw data by its maximum value, so that the highest positive peak gets a Y-value of +1. This normalization, which unfortunately leads to a loss of information, gives in our hands the most reproducible results. In principle, molar extinction differences should also be interesting. However, our opinion is that, unless analytically determined extinction coefficients are available for all oligonucleotides (see Materials and Methods section), it is preferable and more reliable to use the normalization procedure. Analysis of the normalized TDS offers another advantage, as it allows to deal with partial transitions when (i) the Tm is too low or too high (no spectra of the pure folded/unfolded species) (ii) two transitions are overlapping. One may also analyze the shape of the ‘beginning’ of a transition (i.e. the part corresponding to partial dissociation at low temperature, below Tm) and compare it with the global shape or the ‘end’ (see Supplementary Figure S1 for an example).
Figure 2 Normalized differential absorbance signatures: (A) DNA self-complementary duplexes, 100% AT; (B) DNA self-complementary duplexes 100% GC; (C) Z-DNA; (D) Parallel-stranded DNA; (E) GA DNA duplexes; (F) Hoogsteen DNA duplexes; (G) i-DNA; (H) Pyrimidine triplexes; (I) DNA G-quadruplexes in Na+. The complete legend as well as expanded spectra and relevant references may be found as Supplementary Figure S2.
The resulting normalized TDS are concentration independent (data not shown). Except for pyrimidine triplexes (Figure 2H) and to a lesser extent G-quadruplexes (Figure 2I), one cannot but notice that all spectra within a structural family are strikingly similar (e.g. Figure 2A, B and G). Sometimes spectra are essentially identical (Figure D and E). Comparison of these different structures reveals reproducible differences in TDS shapes, allowing an investigator to determine the structural type by the spectrum alone.
The robustness of this method was analyzed by several control experiments (see Figure 3). A simple DNA self-complementary duplex (5'-d-AAATTTAATTAAATTT) was chosen; its differential absorbance spectra was determined at various pH (5.5–7.5), salt concentration, in the presence or absence of magnesium and with different monocations (sodium, lithium or potassium). All spectra were nearly superimposable, despite the fact that the Tm of the duplex may depend on these conditions: as long as the melting transition occurs in the temperature interval defined by the low and high-temperature spectra, very reproducible differences were found (Figure 3A). Similar observations were found for most structures, with a few notable exceptions. For G-quadruplexes, this profile was independent on pH, magnesium, ionic strength, but not on the nature of the monocation . A systematic difference was found between the spectra recorded in sodium and in potassium (Figure 3B). For i-DNA, the spectra appear to be different when working below pH 6: at this pH, a significant proportion of cytosines are protonated in the unfolded form, and C+ has a different absorbance from unprotonated C .
Figure 3 Experimental condition dependency of the differential absorbance shapes. (A) For a 16-base long self-complementary DNA duplex 5'-AAATTTAATTAAATTT; pH 5.5–7.2, with or without 10 mM MgCl2, in a sodium cacodylate buffer with 0.1–1 M NaCl, KCl or LiCl. All curves are superimposable. (B) For the 5'-AGGGTTAGGGTTAGGGTTAGGG intramolecular DNA quadruplex pH 6.0–7.2, in a sodium cacodylate buffer with 0.1–1 M NaCl (in blue), or 0.03–0.6 M KCl (in red). Note that CD spectra of the same sequence are shown as Supplementary Figure S3, and error bars at each wavelength (average of 11 independent measurements) as Supplementary Figure S4.
The first obvious application of this analysis is to help in the choice of appropriate wavelength(s) for a thermal denaturation experiment for more detailed thermodynamic studies. These results have been collected in Table 1. It should be noted that the relation between a differential peak and a ‘good’ wavelength is not always direct. In some cases, the observed hypochromicity does not result from the melting of a structure, but is the consequence of temperature-dependent extinction coefficient of the folded and/or the unfolded forms. For example, the positive peak at 273 nm observed in most quadruplexes (Figure 2I) results from a significant temperature-dependent increase of absorbance of the single strands. This is related to the classical problem of the sloping baseline analysis of melting profiles: all species exhibit temperature-dependent absorbance properties. One would then argue that the most convenient differential spectra should be isothermal; however, this requires the ability to record the absorbance of the folded and the unfolded states at the same temperature. This is possible in only a limited number of cases (21), when association and dissociation of the structure are very slow. Otherwise, T-jump or fast mixing experiments are required. Despite this limitation, this simple dual-temperature measurement is sufficient for most determinations.
Table 1 Useful wavelengths for each nucleic acid structure
These spectra provide information complementary to circular dichroism. The UV absorbance of nucleic acids from 200 to 300 nm is exclusively due to transitions of the planar purine and pyrimidine bases. The phosphate backbone begins to contribute at 190 nm. The CD spectrum is probing the asymmetry of bases. At wavelengths >200 nm, the asymmetry arises from CD induced into the bases as a result of stacking interactions and coupling with backbone transitions. Hypochromicity is mainly the result of base stacking (22). Therefore, one can expect that the relative positioning of the bases will play a crucial role in this differential spectra. The theory of hypochromicity is complex, the changes in light absorption occurring upon denaturation is attributed to the interaction between the dipoles induced in the chromophore by the light (22). It is striking that a number of very different structures involve a negative differential absorbance at 295 nm (Table 1). This might be the result of the significant contribution of n * transition moments at this wavelength: as these transition moments are parallel to the helix axis (23), increased stacking resulting from the formation of a folded structure may result in hyperchromism. However, the absorbance spectra of nucleic acid bases are deceptively simple: they actually correspond to a number of different * and n * transitions. Hypochromism at a given wavelength will not only depend on the intrinsic transition moments of each base, but also on the relative moment of the interacting bases (base-paired or stacked). Such complexity will make an accurate calculation of predicted differential absorbances a difficult task.
Although the dissociation of all G-quadruplexes analyzed so far show an inverted transition at 295 nm (13), this property is not unambiguously indicative of G-quadruplexes. Z-DNA (Figure 2C) i-motif (Figure 2G) Hoogsteen duplexes (Figure 2F) and pyrimidine triplex (Figure 2H) formation also leads to an inverted transition at 295 nm. In Figure 2F–H, this is the result of the differences in absorbance between C and C+ (20). A simple comparison of these spectra indicates that the global shape rather than the wavelength of the peak is specific for each structure (Figure 2).
These observations suggest that the simple recording of two absorbance spectra might be sufficient to obtain qualitative structural information on the identity of the folded form. This was a common accepted practice in the sixties, when the melting of polynucleotides was analyzed . Unfortunately, with a few notable exceptions (27–32), complete spectra are rarely presented in contemporary reports. We demonstrate here that in many cases this difference is very informative. Of course, this requires that the high and low temperatures correspond to the unfolded and folded forms, respectively. In other words, the Tm of the structure must lie between these two values. Although the spectra of different nucleic acids that fold into the same type of structure are not identical, significant homology was found between the spectra. The uniqueness of these spectra is sufficient, in most cases, to suggest a type of folding. However, we could not confidently distinguish the differential spectra of a purine-rich triplex (the so-called GA or GT motif) from the melting of a duplex (data not shown). This is in agreement with the conclusions that the formation of these triplexes occurs without significant hypochromicity (33). Similarly, no significant difference was found between a mismatched and a fully complementary DNA duplex (data not shown).
We have now collected >900 spectra, corresponding to >200 different nucleic acid sequences (mostly DNA), each tested 2–26 times under different experimental conditions. This is, to our knowledge, the first attempt to systematically collect thermal difference spectra as initially proposed 40 years ago by Fresco et al (34). Our goal is to first build an open access database of the spectra of nucleic acids and to propose a predictive algorithm that provides a probability function of a given structure. This should be helpful in the discovery of new nucleic acid folding schemes, and in the crude determination of the global folding of an oligonucleotide (which, for example, might have been found by SELEX). We are now applying this method to the analysis of trinucleotide repeats (35) and to samples which may adopt complex conformations, involving more than one simple structure. Efforts are also being made on RNA structures (and corresponding DNA strands carrying deoxyuracyl bases): RNA duplex and quadruplex structures exhibit related but distinct signatures from the corresponding DNA duplexes (data not shown) and quadruplexes (18,21), respectively. We also wish to analyze the signature of modified bases (such as 5-methyl cytosines) and the impact of modifying the sugar/phosphate backbone (such as 2'-O-methyl phosphorothioates, LNA or PNA).
SUPPLEMENTARY DATA
Supplementary Data are available at NAR Online.
ACKNOWLEDGEMENTS
We thank M. Rougée, T. Garestier, L. Guittat, P. Arimondo, P. Alberti, B. Saccà (MNHN, Paris, France) M. Mills (UCT, Cape Town, Republic of South Africa) for helpful discussions. This work was supported by an ARC grant (#3365 to J.L.M), a EU FP6 MOLCANCERMED grant (to J.L.M.), a fondation Jérome Lejeune research grant (to S.A.) and National Cancer Institute grant CA35635 (to J.B.C). Funding to pay the Open Access publication charges for this article was provided by INSERM.
Conflict of interest statement: None declared.
REFERENCES
Rich, A. (1993) DNA comes in many forms Gene, 135, 99–109 .
Rich, A. and Zhang, S.G. (2003) Z-DNA: the long road to biological function Nature Rev. Genet., 4, 566–572 .
Jovin, T.M., Rippe, K., Ramsing, N.B., Klement, R., Elhorst, W., Vojtisova, M. (1990) Parallel stranded DNA In Sarma, R.H. and Sarma, M.H. (Eds.). Structure & Methods, DNA & RNA, Schenectady, New York Adenine Press Vol 3, pp. 155–174 .
Moser, H.E. and Dervan, P.B. (1987) Sequence specific cleavage of double helical DNA by triple helix formation Science, 238, 645–650 .
Le Doan, T., Perrouault, L., Praseuth, D., Habhoub, N., Decout, J.-L., Thuong, N.T., Lhomme, J., Hélène, C. (1987) Sequence specific recognition, photocrosslinking and cleavage of the DNA double helix by an oligo thymidylate covalently linked to an azidoproflavine derivative Nucleic Acids Res., 15, 7749–7760 .
Gellert, M., Lipsett, M.N., Davies, D.R. (1962) Helix formation by guanylic acid Proc. Natl Acad. Sci. USA, 48, 2013–2018 .
Gehring, K., Leroy, J.L., Guéron, M. (1993) A tetrameric structure with protonated cytosine.cytosine base pairs Nature, 363, 561–565 .
Henderson, E., Hardin, C.C., Walk, S.K., Tinoco, I., Jr, Blackburn, E.H. (1987) Telomeric DNA oligonucleotides form novel intramolecular structures containing guanine.guanine base pairs Cell, 51, 899–908 .
Sundquist, W.I. and Klug, A. (1989) Telomeric DNA dimerizes by formation of guanine tetrads between hairpin loops Nature, 342, 825–829 .
Williamson, J.R., Raghuraman, M.K., Cech, T.R. (1989) Monovalent cation induced structure of telomeric DNA: the G-quartet model Cell, 59, 871–880 .
Kang, C.H., Berger, I., Lockshin, C., Radliff, R., Moyzis, R., Rich, A. (1994) Crystal structure of intercalated four stranded d(C3T) at 1.4 ? resolution Proc. Natl Acad. Sci. USA, 91, 11636–11640 .
Parkinson, G.N., Lee, M.P.H., Neidle, S. (2002) Crystal structure of parallel quadruplexes from human telomeric DNA Nature, 417, 876–880 .
Mergny, J.L., Phan, A.T., Lacroix, L. (1998) Following G-quartet formation by UV-spectroscopy FEBS Lett., 435, 74–78 .
Cantor, C.R., Warshaw, M.M., Shapiro, H. (1970) Oligonucleotide interactions. 3. Circular dichroism studies of the conformation of deoxyoligonucleotides Biopolymers, 9, 1059–1077 .
Plum, G.E. (2000) Determination of oligonucleotide molar extinction coefficients Current Protocols in Nucleic Acid Chemistry, New York Wiley pp. 7.3.1–7.3.17 .
Mergny, J.L. and Lacroix, L. (2003) Analysis of thermal melting curves Oligonucleotides, 13, 515–537 .
Wang, Y. and Patel, D.J. (1993) Solution structure of the human telomeric repeat d G-tetraplex Structure, 1, 263–282 .
Saccà, B., Lacroix, L., Mergny, J.L. (2005) The effect of chemical modifications on the thermal stability of different G-quadruplexes-forming oligonucleotides Nucleic Acids Res., 33, 1182–1192 .
Li, J., Correia, J.J., Wang, L., Trent, J.O., Chaires, J.B. (2005) Not so crystal clear: the structure of the human telomeric G-quadruplex in solution differs from that present in a crystal Nucleic Acids Res., 33, 4649–4659 .
Mergny, J.L., Lacroix, L., Han, X., Leroy, J.L., Hélène, C. (1995) Intramolecular folding of pyrimidine oligodeoxynucleotides into an i-DNA motif J. Am. Chem. Soc., 117, 8887–8898 .
Mergny, J.L., de Cian, A., Ghelab, A., Saccà, B., Lacroix, L. (2005) Kinetics of tetramolecular quadruplexes Nucleic Acids Res., 33, 81–94 .
Tinoco, I., Jr. (1960) Hypochromism in polynucleotides J. Am. Chem. Soc., 82, 4785–4790 .
Rich, A. and Kasha, M. (1960) The n->pi* transition in nucleic acids and polynucleotides J. Am. Chem. Soc., 82, 6197–6199 .
Riley, M., Maling, B., Chamberlin, M.J. (1966) Physical and chemical characterization of two and three stranded adenine thymine and adenine uracil homopolymer complexes J. Mol. Biol., 20, 359–389 .
Blake, R.D., Massoulié, J., Fresco, J.R. (1967) A spectral approach to the equilibria between polyriboadenylate and polyribouridylate and their complexes J. Mol. Biol., 30, 291–308 .
Martin, F.H., Uhlenbeck, O.C., Doty, P. (1971) Self complementary oligoribonucleotides: adenylic acid–uridylic acid block copolylers J. Mol. Biol., 57, 201–215 .
Ramsing, N.B. and Jovin, T.M. (1988) Parallel stranded duplex DNA Nucleic Acids Res., 16, 6659–6676 .
van de Sande, J.H. and Jovin, T.M. (1982) Z* DNA, the left handed helical form of poly d(G–C) in MgCl2-ethanol is biologically active EMBO J., 1, 115–120 .
Lavelle, L. and Fresco, J.R. (2003) Enhanced stabilization of the triplexes d(C+-T)(6):d(A-G)(6).d(C-T)(6), d(T)(21):d(A)21.d(T)(21) and poly r(U:A.U) by water structure-making solutes Biophys. Chem., 105, 701–720 .
Davis, T.M., McFail-Isom, L., Keane, E., Williams, L.D. (1998) Melting of a DNA hairpin without hyperchromism Biochemistry, 37, 6975–6978 .
Kankia, B.I. (2003) Mg2+-induced triplex formation of an equimolar mixture of poly(rA) and poly(rU) Nucleic Acids Res., 31, 5101–5107 .
Haq, I., Chowdhry, B.Z., Chaires, J.B. (1997) Singular value decomposition of 3-D DNA melting curves reveals complexity in the melting process Eur. Biophys. J., 26, 419–426 .
Mills, M., Arimondo, P., Lacroix, L., Garestier, T., Hélène, C., Klump, H.H., Mergny, J.L. (1999) Energetics of strand displacement reactions in triple helices: a spectroscopic study J. Mol. Biol., 291, 1035–1054 .
Fresco, J.R., Klotz, L.C., Richards, E.G. (1963) A new spectroscopic approach to the determination of helical secondary structure in ribonucleic acids Cold Spring Harbor Symp. Quant. Biol., 28, 83–90 .
Amrane, S., Saccà, B., Mills, M., Chauhan, M., Klump, H.H., Mergny, J.L. (2005) Length-dependent energetics of (CTG)n and (CAG)n trinucleotide repeats Nucleic Acids Res., 33, 4065–4077 .
Jose, D. and P?rschke, D. (2004) Dynamics of the B-A transition of DNA double-helices Nucleic Acids Res., 32, 2251–2258 .
Pohl, F.M. and Jovin, T.M. (1972) Salt-induced co-operative conformational change of a synthetic DNA: equilibrium and kinetic studies with poly(dG-dC) J. Mol. Biol., 67, 375–396 .(Jean-Louis Mergny*, Jing Li1, Laurent La)