Hydrogen bonding interactions in noradrenaline-DMSO complexes: DFT and QTAIM studies of structure,properties and topology

The hydrogen bonding interactions between noradrenaline (NA) and DMSO were studied with density functional theory (DFT) regarding their geometries, energies, vibrational frequencies, and topological features of the electron density. The quantum theory of atoms in molecules (QTAIM) and the natural bond orbital (NBO) analyses were employed to elucidate the hydrogen bonding interaction characteristics in noradrenaline-DMSO complexes. The H-bonds involving the hydroxyls hydrogen in NA and the O atom in DMSO are dominant intermolecular H-bonds and are stronger than other H-bonds involving the methyl hydrogen of DMSO as a H-donor. The weak H-bonds also include a π H-bond which involves the benzene ring as a H-donor or H-acceptor. QTAIM identified the weak H-bonds formed between the methyl hydrogen of DMSO and the N atom in NA in some complexes (AB5, AB6 and AB7), which cannot be further confirmed by NBO and other methods, so there are probably no interactions between hydrogen and nitrogen atoms among these complexes. A good linear relationship between logarithmic electron density (lnρ _b ) at the bond critical point (BCP) and structural parameter (δR _H···Y) was found. The formations of new H-bonds in some complexes are helpful to strengthen the original intramolecular H-bond, this is attributed to the cooperativity of H-bonds in complexes and can be learned from the structure results and the NBO and QTAIM analyses. Analysis of various physically meaningful contributions arising from the energy decomposition procedures show that the orbital interactions of H-bond is predominant during the formation of the complex, moreover, both the hydrogen bonding interaction and the structural deformation are responsible for the stability of the complexes.     

A theoretical investigation on the conformation and the interaction of CHF2OCF2CHF2 (desflurane II) with one water molecule

The conformation and the interaction of CHF₂OCF₂CHF₂ (desflurane II) with one water molecule is investigated theoretically using the ab initio MP2/aug-cc-pvdz and DFT-based M062X/6-311++G(d,p) methods. The calculations include the optimized geometries, the harmonic frequencies of relevant vibrational modes along with a natural bond orbital (NBO) analysis including the NBO charges, the hybridization of the C atom and the intra- and intermolecular hyperconjugation energies. In the two most stable conformers, the CH bond of the F₂HCO- group occupies the gauche position. The hyperconjugation energies are about the same for both conformers and the conformational preference depends on the interaction between the non-bonded F and H atoms. The deprotonation enthalpies of the CH bonds are about the same for both conformers, the proton affinity of the less stable conformer being 3 kcal mol^?1 higher. Both conformers of desflurane II interact with water forming cyclic complexes characterized by CH…O and OH…F hydrogen bonds. The binding energies are moderate, ranging from ?2.4 to ?3.2 kcal mol^?1 at the MP2 level. The origin of the blue shifts of the ν(CH) vibrations is analyzed. In three of the complexes, the water molecule acts as an electron donor. Interestingly, in these cases a charge transfer is also directed to the non bonded OH group of the water molecule. This effect seems to be a property of polyfluorinated ethers.     

Theoretical study of the interaction between 5-methylcytosine and acrylamide

The hydrogen-bonded complexes between 5-methylcytosine and acrylamide have been investigated using the density function theory (DFT) method. Five stable complexes have been found with no imaginary frequencies. Complex C3 is the most stable one with interaction energies of -69.01?kJ?mol(-1) corrected for basis set superposition error (BSSE). The charge change in the process of these complexes formation has also been examined. The atoms in molecules (AIM) theory and natural bond orbital (NBO) method have been performed to investigate the hydrogen bonds involved in all the complexes. The electron density and its corresponding Laplacian at the bond and ring critical points have been analyzed. In C3 complex, there is the largest stabilization energy (18.17?kJ?mol(-1)) between N11-H12 antibonding orbital and lone electron pair of O17. It can be seen that the hydrogen bonds play a crucial role in the stability of all the complexes between 5-methylcytosine and acrylamide. The theoretical results could provide helpful information for other researchers in further work.     

