Molecular Tailoring Approach for the Estimation of Intramolecular Hydrogen Bond Energy

Hydrogen bonds (HBs) play a crucial role in many physicochemical and biological processes. Theoretical methods can reliably estimate the intermolecular HB energies. However, the methods for the quantification of intramolecular HB (IHB) energy available in the literature are mostly empirical or indirect and limited only to evaluating the energy of a single HB. During the past decade, the authors have developed a direct procedure for the IHB energy estimation based on the molecular tailoring approach (MTA), a fragmentation method. This MTA-based method can yield a reliable estimate of individual IHB energy in a system containing multiple H-bonds. After explaining and illustrating the methodology of MTA, we present its use for the IHB energy estimation in molecules and clusters. We also discuss the use of this method by other researchers as a standard, state-of-the-art method for estimating IHB energy as well as those of other noncovalent interactions.     

Theoretical Studies on the Stabilities and Hydrogen Bond Actions of （H2O）n Clusters

1 INTRODUCTION In the latest ten years, the structure and function of water clusters have captured the interest of chemists. One of the most important study objects in water cluster is to describe the behavior of water so- lution quantitatively at molecule level, which will pave the way for the solving of some environmental and other scientific problems, such as the formation of acid rain and nucleation mechanism of little water drop. Besides, weak interaction in water clusters could be al…     

Hydrogen and Halogen Bonding in Homogeneous External Electric Fields: Modulating the Bond Strengths

Recent years have witnessed various fascinating phenomena arising from the interactions of noncovalent bonds with homogeneous external electric fields (EEFs). Here we performed a computational study to interpret the sensitivity of intrinsic bond strengths to EEFs in terms of steric effect and orbital interactions. The block-localized wavefunction (BLW) method, which combines the advantages of both ab initio valence bond (VB) theory and molecular orbital (MO) theory, and the subsequent energy decomposition (BLW-ED) approach were adopted. The sensitivity was monitored and analyzed using the induced energy term, which is the variation in each energy component along the EEF strength. Systems with single or multiple hydrogen (H) or halogen (X) bond(s) were also examined. It was found that the X-bond strength change to EEFs mainly stems from the covalency change, while generally the steric effect rules the response of H-bonds to EEFs. Furthermore, X-bonds are more sensitive to EEFs, with the key difference between H- and X-bonds lying in the charge transfer interaction. Since phenylboronic acid has been experimentally used as a smart linker in EEFs, switchable sensitivity was scrutinized with the example of the phenylboronic acid dimer, which exhibits two conformations with either antiparallel or parallel H-bonds, thereby, opposite or consistent responses to EEFs. Among the studied systems, the quadruple X-bonds in molecular capsules exhibit remarkable sensitivity, with its interaction energy increased by −95.2 kJ mol⁻¹ at the EEF strength 0.005 a.u.     

Rapid Prediction of the Hydrogen Bond Cooperativity in N‐methylacetamide Chains

A method is proposed to rapidly predict the hydrogen bond cooperativity in N‐methylacetamide chains. The parameters needed are obtained from the fittings to the hydrogen bonding energies in the formamide chains containing 2 to 8 monomeric units. The scheme is then used to calculate the individual hydrogen bonding energies in N‐methylacetamide chains containing 2 to 7 monomeric units. The cooperativity predicted is in good agreement with those obtained from MP2/6‐31+G** calculations by including the BSSE correction. Our scheme is further employed to predict the individual hydrogen bonding energies in larger N‐methylacetamide chains containing up to 200 monomeric N‐methylacetamide units, to which the MP2 method cannot be applied. Based on our scheme, a cooperative effect of over 170 % of the dimer hydrogen bonding energy in long N‐methylacetamide chains is predicted. The method is also applied to heterogeneous chains containing formamide, acetamide, N‐methylformamide, and N‐methylacetamide. The individual hydrogen bonding energies in these heterogeneous chains are also in good agreement with those obtained from MP2 calculations with the BSSE correction, further demonstrating that our method is reasonable.     

纯水溶液中的氢键交换反应路径

利用分子动力学模拟方法对纯水溶液的氢键转化动力学性质进行了深入的微观探讨,溶液中非氢键构型为寿命较短(0.1～0.2 ps)的过渡态构型,我们发现氢键交换通过两种过渡构型完成,氢键角度扭曲激发后与氢键第一壳层水分子沿路径1交换,氢键径向拉伸激发后与氢键第二壳层水分子沿路径2交换,过渡态路径的选择具有温度依赖性.氢键转化需在旧氢键受体氢键过量和新氢键受体氢键不足,同时满足交换反应空间结构要求下才能完成.氢键交换反应对水分子平动和转动行为起着决定作用.     

