Relationships between Bond Strength and Spectroscopic Quantities in H-Bonds and Related Halogen, Chalcogen, and Pnicogen Bonds

Dynamics of non-covalent interactions during the P-O bond cleavage reaction by ribonuclease A

In this work, an atomistic-scale investigation of the phosphodiester P-O bond cleavage reaction by the enzyme ribonuclease A was carried out using computer simulation techniques. It is shown that during the reaction the network of non-covalent interactions in the active center of the ribonuclease changes significantly, while the role of these non-covalent interactions is different: coordination of the corresponding groups, electron density transfer, and ligand holding in the active center. It is shown that the process of proton transfer from Asp121 to His119 is the first stage of this reaction; at the same time, the hydrogen bond between the phosphate ligand and the imino group of Arg39 is broken, which, although keeping the ligand in the active center, does not allow the ligand to orient itself more conveniently for subsequent proton transfers. Furthermore, the key step of this reaction occurs: proton transfer with the participation of imidazole rings His12 and His119, in which the guiding role is played by several hydrogen bonds with the participation of Phe120, and the role of an electron density carrier is played by the pnictogen bond between the oxygen of the phosphate ligand and the pyridine-like nitrogen of the imidazole ring His119, which was detected for the first time.

Halogen Bonding to the π-Systems of Polycyclic Aromatics

The propensity of the π-electron system lying above a polycyclic aromatic system to engage in a halogen bond is examined by DFT calculations. Prototype Lewis acid CF3I is placed above the planes of benzene, naphthalene, anthracene, phenanthrene, naphthacene, chrysene, triphenyl, pyrene, and coronene. The I atom positions itself some 3.3-3.4 Å above the polycyclic plane, and the associated interaction energy is about 4 kcal/mol. This quantity is a little smaller for benzene, but is roughly equal for the larger polycyclics. The energy only oscillates a little as the Lewis acid slides across the face of the polycyclic, preferring regions of higher π-electron density over minima of the electrostatic potential. The binding is dominated by dispersion which contributes half of the total interaction energy.

Relation between Halogen Bond Strength and IR and NMR Spectroscopic Markers

The relationship between the strength of a halogen bond (XB) and various IR and NMR spectroscopic quantities is assessed through DFT calculations. Three different Lewis acids place a Br or I atom on a phenyl ring; each is paired with a collection of N and O bases of varying electron donor power. The weakest of the XBs display a C-X bond contraction coupled with a blue shift in the associated frequency, whereas the reverse trends occur for the stronger bonds. The best correlations with the XB interaction energy are observed with the NMR shielding of the C atom directly bonded to X and the coupling constants involving the C-X bond and the C-H/F bond that lies ortho to the X substituent, but these correlations are not accurate enough for the quantitative assessment of energy. These correlations tend to improve as the Lewis acid becomes more potent, which makes for a wider range of XB strengths.

Competition between Binding to Various Sites of Substituted Imidazoliums

The imidazolium cation has a number of different sites that can interact with a nucleophile. Adding a halogen atom (X) or a chalcogen (YH) group introduces the possibility of an NX···nuc halogen or NY···nuc chalcogen bond, which competes against the various H-bonds (NH and CH donors) as well as the lone pair···π interaction wherein the nucleophile lies above the plane of the cation. Substituted imidazoliums are paired with the NH₃ base, and the various different complexes are evaluated by density functional theory (DFT) calculations. The strength of XB and YB increases quickly along with the size and polarizability of the X/Y atom, and this sort of bond is the strongest for the heavier Br, I, Se, and Te atoms, followed by the NH···N H-bond, but this order reverses for Cl and S. The various CH···N H-bonds are comparable to one another and to the lone pair···π bond, all with interaction energies of 10-13 kcal/mol, values which show very little dependence upon the substituent placed on the imidazolium.

Exploring the Dynamical Nature of Intermolecular Hydrogen Bonds in Benzamide, Quinoline and Benzoic Acid Derivatives

The hydrogen bonds properties of 2,6-difluorobenzamide, 5-hydroxyquinoline and 4-hydroxybenzoic acid were investigated by Car-Parrinello and path integral molecular dynamics (CPMD and PIMD), respectively. The computations were carried out in vacuo and in the crystalline phase. The studied complexes possess diverse networks of intermolecular hydrogen bonds (N-H…O, O-H…N and O-H…O). The time evolution of hydrogen bridges gave a deeper insight into bonds dynamics, showing that bridged protons are mostly localized on the donor side; however, the proton transfer phenomenon was registered as well. The vibrational features associated with O-H and N-H stretching were analyzed on the basis of the Fourier transform of the atomic velocity autocorrelation function. The spectroscopic effects of hydrogen bond formation were studied. The PIMD revealed quantum effects influencing the hydrogen bridges providing more accurate free energy sampling. It was found that the N…O or O…O interatomic distances decreased (reducing the length of the hydrogen bridge), while the O-H or N-H covalent bond was elongated, which led to the increase in the proton sharing. Furthermore, Quantum Theory of Atoms in Molecules (QTAIM) was used to give insight into electronic structure parameters. Finally, Symmetry-Adapted Perturbation Theory (SAPT) was employed to estimate the energy contributions to the interaction energy of the selected dimers.

Influence of Substituents in the Benzene Ring on the Halogen Bond of Iodobenzene with Ammonia