Cooperativity between fluorine-centered halogen bonds: investigation of substituent effects

This article analyzes the substitution effects on cooperativity between fluorin-centered halogen bonds in NCF?·?·?·?NCF?·?·?·?NCX and CNF?·?·?·?CNF?·?·?·?CNX complexes, where X?=?H, F, Cl, CN, OH, and NH₂. These effects are investigated theoretically in terms of geometric and energetic features of the complexes, which are computed by ab initio methods. The topological analysis, based on the quantum theory of atoms in molecules (QTAIM), is used to characterize the interactions and analyze their enhancement with varying electron density at bond critical points. It is found that the complexes with electron-donating groups exhibit a strong cooperativity, while a much weaker cooperativity occurs in the NCF?·?·?·?NCF?·?·?·?NCCN and CNF?·?·?·?CNF?·?·?·?CNCN trimers. An excellent correlation is found between the cooperative energy in the ternary complexes and the calculated three-body interaction energies. The energy decomposition analysis (EDA) indicates that the electrostatic and dispersion effects play a main role in the cooperativity of fluorine-centered halogen bonding.

Effect of stepwise microhydration on the methylammonium···phenol and ammonium···phenol interaction

A computational study has been performed for studying the characteristics of the interaction of phenol with ammonium and methylammonium cations. The effect of the presence of water molecules has also been considered by microhydrating the clusters with up to three water molecules. Clusters of phenol with ammonium and methylammonium cations present similar characteristics, though ammonium complexes have been found to be more stable than the methylammonium ones. The first water molecule included in the complexes interacts with a N-H group of ammoniun cations and simultaneously with the hydroxyl oxygen atom of phenol (or the aromatic ring). This first water molecule is more tightly bound in the complex, so the stability gain as more water molecules are included drops significantly by 2-3 kcal?mol^?1 with respect to the first one. As more water molecules are included, the differences between favorable coordination sites (the cation, the hydroxyl group or a previous water molecule) decrease. As a consequence, several of the most stable complexes located including three water molecules already exhibit hydrogen bonds between the hydroxyl group and one water molecule. The results indicate that a cyclic pattern formed by a series of hydrogen bonds: π···H-N-H···O-H···O-?, is characteristic of the most stable minima, being kept as more water molecules are included in the system. Therefore, this pattern can be expected to be crucial in ammonium cations···phenol interaction if exposed to the solvent to any degree.     

Enhancing the hydrogen bond between the bridged hydrogen atom of diborane and ammonia

The character of the bridged hydrogen atom (H_b) of B₂H₆ has become a hot issue in recent years. In this work, the complexes B₂H₆?·?·?·?NH₃, B₂H₂X₄?·?·?·?nNH₃ (n?=?1, 2) and 2HF?·?·?·?B₂H₂X₄?·?·?·?2NH₃ (X?=?Cl, Br, I) were constructed and studied based on the M06-2X calculations to investigate how to enhance the H_b?·?·?·?N hydrogen-bonded interaction. When the terminal hydrogen atoms (H_t) of B₂H₆ were replaced by X (X?=?Cl, Br, I) atoms, the H_b?·?·?·?N hydrogen bond were strengthened. According to the electrostatic potentials in B₂H₂X₄, two HF molecules were added to the interspace of the B-H-B-H four-membered ring of the B₂H₂X₄?·?·?·?2NH₃ complexes, and H?·?·?·?X hydrogen bond formed, resulting in further enhancing effect of H_b?·?·?·?N hydrogen bond. As a result, the positive cooperative effect of H_b?·?·?·?N hydrogen bond and H?·?·?·?X hydrogen bond do enhance the interactions of each other. The two measures not only enhance the strength of H_b?·?·?·?N hydrogen bond, but also achieve the goal to make the H_b?·?·?·?N hydrogen bond perpendicular to B?·?·?·?B direction.

Open image in new window

Graphical Abstract Enhancing the hydrogen bond between the bridged hydrogen atom of diborane and ammonia?

     

