Cooperative hydrogen bonds of macromolecules. 3. A model study of the proximity effect

Experimental and theoretical evidence for the proximity effect as a basic mechanism of the hydrogen bond cooperativity was obtained in a model system. Hydrogen bond (HB) interaction between poly(4-vinylpyridine) (PVP) and selected acids as HB donors was studied using PFG NMR self-diffusion measurements, 1H NMR longitudinal relaxation and quantum-mechanical DFT calculations. Bivalent HB donors, such as glutaric (GA) and adipic acid (AA), were compared to univalent donors 4-chlorobutyric acid, 4-acetylbutyric acid and 5-chlorovaleric acid. PFG NMR established substantially larger HB equilibrium constants for AA and, in particular, GA than for univalent donors, thus indicating cooperativity of COOH groups in bivalent donors. According to the values of these constants, the fraction of the transiently bound GA and AA molecules, which are bound by two hydrogen bonds, is 0.70 and 0.63, respectively. This result, which means substantial cooperativity in particular in GA, was then independently verified by a relaxation study comparing longitudinal relaxation rates of univalent and bivalent donors. Analysis of relaxation led to the same probabilities that HB of one COOH group of a bivalent donor will be accompanied by HB of another COOH group of the same molecule, namely 0.70 for GA and 0.63 for AA. Such cooperativity must be due to the proximity effect, i.e., the lowering of the entropy demand of the next binding by the motional restriction imposed by the already existing bonds. This conclusion is in excellent agreement with DFT calculations on the interaction of GA with a model vinylpyridine dimer, which predict preference of double binding of the same GA molecule over that of two GA molecules and show that this preference is due to a substantially lower entropy demand.     

Bifurcated hydrogen-bonding effect on the shielding and coupling constants in trifluoroacetyl pyrroles as studied by 1H, 13C and 15N NMR spectroscopy and DFT calculations

According to the (1)H, (13)C and (15)N NMR spectroscopic data and DFT calculations, bifurcated N--H...N and N--H...O intramolecular hydrogen bond is shown to be present in 2-trifluoroacetyl-5-(2'-pyridyl)-pyrrole. This bifurcated hydrogen bond causes an increase in the absolute size of the (1)J(N,H) coupling constant by about 6 Hz, and the deshielding of the bridge proton by 2 ppm. DFT calculations show that the influence of the N--H...N and N--H...O intramolecular hydrogen bonds on the (1)J(N,H) coupling and proton shielding is almost additive, although the components of the bifurcated hydrogen bond slightly weaken each other. In 2-trifluoroacetyl-5-(2'-pyridyl)-pyrrole, the coupling constants involving the fluorine and the N--H covalent bond nuclei depend dramatically on the spatial position of the pyridine ring. The pyridine ring rotation operates as a quantum switch controlling the spin information transfer between the (19)F and (15)N nuclei, as well as the proton.     

Molecular and crystal structure of 1,8,9,20-tetrahydro-3,6-dimethyl-1,8-diphenyldibenzo[c,m]dipyrazolo[3.4-f:4′, 3′-j][1,2,5,8,9,12]-hexaazacyclotetradecane

An X-ray structural study has been carried out of macroheterocycles constructed from azobispyrazole and azobisbenzene moieties linked by amino nitrogen atoms. The molecule exists as a 14-membered, highly nonplanar hexaaza-heterocycle. The hydrogen atoms of the amino-groups participate in two types of intramolecular hydrogen bonds of the N-H^...H type, one in a forked bond, and the other in a single bond.Translated from Khimiya Geterotsiklicheskikh Soedinenii, No. 11, pp. 1564–1569, November, 1985.     

A Theoretical Study on Hydrogen‐Bonded Complex of Proflavine Cation and Water: The Site‐dependent Feature of Hydrogen Bond Strengthening and Weakening

