Study of the Molecular Architectures of 2-(4-Chlorophenyl)-5-(pyrrolidin-1-yl)-2H-1,2,3-triazole-4-carboxylic Acid Using Their Vibrational Spectra, Quantum Chemical Calculations and Molecular Docking with MMP-2 Receptor

1,2,3-triazole skeleton is a valuable building block for the discovery of new promising anticancer agents. In the present work, the molecular structure of the synthesized anticancer drug 2-(4-chlorophenyl)-5-(pyrrolidin-1-yl)-2H-1,2,3-triazole-4-carboxylic acid (1b) and its anionic form (2b) was characterized by means of the B3LYP, M06-2X and MP2 quantum chemical methods, optimizing their monomer, cyclic dimer and stacking forms using the Gaussian16 program package. The molecular structure was found to be slightly out of plane. The good agreement between the IR and Raman bands experimentally observed in the solid state with those calculated theoretically confirms the synthesized structures. All of the bands were accurately assigned according to functional calculations (DFT) in the monomer and dimer forms, together with the polynomic scaling equation procedure (PSE). Therefore, the effect of the substituents on the triazole ring and the effect of the chlorine atom on the molecular structure and on the vibrational spectra were evaluated through comparison with its non-substituted form. Through molecular docking calculations, it was evaluated as to how molecule 1b interacts with few amino acids of the MMP-2 metalloproteinase receptor, using Sybyl-X 2.0 software. Thus, the relevance of triazole scaffolds in established hydrogen bond-type interactions was demonstrated.

Peculiarities of the Spatial and Electronic Structure of 2-Aryl-1,2,3-Triazol-5-Carboxylic Acids and Their Salts on the Basis of Spectral Studies and DFT Calculations

The molecular structure and vibrational spectra of six 1,2,3-triazoles-containing molecules with possible anticancer activity were investigated. For two of them, the optimized geometry was determined in the monomer, cyclic dimer and stacking forms using the B3LYP, M06-2X and MP2 methods implemented in the GAUSSIAN-16 program package. The effect of the para-substitution on the aryl ring was evaluated based on changes in the molecular structure and atomic charge distribution of the triazole ring. An increment in the positive N4 charge was linearly related to a decrease in both the aryl ring and the carboxylic group rotation, with respect to the triazole ring, and by contrast, to an increment in the pyrrolidine ring rotation. Anionic formation had a larger effect on the triazole ring structure than the electronic nature of the different substituents on the aryl ring. Several relationships were obtained that could facilitate the selection of substituents on the triazole ring for their further synthesis. The observed IR and Raman bands in the solid state of two of these compounds were accurately assigned according to monomer and dimer form calculations, together with the polynomic scaling equation procedure (PSE). The large red-shift of the C=O stretching mode indicates that strong H-bonds in the dimer form appear in the solid state through this group.

Base pairs with 4-amino-3-nitrobenzonitrile: comparison with the natural WC pairs. Dimer and tetramer forms, Infrared and Raman spectra, and several proposed antiviral modified nucleosides

Base pairs of 4-amino-3-nitrobenzonitrile (4A-3NBN) molecule with uracil, thymine and cytosine nucleobases were optimized and compared to natural Watson-Crick (WC) pairs. The slightly greater flexibility of the -NO₂ group of 4A-3NBN than the N3-H group of the natural nucleobases together with a noticeable higher dipole moment of its pairs can facilitate disruption of the DNA/RNA helix formation. Several new mutagenic modified nucleosides with 4A-3NBN and 3-amino-2-nitrobenzonitrile (3A-2NBN) were proposed as antiviral prodrugs and their base pairs optimized. The special characteristics of these prodrugs appear appropriated for their clinical use. The counterpoise (CP) corrected interaction energies of the base pairs were calculated and compared to the natural ones. The M06-2X DFT method was used for this purpose. The molecular structure of 4A-3NBN was analyzed in detail and the crystal unit cell was simulated by a tetramer form and eight dimer forms. The performance of the B3LYP, X3LYP and M06-2X methods was tested on the vibrational wavenumbers in the monomer, dimer and tetramer forms of 4A-3NBN. The observed IR and Raman bands were assigned according to the optimum dimer II form determined by B3LYP and by the tetramer form calculated by M06-2X, which is the expected unit cell that forms the crystal net. The two best scaling procedures were used.Communicated by Ramaswamy H. Sarma.

