Tuning the Competition between Hydrogen and Tetrel Bonds by a Magnesium Bond

A computational study of the complexes formed by TF₃OH (T=C, Si, Ge) with three nitrogen-containing bases NCH, NH₃, and imidazole (IM) is carried out at the MP2/aug-cc-pVTZ level. TF₃OH can participate in two different types of noncovalent interactions: a hydrogen bond (HB) involving the hydroxyl proton and a tetrel bond (TB) with the tetel atom T. The strength of the HB is largely unaffected by the identity of T while the TB is enhanced as T grows larger. The HB is preferred over the TB for most systems, with the exception of GeF₃OH with either NH₃ or IM. MgCl₂ engages in a Mg⋅⋅⋅O Magnesium bond (Mg-bond) with the TF₃OH O atom, which cooperatively enhances both the HB and TB. The HB strengthening is particularly large for the NH₃ or IM bases, and especially for CF₃OH, but is slowly reduced as the T atom grows larger. The TB enhancement, on the other hand, behaves in the opposite fashion, accelerating for the larger T atoms. As a bottom line, the Mg-bond generally reinforces and accentuates the preference for the HB or TB that is already present in the dimer. The Mg-bond is also responsible for a proton transfer in the HB configurations with NH₃ and IM.     

Noncovalent Bonds between Tetrel Atoms

The pairing of TFH₃ with a TH₂CH₃⁻ anion, where T represents tetrel atoms C, Si, Ge, Sn, Pb, results in a strong direct interaction between the two T atoms. The interaction energy is sensitive to the nature of the two T atoms but can be as large as 90 kcal/mol. The noncovalent bond strength rises quickly as the basic T atom of the anion becomes smaller, or as the Lewis acid T grows larger, although there is less sensitivity to the latter atom. The electrostatic component makes up some 55–70 % of the total attraction energy. This term is well accounted for by simple combination of the maximum and minimum values of the molecular electrostatic potential of the Lewis acid and base units, respectively. The complexation induces a rearrangement in the TFH₃ molecule from tetrahedral to trigonal pyramidal. The associated deformation energy reduces the exothermicity of the complexation reaction. Electron density shift patterns reveal a density loss on the basic T atom, along with accompanying increases on the acidic T and its attached F atom.     

The Role of Hydrogen Bonds in Interactions between [PdCl4]2− Dianions in Crystal

[PdCl₄]²⁻ dianions are oriented within a crystal in such a way that a Cl of one unit approaches the Pd of another from directly above. Quantum calculations find this interaction to be highly repulsive with a large positive interaction energy. The placement of neutral ligands in their vicinity reduces the repulsion, but the interaction remains highly endothermic. When the ligands acquire a unit positive charge, the electrostatic component and the full interaction energy become quite negative, signalling an exothermic association. Raising the charge on these counterions to +2 has little further stabilizing effect, and in fact reduces the electrostatic attraction. The ability of the counterions to promote the interaction is attributed in part to the H-bonds which they form with both dianions, acting as a sort of glue.     

Comparison between Hydrogen and Halogen Bonds in Complexes of 6-OX-Fulvene with Pnicogen and Chalcogen Electron Donors

Quantum chemical calculations are applied to complexes of 6-OX-fulvene (X=H, Cl, Br, I) with ZH₃/H₂Y (Z=N, P, As, Sb; Y=O, S, Se, Te) to study the competition between the hydrogen bond and the halogen bond. The H-bond weakens as the base atom grows in size and the associated negative electrostatic potential on the Lewis base atom diminishes. The pattern for the halogen bonds is more complicated. In most cases, the halogen bond is stronger for the heavier halogen atom, and pnicogen electron donors are more strongly bound than chalcogen. Halogen bonds to chalcogen atoms strengthen in the order O<S<Se<Te, whereas the pattern is murkier for the pnicogen donors. In terms of competition, most halogen bonds to pnicogen donors are stronger than their H-bond analogues, but there is no clear pattern with respect to chalcogen donors. O prefers a H-bond, while halogen bonds are favored by Te. For S and Se, I-bonds are strongest, followed Br, H, and Cl-bonds in that order.     

