Bridging QTAIM with vibrational spectroscopy: the energy of intramolecular hydrogen bonds in DNA-related biomolecules

Physical properties of over 8000 intramolecular hydrogen bonds (iHBs), including 2901 ones of the types OH···O, OH···N, NH···O and OH···C, in 4244 conformers of the DNA-related molecules (four canonical 2'-deoxyribonucleotides, 1,2-dideoxyribose-5-phosphate, and 2-deoxy-D-ribose in its furanose, pyranose and linear forms) have been investigated using quantum theory of atoms in molecules (QTAIM) and vibrational analysis. It has been found that for all iHBs with positive red-shift of the proton donating group stretching frequency the shift value correlates with ρ(cp)-the electron charge density at the (3,-1)-type bond critical point. Combining QTAIM and spectroscopic data new relationships for estimation of OH···O, OH···N, NH···O and OH···C iHB enthalpy of formation (kcal mol(-1)) with RMS error below 0.8 kcal mol(-1) have been established: E(OH···O) = -3.09 + 239·ρ(cp), E(OH···N) = 1.72 + 142·ρ(cp), E(NH···O) = -2.03 + 225·ρ(cp), E(OH···C) = -0.29 + 288·ρ(cp), where ρ(cp) is in e a(0)(-3) (a(0)- the Bohr radius). It has been shown that XHY iHBs with red-shift values over 40 cm(-1) are characterized by the following minimal values of the XHY angle, ρ(cp) and nubla(2)ρ(cp): 112°, 0.005 e a(0)(-3) and 0.016 e a(0)(-5), respectively. New relationships have been used to reveal the strongest iHBs in canonical 2'-deoxy- and ribonucleosides and the O(5')H···N(3) H-bond in ribonucleoside guanosine was found to have the maximum energy (8.1 kcal mol(-1)).     

Effect of solvation on ion binding to imidazole and methylimidazole

Quantum chemical [MP2(FULL)/6-311++G-(d,p)] calculations are done on the binding of hydrated Li(+), Na(+), K(+), Mg(2+), Cu(+), and Zn(2+) metal ions with biologically relevant heteroaromatics such as imidazole and methylimidazole. The computed interaction energies are found to be in good agreement with the available experimental data. The effect of hydration on hydrogen bonding has been studied in detail and it shows that the hydrogen bond strength between H(2)O···H-N(1) substantially increases in the presence of metal ions. The present study quantifies the cooperativity between M···imidazole (M = Li(+), Na(+), K(+), Mg(2+), Cu(+), and Zn(2+)) and N(1)-H···OH(2) interactions. Topological atoms in molecules (AIM) analysis and charge analysis support the variation in hydrogen-bonding strength and the variation in M···imidazole binding strength. Effect of hydration on N(1)-H stretching frequency is studied, and it shows a clear shift in the stretching frequency after sequential hydration of metal ion as well as the N(1) of imidazole. The present study provides a detailed account on the biologically important M-histidine motif interaction with metal ions, where histidine is modeled by imidazole and methylimidazole.     

Effects of the biological backbone on stacking interactions at DNA-protein interfaces: the interplay between the backbone···π and π···π components

The (gas-phase) MP2/6-31G*(0.25) π···π stacking interactions between the five natural bases and the aromatic amino acids calculated using (truncated) monomers composed of conjugated rings and/or (extended) monomers containing the biological backbone (either the protein backbone or deoxyribose sugar) were previously compared. Although preliminary energetic results indicated that the protein backbone strengthens, while the deoxyribose sugar either strengthens or weakens, the interaction calculated using truncated models, the reasons for these effects were unknown. The present work explains these observations by dissecting the interaction energy of the extended complexes into individual backbone···π and π···π components. Our calculations reveal that the total interaction energy of the extended complex can be predicted as a sum of the backbone···π and π···π components, which indicates that the biological backbone does not significantly affect the ring system through π-polarization. Instead, we find that the backbone can indirectly affect the magnitude of the π···π contribution by changing the relative ring orientations in extended dimers compared with truncated dimers. Furthermore, the strengths of the individual backbone···π contributions are determined to be significant (up to 18 kJ mol(-1)). Therefore, the origin of the energetic change upon model extension is found to result from a balance between an additional (attractive) backbone···π component and differences in the strength of the π···π interaction. In addition, to understand the effects of the biological backbone on the stacking interactions at DNA-protein interfaces in nature, we analyzed the stacking interactions found in select DNA-protein crystal structures, and verified that an additive approach can be used to examine the strength of these interactions in biological complexes. Interestingly, although the presence of attractive backbone···π contacts is qualitatively confirmed using the quantum theory of atoms in molecules (QTAIM), QTAIM electron density analysis is unable to quantitatively predict the additive relationship of these interactions. Most importantly, this work reveals that both the backbone···π and π···π components must be carefully considered to accurately determine the overall stability of DNA-protein assemblies.     

Electron density characteristics in bond critical point (QTAIM) versus interaction energy components (SAPT): the case of charge-assisted hydrogen bonding

Charge-assisted hydrogen bonds (CAHBs) of N-H···Cl, N-H···Br, and P-H···Cl type were investigated using advanced computational approach (MP2/aug-cc-pVTZ level of theory). The properties of electron density function defined in the framework of Quantum Theory of Atoms in Molecules (QTAIM) were estimated as a function of distance in H-bridges. Additionally, the interaction energy decomposition was performed for H-bonded complexes with different H-bond lengths using the Symmetry-Adapted Perturbation Theory (SAPT). In this way both QTAIM parameters and SAPT energy components could be expressed as a function of the same variable, that is, the distance in H-bridge. A detailed analysis of the changes in QTAIM and SAPT parameters due to the changes in H···A distance revealed that, over some ranges of H···A distances, electrostatic, inductive and dispersive components of the SAPT interaction energy show a linear correlation with the value of the electron density at H-BCP ρ(BCP). The linear relation between the induction component, E(ind), and ρ(BCP) confirms numerically the intuitive expectation that the ρ(BCP) reflects directly the effects connected with the sharing of electron density between interacting centers. These conclusions are important in view of charge density studies performed for crystals in which the distance between atoms results not only from effects connected with the interaction between atomic centers directly involved in bonding, but also from packing effects which may strongly influence the length of the H-bond.     

Ionic liquid based on α-amino acid anion and N7,N9-dimethylguaninium cation ([dMG][AA]): theoretical study on the structure and electronic properties

The interactions between five amino acid based anions ([AA](-) (AA = Gly, Phe, His, Try, and Tyr)) and N7,N9-dimethylguaninium cation ([dMG](+)) have been investigated by the hybrid density functional theory method B3LYP together with the basis set 6-311++G(d,p). The calculated interaction energy was found to decrease in magnitude with increasing side-chain length in the amino acid anion. The interaction between the [dMG](+) cation and [AA](-) anion in the most stable configurations of ion pairs is a hydrogen bonding interaction. These hydrogen bonds (H bonds) were analyzed by the quantum theory of atoms in molecules (QTAIM) and natural bond orbital (NBO) analysis. Finally, several correlations between electron densities in bond critical points of hydrogen bonds and interaction energy as well as vibrational frequencies in the most stable configurations of ion pairs have been checked.     

CH stretching vibration of N-methylformamide as a sensitive probe of its complexation: infrared matrix isolation and computational study

The complexes between trans-N-methylformamide (t-NMF) and Ar, N(2), CO, H(2)O have been studied by infrared matrix isolation spectroscopy and/or ab initio calculations. The infrared spectra of NMF/Ne, NMF/Ar and NMF/N(2)(CO,H(2)O)/Ar matrices have been measured and the effect of the complexation on the perturbation of t-NMF frequencies was analyzed. The geometries of the complexes formed between t-NMF and Ar, N(2), CO and H(2)O were optimized in two steps at the MP2/6-311++G(2d,2p) level of theory. The four structures, found for every system at this level, were reoptimized on the CP-corrected potential energy surface; both normal and CP corrected harmonic frequencies and intensities were calculated. For every optimized structure the interaction energy was partitioned according to the SAPT scheme and the topological distribution of the charge density (AIM theory) was performed. The analysis of the experimental and theoretical results indicates that the t-NMF-N(2) and CO complexes present in the matrices are stabilized by very weak N-H···N and N-H···C hydrogen bonds in which the N-H group of t-NMF serves as a proton donor. In turn, the t-NMF-H(2)O complex present in the matrix is stabilized by O-H···O(C) hydrogen bonding in which the carbonyl group of t-NMF acts as a proton acceptor. Both, the theoretical and experimental results indicate that involvement of the NH group of t-NMF in formation of very weak hydrogen bonds with the N(2) or CO molecules leads to a clearly noticeable red shift of the CH stretching wavenumber whereas engagement of the CO group as a proton acceptor triggers a blue shift of this wavenumber.     

Interaction of coinage metal clusters with chalcogen dihydrides

The interaction of chalcogen dihydrides (H 2E; E = O, S, and Se) with small coinage metal clusters (M n ; M = Cu, Ag, and Au, n = 3 and 4) is studied based on density functional theory, with a focus on the nature of chalcogen-metal bonds. A newly developed pseudopotential-based correlation-consistent basis set is used for metal clusters together with the 6-311++G** basis set for the remaining atoms. Geometrical data identified that no significant deviation has been observed for molecules before and after complexation. For these three metals, binding energy calculations indicate that gold has the highest and silver has the lowest affinities for interaction with H 2E. In comparison with gold and copper, complexation between silver and chalcogen dihydrides is significantly weaker. It is found that interaction of H 2E molecules with the coinage metals have the order of H 2Se > H 2S > H 2O. Therefore, in agreement with experimental works, our calculations confirm that the gold-selenium bond is the most stable. The nature of M-E bonds is also interpreted by means of the quantum theory of atoms in molecules (QTAIM) and natural bond orbital (NBO) analyses. According to the QTAIM results, the bonds are found to be partially ionic and partially covalent. Natural resonance theory (NRT) is used to calculate natural bond order and bond polarity. The NRT result indicates that the percentage of polarity of M-E bonds is affected by coinage metals.     

Interactions of glutathione tripeptide with gold cluster: influence of intramolecular hydrogen bond on complexation behavior

Understanding the nature of the interaction between metal nanoparticles and biomolecules has been important in the development and design of sensors. In this paper, structural, electronic, and bonding properties of the neutral and anionic forms of glutathione tripeptide (GSH) complexes with a Au(3) cluster were studied using the DFT-B3LYP with 6-31+G**-LANL2DZ mixed basis set. Binding of glutathione with the gold cluster is governed by two different kinds of interactions: Au-X (X = N, O, and S) anchoring bond and Au···H-X nonconventional hydrogen bonding. The influence of the intramolecular hydrogen bonding of glutathione on the interaction of this peptide with the gold cluster has been investigated. To gain insight on the role of intramolecular hydrogen bonding on Au-GSH interaction, we compared interaction energies of Au-GSH complexes with those of cystein and glycine components. Our results demonstrated that, in spite of the ability of cystein to form highly stable metal-sulfide interaction, complexation behavior of glutathione is governed by its intramolecular backbone hydrogen bonding. The quantum theory of atom in molecule (QTAIM) and natural bond orbital analysis (NBO) have also been applied to interpret the nature of interactions in Au-GSH complexes. Finally, conformational flexibility of glutathione during complexation with the Au(3) cluster was investigated by means of monitoring Ramachandran angles.     

Analysis of M···H-Si interactions in [{M(CpSiMe2H)Cl3}2], (M = Zr, Hf, Ti and Mo) complexes

For the d(0) complex [{Zr(CpSiMe(2)H)Cl(3)}(2)] which contains a linear Si-H···Zr interaction across the dimer, DFT calculations are in good agreement with X-ray structures. The BP86 functional shows a slightly stronger interaction than B3LYP but for qualitative purposes either functional is sufficient. QTAIM analysis shows a bond critical point (bcp) for the interaction, a small negative value for the total energy density [H((r))] and the H atomic basin decreases in energy, E(H), and atomic volume compared to the free ligand. NBO analysis showed E(2) for Si-H σ to Zr(dz(2)) donation at 42.8 kcal mol(-1) and a 34% spatial overlap for the interaction consistent with an inverse hydrogen bond. The Wiberg bond index for the interaction is 0.1735 (0.7205 for the Si-H bond), ν((Si-H)) and (1)J((Si-H)) at 2060 cm(-1) and 145.4 Hz compared to 2183 cm(-1) and 172.1 Hz in the free ligand. Using a "synthesis by computation" approach to forming like complexes, similar features were found for [{Hf(CpSiMe(2)H)Cl(3)}(2)]. The titanium complex [{Ti(CpSiMe(2)H)Cl(3)}(2)] does not contain any Si-H···Ti interaction as rotation about the C-Si bond of the ligand occurs to place the Si-H bond hydrogen closer to a terminal chloro ligand across the dimer. An increase in electron density on the metal in the d(2) complex [{Mo(CpSiMe(2)H)Cl(3)}(2)] results in a stronger interaction with a distinct QTAIM analysis bcp [ρ((r)) 0.0448 a.u.], a small negative value for H((r)) and a much reduced H atomic volume. NBO analysis shows E(2) for Si-H σ to Mo(dz(2)) donation at 143.1 kcal mol(-1) and a 29% spatial overlap. Mo(dz(2)) to Si-H σ* donation (back donation) is minimal [E(2) 1.3 kcal mol(-1), ~1% spatial overlap]. The Wiberg bond index is 0.3114 (0.5667 for the Si-H bond), ν((Si-H)) 2015 cm(-1) and (1)J((Si-H)) 120.6 Hz.     

QTAIM study of the closed-shell interactions in peptide secondary structures: A cluster treatment of oligo- and polyalanines

The closed-shell interactions in oligo- and polyalanines are studied by the quantum theory of atoms in molecules (QTAIM) using electron densities derived from the B3LYP/6-31+G^** ground-state electronic wave-functions. The QTAIM enabled us to identify a large number of the intraturn closed-shell stabilizing interactions in the β-turns, which were presented by several conformers of the tetrapeptide model compound. We found that only β-turn type IVa exhibits a 10-member pseudocycle. The intrachain H-bonds between the adjacent N–H and CO groups in the antiparallel β-sheet conformation of polyalanine have not been found. At the same time, these interactions do exist in the parallel conformation and are even stronger than the interchain N–H…O bonds. A weak interaction between the CO group at the position i and the side-chain C–H group at the position i + 3 was detected in the -helical conformation of polyalanine.     

Actinium coordination chemistry: A density functional theory study with monodentate and bidentate ligands

In the current study, the coordination chemistry of nine-coordinate Ac(III) complexes with 35 monodentate and bidentate ligands was investigated using density functional theory (DFT) in terms of their geometries, charges, reaction energies, and bonding interactions. The energy decomposition analysis with naturals orbitals for chemical valence (EDA-NOCV) and the quantum theory of atoms in molecules (QTAIM) were employed as analysis methods. Trivalent Ac exhibits the highest affinities toward hard acids (such as charged oxophilic donors, fluoride), so its classification as a hard acid is justified. Natural population analysis quantified the involvement of 5f orbitals on Ac to be about 30% of total valence electron natural configuration indicating that Ac is a member of the actinide series. Pearson correlation coefficients were used to study the pairwise correlations among the bond lengths, ΔG reaction energies, charges on Ac and donor atoms, and data from EDA-NOCV and QTAIM. Strong correlations and anticorrelations were found between Voronoi charges on donor atoms with ΔG, EDA-NOCV interaction energies and QTAIM bond critical point densities.     

Toward a Consistent Interpretation of the QTAIM: Tortuous Link between Chemical Bonds,Interactions, and Bond/Line Paths

Currently, bonding analysis of molecules based on the Quantum Theory of Atoms in Molecules (QTAIM) is popular; however, “misinterpretations” of the QTAIM analysis are also very frequent. In this contribution the chemical relevance of the bond path as one of the key topological entities emerging from the QTAIM’s topological analysis of the one‐electron density is reconsidered. The role of nuclear vibrations on the topological analysis is investigated demonstrating that the bond paths are not indicators of chemical bonds. Also, it is argued that the detection of the bond paths is not necessary for the “interaction” to be present between two atoms in a molecule. The conceptual disentanglement of chemical bonds/interactions from the bonds paths, which are alternatively termed “line paths” in this contribution, dismisses many superficial inconsistencies. Such inconsistencies emerge from the presence/absence of the line paths in places of a molecule in which chemical intuition or alternative bonding analysis does not support the presence/absence of a chemical bond. Moreover, computational QTAIM studies have been performed on some “problematic” molecules, which were considered previously by other authors, and the role of nuclear vibrations on presence/absence of the line paths is studied demonstrating that a bonding pattern consistent with other theoretical schemes appears after a careful QTAIM analysis and a new “interpretation” of data is performed.     

Role of Basic Surface Groups of Activated Carbon in Chlordecone and β-Hexachlorocyclohexane Adsorption: A Molecular Modelling Study

The influence of nitrogen-containing surface groups (SGs) onto activated carbon (AC) over the adsorption of chlordecone (CLD) and β-hexachlorocyclohexane (β-HCH) was characterized by a molecular modelling study, considering pH (single protonated SGs) and hydration effect (up to three water molecules). The interactions of both pollutants with amines and pyridine as basic SGs of AC were studied, applying the multiple minima hypersurface (MMH) methodology and using PM7 semiempirical Hamiltonian. Representative structures from MMH were reoptimized using the M06-2X density functional theory. The quantum theory of atoms in molecules (QTAIM) was used to characterize the interaction types in order understanding the adsorption process. A favorable association of both pesticides with the amines and pyridine SGs onto AC was observed at all pH ranges, both in the absence and presence of water molecules. However, a greater association of both pollutants with the primary amine was found under an acidic pH condition. QTAIM results show that the interactions of CLD and β-HCH with the SGs onto AC are governed by Cl···C interactions of chlorine atoms of both pesticides with the graphitic surface. Electrostatic interactions (H-bonds) were observed when water molecules were added to the systems. A physisorption mechanism is suggested for CLD and β-HCH adsorption on nitrogen-containing SGs of AC.     

Concerted halogen and hydrogen bonding in [RuI2(H2dcbpy)(CO)2]···I2···(CH3OH)···I2···[RuI2(H2dcbpy)(CO)2

A new type of concerted halogen bond-hydrogen bond interaction was found in the solid state structure of [RuI(2)(H(2)dcbpy)(CO)(2)]···I(2)···(MeOH)···I(2)···[RuI(2)(H(2)dcbpy)(CO)(2)]. The iodine atoms of the two I(2) molecules interact simultaneously with each other and with the OH group of methanol of crystallization. The interaction was characterized by single crystal X-ray measurements and by computational charge density analysis based on DFT calculations.     

Tuning the activity of glutathione peroxidase mimics through intramolecular Se···N,O interactions: a DFT study incorporating solvent-assisted proton exchange (SAPE)

Diaryl diselenide mimics of the antioxidant selenoprotein glutathione peroxidase (GPx) often incorporate intramolecular Se···N,O interactions to enhance their GPx-like activity. Although the strength of the interaction is defined by the Lewis basicity of the donating group and the strength of the Se-X bond, there is not a clear relationship between the interaction and the GPx-like activity. Density-functional theory and natural bond orbital (NBO) calculations are used to show the range of Se···N,O interactions for various functional groups. The strongest interactions are found for groups which stabilize the donor-acceptor interaction through aromatic stabilization. The activation barriers for the GPx-like mechanism of activity of several substituted areneselenols are calculated using DFT and solvent-assisted proton exchange (SAPE), a technique that incorporates networks of solvent molecules into the theoretical model to facilitate proton transfer between sites in the reactant and product. DFT-SAPE models show that, in addition to decreasing the barrier to oxidation of the selenol, Se···N,O interactions generally increase the barriers for selenenic acid reduction and selenol regeneration because the Se···N,O interaction must be broken for the reaction to proceed. Calculated activation barriers for the rate-determining step are consistent with the relative experimental GPx-like activities of a series of diaryl diselenides.     

Characterizing Pressure‐Induced Uranium CH Agostic Bonds

The diuranium(III) compound [UN′′₂]₂(μ‐η⁶:η⁶‐C₆H₆) (N′′=N(SiMe₃)₂) has been studied using variable, high‐pressure single‐crystal X‐ray crystallography, and density functional theory . In this compound, the low‐coordinate metal cations are coupled through π‐ and δ‐symmetric arene overlap and show close metal? CH contacts with the flexible methyl CH groups of the sterically encumbered amido ligands. The metal–metal separation decreases with increasing pressure, but the most significant structural changes are to the close contacts between ligand CH bonds and the U centers. Although the interatomic distances are suggestive of agostic‐type interactions between the U and ligand peripheral CH groups, QTAIM (quantum theory of atoms‐in‐molecules) computational analysis suggests that there is no such interaction at ambient pressure. However, QTAIM and NBO analyses indicate that the interaction becomes agostic at 3.2 GPa.     

Theoretical investigation on charge‐assisted halogen bonding interactions in the complexes of bromocarbons with some anions

A series of dimeric complexes formed between bromocarbon molecules and two anions (Br^? and CN^?) have been investigated by using MP2 method. The quantum theory of atoms in molecules (QTAIM) and the second‐order perturbation natural bond orbital (NBO) approaches were applied to analyze the electron density distributions of these complexes and to explore the nature of charge‐assisted halogen bonding interactions. As anticipated, these interactions are significantly stronger relative to the corresponding neutral ones. The results derived from ab initio calculations described herein reveal a major contribution from the electrostatic interaction on the stability of the systems considered. Beside the electrostatic interaction, the charge‐transfer force and the second‐order orbital interaction also play an important role in the formation of the complexes, as a NBO analysis suggested. The presence of halogen bonds in the complexes has been identified in terms of the QTAIM methodology, and several linear relationships have been established to provide more insight into charge‐assisted halogen bonding interactions. © 2007 Wiley Periodicals, Inc. Int J Quantum Chem, 2008     

QTAIM and ETS-NOCV analyses of intramolecular CH···HC interactions in metal complexes

The topological analysis, based on the quantum theory of atoms in molecules (QTAIM) of Bader and the ETS-NOCV charge and energy decomposition method have been used to characterize coordination bonds, chelating rings, and additional intramolecular interactions in the ZnNTA and ZnNTPA complexes in solvent. The QTAIM and ETS-NOCV studies have conclusively demonstrated that the H-clashes (they are observed only in the ZnNTPA complex and classically are interpreted as steric hindrance destabilizing a complex) are characterized by (i) the electron flow channel between the H-atoms involved, as discovered by the ETS-NOCV analysis (on average, ΔE(orb) = -1.35 kcal mol(-1)) and (ii) QTAIM-defined a bond path that indicates the presence of a preferred quantum-mechanical exchange channel, hence, they should be seen as H-H intramolecular bonding interactions. The main reason for the formation of a weaker ZnNTPA complex was attributed to the strain energy (from both QTAIM and ETS-NOCV techniques) and the larger Pauli repulsion contribution found from the ETS-NOCV analysis. An excellent agreement between physical properties controlling the stability of the two complexes was found from the two techniques, QTAIM and ETS-NOCV.