How Strong is a Reverse Dative Bond? Compliance Constants as Unique Bond Strength Descriptors

The concepts of bond strength and bond order are at the very heart of chemistry. Nevertheless, the study of local mechanical bond strengths in large polyatomic molecules is a nontrivial task. In a recent publication, Peters and Lehnert et al. [Inorg. Chem. 2020, 59 (20), 14967-14982] found a quite strong (reverse) dative Fe···B bond for a synthesized low-spin {FeNO}¹⁰ complex supported by an ambiphilic trisphosphineborane ligand. Their conclusion was, among other criteria, based on a surprisingly high Fe-B force constant of 1.56 mdyn/Å. Nevertheless, this high value might be based on numerical ambiguities of (rigid) real-space force constants, a flaw, that may affect many other combined spectroscopic and electronic structure studies.

Interplay of Ring Puckering and Hydrogen Bonding in Deoxyribonucleosides

The Cremer-Pople ring puckering analysis and the Konkoli-Cremer local mode analysis supported by the topological analysis of the electron density were applied for the first comprehensive analysis of the interplay between deoxyribose ring puckering and intramolecular H-bonding in 2'-deoxycytidine, 2'-deoxyadenosine, 2'-deoxythymidine, and 2'-deoxyguanosine. We mapped for each deoxyribonucleoside the complete conformational energy surface and the corresponding pseudorotation path. We found only incomplete pseudorotation cycles, caused by ring inversion, which we coined as pseudolibration paths. On each pseudolibration path a global and a local minimum separated by a transition state were identified. The investigation of H-bond free deoxyribonucleoside analogs revealed that removal of the H-bond does not restore the full conformational flexibility of the sugar ring. Our work showed that ring puckering predominantly determines the conformational energy; the larger the puckering amplitude, the lower the conformational energy. In contrast no direct correlation between conformational energy and H-bond strength was found. The longest and weakest H-bonds are located in the local minimum region, whereas the shortest and strongest H-bonds are located outside the global and local minimum regions at the turning points of the pseudolibration paths, i.e., H-bonding determines the shape and length of the pseudolibration paths. In addition to the H-bond strength, we evaluated the covalent/electrostatic character of the H-bonds applying the Cremer-Kraka criterion of covalent bonding. H-bonding in the puric bases has a more covalent character whereas in the pyrimidic bases the H-bond character is more electrostatic. We investigated how the mutual orientation of the CH₂OH group and the base influences H-bond formation via two geometrical parameters describing the rotation of the substituents perpendicular to the sugar ring and their tilting relative to the ring center. According to our results, rotation is more important for H-bond formation. In addition we assessed the influence of the H-bond acceptor, the lone pair (N, respectively O), via the delocalization energy. We found larger delocalization energies corresponding to stronger H-bonds for the puric bases. The global minimum conformation of 2'-deoxyguanosine has the strongest H-bond of all conformers investigated in this work with a bond strength of 0.436 which is even stronger than the H-bond in the water dimer (0.360). The application of our new analysis to DNA deoxyribonucleotides and to unnatural base pairs, which have recently drawn a lot of attention, is in progress.

The maximum occupancy condition for the localized property-optimized orbitals

It is shown analytically that the Chemist's Localized Property-optimized Orbitals (CLPOs), which are the localized orbitals obtainable from the results of ab initio calculations by the open-source program JANPA (http://janpa.sourceforge.net/) according to the recently proposed optimal property partitioning condition, form the Lewis structure with nearly maximum possible total electron occupancy. The conditions required for this additional optimality to hold are discussed. In particular, when a single-determinant wavefunction is used to describe the molecular system without a noticeable electron delocalization, CLPOs derived from this wavefunction approximately optimize the same target quantity as the Natural Bond Orbitals (NBOs), establishing in this way the link between the two sets of localized orbitals. The performance of CLPO and NBO methods is compared by using a dataset containing 7101 small molecules, and the relevant methodological features of both methods are discussed.

Atomic charges for conformationally rich molecules obtained through a modified principal component regression

