Hydrogen-Bonding Interaction in a Complex of Amino Acid with <Emphasis Type="Italic">N</Emphasis>,<Emphasis Type="Italic">N</Emphasis>-Dimethylformamide Studied by DFT Calculations

Theoretical studies on hydrogen-bonded complexes between amino acids (glycine, alanine and leucine) and N,N-dimethylformamide (DMF) in gas phase have been carried out using density functional theory (DFT) and ab initio calculations at the B3LYP/6-311++G** and MP2/6-311++G** theory levels. The structures, binding energy, stretching frequency and bond characteristics of the mentioned complexes were calculated. The NH₂ and COOH groups of amino acids form different types of hydrogen bonds with the DMF molecule, as well as alkyl side chains. High binding energy suggests multiple hydrogen bonds present in one complex. The nearly linear OH???O and NH???O contacts are stronger than a conventional hydrogen bond interaction with their H???O separation between 1.74 and 2.14 Å. The weaker CH???O H-bond is also discussed as being a crucial interaction in biological systems involving amino acids. The formation of this interaction results in a blue shift in the CH stretching frequency.     

Rapid evaluation of the binding energies between peptide amide and DNA base

The binding energies and the equilibrium hydrogen bond distances as well as the potential energy curves of 20 hydrogen‐bonded amide–base dimers are evaluated from the analytic potential energy function established in our laboratory recently. The analytic potential energy function is used to calculate the N? H···N, N? H···O?C, C? H···N, and C? H···O?C dipole–dipole attractive interaction energies and C?O···O?C, N? H···H? N, and N? H···H? C dipole–dipole repulsive interaction energies in the 20 dimers composed of DNA bases adenine, guanine, cytosine, or thymine and peptide amide. The calculation results show that the potential energy curves obtained from the analytic potential energy function are in good agreement with those obtained from MP2/6‐311+G** calculations by including the basis set superposition error (BSSE) correction. For all the 20 dimers, the analytic potential energy function yields the binding energies of the MP2/6‐311+G** with BSSE correction within the error limits of 0.50 kcal/mol for 19 dimers, only one difference is larger than 0.50 kcal/mol and the difference is only 0.61 kcal/mol. The analytic potential energy function produces the equilibrium hydrogen bond distances of the MP2/6‐311+G** with BSSE correction within the error limits of 0.030 Å for all the 20 dimers. The analytic potential energy function is further applied to four more complicated DNA base‐peptide amide systems involving amino acid side chain and β‐sheet. The values of the binding energies and equilibrium hydrogen bond distances obtained from the analytic potential energy function are also in good agreement with those obtained from MP2 calculations with the BSSE correction. These results demonstrate that the analytic potential energy function can be used to evaluate the binding energies in hydrogen‐bonded peptide amide–DNA base dimers quickly and accurately. © 2011 Wiley Periodicals, Inc. J Comput Chem, 2011     

4-羟甲基吡啶－水红移氢键的密度泛函理论研究

使用密度泛函理论B3LYP方法和6-311++G**基函数对4-羟甲基吡啶与水形成的1:1和1:2（摩尔比）氢键复合物进行了理论计算研究,分别得到稳定的4-羟甲基吡啶－H₂O和4-羟甲基吡啶－(H₂O)₂氢键复合物3个和8个。经基组重叠误差和零点振动能校正后,最稳定的1:1和1:2氢键复合物的相互作用能分别为-20.536和-44.246 kJ/mol。振动分析显示O-H···N(O)氢键的形成使复合物中O-H键对称伸缩振动频率红移(减小)。自然键轨道分析表明,4-羟甲基吡啶与水形成最稳定的1:1和1:2氢键复合物时,分子间电荷转移分别为0.02642 e 和0.03813 e 。含时密度泛函理论TD-B3LYP/ 6-311++G**计算显示,相对于4-羟甲基吡啶单体分子,氢键H-OH···N和H-OH···OH的形成分别使最大吸收光谱波长兰移8~16纳米和红移4~11纳米。     

Density Functional Theory Study of Hydrogen Bonding Dimers with 4-Pyridinecarboxylic Acid Hydrazine