氢键在现代化学中的作用

在各种分子间相互作用中,氢键占有很特殊的地位,被称作为超分子化学中的万能相互作用。讨论了氢键在超分子、自组装、分子识别、晶体工程、材料化学和催化过程等现代化学领域中的作用。     

Hydrogen Bond Cooperativity and the Three‐Dimensional Structures of Water Nonamers and Decamers

Broadband rotational spectroscopy of water clusters produced in a pulsed molecular jet expansion has been used to determine the oxygen atom geometry in three isomers of the nonamer and two isomers of the decamer. The isomers for each cluster size have the same nominal geometry but differ in the arrangement of their hydrogen bond networks. The nearest neighbor O? O distances show a characteristic pattern for each hydrogen bond network isomer that is caused by three‐body effects that produce cooperative hydrogen bonding. The observed structures are the lowest energy cluster geometries identified by quantum chemistry and the experimental and theoretical O? O distances are in good agreement. The cooperativity effects revealed by the hydrogen bond O? O distance variations are shown to be consistent with a simple model for hydrogen bonding in water that takes into account the cooperative and anticooperative bonding effects of nearby water molecules.     

Modulation of hydrogen bonding upon ion binding: Insights into cooperativity

The impact due to the of presence of ions, such as Mg²⁺, Na⁺, H⁺, Cl^?, and OH^? on hydrogen bonded clusters of increasing size (water, formamide, and acetamide [n = 1–10]) in the context of associated cooperativity has been explored using density functinal theory (DFT) calculations. Sequential binding energies (SBE) rise on addition of monomer in case of parent clusters. SBE for ionic clusters are several times higher than that of parent clusters initially. This behavior is more dramatic on addition of either Mg²⁺ or H⁺ compared to other ions. Interestingly, SBE of both parent and ionic clusters approach nearly uniform values beyond n = 6 irrespective of kind of ion present in the cluster with the exception of magnesium. © 2013 Wiley Periodicals, Inc.     

取代基对1-甲基尿嘧啶与N-甲基乙酰胺氢键复合物中氢键强度的影响

使用密度泛函理论B3LYP方法和二阶微扰理论MP2方法对由1-甲基尿嘧啶与N-甲基乙酰胺所形成的氢键复合物中的氢键强度进行了理论研究, 探讨了不同取代基取代氢键受体分子1-甲基尿嘧啶中的氢原子对氢键强度的影响和氢键的协同性. 研究表明: 供电子取代基使N－H…O＝C氢键键长r(H…O)缩短, 氢键强度增强; 吸电子取代基使N－H…O＝C氢键键长r(H…O)伸长, 氢键强度减弱. 自然键轨道(NBO)分析表明: 供电子基团使参与形成氢键的氢原子的正电荷增加, 使氧原子的负电荷增加, 使质子供体和受体分子间的电荷转移量增多; 吸电子基团则相反. 供电子基团使N－H…O＝C氢键中氧原子的孤对电子轨道n(O)对N－H的反键轨道σ^*(N－H)的二阶相互作用稳定化能增强, 吸电子基团使这种二阶相互作用稳定化能减弱. 取代基对与其相近的N－H…O＝C氢键影响更大.     

Making and breaking of small water clusters: A combined quantum chemical and molecular dynamics approach

