Aromaticity-controlled tautomerism and resonance-assisted hydrogen bonding in heterocyclic enaminone-iminoenol systems

On the Relationship between Hydrogen Bond Strength and the Formation Energy in Resonance-Assisted Hydrogen Bonds

Resonance-assisted hydrogen bonds (RAHB) are intramolecular contacts that are characterised by being particularly energetic. This fact is often attributed to the delocalisation of π electrons in the system. In the present article, we assess this thesis via the examination of the effect of electron-withdrawing and electron-donating groups, namely -F, -Cl, -Br, -CF3, -N(CH3)2, -OCH3, -NHCOCH3 on the strength of the RAHB in malondialdehyde by using the Quantum Theory of Atoms in Molecules (QTAIM) and the Interacting Quantum Atoms (IQA) analyses. We show that the influence of the investigated substituents on the strength of the investigated RAHBs depends largely on its position within the π skeleton. We also examine the relationship between the formation energy of the RAHB and the hydrogen bond interaction energy as defined by the IQA method of wave function analysis. We demonstrate that these substituents can have different effects on the formation and interaction energies, casting doubts regarding the use of different parameters as indicators of the RAHB formation energies. Finally, we also demonstrate how the energy density can offer an estimation of the IQA interaction energy, and therefore of the HB strength, at a reduced computational cost for these important interactions. We expected that the results reported herein will provide a valuable understanding in the assessment of the energetics of RAHB and other intramolecular interactions.

Resonance-assisted/impaired anion-π interaction: towards the design of novel anion receptors

Substituents alter the electron density distribution in benzene in various ways, depending on their electron withdrawing and donating capabilities, as summarized by the empirical Hammett equation. The change of the π electron density distribution subsequently impacts the interaction of substituted benzenes or other cyclic conjugated rings with anions. Currently the design and synthesis of conjugated cyclic receptors capable of binding anions is an active field due to their applications in the sensing and removal of environmental contaminants and molecular recognition. By using the block-localized wavefunction (BLW) method, which is a variant of ab initio valence bond (VB) theory and can derive the reference resonance-free state self-consistently, we quantified the resonance-assisted (RA) or resonance-impaired (RI) phenomena in anion–π interactions from both structural and energetic perspectives. The frozen interaction, in which the electrostatic attraction is involved, has been shown to be the governing factor for the RA or RI interactions with anions. Energy analyses based on the empirical point charge (EPC) model indicated that the anion–π interactions can be simplified as the attraction between a negative point charge (anion) and a group of local dipoles, affected by the enriched or diminished π-cloud due to the resonance between the substituents and the conjugated ring. Hence, two strategies for the design of novel anion receptors can be envisioned. One is the enhancement of the magnitudes and/or numbers of local dipoles (polarized σ bonds), and the other is the reduction of π electron density in conjugated rings. For cases with the RI characteristics, “curved” aromatic molecules are preferred to be anion receptors. Indeed, extremely strong binding was found in complexes formed with fluorinated corannulene (F-CDD) and fluorinated [5]cycloparaphenylene (F-[5]CPP). Inspired by the RA phenomenon, complexes of p-, o- and m-benzoquinones with halides were revisited.

Substituents alter the electron density distribution in benzene in various ways, depending on their electron withdrawing and donating capabilities, as summarized by the empirical Hammett equation.

Non-conventional Hydrogen Bonding and Aromaticity: A Systematic Study on Model Nucleobases and Their Solvated Clusters

