This study aims to cast light on the physico-chemical nature and energetic of the non-conventional CH···O/N H-bonds in the biologically important natural nucleobase pairs using a comprehensive quantum-chemical approach. As a whole, the 36 biologically important pairs, involving canonical and rare tautomers of nucleobases, were studied by means of all available up-to-date state-of-the-art quantum-chemical techniques along with quantum theory "Atoms in molecules" (QTAIM), Natural Bond Orbital (NBO) analysis, Grunenberg's compliance constants theory, geometrical and vibrational analyses to identify the CH···O/N interactions, reveal their physico-chemical nature and estimate their strengths as well as contribution to the overall base-pairs stability. It was shown that all the 38 CH···O/N contacts (25 CH···O and 13 CH···N H-bonds) completely satisfy all classical geometrical, electron-topological, in particular Bader's and "two-molecule" Koch and Popelier's, and vibrational criteria of H-bonding. The positive values of Grunenberg's compliance constants prove that the CH···O/N contacts in nucleobase pairs are stabilizing interactions unlike electrostatic repulsion and anti-H-bonds. NBO analysis indicates the electron density transfer from the lone electron pair of the acceptor atom (O/N) to the antibonding orbital corresponding to the donor group σ(∗)(CH). Moreover, significant increase in the frequency of the out-of-plane deformation modes γ (CH) under the formation of the CH···O (by 17.2÷81.3/10.8÷84.7 cm(-1)) and CH···N (by 32.7÷85.9/9.0÷77.9 cm(-1)) H-bonds at the density functional theory (DFT)/second-order Møller-Plesset (MP2) levels of theory, respectively, and concomitant changes of their intensities can be considered as reliable indicators of H-bonding. The strengths of the CH···O/N interactions, evaluated by means of Espinosa-Molins-Lecomte formula, lie within the range 0.45÷3.89/0.62÷4.10 kcal/mol for the CH···O H-bonds and 1.45÷3.17/1.70÷3.43 kcal/mol for the CH···N H-bonds at the DFT/MP2 levels of theory, respectively. We revealed high linear mutual correlations between the H-bond energy and different physico-chemical parameters of the CH···O/N H-bonds. Based on these observations, the authors asserted that the most reliable descriptors of the H-bonding are the electron density ρ at the СН···О/N H-bond critical points and the NBO calculated stabilization energy E((2)). The linear dependence of the H-bond energy ECH···O/N (in kcal/mol) on the electron density ρ (in atomic units) was established (DFT/MP2): ECH···O = 248.501[Formula: see text]ρ-0.367/260.518[Formula: see text]ρ-0.373 and ECH···N = 218.125[Formula: see text]ρ-0.339/243.599[Formula: see text]ρ-0.441. Red-shifted and blue-shifted CH···O/N H-bonds behave in a similar way and can be described with the same fit parameters. It was found that the A-U HH2 and U-U3 nucleobase pairs are stabilized solely by the CH···O/N H-bonds. At the same time, in the A-U HH1, A-U HH2, A-Asyn 1, A-Asyn 2, A-Asyn 3, A-A4, A-G1, A-G2, G-U1, G-U2, G-U3, G-C HH1, U-U1, U-...

The study aimed to cast light on the structure and internal energetics of guanine- and xanthine-based model DNA quadruplexes and the physico-chemical nature of the non-covalent interactions involved. Several independent approaches were used for this purpose: DFT-D3 calculations, Quantum Theory of Atoms in Molecules, Natural Bond Orbital Analysis, Energy Decomposition Analysis, Compliance Constant Theory, and Non-Covalent Interaction Analysis. The results point to an excellent degree of structural and energetic compatibility between the two types of model quadruplexes. This fact stems from both the structural features (close values of van der Waals volumes, pore radii, geometrical parameters of the H-bonds) and the energetic characteristics (comparable values of the energies of formation). It was established that hydrogen bonding makes the greatest (∼50%) contribution to the internal stability of the DNA quadruplexes, whereas the aromatic base stacking and ion coordination terms are commensurable and account for the rest. Energy decomposition analysis performed for guanine (Gua) and xanthine (Xan) quartets B4 and higher-order structures consisting of two or three stacked quartets indicates that whereas Gua structures benefit from a high degree of H-bond cooperativity, Xan models are characterized by a more favorable and cooperative π-π stacking. The results of electron density topological analysis show that Na(+)/K(+) ion coordination deeply affects the network of non-covalent interactions in Gua models due to the change in the twist angle between the stacked tetrads. For Xan models, ion coordination makes tetrads in stacks more planar without changing the twist angle. Therefore, the presence of the ion seems to be essential for the formation of planar stacks in Xan-based DNA quadruplexes. Detailed study of the nature of ion-base coordination suggests that this interaction has a partially covalent character and cannot be considered as purely electrostatic. Investigation of the H-bond and ion-base coordination strengths by various independent approaches agrees well with the results of QTAIM analysis.