Theoretical studies of the interaction between enflurane and water

Increase of the atmospheric concentration of halogenated organic compounds is partially responsible for a change of the global climate. In this work we have investigated the interaction between halogenated ether and water, which is one of the most important constituent of the atmosphere. The structures of the complexes formed by the two most stable conformers of enflurane (a volatile anaesthetic) with one and two water molecules were calculated by means of the counterpoise CP-corrected gradient optimization at the MP2/6–311++G(d,p) level. In these complexes the CH…O_w hydrogen bonds are formed, with the H…O_w distances varying between 2.23 and 2.32 Å. A small contraction of the CH bonds and the blue shifts of the ν(CH) stretching vibrations are predicted. There is also a weak interaction between one of the F atoms and the H atom of water, with the H_w…F distances between 2.41 and 2.87 Å. The CCSD(T)/CBS calculated stabilization energies in these complexes are between ?5.89 and ?4.66 kcal?mol^?1, while the enthalpies of formation are between ?4.35 and ?3.22 kcal?mol^?1. The Cl halogen bonding between enflurane and water has been found in two complexes. The intermolecular (Cl···O) distance is smaller than the sum of the corresponding van der Waals radii. The CCSD(T)/CBS stabilization energies for these complexes are about ?2 kcal?mol^?1.

The effects of tautomerization and protonation on the adenine–cytosine mismatches: a density functional theory study

In the present work, we demonstrate the results of a theoretical study concerned with the question how tautomerization and protonation of adenine affect the various properties of adenine–cytosine mismatches. The calculations, in gas phase and in water, are performed at B3LYP/6-311++G(d,p) level. In gas phase, it is observed that any tautomeric form of investigated mismatches is more stabilized when adenine is protonated. As for the neutral mismatches, the mismatches containing amino form of cytosine and imino form of protonated adenine are more stable. The role of aromaticity on the stability of tautomeric forms of mismatches is investigated by NICS(1)_ZZ index. The stability of mispairs decreases by going from gas phase to water. It can be explained using dipole moment parameter. The influence of hydrogen bonds on the stability of mismatches is examined by atoms in molecules and natural bond orbital analyses. In addition to geometrical parameters and binding energies, the study of the topological properties of electron charge density aids in better understanding of these mispairs.     

Ab Initio Quantum Mechanics Analysis of Imidazole C-H…O Water Hydrogen Bonding and a Molecular Mechanics Forcefield Correction

Abstract

While it is well established that classical hydrogen bonds play an important role in enzyme structure, function and dynamics, the role of weaker, but ‘activated’ C-H donor hydrogen bonds is poorly understood. The most important such case involves histidine which often plays a direct role in enzyme catalysis and possesses the most acidic C-H donor group of the standard amino acids. In the present study, we obtained optimized geometries and hydrogen bond interaction energies for C-H…O hydrogen bonded complexes between methane, ethylene, benzene, acetylene, and imidazole with water at the MP2-FC/6-31++G(2d,2p) and MP2-FC/aug-cc-pVDZ//MP2-FC/6-31++G(2d,2p) levels of theory. A strong linear relationship is obtained between the stability of the various hydrogen bonded complexes and both separation distances for H…0 and C—O. In general, these calculations indicate that C-H…0 interactions can be classified as hydrogen bonding interactions, albeit significantly weaker than the classical hydrogen bonds, but significantly stronger than just van der Waals interactions. For instance, while the electronic energy of stabilization at the MP2-FC/aug-cc-pVDZ//MP2-FC/6-31++G(2d,2p) level of theory of a water C-H…O water hydrogen bond is 4.36 kcal/mol more stable than the methane C-H…O water interaction, the water-water hydrogen bond is only 2.06 kcal/mol more stable than the imidazole C_e?H…O water hydrogen bond. Neglecting this latter hydrogen bonding interaction is obviously unacceptable. We next compare the potential energy surfaces for the imidazole C_e?H…O water and imidazole N_d?H…O hydrogen bonded complexes computed at the MP2/6-31++G(2d,2p) level of theory with the potential energy surface computed using the AMBER molecular mechanics program and forcefields. While the Weiner et al and Cornell et al AMBER forcefields reasonably account for the imidazole N-H…O water interaction, these forcefields do not adequately account for the imidazole C_e?H…O water hydrogen bond. A forcefield modification is offered that results in excellent agreement between the ab initio and molecular mechanics geometry and energy for this C-H…O hydrogen bonded complex.     

