Does electron density in bond critical point reflect the formal charge distribution in H-bridges? The case of charge-assisted hydrogen bonds (CAHBs)

Phosphine oxides as NMR and IR spectroscopic probes for geometry and energy of PO···H–A hydrogen bonds

In this work we evaluate the possibility to use the NMR and IR spectral properties of P=O group to estimate the geometry and strength of hydrogen bonds which it forms...

Resonance-Assisted Hydrogen Bond-Revisiting the Original Concept in the Context of Its Criticism in the Literature

In the presented research, we address the original concept of resonance-assisted hydrogen bonding (RAHB) by means of the many-body interaction approach and electron density delocalization analysis. The investigated molecular patterns of RAHBs are open chains consisting of two to six molecules in which the intermolecular hydrogen bond stabilizes the complex. Non-RAHB counterparts are considered to be reference systems. The results show the influence of the neighbour monomers on the unsaturated chains in terms of the many-body interaction energy contribution. Exploring the relation between the energy parameters and the growing number of molecules in the chain, we give an explicit extrapolation of the interaction energy and its components in the infinite chain. Electron delocalization within chain motifs has been analysed from three different points of view: three-body delocalization between C=C-C, two-body hydrogen bond delocalization indices and also between fragments (monomers). A many-body contribution to the interaction energy as well as electron density helps to establish the assistance of resonance in the strength of hydrogen bonds upon the formation of the present molecular chains. The direct relation between interaction energy and delocalization supports the original concept, and refutes some of the criticisms of the RAHB idea.

Hydration of selenolate moiety: Ab initio investigation of properties of O-H⋯Se(-) hydrogen bonds in CH3 Se(-)⋯(H2 O)n clusters

This work is devoted to investigations of the influence of O-H···Se(-) hydrogen bonds on the electronic shells of selenolate R-Se(-) fragment (R═CH₃ ). The geometric, energetic and nuclear magnetic resonance (NMR) spectral parameters for various conformers of CH₃ Se(-)⋯(H₂ O)_n clusters with n = 0-6 were calculated at CCSD/aug-cc-pVDZ level of theory. For selenolate anion CH₃ Se(-) solvation free energy was calculated, and for water media it is equal to -71.41 kcal/mol. For O-H···Se(-) hydrogen bonds the proportionality coefficients between QTAIM parameters at (3; -1) bond critical point and the strength of an individual hydrogen bond ∆E were proposed. It was shown, that despite a relative weakness of O-H···Se(-) hydrogen bonds, the outer electronic shell of the selenium atom changes significantly upon formation of each hydrogen bond. This, in turn, cause the dramatic change of NMR parameters of selenium nuclei.

Inspiration for revival of old drugs: improving solubility and avoiding hygroscopicity of pipemidic acid by forming two pharmaceutical salts based on charge-assisted hydrogen bond recognitions

Improving solubility and avoiding hygroscopicity of pipemidic acid by forming pharmaceutical salts based on CAHBs.

Theoretical study of hydrogen bonding interactions in substituted nitroxide radicals

Interaction energy (E_int) of hydrogen bonded complexes of nitroxide radicals can be assessed in terms of the deepest minimum of molecular electrostatic potential (V_min).

Correlations of NHN hydrogen bond energy with geometry and 1H NMR chemical shift difference of NH protons for aniline complexes

In this computational work, we propose to use the NMR chemical shift difference of NH₂ protons for 1:1 complexes formed by aniline and nitrogen-containing proton acceptors for the estimation of the hydrogen bond energy and geometry (N⋯H and N⋯N distances). The proposed correlations could be applied to other aromatic amines as well, in a gas phase, a solution, or a solid state, for both inter- and intramolecular hydrogen bonds. We considered a set of 21 complexes with the NHN hydrogen bond without proton transfer, including hydrogen bonds from weak to medium strong ones (2-21 kcal/mol), with neutral or anionic bases and with sp³ and sp² hybridized nitrogen proton acceptors. For each complex apart from direct hydrogen bond energy calculation, we have tested several other ways to estimate the energy: (a) using a correlation between NH stretching band intensity and hydrogen bond energy and (b) using correlations between electron density properties at (3, -1) bond critical point (quantum theory of atoms in molecules analysis) and hydrogen bond energy. Besides for the studied type of complexes, we obtained refined linear correlations linking the local electron kinetic (G) and potential (V) energy densities with the hydrogen bond energy.

Salt Bridge in Ligand-Protein Complexes-Systematic Theoretical and Statistical Investigations

