Substituting CF2 for O4′ in Components of Nucleic Acids: Towards Systems with Reduced Propensity to Form Abasic Lesions

Intrinsic structural features and energetics of nucleotides containing variously fluorinated sugars as potential building blocks of DNA duplexes and quadruplexes are explored systematically using the modern methods of density functional theory (DFT) and quantum chemical topology (QCT). Our results suggest that fluorination at the 2′‐β or 2′‐α,β positions somewhat stabilizes in vacuo the AI relative to the BI conformations. In contrast, substitution of the CF₂ group for the O4′ atom (O4′‐CF₂ modification) leads to a preference of the BI relative to AI DNA‐like conformers. All the studied modifications result in a noticeable increase in the stability of the glycosidic bond [estimated by the relaxed force constants (RFC) approach], with particularly encouraging results for the O4′‐CF₂ derivative. Consequently, the O4′‐CF₂ modified systems are suggested and explored as promising scaffolds for the development of duplex and quadruplex structures with reduced propensity to form abasic lesions and to undergo DNA damage.     

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.     

Designing a New Class of Bases for Nucleic Acid Quadruplexes and Quadruplex‐Active Ligands

A new class of quadruplex nucleobases, derived from 3‐deazaguanine, has been designed for various applications as smart quadruplex ligands as well as quadruplex‐based aptamers, receptors, and sensors. An efficient strategy for modifying the guanine quadruplex core has been developed and tested by using quantum chemistry methods. Several potential guanine derivatives modified at the 3‐ or 8‐position or both are analyzed, and the results compared to reference systems containing natural guanine. Analysis of the formation energies (BLYP‐D3(BJ)/def2‐TZVPP level of theory, in combination with the COSMO model for water) in model systems consisting of two and three stacked tetrads with Na⁺/K⁺ ion(s) inside the internal channel indicates that the formation of structures with 3‐halo‐3‐deazaguanine bases leads to a substantial gain in energy, as compared to the corresponding reference guanine complexes. The results cast light on changes in the noncovalent interactions (hydrogen bonding, stacking, and ion coordination) in a quadruplex stem upon modification of the guanine core. In particular, the enhanced stability of the modified quadruplexes was shown to originate mainly from increased π–π stacking. Our study suggests the 3‐halo‐3‐deazaguanine skeleton as a potential building unit for quadruplex systems and smart G‐quadruplex ligands.     

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.