Intramolecular CH···O hydrogen bonds in the AI and BI DNA-like conformers of canonical nucleosides and their Watson-Crick pairs. Quantum chemical and AIM analysis.

The aim of this work is to cast some light on the H-bonds in double-stranded DNA in its AI and BI forms. For this purpose, we have performed the MP2 and DFT quantum chemical calculations of the canonical nucleoside conformers, relative to the AI and BI DNA forms, and their Watson-Crick pairs, which were regarded as the simplest models of the double-stranded DNA. Based on the atoms-in-molecules analysis (AIM), five types of the CH···O hydrogen bonds, involving bases and sugar, were detected numerically from 1 to 3 per a conformer: C2'H···O5', C1'H···O2, C6H···O5', C8H···O5', and C6H···O4'. The energy values of H-bonds occupy the range of 2.3-5.6 kcal/mol, surely exceeding the kT value (0.62 kcal/mol). The nucleoside CH···O hydrogen bonds appeared to "survive" turns of bases against the sugar, sometimes in rather large ranges of the angle values, pertinent to certain conformations, which points out to the source of the DNA lability, necessary for the conformational adaptation in processes of its functioning. The calculation of the interactions in the dA·T nucleoside pair gives evidence, that additionally to the N6H···O4 and N1···N3H canonical H-bonds, between the bases adenine and thymine the third one (C2H···O2) is formed, which, though being rather weak (about 1 kcal/mol), satisfies the AIM criteria of H-bonding and may be classified as a true H-bond. The total energy of all the CH···O nontraditional intramolecular H-bonds in DNA nucleoside pairs appeared to be commensurable with the energy of H-bonds between the bases in Watson-Crick pairs, which implies their possible important role in the DNA shaping.

Tautomeric equilibrium of uracil and thymine in model protein-nucleic acid contacts. Spectroscopic and quantum chemical approach.

This work deals with tautomeric transformations of uracil (Ura) and thymine (Thy) in their model complexes with the deprotonated carboxylic group. Essential changes in the UV spectra of the bases upon their interaction with NaAc, vanishing signals of both imino protons in (1)H NMR spectra, and a perceptible decrease in intensity of both IR bands, related to the stretching vibrations nu(C=O) of the carbonyl groups, imply involvement of enolic tautomers. Results of quantum chemical calculations of the double complexes of the Ura(Thy) tautomers with CH(3)COO(-) at the MP2/6-311++G(2df,pd)//B3LYP/6-311++G(d,p) level of theory proved to be incompatible with the spectral features: despite the fact that the complexes of the enolic tautomers are much closer in energy to the diketo ones as compared to isolated tautomers, the energy gap between them is such that in tautomeric equilibrium dominate diketo forms. Calculations of triple complexes of the type CH(3)COO(-):Ura(Thy) tautomer:Na(+), taking into account the effect of the Na(+) coordination with tautomers, show that three triple complexes formed by enolic tautomers appeared more stable than those formed by diketo ones. This makes the UV and (1)H NMR data understandable, but the high residual intensity of the nu(C=O) bands in the IR spectra remains unclear. At that ion, Na(+) itself was not able to disturb the tautomeric equilibrium in the coordination complexes of the type Ura(Thy) tautomer:Na(+). To evaluate the DMSO effect, the CPCM solvation model was applied to triple complexes of the Ura tautomers. It appeared that in the solution there is coexistence between the diketo and enolic tautomers in a ratio of 53%:47%. This makes possible reconciliation of our experimental data. The biological significance of high-energy tautomers of nucleotide bases is discussed.

Ab initio comprehensive conformational analysis of 2'-deoxyuridine, the biologically significant DNA minor nucleoside, and reconstruction of its low-temperature matrix infrared spectrum.

A comprehensive conformational analysis of isolated 2'-deoxyuridine (dU), a minor DNA nucleoside, has been performed by means of ab initio calculations at the MP2/6-311++G (d,p)//DFT B3LYP/6-31G (d,p) level of theory. At 298.15 and 420 K, all 94 allowed conformers of dU are within 8.96 and 7.91 kcal/mol Gibbs energy ranges, respectively. Syn orientation for the base and South (S) conformers for the sugar dominate at 298.15 K: syn/anti=62.3%:37.7% and S/N=77.2%:22.8%. At 420 K in the majority of conformers, the base is anti oriented and the population of North (N) sugars increases: syn/anti=39.3%:60.7% and S/N=63.0%:37.0%. Values of all conformational parameters and correlations between them, as well as their correlations with valence bonds, and also correlations between valence bonds and angles were estimated. In general, 14 types of intramolecular H-bonds were detected (1-3 H-bonds per conformer, the total number 175), namely, C1'H...O2 (16 H-bonds), C2'H1...O5' (9), C2'H2...O2 (21), C3'H...O2 (21), C5'H1...O2 (14), C5'H2...O2 (11), C6H...O4' (37), C6H...O5' (22), C3'H...HC6 (4), O5'H...HC6 (2), O3'H...O5' (5), O5'H...O4' (1), O5'H...O3' (4), and O5'H...O2 (8). Geometric, vibrational, structural-topological, and energetic features of the OH...O intramolecular H-bonds in dU conformers were determined. The close similarity between energetic and geometric characteristics of dU and thymidine DNA-like conformers in anti and relevant syn conformations and their transition states of the anti-->syn interconversion implies that mismatch DNA glycosylase discriminates between the two nucleosides, mainly because of the difference in the shapes of their bases. Convolution of calculated IR spectra of all the dU conformers within the limits 3400-3700 cm(-1) appears to be consistent with its low-temperature matrix IR spectrum (Ivanov et al. Spectrochim. Acta, Part A 2003, 59, 1959), wavenumber discrepancy not exceeding 1%. It was concluded that, for a reliable reproduction of the experimental spectrum, the whole set of conformers should be taken into consideration. The suggested method makes reconstruction of the isolated nucleoside IR spectrum at a physiological interval of temperature reasonably possible.

