Density Functional Theory Study on Interaction between Catechin and Thymine

The interacting patterns and mechanism of the catechin and thymine have been investigated with the density functional theory Becke's three-parameter nonlocal exchange functional and the Lee, Yang, and Parr nonlocal correlation functional (B3LYP) method by 6-31+G*basis set. Thirteen stable structures for the catechin-thymine complexes have been found which form two hydrogen bonds at least. The vibrational frequencies are also studied at the same level to analyze these complexes. The results indicated that catechin interactedwith thymine by three different hydrogen bonds as N-H…O、C-H…O、O-H…O and the complexes are mainly stabilized by the hydrogen bonding interactions. Theories of atoms in molecules and natural bond orbital have been adopted to investigate the hydrogen bondsinvolved in all systems. The interaction energies of all complexes have been corrected for basis set superposition error, which are from -18.15 kJ/mol to -32.99 kJ/mol. The results showed that the hydrogen bonding contribute to the interaction energies dominantly. The corresponding bonds stretching motions in all complexes are red-shifted relative to that of the monomer, which is in agreement with experimental results.     

Catalytic peroxide oxidation: The structure of key intermediates in the V<Superscript>V</Superscript>/H<Subscript>2</Subscript>O<Subscript>2</Subscript> system according to quantum chemical data

The inner-sphere isomerization of the peroxo complexes of vanadium(V) with the general formula [VO₆]^? was studied using approximations based on the density functional theory (B3LYP/6-31G**) and the Møller-Plesset perturbation theory (MP2/6-31G**). It was found that the complex [V(=O)(ηO₂)(O₃)]^? containing the O₃ group as a bidentate ligand was the most stable isomer. The transition state region of a rear-rangement of the triperoxo complex [V(ηO₂)₃]^? into [V(=O)(ηO₂)(O₃)]^? was localized. It was found that the activation barrier (～30 kcal/mol) was mainly due to O-O bond cleavage in the peroxo ligand. According to calculations, the reaction proceeds through two intermediate complexes whose structure can be interpreted as that containing coordinated singlet dioxygen (especially in the limiting case) because of noticeably shortened O-O bonds in the ηO₂ ligand. The calculated reaction scheme of the conversion of [V(ηO₂)₃]^? into [V(=O)(ηO₂)(O₃)]^? is qualitatively consistent with the previously found kinetics of the formation of ozone and the oxidation of alkanes, olefins, arenes, and singlet dioxygen traps.     

Conformational Study of α‐N‐Acetyl‐D‐Neuraminic Acid by Density Functional Theory

The stable structures of α‐N‐acetyl‐D‐neuraminic acid (Neu5Acα) in the gas phase were studied at the B3LYP level of theory using 6‐31G(d,p) and 6‐31++G(d,p) basis sets. They are classified into five types according to the patterns of the intramolecular hydrogen bond formations. One of the stable structures had intramolecular hydrogen bond network of O₉H_O9 … O₈H_O8 … O?C₁‐O₁H_O1 and O₇H_O7…O?CHN‐C₅ similar to the crystal structure of Neu5Ac‐α‐methyl glycoside methyl ester. The stable structures of Neu5Acα are reasonable for the following sialooligosaccharide ligand studies with respect to the relationship between OH group orientations and intramolecular hydrogen bond formations. The barrier heights for isomerizations between the stable structures were computed to be 2.8 to 6.7 kcal/mol at the B3LYP/6‐31++G(d,p)//B3LYP/6‐31G(d,p) level, which are basic factors for the conformational behavior of Neu5Acα before its interactions with receptors. We also calculated Neu5Acα–4 or 5‐water complexes to take account of the solvent effect on the intramolecular hydrogen bonds in the stable structures. Consequently, the structures of Neu5Acα in the complexes are similar to each other, which is consistent with the known NMR data. Thus, the optimum Neu5Acα‐water complexes are some of the reasonable pseudohydrous Neu5Acα.     

Density Functional Theory Study on Hydrogen Bonding Interaction of Catechin-(H2O)n

