Symmetry of NHN hydrogen bonds in solution☆

Symmetry of Hydrogen Bonds: Application of NMR Method of Isotopic Perturbation and Relevance of Solvatomers

Short, strong, symmetric, low-barrier hydrogen bonds (H-bonds) are thought to be of special significance. We have been searching for symmetric H-bonds by using the NMR technique of isotopic perturbation. Various dicarboxylate monoanions, aldehyde enols, diamines, enamines, acid-base complexes, and two sterically encumbered enols have been investigated. Among all of these, we have found only one example of a symmetric H-bond, in nitromalonamide enol, and all of the others are equilibrating mixtures of tautomers. The nearly universal lack of symmetry is attributed to the presence of these H-bonded species as a mixture of solvatomers, meaning isomers (or stereoisomers or tautomers) that differ in their solvation environment. The disorder of solvation renders the two donor atoms instantaneously inequivalent, whereupon the hydrogen attaches to the less well solvated donor. We therefore conclude that there is no special significance to short, strong, symmetric, low-barrier H-bonds. Moreover, they have no heightened stability or else they would have been more prevalent.

Theory of Borazine-Derived Nanothreads: Enumeration, Reaction Pathways, and Piezoelectricity

Nanothreads are one-dimensional macromolecules formed by pressure-induced polymerization along stacks of multiply unsaturated (or highly strained) molecules such as benzene (or cubane). Borazine is isoelectronic to benzene yet with substantial bond polarity, thus motivating a theoretical examination of borazine-derived nanothreads with degrees of saturation of 2, 4, and 6 (defined as the number of four-coordinated boron and nitrogen atoms per borazine formula unit). The energy increases upon going from molecular borazine to degree-2 borazine-derived threads and then decreases for degree-4 and degree-6 nanothreads as more σ bonds are formed. With the constraint of no more than two borazine formula units within the repeat unit of the framework's bonding topology, there are only 13 fully saturated (i.e., degree-6) borazine-derived nanothreads that avoid energetically costly homopolar bonds (as compared to more than 50 such candidates for benzene-derived threads). Only two of these are more stable than borazine. Hypothetical pathways from molecular borazine to these two degree-6 borazine-derived nanothreads are discussed. This relative paucity of outcomes may assist in kinetic control of reaction products. Beyond the high mechanical strength also predicted for carbon-based threads, properties such as piezoelectricity and flexoelectricity may be accessible to the polar lattice of borazine-derived nanothreads, with intriguing prospects for expression in these extremely thin yet rigid objects.

N-Methylated 1,8-Diaminonaphthalenes as Bifunctional Nucleo­philes in Reactions with α,ω-Dihalogenoalkanes: A Facile Route to Heterocyclic and Double Proton Sponges

The reaction of 1-dimethylamino-8-(methylamino)naphthalene with 1,3-dibromopropane chemoselectively leads to the product of N,N′-heterocyclization, while in the case of 1,4-dibromobutane and 1,2-bis(bromomethyl)benzene the process results in heterocyclization onto the same nitrogen atom with the formation of previously unknown 1-dimethylamino-8-pyrrolidino- and 1-dimethylamino-8-isoindolino-naphthalenes. The same reactions conducted without adding any auxiliary base lead to the formation of N,N′-linked double proton sponges as a new type of polynitrogen organic receptor. Proceeding as a sequence of quaternization–demethylation–cyclization steps, this heterocyclization process can also be used to construct six-membered rings (piperidino, morpholino), albeit in lower yields. The ability of 1,2-dibromoethane to brominate N-alkylated 1,8-diaminonaphthalenes is also described. It is shown for the first time that a commercially available 1,8-bis(dimethylamino)naphthalene (DMAN) can be used as a starting material in a heterocyclization reaction, which via a one-pot approach and in a short time can be converted into 1,5-dimethylnaphtho[1,8-bc]-1,5-diazacyclooctane or 1-dimethylamino-8-(pyrrolidin-1-yl)naphthalene.

Neutral Pyrrole Nitrogen Atom as a π- and Mixed n,π-Donor in Hydrogen Bonding