The intermolecular hydrogen‐bonds between proflavine cation (PC) and water molecules are investigated by density functional theory (DFT) and time‐dependent density functional theory (TDDFT) methods. The ground‐state geometry optimizations, electronic excitation energies and corresponding oscillation strengths of the low‐lying electronically excited states for the isolated proflavine cation, the hydrogen‐bonded PC–H2O dimer and PC–(H2O)2 trimer are calculated. Intermolecular hydrogen bonds at the central site of proflavine molecule are found to be stronger than the peripheral site. The hydrogen bond N–H???O for the hydrogen‐bonded dimer are indicated to be weakened in the excited states, since the excitation energy is increased slightly comparing to the monomer. Hydrogen bonds of PC–(H2O)2 trimer with the same type as the dimer are strengthened in the excited state, which is demonstrated by the decrease of the excited energies. Thus, hydrogen bond strengthening and weakening are observed to reveal site dependent feature in proflavine molecule. Furthermore, the hydrogen bond at central site induces the blue‐shift of the absorption spectrum, while the ones at peripheral site induce red‐shift. Hydrogen bonds with the same type at peripheral and central sites of proflavine molecule provide different effects on the photochemical and photophysical properties of proflavine.     

New thermochemical parameter for describing solvent effects on IR stretching vibration frequencies. Communication 2. Assessment of cooperativity effects

Solvent effects on O-H stretching vibration frequency of methanol in hydrogen bond complexes with different bases, CH3OH...B, have been investigated by FTIR spectroscopy. Using chloroform as a solvent results in strengthening of CH3OH...B hydrogen bonding due to cooperativity between CH3OH...B and Cl3CH...CH3OH bonds. A method is proposed for quantifying the hydrogen bond cooperativity effect. The determined cooperativity factors take into account all specific interactions of the solute in proton-donor solvents. In addition, a method of estimation of cooperativity factors Ab and AOX in system (CH3OH)2...B is proposed. It is demonstrated that in such systems, the cooperativity factor of the OH...B bond decreases and that of the OH...O bond increases with increasing the acceptor strength of the base B. The obtained results are in a good agreement with the data obtained previously from matrix-isolation FTIR spectroscopy.     

Theoretical Study on the Hydrogen Bonding Interactions in Complexes of 5‐Hydroxytryptamine with Water

The energies, geometries and harmonic vibrational frequencies of 1:1 5‐hydroxytryptamine‐water (5‐HT‐H₂O) complexes are studied at the MP2/6‐311++G(d,p) level. Natural bond orbital (NBO), quantum theory of atoms in molecules (QTAIM) analyses and the localized molecular orbital energy decomposition analysis (LMO‐EDA) were performed to explore the nature of the hydrogen‐bonding interactions in these complexes. Various types of hydrogen bonds (H‐bonds) are formed in these 5‐HT‐H₂O complexes. The intermolecular C4H5^5‐HT···O^w H‐bond in HTW3 is strengthened due to the cooperativity, whereas no such cooperativity is found in the other 5‐HT‐H₂O complexes. H‐bond in which nitrogen atom of amino in 5‐HT acted as proton donors was stronger than other H‐bonds. Our researches show that the hydrogen bonding interaction plays a vital role on the relative stabilities of 5‐HT‐H₂O complexes.     

Hydrogen‐bonded structures of the isomeric 2‐, 3‐ and 4‐carbamoylpyridinium hydrogen chloranilates

In the three isomeric salts, all C₆H₇N₂O⁺·C₆HCl₂O₄⁻, of chloranilic acid (2,5‐dichloro‐3,6‐dihydroxy‐1,4‐benzoquinone) with 2‐, 3‐ and 4‐carbamoylpyridine, namely, 2‐carbamoylpyridinium hydrogen chloranilate (systematic name: 2‐carbamoylpyridinium 2,5‐dichloro‐4‐hydroxy‐3,6‐dioxocyclohexa‐1,4‐dienolate), (I), 3‐carbamoylpyridinium hydrogen chloranilate, (II), and 4‐carbamoylpyridinium hydrogen chloranilate, (III), acid–base interactions involving H‐atom transfer are observed. The shortest interactions between the cation and the anion in (I) and (II) are pyridinium N—H...(O,O) bifurcated hydrogen bonds, which act as the primary intermolecular interaction in each crystal structure. In (III), an amide N—H...(O,O) bifurcated hydrogen bond, which is much weaker than the bifurcated hydrogen bonds in (I) and (II), connects the cation and the anion.     