Density functional theory B3LYP method with 6‐31G* basis set has been used to optimize the geometries of the catechin, water and catechin‐(H₂O)_n complexes. The vibrational frequencies have been studied at the same level to analyze these complexes. Six and eleven stable structures for the catechin‐H₂O and catechin‐(H₂O)₂ have been found, respectively. Theories of atoms in molecules (AIM) and natural bond orbital (NBO) have been utilized to investigate the hydrogen bonds involved in all the systems. The interaction energies of all the complexes corrected by basis set superposition error, are from ?13.27 to ?83.56 kJ/mol. All calculations also indicate that there are strong hydrogen‐bonding interactions in catechin‐water complexes. The strong hydrogen‐bonding contributes to the interaction energies dominantly. The O–H stretching motions in all the complexes are red‐shifted relative to that of the monomer.     

Theoretical study on stabilities of multiple hydrogen bonded dimers

A series of self-constituted multiple hydrogen bonded (MHB) complexes has been investigated systematically by density functional theory (PBE1PBE /6-31G**), the Morokuma energy decomposition method (HF/6-31G**) and MP2 (6-31G** and 6-311++G**) calculation. We have discovered that (i) for doubly hydrogen bonded (DHB) complexes, both the interaction energy and stability increase with the charge transfer energy; (ii) for quadruple hydrogen bonded (QHB) complexes, cooperativity is the most important factor determining stability of the complex: stronger cooperative energy correlates well with larger interaction energy and thus more stable complex and vice versa; (iii) correlation energy plays an important role in intermolecular interactions. The correlation energy, mainly consisting of dispersive energy, also exhibits cooperativity in MHB dimers: positive for M-aadd and generally negative for other complexes.     

Theoretical study of the N-H…O red-shifted and blue-shifted hydrogen bonds

Theoretical calculations are performed to study the nature of the hydrogen bonds in complexes HCHO…HNO, HCOOH…HNO, HCHO…NH3, HCOOH…NH3, HCHO…NH2F and HCOOH…NH2F. The geometric structures and vibrational frequencies of these six complexes at the MP2/6-31+G(d,p), MP2/6-311++G(d,p), B3LYP/6-31+G(d,p) and B3LYP/6-311++G(d,p) levels are calculated by standard and counterpoise-corrected methods, respectively. The results indicate that in complexes HCHO…HNO and HCOOH…HNO the N-H bond is strongly contracted and N-H…O blue-shifted hydrogen bonds are observed. While in complexes HCHO…NH3, HCOOH…NH3, HCHO…NH2F and HCOOH…NH2F, the N-H bond is elongated and N-H…O red-shifted hydrogen bonds are found. From the natural bond orbital analysis it can be seen that the X-H bond length in the X-H…Y hydrogen bond is controlled by a balance of four main factors in the opposite directions: hyperconjugation, electron density redistribution, rehybridization and structural reorganization. Among them hyperconjugation has the effect of elongating the X-H bond, and the other three factors belong to the bond shortening effects. In complexes HCHO…HNO and HCOOH…HNO, the shortening effects dominate which lead to the blue shift of the N-H stretching frequencies. In complexes HCHO…NH3, HCOOH…NH3, HCHO…NH2F and HCOOH…NH2F where elongating effects are dominant, the N-H…O hydrogen bonds are red-shifted.     

Modelling protein-RNA interactions: an electron density study of the guanidinium and formate complexes with RNA bases

The complexes formed by the double interaction established between RNA bases and guanidinium and formate ions, as a model for the interacting groups of arginine and glutamic or aspartic amino acid side chains, have been theoretically studied. A density functional theory method (B3LYP/6-31 + G**) has been used for this study. The range of interaction energies obtained allowed for a distinction between bidentate and bifurcate hydrogen bond interactions. The analysis of the electron density and the natural bond orbital analysis shows that these complexes are bound by double hydrogen bonds established between the donor and acceptor groups of guanidinium and formate respectively and those of the RNA bases. Comparisons are made with the results obtained in some previous theoretical and experimental studies.     

Modeling of hydrogen bonds in monohydrated 2,4-dithiothymine: an ab initio and AIM study

Twelve tautomers of 2,4-dithiothymine are calculated at the MP2/6-31+G(d) level, and the most stable one is referred to the di-keto form (P12). Then four H-bonded complexes between P12 and water are optimized at the MP2/6-31+G(d) level of theory. The calculation of vibrational frequencies and natural bond orbital analysis are also carried out at the same level to investigate the hydrogen bonds involved in all the systems. Within all the four complexes, three types of hydrogen bonds are formed, in which the O-H...S and N-H...O bonds are the normal bonds with the X-H bond elongation and red shift of the corresponding stretch frequencies, while the C-H...O interaction is an improper, blue-shifting hydrogen bond accompanied with the contraction of the C-H bond and a blue shift of the C-H stretch frequency. The topological properties are investigated with the atoms-in-molecules (AIM) theory. The NMR chemical shielding for the isolated and the four monohydrated 2,4-dithiothymine are calculated using the "gauge-including atomic orbital" (GIAO) method. The 1H chemical shifts are influenced by the formation of hydrogen bonds.     

