On Integral INICS Aromaticity of Pyridodiazepine Constitutional Isomers and Tautomers

The structure, energetics, and aromaticity of c.a. 100 constitutional isomers and tautomers of pyrido[m,n]diazepines (m = 1, 2; n = 2, 3, 4, 5; m ≠ n) were studied at the B3LYP/cc-pVTZ level. The pyrido[1,3]diazepines appear the most, while pyrido[2,4]diazepines are the least stable (ca. 26 kcal/mol). In the pyrido[1,n]diazepine group (n = 2-5), the [1,5] isomers are higher in energy by ca. 4.5 kcal/mol and the [1,4] ones by ca. 7 kcal/mol, and the pyrido[1,2]diazepines are the least stable (ca. 20 kcal/mol). All the most stable pyrido[1,n]diazepines have N-atoms near the ring's junction bond but on opposite sites. The most stable [2,n]-forms are also those with the pyridine ring N6-atom near the junction bond. Surprisingly, for the [1,2]-, [1,3]-, and [1,4]-isomer condensation types of pyridine and diazepine rings, the same N9 > N7 > N6 > N8 stability pattern obeys. The stability remains similar in a water medium simulated with the Polarizable Continuum Model of the solvent and is conserved when calculated using the CAM-B3LYP or BHandHlyp functionals. The ring's aromaticity in the pyridine[m,n]diazepines was established based on the integral INICS index resulting from the NICSzz-scan curves' integration. The integral INICS index is physically justified through its relation to the ringcurrent as demonstrated by Berger, R.J.F., et al. Phys. Chem. Chem. Phys. 2022, 24, 624. The six-membered pyrido rings have negative INICSZZ indices and can be aromatic only if they are not protonated at the N-atom. All protonated pyrido and seven-membered rings exhibit meaningful positive INICSZZ values and can be assigned as antiaromatic. However, some non-protonated pyrido rings also have substantial positive INICSZZ indices and are antiaromatic. A weak linear correlation (R2 = 0.72) between the INICSZZ values of the pyridine I(6) and diazepine I(7) rings exists and is a consequence of the communication between the π-electron systems of the two rings. The juxtaposition of the INICS descriptor of the six- and seven-membered rings and diverse electron density parameters at the Ring Critical Points (RCP) revealed good correlations only with the Electrostatic Potentials from the electrons and nuclei (ESPe and ESPn). The relationships with other RCP parameters like electron density and its Laplacian, total energy, and the Hamiltonian form of kinetic energy density were split into two parts: one nearly constant for the six-membered rings and one linearly correlating for the seven-membered rings. Thus, most of the electron density parameters at the RCP of the six-membered rings of pyridodiazepines practically do not change with the diazepine type and the labile proton position. In contrast, those of the seven-membered rings display aromaticity changes in the antiaromatic diazepine with its ring structural modifications.

Three decades of unveiling the complex chemistry of C-nitroso species with computational chemistry

This review revisits the complex reactivity of C-nitroso derivatives through the synergistic combination of computational and synthetic organic chemistry, with an emphasis on the rationalization of mechanisms and selectivities. 

Application of the Extended HOMED (Harmonic Oscillator Model of Aromaticity) Index to Simple and Tautomeric Five-Membered Heteroaromatic Cycles with C, N, O, P, and S Atoms