The H‐bonding dimers of 4‐pyridinecarboxylic acid hydrazine were studied using density functional theory (DFT) at B3LYP/6‐311++G** level. The results showed that the most stable dimer D1 had two same linear N H···O hydrogen bonds, and the interaction energy between them was 51.038 kJ·mol⁻¹ which was corrected by the basis set superposition error and zero‐point. The stretching vibration frequency of N H bond had a red shift because of the hydrogen bonds. The natural bond orbital analysis showed that each N H···O hydrogen bond in D1 had the biggest interaction stabilization energy of 69.078 kJ·mol⁻¹. Thermodynamic analysis indicated that the formation process of D1 was exothermic and spontaneous at low and room temperatures.     

Intramolecular hydrogen bonding in 3‐imino‐propenylamine: Theoretical investigations

Intramolecular H‐bonds existing for derivatives of 3‐imino‐propenylamine have been studied using the B3LYP/6‐311++G** level of theory. The nature of these interactions, known as resonance‐assisted hydrogen bonds, has been discussed. Vibrational frequencies for α‐derivatives were calculated at the same level of theory. The topological properties of the electron density distributions for N? H···N intramolecular bridges have been analyzed in terms of the Bader theory of atoms in molecules (AIM). Calculation for 3‐imino‐propenylamine derivatives in water solution were also carried out at B3LYP/6‐311++G** level of theory. Finally, the analysis of hydrogen bond in this molecule and their derivatives by quantum theory of natural bond orbital methods fairly support the ab initio results. © 2009 Wiley Periodicals, Inc. Int J Quantum Chem, 2009     

Intramolecular hydrogen bond in 3‐imino‐propenylamine isomers: AIM and NBO studies

The molecular structure and intramolecular hydrogen bond energy of 18 conformers of 3‐imino‐propenyl‐amine were investigated at MP2 and B3LYP levels of theory using the standard 6‐311++G** basis set. The atom in molecules or AIM theory of Bader, which is based on the topological properties of the electron density (ρ), was used additionally and the natural bond orbital (NBO) analysis was also carried out. Furthermore calculations for all possible conformations of 3‐imino‐propenyl‐amin in water solution were also carried out at B3LYP/6‐311++G** and MP2/6‐311++G** levels of theory. The calculated geometrical parameters and conformational analyses in gas phase and water solution show that the imine–amine conformers of this compound are more stable than the other conformers. B3LYP method predicts the IMA‐1 as global minimum. This stability is mainly due to the formation of a strong N? H···N intramolecular hydrogen bond, which is assisted by π‐electrons resonance, and this π‐electrons are established by NH2 functional group. Hydrogen bond energies for all conformers of 3‐imino‐propenyl‐amine were obtained from the related rotamers methods. © 2009 Wiley Periodicals, Inc. Int J Quantum Chem, 2010     

Characterization of the N?H···ON and N?H···NO H‐bonds in nitrosamine dimers

Ab initio molecular orbital and DFT calculations have been carried out for three most stable dimers of parent nitrosamine (NA) in order to elucidate the structures and energetics of the dimers. The structures were optimized using HF, B3LYP, and MP2 methods with 6‐311+G(d,p) and 6‐311++G(2d,2p) basis sets. At the optimized geometries obtained at MP2/6‐311++G(2d,2p) level of theory, the energies were evaluated at QCISD/aug‐cc‐pVDZ and CCSD/aug‐cc‐pVDZ levels. The most stable dimer has two N? H···O?N hydrogen bonds and the least stable dimer has two N? H···N?O hydrogen bonds. The natural bond orbital analysis showed that the lpO(N) → BD*(N? N) and lpO(N) → BD*(N? H_b) interactions play a decisive role in the stabilization of the NH···O(N) hydrogen bonds in dimers. The atoms in molecules results reveal that the intermolecular N? H···O(N) H‐bonds in dimers have electrostatic character. © 2007 Wiley Periodicals, Inc. Int J Quantum Chem, 2008     

Basicity of the nitrogen atom and of the carbonyl function of methyl isonicotinate and 4-acetylpyridine. An infrared study of the interaction with proton donors

