|Home | About | Journals | Submit | Contact Us | Français|
The online version of this article has been published under an open access model. Users are entitled to use, reproduce, disseminate, or display the open access version of this article for non-commercial purposes provided that: the original authorship is properly and fully attributed; the Journal and Oxford University Press are attributed as the original place of publication with the correct citation details given; if an article is subsequently reproduced or disseminated not in its entirety but only in part or as a derivative work this must be clearly indicated. For commercial re-use, please contact email@example.com
Two factors are mainly responsible for the stability of the DNA double helix: base pairing between complementary strands and stacking between adjacent bases. By studying DNA molecules with solitary nicks and gaps we measure temperature and salt dependence of the stacking free energy of the DNA double helix. For the first time, DNA stacking parameters are obtained directly (without extrapolation) for temperatures from below room temperature to close to melting temperature. We also obtain DNA stacking parameters for different salt concentrations ranging from 15 to 100 mM Na+. From stacking parameters of individual contacts, we calculate base-stacking contribution to the stability of A•T- and G•C-containing DNA polymers. We find that temperature and salt dependences of the stacking term fully determine the temperature and the salt dependence of DNA stability parameters. For all temperatures and salt concentrations employed in present study, base-stacking is the main stabilizing factor in the DNA double helix. A•T pairing is always destabilizing and G•C pairing contributes almost no stabilization. Base-stacking interaction dominates not only in the duplex overall stability but also significantly contributes into the dependence of the duplex stability on its sequence.
Stability of the DNA double helix with respect to separation of complementary strands is known to depend on the base composition of the duplex (1–4). A classical study of Marmur and Doty (1) on DNA polymer stability gives a linear relationship between the G•C content of the polymer and its melting temperature. The simplest explanation of the linear dependence is that A•T- and G•C-pairs differ in stability independently of their neighbors. Accordingly, base-stacking interactions have been thought to constitute only a small correction to the major effect of the differences in G•C and A•T pair stabilities (5). Nearest-neighbor stability parameters have been introduced to account for sequence effects in DNA stability. These parameters are obtained from the analysis of the melting data for DNA polymers (5), DNA oligomers (6–8) and DNA dumbbells (9,10), and present the cumulative (base pairing and stacking) contribution of each dinucleotide stack to the overall stability of the molecule. In fact, DNA melting experiments do not allow separation of the two contributions.
Partitioning of base pairing and stacking contributions to DNA stability not only delivers a new aspect in the fundamental understanding of DNA structure and energetics, but also it has significant implications in a number of biological processes. Fluctuations in local helical conformation of DNA, the phenomenon known as DNA breathing, lead to infrequent events of base pair opening thus making normally buried groups available for modification and interaction with proteins (11,12). Fluctuational base pair opening implies disruption of hydrogen bonds between the complementary bases and flipping the base out of the helical stack disrupting two contacts. Heterogeneous stacking at these contacts determine sequence dependence of the base pair fluctuational motility. Moreover, single-stranded break (a nick) in the DNA double helix is stabilized by stacking interactions between base pairs flanking the lesion; these interactions are sequence-dependent (13). In the cell, DNA nicks are substrates for DNA damage-detecting and DNA-repair proteins (14–17). In vitro, DNA nick-based approaches are often used in DNA detection and amplification protocols (18,19). Also, stabilization achieved through coaxial stacking interactions has been put into effect to improve the efficiency of short primer hybridization for standard sequencing protocols (20) and within the format of sequencing by hybridization approach (21). By using partially double-stranded probe with a 5 nt overhang, the latter method takes advantage of the enhancement of the stability of the short duplex formed between the probe overhang and single-stranded target by coaxial stacking with the duplex interface of the probe. Another biotechnology application involving nicked DNA is the formation of padlock probes via ligation of a nick formed by two termini of artificial oligonucleotide hybridized to single-stranded or locally opened duplex DNA (22–24).
A number of theoretical studies addressed characterization of stacking and base pairing contributions to the duplex stability (25–30). Experimentally, energetics of base-stacking interactions in nucleic acids has been evaluated by studying the effect of dangling (unpaired) terminal bases on the overall stability of duplexes (31–35) and in the coaxial stacking hybridization experiments where binding of a short oligonucleotide to the single-stranded DNA template is assisted by stacking interaction with the nearby duplex interface (36–41). Both of these approaches rely on thermal denaturation measurements of short duplex molecules.
We have recently introduced an entirely new approach for characterization of stacking interactions in the DNA double helix (13). For this purpose we study DNA molecules with solitary nicks positioned in a strictly defined sequence context. We subject these molecules to PAGE. Structurally, introduction of a break to one of the strands of DNA perturbs only slightly the double helix conformation with stacking and base pairing at the nick site remaining intact (42–49). Nicked molecules have been shown to move somewhat slower during PAGE than intact molecules of the same size; this retardation is enhanced at higher temperatures (14,50,51).
The model of the DNA nick underlying our approach involves an equilibrium between the two conformations as shown in Figure 1a (13). One conformation is very close to that of the intact double helix where stacking between the base pairs flanking the nick is conserved. We assume that no optimization of stacking interactions—like in case of duplexes with dangling nucleotides (52)—occurs at the site of the nick in DNA. The other conformation corresponds to complete loss of stacking at the nick site thus inducing a kink in DNA. The fast equilibration between stacked/straight and unstacked/bent conformations of the nick directly affects the mobility of DNA molecule during PAGE leading to a differential retardation characteristic to a particular dinucleotide carrying the nick, i.e. KL. For each contact, stacking free energy parameter, , is ascribed governing the distribution of the molecules between the two conformations. Using PAGE mobility data of DNA molecules with solitary nicks flanked by all possible combinations of the base pairs we extract stacking parameters of the DNA double helix (13). Using heterogeneous stacking parameters and a well-known Marmur–Doty (1) dependence of DNA stability on its content, we have previously arrived at a revised set of the nearest-neighbor stability parameters which is in full agreement with the DNA melting data (13).
Here we measure, for the first time, temperature dependence and salt dependence of base-stacking contribution to the DNA duplex stability. Contributions of A•T and G•C pairing are estimated from the comparison of A•T- and G•C-containing polymer stability parameters with the stacking terms. We find that throughout the temperature range employed, base-stacking interactions stabilize DNA double helix. Temperature dependence of the base-stacking term fully determines the temperature dependence of the DNA stability parameter. Base pairing term is destabilizing in case of A•T and somewhat stabilizing in case of G•C pairs. Differential contribution of base-stacking in A•T- and G•C-containing contacts is responsible for 50% of the dependence of DNA stability on its G•C content. Salt conditions affect stacking parameters leading to stabilization upon an increase in sodium concentration. In this case, we also find that salt dependence of stacking term determines the salt dependence of the DNA stability term.
The 300 bp-long DNA molecules used in this study are PvuII/PvuII (all enzymes used in this study were purchased from New England Biolabs, Ipswich MA) restriction fragments of pUC19 derivatives described before (13). Briefly, pKL (were K and L are DNA bases) carries recognition sites for nicking endonucleases each of which introduces a single nick to the forward (N.BstNB I) or reverse (N.Alw I) strand of KL/K′L′ dinucleotide stack (or simply KL; K′ and L′ are complementary to K and L) yielding KL/F and KL/R fragments, respectively (Figure 1b). The nick is located 104 bp apart from the terminus of PvuII/PvuII fragment. Molecules with a single 2 nt gap in place of a nick are obtained by digesting pG2F and pG2R by N.BbvC IA followed by N.BstNBI digestion. As shown in Figure 1b, nicking sites of these nicking endonucleases are 2 nt apart and located on the forward strand of G2/F and reverse strand of G2/R. To remove 5′ phosphoryl groups, nicked and gapped molecules were incubated with Alkaline Phosphatase (CIP). Following treatment with enzymes, all DNA samples were purified by phenol extraction, precipitated with ethanol and resuspended in TE buffer [10 mM Tris (pH 7.4) and 1 mM EDTA].
Nicked, gapped and intact DNA fragments 300 bp in length were concurrently subjected to electrophoresis through 6.7% (w/v) polyacrylamide gels run at 17 V/cm with 1× TBE [90 mM Tris (pH 8.0), 90 mM boric acid and 1 mM EDTA] as a running buffer. We have shown before that under this set of PAGE conditions stacked-unstacked equilibrium at the nick site of DNA is not affected by specific experimental design (13). Specifically, Figure 5 of (13) demonstrates that the presence of the gel matrix and force applied to the molecule during PAGE do not affect apparent free energy parameters of stacked-unstacked equilibrium at the DNA nick. We also conducted PAGE in the presence of increasing concentrations of sodium by adding NaCl up to 85 mM to the gel and to the running buffer; the strength of electric field for these experiments was reduced to 12 V/cm. The gel electrophoresis was performed using Protean II xi Cell (BioRad, Hercules, CA) where the gel-plate sandwich is mounted against the temperature-controlling core connected to a circulating water bath set to a desired temperature. The temperature during PAGE was measured by a thermocouple (Oakton, Vernon Hills, IL) with a type T wire probe inserted directly into the gel. DNA band patterns were visualized by ethidium bromide staining and detected by the CCD camera.
Relative mobilities of DNA fragments with a nick, μ, were calculated as a ratio of the distance migrated by the nicked fragments to the distance migrated by the intact fragment during PAGE. All experiments were conducted at least in triplicates. To enhance the separation of the molecules with solitary nicks from the intact molecules we added various concentrations of urea (0.4 to 8 M) to the gel thus allowing subsequent extrapolation of measured parameters to zero urea (13). We also conducted PAGE in the presence of dimethylformamide (DMF) in place of urea. In this case, a smaller range of concentrations was used since the gels became fragile at concentrations of DMF above 1.5 M.
Two synthetic 65 nt-long complementary DNA strands (MWG Biotech Inc., High Point, NC) were annealed in TE buffer appended with a desired concentration of NaCl (from 20 to 200 mM) at 0.5 µM of each strand. The same duplex was also prepared in the 1×TBE buffer. The purity of the duplex was checked by PAGE. Melting curves of DNA duplexes were collected using a CARY Bio100 spectrophotometer equipped with the peltier thermocontroller (Varian Instruments, Walnut Creek, CA). Both denaturation and renaturation curves were collected at slow heating/cooling rate of at 0.08 °/min; no hysteresis between heating and cooling curves has been observed. Melting temperature, TM, was estimated as a temperature at the midpoint of the duplex-to-coil transition.
In melting studies we determined stabilities of four DNA duplexes, which differed in G•C content, χ (DNA sequences are listed in Supplementary Table 1). Salt dependence of A•T-containing duplex stability (χ = 0%) was measured directly for 22 mM < [Na+] < 202 mM (i.e. 20–200 mM NaCl plus 1 mM EDTA). Salt dependence of more stable G•C-containing duplex was obtained by extrapolating TM's measured for four duplexes with χ = 0%, χ = 14%, 32% and 41.5% to χ = 100%. Corresponding plots are shown in Figure 4a.
During gel electrophoresis, a 300 bp-long DNA molecule with a solitary nick migrates slower with respect to the DNA fragment without the lesion (Figure 1c). Retardation of nicked fragments depends on the identity of the base pairs flanking the nick like in case of the molecules with nicks flanked by As and Ts shown in Figure 1c. Differential retardation of nicked DNA is due to specific interactions characteristic to each nicked dinucleotide stack. Quantitatively, equilibrium between stacked/closed and unstacked/open conformations at the nick site is governed by stacking free energy, ΔGST, so that:
where Nclosed and Nopen are occupancies of stacked and unstacked conformations at the DNA nick, respectively, R is the universal gas constant and T is the absolute temperature.
At a given temperature, due to fast equilibration between the two states, mobility of nicked DNA is a weighted average of the mobility of molecules in closed state and mobility of molecules in open state (13). Since stacked conformation of the DNA nick is very close to the DNA molecule without the lesion μclosed is given by the mobility of the intact DNA fragment. The presence of a single-stranded stretch in a DNA molecule with a short gap prevents stacking interactions between the base pairs flanking the gap thus forcing the molecule into the open state. We use the mobility of molecules with a single 2 nt-long gap to estimate μopen. Thus stacking parameter ΔGST describing stacked-unstacked equilibrium at a nick site can be calculated directly from the mobility data using Equations 1 and 2 (13):
Since the perturbations brought about by the nick are often minor we introduce urea to the gel in order to enhance the separation of nicked molecules from intact molecules (Figure 2c) (13,51). Data obtained in urea-enhanced PAGE is then extrapolated to zero urea concentration resulting in stacking parameters for each nicked dinucleotide stack in the absence of urea, (Figure 2d). The use of denaturant during PAGE is an imperative part of our experimental approach since it allows characterization of even subtle differences in nicked molecule mobilities. To ensure that the nature of the denaturant does not affect the resulting stacking parameters we also conducted PAGE in the presence of varying concentrations of DMF (Figure 2d). Additionally, for the least stable stacks, we can measure stacking parameters directly at zero urea without extrapolating the mobility data collected in the presence of denaturant; this parameter is also plotted in Figure 2d. DMF-enhanced and urea-enhanced PAGE yield quantitatively similar dependencies of apparent stacking parameters on the concentration of the denaturant. This finding is quite unexpected and is most likely due to a narrow range of DMF concentrations used in PAGE; an increase in DMF concentration could yield more significant differences in apparent values measured in the presence of the two denaturants. We are unable to conduct these experiments since the gels containing more than 1.5 M DMF are too fragile to handle. Extrapolation of the data obtained in DMF-enhanced and urea-enhanced PAGE experiments to zero denaturant concentration is consistent with values measured directly in the absence of denaturant.
Our further experiments involve conducting urea-enhanced PAGE at different temperatures to obtain temperature dependence of DNA stacking parameters. Molecules with solitary nicks flanked by A•T pairs shown in Figure 1c, display a defined dependence of the electrophoretic shift towards the molecules with a 2 nt gap when temperature increases from 22 to 52°C. Relative electrophoretic mobilities of nicked AA/F fragment and gapped G2/F fragment are plotted as a function of temperature in Figure 2a. A 30 degrees increase in temperature leads to gradually decreasing PAGE mobility of the nicked fragment while mobility of the gapped fragment remains largely unaffected. In terms of occupancies of the two states of the nicked site, the shift in mobility is equivalent to the shift in the equilibrium towards the open conformation of the nicked stack. Similar results were obtained for the G•C-containing contacts. More narrow temperature range has been employed in this case since molecules with nicks located in more stable G•C-containing stacks do not resolve well enough from intact molecules at temperatures under 30°C.
Stacking free energies were calculated from the mobility data for all four A•T- and four G•C-containing nicked stacks using Equations 1 and 2. For each stack, molecules with solitary nicks located on the forward strand, KL/F and on the reverse strand, KL/R, were examined. Our further analysis is grounded on the assumption that the stacking free energies of the stacked-unstacked equilibrium at the site of the DNA nick are close to the stacking parameters of the intact DNA double helix (13). In this case we calculate the by averaging (i) the parameters obtained for KL/F and KL/R and (ii) the parameters obtained for the stacks equivalent due to dyad symmetry of the double helix (in our case there are two pairs of stacks: AA/TT and TT/AA pair and GG/CC and CC/GG pair). Individual stacking parameters are listed in Supplementary Table 2 and are plotted as a function of temperature in Figure 3a. Error bars on this plot correspond to the scatter of individual stacking parameters around the mean. As expected, throughout the temperature range employed, values for both A•T- and G•C-containing contacts increase (become less stable) gradually with increasing temperature.
To address the effect of ionic strength on stacked-unstacked equilibrium at the site of DNA nick we conducted urea-enhanced PAGE at different salt concentrations by adding NaCl to the gel and to the running buffer. Salt dependence of relative PAGE mobilities of AA/F nicked fragment and corresponding G2/F gapped fragment is shown in Figure 2b. Unlike in case of increasing temperature (Figure 2a), opposite tendencies are displayed by the two molecules with increasing salt concentration. Mobility of the fragment with a gap decreases. The decrease is likely due to more efficient screening of phosphate charges reducing electrostatic repulsion between the two double-stranded segments on either side of the gap thus making the gap more flexible (53). Similar effect of the salt concentration is expected at the nick site in open conformation. Mobility of the nicked fragment increases corresponding to the salt-assisted shift towards stacked state of the nick.
To facilitate the comparison of DNA stacking parameters we determine experimentally, with the DNA stability parameters we have to account for the difference in salt conditions. Our PAGE-based experiments were conducted in a standard gel electrophoresis buffer 1× TBE. In order to find concentration of NaCl equivalent to PAGE conditions with respect to DNA stability, we conducted melting experiments, i.e. we determined TM of duplex DNA in 1× TBE and in the presence of various concentrations of NaCl (Figure 4). For A•T-containing DNA, melting temperatures under different salt conditions were determined directly. A DNA molecule consisting only of guanines and cytosines is much more stable and melts around or above 100°C depending on the concentration of NaCl. Melting temperatures of G•C-containing DNA were obtained by extrapolating TM's measured for four DNA duplexes with different G•C-contents (χ = 0, 14, 32 and 41.5%) to χ = 100% (Figure 4a). Thermal stabilities of A•T- and G•C-containing DNA duplexes in the presence of NaCl and in 1× TBE are compared in Figure 4b. Melting temperatures of both molecules in 1× TBE correspond to the TM of these molecules in the presence of 15 mM [Na+].
We studied eight A•T-containing and eight G•C-containing nicked dinucleotide stacks. Like in the experiments on the temperature dependence of stacking parameters described above, for each KL contact we examined KL/F and KL/R molecules and the stacking parameter was obtained as their average in case of unique contacts (AT, TA, GC and CG). In case of stacks with dyad symmetry, we also average values of equivalent stacks. Individual stacking parameters are plotted in Figure 3b and listed in Supplementary Table 3. Note that all values are plotted as a function of the ‘effective’ concentration of sodium, i.e. the fist point at [Na+] = 15 mM corresponds to the experiment conducted in 1× TBE, the last point at [Na+] = 100 mM corresponds to 1× TBE appended with 85 mM of NaCl. All values become more negative gradually as concentration of sodium increases.
Site-specific nicks in KL/F and KL/R fragments were introduced enzymatically. Nicking enzymes, like any restriction enzymes, cleave DNA strand producing a 5′ phosphate and a 3′ hydroxyl termini, which translate into two negative charges at the nick site of the double helix. To estimate the effect of the phosphate negative charges on the stacked-unstacked equilibrium at the DNA nick we considered DNA molecules where the phosphate has been removed from the nick site. Dephosphorylated nicks carry no charge.
The molecules with dephosphorylated nicks and gaps were studied by urea-enhanced PAGE in the same manner as the molecules with phosphates. As expected, treatment of DNA fragments with Alkaline Phosphatase, which results in the removal of terminal phosphates as well as dephosphorylation of the gap, had no effect on PAGE mobility of intact and gapped fragments (data not shown). ΔGST values obtained for molecules with nicked stacks before and after dephosphorylation are compared in Figure 5a. A defined trend is observed for both A•T- and G•C-containing contacts—there is a significant ~0.4–0.6 kcal/mol shift towards more negative energies upon removal of a phosphate from the nick site only if phosphorylated 5′ nt is a purine nucleotide, like in case of AA/F and TT/R, but not in case of AA/R or TT/F. When the phosphate is located on a pyrimidine nucleotide, its removal has no effect on ΔGST values. This trend holds for 15 out of 16 nicked fragments shown in Figure 5a; the only exception is GC/F. The origin of such diverse behavior cannot be addressed by our approach and it deserves a separate detailed study. We can only suppose that removal of the phosphate from purine nucleotides allows rearrangement—in a sequence-independent manner—of the local structure benefiting stacked conformation. Alternatively, differential effect of phosphate positioning on purines or pyrimidines might be entropic in origin due to larger conformational flexibility of purine nucleotides without 5′ phosphates (54).
Note that ΔGST values of nicked stacks with and without phosphate display similar dependences on salt concentration (Figure 5b). For nicks flanked by A•T pairs, this is the case for AA/R, TT/F, AT/F and AT/R fragments which display no change in ΔGST values upon removal of the phosphate from 5′-T nt and for AA/F, TT/R, TA/F and TA/R fragments which shift towards more negative ΔGST values upon dephosphorylation of 5′-A nt. This result indicates that the phosphate charge at the nick does not contribute to the salt dependence of the measured stacking component; differences seen in Figure 5a between the nicked stacks with or without phosphate on 5′-purine nucleotides are not due to electrostatic interactions. In our further analysis we consider only molecules carrying phosphates at the nicks.
We obtain stacking free energy parameters by studying the equilibrium between stacked and unstacked form of DNA nick schematically presented in Figure 1a. The DNA molecule in the stacked state is very close the molecule without the lesion (42–49). Strictly speaking, the conformation of the unstacked state is unknown and our analysis is based on the assumption that unstacked state is close to the molecule with a short gap in place of a nick. Indeed, presence of a single-stranded gap precludes stacking interactions between two helical interfaces. Note, that this is not the case for 1 nt gaps—it appears that stacking between the base pairs flanking the gap is restored to some degree leading to anisotropic, directional bending of the molecule reducing the size of the gapped cavity (48,55). Molecules with longer gaps however, have been shown to possess isotropic bending flexibility which manifests itself in the absence of helical periodicity in electrophoretic mobility and cyclization kinetics measurements (50,53). Molecules with gaps 2, 3 and 4 nt in length migrate very closely during PAGE revealing similarity of their effective conformations (13,50). Additional factors [i.e. sequence of single-stranded linker (56)] come into effect once the gap size is longer than persistent length of single chains—in this case, the gap is likely to act as a hinge.
Thus, we use molecules with 2 nt-long gaps positioned in forward or reverse strand to approximate the open state of the nick; the site of the gap is thought of as a point of increased isotropic bending flexibility, which leads to retardation of the molecule during PAGE (57). Perturbations of DNA fragment at the gap sites are largely unaffected by ambient conditions—temperature, salt and denaturant concentration (Figure 2a–c).
Nicked DNA undergoes transition between the stacked (closed or intact-like) and unstacked (open or gapped-like) states (13). Note, that like in case of short gaps, perturbation at the nick site has been shown to be isotropic with no preferential directionality (50,58). The partitioning of molecules between these conformations determines the electrophoretic mobility of the nicked fragment and is the basis for measuring stacking parameters using Equations 1 and 2.
For further analysis we assume that the parameters describing the equilibrium at the nick site of DNA that we measure are in fact the same stacking parameters that stabilize intact double helix and determine the first term on the right-hand side of Equation 3:
where ΔGKL is the total melting free energy parameter (per 1 base pair) of a long DNA molecule (neglecting the end effects) and ΔGBP terms are the free energy contributions due to pairing between complementary bases (5). There are two ΔGBP parameters and for A•T and G•C pairs, respectively.
How accurate is the above assumption? Experimentally, the assumption proves to stand a critical test when we compare the values obtained for KL/F and KL/R series: these values indeed are very close [see Figure 4b in ref. (13)]. The values for dinucleotide stacks equivalent due to dyad axis of symmetry of the double helix (e.g. AA/TT and TT/AA) are close as well.
Let us examine the stacked-unstacked equilibrium at the nick of the dinucleotide stack and compare it to the helix-to-coil transition of a DNA doublet. The schematic in Figure 6 depicts these two transitions. Hydrogen bonding between complementary bases and stacking of pairs along the helical axis are disrupted upon melting; duplex melting to single strands is accompanied by the gain in conformational entropy, release of counterions and change in interactions with the solvent. Strictly speaking, partitioning of these effects between ΔGST and ΔGBP terms is not known. For example, base pairing terms in Equation 3 include breakage of hydrogen bonds between bases and formation of hydrogen bonds with water molecules rearranging the water structure around DNA chains. Likewise, base-stacking interactions have electrostatic and hydrophobic components.
The stacked-unstacked transition of a nicked DNA doublet involves loss of stacking interactions while hydrogen bonding between the complementary bases is preserved (Figure 6). It also involves removal of structural constraints imposed in the double helix making most of the conformational space of the intact sugar-phosphate backbone accessible. This part of backbone conformational entropy gain is close to the gain per one strand upon duplex melting. The conformation about the glycosidic bond in an open stack, however, is preserved as in duplex form. Nicked strand is lacking O3′-P linkage; the absence of this constraint allows for additional conformations when second-order interactions (interactions depending on two torsional angles) are repealed (59). Note that the backbone conformational entropy change is independent of the particular sequence and does not affect heterogeneity of values. We assume that conformation of base pairs flanking the nick is preserved as in stacked duplex. Like in case of helix-to-coil transition, other factors contributing to transition parameters come from interactions with solvent and counterions.
We assign free energy difference to stacked-unstacked equilibrium at the nick site of DNA—this parameter describes interaction of two neighboring base pairs and is sequence-dependent. As noted above, it also includes a part of conformational entropy change and other relevant interactions. These factors are universal for all neighboring pairs and do not contribute to heterogeneity of measured . We assume further that this parameter represents the stacking term of intact DNA duplex in Equation 3. Note that this assumption explicitly implies that stacked/closed state of the nicked stack is equivalent to the stack in duplex without the lesion and no optimization of the base pair stacking occurs upon introduction of the nick like in case of dangling unpaired nucleotides (52). If this assumption is valid, then parameters of stacked-unstacked equilibrium at the nicks positioned on one or the other strand of the each dinucleotide stack will be close. Indeed, we have shown previously that stacking parameters for KL/F and KL/R molecules are similar (13). Parameters obtained for stacks equivalent due to dyad axis of symmetry, i.e. AA and TT, GG and CC, are also similar. Note that marked difference is observed in case of coaxial stacking parameters determined in melting experiments (39).
Let us consider a DNA polymer comprised of adenines and thymines. The set of nearest-neighbor stability parameters of all A•T-containing dinucleotides, i.e. ΔGAT, ΔGAA = ΔGTT, ΔGTA, fully describes the contacts found in this polymer. According to Equation 3, stability parameter per 1 bp, ΔGA•T, of the A•T-polymer with a random sequence is given by [see ref. (13)]:
Stacking contribution to the polymer stability is the first term of the right-hand side of Equation 4; the second term presents the base pairing contribution. Experimentally determined stacking term is the mean of stacking parameters we obtained for all A•T-containing contacts. Stability parameter for G•C-containing polymer, ΔGG•C, is calculated analogously from stacking, and base pairing contributions:
DNA stability parameters in terms of melting temperatures, and , are salt-dependent and are given by the following empirical equations for A•T- and G•C-containing polymers, respectively (60):
Stability parameters of A•T-containing and G•C-containing polymers calculated for [Na+] = 15 mM using Equations 6 and 7 are plotted in Figure 7a. Salt dependences of these parameters calculated at 37°C are presented in Figure 8a.
Stability of a DNA polymer comprised only of A's and T's arranged in a random order is given by two terms—the stacking term and the base pairing term—that appear on the right-hand side of Equation 4. Stacking term is calculated as an average of stacking parameters of A•T-containing contacts, which we determine in PAGE experiments conducted under certain set of ambient conditions (temperature and salt concentration). Stacking term of the G•C-containing polymer is obtained analogously.
Stacking contributions to the stability of A•T- and G•C-containing polymers for different temperatures are plotted in Figure 7a. Comparison of the temperature dependence of stacking terms with the temperature dependence of stability terms for both polymers reveals remarkable similarity between the two. Linear regression analysis of the temperature dependences of ΔGST values gives mean slope of kcal/molK, which is very close to the slope of temperature dependence of stability parameters described by Equations 6 and 7. It is evident that temperature dependence of stacking component fully determines the temperature dependence of DNA stability parameters throughout the temperature range employed in our experiments.
According to Equations 4 and 5, the difference between the melting stability of the polymer and the stacking term gives the contribution of A•T or G•C base pairing. The plots shown in Figure 7b present the base pairing term calculated for different temperatures by subtracting the experimentally determined stacking term from the polymer stability parameter. With good accuracy, the base pairing term is independent of temperature within the temperature range used in our experiments. The horizontal lines give the mean values of base pairing parameters at 15 mM Na+ of and . We conclude that stability of the DNA polymer with respect to strand separation is mainly determined by stacking interactions with base pairing being destabilizing (A•T pairs) or contributing almost nothing (G•C-pairs) (Table 1).
Furthermore, throughout the temperature range used in our experiments, heterogeneity of base-stacking interactions of A•T- and G•C-containing contacts, , is responsible for a substantial part of the difference in the stabilities of A•T- and G•C-polymers, (Figure 7a). In the presence of 15 mM Na+, heterogeneous stacking accounts for δΔGST = 0.6 kcal/mol out of δΔG = 1.2 kcal/mol stability difference of A•T- and G•C-containing polymers. Thus, the dependence of DNA stability on the G•C content is not only due to the fact that G•C pairs are stronger than A•T pairs (Figure 7b) but also is due to the stronger stacking component of G•C-containing contacts (Figure 7a). This understanding of stacking domination constitutes a paradigm shift in the view on DNA stability.
Stacking terms of A•T- and G•C-containing polymers are calculated from stacking parameters obtained in PAGE experiments that have been carried out under different salt conditions (Figure 8a). Salt dependences of the stability parameters of both polymers at 37°C calculated using Equation 6 and 7, are also shown in Figure 8a. The two dependencies are very similar indicating that salt dependence of stacking term determines salt dependence of DNA stability parameters throughout the range of salt concentration used in our experiments. In case of A•T-containing contacts, salt dependence of stacking component reproduces salt dependence of DNA stability parameters with a great degree of accuracy. Salt dependences of stacking terms of A•T- and G•C-containing contacts are similar to . The dependence of stacking parameters on salt concentration was unexpected. In fact, in our original publication we assumed that base pairing term is fully responsible for ionic strength dependence of DNA melting stability (13). Figure 8a shows that this assumption is incorrect and base pairing term plotted in Figure 8b is virtually independent of salt concentration. Horizontal lines correspond to the values of base pairing parameters averaged over the whole range of salt concentrations: and (Table 1).
What is the origin of the ionic strength dependence of stacking parameters measured in PAGE experiments? It is clear that stacking is a net phenomenon involving hydrophobic, electrostatic and dispersion components (28,29,62–65). There is no agreement between researchers, however, what is the dominant force in this interaction. Luo et al. (29) emphasize a dominant role of nonelectrostatic interactions. Gellman and co-workers (65,66) report that neither solvophobic effects nor dispersion forces are important in aromatic stacking; stabilization is achieved through attractive interactions between positive and negative partial charges on bases. Guckian et al. found that hydrophobic effects dominate stacking but dispersion forces and electrostatic interaction between partial charges contribute substantially to the total base-stacking as well. This is particularly true for natural bases, which have significant charge localization (28,67). In either case, electrostatic component of base-stacking interaction (substantial or not) will depend on the ambient ionic strength through local screening of partial charges at the solvent accessible positions. In the context of the DNA double helix, however, this effect is likely to be overwhelmed by the salt-dependent interactions between highly negatively charged DNA backbones as described by the polyelectrolyte theory of DNA stability (68,69). Similarity between salt dependences of DNA stacking parameters measured in PAGE experiments and DNA stability parameters indicates that electrostatic component contributes exclusively to of Equation 3. Note, that charge of the phosphate at the 5′ nt of the nicked stack does not affect salt dependence of measured parameters.
Separation of the two contributions to thermal stability of the DNA double helix is achieved by studying PAGE of DNA molecules with solitary nicks and gaps. For the first time, the dependence of DNA stacking parameters on ambient conditions (temperature and salt concentration) is determined. Throughout the temperature and salt concentration range of our experiments, base-stacking interactions are always stabilizing for both A•T- and G•C-containing contacts in the DNA double helix. In fact, DNA stability is mainly determined by base-stacking interactions. G•C pairing does not contribute to stabilization of DNA duplex, while A•T pairing is always destabilizing. This finding presents a paradigm shift in the understanding of the interplay of the forces stabilizing DNA double helix. For all temperatures heterogeneity of stacking interactions in A•T- and G•C-containing contacts accounts for at least half of heterogeneity in the stability of A•T- and G•C-polymers; the other half is due to the difference in the energetics of A•T and G•C base pairing. The data on separation of stacking and base pairing contributions have made it possible to describe sequence-dependent fluctuational opening of the DNA double helix (70).
Supplementary Data are available at NAR Online.
The authors thank Dr R.M. Georgiadis for providing access to spectrophotometer and Ms J. Ruemmele for technical assistance with melting experiments. The authors thank Dr Y. Mamasakhlisov and Dr P. E. Nielsen for fruitful discussions. Funding to pay the Open Access publication charges for this article was provided by NIH.
Conflict of interest statement. None declared.