Theoretical study of the symmetry of the (OH...O)- hydrogen bonds in vinyl alcohol-vinyl alcoholate systems

The interactions between substituted vinyl alcohols and vinyl alcoholates (X = NH(2), H, F, Cl, CN) are studied at the B3LYP/6-311++G(d,p) level of theory. In a first step, the conformation of the monomers is investigated and the proton affinities (PA(A(-))) of the enolates are calculated. The enols and enolates are held together by strong (OH...O)(-) hydrogen bonds, the hydrogen bond energies ranging from 19.1 to 34.6 kcal mol(-1). The optimized O...O distances are between 2.414 and 2.549 A and the corresponding OH distances from 1.134 and 1.023 A. The other geometry parameters such as C[double bond]C or CO distances also indicate that, in the minimum energy configuration, the hydrogen bonds are characterized by a double well potential. The Mulliken charges on the different atoms of the proton donors and proton acceptors and the frequencies of the nu(OH) stretching vibrations agree with this statement. All the data indicate that the hydrogen bonds are the strongest in the homomolecular complexes. The transition state for hydrogen transfer is located with the transition barrier estimated to be about zero. Upon addition of the zero-point vibration energies to the total potential energy, the barrier vanishes. This is a characteristic feature of low-barrier hydrogen bonds (LBHBs). The hydrogen bond energies are correlated to the difference 1.5 PA(AH) - PA(A(-)). The correlation predicts different energies for homomolecular hydrogen bonds, in agreement with the theoretical calculations. Our results suggest that a PA (or pK(a)) match is not a necessary condition for forming LBHBs in agreement with recent data on the intramolecular hydrogen bond in the enol form of benzoylacetone (J. Am. Chem. Soc. 1998, 120, 12117).     

Blue-shifted A-H stretching modes and cooperative hydrogen bonding. 1. Complexes of substituted formaldehyde with cyclic hydrogen fluoride and water clusters

The structures and vibrational spectra of the intermolecular complexes formed by insertion of substituted formaldehyde molecules HRCO (R = H, Li, F, Cl) into cyclic hydrogen fluoride and water clusters are studied at the MP2/aug-cc-pVTZ computational level. Depending on the nature of the substituent R, the cluster type, and its size, the C-H stretching modes of HRCO undergo large blue and partly red shifts, whereas all the F-H and O-H stretching modes of the conventional hydrogen bonds are strongly red-shifted. It is shown that (i) the mechanism of blue shifting can be explained within the concept of the negative intramolecular coupling between C-H and C=O bonds that is inherent to the HRCO monomers, (ii) the blue shifts also occur even if no hydrogen bond is formed, and (iii) variation of the acceptor X or the strength of the C-H...X hydrogen bond may either amplify the blue shift or cause a transition from blue shift to red shift. These findings are illustrated by means of intra- and intermolecular scans of the potential energy surfaces. The performance of the negative intramolecular coupling between C-H and C=O bonds of H(2)CO is interpreted in terms of the NBO analysis of the isolated H(2)CO molecule and H(2)CO interacting with (H2O)n and (HF)n clusters.     

The progression of strong and weak hydrogen bonds in a series of ethylenediammonium dithiocyanate derivatives--a new bonding protocol for macromolecules?

