The role of solvent on the base stacking properties of the stacked cytosine dimer

Nature and magnitude of aromatic stacking of nucleic acid bases

This review summarises recent advances in quantum chemical calculations of base-stacking forces in nucleic acids. We explain in detail the very complex relationship between the gas-phase base-stacking energies, as revealed by quantum chemical (QM) calculations, and the highly variable roles of these interactions in nucleic acids. This issue is rarely discussed in quantum chemical and physical chemistry literature. We further extensively discuss methods that are available for base-stacking studies, complexity of comparison of stacking calculations with gas phase experiments, balance of forces in stacked complexes of nucleic acid bases, and the relation between QM and force field descriptions. We also review all recent calculations on base-stacking systems, including details analysis of the B-DNA stacking. Specific attention is paid to the highest accuracy QM calculations, to the decomposition of the interactions, and development of dispersion-balanced DFT methods. Future prospects of computational studies of base stacking are discussed.

Stability of nucleic acid base pairs in organic solvents: molecular dynamics, molecular dynamics/quenching, and correlated ab initio study

The dynamic structure and potential energy surface of adenine...thymine and guanine...cytosine base pairs and their methylated analogues interacting with a small number (from 1 to 16 molecules) of organic solvents (methanol, dimethylsulfoxide, and chloroform) were investigated by various theoretical approaches starting from simple empirical methods employing the Cornell et al. force field to highly accurate ab initio quantum chemical calculations (MP2 and particularly CCSD(T) methods). After the simple molecular dynamics simulation, the molecular dynamics in combination with quenching technique was also used. The molecular dynamics simulations presented here have confirmed previous experimental and theoretical results from the bulk solvents showing that, whereas in chloroform the base pairs create hydrogen-bonded structures, in methanol, stacked structures are preferred. While methanol (like water) can stabilize the stacked structures of the base pairs by a higher number of hydrogen bonds than is possible in hydrogen-bonded pairs, the chloroform molecule lacks such a property, and the hydrogen-bonded structures are preferred in this solvent. The large volume of the dimethylsulfoxide molecule is an obstacle for the creation of very stable hydrogen-bonded and stacked systems, and a preference for T-shaped structures, especially for complexes of methylated adenine...thymine base pairs, was observed. These results provide clear evidence that the preference of either the stacked or the hydrogen-bonded structures of the base pairs in the solvent is not determined only by bulk properties or the solvent polarity but rather by specific interactions of the base pair with a small number of the solvent molecules. These conclusions obtained at the empirical level were verified also by high-level ab initio correlated calculations.

Role of nitrogen atom in aromatic stacking

Aromatic stacking of 9,9'-(alpha,omega-alkanediyl)bis[adenine] (1), 1,1'-(alpha,omega-alkanediyl)bis[benzimidazole] (2), and 9-[omega-(benzimidazol-1-yl)alkyl]adenine (3) were studied at low concentrations of these compounds by means of UV and NMR spectroscopies. The UV hypochromic effect at T degrees C was determined as the ratio of the integration strength B at T degrees C (T = 27, 40, and 50) to that at 60 degrees C. The UV hypochromic effects of 1 and 3 were remarkable in water, suggesting a formation of intramolecular aromatic stacking, while the UV data of 2 did not present unambiguous evidence supporting aromatic stacking. A difference of chemical shift of each aromatic ring proton between 27 and 80 degrees C, that is Deltadelta = delta(80 degrees C) - delta(27 degrees C), was given as an indication of the aromatic stacking in the NMR study. On the basis the data of Deltadelta, 1 and 3 were stabilized by a stacking interaction in the buffer solution at pD 7.0 but not in the organic solvents. On the other hand, the NMR data did not indicate the formation of aromatic stacking of 2 either in the organic solvents or in the aqueous solution. The thermodynamic parameters of the intramolecular aromatic stacking of 3 were determined by means of NMR spectroscopy.

Hydration and stability of nucleic acid bases and base pairs

Empirical, quantum chemical calculations and molecular dynamics simulations of the role of a solvent on tautomerism of nucleic acid bases and structure and properties of nucleic acid base pairs are summarized. Attention was paid to microhydrated (by one and two water molecules) complexes, for which structures found by scanning of empirical potential surfaces were recalculated at a correlated ab initio level. Additionally, isolated as well as mono- and dihydrated H-bonded, T-shaped and stacked structures of all possible nucleic acid base pairs were studied at the same theoretical levels. We demonstrate the strong influence of a solvent on the tautomeric equilibrium between the tautomers of bases and on the spatial arrangement of the bases in a base pair. The results provide clear evidence that the prevalence of either the stacked or hydrogen-bonded structures of the base pairs in the solvent is not determined only by its bulk properties, but rather by specific hydrophilic interactions of the base pair with a small number of solvent molecules.