9-Dimethylaminobenzo[ g]indoles 3-6 and 1-dimethylamino-8-(pyrrolyl-1)naphthalene 7 were examined as possible models for establishing the ability of the pyrrole nitrogen atom to participate in [NHN]⁺ hydrogen bonding as a proton acceptor. Indoles 3-5 (to a lesser extent 6) form rather stable tetrafluoroborates, with the proton mostly located on the NMe₂ group but simultaneously engaged in the formation of a charged intramolecular [NHN]⁺ hydrogen bond (IHB) with the pyrrole N atom. The theoretically estimated energies of IHB in salts 3H⁺BF₄^--6H⁺BF₄^- vary between 7.0-10.7 and 6.2-7.0 kcal mol^-1 in vapor and MeCN, respectively. The pyrrole N atom undergoes a perceptible pyramidalization but still remains involved in the 6π-electron aromatic system, suggesting that the hydrogen bonding in salts 3H⁺BF₄^--6H⁺BF₄^- represents a previously unknown mixed NH···N(n,π) interaction. Despite the favorable orientation of the N-H bond and the pyrrole ring in salt 7H⁺BF₄^-, no signs of NH···N(n) bonding in it were noticed, and the existing interaction was classified as pure NH···N(π). The results obtained may be useful in studies of secondary protein structures, especially those α-helix sections which contain tryptophan residues.

Spectral Signatures of Proton-Transfer Dynamics at the Cusp of Low-Barrier Hydrogen Bonding

Despite their importance in diverse chemical and biochemical processes, low-barrier hydrogen bonds remain elusive targets to classify and interpret spectroscopically. Here the correlated nature of hydrogen bonding and proton transfer in the low-barrier regime has been probed for the ground and excited electronic states of 6-hydroxy-2-formylfulvene by acquiring jet-cooled fluorescence spectra of the parent and monodeuterated isotopologs. While excited-state profiles reveal regular vibronic patterns devoid of obvious dynamical signatures, their ground-state counterparts display a radically altered energy landscape characterized by spectral bifurcations comparable in magnitude to typical vibrational spacings (>100 cm^-1). Quantitative analyses yield unusual deuterium kinetic isotope effects that straddle limiting values attributed to above-barrier vibration and below-barrier tunneling of the proton adjoining donor/acceptor sites. Our findings provide compelling experimental evidence for ultrafast hydron-migration events commensurate with the onset of low-barrier hydrogen bonding and afford a trenchant glimpse of molecular phenomena taking place at the "tipping point" between disparate dynamical regimes.

pH-Free Measurement of Relative Acidities, Including Isotope Effects

A powerful pH-free multicomponent NMR titration method can measure relative acidities, even of closely related compounds, with excellent accuracy. The history of the method is presented, along with details of its implementation and a comparison with earlier NMR titrations using a pH electrode. Many of its areas of applicability are described, especially equilibrium isotope effects. The advantages of the method, some practical considerations, and potential pitfalls are considered.

Symmetry and dynamics of FHF− anion in vacuum, in CD2Cl2 and in CCl4. Ab initio MD study of fluctuating solvent–solute hydrogen and halogen bonds

Asymmetric solvation of FHF⁻ by halogen- and hydrogen-bonding solvents breaks the symmetry of the anion.

Hydrogen-bond symmetry in difluoromaleate monoanion

The symmetry of the hydrogen bond in hydrogen difluoromaleate monoanion is probed by X-ray crystallography and by the NMR method of isotopic perturbation in water, in two aprotic organic solvents, and in an isotropic liquid crystal. The X-ray crystal structure of potassium hydrogen difluoromaleate shows a remarkably short O-O distance of 2.41 Å and equal O-H distances of 1.206 Å, consistent with a strong and symmetric hydrogen bond. Incorporation of (18)O into one carboxyl group allows investigation of the symmetry of the H-bond in solution by the method of isotopic perturbation. The (19)F NMR spectra of the mono-(18)O-substituted monoanion in water, CD(2)Cl(2), and CD(3)CN show an AB spin system, corresponding to fluorines in different environments. The difference is attributed to the perturbation of the acidity of a carboxylic acid by (18)O, not to the mere presence of the (18)O, because the mono-(18)O dianion shows equivalent fluorines. Therefore, it is concluded that the monoanion exists as an equilibrating pair of interconverting tautomers and not as a single symmetric structure not only in water but also in organic solvents. However, in the isotropic liquid crystal phase of 4-cyanophenyl 4-heptylbenzoate, tetrabutylammonium hydrogen difluoromaleate-(18)O shows equivalent fluorines, consistent with a single symmetric structure. These results support earlier studies, which suggested that the symmetry of hydrogen bonds can be determined by the local environment.

On the Intramolecular Hydrogen Bond in Solution: Car-Parrinello and Path Integral Molecular Dynamics Perspective