The geometry-based HOMA (Harmonic Oscillator Model of Aromaticity) descriptor, based on the reference compounds of different delocalizations of n- and π-electrons, can be applied to molecules possessing analogous bonds, e.g., only CC, only CN, only CO, etc. For compounds with different heteroatoms and a different number of CC, CX, XX, and XY bonds, its application leads to some discrepancies. For this reason, the structural descriptor was modified and the HOMED (Harmonic Oscillator Model of Electron Delocalization) index defined. In 2010, the HOMED index was parameterized for compounds with C, N and O atoms. For parametrization, the reference molecules of similar delocalizations of n- and π-electrons were employed. In this paper, the HOMED index was extended to compounds containing the CP, CS, NN, NP, PP, NO, NS, PO, and PS bonds. For geometrical optimization of all reference molecules and of all investigated heterocompounds, the same quantum–chemical method {B3LYP/6-311+G(d,p)} was used to eliminate errors of the HOMED estimation. For some tautomeric systems, the Gn methods were also employed to confirm tautomeric preferences. The extended HOMED index was applied to five-membered heterocycles, simple furan and thiophene, and their N and P derivatives as well as for tautomeric pyrrole and phosphole and their N and P derivatives. The effects of additional heteroatom(s) in the ring on the HOMED values for furan are parallel to those for thiophene. For pyrroles, aromaticity dictates the tautomeric preferences. An additional N atom in the ring only slightly affects the HOMED values for the favored and well delocalized NH tautomers. Significant changes take place for their rare CH forms. When intramolecular proton-transfer is considered for phosphole and its P derivatives, the PH tautomers seem to be favored only for 1,2,3-triphosphole/1,2,5-triphosphole and for 1,2,3,5-tetraphosphole. For other phospholes, the CH forms have smaller Gibbs energies than the PH isomers. For phosphazoles, the labile proton in the favored form is linked to the N atom. The PH forms have smaller HOMED indices than the NH tautomers but higher than the CH ones.

NMR structural study of the prototropic equilibrium in solution of Schiff bases as model compounds

An NMR titration method has been used to simultaneously measure the acid dissociation constant (pKa) and the intramolecular NHO prototropic constant ΔKNHO on a set of Schiff bases. The model compounds were synthesized from benzylamine and substituted ortho-hydroxyaldehydes, appropriately substituted with electron-donating and electron-withdrawing groups to modulate the acidity of the intramolecular NHO hydrogen bond. The structure in solution was established by 1H-, 13C- and 15N-NMR spectroscopy. The physicochemical parameters of the intramolecular NHO hydrogen bond (pKa, ΔKNHO and ΔΔG°) were obtained from 1H-NMR titration data and pH measurements. The Henderson-Hasselbalch data analysis indicated that the systems are weakly acidic, and the predominant NHO equilibrium was established using Polster-Lachmann δ-diagram analysis and Perrin model data linearization.

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.

Stability, polarity, intramolecular interactions and π-electron delocalization for all eighteen tautomers rotamers of uracil. DFT studies in the gas phase

Complete tautomeric equilibria were investigated for uracil at the DFT(B3LYP)/6-311+G(d,p) level to establish the stability order of all possible 18 tautomers-rotamers in the gas phase and to characterize their internal effects, polarity and aromatic character. Although the di-NH form strongly predominates (100%) in the mixture, the NH–OH, di-OH and CH–NH forms can be also considered. The favored tautomer is moderately stabilized by intramolecular interactions (attractions of the NH and C=O groups); it is also moderately polar and moderately delocalized. Stability of the functional groups (both amide functions) seems to be more important than intramolecular interactions, polarity and aromaticity. This is probably the main factor that dictates the tautomeric preferences in the uracil molecule.

Non-resonance-assisted hydrogen bonding in hydroxymethylene and aminomethylene cyclobutanones and cyclobutenones and their nitrogen counterparts

The characteristics of intramolecular hydrogen bonds (IMHB) have been systematically analyzed for a series of 32 different enols of derivatives of cyclobutane, cyclobutene, and cyclobutadiene bearing oxygen and nitrogen functionalities, at the B3LYP/6-311+G(3df,2p)//B3LYP/6-311+G(d,p) level of theory. In those cases where two tautomers (interconnected by a hydrogen shift through the IMHB) exist, tautomer a, in which the HB-donor group (YH) is attached to the four-membered ring, is less stable than tautomer b, in which is the HB-acceptor (X) is the one attached to the four-membered ring. As expected the OH group behaves as a better HB-donor than the NH(2) group and the C==NH group as a better HB-acceptor than the C==O group, although the first effect clearly dominates. Accordingly, the expected IMHB strength follows the [donor, acceptor] trend: [OH,C==NH]>[OH,C==O]>[NH(2),C==NH]>[NH(2),C==O]. Quite unexpectedly, in all those compounds in which the functionality exhibiting the IMHB is unsaturated, the IMHB is weaker than in their saturated counterparts, verifying that the primary effect on the strength of the IMHB is the structure of the sigma-skeleton of the system. In the conjugated systems investigated here, the severe constraints imposed by the four-membered ring force the HB donor and acceptor to be far apart and the IMHB is rather weak. These geometrical constraints are less severe for the saturated analogues and the IMHB becomes stronger, confirming that the characteristics of the sigma-skeleton, and not the resonance-assisted hydrogen bond (RAHB) phenomenon, are the primary contributors to the strength of the IMHB in conjugated compounds.

