Non-Covalent Interactions, Conjugation, and DNA Duplex Stability

Chemical Nature of Nucleic Acids and Intermolecular Interactions

• Electronegativity and Bond Polarity:

  • Polar covalent bonds form between atoms possessing significantly different electronegativity values.

  • Relative electronegativity scale:     $P \approx H < C \approx S < N < O$

  • The vast majority of covalent bonds in nucleic acid backbones (e.g., $\text{P}=\text{O}$, $\text{P}-\text{O}$, $\text{C}-\text{O}$, $\text{C}-\text{N}$) are strongly polar, giving nucleic acids exceptional aqueous solubility and high hydrophilicity.

  • Although nucleobases contain aromatic rings, their constituent polar covalent bonds prevent them from being nonpolar or apolar.

• Phosphodiester Backbone Charges and Electrostatic Repulsion:

  • The phosphodiester group in nucleic acid backbones has a $pK_a \approx 1$.

  • At physiological pH (near $7.0$), phosphodiester groups are essentially $100\%$ deprotonated, carrying a full formal negative charge ($-1$) per phosphate linkage.

  • In a double-stranded DNA helix, the two opposed sugar-phosphate backbones exert strong electrostatic repulsion against each other (negative charge vs. negative charge).

  • The two strands are held together in a stable duplex exclusively by attractive interactions between opposing nucleobases.

Non-Covalent Forces Governing DNA Duplex Stability

• Summary of Intermolecular Forces in DNA:

  • Ionic Interactions: Repulsive electrostatic forces between negatively charged backbone phosphodiester groups.

  • Hydrogen Bonds ($H$-bonds): Specific, complementary interactions occurring between opposing base pairs on opposite strands.

  • van der Waals Interactions (Three Distinct Types):

  1. London Dispersion Attraction: Attractive forces between transient dipoles of closely packed atoms along the duplex.

  2. Aromatic $\pi$-Stacking: Attractive interactions between overlapping $\pi$-electron clouds of vertically stacked nucleobases.

  3. van der Waals Repulsion: Steric repulsion that occurs when electron clouds of adjacent atoms overlap too closely.

  • Note: The hydrophobic effect plays a minimal role in nucleic acid duplex stability compared to its dominant role in protein folding.

• Directionality and Geometry of Hydrogen Bonds:

  • Hydrogen bonding is a strongly direction-dependent interaction.

  • An $H$-bond achieves maximum strength when the donor atom, the shared hydrogen atom, and the acceptor atom lie in a straight line ($180^\circ$).

  • Deviations from $180^\circ$ reduce hydrogen bond strength.

  • In contrast, ionic electrostatic interactions are non-directional and depend solely on inter-charge distance.

  • Because $\pi$-conjugation forces nucleobases to remain planar, all $H$-bonding functional groups reside in a single plane, forming optimal, linear $180^\circ$ hydrogen bonds across complementary strands.

  • Guanine-Cytosine (G:C) base pairs form $3$ inter-strand hydrogen bonds, whereas Adenine-Thymine (A:T) base pairs form $2$.

• Distance Dependence of van der Waals Interactions:

  • Any pair of atoms experiences opposing van der Waals forces:

  • At short distances, Pauli overlap causes sharp van der Waals repulsion.

  • At larger distances, weak London dispersion attraction dominates.

  • The interatomic distance $r$ corresponding to the potential energy minimum defines the van der Waals radius (the distance depicted in space-filling molecular models).

  • The strength of van der Waals interactions decays far more steeply with increasing distance than the strength of long-range ionic interactions.

• London Dispersion Forces in Macromolecules:

  • London attraction stems from temporary, transient dipoles caused by brief asymmetries in electron distribution around atomic nuclei.

  • Individual atomic London forces are non-directional and weak ($<4\,kJ\,mol^{-1}$).

  • In large, tightly packed biological macromolecules, thousands of simultaneous London attractions accumulate to yield substantial total thermodynamic stability.

• Planar Base Packing Maximizes Dispersion and Stacking:

  • Planar nucleobases stack on top of one another along the central axis of the DNA duplex.

  • This parallel alignment maximizes contact surface area and minimizes interatomic distance, optimizing London dispersion forces.

  • G:C base pairs contain more co-planar heavy atoms than A:T base pairs, generating stronger overall London attractions.

  • Aromaticity permits overlapping $\pi$-electron cloud interactions ($\pi$-stacking) between stacked base pairs, which are also intrinsically stronger in G:C pairs than in A:T pairs.
    
    

• Melting Temperature ($T_m$): The specific temperature at which $50\%$ of the DNA molecules in a sample are denatured into single strands.

Dissecting $T_m$ Determinants: Base Stacking vs. Hydrogen Bonding

• Effect of Ionic Strength (Charge Screening):

  • Increasing cation concentration (e.g., $50\,mM \rightarrow 75\,mM \rightarrow 100\,mM\text{ NaCl}$) shifts the DNA melting curve toward higher temperatures, increasing $T_m$.

  • Mechanism: Positively charged sodium ions ($Na^+$) cluster around the negatively charged phosphodiester backbone, screening the negative charges and suppressing inter-strand electrostatic repulsion.

• Effect of Duplex Length and Mismatches:

  • Duplex Length: Longer DNA duplexes ($20\,bp \rightarrow 40\,bp \rightarrow 60\,bp$) exhibit higher $T_m$ values. Adding a base pair introduces two additional repulsive backbone phosphate groups, but the attractive base-pairing and base-stacking forces outweigh phosphate repulsion, increasing overall duplex stability.

  • Mismatches: Introducing a base mutation disrupts complementary hydrogen bonding and structural packing, lowering $T_m$.

• Evaluating the G:C Content Stability Hypothesis:

  • Initial Hypothesis: G:C-rich duplexes are more stable because G:C pairs form $3$ hydrogen bonds while A:T pairs form only $2$.

  • Experimental Test: Compare two DNA duplexes engineered to contain the exact same total number of hydrogen bonds ($48\text{ H-bonds}$):

  • Duplex 1: $24\,bp$ duplex, $100\%\text{ A:T}$ content ($48\text{ total H-bonds}$) $\rightarrow T_m = 37\,^{\circ}C$

  • Duplex 2: $16\,bp$ duplex, $100\%\text{ G:C}$ content ($48\text{ total H-bonds}$) $\rightarrow T_m = 78\,^{\circ}C$

  • Conclusion: Despite having identical total numbers of hydrogen bonds, the $16\,bp$ G:C duplex is dramatically more stable ($T_m = 78\,^{\circ}C$) than the $24\,bp$ A:T duplex ($T_m = 37\,^{\circ}C$).

  • Key Insight: The enhanced stability of G:C-rich DNA is NOT primarily due to forming additional hydrogen bonds. Instead, it is predominantly driven by superior London dispersion forces and stronger aromatic $\frac{\pi}{\pi}$-stacking interactions between adjacent G:C base pairs!