Unraveling Base Stacking Driving Forces in DNA.

Base stacking is a key determinant of nucleic acid structures, but the precise origin of the thermodynamic driving force behind the stacking of nucleobases remains open. The rather mild stacking free energy measured experimentally, roughly a kcal/mol depending on the identity of the bases, is physiologically significant because while base stacking confers stability to the genome in its double helix form, the duplex also has to be unwound in order to be replicated or transcribed. A stacking free energy that is either too high or too low will over- or understabilize the genome, impacting the storage of genetic information and also its retrieval. While the molecular origin of stacking driving force has been attributed to many different sources including dispersion, electrostatics, and solvent hydrogen bonding, here we show via a systematic decomposition of the stacking free energy using large-scale computer simulations that the dominant driving force stabilizing base stacking is nonhydrophobic solvent entropy. Counteracting this is the conformational entropic penalty on the sugar-phosphate backbone against stacking, while solvent hydrogen-bonding, charge-charge interactions, and dispersive forces produce only secondary perturbations. Solvent entropic forces and DNA backbone conformational strains therefore work against each other, leading to a very mild composite stacking free energy in agreement with experiments.