The conceptual development of aromaticity is essential to rationalize and understand the structure and behavior of aromatic heterocycles. This work addresses for the first time, the interconnection between aromaticity and sulfur/selenium centered hydrogen bonds (S/SeCHBs) involved in representative heterocycle models of canonical nucleobases (2-Pyridone; 2PY) and its sulfur (2-Thiopyridone; 2TPY) and selenium (2-Selenopyridone; 2SePY) analogs. The nucleus-independent chemical shift (NICS) and gauge induced magnetic current density (GIMIC) values suggested significant reduction of aromaticity upon replacement of exocyclic carbonyl oxygen with sulfur and selenium. However, we observed two-fold (57 %) and three-fold (80 %) enhancement in the aromaticity for 2TPY dimer, and 2SePY dimer, respectively which are connected through S/SeCHBs. Aromaticity enhancement was also noticed in 1 : 1 H-bonded complexes (heterodimers), micro hydrated clusters and for bulk hydration. It is expected that exocyclic S and Se incorporation into heterocycles without compromising aromatic loss would definitely reinforce to design new supramolecular building blocks via S/SeCH-bonded complexes.

Using Chou’s 5-Step Rule to Evaluate the Stability of Tautomers: Susceptibility of 2-[(Phenylimino)-methyl]-cyclohexane-1,3-diones to Tautomerization Based on the Calculated Gibbs Free Energies

Gibbs free energies, based on DFT (Density Functional Theory) calculations, prove that enaminone (2-(anilinemethylidene)cyclohexane-1,3-dione) and ketamine (2-[(phenylimino)-methyl]cyclohexane-1,3-dione) are the most and least stable tautomeric forms of the studied systems, respectively. 1H and 13C NMR spectra prove that 2-(anilinemethylidene)cyclohexane-1,3-diones are the only tautomeric species present in dimethylsulfoxide solution (a very weak signal can be seen only for the p-methoxy derivatives). The zwitterionic character of these enaminones is strengthened by naphthoannulation and by the insertion of the electron-withdrawing substituent into the benzene ring (the latter weakens the intramolecular hydrogen bond in the compound). Substituent and naphtoannulation have no effect on the stability of the studied tautomers. Slight twisting of the benzene ring, with respect to the CArNC plane (seen in the crystalline state), was proven to also take place in vacuum and in solution.

Intramolecular Hydrogen Bonds in Selected Aromatic Compounds: Recent Developments

A review of intramolecular hydrogen bonding in ortho-hydroxyaryl Schiff bases, ortho-hydroxyaryl Mannich bases, dipyrrins, ortho-hydroxyaryl ketones, ortho-hydroxyaryl amides, and 4-Bora-3a,4a-diaza-s-indacene (BODIPY) dyes with tautomeric sensors as substituents is presented in this paper. Ortho-hydroxy Schiff and Mannich base derivatives are known as model molecules for analysing the properties of intramolecular hydrogen bonding. The compounds under discussion possess physicochemical features modulated by the presence of strong intramolecular hydrogen bonds. The equilibrium between intra- and inter-molecular hydrogen bonds in BODIPY is discussed. Therefore, the summary can serve as a knowledge compendium of the influence of the hydrogen bond on the molecular properties of aromatic compounds.

Enol-imino-Keto-enamine Tautomerism in a Diazepine Derivative: How Decisive Are the Intermolecular Interactions in the Equilibrium?

The strange tautomeric equilibrium behavior exhibited by a new o-hydroxyphenyl diazepine derivative when the compound is analyzed both in solution and solid state opens the structural study of the enol-imino-keto-enamine forms and the influence of the intermolecular interactions in their equilibrium. The expected enol-imino form, in which the enol is part of a phenyl system and a strong O-H···N intramolecular hydrogen bond is established, results the most stable tautomer in gas phase (theoretical calculations) and was detected by NMR spectroscopy when the compound was dissolved in aprotic solvents. Nevertheless, the keto-enamine form ,in which the keto group integrates a cyclohexadienone moiety and the aromaticity of the phenol is lost, was the only tautomer in the crystal lattice according to single-crystal X-ray diffraction, vibrational spectroscopy, and diffuse reflectance results. The last form was also found as the main tautomer in UV-vis and NMR spectroscopy when a protic solvent was employed.

Open-Resonance-Assisted Hydrogen Bonds and Competing Quasiaromaticity