A modification of the principal component regression model is proposed for obtaining a fixed set of atomic charges (referred to as dipole-derived charges) optimized for reproducing the dipole moment of a conformationally rich molecule, i.e., a molecule with multiple local minima on the potential energy surface. The method does not require any adjustable parameters and requires the geometries of conformers, their dipole moments and atomic polar tensor (APT) charges as the only input data. The fixed atomic charges generated by the method not only reproduce the molecular dipole moment in all the conformers accurately, but are also numerically close to the APT charges, thereby ensuring accurate reproduction of the dipole moment variations caused by small geometrical distortions (e.g., by vibrations) of the conformers. The proposed method has been applied to canonical 2'-deoxyribonucleotides, the model DNA monomers, and the dipole-derived charges have been shown to outperform both the averaged APT and RESP charges in reproducing the dipole moments of large sets of conformers, thus demonstrating a potential usefulness of the dipole-derived charges as a 'reference point' for modeling polarization effects in conformationally rich molecules.

Estimating the energy of intramolecular hydrogen bonds from 1H NMR and QTAIM calculations

The values of the downfield chemical shift of the bridge hydrogen atom were estimated for a series of compounds containing an intramolecular hydrogen bond O-HO, O-HN, O-HHal, N-HO, N-HN, C-HO, C-HN and C-HHal. Based on these values, the empirical estimation of the hydrogen bond energy was obtained by using known relationships. For the compounds containing an intramolecular hydrogen bond, the DFT B3LYP/6-311++G(d,p) method was used both for geometry optimization and for QTAIM calculations of the topological parameters (electron density ρ^BCP and the density of potential energy V in the critical point of the hydrogen bond). The calculated geometric and topological parameters of hydrogen bonds were also used to evaluate the energy of the hydrogen bond based on the equations from the literature. Comparison of calibrating energies from the ¹H NMR data with the energies predicted by calculations showed that the most reliable are the linear dependence on the topological ρ^BCP and V parameters. However, the correct prediction of the hydrogen bond energy is determined by proper fitting of the linear regression coefficients. To obtain them, new linear relationships were found between the calculated ρ^BCP and V parameters and the hydrogen bond energies obtained from empirical ¹H NMR data. These relationships allow the comparison of the energies of different types of hydrogen bonds for various molecules and biological ensembles.

Effective in silico prediction of new oxazolidinone antibiotics: force field simulations of the antibiotic-ribosome complex supervised by experiment and electronic structure methods

We propose several new and promising antibacterial agents for the treatment of serious Gram-positive infections. Our predictions rely on force field simulations, supervised by first principle calculations and available experimental data. Different force fields were tested in order to reproduce linezolid's conformational space in terms of a) the isolated and b) the ribosomal bound state. In a first step, an all-atom model of the bacterial ribosome consisting of nearly 1600 atoms was constructed and evaluated. The conformational space of 30 different ribosomal/oxazolidinone complexes was scanned by stochastic methods, followed by an evaluation of their enthalpic penalties or rewards and the mechanical strengths of the relevant hydrogen bonds (relaxed force constants; compliance constants). The protocol was able to reproduce the experimentally known enantioselectivity favoring the S-enantiomer. In a second step, the experimentally known MIC values of eight linezolid analogues were used in order to crosscheck the robustness of our model. In a final step, this benchmarking led to the prediction of several new and promising lead compounds. Synthesis and biological evaluation of the new compounds are on the way.

Behavior of Halogen Bonds of the Y-X⋅⋅⋅π Type (X, Y=F, Cl, Br, I) in the Benzene π System, Elucidated by Using a Quantum Theory of Atoms in Molecules Dual-Functional Analysis

The nature of halogen bonds of the Y-X-✶-π(C6 H6 ) type (X, Y=F, Cl, Br, and I) have been elucidated by using the quantum theory of atoms in molecules (QTAIM) dual-functional analysis (QTAIM-DFA), which we proposed recently. Asterisks (✶) emphasize the presence of bond-critical points (BCPs) in the interactions in question. Total electron energy densities, Hb (rc ), are plotted versus Hb (rc )-Vb (rc )/2 [=(ħ(2) /8m)∇(2) ρb (rc )] for the interactions in QTAIM-DFA, in which Vb (rc ) are potential energy densities at the BCPs. Data for perturbed structures around fully optimized structures were used for the plots, in addition to those of the fully optimized ones. The plots were analyzed by using the polar (R, θ) coordinate for the data of fully optimized structures with (θp , κp ) for those that contained the perturbed structures; θp corresponds to the tangent line of the plot and κp is the curvature. Whereas (R, θ) corresponds to the static nature, (θp , κp ) represents the dynamic nature of the interactions. All interactions in Y-X-✶-π(C6 H6 ) are classified by pure closed-shell interactions and characterized to have vdW nature, except for Y-I-✶-π(C6 H6 ) (Y=F, Cl, Br) and F-Br-✶-π(C6 H6 ), which have typical hydrogen-bond nature without covalency. I-I-✶-π(C6 H6 ) has a borderline nature between the two. Y-F-✶-π(C6 H6 ) (Y=Br, I) were optimized as bent forms, in which Y-✶-π interactions were detected. The Y-✶-π interactions in the bent forms are predicted to be substantially weaker than those in the linear F-Y-✶-π(C6 H6 ) forms.