On the Stability of Interactions between Pairs of Anions – Complexes of MCl3− (M=Be,Mg, Ca,Sr, Ba) with Pyridine and CN−

The ability of the central M atom of the MCl₃⁻ anion, with M=Be, Mg, Ca, Sr, Ba, to engage in a noncovalent bond with an approaching nucleophile is gauged by ab initio methods. The N atom of pyridine forms a M⋅⋅⋅N bond with an interaction energy between 12 and 21 kcal mol⁻¹, even though the π-hole above the M atom is not necessarily positive in sign. Despite a strong Coulombic repulsion between two anions, CN⁻ is also able to approach the M atom so as to engage in a metastable complex that is higher in energy than the individual anions. The energy barrier separating this complex from its constituent anion pair is roughly 20 kcal mol⁻¹. Despite the endothermic formation reaction energy of the CN⁻⋅⋅⋅MCl₃⁻ complex, the electron topology signals a strong interaction, more so than in pyridine⋅⋅⋅MCl₃⁻ with its exothermic binding energy. The dianionic complex is held together largely on the strength of interorbital interactions, thereby overcoming a repulsive electrostatic component. The latter is partially alleviated by the pyramidalization of the MCl₃ unit which makes its π-hole more positive. The complex sinks below the separate monomers in energy when the system is immersed in an aqueous medium, with a binding energy that varies from as much as 20 kcal mol⁻¹ for Be down to 1.2 kcal mol⁻¹ for Ba.     

Anion–Anion Interactions in Aerogen-Bonded Complexes. Influence of Solvent Environment

Ab initio calculations are applied to the question as to whether a AeX₅⁻ anion (Ae = Kr, Xe) can engage in a stable complex with another anion: F⁻, Cl⁻, or CN⁻. The latter approaches the central Ae atom from above the molecular plane, along its C₅ axis. While the electrostatic repulsion between the two anions prevents their association in the gas phase, immersion of the system in a polar medium allows dimerization to proceed. The aerogen bond is a weak one, with binding energies less than 2 kcal/mol, even in highly polar aqueous solvent. The complexes are metastable in the less polar solvents THF and DMF, with dissociation opposed by a small energy barrier.     

Attraction between electrophilic caps: A counterintuitive case of noncovalent interactions

Intermolecular attractive interaction between electrophilic sites is a counterintuitive phenomenon, as the electrostatic interaction therein is repulsive and destabilizing. Here, we confirm this phenomenon in four representative complexes, using state-of-the-art quantum mechanical methods. By employing the block-localized wavefunction (BLW) method, which can turn off intermolecular charge transfer interactions, we profoundly demonstrated the significance of charge transfer interactions in these seemingly counterintuitive complexes. Indeed, after being “turned off” the intermolecular charge transfer interaction in, for example, the FNSi···BrF complex, the originally attractive intermolecular interaction turns to be repulsive. The energy decomposition approach based on the BLW method (BLW-ED) can partition the overall stability gained on the formation of intermolecular noncovalent interaction into several physically meaningful components. According to the BLW-ED analysis, the electrostatic repulsion in these counterintuitive cases is overwhelmed by the stabilizing polarization, dispersion interaction, and most importantly, the charge transfer interaction, resulting in the eventual counterintuitive overall attraction. The present study suggests that, predicting bonding sites of noncovalent interactions using only the “hole” concept may be not universally sufficient, because other significant stabilizing factors will contribute to the stability and sometimes, play even bigger roles than the electrostatic interaction and consequently govern the complex structures. © 2018 Wiley Periodicals, Inc.     

Competition between Intra and Intermolecular Pnicogen Bonds: Complexes between Naphthalene Derivatives and Neutral or Anionic Bases

