DNA duplex stability

Contents

Introduction

The stability of the DNA double helix depends on a fine balance of interactions including hydrogen bonds between bases, hydrogen bonds between bases and surrounding water molecules, and base-stacking interactions between adjacent bases.

Slight variations in the DNA sequence can have profound implications on the stability of the DNA duplex. For example, mutations in the base sequence that result from errors that occur during DNA replication can result in mismatches that lead to relatively unstable duplexes. This instability is exploited by proofreading enzymes which recognize the mutation and replace it with the correct nucleotide (see Mutagenesis and DNA repair).

To gain an insight into DNA duplex stability, and how it is affected by changes in primary structure, scientists have studied the structure and thermodynamic stability of a variety of DNA duplexes by using a combination of physical methods including X-ray crystallography, ultraviolet (UV) melting and NMR.

Factors influencing DNA duplex stability

DNA duplex stability is determined primarily by hydrogen bonding, but base stacking also plays an important role.

Hydrogen bonding

The heterocyclic bases of single-stranded DNA have polar amido, amidino, guanidino and carbonyl groups that form a complex network of hydrogen bonds with the surrounding water molecules. Some of these bonds must be broken during duplex formation as the inter-base hydrogen bonds are formed. The overall process is one of "hydrogen bond exchange" and the net change in enthalpy upon duplex formation is partly due to ∆H(H-bonds formed) − ∆H(H-bonds broken). For duplexes of any significant length this is an exothermic process at ambient temperature. Not surprisingly the coming together of two large oligomeric molecules is entropically unfavourable (∆S is negative).

Base stacking

Inter-strand hydrogen bonding is clearly important in driving the formation of DNA duplexes, but it is by no means the only contributing factor. The individual bases form strong stacking interactions which are major contributors to duplex stability, as base stacking is much more prevalent in duplexes than in single strands (Figure 1). Base-stacking interactions are hydrophobic and electrostatic in nature, and depend on the aromaticity of the bases and their dipole moments. Base-stacking interactions in nucleic acid duplexes are partly inter-strand and partly intra-strand in nature. However, it is probably more informative to consider base pairs rather than individual bases as discrete units in order to visualize the stabilising effects of base stacking.

Cartoon representation of base stacking in a DNA double helix

Figure 1 | Cartoon representation of base stacking in a DNA double helix

The degree of stabilization afforded by base stacking depends on the DNA sequence. Some combinations of base pairs form more stable interactions than others, so nearest neighbour base-stacking interactions are important determinants of duplex stability.

Base-stacking interactions increase with increasing salt concentration, as high salt concentrations mask the destabilising charge repulsion between the two negatively charged phosphodiester backbones. DNA duplex stability therefore increases with increasing salt concentration. Divalent cations such as Mg2+ are more stabilising than Na+ ions, and some metal ions bind to specific loci on the DNA duplex.

Thermodynamics of DNA duplex formation

Non self-complementary DNA

When two complementary (but not self-complementary) strands of DNA are mixed in solution, an equilibrium is rapidly established between the single strands, Ass and Bss, and the (double-stranded) duplex, Cds:

Ass + Bss ⇌ Cds

Eliminating ∆rG, the change in Gibbs energy for the forward reaction (duplex formation), from the thermodynamic equations

rG = −RT ln Keq

and

rG = ∆rHTrS

gives the van't Hoff equation:

− ln Keq = (∆rH / RT) − (∆rS / R)

This can be rearranged to give

1 / T = − (R / ∆rH) ln Keq + ∆rS / ∆rH)

in which ∆rH is the enthalpy change, ∆rS is the entropy change, R is the gas constant and Keq is the equilibrium constant.

The equilibrium constant Keq can be written as a function of the extent of association, α, and of the initial total concentration of single-stranded DNA, CT (CT = [Ass]initial + [Bss]initial). At equilibrium there are 1−α moles of Ass, 1−α moles of Bss and α moles of Cds, so

