Microsolvation of aminoethanol: a study using DFT combined with QTAIM

The microsolvation of aminoethanol (AE) with one, two, three or four water molecules was investigated using a density functional theory (DFT) approach. Quantum theory of atoms in molecules (QTAIM) analyses were employed to elucidate the hydrogen-bonding characteristics of AE–(H₂O) _n (n = 1–4) complexes. The results showed that AE tends to break its intramolecular OH^AE···N^AE hydrogen bond (H-bond) upon microsolvation and form intermolecular H-bonds with water molecules, while complexes that retain the intramolecular OH^AE···N^AE H-bond show reduced stabilities. The intermolecular H-bond that forms between the nitrogen atom of AE and the hydroxyl of a water molecule is the strongest one for the most stable AE–(H₂O) _n (n = 1–4) complexes, and as n increases from 1 to 4 they grow stronger. The partial covalent character of this H-bond was confirmed by QTAIM analyses. Many-body interaction analysis showed that the relaxation energies and two- and three-body energies make significant contributions to the binding energies of the complexes.     

Hydrogen-bonding interactions in adrenaline–water complexes: DFT and QTAIM studies of structures, properties, and topologies

ωB97XD/6-311++G(d,p) calculations were carried out to investigate the hydrogen-bonding interactions between adrenaline (Ad) and water. Six Ad-H(2)O complexes possessing various types of hydrogen bonds (H-bonds) were characterized in terms of their geometries, energies, vibrational frequencies, and electron-density topology. Natural bond orbital (NBO) and quantum theory of atoms in molecules (QTAIM) analyses were performed to elucidate the nature of the hydrogen-bonding interactions in these complexes. The intramolecular H-bond between the amino and carboxyl oxygen atom of Ad was retained in most of the complexes, and cooperativity between the intra- and intermolecular H-bonds was present in some of the complexes. H-bonds in which hydroxyls of Ad/water acted as proton donors were stronger than other H-bonds. Both hydrogen-bonding interactions and structural deformation play important roles in the relative stabilities of the complexes. The intramolecular H-bond was broken during the formation of the most stable complex, which indicates that Ad tends to break the intramolecular H-bond and form two new intermolecular H-bonds with the first water molecule.     

Microsolvation effect and hydrogen-bonding pattern of taurine-water TA-(H2O)n (n = 1-3) complexes

The microsolvation of taurine (TA) with one, two or three water molecules was investigated by a density functional theory (DFT) approach. Quantum theory of atoms in molecules (QTAIM) analyses were employed to elucidate the hydrogen bond (H-bond) interaction characteristics in TA-(H₂O)_n (n = 1–3) complexes. The results showed that the intramolecular H-bond formed between the hydroxyl and the N atom of TA are retained in most TA-(H₂O)_n (n = 1–3) complexes, and are strengthened via cooperative effects among multiple H-bonds from n = 1–3. A trend of proton transformation exists from the hydroxyl to the N atom, which finally results in the cleavage of the origin intramolecular H-bond and the formation of a new intramolecular H-bond between the amino and the O atom of TA. Therefore, the most stable TA-(H₂O)₃ complex becomes a zwitterionic complex rather than a neutral type. A many-body interaction analysis showed that the major contributors to the binding energies for complexes are the two-body energies, while three-body energies and relaxation energies make significant contributions to the binding energies for some complexes, whereas the four-body energies are too small to be significant.     

Quantum chemical topological analysis of hydrogen bonding in HX…HX and CH3X…HX dimers (X = Br,Cl, F)

We present a systematic investigation of the nature and strength of the hydrogen bonding in HX···HX and CH₃X…HX (X = Br, Cl and F) dimers using ab initio MP2/aug-cc-pVTZ calculations in the framework of the quantum theory of atoms in molecules (QTAIM) and electron localisation functions (ELFs) methods. The electron density of the complexes has been characterised, and the hydrogen bonding energy, as well as the QTAIM and ELF parameters, is consistent, providing deep insight into the origin of the hydrogen bonding in these complexes. It was found that in both linear and angular HX…HX and CH₃X…HX dimers, F atoms form stronger HB than Br and Cl, but they need short (～2 Å) X…HX contacts.     