The PnF₂ (Pn=P,As,Sb,Bi) on a naphthalene scaffold can engage in an internal pnicogen Pn⋅⋅⋅N bond (PnB) with an NH₂ group placed close to it on the adjoining ring. An approaching neutral NH₃ molecule can engage in a second PnB with the central Pn, which tends to weaken the intramolecular bond. The presence of the latter in turn weakens the intermolecular PnB with respect to that formed in its absence. Replacement of the external NH₃ by a CN⁻ anion causes a fundamental change in the bonding pattern, producing a fourth covalent bond with Pn, which rearranges into a trigonal bipyramidal motif. This addition disrupts the internal Pn⋅⋅⋅N pnicogen bond, recasting the PnF₂⋅⋅⋅NH₂ interaction into an NH⋅⋅⋅F H-bond.     

Anion⋅⋅⋅Anion Attraction in Complexes of MCl3− (M=Zn,Cd, Hg) with CN−

High-level ab initio calculations show that the MCl₃⁻ anions comprising Group 2B M atoms Zn, Cd, and Hg form a stable complex with the CN⁻ anion, despite the like charge of the two ions. The complexation occurs despite a negative π-hole region above the M atom of MCl₃⁻. The dimerization distorts the planar geometry of MCl₃⁻ into a pyramidal shape which reduces the negative potential above the M atom, facilitating a close approach of the two anions, with R(M⋅⋅⋅C)∼2 Å, and an overall attractive electrostatic attraction within the dimer. In the gas phase, this dimer is less stable than the pair of separated ions by some 30 kcal/mol. However, the dissociation must surmount an energy barrier of roughly 25 kcal/mol which occurs at an intermolecular distance of 4 Å. In aqueous solution, the dimerization process is exothermic and barrier-free, with a binding energy in the 11–18 kcal/mol range.     

Complexes of HArF and AuX (X = F,Cl, Br,I). Comparison of H-bonds,halogen bonds,F-shared bonds and covalent bonds

The various sorts of complexes in which HArF and AuX (X = F, Cl, Br, I) can engage are probed by MP2/aug-cc-pVTZ calculations. The most weakly bound are those containing a halogen bond (XB) of the AuX⋯FArH sort, with binding energies less than 8 kcal/mol. H-bonded dimers FArH⋯XAu are a little stronger, held together by some 12 kcal/mol. Being the most strongly bound places the F atom of HArF roughly midway between Ar and Au in an F-shaped structure, bound by some 43–54 kcal/mol. The last sort of product involves atomic rearrangements wherein the H atom migrates from Ar to Au, followed by formation of a covalent Ar–Au bond. The resulting molecular unit is stabilized by 30–40 kcal/mol relative to the original HArF and AuX reactants. The H-bonded dimers are held together by an unusually large polarization component, surpassing electrostatic attraction, while dispersion predominates for the halogen bonds. Perturbations of the geometries and stretching frequencies offer a ready means of distinguishing the different types of complexes by spectroscopic techniques.     

Enhancement of the Tetrel Bond by the Effects of Substituents,Cooperativity, and Electric Field: Transition from Noncovalent to Covalent Bond

The T⋅⋅⋅N tetrel bond (TB) formed between TX₃OH (T=C, Si, Ge; X=H, F) and the Lewis base N≡CM (M=H, Li, Na) is studied by ab initio calculations at the MP2/aug-cc-pVTZ level. Complexes involving TH₃OH contain a conventional TB with interaction energy less than 10 kcal/mol. This bond is substantially strengthened, approaching 35 kcal/mol and covalent character, when fluorosubstituted TF₃OH is combined with NCLi or NCNa. Along with this enhanced binding comes a near equalization of the TB T⋅⋅⋅N and the internal T−O bond lengths, and the associated structure acquires a trigonal bipyramidal shape, despite a high internal deformation energy. This structural transformation becomes more complete, and the TB is further strengthened upon adding an electron acceptor BeCl₂ to the Lewis acid and a base to the NCM unit. This same TB strengthening can be accomplished also by imposition of an external electric field.     

Theoretical Studies on the Hydrogen Bonds of Different Position Action Mechanisms of Thymine with Uracil