Although the salt bridge is the strongest among all known noncovalent molecular interactions, no comprehensive studies have been conducted to date to examine its role and significance in drug design. Thus, a systematic study of the salt bridge in biological systems is reported herein, with a broad analysis of publicly available data from Protein Data Bank, DrugBank, ChEMBL, and GPCRdb. The results revealed the distance and angular preferences as well as privileged molecular motifs of salt bridges in ligand-receptor complexes, which could be used to design the strongest interactions. Moreover, using quantum chemical calculations at the MP2 level, the energetic, directionality, and spatial variabilities of salt bridges were investigated using simple model systems mimicking salt bridges in a biological environment. Additionally, natural orbitals for chemical valence (NOCV) combined with the extended-transition-state (ETS) bond-energy decomposition method (ETS-NOCV) were analyzed and indicated a strong covalent contribution to the salt bridge interaction. The present results could be useful for implementation in rational drug design protocols.

Extremely Strong Halogen Bond. The Case of a Double-Charge-Assisted Halogen Bridge

The stable model of a double (±)charge-assisted halogen bridge has been built on the basis of searches of the Crystal Structure Database. The model, investigated by DFT theory, consists of quinuclidine-like cation derivatives and a set of simple anions. These charged fragments form halogen-bonded complexes of which the energy of complexation in some cases reaches 100 kcal/mol. Even for such strong interactions, the QTAIM characteristics are similar to those of the more classic, relatively weak halogen bonds. An important effect of complexation is the charge transfer measured by means of QTAIM and NBO. It can also be supposed, on the basis of detailed structural and QTAIM analysis, that the delocalization of the charge in a quinuclidine moiety occurs through space and not necessarily along formal bonds. The analysis of only partially charged and fully neutral counterparts of a double (±)charge-assisted halogen bridge shows significantly weaker bonding, being less than 10 kcal/mol.

Unravelling hydrogen bonding interactions of tryptamine–water dimer from neutral to cation

The physical origin for the three intermolecular hydrogen bonds in the neutral and cationic forms of the tryptamine–water dimer is explored.

Quantification of Conventional and Nonconventional Charge-Assisted Hydrogen Bonds in the Condensed and Gas Phases