Nucleotides containing variously modified sugars: energetics, structure, and mechanical properties

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.

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

Structure and conformational analysis of the anti-HIV reverse transcriptase inhibitor AZT using MP2 and DFT methods. Differences with the natural nucleoside thymidine. Simulation of the 1st phosphorylation step with ATP

A comprehensive quantum-chemical investigation of the conformational landscape of the HIV-1 reverse transcriptase inhibitor AZT (3'-azido-3'-deoxythymidine) nucleoside analogue was carried out. The whole conformational parameters (χ, γ, β, δ, ϕ, P, νmax) were analysed as well as the NBO charges. The search located at least 55 stable structures, 9 of which were by MP2 within a 1 kcal mol(-1) electronic energy range of the global minimum. Most conformers were anti or high-anti around the glycoside bond and with North sugar ring puckering angles. The distribution of all the conformers according to the ranges of stability of the characteristic torsional angles was established. The results obtained were in accordance with those found in related anti-HIV nucleoside analogues. The best conformer in the anti form corresponded to the calculated values by MP2 of χ = -126.9°, β = 176.4° and γ = 49.1°. An analysis of the lowest vibrations in conformer C1 was carried out. The first hydration shell was simulated and the structural differences with the natural nucleoside deoxythymidine (dT) were determined. The first phosphorylation step was simulated by interacting ATP with the best hydrated clusters of AZT and dT. The Na cations act as a bridge between the phosphate moieties of ATP making it easy for -P3O3 to receive the H5' proton from AZT or dT. A proton-transfer mechanism is proposed through the water molecules. When the number of the water molecules surrounding AZT is lower than 8, the first phosphorylation step of AZT can be carried out. However, the appropriate orientation of the O5'-H in dT avoids this limitation and it can be performed with large numbers of water molecules.

III-defined concepts in chemistry: rigid force constants vs. compliance constants as bond strength descriptors for the triple bond in diboryne

In a recent publication, the interpretation of Braunschweig's diboryne NHC-BB-NHC as a true triple bond is questioned. The analysis by Köppe and Schnöckel is based, inter alia, on the calculation of rigid coupling force constants. Nevertheless, since it is known for a long time that the use of rigid force constants as bond strength descriptors is by no means straightforward, we recomputed the rigid force constants for a model diboryne, applying different coordinate systems and compared the values with the relaxed force constants (generalized compliance constants, GCC). In contrast with the results by Schnöckel and Köppe, the true coupling between the boron-boron bond and the boron-carbon bond, that is, after the elimination of all numerical artifacts, is negligible (fBB/BC = -0.003).

How does the long G·G* Watson-Crick DNA base mispair comprising keto and enol tautomers of the guanine tautomerise? The results of a QM/QTAIM investigation

The double proton transfer (DPT) in the long G·G* Watson-Crick base mispair (|C6N1(G*)N1C6(G)| = 36.4°; C1 symmetry), involving keto and enol tautomers of the guanine (G) nucleobase, along two intermolecular neighboring O6H···O6 (8.39) and N1···HN1 (6.14 kcal mol(-1)) H-bonds that were established to be slightly anti-cooperative, leads to its transformation into the G*·G base mispair through a single transition state (|C6N1N1C6| = 37.1°; C1), namely to the interconversion into itself. It was shown that the G·G* ↔ G*·G tautomerisation via the DPT is assisted by the third specific contact, that sequentially switches along the intrinsic reaction coordinate (IRC) in an original way: (G)N2H···N2(G*) H-bond (-25.13 to -10.37) → N2···N2 van der Waals contact (-10.37 to -9.23) → (G)N2···HN2(G*) H-bond (-9.23 to 0.79) → (G*)N2···HN2(G) H-bond (0.79 to 7.35 Bohr). The DPT tautomerisation was found to proceed through the asynchronous concerted mechanism by employing the QM/QTAIM approach and the methodology of the scans of the geometric, electron-topological, energetic, polar and NBO properties along the IRC. Nine key points, that can be considered as part of the tautomerisation repertoire, have been established and analyzed in detail. Furthermore, it was shown that the G·G* or G*·G base mispair is a thermodynamically and dynamically stable structure with a lifetime of 8.22 × 10(-10) s and all 6 low-frequency intermolecular vibrations are able to develop during this time span. Lastly, our results highlight the importance of the G·G* ↔ G*·G DPT tautomerisation, which can have implications for biological and chemical sensing applications.

