Thermodynamic Parameters for Stacking and Hydrogen Bonding of Nucleic Acid Bases in Aqueous Solution: Ab Initio/Langevin Dipoles Study

At nonzero temperatures, stacked structures of methylated nucleic acid base pairs and microhydrated nonmethylated nucleic acid base pairs are favored over planar hydrogen-bonded structures: a molecular dynamics simulations study

The dynamic structure of all ten possible nucleic acid (NA) base pairs and methylated NA base pairs hydrated by a small number of water molecules (from 1 to 16) was determined by using molecular dynamics simulations in the NVE microcanonical and NVT canonical ensembles with the Cornell force field (W. D. Cornell, P. Cieplak, C. I. Bayly, I. R. Gould, K. M. Merz, D. M. Ferguson, D. C. Spellmeyer, T. Fox, J. E. Caldwell, P. Kollman, J. Am. Chem. Soc. 1995, 117, 5179). The presence of one water molecule does not affect the structure of any hydrogen-bonded (H-bonded) nonmethylated base pair. An equal population of H-bonded and stacked structures of adenine...adenine, adenine...guanine and adenine... thymine pairs is reached if as few as two water molecules are present, while obtaining equal populations of these structures in the case of adenine...cytosine, cytosine...thymine, guanine... guanine and guanine...thymine required the presence of four water molecules, and in the case of guanine...cytosine, six. A comparable population of planar, H-bonded and stacked structures for cytosine...cytosine and thymine... thymine base pairs was only obtained if at least eight water molecules hydrated a pair. Methylation of bases changed the situation dramatically and stacked structures were favoured over H-bonded ones even in the absence of water molecules in most cases. Only in the case of methyl cytosine...methyl cytosine, methyl guanine...methyl guanine and methyl guanine...methyl cytosine pairs were two, two or six water molecules, respectively, needed in order to obtain a comparable population of planar, H-bonded and stacked structures. We believe that these results give clear evidence that the preferred stacked structure of NA base pairs in the microhydrated environment, and also apparently in a regular solvent, is due to the hydrophilic interaction of a small number of water molecules. In the case of methylated bases, it is also due to the fact that the hydrogen atoms most suitable for the formation of H-bonds have been replaced by a methyl group. A preferred stacked structure is, thus, not due to a hydrophobic interaction between a large bulk of water molecules and the base pair, as believed.

The physical basis of nucleic acid base stacking in water

It has been argued that the stacking of adenyl groups in water must be driven primarily by electrostatic interactions, based upon NMR data showing stacking for two adenyl groups joined by a 3-atom linker but not for two naphthyl groups joined by the same linker. In contrast, theoretical work has suggested that adenine stacking is driven primarily by nonelectrostatic forces, and that electrostatic interactions actually produce a net repulsion between adenines stacking in water. The present study provides evidence that the experimental data for the 3-atom-linked bis-adenyl and bis-naphthyl compounds are consistent with the theory indicating that nonelectrostatic interactions drive adenine stacking. First, a theoretical conformational analysis is found to reproduce the observed ranking of the stacking tendencies of the compounds studied experimentally. A geometric analysis identifies two possible reasons, other than stronger electrostatic interactions, why the 3-atom-linked bis-adenyl compounds should stack more than the bis-naphthyl compounds. First, stacked naphthyl groups tend to lie further apart than stacked adenyl groups, based upon both quantum calculations and crystal structures. This may prevent the bis-naphthyl compound from stacking as extensively as the bis-adenyl compound. Second, geometric analysis shows that more stacked conformations are sterically accessible to the bis-adenyl compound than to the bis-naphthyl compound because the linker is attached to the sides of the adenyl groups, but to the ends of the naphthyl groups. Finally, ab initio quantum mechanics calculations and energy decompositions for relevant conformations of adenine and naphthalene dimers support the view that stacking in these compounds is driven primarily by nonelectrostatic interactions. The present analysis illustrates the importance of considering all aspects of a molecular system when interpreting experimental data, and the value of computer models as an adjunct to chemical intuition.

Conformational effects on photophysical characteristics of C5-C5'-linked dihydrothymine dimers in solution

Photophysical characteristics of N-substituted C5-C5'-linked dihydrothymine dimers (1a,b[meso], meso compounds of [5R,5'S]-bi-5,6-dihydrothymines; 1a,b[rac], racemic compounds of [5R,5'R]-bi-5,6-dihydrothymines and [5S,5'S]-bi-5,6-dihydrothymines) in aqueous solution with varying contents of less-polar aprotic solvent such as tetrahydrofuran or dioxane have been investigated by UV-absorption, and steady-state and time-resolved fluorescence spectroscopies. Among the C5-C5'-linked dimers, (5R,5'S)-bi-5,6-dihydro-1-methylthymine (1a[meso]) showed a red-shifted weak UV-absorption band at 270-350 nm and excimer fluorescence emission at lambda max = 370 nm with a quantum yield (phi F) of approximately 0.1 in phosphate buffer (pH < 10) at 293 K. Racemic compound of 5,6-dihydro-1-methylthymine dimer (1a[rac]), meso and racemic compounds of 5,6-dihydro-1,3-dimethylthymine dimers (1b[meso] and 1b[rac]) in phosphate buffer were nonfluorescent under similar conditions. The UV-absorption and fluorescence spectral characteristics of 1a[meso] in aqueous solution were interpreted in terms of intramolecular stacking interactions between the dihydropyrimidine chromophores leading to a preferential "closed-shell" conformation both in the ground state and the excited singlet state. In basic solutions at pH > pKa = 11.7, the fluorescence quantum yield of 1a[meso] decreased due to a dominant "open-shell" conformation resulting from the electrostatic repulsion between the deprotonated dihydrothymine chromophores of 1a[meso] in a dianion form.

Thermodynamic parameters for DNA sequences with dangling ends

The thermodynamic contributions to duplex formation of all 32 possible single-nucleotide dangling ends on a Watson-Crick pair are reported. In most instances, dangling ends are stabilizing with free energy contributions ranging from +0.48 (GT(A)) to-0.96 kcal/mol (). In comparison, Watson-Crick nearest-neighbor increments range from -0. 58 (TA/AT) to -2.24 (GC/CG) kcal/mol. Hence, in some cases, a dangling end contributes as much to duplex stability as a Watson-Crick A-T base pair. The implications of these results for DNA probe design are discussed. Analysis of the sequence dependence of dangling-end stabilities show that the nature of the closing base pair largely determines the stabilization. For a given closing base pair, however, adenine dangling ends are always more or equally as stable as the other dangling nucleotides. Moreover, 5' dangling ends are more or equally as stabilizing as their 3' counterparts. Comparison of DNA with RNA dangling-end motifs shows that DNA motifs with 5' dangling ends contribute to stability equally or more than their RNA counterparts. Conversely, RNA 3' dangling ends contribute to stability equally or more than their DNA counterparts. This data set has been incorporated into a DNA secondary structure prediction algorithm (DNA MFOLD) (http://mfold2.wustl.edu/mfold/dna/for m1.cgi) as well as a DNA hybridization prediction algorithm (HYTHERtrade mark) (http://jsl1.chem.wayne.edu/Hyther/hythermenu .html).