A density functional theory study on the hydrogen bonding interactions between luteolin and ethanol

Ethanol is one of the most commonly used solvents to extract flavonoids from propolis. Hydrogen bonding interactions play an important role in the properties of liquid system. The main objective of the work is to study the hydrogen bonding interactions between flavonoid and ethanol. Luteolin is a very common flavonoid that has been found in different geographical and botanical propolis. In this work, it was selected as the representative flavonoid to do detailed research. The study was performed from a theoretical perspective using density functional theory (DFT) method. After careful optimization, there exist nine optimized geometries for the luteolin ? CH₃CH₂OH complex. The binding distance of X ? H···O, and the bond length, vibrational frequency, and electron density changes of X ? H all indicate the formation of the hydrogen bond in the optimized geometries. In the optimized geometries, it is found that: (1) except for the H2’, H5’, and H6’, CH₃CH₂OH has formed hydrogen bonds with all the hydrogen and oxygen atoms in luteolin. The hydrogen atoms in the hydroxyl groups of luteolin form the strongest hydrogen bonds with CH₃CH₂OH; (2) all of the hydrogen bonds are closed-shell interactions; (3) the strongest hydrogen bond is the O3’ ? H3’···O in structure A, while the weakest one is the C3 ? H3···O in structure E; (4) the hydrogen bonds of O3’ ? H3’···O, O ? H···O4, O ? H···O3’ and O ? H···O7 are medium strength and covalent dominant in nature. While the other hydrogen bonds are weak strength and possess a dominant character of the electrostatic interactions in nature.     

DFT study of water adsorption on lignite molecule surface

High moisture content is a main characteristic of low-rank coal, such as lignite. Numerous oxygen containing functional groups in lignite make it represent some special properties, and these functional groups affect the adsorption mechanisms of water molecules on lignite surface. This study reports some typical water?·?·?·?lignite conformations, along with a detailed analysis of the geometry, electrostatic potential distribution, reduced density gradient of interaction, and interaction energy decomposition. The results show that water molecules tend to aggregate around functional groups, and hydrogen bonds play a dominant role in the interaction. The adsorption energy of water cluster on lignite surface is larger than that of isolated water molecule, a good linear relationship between the interaction distance and adsorption energy of layers has been found. Since water is a polar molecule, the local minima and maxima of electrostatic potential in conformations increase along with more water adsorbing on lignite surface. Reduced density gradient analysis shows that H-bonds, van der Waals interaction, and a little steric make up the interaction between water cluster and lignite molecule. In these studied conformations which mainly are H-bond complexes, electrostatic and exchange repulsion play a dominant role, whereas polarization and dispersion make relatively small contribution to the interaction. Attractive and repulsive interaction both affect the stability of water?·?·?·?lignite conformations.     

Calixarene building block bis(2-hydroxyphenyl)methane (2HDPM) and hydrogen-bonded 2HDPM-H2O complex in electronic excited state