In this research, the hydrogen bonds Y···H-X(X = C, N; Y = N, O) of thymine and uracil have been theoretically studied. The results show that hydrogen bond leads to bond length elongation and stretches the frequency red-shift of N-H···Y. Meanwhile, the C-H···O bonds shorten and stretch the frequency blue-shift. They all belong to traditional hydrogen bonds. The intermolecular charge transfer caused by the intermolecular hyperconjugation ρ*(N–H) →n(Y) and intramolecular charge redistribution by intramolecular hyperconjugation ρ(C-H)→ρ*(C-N) play important roles in the formation of hydrogen bonds. According to the judgment standards proposed by Bader and Popelier, these hydrogen bonds have typical electron density topological properties. Electrostatic surface potential(ESP) is a useful physicochemical property of a molecule that provides insights into inter- and intramolecular associations, as well as the prediction of likely sites of electrophilic and nucleophilic metabolic attack.     

Tetrel Bonds between Phenyltrifluorosilane and Dimethyl Sulfoxide: Influence of Basis Sets,Substitution and Competition

Ab initio calculations have been performed for the complexes of DMSO and phenyltrifluorosilane (PTS) and its derivatives with a substituent of NH₃, OCH₃, CH₃, OH, F, CHO, CN, NO₂, and SO₃H. It is necessary to use sufficiently flexible basis sets, such as aug’-cc-pVTZ, to get reliable results for the Si···O tetrel bonds. The tetrel bond in these complexes has been characterized in views of geometries, interaction energies, orbital interactions and topological parameters. The electron-donating group in PTS weakens this interaction and the electron-withdrawing group prominently strengthens it to the point where it exceeds that of the majority of hydrogen bonds. The largest interaction energy occurs in the p-HO₃S-PhSiF₃···DMSO complex, amounting to −122 kJ/mol. The strong Si···O tetrel bond depends to a large extent on the charge transfer from the O lone pair into the empty p orbital of Si, although it has a dominant electrostatic character. For the PTS derivatives of NH₂, OH, CHO and NO₂, the hydrogen bonded complex is favorable to the tetrel bonded complex for the NH₂ and OH derivatives, while the σ-hole interaction prefers the π-hole interaction for the CHO and NO₂ derivatives.     

二元类水滑石层板组成、结构与性能的理论研究

采用晶体学理论建立二元类水滑石(LDHs)微观结构模型与静电势能模型，将层板金属离子间距、层板电荷密度、层间阴离子间距等微观结构参数定量化，并将层间阴离子的静电势能表示成层板金属离子半径和物质的量之比、插层阴离子尺寸和电荷的函数。研究结果表明：LDHs层板金属离子间距应用离子紧密堆积来估算和孔径按阴离子平面六方点阵分布来计算是可行的；调变层板金属离子种类与物质的量之比影响层间阴离子的稳定性，势能计算值与文献报道的LDHs热稳定性次序一致。所以该模型可用于预测LDHs的微观结构参数以及热稳定性，为新型层状双羟基材料的定向合成提供思路。     

结合三维静电势参数研究化合物溶解度与溶剂分子结构间的定量关系

对一系列溶剂分子进行了结构优化和静电势及其导出参数的计算,运用多元线 性回归方法对10种化合物萘、菲、蒽、联苯、苊、六氯苯、苯偶酰、噻吨-9-酮、 二苯砜和敌草隆的溶解性能与溶剂分子的结构参数进行了关联。结果表明：分子表 面最正和最负的静电势V_(s,max)和V_(s,min)、电荷分离度Π以及分子的静电相互 作用趋势量τ这四个三维静电势参数,加上分子的前线轨道能级ε_(HOMO)和ε_ (LUMO)能很好地用于表达这些化合物在不同纯溶剂中的溶解度与溶剂分子结构间的 定量关系。     

Investigation of binding features: Effects on the interaction between CYP2A6 and inhibitors