The experimental charge densities for a series of sym-N-methyl-substituted ethylenediammonium dithiocyanate salts have been investigated based on low-temperature and high-resolution X-ray diffraction data. This series of organic dications provides both strong and weak hydrogen bonding networks that vary depending on the N-H : SCN(-) (donor/acceptor) ratios. The number of N-HN hydrogen bonds connected to each cation increases (linear to bifurcated) as the number of N-H donor groups increases. The bifurcated thiocyanate anions also form a less energetic N-HS hydrogen bond. The presence of more than one hydrogen bond acceptor on each thiocyanate anion results in a competition between the sulfur and nitrogen atoms in forming both strong and weak hydrogen bonds. The formation of a significant number of weak hydrogen bonds is shown to play a crucial role in stabilizing these organic ionic crystals. The progression of these organic dications (smaller to larger N-H : SCN(-) ratios) results in the weaker hydrogen bonds playing a smaller role in stabilizing the crystalline structures. In addition, the electron density along the saddle point has been shown to vary significantly from weak hydrogen bonds to van der Waals interactions. This has led to a better understanding of the progression of hydrogen bonding in the crystalline states of sym-N-methyl substituted ethylenediammonium dithiocyanate salts and provides insight into the relationship between strong and weak hydrogen bonds in organic ionic crystals.     

Competition between nonclassical hydrogen-bonded acceptor sites in complexes of neutral AH2 Radicals (A = B, Al, and Ga): A theoretical investigation

An ab initio computational study of the properties of the neutral AH2 radicals (A = B, Al, Ga) as hydrogen-bond (HB) acceptors, with H-X (X = F, Cl, Br, CN, and CCH) as HB donors, is carried out at the UMP2/6-311++G(2d,2p) level. Two different minima have been found for each of the 15 possible dimers. One structure corresponds to a single-electron hydrogen-bonded complex (SEHB), with the A atom acting as an HB acceptor. The second corresponds to a dihydrogen bond complex between one of the hydrogen atoms of AH2 and the H-X molecule. Thus, all the atoms of the neutral AH2 molecule can act as HB acceptors and none as donors. The stability of the SEHB complexes decreases as BH2 > AlH2 > GaH2, while for the dihydrogen-bonded complexes the order is AlH2 > GaH2 > BH2. For the BH2 radical the SEHB complexes are stronger than the dihydrogen bonded ones, while the opposite is found for the AlH2 and GaH2 systems. Regarding the HB donors, the order found for the binding energy in the two types of complexes is H2A...HF > H2A...HCl > H2A...HBr > H2A...HCN > H2A...HCCH.     

Strong Short‐Range Cooperativity in Hydrogen‐Bond Chains

Chains of hydrogen bonds such as those found in water and proteins are often presumed to be more stable than the sum of the individual H bonds. However, the energetics of cooperativity are complicated by solvent effects and the dynamics of intermolecular interactions, meaning that information on cooperativity typically is derived from theory or indirect structural data. Herein, we present direct measurements of energetic cooperativity in an experimental system in which the geometry and the number of H bonds in a chain were systematically controlled. Strikingly, we found that adding a second H‐bond donor to form a chain can almost double the strength of the terminal H bond, while further extensions have little effect. The experimental observations add weight to computations which have suggested that strong, but short‐range cooperative effects may occur in H‐bond chains.     

Quantifying weak hydrogen bonding in uracil and 4-cyano-4'-ethynylbiphenyl: a combined computational and experimental investigation of NMR chemical shifts in the solid state