Intramolecular and intermolecular hydrogen bonding in electronic excited states of calixarene building blocks bis(2-hydroxyphenyl)methane (2HDPM) monomer and hydrogen-bonded 2HDPM-H₂O complex were studied theoretically using the time-dependent density functional theory (TDDFT). Twenty-four stable conformations (12 pairs of enantiomers) of 2HDPM monomer have been found in the ground state. From the calculation results, the conformations 1a and 1b which both have an intramolecular hydrogen bond are the most stable ones. The infrared spectra of 2HDPM monomer and 2HDPM-H₂O complex in ground state and S₁ state were calculated. The stretching vibrational absorption band of O₂???H₃ group in the monomer and complex disappeared in the S₁ state. At the same time, a new strong absorption band appeared at the C=O stretching region. From the calculation of bond lengths, it indicates that the O₂???H₃ bond is significantly lengthened in the S₁ state. However, the C₁???O₂ bond is drastically shortened upon electronic excitation to the S₁ state and has the characteristics of C=O band. Furthermore, the intramolecular hydrogen bond O₂???H₃?·?·?·?O₄ of the 2HDPM monomer and the intermolecular hydrogen bonds O₂???H₃?·?·?·?O₇ and O₇???H₉?·?·?·?O₄ of 2HDPM-H₂O complex are all shortened and strengthened in the S₁ state.

Competition and interplay between the lithium bonding and hydrogen bonding: R3C···HY···LiY and R3C···LiY···HY triads as a working model (R=H, CH3; Y=CN, NC)

UMP2 calculations with aug-cc-pVDZ basis set were used to analyze intermolecular interactions in R₃C···HY···LiY and R₃C···LiY···HY triads (R=H, CH₃; Y=CN, NC), which are connected via lithium and hydrogen bonds. To better understand the properties of these systems, the corresponding dyads were also studied. Molecular geometries and binding energies of dyads, and triads were investigated at the UMP2/aug-cc-pVDZ computational level. Particular attention was paid to parameters such as cooperative energies, and many-body interaction energies. All studied complexes, with the simultaneous presence of a lithium bond and a hydrogen bond, showed cooperativity with energy values ranging between ?1.71 and ?9.03 kJ mol^?1. The electronic properties of the complexes were analyzed using parameters derived from atoms in molecules (AIM) methodology. Energy decomposition analysis revealed that the electrostatic interactions are the major source of the attraction in the title complexes.     

Intermolecular interactions between gold clusters and selected amino acids cysteine and glycine: a DFT study

The intermolecular interactions between Au_n (n = 3–4) clusters and selected amino acids cysteine and glycine have been investigated by means of density functional theory (DFT). Present calculations show that the complexes possessing Au-NH₂ anchoring bond are found to be energetically favored. The results of NBO and frontier molecular orbitals analysis indicate that for the complex with anchoring bonds, lone pair electrons of sulfur, oxygen, and nitrogen atoms are transferred to the antibonding orbitals of gold, while for the complex with the nonconventional hydrogen bonds (Au···H–O), the lone pair electrons of gold are transferred to the antibonding orbitals of O-H bonds during the interaction. Furthermore, the interaction energy calculations show that the complexes with Au-NH₂ anchoring bond have relatively high intermolecular interaction energy, which is consistent with previous computational studies.     

Theoretical investigations on the hydrogen bonding of nitrile isomers with H2O,HF, NH3 and H2S

The hydrogen bonds formed by the interaction of nitriles with water, hydrogen fluoride, ammonia and hydrogen sulphide have been studied using B3LYP and second-order Møller–Plesset perturbation (MP2) methods and 6-311+ + G(d,p) basis set. The energies and structures of 80 hydrogen-bonded complexes between nitriles and small molecules were examined systematically using B3LYP and MP2 procedure. Categorisation of the hydrogen bonds involved in the various complexes led to an ordering of hydrogen bond donor and acceptor abilities for some functional groups. The interaction energies have been corrected for the basis set superposition error using Boy's counterpoise correction method. The Morokuma energy decomposition analysis reveals that the strong interactions are due to the attractive contributions from the electrostatic (ES), polarisation (PL) and charge transfer (CT) components. The topological parameters, electron density and Laplacian of electron density show excellent correlation with the hydrogen bond length. Natural bond orbital (NBO) analysis has also been performed to study the CT from proton acceptor to the antibonding orbital of the H–Y bond in the proton donor part of complexes. The frequency analysis of C–H…Y bond in the complexes indicates the blue-shifting nature largely in case of sp² hybridised carbon atom.     

Effects of ionization on stability of 1-methylcytosine — DFT and PCM studies

