Cooperativity in Intramolecular Bifurcated Hydrogen Bonds: An Ab Initio Study

A computational study on acetaminophen drug complexed with Mn+, Fe2+, Co+, Ni2+, and Cu+ ions: structural analysis, electronic properties, and solvent effects

In the present research, the cation-π interactions in acetaminophen-M complexes (M = Mn⁺, Fe²⁺, Co⁺, Ni²⁺, and Cu⁺) are investigated using density functional theory (DFT/ωB97XD) in the gas phase and solution. The results show that the absolute values of energy are reduced in going from the gas phase to the solution. Based on the obtained data, the complexes in water are the most stable. The natural bond orbital (NBO) and the atoms in molecules (AIM) analyses are also applied to achieve more details about the nature of interactions. These results are useful for understanding the role of the drug-receptor interactions in the complexes. According to AIM outcomes, the cation-π interactions are the closed-shell and may indicate the partially covalent nature in the complexes. A comprehensive analysis is also performed on the conceptual DFT parameters of the complexes to evaluate their electronic properties. Our findings show increasing the stability and decreasing the reactivity of the complexes in the solution phase with respect to the gas phase. These interactions are ubiquitous in biological systems, and their importance in theoretical models led us to study such important interactions. The results of this study may be useful for the design and synthesis of a variety of supramolecular complexes with the desired properties.

New benzyltriethylammonium/urea deep eutectic solvent: Quantum calculation and application to hyrdoxylethylcellulose modification

In this paper, a new deep eutectic solvent (DES) has been successfully synthesized that is based on benzyltriethylammonium bromide as a hydrogen bond acceptor (HBA) and urea as a hydrogen bond donor (HBD). However, its usability in modifying cellulose derivatives, especially acylating hydroxyethylcellulose (HEC) was investigated. The chemical modification (acetylation) of HEC was carried out in BTEAB/urea DES system without any additional conventional solvent or catalyst. However, the proposed structure of acetylated HEC (HECA) was confirmed according to the structural spectra analyses FTIR-ATR, ¹H, ¹³C, and APT-NMR. The crystalline behavior of acetylated and unmodified HEC in the DES system has been evaluated using XRD patterns, where the thermal stability was evaluated basing on the TD-TGA thermograms. Hence, SEM images and EDX spectra were recorded to prove the changes that are expected at the morphological level and elemental profile. Yet, the nanometric sheets aspect was observed. The Functional Density Theory (DFT) was investigated as a useful computational tool to understand mechanism and donor-acceptor interactions. The topological parameters (electron density Laplacian, kinetic energy density, potential energy density, and energy density) at the bond critical points (BCP), between TBEAB and urea, are deducted according to Quantum Bader's theory, and Atoms-in-molecules (AIM). The non-covalent interactions and steric effect in the DES system were studied using the reduced density gradient isosurface (RDG). Theoretical and computational calculations revealed that the H-bonds and the electrostatic coexist, as predominant interactions in the BTEAB-based DES resulting chemical structure, and mechanism formation. The physical interactions between the component entities of DES lead to a new equilibrium that is more stable than that of HBA and HBD in their separate states.

Folded network and structural transition in molten tin

The fundamental relationships between the structure and properties of liquids are far from being well understood. For instance, the structural origins of many liquid anomalies still remain unclear, but liquid-liquid transitions (LLT) are believed to hold a key. However, experimental demonstrations of LLTs have been rather challenging. Here, we report experimental and theoretical evidence of a second-order-like LLT in molten tin, one which favors a percolating covalent bond network at high temperatures. The observed structural transition originates from the fluctuating metallic/covalent behavior of atomic bonding, and consequently a new paradigm of liquid structure emerges. The liquid structure, described in the form of a folded network, bridges two well-established structural models for disordered systems, i.e., the random packing of hard-spheres and a continuous random network, offering a large structural midground for liquids and glasses. Our findings provide an unparalleled physical picture of the atomic arrangement for a plethora of liquids, shedding light on the thermodynamic and dynamic anomalies of liquids but also entailing far-reaching implications for studying liquid polyamorphism and dynamical transitions in liquids.

Unraveling the structural origin of liquid anomalies remains a challenging topic. Xu et al. propose a folded-network structural model for molten tin and provide insights into the observed second-order-like structural transition.

Hybrid Molecular Dynamics for Elucidating Cooperativity Between Halogen Bond and Water Molecules During the Interaction of p53-Y220C and the PhiKan5196 Complex