Weak hydrogen bonding in uracil and 4-cyano-4'-ethynylbiphenyl, for which single-crystal diffraction structures reveal close CH...O=C and C[triple bond]CH...N[triple bond]C distances, is investigated in a study that combines the experimental determination of 1H, 13C, and 15N chemical shifts by magic-angle spinning (MAS) solid-state NMR with first-principles calculations using plane-wave basis sets. An optimized synthetic route, including the isolation and characterization of intermediates, to 4-cyano-4'-ethynylbiphenyl at natural abundance and with 13C[triple bond]13CH and 15N[triple bond]C labeling is described. The difference in chemical shifts calculated, on the one hand, for the full crystal structure and, on the other hand, for an isolated molecule depends on both intermolecular hydrogen bonding interactions and aromatic ring current effects. In this study, the two effects are separated computationally by, first, determining the difference in chemical shift between that calculated for a plane (uracil) or an isolated chain (4-cyano-4'-ethynylbiphenyl) and that calculated for an isolated molecule and by, second, calculating intraplane or intrachain nucleus-independent chemical shifts that quantify the ring current effects caused by neighboring molecules. For uracil, isolated molecule to plane changes in the 1H chemical shift of 2.0 and 2.2 ppm are determined for the CH protons involved in CH...O weak hydrogen bonding; this compares to changes of 5.1 and 5.4 ppm for the NH protons involved in conventional NH...O hydrogen bonding. A comparison of CH bond lengths for geometrically relaxed uracil molecules in the crystal structure and for geometrically relaxed isolated molecules reveals differences of no more than 0.002 A, which corresponds to changes in the calculated 1H chemical shifts of at most 0.1 ppm. For the C[triple bond]CH...N[triple bond]C weak hydrogen bonds in 4-cyano-4'-ethynylbiphenyl, the calculated molecule to chain changes are of similar magnitude but opposite sign for the donor 13C and acceptor 15N nuclei. In uracil and 4-cyano-4'-ethynylbiphenyl, the CH hydrogen-bonding donors are sp2 and sp hybridized, respectively; a comparison of the calculated changes in 1H chemical shift with those for the sp3 hybridized CH donors in maltose (Yates et al. J. Am. Chem. Soc. 2005, 127, 10216) reveals no marked dependence on hybridization for weak hydrogen-bonding strength.     

2‐{[(3‐Fluorophenyl)amino]methylidene}‐3‐oxobutanenitrile and 5‐{[(3‐fluorophenyl)amino]methylidene}‐2,2‐dimethyl‐1,3‐dioxane‐4,6‐dione: X‐ray and DFT studies

In the crystal structures of the title compounds, C₁₁H₉FN₂O, (I), and C₁₃H₁₂FNO₄, (II), the molecules are joined pairwise via different hydrogen bonds and the constituent pairs are crosslinked by weak C—H...O hydrogen bonds. The basic structural motif in (I), which is partially disordered, comprises pairs of molecules arranged in an antiparallel fashion which enables C—H...N[triple‐bond]C interactions. The pairs of molecules are crosslinked by two weak C—H...O hydrogen bonds. The constituent pair in (II) is formed by intramolecular bifurcated C—H...O/O′ and combined inter‐ and intramolecular N—H...O hydrogen bonds. In both structures, F atoms form weak C—F...H—C interactions with the H atoms of the two neighbouring methyl groups, the H...F separations being 2.59/2.80 and 2.63/2.71 Å in (I) and (II), respectively. The bond orders in the molecules, estimated using the natural bond orbitals (NBO) formalism, correlate with the changes in bond lengths. Deviations from the ideal molecular geometry are explained by the concept of non‐equivalent hybrid orbitals. The existence of possible conformers of (I) and (II) is analysed by molecular calculations at the B3LYP/6–31+G** level of theory.     

A whole zoo of hydrogen bonds in one crystal structure: tris(isonicotinium) hydrogensulfate sulfate monohydrate

Depending on the reaction partner, the organic ditopic molecule isonicotinic acid (Hina) can act either as a Brønsted acid or base. With sulfuric acid, the pyridine ring is protonated to become a pyridinium cation. Crystallization from ethanol affords the title compound tris(4‐carboxypyridinium) hydrogensulfate sulfate monohydrate, 3C₆H₆NO₂⁺·HSO₄⁻·SO₄²⁻·H₂O or [(H₂ina)₃(HSO₄)(SO₄)(H₂O)]. This solid contains 11 classical hydrogen bonds of very different flavour and nonclassical C—H…O contacts. All N—H and O—H donors find at least one acceptor within a suitable distance range, with one of the three pyridinium H atoms engaged in bifurcated N—H…O hydrogen bonds. The shortest hydrogen‐bonding O…O distance is subtended by hydrogensulfate and sulfate anions, viz. 2.4752 (19) Å, and represents one of the shortest hydrogen bonds ever reported between these residues.     

Factors controlling extremely strong AAA-DDD triply hydrogen-bonded complexes

Triply hydrogen-bonded complexes of the form AAA-DDD are shown to have the strongest interaction when the complex is substituted with electron withdrawing groups on the donor molecule (DDD) and electron donating groups on the acceptor molecule (AAA). In particular, the largest effects are observed when the withdrawing groups act through resonance. This serves to flatten the entire system resulting in more linear, and consequently stronger, hydrogen bonds. Furthermore, the present calculations show that the binding energy correlates with the electron density at the bond critical points and inversely with the hydrogen bond lengths.     

Microsolvation of formamide: a rotational study

Microsolvated formamide clusters have been generated in a supersonic jet expansion and characterized using Fourier transform microwave spectroscopy. Three conformers of the monohydrated cluster and one of the dihydrated complex have been observed. Seven monosubstituted isotopic species have been measured for the most stable conformer of formamide...H(2)O, which adopts a closed planar ring structure stabilized by two intermolecular hydrogen bonds (N-H...O(H)-H...O=C). The two higher energy forms of formamide...H(2)O have been observed for the first time. The second most stable conformer is stabilized by a O-H...O=C and a weak C-H...O hydrogen bond, while, in the less stable form, water accepts a hydrogen bond from the anti hydrogen of the amino group. For formamide...(H(2)O)(2), the parent and nine monosubstituted isotopic species have been observed. In this cluster the two water molecules close a cycle with the amide group through three intermolecular hydrogen bonds (N-H...O(H)-H...O(H)-H...O=C), the nonbonded hydrogen atoms of water adopting an up-down configuration. Substitution (r(s)) and effective (r(0)) structures have been determined for formamide, the most stable form of formamide...H(2)O and formamide...(H(2)O)(2). The results on monohydrated formamide clusters can help to explain the observed preferences of bound water in proteins. Clear evidence of sigma-bond cooperativity effects emerges when comparing the structures of the mono- and dihydrated formamide clusters. No detectable structural changes due to pi-bond cooperativity are observed on formamide upon hydration.     

Cooperativity between the Dihydrogen Bond and the N⋅⋅⋅HC Hydrogen Bond in LiH–(HCN)n Complexes

The cooperativity between the dihydrogen bond and the N???HC hydrogen bond in LiH–(HCN)_n (n=2 and 3) complexes is investigated at the MP2 level of theory. The bond lengths, dipole moments, and energies are analyzed. It is demonstrated that synergetic effects are present in the complexes. The cooperativity contribution of the dihydrogen bond is smaller than that of the N???HC hydrogen bond. The three‐body energy in systems involving different types of hydrogen bonds is larger than that in the same hydrogen‐bonded systems. NBO analyses indicate that orbital interaction, charge transfer, and bond polarization are mainly responsible for the cooperativity between the two types of hydrogen bonds.     

Weak C—H...X (X = O,N) hydrogen bonds in the crystal structure of dihydroberberine

Dihydroberberine (systematic name: 9,10‐dimethoxy‐6,8‐dihydro‐5H‐1,3‐dioxolo[4,5‐g]isoquinolino[3,2‐a]isoquinoline), C₂₀H₁₉NO₄, a reduced form of pharmacologically important berberine, crystallizes from ethanol without interstitial solvent. The molecule shows a dihedral angle of 27.94 (5)° between the two arene rings at the ends of the molecule, owing to the partial saturation of the inner quinolizine ring system. Although lacking classical O—H or N—H donors, the packing in the crystalline state is clearly governed by C—H...N and C—H...O hydrogen bonds involving the two acetal‐type C—H bonds of the 1,3‐dioxole ring. Each dihydroberberine molecule is engaged in four hydrogen bonds with neighbouring molecules, twice as donor and twice as acceptor, thus forming a two‐dimensional sheet network that lies parallel to the (100) plane.