A computational investigation has been carried out on CYP2A6 and its naphthalene inhibitors to explore the crucial molecular features contributing to binding specificity. The molecular bioactive orientations were obtained by docking (FlexX) these compounds into the active site of the enzyme. And the density functional theory method was further used to optimize the molecular structures with the subsequent analysis of molecular lipophilic potential (MLP) and molecular electrostatic potential (MEP). The minimal MLPs, minimal MEPs, and the band gap energies (the energy difference between the highest occupied molecular orbital and lowest unoccupied molecular orbital) showed high correlations with the inhibition activities (pIC₅₀s), illustrating their significant roles in driving the inhibitor to adopt an appropriate bioactive conformation oriented in the active site of CYP2A6 enzyme. The differences in MLPs, MEPs, and the orbital energies have been identified as key features in determining the binding specificity of this series of compounds to CYP2A6 and the consequent inhibitory effects. In addition, the combinational use of the docking, MLP and MEP analysis is also demonstrated as a good attempt to gain an insight into the interaction between CYP2A6 and its inhibitors. © 2010 Wiley Periodicals, Inc. J Comput Chem, 2010     

Interactions between chlorinated dioxins and a positively charged molecular probe: New molecular interaction potential

Interaction with the ligand binding domain of receptors for natural chemicals present one potential mechanism for the biological effects of environmental chemicals. Evidence suggests that the electrostatic interaction between the ligand and the receptor is an important component for binding to some of the relevant receptors. The presence of charged residues near the binding site suggests that the charge distribution of the free ligand may be different from the charge distribution of the ligand as it approaches the binding domain of the protein. In this study a new type of potential is computed for a series of dibenzo-p-dioxin (dioxin) ligands. This quantum mechanically computed potential results from interaction between the ligand and a trimethyl ammonium probe at a set of grid points. This interaction potential is compared with the molecular electrostatic potential computed from the wave function of the isolated ligands. Three types of local minima are found: (1) above the oxygen; (2) above the conjugated ring; and (3) above the chlorine(s). The molecular electrostatic potential emphasizes the minima associated with the chlorine atoms and, in that potential, the minima associated with the oxygen atoms disappear with chlorination. In the new potential, the minima over the oxygen atoms are maintained even in tetrachlorodioxin. As chlorination is increased the differences between the two potentials increases. The new potential shows the influence of the π-cation interaction, which is largest when there is little substitution on the ring. The presence of the probe induces a dipole component of 1 debye perpendicular to the plane of the ligand. Local minima in the interaction potential are then used as starting structures for the determination of the most stable ligand–probe complexes. The most stable structures are obtained from the minima associated with the oxygen atoms. These structures are stabilized by a hydrogen bond formation between the probe and the oxygen and the molecule is bent by 30° about the O(SINGLE BOND)O axis. For this series of molecules, the new potential retains some of the features that determine the hydrogen bond whereas the molecular electrostatic potential does not. © 1998 John Wiley & Sons, Inc. J Comput Chem 19: 673–684, 1998     

Anionic CH⋅⋅⋅X− Hydrogen Bonds: Origin of Their Strength,Geometry, and Other Properties

CF₃H as a proton donor was paired with a variety of anions, and its properties were assessed by MP2/aug‐cc‐pVDZ calculations. The binding energy of monoanions halide, NO₃^?, formate, acetate, HSO₄^?, and H₂PO₄^? lie in the 12–17 kcal mol^?1 range, although F^? is more strongly bound, by 26 kcal mol^?1. Dianions SO₄^2? and HPO₄^2? are bound by 27 kcal mol^?1, and trianion PO₄^3? by 45 kcal mol^?1. When two O atoms are available on the anion, the CH???O^? H‐bond (HB) is usually bifurcated, although asymmetrically. The CH bond is elongated and its stretching frequency redshifted in these ionic HBs, but the shift is reduced in the bifurcated structures. Slightly more than half of the binding energy is attributed to Coulombic attraction, with smaller contributions from induction and dispersion. The amount of charge transfer from the anions to the σ*(CH) orbital correlates with many of the other indicators of bond strength, such as binding energy, CH bond stretch, CH redshift, downfield NMR spectroscopic chemical shift of the bridging proton, and density at bond critical points.