How many conformers determine the thymidine low-temperature matrix infrared spectrum? DFT and MP2 quantum chemical study.

A comprehensive conformational analysis of isolated 2'-beta-deoxy-thymidine (T), canonical DNA nucleoside, has been performed by means of ab initio calculations at the MP2/6-311++G(d,p)//DFT B3LYP/6-31G(d,p) level of theory. At 298.15 K, all 92 conformers of isolated dT are within a 7.49 kcal/mol Gibbs energy range. Syn orientation for the base and South (S) conformers for the sugar dominate at this temperature: syn/anti = 61.6%:38.4% and S/N = 74.5%:25.5%. However, at 420 K, the majority of conformers contain anti base and the population of North (N) sugars increases: syn/anti = 38.0%:62.0% and S/N = 59.5%:40.5%. The whole conformational parameters (P, chi, gamma, delta, beta, epsilon, nu max) were analyzed as well as the energies of the OH...O intramolecular H-bonds on the basis of nu(OH) stretching vibrations. Convolution of calculated IR spectra of all of the T conformers appears consistent with its low-temperature matrix spectrum (Ivanov et al. Low Temp. Phys. 2003, 29, 809). The maximal discrepancy in frequencies between calculated and experimental spectra is less than 1%. A conclusion was made that for reliable reconstruction of the isolated nucleoside IR spectrum the quasi whole set of conformers should be taken into consideration. In essence, this result opens up a possibility to reconstruct IR spectra of isolated nucleosides at physiological temperatures with rather satisfactory probability.

Comprehensive conformational analysis of the nucleoside analogue 2'-beta-deoxy-6-azacytidine by DFT and MP2 calculations.

A comprehensive conformational analysis of isolated 2'-beta-deoxy-6-azacytidine (d6AC), an analogue of therapeutically active 6-azacytidine (6AC), has been performed by means of ab initio calculations at the MP2/6-311++G(2df,pd)//DFT B3LYP/6-31G(d,p) level of theory. Among the 81 conformers located within a 7.83 kcal/mol Gibbs energy range at T = 298.15 K, 38 contain syn-oriented bases with respect to 2'-deoxyribose; the other conformers include anti-oriented bases. Energetic analysis of these conformers shows that conformational equilibrium of isolated d6AC at T = 298.15 K is shifted to syn conformation with a syn/anti ratio estimated as 61.4%:38.6%. As far as the sugar conformation is concerned, 40 conformers contain north (N) (with 0.3 degrees < or = P < or = 40.1 degrees), and the rest possess south (S) (with 157.1 degrees < or = P < or = 207.0 degrees) puckers, where P is the pseudorotational angle of the furanose ring. The S/N occupancy ratio is estimated as 80.2%:19.8% (T = 298.15 K). The two most stable conformers are energetically quasidegenerate and correspond to both C2'-endo/syn conformers differing only by orientation of the O3'H hydroxyl group. They are both stabilized by means of similar intramolecular H-bonds, i.e., O5'H...O2, C2'H2...O2, and C2'H2...O5'. As examined by AIM criteria, from 1 to 3 H-bonds per conformer were identified among 13 possible interactions: O5'H...O2, O5'H...N6, O3'H...O5', O5'H...O3', C1'H...O2, C2'H2...O2, C2'H2...O5', C3'H...O2, C3'H...N6, C5'H1...O2, C5'H2...O2, C5'H1...N6, and C5'H2...N6. The biological effect of d6AC is conceived as an inhibition of replicative DNA polymerase caused by an unusual orientation of the sugar residue against the base in the only A form DNA-like conformer.

ATR-IR spectroscopy as applied to nucleic acid films.

For the first time the ATR technique was applied to obtain IR absorption spectra of DNA and RNA dry films. There was worked out procedure of the nucleic acid removal from germanium plate, which obviously was a main obstacle to application of ATR-IR spectroscopy to nucleic acids. This technique of IR spectroscopy was applied to confirmation of RNA tropism of aurin tricarboxylic acid observed by molecular biological methods.

Structural transitions in polycytidylic acid: proton buffer capacity data.

The pH-dependences of proton buffer capacity of poly(C) were computed on the basis of the literature data. In these curves there were observed four peaks: two narrow and two wide ones. The first narrow peak reflects the process of cooperative formation of double helices, which is induced by protonation of the N3 atom of nucleotide bases. The first wide peak is assigned to noncooperative process of poly(C) double helices protonation at the N3 nitrogen atom. It is proposed that the second wide peak corresponds to noncooperative protonation of the neutral cytosine bases at the oxygen atom. This reaction causes cooperative dissociation of the poly(C) double helices. The second narrow peak reflects the dissociation process.