We present a combined quantum chemical and molecular dynamics study of cyclic and noncyclic water n-mers ([(H₂O]_n, n = 2–6) at four different temperatures and showcase that the dynamics of small water clusters can reproduce the known properties of bulk water reasonably well. We investigate the making and breaking of the water clusters by computing the hydrogen bond strengths, average lifetimes, and relative stabilities, which are important to understand the complex solution dynamics. We compare the behavior of water clusters in the gas phase and in the solution phase as well as the variation in the properties as a function of cluster size and highlight the notably more interesting cluster dynamics of the water trimer when compared to the other water clusters. © 2019 Wiley Periodicals, Inc.     

The Role of Hydrogen Bond in Catalytic Triad of Serine Proteases

In order to investigate the origin of catalytic power for serine proteases, the role of the hydrogen bond in the catalytic triad was studied in the proteolysis process of the peptides chymotrypsin inhibitor 2 (CI2), MCTI-A, and a hexapeptide (SUB), respectively. We first calculated the free energy profile of the proton transfer between His and Asp residues of the catalytic triad in the enzyme-substrate state and transition state by employing QM/MM molecular dynamics simulations. The results show that a low-barrier hydrogen bond (LBHB) only forms in the transition state of the acylation of CI2, while it is a normal hydrogen bond in the acylation of MCTI-A or SUB. In addition, the change of the hydrogen bond strength is much larger in CI2 and SUB systems than in MCTI-A system, which decreases the acylation energy barrier significantly for CI2 and SUB. Clearly, a LBHB formed in the transition state region helps accelerate the acylation reaction. But to our surprise, a normal hydrogen bond can also help to decrease the energy barrier. The key to reducing the reaction barrier is the increment of hydrogen bond strength in the transition state state, whether it is a LBHB or not. Our studies cast new light on the role of the hydrogen bond in the catalytic triad, and help to understand the catalytic triad of serine proteases.     

Elucidating Solvent Effects on Strong Intramolecular Hydrogen Bond: DFT-MD Study of Dibenzoylmethane in Methanol Solution

The dynamic aspect of solvation plays a crucial role in determining properties of strong intramolecular hydrogen bonds since solvent fluctuations modify instantaneous hydrogen-bonded proton transfer barriers. Previous studies pointed out that solvent-solute interactions in the first solvation shell govern the position of the proton but the ability of the electric field due to other solvent molecules to localize the proton remains an important issue. In this work, we examine the structure of the O−H⋅⋅⋅O intramolecular hydrogen bond of dibenzoylmethane in methanol solution by employing density functional theory-based molecular dynamics and quantum chemical calculations. Our computations showed that homogeneous electric fields with intensities corresponding to those found in polar solvents are able to considerably alter the proton transfer barrier height in the gas phase. In methanol solution, the proton position is correlated with the difference in electrostatic potentials on the oxygen atoms of dibenzoylmethane even when dibenzoylmethane-methanol hydrogen bonding is lacking. On a timescale of our simulation, the hydrogen bonding and solvent electrostatics tend to localize the proton on different oxygen atoms. These findings provide an insight into the importance of the solvent electric field on the structure of a strong intramolecular hydrogen bond.     

Theoretical Study of the Interplay between Lithium Bond and Hydrogen Bond in Complexes Involved with HLi and HCN

The lithium‐ and hydrogen‐bonded complex of HLi? NCH? NCH is studied with ab initio calculations. The optimized structure, vibrational frequencies, and binding energy are calculated at the MP2 level with 6‐311++G(2d,2p) basis set. The interplay between lithium bonding and hydrogen bonding in the complex is investigated with these properties. The effect of lithium bonding on the properties of hydrogen bonding is larger than that of hydrogen bonding on the properties of lithium bonding. In the trimer, the binding energies are increased by about 19 % and 61 % for the lithium and hydrogen bonds, respectively. A big cooperative energy (?5.50 kcal mol^?1) is observed in the complex. Both the charge transfer and induction effect due to the electrostatic interaction are responsible for the cooperativity in the trimer. The effect of HCN chain length on the lithium bonding has been considered. The natural bond orbital and atoms in molecules analyses indicate that the electrostatic force plays a main role in the lithium bonding. A many‐body interaction analysis has also been performed for HLi? (NCH)_N (N=2–5) systems.     

含有氢键供体大环化合物的构筑及其功能研究进展

含有氢键供体基团构筑的大环化合物因其结构中具有可以提供氢键供体的N—H基团,可以为大环化合物的主客体化学提供额外的分子间作用力,在分子识别、自组装以及超分子催化等领域被广泛应用.综述了近十年基于(硫)脲键、酰胺键构筑的大环化合物的合成方法及其在分子识别中的最新研究进展.为今后此类大环化合物的合成及应用提供参考.     

Impact of Covalent Modifications on the Hydrogen Bond Strengths in Diaminotriazine Supramolecules

Melamine (M) is a popular triamine triazine compound in the field of supramolecular materials. In this work, we have computationally investigated how substituents can be exploited to improve the binding strength of M supramolecules. Two types of covalent modifications were studied: the substitution of an H atom within an amine group −NHR, and the replacement of the whole −NH₂ group (R=H, F, CH₃ and COCH₃). Through our dispersion-corrected density functional theory computations, we explain which covalent modification will show the best self-assembling capabilities, and why the binding energy is enhanced. Our charge density and molecular orbital analyses indicate that the best substituents are those that generate a charge accumulation on the endocyclic N atom, providing an improvement of the electrostatic attraction. At the same time the substituent assists the main N−H⋅⋅⋅N hydrogen bonds by interacting with the amino group of the other monomer. We also show how the selected group notably boosts the strength of hexameric rosettes. This research, therefore, provides molecular tools for the rational design of emerging materials based on uneven hydrogen-bonded arrangements.     

Molecular processes in an aqueous cluster

A cluster of 216 water molecules in a vacuum is simulated (molecular-dynamic trajectory is 40 ps); the average kinetic energy of water molecule movements corresponds to 300 K. The HLC model of data representation in time intervals of 0.082 ps is used to investigate the molecular scale processes. Volumetric and surface movements of water molecules are distinguished by the rank distributions based on the molecular movement character. The method for calculating volumetric and surface self-diffusion coefficients of molecules in a cluster is proposed. The cooperative character of movements of water molecules in a 5.4 Å-thick surface layer of the cluster is revealed.     

A Mixed Quantum-Classical Approach for Computing Effects of Intramolecular Motion on the Low-Barrier Hydrogen Bond

The fractionation factor is defined as the equilibrium constant for the reaction: R – H + DOH R – D + HOH. Of interest are values of fractionation factors for reactions where reactants and/or products form intramolecular low-barrier hydrogen bonds. Experimentally measured isotopic fractionation factors are usually interpreted via a one-dimensional potential energy surface along the intrinsic proton hydrogen bond coordinate. Such a one-dimensional picture cannot be completely correct. Intramolecular motions, such as vibrations and librations, can modulate the underlying potential energy surface along the hydrogen bond coordinate and thus affect the isotopic fractionation factor. We have recently generated a picture of the motion of the proton in a low-barrier hydrogen bond as taking place in an effective single-dimensional potential, which we term the potential of mean force (PMF). In this paper, we compute the PMF for a molecule with an intramolecular hydrogen bond in order to quantify the effect of intramolecular motions on the fractionation factor. The PMF and isotopic fractionation factor are computed with a combination of high-level density functional theory and molecular dynamics simulations.     

Clusters of Ammonium Cation—Hydrogen Bond versus σ‐Hole Bond

MP2/6‐311++G(d,p) calculations were performed on the NH₄⁺ ??? (HCN)_n and NH₄⁺ ??? (N₂)_n clusters (n=1–8), and interactions within them were analyzed. It was found that for molecules of N₂ and HCN, the N centers play the role of the Lewis bases, whereas the ammonium cation acts as the Lewis acid, as it is characterized by sites of positive electrostatic potential, that is, H atoms and the sites located at the N atom in the extension of the H?N bonds. Hence, the coordination number for the ammonium cation is eight, and two types of interactions of this cation with the Lewis base centers are possible: N?H ??? N hydrogen bonds and H?N ??? N interactions that are classified as σ‐hole bonds. Redistribution of the electronic charge resulting from complexation of the ammonium cation was analyzed. On the one hand, the interactions are similar, as they lead to electronic charge transfer from the Lewis base (HCN or N₂ in this study) to NH₄⁺. On the other hand, the hydrogen bond results in the accumulation of electronic charge on the N atom of the NH₄⁺ ion, whereas the σ‐hole bond results in the depletion of the electronic charge on this atom. Quantum theory of “atoms in molecules” and the natural bond orbital method were applied to deepen the understanding of the nature of the interactions analyzed. Density functional theory/natural energy decomposition analysis was used to analyze the interactions of the ammonium ion with various types of Lewis bases. Different correlations between the geometrical, energetic, and topological parameters were found and discussed.     

ab-Initio Study of Hydrogen Bond Networks in 1,2,3-Triazole Phases

The research in storage and conversion of energy is an everlasting process. The use of fuel cells is very tempting but up to now there are still several conceptual challenges to overcome. Especially, the requirement of liquid water causes difficulties due to the temperature limit. Therefore, imidazoles and triazoles are increasingly investigated in a manifold of experimental and theoretical publications as they are both very promising in overcoming this problem. Recently, triazoles were found to be superior to imidazoles in proton conduction. An ab-initio molecular dynamics simulation of pure triazole phases for investigating the behavior of both tautomer species of the triazole molecule has never been done. In this work, we investigate the structural and dynamical properties of two different solid phases and the liquid phase at two different temperatures. We are able to show how the distinct tautomers contribute to the mechanism of proton conduction, to compute dynamical properties of the four systems and to suggest a mechanism of reorientation in solid phase.