Theoretical investigations of the H···π and X (X = F,Cl, Br,I)···π complexes between hypohalous acids and benzene

The H···π and X (X = F, Cl, Br, I)···π interactions between hypohalous acids and benzene are investigated at the MP2/6-311++G(2d,2p) level. Four hydrogen-bonded and three halogen-bonded complexes were obtained. Ab initio calculations indicate that the X···π interaction between HOX and C₆H₆ is mainly electrostatically driven, and there is nearly an equal contribution from both electrostatic and dispersive energies in the case of XOH–C₆H₆ complexes. Natural bond orbital (NBO) analysis reveals that there exists charge transfer from benzene to hypohalous acids. Atom in molecules (AIM) analysis locates bond critical points (BCP) linking the hydrogen or halogen atom and carbon atom in benzene.     

Density functional theory study of the potassium complexation of an unsymmetrical 1,3-alternate calix[4]-crown-5-N-azacrown-5 bearing two different crown rings

Theoretical studies of an unsymmetrical calix[4]-crown-5-N-azacrown-5 (1) in a fixed 1,3-alternate conformation and the complexes 1·K⁺(a), 1·K⁺(b), 1·K⁺(c) and 1·K⁺K⁺ were performed using density functional theory (DFT) at the B3LYP/6-31G* level. The fully optimized geometric structures of the free macroligand and its 1:1 and 1:2 complexes, as obtained from DFT calculations, were used to perform natural bond orbital (NBO) analysis. The two main types of driving force metal–ligand and cation–π interactions were investigated. NBO analysis indicated that the stabilization interaction energies (E ₂) for O…K⁺ and N…K⁺ are larger than the other intermolecular interactions in each complex. The significant increase in electron density in the RY* or LP* orbitals of K⁺ results in strong host–guest interactions. In addition, the intermolecular interaction thermal energies (ΔE, ΔH, ΔG) were calculated by frequency analysis at the B3LYP/6-31G* level. For all structures, the most pronounced changes in the geometric parameters upon interaction are observed in the calix[4]arene molecule. The results indicate that both the intermolecular electrostatic interactions and the cation–π interactions between the metal ion and π orbitals of the two pairs that face the inverted benzene rings play a significant role.     

Insight into the nature of the interactions of furan and thiophene with hydrogen halides and lithium halides: ab initio and QTAIM studies

The nature of the interactions of furan and thiophene with hydrogen halides and lithium halides has been investigated using ab initio calculations and QTAIM analysis. The concept of molecule formation density difference (MFDD) is introduced to study weak hydrogen bond (HB) and lithium bond (LB) interactions. The results have shown the molecular electrostatic potentials of furan and thiophene, as well as of the hydrogen halides and lithium halides, determine the geometries of the complexes. Both the studied HB and LB interactions can be classified as “closed-shell” weak interactions. The topological properties and energy properties at the bond critical points of HB and LB have been shown to be exponentially dependent on intermolecular distances d(H-bond) and d(Li-bond), which enables interpretation of the strength of the HB and LB interactions in terms of these ρ(r) properties. Electron transfer plays a more important role in the formation of HB than in that of LB, while electrostatic interaction in LB is more dominant than that in HB.     

Hydrogen bonding interactions in PN···HX complexes: DFT and ab initio studies of structure, properties and topology