The cooperativity between hydrogen and halogen bonds plays an important role in rational drug design. However, mimicking the dynamic cooperation between these bonds is a challenging issue, which has impeded the development of the halogen bond force field. In this study, the Y220C–PhiKan5196 complex of p53 protein was adopted as a model, and the functions of three water molecules that formed hydrogen bonds with halogen atoms were analyzed by the simulation method governed by the hybrid quantum mechanical/molecular mechanical molecular dynamics. A comparison with the water-free model revealed that the strength of the halogen bond in the complex was consistently stronger. This confirmed that the water molecules formed weak hydrogen bonds with the halogen atom and cooperated with the halogen atom to enhance the halogen bond. Further, it was discovered that the roles of the three water molecules were not the same. Therefore, the results obtained herein can facilitate a rational drug design. Further, this work emphasizes on the fact that, in addition to protein pockets and ligands, the role of voids should also be considered with regard to the water molecules surrounding them.

Energy of Intramolecular Hydrogen Bonding in ortho-Hydroxybenzaldehydes, Phenones and Quinones. Transfer of Aromaticity from ipso-Benzene Ring to the Enol System(s)

Intramolecular hydrogen bonding (HB) is one of the most studied noncovalent interactions of molecules. Many physical, spectral, and topological properties of compounds are under the influence of HB, and there are many parameters used to notice and to describe these changes. Hitherto, no general method of measurement of the energy of intramolecular hydrogen bond (EHB) has been put into effect. We propose the molecular tailoring approach (MTA) for EHB calculation, modified to apply it to Ar-O-H∙∙∙O=C systems. The method, based on quantum calculations, was checked earlier for hydroxycarbonyl-saturated compounds, and for structures with resonance-assisted hydrogen bonding (RAHB). For phenolic compounds, the accuracy, repeatability, and applicability of the method is now confirmed for nearly 140 structures. For each structure its aromaticity HOMA indices were calculated for the central (ipso) ring and for the quasiaromatic rings given by intramolecular HB. The comparison of calculated HB energies and values of estimated aromaticity indices allowed us to observe, in some substituted phenols and quinones, the phenomenon of transfer of aromaticity from the ipso-ring to the H-bonded ring via the effect of electron delocalization.

What causes the weakest host to act as the strongest one? A theoretical study on the host-guest chemistry of five azacryptands and fluoride anions

In this work we have attempted to computationally analyze the important parameters affecting the selectivity in host-guest systems in order to show that the solvent effect in some host-guest systems is even more important than the hole-size fitting and/or host-guest interaction energy. For this purpose, the fluoride anion selectivities of the five most studied azacryptands with different affinities, moieties, cavity sizes and degrees of preorganization, in their fully protonated forms, were studied at PBE/TZVP and B3P86/TZVP levels of theory. The factors affecting the selectivity such as hole size fitting, steric effects, electronic properties, preorganization of hosts, as well as solvent effects were investigated. The results showed that among the five studied azacryptands, one of them, which ranks fifth for selectivity according to its weak interaction energy with fluoride anions, ranks fourth due to its good preorganization. However, a surprising improvement occurs when its selectivity for fluoride anions ranks first due to the existence of less solvent hindrance. The results of the calculations were further confirmed by a good correlation between the calculated and experimental formation constants.

Solid state structure and solution thermodynamics of three-centered hydrogen bonds (O∙∙∙H∙∙∙O) using N-(2-benzoyl-phenyl) oxalyl derivatives as model compounds

Intramolecular hydrogen bond (HB) formation was analyzed in the model compounds N-(2-benzoylphenyl)acetamide, N-(2-benzoylphenyl)oxalamate and N¹,N²-bis(2-benzoylphenyl)oxalamide. The formation of three-center hydrogen bonds in oxalyl derivatives was demonstrated in the solid state by the X-ray diffraction analysis of the geometric parameters associated with the molecular structures. The solvent effect on the chemical shift of H6 [δH6(DMSO-d₆)–δH6(CDCl₃)] and Δδ(ΝΗ)/ΔT measurements, in DMSO-d₆ as solvent, have been used to establish the energetics associated with intramolecular hydrogen bonding. Two center intramolecular HB is not allowed in N-(2-benzoylphenyl)acetamide either in the solid state or in DMSO-d₆ solution because of the unfavorable steric effects of the o-benzoyl group. The estimated ΔHº and ΔSº values for the hydrogen bonding disruption by DMSO-d₆ of 28.3(0.1) kJ·mol⁻¹ and 69.1(0.4) J·mol⁻¹·K⁻¹ for oxalamide, are in agreement with intramolecular three-center hydrogen bonding in solution. In the solid, the benzoyl group contributes to develop 1-D and 2-D crystal networks, through C–H∙∙∙A (A = O, π) and dipolar C=O∙∙∙A (A = CO, π) interactions, in oxalyl derivatives. To the best of our knowledge, this is the first example where three-center hydrogen bond is claimed to overcome steric constraints.

Ion-pairing of phosphonium salts in solution: C-H⋅⋅⋅halogen and C-H⋅⋅⋅π hydrogen bonds