The hydrogen bonded complexes between methyl isonicotinate and 4-acetylpyridine and phenol derivatives acting as proton donors have been investigated by ir spectroscopy. The thermodynamic parameters (K, — ΔH°, — ΔS°) have been determined. The ir data show that the hydrogen bond interaction occurs at the carbonyl function and at the nitrogen atom of the ring. When the acidity of the proton donor increases, N-complexation is favoured over carbonyl complexation and with the stronger acid hydrochloric acid, only N-protonated species are observed. The data are compared with those obtained for closely related pyridine derivatives, bearing a X-C=O substituent and it is shown that the proportion of OH···N and OH···O=C species is related to the basicity of the nitrogen atom of the heterocyclic ring and to inductive and mesomeric effects depending on the nature of X.     

Density functional theory and topological analysis on the hydrogen bonds in cysteine–propanoic acid complexes

The complexes formed via hydrogen bonding interactions between cysteine and propanoic acid have been studied at the density three-parameter hybrid functional DFT-B3LYP/6-311++G(d,p) level regarding their geometries, energies, vibrational frequencies, and topological features of the electron density. The quantum theory of atoms in molecules (QTAIM) and natural bond orbital (NBO) analysis was employed to elucidate the interaction characteristics in cysteine–propanoic acid (Cys–Prop) complexes. More than 10 kinds of hydrogen bonds (H-bonds) including intra- and inter-molecular H-bonds have been found in Cys–Prop complexes. The results show that both the strength of H-bonds and the deformation are important factors for the stability of Cys–Prop complexes. The strongest H-bonds (O2H^A···O1^B and O2H^A···O1^B) exist in the most stable Cys–Prop complex. The stronger H-bonds formed between hydroxyl and O (or N) atom usually stronger than those involve C (or S) atom. Relationships between the electron density (ρ) of BCP and H-bond length as well as the Fock matrix element (F _ij) has also been investigated and used to study the nature of H-bonds. Moreover, the results show that the change of the bond length linearly correlates with the corresponding frequency shift.     

Hydrogen Bonding Interaction of Formic Acid-, Formaldehyde-, Formylfluoride-Nitrosyl Hydride： Theoretical Study on the Geometries, Interaction Energies and Blue- or Red-Shifted Hydrogen Bonds

The hydrogen bonding interaction of formic acid-, formaldehyde-, formylfluoride-nitrosyl hydride complexes was investigated by the density functional theory （DFT） and ab inito method in conjunction with 6-311＋＋G（2d,2p） basis set. The geometries, vibrational frequencies and interaction energies of the complexes were calculated by both standard and CP-corrected methods respectively. Moreover, G3B3 method was employed to estimate the interaction energies. There are C--H…O, N--H…O, N--H…F blue-shifted H-bonds and red-shifted O----H…O H-bond in the complexes. Electron density redistribution and rehybridization contribute to the N--H and C--H blue shifts. All geometric reorganizations contribute to the N--H blue shifts and partial geometric reorganizations contribute to the C--H blue shifts. The geometric reorganizations of the complex C except ZH（5）-O（4）-C（1） contribute to the O----H red shift. For the N--H blue shifts, the effect of r（N--O） variation on the N--H blue shifts is larger than that of ZH-N-O variation. Rehybridization plays a dominant role in the degree of N--H blue shifts, whereas the electron density redistribution contributes more to the degree of C--H blue shifts than the other effects do.     

Theoretical study of molecular interactions of sulfur ylide with HF, HCN, and HN3

The molecular interactions between sulfur ylide (SY) and HX molecules (X = F, CN, and N₃) were investigated using the MP2 method at 6-311++G(2d,2p) basis set. Three different patterns including non-classical hydrogen bond (HB) H···C and classical HB H···X were found for complex formation between SY and HX molecules. Stability of the H···C type complexes are greater than H···X complexes. Quantum theories of atoms in molecules, natural bond orbitals, and energy decomposition analysis methods have been applied to analyze the intermolecular interactions. Good correlations have been found between the interaction energies (SE), the second-order perturbation energy E ⁽²⁾ and the charge transfer qCT in the studied systems.     

复合物HNO···H2O2内N-H···O蓝移氢键的理论研究