Spin-restricted DFT (X3LYP and B3LYP) and ab initio (MP2(fc) and CCSD(fc)) calculations in conjunction with the Aug-CC-pVDZ and Aug-CC-pVTZ basis sets were performed on a series of hydrogen bonded complexes PN···HX (X = F, Cl, Br) to examine the variations of their equilibrium gas phase structures, energetic stabilities, electronic properties, and vibrational characteristics in their electronic ground states. In all cases the complexes were predicted to be stable with respect to the constituent monomers. The interaction energy (ΔE) calculated using a super-molecular model is found to be in this order: PN···HF > PN···HCl > PN···HBr in the series examined. Analysis of various physically meaningful contributions arising from the Kitaura-Morokuma (KM) and reduced variational space self-consistent-field (RVS-SCF) energy decomposition procedures shows that the electrostatic energy has significant contribution to the over-all interaction energy. Dipole moment enhancement (Δμ) was observed in these complexes expected of predominant dipole-dipole electrostatic interaction and was found to follow the trend PN···HF > PN···HCl > PN···HBr at the CCSD level. However, the DFT (X3LYP and B3LYP) and MP2 levels less accurately determined these values (in this order HF < HCl < HBr). Examination of the harmonic vibrational modes reveals that the PN and HX bands exhibit characteristic blue- and red shifts with concomitant bond contraction and elongation, respectively, on hydrogen bond formation. The topological or critical point (CP) analysis using the static quantum theory of atoms in molecules (QTAIM) of Bader was considered to classify and to gain further insight into the nature of interaction existing in the monomers PN and HX, and between them on H-bond formation. It is found from the analysis of the electron density ρ _c , the Laplacian of electron charge density ∇²ρ_c, and the total energy density (H _c ) at the critical points between the interatomic regions that the interaction N···H is indeed electrostatic in origin (ρ_c > 0, ∇²ρ_c > 0 and H_c > 0 at the BCP) whilst the bonds in PN (ρ_c > 0, ∇²ρ_c > 0 and H_c < 0) and HX ((ρ_c > 0, ∇²ρ_c < 0 and H_c < 0)) are predominantly covalent. A natural bond orbital (NBO) analysis of the second order perturbation energy lowering, E⁽²⁾, caused by charge transfer mechanism shows that the interaction N···H is n(N) → BD*(HX) delocalization.     

Exploring the nature of the H-bonds between the human class II MHC protein,HLA-DR1 (DRB*0101) and the influenza virus hemagglutinin peptide,HA306-318, using the quantum theory of atoms in molecules

The nature of the H-bonds between the human protein HLA-DR1 (DRB*0101) and the hemagglutinin peptide HA306-318 has been studied using the Quantum Theory of Atoms in Molecules for the first time. We have found four H-bond groups: one conventional CO··HN bond group and three nonconventional CO··HC, π··HC involving aromatic rings and HN··HC_aliphatic groups. The calculated electron density at the determined H-bond critical points suggests the follow protein pocket binding trend: P1 (2,311) >> P9 (1.109) > P4 (0.950) > P6 (0.553) > P7 (0.213) which agrees and reveal the nature of experimental findings, showing that P1 produces by a long way the strongest binding of the HLA-DR1 human protein molecule with the peptide backbone as consequence of the vast number of H-bonds in the P1 area and at the same time the largest specific binding of the peptide Tyr308 residue with aromatic residues located at the binding groove floor. The present results suggest the topological analysis of the electronic density as a valuable tool that allows a non-arbitrary partition of the pockets binding energy via the calculated electron density at the determined critical points.     

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.

     

Aromaticity balance, π-electron cooperativity and H-bonding properties in tautomerism of salicylideneaniline: the quantum theory of atoms in molecules (QTAIM) approach