Charge-assisted hydrogen bonds (CAHBs) play critical roles in many systems from biology through to materials. In none of these areas has the role and function of CAHBs been explored satisfactorily because of the lack of data on the energy of CAHBs in the condensed phases. We have, for the first time, quantified three types of CAHBs in both the condensed and gas phases for 1-(2'-hydroxylethyl)-3-methylimidazolium acetate ([C2OHmim][OAc]). The energy of conventional OH···[OAc](-) CAHBs is ∼10 kcal·mol(-1), whereas nonconventional C(sp2)H···[OAc](-) and C(sp3)H···[OAc](-) CAHBs are weaker by ∼5-7 kcal·mol(-1). In the gas phase, the strength of the nonconventional CAHBs is doubled, whereas the conventional CAHBs are strengthened by <20%. The influence of cooperativity effects on the ability of the [OAc](-) anion to deprotonate the imidazolium cation is evaluated. The ability to quantify CAHBs in the condensed phase on the basis of easier accessible gas-phase estimates is highlighted.

Hydrogen bonding in ionic liquids

Ionic liquids (IL) and hydrogen bonding (H-bonding) are two diverse fields for which there is a developing recognition of significant overlap. Doubly ionic H-bonds occur when a H-bond forms between a cation and anion, and are a key feature of ILs. Doubly ionic H-bonds represent a wide area of H-bonding which has yet to be fully recognised, characterised or explored. H-bonds in ILs (both protic and aprotic) are bifurcated and chelating, and unlike many molecular liquids a significant variety of distinct H-bonds are formed between different types and numbers of donor and acceptor sites within a given IL. Traditional more neutral H-bonds can also be formed in functionalised ILs, adding a further level of complexity. Ab initio computed parameters; association energies, partial charges, density descriptors as encompassed by the QTAIM methodology (ρBCP), qualitative molecular orbital theory and NBO analysis provide established and robust mechanisms for understanding and interpreting traditional neutral and ionic H-bonds. In this review the applicability and extension of these parameters to describe and quantify the doubly ionic H-bond has been explored. Estimating the H-bonding energy is difficult because at a fundamental level the H-bond and ionic interaction are coupled. The NBO and QTAIM methodologies, unlike the total energy, are local descriptors and therefore can be used to directly compare neutral, ionic and doubly ionic H-bonds. The charged nature of the ions influences the ionic characteristics of the H-bond and vice versa, in addition the close association of the ions leads to enhanced orbital overlap and covalent contributions. The charge on the ions raises the energy of the Ylp and lowers the energy of the X-H σ* NBOs resulting in greater charge transfer, strengthening the H-bond. Using this range of parameters and comparing doubly ionic H-bonds to more traditional neutral and ionic H-bonds it is clear that doubly ionic H-bonds cover the full range of weak through to very strong H-bonds.

Effect of the H-bonding on aromaticity of purine tautomers

Four tautomers of purine (1-H, 3-H, 7-H, and 9-H) and their equilibrium H-bonded complexes with F(-) and HF for acidic and basic centers, respectively, were optimized by means of the B3LYP/6-311++G(d,p) level of theory. Purine tautomer stability increases in the following series: 1-H < 3-H < 7-H < 9-H, consistent with increasing aromaticity. Furthermore, the presence of a hydrogen bond with HF does not change this order. For neutral H-bonded complexes, the strongest and the weakest intermolecular interactions occur (-14.12 and -10.49 kcal/mol) for less stable purine tautomers when the proton acceptor is located in the five- and six-membered rings, respectively. For 9-H and 7-H tautomers the order is reversed. The H-bond energy for the imidazole complex with HF amounts to -14.03 kcal/mol; hence, in the latter case, the fusion of imidazole to pyrimidine decreases its basicity. The ionic H-bonds of N(-)···HF type are stronger by ~10 kcal/mol than the neutral N···HF intermolecular interactions. The hydrogen bond N(-)···HF energies in pyrrole and imidazole are -32.28 and -30.03 kcal/mol, respectively, and are substantially stronger than those observed in purine complexes. The aromaticity of each individual ring and of the whole molecule for all tautomers in ionic complexes is very similar to that observed for the anion of purine. This is not the case for neutral complexes and purine as a reference. The N···HF bonds perturb much more the π-electron structure of five-membered rings than that of the six-membered ones. The H-bonding complexes for 7-H and 9-H tautomers are characterized by higher aromaticity and a much lower range of HOMA variability.

Nature of a hydride-halogen bond. A SAPT-, QTAIM-, and NBO-based study

The nature of a hydride-halogen bond is investigated by means of high-level quantum mechanical calculations expended with symmetry-adapted perturbation theory (SAPT), quantum theory of atoms in molecules (QTAIM), and natural bond orbital (NBO) methods. As model hydride-halogen bonded systems complexes between either LiH or HBeH and either XCF(3) or XCCH (X = F, Cl, Br, I) are used. It is shown that the formation of a hydride-halogen bond leads to the elongation of the R(δ+)-H(δ-) hydride bond, which is accompanied by the blue shift of the ν(R-H) stretching vibration frequency and the increase of the IR intensity of this mode. All these effects, although untypical in the case of, e.g., hydrogen bonds, can be considered as rather typical for hydride-halogen bonded systems. The decomposition of the interaction energy based on the SAPT method clearly indicates the dominant role of the induction term, thus the inductive nature of a hydride-halogen bond in opposition to previous findings. NBO-based analysis indicates the charge transfer from the hydride molecule to the more remote parts of the halogen donor and that the elongation of the R-H bond is caused by the charge outflow from the σ(RH) bonding orbital.

Long, directional interactions in cofacial silicon phthalocyanine oligomers

Single crystal structures have been determined for the three cofacial, oxygen-bridged, silicon phthalocyanine oligomers, [((CH(3))(3)SiO)(2)(CH(3))SiO](SiPcO)(2-4)[Si(CH(3))(OSi(CH(3))(3))(2)], and for the corresponding monomer. The data for the oligomers give structural parameters for a matching set of three cofacial, oxygen-bridged silicon phthalocyanine oligomers for the first time. The staggering angles between the six adjacent cofacial ring pairs in the three oligomers are not in a random distribution nor in a cluster at the intuitively expected angle of 45° but rather are in two clusters, one at an angle of 15° and the other at an angle of 41°. These two clusters lead to the conclusion that long, directional interactions (LDI) exist between the adjacent ring pairs. An understanding of these interactions is provided by atoms-in-molecules (AIM) and reduced-density-gradient (RDG) studies. A survey of the staggering angles in other single-atom-bridged, cofacial phthalocyanine oligomers provides further evidence for the existence of LDI between cofacial phthalocyanine ring pairs in single-atom-bridged phthalocyanine oligomers.