Keq = ([Cds]eq / [Ass]eq[Bss]eq) = (xCCT / xACTxBCT) = 2α / [(1−α)2CT]

At the midpoint of the melting curve, α = ½, so

Keq = (1 / ¼CT) = (4 / CT)

therefore

− ln Keq = − ln (4 / CT) = ln (CT / 4)

Self-complementary DNA

For self-complementary DNA, an equilbrium is established in solution between single-stranded and double-stranded DNA:

2 Ass ⇌ Cds

As with non self-complementary DNA, the equilibrium constant Keq can be written as a function of the extent of association, α, and of the initial concentration of single-stranded DNA, CT (CT = [Ass]initial). At equilibrium there are 1−α moles of Ass and α moles of Cds, so

Keq = ([Cds]eq / [Ass]eq2) = (xCCT / xA2CT2) = α / [2(1−α)2CT]

At the midpoint of the melting curve, α = ½, so

Keq = (1 / CT)

therefore

− ln Keq = − ln (1 / CT) = ln (CT)

Experimental determination of DNA duplex stability

UV melting

Heating double stranded nucleic acids causes the strands to unwind (denature) by disrupting the ordered stacking of the bases and breaking hydrogen bonds. The process can be conveniently monitored by an increase in UV absorbance as the double strands unwind to single strands (owing to hypochromicity). The thermal denaturation of double-stranded DNA is progressive and the concerted melting of the whole structure occurs at a well-defined temperature, corresponding to the mid-point of a smooth transition. This temperature is known as the melting temperature (Tm).

The example in Figure 2 shows a melting curve (UV absorption as a function of temperature). Slow heating of double-stranded DNA causes the unwinding of the ordered helical structure into the two single-stranded constituents. This can be seen as a sigmoidal curve of increasing UV absorption. The mid-point, corresponding to the precise melting temperature (Tm) of the duplex, is indicated.

A typical UV melting curve

Figure 2 | A typical UV melting curve

To determine thermodynamic parameters, samples are made up at different concentrations, each containing equal amounts of two complementary sequences of DNA (Ass and Bss) in a salt/buffer solution. On mixing at ambient temperature the two strands of DNA spontaneously form a duplex (Cds) (i.e. α = 1) (this is true for DNA strands of any reasonable length). Each sample is heated in a UV spectrometer from 20 ºC to (say) 80 ºC at a rate of one degree per minute and the absorbance is monitored as shown in Figure 2. The melting temperature, Tm, is the temperature at the maximum absorbance gradient (dA/dT) and corresponds to 50% association (α = ½). The graph (Figure 3) shows the relationship between absorbance and concentration and illustrates the slight increase in Tm with increasing DNA concentration.

UV melting curves and derivatives

Figure 3 | UV melting curves and derivatives

Since α = ½,

1 / Tm = (R / ∆rH) ln (CT / 4) + (∆rS / ∆rH)

A graph of 1/Tm against ln(CT/4) therefore has a gradient R/∆rH and an intercept ∆rS/∆rH (neglecting the temperature-dependence of ∆rS and ∆rH ). The reciprocals of the recorded Tm values are plotted against the ln(CT/4) values and linear regression is used to determine ∆rH and ∆rS. Finally, ∆rG is evaluated. Figure 4 shows such a van't Hoff plot:

Van’t Hoff plot for the duplex dGGCTATCTATCTGC-dGCAGATATATAGCC (containing a C-T mismatch) giving ∆<sub>r</sub><em>H</em>= −453.5 kJ mol-1, ∆<sub>r</sub><em>S</em> = −1295 J K<sup>-1</sup> mol<sup>-1</sup> and ∆<sub>r</sub><em>G</em> = −52.01 kJ mol<sup>-1</sup> (at 298K)