The (1) H NMR chemical shifts of the C(α)H protons of arylmethyl triphenylphosphonium ions in CD2 Cl2 solution strongly depend on the counteranions X(-) . The values for the benzhydryl derivatives Ph2 CHPPh3 (+)  X(-) , for example, range from δH =8.25 (X(-) =Cl(-) ) over 6.23 (X(-) =BF4 (-) ) to 5.72 ppm (X(-) =BPh4 (-) ). Similar, albeit weaker, counterion-induced shifts are observed for the ortho-protons of all aryl groups. Concentration-dependent NMR studies show that the large shifts result from the deshielding of the protons by the anions, which decreases in the order Cl(-) > Br(-) ≫ BF4 (-) > SbF6 (-) . For the less bulky derivatives PhCH2 PPh3 (+)  X(-) , we also find CH⋅⋅⋅Ph interactions between C(α)H and a phenyl group of the BPh4 (-) anion, which result in upfield NMR chemical shifts of the C(α)H protons. These interactions could also be observed in crystals of (p-CF3 -C6 H4 )CH2 PPh3 (+)  BPh4 (-) . However, the dominant effects causing the counterion-induced shifts in the NMR spectra are the CH⋅⋅⋅X(-) hydrogen bonds between the phosphonium ion and anions, in particular Cl(-) or Br(-) . This observation contradicts earlier interpretations which assigned these shifts predominantly to the ring current of the BPh4 (-) anions. The concentration dependence of the (1) H NMR chemical shifts allowed us to determine the dissociation constants of the phosphonium salts in CD2 Cl2 solution. The cation-anion interactions increase with the acidity of the C(α)H protons and the basicity of the anion. The existence of CH⋅⋅⋅X(-) hydrogen bonds between the cations and anions is confirmed by quantum chemical calculations of the ion pair structures, as well as by X-ray analyses of the crystals. The IR spectra of the Cl(-) and Br(-) salts in CD2 Cl2 solution show strong red-shifts of the CH stretch bands. The CH stretch bands of the tetrafluoroborate salt PhCH2 PPh3 (+)  BF4 (-) in CD2 Cl2 , however, show a blue-shift compared to the corresponding BPh4 (-) salt.

"Additive" cooperativity of hydrogen bonds in complexes of catechol with proton acceptors in the gas phase: FTIR spectroscopy and quantum chemical calculations

Experimental study of hydrogen bond cooperativity in hetero-complexes in the gas phase was carried out by IR-spectroscopy method. Stretching vibration frequencies of O-H groups in phenol and catechol molecules as well as of their complexes with nitriles and ethers were determined in the gas phase using a specially designed cell. O-H groups experimental frequency shifts in the complexes of catechol induced by the formation of intermolecular hydrogen bonds are significantly higher than in the complexes of phenol due to the hydrogen bond cooperativity. It was shown that the cooperativity factors of hydrogen bonds in the complexes of catechol with nitriles and ethers in the gas phase are approximately the same. Quantum chemical calculations of the studied systems have been performed using density functional theory (DFT) methods. It was shown, that theoretically obtained cooperativity factors of hydrogen bonds in the complexes of catechol with proton acceptors are in good agreement with experimental values. Cooperative effects lead to a strengthening of intermolecular hydrogen bonds in the complexes of catechol on about 30%, despite the significant difference in the proton acceptor ability of the bases. The analysis within quantum theory of atoms in molecules was carried out for the explanation of this fact.

Bifurcated hydrogen bonds stabilize fibrils of poly(L-glutamic) acid

Model fibrillating homopolypeptides have been providing many insightful analogies to the clinically important phenomena of protein misfolding and amyloidogenesis. Here we show that the beta(2) structural variant of poly(l-glutamic) acid forms fibrils with an amyloid-like morphology, ability to enhance fluorescence of thioflavin T, and seeding properties. The beta(2) fibrils are formed upon heating of aqueous solutions of alpha-helical poly(l-glutamic) acid, which leads to a significant increase of pD (pH) of unbuffered samples and a concomitant precipitation of fibrils with unusual infrared traits: amide I' band being dramatically red-shifted to 1596 cm(-1), and the -COOD stretching band split into two peaks around 1730 and 1719 cm(-1). We are proposing that formation of three-center hydrogen bonds involving bifurcated peptide carbonyl acceptors (>C=O) and main chains' NH, as well as side chains' -COOH proton donors is likely to underlie the observed infrared characteristics of beta(2) fibrils. Such bonds provide additional conformational constraints in a tightly packed environment around glutamate side chains resulting in the decreased overall acidity of the polypeptide. The presence of bifurcated hydrogen bonds in amyloid fibrils may be an overlooked factor in fibrils' robustness, thermodynamic stability and the ability to propagate their own growth.