This report present the results of natural energy decomposition analysis (NEDA), natural bond orbital (NBO), and quantum theory of atoms in molecules (QTAIM) calculations of three derivatives of biphenyl-1-aza-18-crown-6 ether and their 1:1 complexes with Cd²⁺. All calculations used the B3LYP density functional theory in combination with the 6-311G and WTBS basis sets for ligands and Cd²⁺ ion, respectively. Ligands 1 and 3 have a single 1-aza-18-crown-6, substituent; ligand 2 has two such substituents. The results show that, in the optimized geometries of the complexes, the distance between N and Cd²⁺ is greater than the distance between O and Cd²⁺. NBO and QTAIM data confirm these results. There was no stabilization energy or bond critical point for N · · · Cd²⁺ in NBO or QTAIM, respectively. Data show that the O · · · Cd²⁺ interaction is a kind of closed shell interaction. The trend of the calculated stabilization energy was similar to the experimental data. Different contributions of interaction energies for complex formation were analyzed by NEDA, and the results show that the main component of the interactions is accounted for by polarization.

     

Theoretical investigation on the structure and thermodynamic properties of the 2,4-dinitroimidazole complex with methanol

The structure and thermodynamic properties of the 2, 4-dinitroimidazole complex with methanol were investigated using the B3LYP and MP2(full) methods with the 6-31++G(2d,p) and 6-311++G(3df,2p) basis sets. Four types of hydrogen bonds [N–H?O, C–H?O, O–H?O (nitro oxygen) and O–H?π] were found. The hydrogen-bonded complex having the highest binding energy had a N–H?O hydrogen bond. Analyses of natural bond orbital (NBO) and atoms-in-molecules (AIM) revealed the nature of the intermolecular hydrogen-binding interaction. The changes in thermodynamic properties from monomers to complexes with temperatures ranging from 200.0 to 800.0 K were investigated using the statistical thermodynamic method. Hydrogen-bonded complexes of 2,4-dinitroimidazole with methanol are fostered by low temperatures.

Interplay between halogen bonds and hydrogen bonds in OH/SH···HOX···HY (X = Cl, Br; Y = F, Cl, Br) complexes

The character of the cooperativity between the HOX···OH/SH halogen bond (XB) and the Y―H···(H)OX hydrogen bond (HB) in OH/SH···HOX···HY (X = Cl, Br; Y = F, Cl, Br) complexes has been investigated by means of second-order Møller?Plesset perturbation theory (MP2) calculations and “quantum theory of atoms in molecules” (QTAIM) studies. The geometries of the complexes have been determined from the most negative electrostatic potentials (V _S,min) and the most positive electrostatic potentials (V _S,max) on the electron density contours of the individual species. The greater the V _S,max values of HY, the larger the interaction energies of halogen-bonded HOX···OH/SH in the termolecular complexes, indicating that the ability of cooperative effect of hydrogen bond on halogen bond are determined by V _S,max of HY. The interaction energies, binding distances, infrared vibrational frequencies, and electron densities ρ at the BCPs of the hydrogen bonds and halogen bonds prove that there is positive cooperativity between these bonds. The potentiation of hydrogen bonds on halogen bonds is greater than that of halogen bonds on hydrogen bonds. QTAIM studies have shown that the halogen bonds and hydrogen bonds are closed-shell noncovalent interactions, and both have greater electrostatic character in the termolecular species compared with the bimolecular species.

Synthesis,characterization, DNA interaction,and cytotoxicity of novel Pd(II) and Pt(II) complexes