Guanine DNA quadruplexes are interesting and important biological objects because they represent potential targets for regulatory drugs. Their use as building blocks for biomaterial applications is also being investigated. This contribution reports the in silico design of artificial building blocks derived from xanthine. Methods of quantum chemistry were used to evaluate the properties of xanthine structures relative to those containing guanine, the natural reference used. Tailoring the xanthine core showed that the base stacking and the ion coordination were significantly enhanced in the designed systems, while the ion-transport properties were not affected. Our study suggests that the 9-deaza-8-haloxanthine bases (where the halogen is fluorine or chlorine) are highly promising candidates for the development of artificial quadruplexes and quadruplex-active ligands.

Ab initio quantum-chemical study of specific point contacts of replisome proteins with DNA modeled by acetic acid with canonical and mutagenic tautomers of DNA bases methylated at the glycosidic nitrogen atoms was performed in vacuo and continuum with a low dielectric constant (ϵ ∼ 4) corresponding to a hydrophobic interface of protein-nucleic acid interaction. All tautomerized complexes were found to be dynamically unstable, because the electronic energies of their back-reaction barriers do not exceed zero-point vibrational energies associated with the vibrational modes whose harmonic vibrational frequencies become imaginary in the transition states of the tautomerization reaction. Additionally, based on the physicochemical arguments, it was demonstrated that the effects of biomolecular environment cannot ensure dynamic stabilization. This result allows suggesting that hypothetically generated by DNA-binding proteins of replisome rare tautomers will have no impact on the total spontaneous mutation due to the low reverse barrier allowing a quick return to the canonical form.

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.

A comprehensive quantum-chemical conformational analysis of two nucleoside analogues, 2',3'-didehydro-2',3'-dideoxyuridine (d4U) and 2',3'-didehydro-2',3'-dideoxycytidine (d4C), is reported. The electronic structure calculations were performed at the MP2/6-311++G(d,p)//B3LYP/6-31++G(d,p) level of theory. It was found that d4U and d4C adopt 20 conformers and 19 conformers, respectively, which correspond to local minima on the respective potential energy landscapes. QTAIM and NBO analyses show that the d4U and d4C conformers are stabilised by a complicated network of specific intramolecular interactions, which includes conventional (OHO) and non-conventional (CHO, CHHC) H-bonds as well as closed-shell van der Waals (CO) contacts. A satisfactory linear correlation was found between Grunenberg's compliance constants for closed-shell intramolecular interactions and their energy. It is shown that there are no conformational obstacles for incorporation of d4U and d4C into the double helical A and B forms of DNA. The less pronounced biological activity of d4U as compared to 2',3'-didehydro-2',3'-dideoxythymidine (d4T) is most likely due to the presence of the bulky methyl group at the 5-position of d4T, which can be recognised by target enzymes.

This paper is a logical continuation of the theoretical survey of the CH⋯O/N specific contacts in the nucleobase pairs using a wide arsenal of the modern methods, which was initiated in our previous study [J. Biomol. Struct. & Dynam., 2014, 32, 993-1022]. It was established that 34 CH⋯O and 7 CH⋯N interactions, that were detected by quantum-chemical calculations in the 39 biologically important pairs involving modified nucleobases, completely satisfy all geometrical, vibrational, electron-topological, in particular Bader's and "two-molecule" Koch and Popelier's, Grunenberg's compliance constants theory and natural bond orbital criteria indicating that they can be identified as true H-bonds. The geometrical criteria of the H-bond formation are fulfilled for all considered CH⋯O/N H-bonds without any exception. It was shown that the classical rule of the stretching vibration shifts does not work in the ~95% cases of the CH⋯O/N H-bonds. Furthermore, significant increase in the frequency of the out-of-plane deformation modes γ(CH) under the formation of CH⋯O/N H-bonds and corresponding changes of their intensities can be also considered as reliable indicators of the H-bonding. We revealed high linear mutual correlations between the electron density, Laplacian of the electron density, H-bond energy at the (3, -1) bond critical points of the CH⋯O/N H-bonds, and different physico-chemical parameters of the CH⋯O/N H-bonds. We suggested that the electron density ρ and the interaction energy E((2)) of the lone orbital pairs are the most reliable descriptors of the H-bonding. The linear dependence of the H-bond energy ECH⋯O/N on the electron density ρ was established: ECH⋯O = 250.263∙ρ - .380/258.255∙ρ - .396 and ECH⋯N = 196.800∙ρ - .172/268.559∙ρ - .703 obtained at the density functional theory (DFT)/Møller-Plesset (MP2) levels of theory, respectively. The studies of the interaction energies show that the contribution of the CH⋯O and CH⋯N H-bonds into the base pairs stability varies from 3.0/4.2 to 35.1/31.2% and from 3.0/4.3 to 44.4/46.5% at the DFT/MP2 levels of theory, accordingly. Energy decomposition analysis performed for all base pairs involving canonical and modified nucleobases defines the electrostatic attraction and Pauli repulsion as dominant stabilizing forces in all complexes. This observation was additionally confirmed by the results of the QTAIM delocalization indexes analysis. The studies reported here advance our understanding of the biological role of the weak CH⋯O/N H-bonds, that dictates the requirements for the structural and dynamical similarity of the canonical and mismatched pairs with Watson-Crick (WC) geometry, which facilitates their enzymatic incorporation into the DNA double helix during DNA replication. Thus, these H-bonds in the base pairs with WC geometry may be also considered as "the last drop" at the transmission of the electronic signal that launches the chemical incorporation of the incoming nucleoside triphosphate into DNA.

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.