Intermolecular interactions in uracil–nitrous acid complexes: structures,binding energy,topological properties,and nuclear magnetic resonance study

Molecular interactions between uracil and nitrous acid (U–NA) [C₄N₂O₂H₄? NO₂H] have been studied using B3LYP, B3PW91, and MP2 methods with different basis sets. The optimized geometries, harmonic vibrational frequencies, charge transfer, topological properties of electron density, nucleus‐independent chemical shift (NICS), and nuclear magnetic resonance one‐ and two‐bonds spin–spin coupling constants were calculated for U–NA complexes. In interaction between U and NA, eight cyclic complexes were obtained with two intermolecular hydrogen bonds N(C)HU…N(O) and OH_NA…O_U. In these complexes, uracil (U) simultaneously acts as proton acceptor and proton donor. The most stable complexes labeled, UNA1 and UNA2, are formed via NH bond of U with highest acidity and CO group of U with lowest proton affinity. There is a relationship between hydrogen bond distances and the corresponding frequency shifts. The solvent effect on complexes stability was examined using B3LYP method with the aug‐cc‐pVDZ basis set by applying the polarizable continuum model (PCM). The binding energies in the gas phase have also been compared with solvation energies computed using the PCM. Natural bond orbital analysis shows that in all complexes, the charge transfer takes place from U to NA. The results predict that the Lone Pair (LP)(O)_U → σ*(O? H) and LP(N(O)_NA → σ*(N(C)? H)_U donor–acceptor interactions are most important interactions in these complexes. Atom in molecule analysis confirms that hydrogen bond contacts are electrostatic in nature and covalent nature of proton donor groups decreases upon complexation. The relationship between spin–spin coupling constant (^1hJ_H…_Y and ^2hJ_H…_Y) with interaction energy and electronic density at corresponding hydrogen bond critical points and H‐bonds distances are investigated. NICS used for indicating of aromaticity of U ring upon complexation. © 2013 Wiley Periodicals, Inc.     

乙醛二聚体的从头计算

用ab initio分子轨道法在STO-3G和6-31G水平上, 全构型优化, 对二聚乙醛的各种稳定构型进行了研究。结果表明, 在STO/6-31G水平上, 最稳定的乙醛二聚体为具有对称中心的环状结构, 包含由醛基氢和醛基氧组成的二个C—H…O氢键。结合能为-20.17 kJ·mol~(-1), 与实验估计的结合能-22.39±0.15 kJ·mol~(-1)比较接近。STO-3G过低估计了乙醛二聚体的结合能, 在不同构型的相对稳定性方面也与6-31G不一致。     

Structure and proton-donor ability of N-(1-trifluoromethylsulfonylamino-2,2,2-trichloroethyl)-acrylamide

Quantum-chemical calculations (B3LYP/6-311G**) of N-(1-trifluoromethylsulfonylamino-2,2,2-trichloroethyl)acrylamide CF₃SO₂NHCH(CCl₃)NHC(O)CH=CH₂ (I) in the isolated state revealed four local minima corresponding to the conformers with the syn- and antiperiplanar orientation of the C=O and N-H bonds in the amide fragment, two of which containing the intramolecular C(O)NH…O=S or SO₂NH…O=C hydrogen bonds. Judged from the data of IR spectroscopy and dielectrometry, compound I in inert media exists predominantly in the form of conformer with antiperiplanar amide fragment and free NH group. Its self-associates in molecular crystals and solutions are formed by hydrogen bonds SO₂NH…O=C. Spectroscopic acidity of compound I determined as the value of Δν(NH) upon interactions with DMF in CCl₄ is higher than that of N-methyltrifluoromethanesulfonamide.     

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纳米。     

Theoretical study on coordination of methanol clusters to 3-methyl-4-pyrimidone

Density functional theory (DFT), MP2, and couple cluster ab initio methods were employed to investigate the microsolvation of 3-methyl-4-pyrimidone (3M4P) surrounded by methanol (MeOH) molecules. Structures are analyzed based on hydrogen bonds with a focus on relative energies, interaction energies, hydrogen bond cooperativity, hydrogen bonding geometries, and redshifts in the frequencies of O–H and C=O stretching modes. Our results show that there is no preferential orientation of MeOH attacks on the carbonyle site of 3M4P; both trans and cis 3M4P-MeOH complexes have same chance to be observed. cis 3M4P-MeOH and 3M4P-MeOH complex in which MeOH is located on N lie 0.56 and 3.11 kJ/mol at CCSD(T)/6-31+G(d,p) (0.63 and 1.67 kJ/mol at MP2/6-311++G(d,p)) above trans 3M4P-MeOH. MeOH dimers form more stable 3M4P-(MeOH)₂ complexes compare to 3M4P-(MeOH)₂ complexes in which individual MeOH molecules bind to carbonyl and N. Relative energies of 3M4P-(MeOH)₃ computed using various DFT methods point out the complex formed by linear MeOH trimer along methyl group of 3M4P (cis 3M4P-(MeOH)₃) as lowest. Carbonyl group is predicted as preferential site for hydrogen bond interaction. Besides O–H…O and O–H…N hydrogen bonds, 3M4P-(MeOH)₂ and 3M4P-(MeOH)₃ complexes are also stabilized by H–O…H–C weak interactions.     

DFT study on the intermolecular interactions between Aun (n = 2–4) and thymine

Density functional B3LYP method with 6-31++G** basis set is applied to optimize the geometries of the luteolin, water and luteolin–(H₂O)_n complexes. The vibrational frequencies are also studied at the same level to analyze these complexes. We obtained four steady luteolin–H₂O, nine steady luteolin–(H₂O)₂ and ten steady luteolin–(H₂O)₃, respectively. Theories of atoms in molecules (AIM) and natural bond orbital (NBO) are used to investigate the hydrogen bonds involved in all the systems. The interaction energies of all the complexes corrected by basis set superposition error, are within −13.7 to −82.5 kJ/mol. The strong hydrogen bonding mainly contribute to the interaction energies, Natural bond orbital analysis is performed to reveal the origin of the interaction. All calculations also indicate that there are strong hydrogen bonding interactions in luteolin–(H₂O)_n complexes. The OH stretching modes of complexes are red-shifted relative to those of the monomer.     

Theoretical study of the red‐ and blue‐shifted hydrogen bonds of nitroxyl and acetylene dimers

Ab initio molecular orbital and density functional theory (DFT) in conjunction with different basis sets calculations were performed to study the C? H…O red‐shifted and N? H…π blue‐shifted hydrogen bonds in HNO? C₂H₂ dimers. The geometric structures, vibrational frequencies and interaction energies were calculated by both standard and counterpoise (CP)‐corrected methods. In addition, the G3B3 method was employed to calculate the interaction energies. The topological and natural bond orbital (NBO) analysis were investigated the origin of N? H…π blue‐shifted hydrogen bond. From the NBO analysis, the electron density decrease in the σ* (N? H) is due to the significant electron density redistribution effect. The blue shifts of the N? H stretching frequency are attributed to a cooperative effect between the rehybridization and electron density redistribution. © 2006 Wiley Periodicals, Inc. Int J Quantum Chem, 2006     

水溶液中氨基酸侧链对G∶C碱基对间氢键影响的理论研究

利用量子化学方法研究了气相和水溶液下,氨基酸侧链与鸟嘌呤和胞嘧啶间的氢键作用.应用B3LYP/6-31+G(d,p)方法优化复合物几何结构,使用MP2/aug-cc-p VDZ方法进行复合物能量、自然键轨道(NBO)电荷和二阶稳定化能的计算.结果表明,水溶液可使氨基酸侧链与碱基或碱基对之间氢键键能显著减小;带电复合物气相和水溶液氢键键能之差范围为50.63~146.48 k J/mol,中性为0.17~24.94 k J/mol;电荷的转移量与氢键键能成正比,电荷转移量越多,复合物越稳定;二阶稳定化能与氢键键长成反比,与电荷转移量成正比,且气相与水溶液氢键二阶稳定化能之比约为两相的电荷转移量之比.水溶液对该类体系中氢键作用具有明显影响.     

Study of hydrogen bonding interactions relevant to biomolecular structure and function

Ab initio molecular orbital calculations were used to study hydrogen bonding interactions and interatomic distances of a number of hydrogen bonded complexes that are germane to biomolecular structure and function. The calculations were carried out at the STO-3G, 3-21G, 6-31G*, and MP2/6-31G* levels (geometries were fully optimized at each level). For anionic species, 6-31 + G* and MP2/6-31 + G* were also used. In some cases, more sophisticated calculations were also carried out. Whenever possible, the corresponding enthalpy, entropy, and free energy of complexation were calculated. The agreement with the limited quantity of experimental data is good. For comparison, we also carried out semiempirical molecular orbital calculations. In general, AM1 and PM3 give lower interaction enthalpies than the best ab initio results. With regard to structural results, AM1 tends to favor bifurcated structures for O? H-O and N? HO types of hydrogen bonds, but not for hydrogen bonds involving O-H? S and S-H? O, where the usual hydrogen bond patterns are observed. Overall, AM1 geometries are in general in poor agreement with ab initio structural results. On the other hand, PM3 gives geometries similar to the ab initio ones. Hence, from the structural point of view PM3 does show some improvement over AM1. Finally, insights into the formation of cyclic or open formate–water hydrogen bonded complexes are presented. © 1992 by John Wiley & Sons, Inc.     

Theoretical Study on the Hydrogen Bonding Interactions in Complexes of 5‐Hydroxytryptamine with Water

The energies, geometries and harmonic vibrational frequencies of 1:1 5‐hydroxytryptamine‐water (5‐HT‐H₂O) complexes are studied at the MP2/6‐311++G(d,p) level. Natural bond orbital (NBO), quantum theory of atoms in molecules (QTAIM) analyses and the localized molecular orbital energy decomposition analysis (LMO‐EDA) were performed to explore the nature of the hydrogen‐bonding interactions in these complexes. Various types of hydrogen bonds (H‐bonds) are formed in these 5‐HT‐H₂O complexes. The intermolecular C4H5^5‐HT···O^w H‐bond in HTW3 is strengthened due to the cooperativity, whereas no such cooperativity is found in the other 5‐HT‐H₂O complexes. H‐bond in which nitrogen atom of amino in 5‐HT acted as proton donors was stronger than other H‐bonds. Our researches show that the hydrogen bonding interaction plays a vital role on the relative stabilities of 5‐HT‐H₂O complexes.     

Theoretical study of the interaction of fluorinated dimethyl ethers and the ClF and HF molecules. Comparison between halogen and hydrogen bonds

A theoretical study of the halogen‐bonded complexes formed between fluorinated dimethyl ethers (n_F = 0–4) and ClF is carried out using the wB97XD method combined with the 6‐311++G(d,p) basis set. The properties of the complexes are compared with the corresponding properties of the hydrogen‐bonded complexes formed between the same electron donors and HF. The optimized geometries, the interaction energies, relevant natural bonding orbital characteristics along with some vibrational data are calculated. The analyzed properties also include the symmetry adapted perturbation theory decomposition of the energies along with the atoms‐in molecule analysis. For both the halogen and hydrogen bonds, the interaction energies are ruled by the intermolecular hyperconjugation energies. In contrast, the correlations between the binding energies and the basic properties of the ethers or the charge transfer are different for the halogen and hydrogen bonds. The applicability of the Bent's rule to these systems is discussed. © 2016 Wiley Periodicals, Inc.     

Density Functional Theory Study of the Interaction between Guanine and Catechin

The interacting patterns and mechanism of the catechin and guanine have been investigated with the density functional theory B3LYP method by 6‐31G* basis set. Fourteen stable structures for the catechin‐guanine complexes have been found which form two hydrogen bonds at least. The results indicate that the complexes are mainly stabilized by the hydrogen bonding interactions. At the same time, the number and strength of hydrogen bond play a co‐determinant parts in the stability of the complexes which can form two or more hydrogen bonds. Theories of atoms in molecules (AIM) and natural bond orbital (NBO) have been adopted to investigate the hydrogen bonds involved in all systems. The interaction energies of all complexes have been corrected for basis set superposition error (BSSE), ranging from ?38.86 to ?14.56 kJ/mol. The results showed that the hydrogen bonding contributes to the interaction energies dominantly. The corresponding bonds stretching motions in all complexes are red‐shifted relative to that of the monomer, which is in agreement with experimental results.