The delocalization of electron density upon tautomerization of a proton across a conjugated bridge can alter the strength of hydrogen bonds. This effect has been dubbed resonance-assisted hydrogen bonding (RAHB) and plays a major role in the energetics of the tautomeric equilibrium. The goal of this work was to investigate the role that π-delocalization plays in the stability of RAHBs by engaging other isomerization processes. Similarly, acid-base chemistry has received little experimental attention in studies of RAHB, and we address the role that acid-base effects play in the tautomeric equilibrium. We find that π-delocalization and the disruption of adjacent aromatic rings is the dominant effect in determining the stability of a RAHB.

Deciphering Stability of Five-Membered Heterocyclic Radicals: Balancing Act Between Delocalization and Ring Strain

Computational studies on five-membered heterocycles with single heteroatom and their isomeric dehydro-borole 1a-1c, cyclopentadiene 2a-2c, pyrrole 3a-3c, furan 4b-4c, phosphole 5a-5c, and thiophene 6b-6c radicals have been carried out. Geometrical aspects through ground state electronic structures and stability aspects using bond dissociation energies (BDE) and radical stabilization energies (RSE) have been envisaged in this regard. Spin densities, electrostatic potentials (ESP), and natural bond orbital (NBO) analysis unveiled the extent of spin delocalization. The estimated nucleus-independent chemical shifts (NICS) values revealed the difference in aromaticity characteristics of radicals. Particularly the heteroatom centered radicals exhibit odd electron π-delocalized systems with a quasi-antiaromatic character. Various factors such as, the relative position of the radical center with respect to heteroatoms, resonance, ring strain and orbital interactions influence the stability that follows the order: heteroatom centered > β-centered > α-centered radicals. Among the influences of various factors, we confirmed the existence of a competition between delocalization and the ring strain, and the interplay of both decides the overall stability order.

A Direct Proof of the Resonance-Impaired Hydrogen Bond (RIHB) Concept

The concept of resonance-enhanced hydrogen bond (RAHB) has been widely accepted and applied as it highlights the positive impact of π-conjugation on intramolecular H-bonds. However, electron delocalization is directional and there is a possibility that π-resonance goes from the H-bond acceptor to the H-bond donor, leading to a negative impact on H-bonds. Here we used the block-localized wavefunction (BLW) method which is a variant of ab initio valence bond (VB) theory and able to derive strictly electron-localized structures self-consistently, to quantify the interplay between H-bond and π-resonance in the terms of geometry, energetics and spectral properties. The comparison of geometrical optimizations with and without π-resonance shows that conjugation can indeed either enhance or weaken intramolecular H-bonds. We further experimented with various substituents attached to either the H-bond acceptor and/or H-bond donor side(s) to tune the H-bonding strength in both directions.

Some Brief Notes on Theoretical and Experimental Investigations of Intramolecular Hydrogen Bonding

A review of selected literature data related to intramolecular hydrogen bonding in ortho-hydroxyaryl Schiff bases, ortho-hydroxyaryl ketones, ortho-hydroxyaryl amides, proton sponges and ortho-hydroxyaryl Mannich bases is presented. The paper reports on the application of experimental spectroscopic measurements (IR and NMR) and quantum-mechanical calculations for investigations of the proton transfer processes, the potential energy curves, tautomeric equilibrium, aromaticity etc. Finally, the equilibrium between the intra- and inter-molecular hydrogen bonds in amides is discussed.

H-Bond: Τhe Chemistry-Biology H-Bridge