DPT tautomerisation of the G·A(syn) and A*·G*(syn) DNA mismatches: a QM/QTAIM combined atomistic investigation

By applying a combined QM and QTAIM atomistic computational approach we have established for the first time that the G·A(syn) and A*·G*(syn) DNA mismatches (rare tautomers are marked with an asterisk), causing spontaneous transversions with substantially various probabilities, radically differ from each other in their ability to tautomerise through the double proton transfer (DPT). The A*·G*(syn) mismatch tautomerises quite easily (ΔΔG(TS) ≈ 4·kT at room temperature) into the A·G*(syn) mismatch through the asynchronous concerted mechanism, whereas the G·A(syn) base mispair does not tautomerise via the DPT at all, since there is no local minimum corresponding to the tautomerised G*·A*(syn) mismatch on the potential energy surface. It was established that the A·G*(syn) base mispair is a dynamically unstable H-bonded complex with an extremely short lifetime of 2.17 × 10(-13) s. Consequently, the obtained results allow us to believe that spontaneous or forced dissociation of both the G·A(syn) and A*·G*(syn) DNA mismatches by the DNA-polymerase occurs with the preservation of the tautomeric status of the bases.

Direct and solvent-assisted keto-enol tautomerism and hydrogen-bonding interactions in 4-(m-chlorobenzylamino)-3-phenyl-4,5-dihydro-1H-1,2,4-triazol-5-one: a quantum-chemical study

The tautomeric equilibrium of the title triazole compound was computationally analyzed at the B3LYP/6-311++G(d,p) and MP2/6-311++G(d,p) levels of theory. The solvent effect was considered for three solvents (chloroform, methanol, and water). Two distinct mechanisms were applied: a direct intramolecular transfer using the polarizable continuum model (PCM) and a solvent-assisted mechanism. The calculations indicated that the keto form is more stable in all cases. It was found that the barrier heights for the tautomerization reaction are very high, indicating a relatively disfavored process. Although the barrier heights for solvent-assisted reactions are significantly lower than those for the unassisted tautomerization reaction, implying the importance of the superior catalytic effect of the solvents, monosolvation was not found to be sufficient for the reaction to occur. Finally, the two intermolecular hydrogen-bonding interactions in the crystal structure were investigated in the gas phase; according to the calculated energies and structural parameters, the order of stability is N3-H3···O1 > N1-H1···O1.

The nature of the transition mismatches with Watson-Crick architecture: the G*·T or G·T* DNA base mispair or both? A QM/QTAIM perspective for the biological problem

This study provides the first accurate investigation of the tautomerization of the biologically important guanine*·thymine (G*·T) DNA base mispair with Watson-Crick geometry, involving the enol mutagenic tautomer of the G and the keto tautomer of the T, into the G·T* mispair (∆G = .99 kcal mol(-1), population = 15.8% obtained at the MP2 level of quantum-mechanical theory in the continuum with ε = 4), formed by the keto tautomer of the G and the enol mutagenic tautomer of the T base, using DFT and MP2 methods in vacuum and in the weakly polar medium (ε = 4), characteristic for the hydrophobic interfaces of specific protein-nucleic acid interactions. We were first able to show that the G*·T↔G·T* tautomerization occurs through the asynchronous concerted double proton transfer along two antiparallel O6H···O4 and N1···HN3 H-bonds and is assisted by the third N2H···O2 H-bond, that exists along the entire reaction pathway. The obtained results indicate that the G·T* base mispair is stable from the thermodynamic point of view complex, while it is dynamically unstable structure in vacuum and dynamically stable structure in the continuum with ε = 4 with lifetime of 6.4·10(-12) s, that, on the one side, makes it possible to develop all six low-frequency intermolecular vibrations, but, on the other side, it is by three orders less than the time (several ns) required for the replication machinery to forcibly dissociate a base pair into the monomers during DNA replication. One of the more significant findings to emerge from this study is that the short-lived G·T* base mispair, which electronic interaction energy between the bases (-23.76 kcal mol(-1)) exceeds the analogical value for the G·C Watson-Crick nucleobase pair (-20.38 kcal mol(-1)), "escapes from the hands" of the DNA replication machinery by fast transforming into the G*·T mismatch playing an indirect role of its supplier during the DNA replication. So, exactly the G*·T mismatch was established to play the crucial role in the spontaneous point mutagenesis.