The issue of the symmetry of short, low-barrier hydrogen bonds in solution is addressed here with advanced ab initio simulations of a hydrogen maleate anion in different environments, starting with the isolated anion, going through two crystal structures (sodium and potassium salts), then to an aqueous solution, and finally in the presence of counterions. By Car-Parrinello and path integral molecular dynamics simulations, it is demonstrated that the position of the proton in the intramolecular hydrogen bond of an aqueous hydrogen maleate anion is entirely related to the solvation pattern around the oxygen atoms of the intramolecular hydrogen bond. In particular, this anion has an asymmetric hydrogen bond, with the proton always located on the oxygen atom that is less solvated, owing to the instantaneous solvation environment. Simulations of water solutions of hydrogen maleate ion with two different counterions, K(+) and Na(+), surprisingly show that the intramolecular hydrogen-bond potential in the case of the Na(+) salt is always asymmetric, regardless of the hydrogen bonds to water, whereas for the K(+) salt, the potential for H motion depends on the location of the K(+). It is proposed that repulsion by the larger and more hydrated K(+) is weaker than that by Na(+) and competitive with solvation by water.

Are short, low-barrier hydrogen bonds unusually strong?

In a symmetric hydrogen bond (H-bond), the hydrogen atom is perfectly centered between the two donor atoms. The energy diagram for hydrogen motion is thus a single-well potential, rather than the double-well potential of a more typical H-bond, in which the hydrogen is covalently bonded to one atom and H-bonded to the other. Examples of symmetric H-bonds are often found in crystal structures, and they exhibit the distinctive feature of unusually short length: for example, the O-O distance in symmetric OHO H-bonds is found to be less than 2.5 Å. In comparison, the O-O distance in a typical asymmetric H-bond, such as ROH···OR(2), ranges from about 2.7 to 3.0 Å. In this Account, we briefly review and update our use of the method of isotopic perturbation to search for a symmetric, centered, or single-well-potential H-bond in solution. Such low-barrier H-bonds are thought to be unusually strong, owing perhaps to the resonance stabilization of two identical resonance forms [A-H···B ↔ A···H-B]. This presumptive bond strength has been invoked to explain some enzyme-catalyzed reactions. Yet in solution, a wide variety of OHO, OHN, and NHN H-bonds have all been found to be asymmetric, in double-well potentials. Examples include the monoanion of (±)-2,3-di-tert-butylsuccinic acid and a protonated tetramethylnaphthalenediamine, even though these two ions are often considered prototypes of species with strong H-bonds. In fact, all of the purported examples of strong, symmetric H-bonds have been found to exist in solution as pairs of asymmetric tautomers, in contrast to their symmetry in some crystals. The asymmetry can be attributed to the disorder of the local solvation environment, which leads to an equilibrium among solvatomers (that is, isomers that differ in solvation). If the disorder of the local environment is sufficient to break symmetry, then symmetry itself is not sufficient to stabilize the H-bond, and symmetric H-bonds do not have an enhanced stability or an unusual strength. Nor are short H-bonds unusually strong. We discuss previous evidence for "short, strong, low-barrier" H-bonds and show it to be based on ambiguous comparisons. The role of such H-bonds in enzyme-catalyzed reactions is then ascribed not to any unusual strength of the H-bond itself but to relief of "strain."

Symmetry of hydrogen bonds in solution

A classic question regarding hydrogen bonds (H-bonds) concerns their symmetry. Is the hydrogen centered or is it closer to one donor and jumping between them? These possibilities correspond to single- and double-well potentials, respectively. The NMR method of isotopic perturbation can answer this question. It is illustrated with 3-hydroxy-2-phenylpropenal and then applied to dicarboxylate monoanions. The 18O-induced 13C NMR splittings signify that their intramolecular H-bonds are asymmetric and that each species is a pair of tautomers, not a single symmetric structure, even though maleate and phthalate are symmetric in crystals. The asymmetry is seen across a wide range of solvents and a wide variety of monoanions, including 2,3-di-tert-butylsuccinate and zwitterionic phthalates. Asymmetry is also seen in monoprotonated 1,8-bis(dimethylamino)naphthalenediamines, N,N'-diaryl-6-aminofulvene-2-aldimines, and 6-hydroxy-2-formylfulvene. The asymmetry is attributed to the disorder of the local environment, establishing an equilibrium between solvatomers. The broader implications of these results regarding the role of solvation in breaking symmetry are discussed. It was prudent to confirm a secondary deuterium isotope effect (IE) on amine basicity by NMR titration of a mixture of PhCH2NH2 and PhCHDNH2. The IE is of stereoelectronic origin. It is proposed that symmetric H-bonds can be observed in crystals but not in solution because a disordered environment induces asymmetry, whereas a crystal can guarantee a symmetric environment. The implications for the controversial role of low-barrier H-bonds in enzyme-catalyzed reactions are discussed.

Substituent Effects on the Basicity of 3,7-Diazabicyclo[3.3.1]nonanes