H-bonding, as a non covalent stabilizing interaction of diverse nature, has a central role in the structure, function and dynamics of chemical and biological processes, pivotal to molecular recognition and eventually to drug design. Types of conventional and non conventional (H-H, dihydrogen, H- π, CH- π, anti- , proton coordination and H-S) H-bonding interactions are discussed as well as features emerging from their interplay, such as cooperativity (σ- and π-) effects and allostery. Its utility in many applications is described. Catalysis, proton and electron transfer processes in various materials or supramolecular architectures of preorganized hosts for guest binding, are front-line technology. The H-bond-related concept of proton transfer (PT) addresses energy issues or deciphering the mechanism of many natural and synthetic processes. PT is also of paramount importance in the functions of cells and is assisted by large complex proteins embedded in membranes. Both intermolecular and intramolecular PT in H-bonded systems has received attention, theoretically and experimentally, using prototype molecules. It is found in rearrangement reactions, protein functions, and enzyme reactions or across proton channels and pumps. Investigations on the competition between intra- and intermolecular H bonding are discussed. Of particular interest is the H-bond furcation, a common phenomenon in protein-ligand binding. Multiple H-bonding (H-bond furcation) is observed in supramolecular structures.

Resonance-Assisted Hydrogen Bonding as a Driving Force in Synthesis and a Synthon in the Design of Materials

Resonance-assisted hydrogen bonding (RAHB), a concept introduced by Gilli and co-workers in 1989, concerns a kind of intramolecular H-bonding strengthened by a conjugated π-system, usually in 6-, 8-, or 10-membered rings. This Review highlights the involvement of RAHB as a driving force in the synthesis of organic, coordination, and organometallic compounds, as a handy tool in the activation of covalent bonds, and in starting moieties for synthetic transformations. The unique roles of RAHB in molecular recognition and switches, E/Z isomeric resolution, racemization and epimerization of amino acids and chiral amino alcohols, solvatochromism, liquid-crystalline compounds, and in synthons for crystal engineering and polymer materials are also discussed. The Review can provide practical guidance for synthetic chemists that are interested in exploring and further developing RAHB-assisted synthesis and design of materials.

An intramolecular hydrogen bond study in some Schiff bases of fulvene: a challenge between the RAHB concept and the σ-skeleton influence

DFT, NBO and AIM analyses have been employed to investigate which one, the resonance-assisted hydrogen bond concept or the σ-skeleton of the system, has more influence on making intramolecular hydrogen bonds stronger.

Tautomerism and metal complexation of 2-acylmethyl-2-oxazolines: a combined synthetic, spectroscopic, crystallographic and theoretical treatment

A synthetic, structural and theoretical investigation into the solid-state, solution and gas phase structure(s) of six 2-acylmethyl-4,4-dimethyl-2-oxazolines is reported. Four of these materials, viz.α-[(4,5-dihydro-4,4-dimethyl-2-oxazolyl)methylene]benzenemethanol (3a), α-[(4,5-dihydro-4,4-dimethyl-2-oxazolyl)methylene]-(4-nitrobenzene)methanol (3b), 1-(4,5-dihydro-4,4-dimethyl-2-oxazolyl)-3,3-dimethyl-1-buten-2-ol (3d) and (E)-1-phenyl-2-((3aR)-3,3a,8,8a-tetrahydro-2H-indeno[1,2-d]oxazol-2-ylidene)ethanone (3f) have been characterised in the solid-state by single crystal X-ray diffraction studies. These data represent the first solid-state structural studies of this class of compounds and details the first synthesis and full characterisation of chiral derivative 3f. All four of these materials are shown to exist in the solid phase in the enamine tautomeric form (e.g., 3a is best described as 2-[4,4-dimethyl-2-oxazolidinylidene]-1-phenylethanone) and it is suggested (NMR, IR) that this isomeric form is likely also retained in solution (e.g., CDCl3) as the more stable isomer. An investigation of the relative gas phase stabilities of the three possible (i.e., the (Z)-enol, keto and enamine) isomers of all five compounds by DFT at the B3LYP/6-311G(d) level of theory confirms the latter as the most stable form. The energy differences between the enamine and keto tautomers have been calculated to be the lowest for derivative 3d. These results are compared and contrasted with the previously reported NMR studies of such compounds which have identified the keto form as being a minor (albeit solution) tautomer. Equilibrium solution tautomer distributions for 3d are found to be solvent dependent. The protonated form of 3a, isolated as the HSO4(-) salt (i.e.4a), has been further characterised in the solid state by single crystal X-ray diffraction. These data represent the first example of a protonated oxazoline to be structurally elucidated and confirms that upon protonation, the keto (oxazoline) tautomer is the energetically favoured form in the solid-state. This observation is further supported by DFT studies for the gas phase protonated forms of such materials. Further DFT (B3LYP/6-311G(d)) calculations employing the SM8 or SMD solvation models were then applied to address the observed solution isomeric distribution for 3d; these results corroborate the gas phase theoretical treatment and also yield values that predict the higher solution stability of the enamine form as observed, although they fail to account for the existence of the keto form as a minor solution state tautomer. To access the availability of an enol-form, via hypothetical de-protonation to the enolate, compound 3a was treated with hydrated Cu(NO3)2 in EtOH solution. The resulting isolated green-coloured product (5), the first metal derivative of this entire class of ligands, is best described (IR, X-ray diffraction) as a coordinated enolate complex, i.e., Cu(3a-H)2. Complex 5 crystallizes in the P21/c space group with four molecules in the unit cell. The coordination geometry around the formal Cu(2+) metal centre is determined to be highly distorted square planar in nature (τ4 = 0.442). TD-DFT is used to give a reasonable explanation for the intensity of the absorbance band observed in the visible region for solutions of 5. These latter experiments strongly suggest that the title class of compounds may have considerable potential as ligands in coordination chemistry and/or metal-mediated catalysis.