Atomistic nature of the DPT tautomerisation of the biologically important C·C* DNA base mispair containing amino and imino tautomers of cytosine: a QM and QTAIM approach

A theoretical study of tautomerisation of the biologically important cytosine·cytosine* (C·C*) DNA mismatch with a propeller-like structure (|C4N3N3C4| = 32.4°; C1 symmetry) and cis-oriented N1H glycosidic bonds, formed by the amino and imino tautomers of the C nucleobase, via the asynchronous concerted double proton transfer (DPT) along two H-bonds through the transition state (TSC·C*↔C*·C) (|C4N3N3C4| = 48.5°; C1 symmetry) into the C*·C mispair was carried out for the first time. It was established that the C·C*/C*·C DNA base mispair is associated by the antiparallel N4H···N4 (6.66 kcal mol(-1)), N3H···N3 (6.47 kcal mol(-1)) H-bonds and the O2···O2 van der Waals (vdW) contact (0.33 kcal mol(-1)), while the zwitterionic TSC·C*↔C*·C is stabilized by the parallel N4(+)H···N4(-) (13.55 kcal mol(-1)), N3(+)H···N3(-) (13.20 kcal mol(-1)) H-bonds and the O2(+)···O2(-) vdW contact (0.60 kcal mol(-1)). It was shown that the C·C* ↔ C*·C tautomerisation via the DPT is assisted by the O2···O2 vdW contact, that in contrast to the two others N4H···N4 and N3H···N3 H-bonds exists along the entire intrinsic reaction coordinate (IRC) range. The positive values of the Grunenberg's compliance constants (30.919 and 21.384 Å mdyn(-1) for C·C*/C*·C and TSC·C*↔C*·C, respectively) indicate that the O2···O2 vdW contact is a stabilizing closed-shell interaction. It was found that the middle N3H···N3 H-bond is anti-cooperative with the upper N4H···N4 H-bond and cooperative with the lower O2···O2 vdW contact. The 9 key points, which can be considered as electron-topological "fingerprints" of the asynchronous concerted C·C* ↔ C*·C tautomerisation process via the DPT were revealed along the IRC and examined in detail. It was shown that the C·C*/C*·C base mispair is a thermodynamically and dynamically stable structure. Its lifetime is equal to 1.53 × 10(-7) s at the MP2/cc-pVQZ//B3LYP/6-311++G(d,p) level of theory in vacuum. All 6 low-frequency intermolecular vibrations are able to develop during this time span.

DPT tautomerisation of the wobble guanine·thymine DNA base mispair is not mutagenic: QM and QTAIM arguments

We have shown for the first time, connecting QM methods with QTAIM analysis and using the methodology of the sweeps of the energetical, electron-topological and geometrical parameters, that the tautomerisation of the wobble guanine·thymine (wG·T) DNA base mispair into the wG(*)·T(*) base mispair induced by the double proton transfer (DPT), which undergoes a concerted asynchronous pathway, is not mutagenic. The wG·T → wG(*)·T(*) DPT tautomerisation does not result in the transition of the G base into its mutagenic tautomeric form G(*) able to mispair with the T base within the Watson-Crick base pairing scheme. This observation is explained by the so-called quantum protection of the wG·T DNA base mispair from its mutagenic tautomerisation - the dynamical non-stability of the tautomerised wG(*)·T(*) base mispair and significantly negative value of the Gibbs free energy of activation for the reverse reaction of the wG·T → wG(*)·T(*) DPT tautomerisation.

Is the DPT tautomerization of the long A·G Watson-Crick DNA base mispair a source of the adenine and guanine mutagenic tautomers? A QM and QTAIM response to the biologically important question