Basicity constants for a series of 3,7-diazabicyclo[3.3.1]nonane derivatives in acetonitrile with a variation over 13 orders of magnitude have been determined using a spectrophotometric titration technique. An excellent correlation between basicity and calculated proton affinities obtained at PCM-B3LYP/6-31+G(d)//B3LYP/6-31G(d) level was found. The results are discussed in terms of substituent effects and compared to (15)N NMR chemical shifts.

Hydrogen-Bond Symmetry in Zwitterionic Phthalate Anions: Symmetry Breaking by Solvation

The cationic nitrogen of zwitterion 1 is located symmetrically with respect to its intramolecular OHO hydrogen bond. Incorporation of one (18)O allows investigation of the H-bond symmetry by the NMR method of isotopic perturbation. In both CD(3)OD and CD(2)Cl(2) equilibrium isotope shifts are detected at the carboxyl and ipso carbons. Therefore, 1 exists as a pair of interconverting tautomers, not as a single symmetric structure with its hydrogen centered between the two oxygens. The H-bond is instantaneously asymmetric, and there is an equilibrium between solvatomers (isomers or stereoisomers that differ in solvation). The broader implications of this result regarding the role of the local environment ("solvation") in breaking symmetry are discussed.

Design of C2-Chiral Diamines That Are Computationally Predicted To Be a Million-fold More Basic than the Original Proton Sponges

A set of C2-chiral diamines 18-21 based on 1,6-diazacyclodecane have been identified whose conjugate acids are predicted by B3LYP/6-31G calculations to have pKa values of approximately 23-6 on the water scale (pKa = 30-33 in MeCN); they are also expected to be kinetically active, but essentially nonnucleophilic. Strain relief on protonation largely determines the basicity of these compounds, and the key to the design of stronger bases is limiting conformational freedom, especially by preventing nitrogen inversion, through the introduction of additional ring fusions. 15,16-Dimethyl-15,16-diazatricyclo[9.3.1.1(4,8)]hexadecane (20) is examined in detail and shown to exist in 10 diastereomeric forms as a result of in-/out-isomerism. The predicted pKa values for these diastereomers range over 14 log units.

Experimental and Theoretical Evidence of Basic Site Preference in Polyfunctional Superbasic Amidinazine: N1,N1-Dimethyl-N2-β-(2-pyridylethyl)formamidine

The gas-phase basicity (GB) of the flexible polyfunctional N(1),N(1)-dimethyl-N(2)-beta-(2-pyridylethyl)formamidine (1) containing two potential basic sites (the ring N-aza and the chain N-imino) is obtained from proton-transfer equilibrium constant measurements, using Fourier-transform ion-cyclotron resonance mass spectrometry. Comparison of the experimental GB obtained for 1 with those reported for model amidines and azines indicates that the chain N-imino in the amidine group is the favored site of protonation. Semiempirical (AM1) and ab initio calculations (HF, MP2, and DFT), performed for 1 and its protonated forms, confirm this interpretation. These results are in contrast to those found previously for N(1),N(1)-dimethyl-N(2)-azinylformamidines (containing the amidine function directly linked to the azinyl ring), in which the ring N-aza is the most basic site in the gas phase. The separation of the two potential basic sites in 1 by the ethylene chain interrupts the resonance conjugation between the two functions and changes their relative basicities and, thus, the preferable site of protonation. It also increases the chelation effect against the proton and the gas-phase basicity of 1 in such a magnitude that consequently 1 may be classified as a superbase (GB = 241.1 kcal mol(-)(1)). A transition state corresponding to the internal transfer of the proton (ITP) between the ring N-aza and the chain N-imino in 1 is investigated at the DFT(B3LYP)/6-31G level. The energy barrier calculated for the ITP between the two basic sites is small and vanishes when zero-point vibrational terms and thermal corrections are applied to obtain the enthalpy or Gibbs energy of activation for the proton transfer. Additional calculations at the DFT(MPW1K)/6-31G level confirm this behavior. This indicates that the quantum-chemical ITP in 1 has a single-well character. The proton is located on the N-imino site, and the H-bond is formed with the N-aza site.

β-Deuterium Isotope Effects on Amine Basicity, “Inductive” and Stereochemical

Secondary beta deuterium isotope effects on acidity constants of ammonium ions are measured using a remarkably precise NMR titration method. Deuteration is found to increase the basicity of methylamine, dimethylamine, benzylamine, and N,N-dimethylaniline. The effect is attributed to a lowered zero-point energy of a CH bond adjacent to an amine nitrogen. The method permits a determination of the stereochemical dependence of the isotope effect in a locked piperidine, and it is found that deuteration is more effective when antiperiplanar to a lone pair. The values are consistent with a cos(2) dependence on dihedral angle, with no detectable angle-independent inductive effect.