Electron-topological, energetic and π-electron delocalization analysis of ketoenamine-enolimine tautomeric equilibrium

The ketoenamine-enolimine tautometic equilibrium has been studied by the analysis of aromaticity and electron-topological parameters. The influence of substituents on the energy of the transition state and of the tautomeric forms has been investigated for different positions of chelate chain. The quantum theory of atoms in molecules method (QTAIM) has been applied to study changes in the electron-topological parameters of the molecule with respect to the tautomeric equilibrium in intramolecular hydrogen bond. Dependencies of the HOMA aromaticity index and electron density at the critical points defining aromaticity and electronic state of the chelate chain on the transition state (TS), OH and HN tautomeric forms have been obtained.

Influence of bond fixation in benzo-annulated N-salicylideneanilines and their ortho-C(=O)X derivatives (X = CH3, NH2, OCH3) on tautomeric equilibria in solution

1H, 13C, and 15N NMR spectra show that an ortho-C(=O)X group present in the molecules of N-salicylideneanthranilamide (X = NH2), methyl N-salicylideneanthranilate (X = OCH3), N-salicylidene-o-aminoacetophenone (X = CH3), and their benzo analogues have only a minor effect on the tautomeric OH/NH-equilibrium in solution. Only two of three possible tautomers were detected. Lability of the absent form was proved by theoretical calculations. Calculated energies show that the enolimino form (OH) is less stable than the enaminone (NH) form only for dibenzo-annulated N-salicylideneanilines. The population of each species in the tautomeric mixture was found to be inversely proportional to its energy. Application of the geometry-based aromaticity index HOMA shows that the effectiveness of the pi-electron delocalization in different rings in the molecule depends mostly on the position of benzo-annulation. Both the NH...O and N...HO hydrogen bonds present in the NH and OH tautomers, respectively, increase the aromaticity of the quasirings H-O-C=C-C=N and O=C-C=C-N-H and decrease the aromatic character of the fused benzene ring. These results seem to be reliable when N-salicylideneanilines studied are compared with naphthalene and their benzo-annulated derivatives, i.e., phenanthrene, anthracene, and triphenylene. An analysis of the effectiveness of pi-electron delocalization confirms that in all cases studied, the OH form is more stable. Although the HOMA values and calculated energies are not a criterion that allows determination of the dominating tautomer, both of these parameters correctly show the effect of changes in the molecular topology on tautomeric preferences.