Four complexes [Pd(L)(bipy)Cl]·4H₂O (1), [Pd(L)(phen)Cl]·4H₂O (2), [Pt(L)(bipy)Cl]·4H₂O (3), and [Pt(L)(phen)Cl]·4H₂O (4), where L = quinolinic acid, bipy = 2,2’-bipyridyl, and phen = 1,10-phenanthroline, have been synthesized and characterized using IR, ¹H NMR, elemental analysis, and single-crystal X-ray diffractometry. The binding of the complexes to FS-DNA was investigated by electronic absorption titration and fluorescence spectroscopy. The results indicate that the complexes bind to FS-DNA in an intercalative mode and the intrinsic binding constants K of the title complexes with FS-DNA are about 3.5?×?10⁴ M^?1, 3.9?×?10⁴ M^?1, 6.1?×?10⁴ M^?1, and 1.4?×?10⁵ M^?1, respectively. Also, the four complexes bind to DNA with different binding affinities, in descending order: complex 4, complex 3, complex 2, complex 1. Gel electrophoresis assay demonstrated the ability of the Pt(II) complexes to cleave pBR322 plasmid DNA.     

Prototropic tautomerism and basic molecular principles of hypoxanthine mutagenicity: an exhaustive quantum-chemical analysis

The molecular structures, relative stability order, and dipole moments of a complete family of 21 planar hypoxanthine (Hyp) prototropic molecular–zwitterionic tautomers including ylidic forms were computationally investigated at the MP2/6–311++G(2df,pd)//B3LYP/6–311++G(d,p) level of theory in vacuum and in three different surrounding environments: continuum with a low dielectric constant (??=?4) corresponding to a hydrophobic interface of protein–nucleic acid interactions, dimethylsulfoxide (DMSO), and water. The keto-N1HN7H tautomer was established to be the global minimum in vacuum and in continuum with ??=?4, while Hyp molecule exists as a mixture of the keto-N1HN9H and keto-N1HN7H tautomers in approximately equal amounts in DMSO and in water at T?=?298.15?K. We found out that neither intramolecular tautomerization by single proton transfer in the Hyp base, nor intermolecular tautomerization by double proton transfer in the most energetically favorable Hyp·Hyp homodimer (symmetry C _2h ), stabilized by two equivalent N1H…O6 H-bonds, induces the formation of the enol tautomer (marked with an asterisk) of Hyp with cis-oriented O6H hydroxyl group relative to neighboring N1C6 bond. We first discovered a new scenario of the keto–enol tautomerization of Hyp?·?Hyp homodimer (C _2h ) via zwitterionic near-orthogonal transition state (TS), stabilized by N1⁺H…N1^? and O6⁺H…N1^? H-bonds, to heterodimer Hyp^??·?Hyp (C _s ), stabilized by O6H…O6 and N1H…N1 H-bonds. We first showed that Hyp^??·?Thy mispair (C _s ), stabilized by O6H…O4, N3H…N1, and C2H…O2 H-bonds, mimicking Watson–Crick base pairing, converts to the wobble Hyp?·?Thy base pair (C _s ), stabilized by N3H…O6 and N1H…O2 H-bonds, via high- and low-energy TSs and intermediate Hyp?·?Thy^?, stabilized by O4H…O6, N1H…N3, and C2H…O2 H-bonds. The most energetically favorable TS is the zwitterionic pair Hyp⁺?·?Thy^? (C _s ), stabilized by O6⁺H…O4^?, O6⁺H…N3^?, N1⁺H…N3^?, and N1⁺H…O2^? H-bonds. The authors expressed and substantiated the hypothesis, that the keto tautomer of Hyp is a mutagenic compound, while enol tautomer Hyp^? does not possess mutagenic properties. The lifetime of the nonmutagenic tautomer Hyp^? exceeds by many orders the time needed to complete a round of DNA replication in the cell. For the first time purine–purine planar H-bonded mispairs containing Hyp in the anti-orientation with respect to the sugar moiety – Hyp?·?Ade _syn , Hyp?·?Gua^? _syn , and Hyp?·?Gua _syn , that closely resembles the geometry of the Watson–Crick base pairs, have been suggested as the source of transversions. An influence of the surrounding environment (??=?4) on the stability of studied complexes and corresponding TSs was estimated by means of the conductor-like polarizable continuum model. Electron-topological, structural, vibrational, and energetic characterictics of all conventional and nonconventional H-bonds in the investigated structures are presented. Presented data are key to understanding elementary molecular mechanisms of mutagenic action of Hyp as a product of the adenine deamination in DNA.