Figure 4 | Van’t Hoff plot for the duplex dGGCTATCTATCTGC-dGCAGATATATAGCC (containing a C-T mismatch) giving ∆rH= −453.5 kJ mol-1, ∆rS = −1295 J K-1 mol-1 and ∆rG = −52.01 kJ mol-1 (at 298K)

Calorimetry

Calorometry is an classically important method for measuring DNA duplex stability, although UV melting has proved more popular owing the available of UV spectrometers with temperature-controlled sample chambers, and the relative ease of recording Tm values. Whereas UV melting provides an indirect method for measuring ∆rH, ∆rS and ∆rG for melting, calorimetry allows more direct measurement of these thermodynamic parameters.

Prediction of DNA duplex stability

The melting temperature (Tm) of a DNA duplex and the effect of salt can be predicted using the formula

Tm = 81.5 + 16.6 × (−1) + 41 × (nG + nC) / (nA + nG + nC + nT) − 675 / (nA + nG + nC + nT)

However, this formula only applies to DNA duplexes containing exclusively Watson-Crick A.T and G.C base pairs under a well-defined set of conditions, and cannot be used to determine the melting temperatures of mismatch-containing duplexes.

Tuning DNA duplex stability

When designing and synthesizing synthetic DNA oligonucleotides, for example for use in antisense therapy or as PCR probes, it is often desirable to either increase or decrease the stability of the duplex formed on binding to complementary DNA.

Nucleic acid analogues

Artificial nucleic acid analogues bind to complementary DNA. Peptide nucleic acid (PNA), a neutral molecule, forms particularly strong and specific interactions with DNA. Stretches of PNA can be incorporated into DNA oligonucleotides to increase duplex stability; but care must be taken in desigining the PNA-DNA junction, and the stablizing effect of the PNA can be offset by the destabilizing impact of the junction. Locked nucleic acid (LNA) and unlocked nucleic acid (UNA) are analogues of RNA that can be easily incorporated into DNA oligonucleotides during solid-phase oligonucleotide synthesis, and respectively increase and decrease duplex stability. LNA and UNA phosphoramidite monomers are not readily available, however, and it is not always easy to predict the effect of incorporating LNA and UNA units into DNA oligonucleotides before they are synthesized.

Modified bases

Modified DNA bases that increase the stability of base pairs (and therefore the duplex as a whole) have been developed. These modified bases can be incorporated into oligonucleotides during solid-phase synthesis and offer a more predictable method of increasing DNA duplex stability.

AP-dC (G-clamp)

AP-dC (G-clamp) is a tricyclic analogue of cytosine. Like cytosine, G-clamp forms base pairs with guanine; but, while G·C base pairs are stabilized by three hydrogen bonds, G·G-clamp base pairs have four hydrogen bonds (Figure 5). The extra hydrogen bond can lead to an increase in Tm of 18 °C. The incorporation of multiple G-clamp bases can result in further increases in duplex stability.

Structure of the G·G-clamp base pair

Figure 5 | Structure of the G·G-clamp base pairG·G-clamp base pairs (right) increase the stability of a DNA duplex by up to 18 °C relative to G·C base pairs (left).

2-Aminoadenine

2-Aminoadenine forms base pairs with thymine that are stabilized by three hydrogen bonds (one more than the two hydrogen bonds of A·T base pairs) (Figure 6). Each substitution of an A·T base pair by a 2-aminoadenine·thymine base pair increases the stability of the duplex by about 3 °C.

Structure of the 2-aminoadenine·thymine base pair

Figure 6 | Structure of the 2-aminoadenine·thymine base pair2-aminoadenine·thymine base pairs (right) increase the stability of a DNA duplex by around 3 °C (per substitution) relative to A·T base pairs (left).

Other modified bases

Modest increases in duplex stability have been observed with other modified bases, including 5-methylcytosine and C(5)-propynylcytosine (replacing cytosine), and C(5)-propynyluracil (replacing thymine).

See also

Our free online Nucleic Acids Book contains information on all aspects of nucleic acids chemistry and biology.