The effects on the CI··N halogen bond between iodobenzene and NH3 of placing various substituents on the phenyl ring are monitored by quantum calculations. Substituents R = N(CH3)2, NH2, CH3, OCH3, COCH3, Cl, F, COH, CN, and NO2 were each placed ortho, meta, and para to the I. The depth of the σ-hole on I is deepened as R became more electron-withdrawing which is reflected in a strengthening of the halogen bond, which varied between 3.3 and 5.5 kcal/mol. In most cases, the ortho placement yields the largest perturbation, followed by meta and then para, but this trend is not universal. Parallel to these substituent effects is a progressive lengthening of the covalent C-I bond. Formation of the halogen bond reduces the NMR chemical shielding of all three nuclei directly involved in the C-I··N interaction. The deshielding of the electron donor N is most closely correlated with the strength of the bond, as is the coupling constant between I and N, so both have potential use as spectroscopic measures of halogen bond strength.

Dissection of the Origin of π-Holes and the Noncovalent Bonds in Which They Engage

Accompanying the rapidly growing list of σ-hole bonds has come the acknowledgment of parallel sorts of noncovalent bonds which owe their stability in large part to a deficiency of electron density in the area above the molecular plane, known as a π-hole. The origins of these π-holes are probed for a wide series of molecules, comprising halogen, chalcogen, pnicogen, tetrel, aerogen, and spodium bonds. Much like in the case of their σ-hole counterparts, formation of the internal covalent π-bond in the Lewis acid molecule pulls density toward the bond midpoint and away from its extremities. This depletion of density above the central atom is amplified by an electron-withdrawing substituent. At the same time, the amplitude of the π*-orbital is enhanced in the region of the density-depleted π-hole, facilitating a better overlap with the nucleophile's lone pair orbital and a stabilizing n → π* charge transfer. The presence of lone pairs on the central atom acts to attenuate the π-hole and shift its position somewhat, resulting in an overall weakening of the π-hole bond. There is a tendency for π-hole bonds to include a higher fraction of induction energy than σ-bonds with proportionately smaller electrostatic and dispersion components, but this distinction is less a product of the σ- or π-character and more a function of the overall bond strength.

Spectroscopic Investigation of Chalcogen Bonding: Halide-Carbon Disulfide Complexes

A combined experimental and theoretical approach has been used to study intermolecular chalcogen bonding. Specifically, the chalcogen bonding occurring between halide anions and CS₂ molecules has been investigated using both anion photoelectron spectroscopy and high-level CCSD(T) calculations. The relative strength of the chalcogen bond has been determined computationally using the complex dissociation energies as well as experimentally using the electron stabilisation energies. The anion complexes featured dissociation energies on the order of 47 kJ/mol to 37 kJ/mol, decreasing with increasing halide size. Additionally, the corresponding neutral complexes have been examined computationally, and show three loosely-bound structural motifs and a molecular radical.

Chalcogen bonds: Hierarchical ab initio benchmark and density functional theory performance study

We have performed a hierarchical ab initio benchmark and DFT performance study of D2 Ch•••A- chalcogen bonds (Ch = S, Se; D, A = F, Cl). The ab initio benchmark study is based on a series of ZORA-relativistic quantum chemical methods [HF, MP2, CCSD, CCSD(T)], and all-electron relativistically contracted variants of Karlsruhe basis sets (ZORA-def2-SVP, ZORA-def2-TZVPP, ZORA-def2-QZVPP) with and without diffuse functions. The highest-level ZORA-CCSD(T)/ma-ZORA-def2-QZVPP counterpoise-corrected complexation energies (ΔECPC ) are converged within 1.1-3.4 kcal mol-1 and 1.5-3.1 kcal mol-1 with respect to the method and basis set, respectively. Next, we used the ZORA-CCSD(T)/ma-ZORA-def2-QZVPP (ΔECPC ) as reference data for analyzing the performance of 13 different ZORA-relativistic DFT approaches in combination with the Slater-type QZ4P basis set. We find that the three-best performing functionals are M06-2X, B3LYP, and M06, with mean absolute errors (MAE) of 4.1, 4.2, and 4.3 kcal mol-1 , respectively. The MAE for BLYP-D3(BJ) and PBE amount to 8.5 and 9.3 kcal mol-1 , respectively.

Participation of S and Se in hydrogen and chalcogen bonds

The heavier chalcogen atoms S, Se, and Te can each participate in a range of different noncovalent interactions. They can serve as both proton donor and acceptor in H-bonds. Each atom can also act as electron acceptor in a chalcogen bond.

Origins and properties of the tetrel bond

The tetrel bond (TB) recruits an element drawn from the C, Si, Ge, Sn, Pb family as electron acceptor in an interaction with a partner Lewis base. The underlying principles that explain this attractive interaction are described in terms of occupied and vacant orbitals, total electron density, and electrostatic potential. These principles facilitate a delineation of the factors that feed into a strong TB. The geometric deformation that occurs within the tetrel-bearing Lewis acid monomer is a particularly important issue, with both primary and secondary effects. As a first-row atom of low polarizability, C is a reluctant participant in TBs, but its preponderance in organic and biochemistry make it extremely important that its potential in this regard be thoroughly understood. The IR and NMR manifestations of tetrel bonding are explored as spectroscopy offers a bridge to experimental examination of this phenomenon. In addition to the most common σ-hole type TBs, discussion is provided of π-hole interactions which are a result of a common alternate covalent bonding pattern of tetrel atoms.