Topological analysis based on DFT calculations regarding proton transfer reaction in salicylideneaniline (SA) was performed to scrutinize possible changes in the intramolecular H-bond, π-electron delocalization and aromaticity levels of certain fragments. Quantum chemical calculations and natural bond orbital (NBO) analyses were carried out over a tautomeric ensemble whose members correspond to the molecules at different stages in tautomeric interconversion of SA. The elaboration of intramolecular hydrogen bonding in terms of the relevant topological parameters and the interpretation of certain dependencies regarding its strength were examined. The results show that delocalization index (DI) between donor and acceptor atom δ(O,N) is a useful topological parameter for describing H-bond strength, which is influenced by π-delocalization level within quasiaromatic chelate ring, indicating its resonance-assisted character. NBO analyses reveal that lone-pair (LP) population on N center also affects the strength of intramolecular H-bond in SA. Furthermore, π-electron transfer accompanying intramolecular proton migration in SA is brought into being through formally vacant non-Lewis type LP* orbital on the tautomeric proton. As a result of this, tautomeric protons in molecular entities near TS have hypovalent character due to the lack of electron population in the bonding orbital relative to that in LP* orbital. While H-bonds in the tautomeric ensemble of SA are predominantly partial covalent, molecular entities close to transition state have the strongest covalent H-bonds. The most important result is also that there are linear correlations between the orders of bonds (hydroxyl and amine) involving intramolecular H-bond and electron density values at the relevant BCPs due to partially covalent character of these bonds, contrary to exponential behavior as for purely covalent bonds. Quasiaromatic chelate ring formation is established not only to compel a reduced aromaticity of salicylidene ring but also to decrease in LP-population on N.     

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.     

Non-covalent interactions – QTAIM and NBO analysis

MP2(full)/6-311++G(3df,3pd) calculations were carried out on complexes linked through various non-covalent Lewis acid – Lewis base interactions. These are: hydrogen bond, dihydrogen bond, hydride bond and halogen bond. The quantum theory of ´atoms in molecules´ (QTAIM) as well as the natural bond orbitals (NBO) method were applied to analyze properties of these interactions. It was found that for the A-H…B hydrogen bond as well as for the A-X…B halogen bond (X designates halogen) the complex formation leads to the increase of s-character in the A-atom hybrid orbital aimed toward the H or X atom. In opposite, for the A…H-B hydride bond, where the H-atom possesses negative charge, the decrease of s-character in the B-atom orbital is observed. All these changes connected with the redistribution of the electron charge being the effect of the complex formation are in line with Bent´s rule. The numerous correlations between energetic, geometrical, NBO and QTAIM parameters were also found.

Quantum-Chemical Analysis of C-H…O and C-H…N Interactions in RNA Base Pairs—H-Bond Versus Anti-H-Bond Pattern

Abstract

Geometries and interaction energies of unusual UU and AA base pairs with one standard hydrogen bond (H-bond) and additional C-H…O or C-H…N contacts have been determined by quantum-chemical methods taking into account electron correlation. Whereas the C-H bond length in the UU C-H…O contact increases upon complex formation (H-bond pattern), the C-H bond of the AA C-H….N interaction is shortened (anti-H-bond pattern). The same properties are found for model complexes between U or A and formaldehyde that have intermolecular C-H…acceptor contacts but no standard H-bonds. Both the C-H…acceptor H-bond and anti-H-bond interactions are attractive. A possible influence of the donor CH group charge distribution on the interaction pattern is discussed.     

Density functional theory study on the interaction between keto-9H guanine and aspartic acid

A theoretical study was performed using density functional theory (DFT) to investigate hydrogen bonding interactions in signature complexes formed between keto-9H guanine (Gua) and aspartic acid (Asp) at neutral pH. Optimized geometries, binding energies and the theoretical IR spectra of guanine, aspartic acid and their corresponding complexes (Gua-Asp) were calculated using the B3LYP method and the 6-31+G(d) basis set. Stationary points found to be at local minima on the potential energy surface were verified by second derivative harmonic vibrational frequency calculations at the same level of theory. AIM theory was used to analyze the hydrogen bonding characteristics of these DNA base complex systems. Our results show that the binding motif for the most stable complex is strikingly similar to a Watson-Crick motif observed in the guanine-cytosine base pair. We have found a range of hydrogen bonding interactions between guanine and aspartic acid in the six complexes. This was further verified by theoretical IR spectra of ω(C-H---O-H) cm⁻¹ stretches for the Gua-Asp complexes. The electron density plot indicates strong hydrogen bonding as shown by the 2p _z dominant HOMO orbital character.     