A theoretical study on the blue-shifted H-bond N-H…O and red-shifted H-bond O-H…O in the complexHNO…H_2O_2 was conducted by employment of both standard and counterpoise-corrected methods to calculate thegeometric structures and vibrational frequencies at the MP2/6-31G(d),MP2/6-31 G(d,p),MP2/6-311 q G(d,p),B3LYP/6-31G(d),B3LYP/6-31 G(d,p) and B3LYP/6-311 G(d,p) levels.In the H-bond N-H…O,the calcu-lated blue shift of N-H stretching frequency is in the vicinity of 120 cm~(-1) and this is indeed the largest theoreticalestimate of a blue shift in the X-H…Y H-bond ever reported in the literature.From the natural bond orbital analy-sis,the red-shifted H-bond O-H…O can be explained on the basis of the dominant role of the hyperconjugation.For the blue-shifted H-bond N-H…O,the hyperconjugation was inhibited due to the existence of significant elec-tron density redistribution effect,and the large blue shift of the N-H stretching frequency was prominently due tothe rehybridization of sp~n N-H hybrid orbital.     

Intramolecular hydrogen bonding in chemoselective synthesized 2-substituted pyrrole stable phosphorus ylide: GIAO, AIM, and NBO approaches

The chemoselectivity of geometrically ylide compounds is often hard to assign from experimental techniques, particular system with intramolecular hydrogen bonding (IHB) are even more challenging. Herein, theoretical calculations were performed to investigate whether theoretical results would provide consistent evidence for the existence of IHB to confirm experimental data and to evaluate strength of the N–H···O IHB from geometrical synthesized 2-substituted pyrrole stable phosphorus ylide (dimethyl 2-(1H-pyrrol-2-yl)-3-(triphenylphosphoranylidene) butanedioate in a single chemoselective compound. Topological parameters at the bond critical points (BCP) of intramolecular hydrogen bonds from Bader’s atoms in molecules (AIM) theory and Winhold’s natural bond orbital (NBO) calculations were analyzed at the B3LYP/6-311++g** level in details. A series of gage-including atomic orbital chemical shift (GIAO c.s.) calculations at the HF and DFT levels of theory were carried out to assign the ¹H NMR chemical shifts. The best prediction of the experimental ¹H NMR values was obtained at the mPW1PW91 levels using the 6-31G** basis set. Theoretical results, in agreement with the experimental data, were confirmed the N–H···O IHB was caused the deshielding of the proton to lower field. The barriers in double bond and single bond rotation were theoretically estimated in detailed and the AIM and NPA approaches were confirmed the loss of charge of the hydrogen atom involving in intramolecular N–H···O hydrogen bonding. The geometrical and topological parameters from AIM and NBO analyses were indicated the medium N–H···O IHB.     

Understanding hydrogen bonding of hydroxamic acids with some amino acid side chain model molecules

The hydrogen-bonding abilities of a few amino acid side chains have been studied through aggregation of methylamine, methanol, and acetic acid (as model molecules) with formo- and thioformo- hydroxamic acids using ab initio calculations. Forty six aggregates representing all possible H-bond interactions between these amino acid side chain groups and two most stable keto and enol tautomeric forms of both hydroxamic acids have been optimized. Although participation of conventional H-bond donors and acceptors leads to significant stabilization energies, yet C–H···O, C–H···N, S–H···O, and S–H···N etc. unconventional H-bonds also contribute to stabilize interactions in many aggregates. Strength of H-bonds of the molecules with formo- and thioformo- hydroxamic acid studied follows the order acetic acid > methylamine > methanol. A comparative study of atomic charges and orbital interactions employing NBO analysis has been carried out to explore the role of bond polarizations, charge transfer, and electron delocalizations as contributors to stabilization energy.     

气相中O3与HSO自由基间的氢键复合物

气相中O3与HSO自由基之间的相互作用及其反应在大气化学中非常重要. 在DFT-B3LYP/6-311++G**和MP2/6-311++G**水平上求得O3+HSO复合物势能面上的稳定构型, B3LYP方法得到了三种构型(复合物I, II和III), 而MP2方法只能得到一种构型(复合物II). 在复合物I和III中, HSO单元中的1H原子作为质子供体, 与O3分子中的端基O原子作为质子受体相互作用, 形成红移氢键复合物; 而在复合物II中, 虽与复合物I和III中具有相同的质子供体和质子受体, 却形成了蓝移氢键复合物. B3LYP/6-311++G**水平上计算的单体间相互作用能的计算考虑了基组重叠误差(BSSE)和零点振动能(ZPVE)校正, 其值在-3.37到-4.55 kJ·mol-1之间. 采用自然键轨道理论(NBO)对单体间相互作用的本质进行了考查, 并通过分子中原子理论(AIM)分析了三种复合物中氢键的电子密度拓扑性质.     

气相中O3与HSO自由基间的氢键复合物

气相中O3与HSO自由基之间的相互作用及其反应在大气化学中非常重要.在DFT-B3LYP/6-311++G**和MP2/6-311++G**水平上求得O3+HSO复合物势能面上的稳定构型,B3LYP方法得到了三种构型(复合物Ⅰ,Ⅱ和Ⅲ),而MP2方法只能得到一种构犁(复合物Ⅱ).在复合物Ⅰ和Ⅲ中,HSO单元中的1H原子作为质子供体.与O3分子中的端基O原子作为质子受体相互作用,形成红移氢键复合物;而在复合物Ⅱ中,虽与复合物Ⅰ和Ⅲ中具有相间的质子供体和质子受体,却形成了蓝移氢键复合物.B3LYP/6-311++G**水平上计算的单体间相互作用能的计算考虑了基组重甍误差(BSSE)和零点振动能(ZPVE)校正,其值在-3.37到-4.55 kJ·mol-1之间.采用自然键轨道理论(NBO)对单体间相互作用的本质进行了考查,并通过分子中原子理论(AIM)分析了三种复合物中氢键的电子密度拓扑性质.     

Characterization of intermolecular interactions in crystalline aspirin: A computational NQR study

A computational study at the level of density functional theory was carried out to characterize the ¹⁷O and ²H nuclear quadrupole resonance (NQR) spectroscopy parameters in crystalline aspirin. To include O? H···O and C? H···O hydrogen bonding effects in the calculations, the most probable interacting molecules with the target molecule in the crystalline phase were considered through a pentamer cluster. The NQR calculations were performed with BLYP, B3LYP, and M06 functionals employing 6‐311++G** and Jensen's polarization‐consistent pcJ‐1 basis sets. Linear correlations are observed between the calculated ¹⁷O and ²H NQR parameters and the hydrogen bond strengths, suggesting the possibility of estimating hydrogen bonding information from calculated NQR data. Different contributions of various nuclei to hydrogen bonding interactions and observed trends of calculated NQR parameters are well justified by atoms in molecules analyses at the BCPs of these interactions. © 2011 Wiley Periodicals, Inc. Int J Quantum Chem, 2012     

Ab initio study of water clustering in the presence of a methyl radical

The geometrical structure and binding energy of small clusters of methyl radical and water molecules (up to five water molecules) in gas phase and water media have been investigated at the MP2 level of theory using 6-311++G(2df,2p) basis set. The complexes characterized contain OH···O, CH···O, and OH···C attractive interactions with stabilization energies in the range 6–143 kJ mol^?1. The solvent has an enhancing influence on the stabilities of studied clusters. The atoms in molecules theory were also applied to explain the nature of the complexes. The interaction energies have been partitioned with the natural energy decomposition analysis showing that the most important attractive term corresponds to the charge transfer one.     

Theoretical Studies on Energetic Property and Stability of Pyrazine and Pyridine Derivatives

Density functional calculations at the B3LYP level with 6‐311G** and aug‐cc‐pVDZ basis sets were performed to predict the heats of formation (HOFs) for two pyrazine derivatives and eight pyridine derivatives. In the isodesmic reactions designed for the computation of heats of formation (HOFs), pyrazine and pyridine were chosen as reference compounds. The N‐oxidations for the ring nitrogen of pyrazine and pyridine derivatives decrease the HOF values when N‐oxide oxygen is neighboring with amino groups, but increase when it neighbors with nitro groups. Thermal stability was evaluated via bond dissociation energies (BDE) at the UB3LYP/6‐311G** level. As a whole, the homolysis of C–NO₂ bonds is the main step for bond dissociation of the title compounds. The BDE values of title compounds are influenced by intramolecular hydrogen bonds. The hydrogen bond effects associated with the length of the H···O bonds were analyzed by the electron density at the critical points and natural bond orbital.