Study on the kinetics of keto-enol tautomerism of p-hydroxyphenylpyruvic acid using capillary electrophoresis

The kinetics of keto-enol tautomerism of p-hydroxyphenylpyruvic acid (pHPP) as a model of alpha-carbonyl compounds in aqueous solution at room temperature (25 degrees C) was first investigated by capillary electrophoresis with UV detection at 200 nm. The two tautomers could be separated and detected within 3 min. Since the ketonization of enolic pHPP varied with the buffer composition and buffer pH, the kinetics of pHPP was studied under different conditions, and relevant distributing fractions of enolic pHPP, ketonization rate constants and half-life were determined. In addition, beta-CD played an important part in the separation of the two tautomers, thus, the interaction between pHPP and beta-CD was also investigated by electrochemical techniques.

Tautomeric equilibria and pi electron delocalization for some monohydroxyarenes--quantum chemical studies

Keto-enol tautomeric interconversions and variations of the pi-electron distribution were studied for 11 isolated monohydroxyarenes at the DFT(B3LYP)/6-311++G(2df,2p) level. For two monohydroxyarenes (phenol and 9-anthrol), the PCM model of solvation (water) was also applied to the DFT geometries. The geometry-based HOMA index was applied to estimate pi-electron delocalization in the keto and enol tautomeric forms. Thermodynamic parameters of tautomeric interconversions (DeltaET, DeltaGT, TDeltaST, pKT) were calculated to estimate relative stabilities of individual tautomers and their percentage contents in the tautomeric mixtures. In almost all cases, the aromatic enol forms are strongly favored. An exception is 9-anthrol, which prefers its keto form. The resonance stabilization of this form comes from the central ring. Generally, aromaticity is the main factor that influences tautomeric equilibria in monohydroxyarenes. Hydration effect is considerably smaller and it does not change the tautomeric preference.

Resonance-assisted hydrogen bonds: a critical examination. Structure and stability of the enols of beta-diketones and beta-enaminones

The characteristics of the intramolecular hydrogen bond (IMHB) for a series of 40 different enols of beta-diketones and their nitrogen counterparts have been systematically analyzed at the B3LYP/6-311+G(3df,2p)//B3LYP/6-311+G(d,p) level of theory. In some cases, two tautomers may exist which are interconnected by a hydrogen shift through the IMHB. In tautomer a the HB donor group (YH) is attached to the six-membered ring, while in tautomer b the HB acceptor (X) is the one that is attached to the six-membered ring. We found that changing an O to a N favors the a tautomer when the atom is endo and the contrary when it is exo, while the presence of a double bond favors the a tautomers. As expected, the OH group behaves as a better HB donor than the NH2 group and the C=NH group as a better HB acceptor than the C=O group, although the first effect clearly dominates. Accordingly, the expected IMHB strength follows the [donor, acceptor] trend: [OH, C=NH] > [OH, C=O] > [NH2, C=NH] > [NH2, C=O]. For all those compounds in which the functionality exhibiting the IMHB is unsaturated (I-type), the IMHB is much stronger than in their saturated counterparts (II-type). However, when the systems of the II-type subset, which are saturated, are constrained to have the HB donor and the HB acceptor lying in the same plane and at the same distance as in the corresponding unsaturated analogue, the IMHB is of similar or even larger strength. Hence, we conclude that, at least for this series of unsaturated compounds, the resonance-assisted hydrogen bond effect is not the primary reason behind the strength of their IMHBs, which is simply a consequence of the structure of the sigma-skeleton of the system that keeps the HB donor and the HB acceptor coplanar and closer to each other.