Hydrogen bonds in Zif268 proteins – a theoretical perspective

The aim of the work was to elucidate the presence of different hydrogen bond (H-bond) in five Zif268 proteins (1A1F, 1A1G, 1A1H, 1A1I and 1A1K). For this purpose, we have performed the QM/MM and molecular dynamics (MD) studies, the results of which reveal that H-bonds depend on the amino acid sequence and orientation of the H-bond donor atoms. Further, high specificity of Arg and Asn is observed for guanine and adenine, respectively. Furthermore, both conventional and non-conventional hydrogen bond also exists in the proteins, among them N–H?O H-bonds are the strongest. Besides, the non-conventional bonds play a role in the protein folding and DNA stacking. From the QSAR properties, amino acids such as asparagine and aspartic acids are the major reactive sites in the Zif268 protein. The electron affinities of Zif268 proteins are high, so the charge transfer occurs from the DNA to the protein molecules. NBO analysis indicates the majority of charge transfer occurs from DNA to the corresponding anti-bonding orbital of the peptides. Root mean square deviation and R_g (radius of gyration) show that 1A1F is more compact and in native state during MD simulation. The minimum R_g leads to the large number of hydrogen bonds formation in 1A1F. Higher solvent accessible surface area in 1A1I indicates that the cavity inside the protein is large.     

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.

Abundance and arrangement of single,double and bifurcated hydrogen bonds in protein α-helices: Statistical analysis of PDB-select

Statistics are collected and analyzed for the possibility of hydrogen bonding in the secondary structures of globular proteins, based on geometric criteria. Double and bifurcated bonds are considered as pairs of admissible H-bonds with two proton donors or two proton acceptors, respectively. Most of such bonds belong to peptide groups in α-helices, with O _i …N _{i + 3} nearly as frequent as O _i …N _{i + 4}; in contrast, most of the 3/10-helical segments are too short to have any. Alternating double and bifurcated bonds in α-helices form an apparently cooperative network structure. A typical α-helical segment perhaps carries two stretches of the H-bond network broken in the middle. The constituent H-bonds are nonlinear: the hydrogen atom is off the straight line connecting the proton donor and proton acceptor atoms. This deflection is larger for H _{i + 3} vs. bond line O _i −N _{i + 3} than for H _{i + 4} vs. O _i −N _{i + 4}, and though the two kinds of bond have about the same length (exceeding those typical of low-molecular compounds), O _i …N _{i + 4} must be stronger than O _i …N _{i + 3}. Double/bifurcated bonds are also not coplanar, i.e., hydrogen atoms are beyond the N…O…N (or O…N…O) plane. The text was submitted by the authors in English.     

CH...O hydrogen bonds at protein-protein interfaces

For the first time, a statistical potential has been developed to quantitatively describe the CH.O hydrogen bonding interaction at the protein-protein interface. The calculated energies of the CH.O pair interaction show a favorable valley at approximately 3.3 A, exhibiting a feature typical of an H-bond and similar to the ab initio quantum calculation result (Scheiner, S., Kar, T., and Gu, Y. (2001) J. Biol. Chem. 276, 9832-9837). The potentials have been applied to a set of 469 protein-protein complexes to calculate the contribution of different types of interactions to each protein complex: the average energy contribution of a conventional H-bond is approximately 30%; that of a CH.O H-bond is 17%; and that of a hydrophobic interaction is 50%. In some protein-protein complexes, the contribution of the CH.O H-bond can reach as high as approximately 40-50%, indicating the importance of the CH.O H-bond at the protein interface. At the interfaces of these complexes, C(alpha)H.O H-bonds frequently occur between adjacent strands in both parallel and antiparallel orientations, having the obvious structural motif of bifurcated H-bonds. Our study suggests that the weak CH.O H-bond makes an important contribution to the association and stability of protein complexes and needs more attention in protein-protein interaction studies.