Herein, we first address the question posed in the title by establishing the tautomerization trajectory via the double proton transfer of the adenine·guanine (A·G) DNA base mispair formed by the canonical tautomers of the A and G bases into the A*·G* DNA base mispair, involving mutagenic tautomers, with the use of the quantum-mechanical calculations and quantum theory of atoms in molecules (QTAIM). It was detected that the A·G ↔ A*·G* tautomerization proceeds through the asynchronous concerted mechanism. It was revealed that the A·G base mispair is stabilized by the N6H···O6 (5.68) and N1H···N1 (6.51) hydrogen bonds (H-bonds) and the N2H···HC2 dihydrogen bond (DH-bond) (0.68 kcal·mol(-1) ), whereas the A*·G* base mispair-by the O6H···N6 (10.88), N1H···N1 (7.01) and C2H···N2 H-bonds (0.42 kcal·mol(-1) ). The N2H···HC2 DH-bond smoothly and without bifurcation transforms into the C2H···N2 H-bond at the IRC = -10.07 Bohr in the course of the A·G ↔ A*·G* tautomerization. Using the sweeps of the energies of the intermolecular H-bonds, it was observed that the N6H···O6 H-bond is anticooperative to the two others-N1H···N1 and N2H···HC2 in the A·G base mispair, while the latters are significantly cooperative, mutually strengthening each other. In opposite, all three O6H···N6, N1H···N1, and C2H···N2 H-bonds are cooperative in the A*·G* base mispair. All in all, we established the dynamical instability of the А*·G* base mispair with a short lifetime (4.83·10(-14) s), enabling it not to be deemed feasible source of the A* and G* mutagenic tautomers of the DNA bases. The small lifetime of the А*·G* base mispair is predetermined by the negative value of the Gibbs free energy for the A*·G* → A·G transition. Moreover, all of the six low-frequency intermolecular vibrations cannot develop during this lifetime that additionally confirms the aforementioned results. Thus, the A*·G* base mispair cannot be considered as a source of the mutagenic tautomers of the DNA bases, as the A·G base mispair dissociates during DNA replication exceptionally into the A and G monomers in the canonical tautomeric form.

Atomistic understanding of the C·T mismatched DNA base pair tautomerization via the DPT: QM and QTAIM computational approaches

It was established that the cytosine·thymine (C·T) mismatched DNA base pair with cis-oriented N1H glycosidic bonds has propeller-like structure (|N3C4C4N3| = 38.4°), which is stabilized by three specific intermolecular interactions-two antiparallel N4H…O4 (5.19 kcal mol(-1)) and N3H…N3 (6.33 kcal mol(-1)) H-bonds and a van der Waals (vdW) contact O2…O2 (0.32 kcal mol(-1)). The C·T base mispair is thermodynamically stable structure (ΔG(int) = -1.54 kcal mol(-1) ) and even slightly more stable than the A·T Watson-Crick DNA base pair (ΔG(int) = -1.43 kcal mol(-1)) at the room temperature. It was shown that the C·T ↔ C*·T* tautomerization via the double proton transfer (DPT) is assisted by the O2…O2 vdW contact along the entire range of the intrinsic reaction coordinate (IRC). The positive value of the Grunenberg's compliance constants (31.186, 30.265, and 22.166 Å/mdyn for the C·T, C*·T*, and TS(C·T ↔ C*·T*), respectively) proves that the O2…O2 vdW contact is a stabilizing interaction. Based on the sweeps of the H-bond energies, it was found that the N4H…O4/O4H…N4, and N3H…N3 H-bonds in the C·T and C*·T* base pairs are anticooperative and weaken each other, whereas the middle N3H…N3 H-bond and the O2…O2 vdW contact are cooperative and mutually reinforce each other. It was found that the tautomerization of the C·T base mispair through the DPT is concerted and asynchronous reaction that proceeds via the TS(C·T ↔ C*·T*) stabilized by the loosened N4-H-O4 covalent bridge, N3H…N3 H-bond (9.67 kcal mol(-1) ) and O2…O2 vdW contact (0.41 kcal mol(-1) ). The nine key points, describing the evolution of the C·T ↔ C*·T* tautomerization via the DPT, were detected and completely investigated along the IRC. The C*·T* mispair was revealed to be the dynamically unstable structure with a lifetime 2.13·× 10(-13) s. In this case, as for the A·T Watson-Crick DNA base pair, activates the mechanism of the quantum protection of the C·T DNA base mispair from its spontaneous mutagenic tautomerization through the DPT.