Search

Information on how insulin and insulin-like growth factors 1 and 2 (IGF-1 and -2) activate insulin receptors (IR-A and -B) and the IGF-1 receptor (IGF-1R) is crucial for understanding the difference in the biological activities of these peptide hormones. Cryo-EM studies have revealed that insulin uses its binding sites 1 and 2 to interact with IR-A and have identified several critical residues in binding site 2. However, mutagenesis studies suggest that Ile-A10, Ser-A12, Leu-A13, and Glu-A17 also belong to insulin's site 2. Here, to resolve this discrepancy, we mutated these insulin residues and the equivalent residues in IGFs. Our findings revealed that equivalent mutations in the hormones can result in differential biological effects and that these effects can be receptor-specific. We noted that the insulin positions A10 and A17 are important for its binding to IR-A and IR-B and IGF-1R and that A13 is important only for IR-A and IR-B binding. The IGF-1/IGF-2 positions 51/50 and 54/53 did not appear to play critical roles in receptor binding, but mutations at IGF-1 position 58 and IGF-2 position 57 affected the binding. We propose that IGF-1 Glu-58 interacts with IGF-1R Arg-704 and belongs to IGF-1 site 1, a finding supported by the NMR structure of the less active Asp-58-IGF-1 variant. Computational analyses indicated that the aforementioned mutations can affect internal insulin dynamics and inhibit adoption of a receptor-bound conformation, important for binding to receptor site 1. We provide a molecular model and alternative hypotheses for how the mutated insulin residues affect activity.

The influence of various sugar residue modifications on intrinsic energetic, conformational, and mechanical properties of 2'-deoxyribonucleotide-5'-monophosphates (dNs) was comprehensively investigated using modern quantum chemical approaches. In total, fourteen sugar modifications, including double bonds and heteroatoms (S and N) inside the sugar ring, as well as fluorination in various positions, were analyzed. Among hundreds of possible conformational states of dNs, only two - AI and BI, corresponding to the most biologically significant forms of a double-helical DNA, were considered for each dN. It was established that the most of the studied modifications tend to strongly stabilize either AI or BI conformation of dNs both in the gas phase and in aqueous solution (modelled by implicit solvent models). Therefore, some of these modifications can be used as a tool for reducing structural polymorphism of nucleic acids in solution as well as for designing oligonucleotides with specific structural features. The evaluation of relaxed force constants (RFC) for glycosidic bonds suggests that the majority of the studied modifications of the sugar residue yield increased strengths of glycosidic bonds in dNs, and can therefore be used for designing modified nucleic acids with an increased resistance to abasic lesions. The most significant reinforcement of the glycosidic bond occurs in dNs containing the CF2 group instead of the O4' oxygen and the fluorine atom at the 2'-α-position. The calculation of the RFC and vibrational root-mean-square (VRMS) deviations for conformational degrees of freedom revealed a strong dependence between mechanical properties of dNs and their energetic characteristics. In particular, electronic energies of AI and BI conformers of dNs calculated in vacuo are closely connected with the values of relaxed force constants (RFC) for the δ angle: the higher RFC(δ) values correspond to more energetically favorable conformers.

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.

Anion-π interactions have been shown to stabilize flavoproteins and to regulate the redox potential of the flavin cofactor. They are commonly attributed to electrostatic forces. Herein we show that anion-flavin interactions can have a substantial charge-transfer component. Our conclusion emanates from a multi-approach theoretical analysis and is backed by previously reported observations of absorption bands, originating from charge transfer between oxidized flavin and proximate cysteine thiolate groups. This partial covalency of anion-flavin contacts renders classical simulations of flavoproteins questionable.

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 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.

Based on comprehensive conformational analysis at the DFT B3LYP/6-31G(d,p) level of theory [Ye.P. Yurenko et al., J. Phys. Chem. B, 111 (2007) 9655] 13 types of intramolecular hydrogen bonds (171 in total number, from 1 to 3 H-bonds per conformer) were identified in the 92 allowed conformers of isolated thymidine, a DNA canonical nucleoside. Weak interactions CH ⋯ O , OH ⋯ HC and CH ⋯ HC were analyzed by all the Koch and Popelier criteria of H-bonding. Quite satisfactory linear correlation was found between Grabowski complex measure of H-bond strength Δ com and Iogansen H-bond energy −Δ H for 17 OH⋯O hydrogen bonds.

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 degreesor = Por = 40.1 degrees), and the rest possess south (S) (with 157.1 degreesor = Por = 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.

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.

Structural and energetic features of arti fi cial DNA quadruplexes consisting of base tetrads and their stacks with Na + /K + ion(s) inside the central pore and incorporating halogenated derivatives of xanthine, 8- fl uoro-9-deazaxanthine (FdaX), 8- chloro-9-deazaxanthine (CldaX), 8-bromo-9-deazaxanthine (BrdaX), or 8-iodo-9-deazaxanthine (IdaX), have been investigated by modern state-of-the-art computational tools. The DNA (or RNA) quadruplex models based on 8-halo-9-deazaxanthines are predicted to be more stable relative to those with unmodi fi ed xanthine due to the increased stabilizing contributions coming from all three main types of weak interactions (H-bonding, stacking, and ion coordination). Methods for analyzing the electron density are used to understand the nature of forces determining the stability of the system and to gain a predictive potential. Quadruplex systems incorporating polarizable halogen atoms (chlorine, bromine, or iodine) bene fi t signi fi cantly from the stabilizing stacking between the individual tetrads due to an increased dispersion contribution as compared to xanthine and guanine, natural references used. Ion coordination induces a signi fi cant rearrangement of electron density in the quadruplex stem as visualized by electron deformation density (EDD) and analyzed by ETS-NOCV and Voronoi charges. Na + induces larger electron polarization from the quadruplex toward the ion, whereas K + has a higher propensity to electron sharing (identi fi ed by QTAIM delocalization index). We expect that our results will contribute to the development of novel strategies to further modify and analyze the natural G-quadruplex core.

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.

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-U2, U-U3 and A-C nucleobase pairs the CH···O/N H-bonds play a prominent role (30%) in their stabilization. We suppose that unconventional CH···O/N H-bond plays the role of the third "fulcrum", ensuring structurally dynamic similarity of the isomorphic base pairs of different origin, which are incorporated equally well into the structure of the DNA double helix.

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.

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.

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 (OH⋯O) and non-conventional (CH⋯O, CH⋯HC) H-bonds as well as closed-shell van der Waals (C⋯O) 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.

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.

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'-B or 2'-A,B positions somewhat stabilizes in vacuo the AI relative to the BI conformations. In contrast, substitution of the CF2 group for the O4' atom (O4'-CF2 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'-CF2 derivative. Consequently, the O4'-CF2 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.