First

Scientific Reports volume 13, Article number: 8909 (2023) Cite this article

394 Accesses

2 Altmetric

Metrics details

The tautomerizations mechanism of 4-(methylsulfanyl)-3[(1Z)-1-(2-phenylhydrazinylidene) ethyl] quinoline-2(1H)-one were inspected in the gas phase and ethanol using density function theory (DFT) M06-2X and B3LYP methods. Thermo-kinetic features of different conversion processes were estimated in temperature range 273–333 K using the Transition state theory (TST) accompanied with one dimensional Eckert tunneling correction (1D-Eck). Acidity and basicity were computed as well, and the computational results were compared against the experimental ones. Additionally, NMR, global descriptors, Fukui functions, NBO charges, and electrostatic potential (ESP) were discussed. From thermodynamics analysis, the keto form of 4-(methylsulfanyl)-3-[(1Z)-1-(2 phenylhydrazinylidene) quinoline-2(1H)-one is the most stable form in the gas phase and ethanol and the barrier heights required for tautomerization process were found to be high in the gas phase and ethanol ~ 38.80 and 37.35 kcal/mol, respectively. DFT methods were used for UV–Vis electronic spectra simulation and the time-dependent density functional theory solvation model (TDDFT-SMD) in acetonitrile compounds.

Quinolinones, a quinoline heterocyclic analogue, have gotten a lot of attention for their physical, chemical, and biological activity as a treatment for sexually transmitted diseases, genitourinary infections, respiratory diseases, skin and soft tissue infections, prostate, and gastroenteritis1. Quinolinones have molluscicidal, fungicidal, and bactericidal activity; anti-HSV; anti-convulsion; anti-tumor; anti-oxidation; and anti-inflammatory activity2,3,4,5,6,7.

The position of the heterocyclic category depicted our consideration of the treatment of 3-acetyl-4-(methylsulfanyl)-quinolin-2(1H)-one with phenylhydrazine in boiling ethanol, which afforded a pale brown product that was identified as phenylhydrazone with a Z/E isomer ratio of 65:35. According to the literature8 elemental analysis revealed that compound 4-(methylsulfanyl)-3[(1Z)-1-(2-phenylhydrazinylidene) ethyl] quinoline-2(1H)-one lacked sulphur all at once.

Tautomeric phenomena occur when a molecule has multiple isomers. In organic chemistry, molecular biology, medicinal chemistry, and pharmacology, intramolecular H-atom transfer between two mixed equilibrium structures is extremely important. Tautomerism occurs in this work through the migration of an H-atom from one location to another. H-atom transfer produces keto-enol tautomers in the system under discussion. Tautomers are created when hydrogen atoms are exchanged between the N and O atoms of the heterocyclic ring9,10. Proton transport and H-bonding are important properties of the H-atom in chemistry. Theoretical research using the DFT approach, computational studies, and electrical, optical, and photoelectrical characteristics is very important to finding new drug candidates and understanding the electrical properties of different molecular structures11,12,13,14.

Due to lack of data of 4-(methylsulfanyl)-3[(1Z)-1-(2-phenylhydrazinylidene) ethyl] quinoline-2(1H)-one and the compound is partially soluble in water and completely soluble in ethanol8. We present here a computational study on the structures, relative stabilities, and electronic absorption spectra of tautomeric forms of 4-(methylsulfanyl)-3[(1Z)-1-(2-phenylhydrazinylidene) ethyl] quinoline-2(1H)-one (Scheme 1) in the gas phase and ethanol at B3LYP with 6-31G(d,p), 6-311++G(2d,2p) basis sets, and M06-2X/6-311++G(2d,2p) level. Proton transport and hydrogen bonding are two crucial factors that will be discussed.

The tautomeric structure keto and enol-form of 4-(methylsulfanyl)-3[(1Z)-1-(2 phenylhydrazinylidene) ethyl] quinoline-2(1H)-one.

There is a strong link between structure and stability. Understanding the chemical and physical features of these tautomers could aid future experimental studies on their expected applications, especially their ability to form metal complexes for analytical and biological purposes15,16. By comparing the resulting spectra against the experimental data, it is possible to understand the nature of the observed spectra and all its features17.

Becke's three-parameter and Lee–Yang–Parr hybrid functional (B3LYP) density functional theory (DFT) was used in conjunction with the 6-31G(d,p) basis set18,19,20 to fully optimize tautomer's and rotamers for their interconversions (detailed structures are given in Supplementary Table S1 in the Supporting information (SI) file). To characterize the nature of each stationary point on the potential energy surface, vibrational mode calculations were performed at the same level of optimization (B3LYP/6-31G(d,p) level). For more accurate energies, the B3LYP/6-31G (d, p) optimized structures were refined at the B3LYP/6-311++G (d, p) and meta-hybrid generalized gradient approximation M06-2X/6-311++G (d, p) levels of theories21,22.

For the transformation of keto—enol and enol—rotomar processes, rate coefficients were estimated in temperature range 273–333 K and calculated for the unimolecular reactions (kuni, in s−1) using the transition state theory (TST)23,24,25 at the high pressure (HP) limits. TST was recently used in many advanced studies26,27,28,29,30,31,32 and can be obtained from Eq. (1).

where σ, χ(T), kB, T, h, R, P°, and Go(T) refer to the reaction path degeneracy, tunneling correction, Boltzmann constant, temperature in Kelvin, Planck constant, ideal gas constant, the standard pressure and standard Gibbs free energy of activation for reaction, respectively. For unimolecular reactions, Δn = 0.

Tunneling coefficient χ(T) was estimated using the 1D-Eckart tunneling correction (Eck)33 that previously mentioned in many studies29,30,31,32. 1D-Eckart tunneling correction can be obtained numerically by integration of probability of transmission (p(E)) over a Boltzmann distribution of energy as given in Eq. (2):

where ΔHf ≠,0K and is the zero-point corrected activation enthalpy in the forward direction.

The orbital interactions, atomic charges, and their effects on the structure and stability of the examined structures were calculated using the natural bond orbital (NBO) approach34 using NBO program version 3.135.

Optimization was also performed in ethanol at the B3LYP/6-31G (d, p) level using the solvation model based on density (SMD)36. The solvation effect in ethanol was computed using B3LYP/6-311++G (2d, 2p) and M06-2X/6-311++G (2d, 2p) levels. The Gaussian 09W program37 was used to do all calculations. The study of tautomeric equilibria can be supported by nuclear magnetic resonance (NMR) spectroscopy. The gauge independent atomic orbital (GIAO) method was used to obtain NMR shielding relative to the 13C and 1H isotropic chemical shielding of tetramethyl silane (TMS) at the B3LYP/6-31G(d,p) optimized gas phase geometry in chloroform38,39,40. To recognize the stability and reactivity of the investigated structures, the global chemical reactivity descriptors were established from the highest occupied molecular orbital (HOMO) and lowest unoccupied molecular orbital (LUMO) energies41,42,43,44 consequently, electron affinity (EA), ionization potential (IP), absolute hardness (η), softness (S), electronegativity (X), chemical potential (\(\mu \)), electrophilicity index (\(\omega )\) have been estimated in the gas phase at the B3LYP/6-31G (d, p) level of theory. The global chemical reactivity descriptors have been obtained from the following relations, Eqs. (3–9):

The Fukui function is one of the main density functional descriptors that are used to model chemical reactivity and intramolecular site selectivity45. Fukui functions condensed (CFF) f+(r), f−(r), f0(r) of investigated structures are computed using individual atomic charges for natural population analysis (NPA).

Fukui functions can be extracted from Eqs. (10–12):

where f+(r), f0(r), and f-(r) represent nucleophilic, radical, and electrophilic attacks. The dual descriptor ∆f(r) resembles the difference between nucleophilic (Eq. (10)) and electrophilic attacks (Eq. (12)). For different sites, if ∆f(r) > 0, this site tends to undergo a nucleophilic attack. Chemical reactivity towards negative and positive charges could also be expected through mapping electrostatic potential (ESP).

For the investigated structures, the electronic absorption spectra (EAS) have been inspected using time-dependent density functional theory (TD-DFT) and the Perdew, Burke, and Ernzerhof (PBE) method (abbreviated as TD-PBE) in acetonitrile via the SMD approach at the optimized gas phase geometry of B3LYP/6-31G(d,p), level46,47. The accurate PBE0 method can be used to estimate electron excitations for different dye compounds48,49. Though, for the inspected molecule compared to experimental results, it gives an underestimation of λmax by 59 nm, whereas TD-PBE functionally overestimates and yields an 11–13 nm difference depending on basis sets. In order to draw the ultraviolet–visible (UV–Vis) spectra, Gauss Sum program50 was used. For more accurate results, the natural transition orbitals (NTOs)51 was investigated for different electrons excitations instead of the canonical orbitals. The NTOs were sketched using Gaussview52.

Previous studies demonstrated the accuracy of M06-2X functional in predicting the stability of tautomers and conformers21,22,53,54. Hence, the structures will be considered at B3LYP/6-31G (d, p) and energies at B3LYP/6-31G (d, p), B3LYP/6-311++G (d, p), and M06-2X/6-311++G (d, p). Three isomers of 4-(methylsulfanyl)-3[(1Z)-1-(2 phenylhydrazinylidene) ethyl] quinoline-2(1H)-one are discussed. Figure 1 shows the optimized structures of different structures at B3LYP/6-31G (d, p) level. To estimate and understand the stability order of the investigated system, reliable structures are necessary. A comparison between theory and experiment must be made to explain the reliability of the obtained data. The B3LYP/6-31G (d, p) level nearly reproduces the structure as obtained from X-ray17 and this supports the reliability of this level for structure optimization.

The optimized structures of keto-enol tautomers, transition state (TS) and rotamer (Rot) of 4-(methylsulfanyl)-3[(1Z)-1-(2 phenylhydrazinylidene) ethyl] quinoline-2(1H)-one in gas phase and ethanol (in parentheses) at B3LYP/6-31G(d,p) level of theory.

Scheme 1 depicts two tautomeric forms of 4-(methylsulfanyl)-3[(1Z)-1-(2 phenylhydrazinylidene) ethyl] quinoline-2(1H)-one with an intramolecular hydrogen bond (HB). The enol (Quinolinol, E-form) form indicates a phenol-imine HB (O–H···N, 1.816 Å) between the H- atom of the phenolic group and the N atom of the quinoline moiety. It is significant that the short donor–acceptor atom distance in both enol and keto forms provides a suggestion about the presence of a low-barrier hydrogen bond (LBHB), which is between 2.6 and 2.8 Å and has a noted biological importance55,56,57,58,59,60,61. The strength of HB in the tautomeric structures can be estimated from the donor (O15)-acceptor (N17) distance, which is 1.816 Å in enol form and 2.004 Å in keto form. Accordingly, the enol form has the strongest H-bond. In contrast, rotating the OH hydrogen atom to be far from the acceptor atom (nitrogen atom) gives the corresponding rotamer with no hydrogen bond and increases the donor–acceptor distance to 2.99 Å. Therefore, the keto-enol tautomers will be more stable than the rotamer.

Table 1 collects the barrier height and reaction energy for Keto–Enol reaction in gas phase and ethanol (in parentheses) at the B3LYP/6-31G (d, p) and M06-2X/6-311++G (2d, 2p)//B3LYP/6-31G (d, p) levels, while Fig. 2 shows the potential energy diagram using M06-2X energies. From M06-2X energies, the keto form is the most stable structure followed by rotamer in gas phase and in ethanol. The barrier heights required for the transformation process to the enolare 38.80 and 37.35 kcal/mol, respectively in gas phase and in ethanol, while the transformation of enol to rotamer accompanied barrier heights 12.54 and 11.54 kcal/mol and reaction energies 2.10 and 2.13 kcal/mol in gas phase and ethanol, respectively relative to keto form. The intrinsic reaction coordinates (IRC) and the potential energy changes during keto-enol conversion were drawn in Supplementary Figs. S1 and S2 (SI), respectively. From Supplementary Fig. S1, O–H bond is formed gradually with the breaking of the N–H bond and the two curves cross each other at s = 0 amu1/2 bohr. The formed N–C bond and the broken C–O bond are gently formed during the progress of the reaction.

Potential energy diagram (∆E0K, ΔE0K†, in kcal/mol) for Keto → Enol (TS1) and Enol → Rotamer (TS2) conversions in gas phase and ethanol (in parentheses) using M06-2X energies.

The calculated rate coefficients for the selected transitions at TST and 1D-Eck tunneling are given in Table 2. The results show that, the rate of the transformation of keto to enol in ethanol is higher about 12–20 times than in gas phase and a high effect for tunneling correction during the applied range of temperature especially for Keto–enol reaction compared to the enol-rotomar conversion.

In the NMR spectrum, the development of a low-field proton signal (high chemical shifts) is a well-known effect of forming a hydrogen bond with a sign for an LBHB62,63,64. Table 3 contains the complete data on NMR in enol, keto, and rotamer at B3LYP/6-31G(d,p) in comparison to experimental NMR. Figure 3 shows the calculated 13C and 1H NMR chemical shifts for 4-(methylsulfanyl)-3[(1Z)-1-(2 phenylhydrazinylidene) ethyl] quinoline-2(1H)-one in CHCl3 agree well with the experimental results8. Table 3 shows that C18 in the keto form has the highest chemical shift (171.33 ppm) that can be attributed to its closeness to the oxygen atom of the carbonyl group. Obviously, H1 of the keto form has the highest chemical shift (19.40 ppm) that attends the formation of strong LBHB (N1–H1…O1=C2), followed by H1 of the enol form (14.76 ppm). On the contrary, H1 of the rotamer has the lowest chemical shift (4.51 ppm). Hibbert and Emsley calculated 1H NMR chemical shifts to help predict the presence or absence of HB and distinguish LBHB for a proton chemical shift of up to 20 ppm (relative to TMS)65.

Sketches the experimental chemical shifts verses the estimated 13C and 1H NMR chemical shifts at the B3LYP/6-31G(d,p) level 4-(methylsulfanyl)-3[(1Z)-1-(2 phenylhydrazinylidene) ethyl] quinoline-2(1H)-one using the GIAO method in CHCl3.

Table 2 shows the energies of LUMO, HOMO, the energy gap (Eg), vertical IP, EA, \(\eta \),\(S\), \(\chi \), \(\mu \) and \(\omega \). The IP represents the ability to donate an electron, while the EA illustrates the ability to accept an electron. According to Table 4, the keto form has the smallest HOMO–LUMO gap (3.00 eV), followed by rotamer (3.53 eV). The enol form has the highest gap (3.78 eV). The smaller the energy gaps, the higher the reactivity of the molecule66,67. So, the keto form is expected to have high chemical reactivity, low hardness, and high softness compared to enol and their rotamers. The keto form can act as an electron donor and acceptor by having the highest HOMO energy (EHOMO = − 5.36 eV, the lowest value of the IP) and the lowest LUMO energy (ELUMO = − 2.36 eV, the largest value of the affinity). However, based on the calculations, ω and \(\eta \) they are better suited to act as strong electrophiles. A good electrophile has a high chemical potential as well as a low hardness68. Because the keto form has the highest electronegativity (χ = 3.86 eV) and the most negative chemical potential, it is the best electron acceptor.

To determine the most favorable place to add or remove an electron from a molecule, one looks at the chemical reactivity, which is one of the most fundamental questions. Electron density distributions are basic to understanding electrophilic and nucleophilic attacks. Table 5 shows the calculated condensed Fukui function value for enol and keto forms, and values for rotamers are given in the supporting information.

From Table 5 values, the electrophilic attack order for enol in the gas phase is C6 > N17 > C4 > C5 > C9 > C14 > C13 > C12. The C6 and N17 atoms have a higher f– value, indicating possible electron acceptor sites. There is some evidence that reactive electrophilic sites are primarily found on the hydrazinylidene ring. Conversely, for the nucleophilic attack, the reactivity order is C6 > C4 > C14 > C33 > C34. According to the dual descriptor (∆f) > 0 value for N17 and O15, these sites are favored for nucleophilic attack. As well as all hydrogen atoms, H is highly nucleophilic and H7 is highly electrophilic.

The C6, C4, C14, C12, N29, C5, C13, and C9 atoms in keto form are more sensitive sites for accepting electrons, and the C6, C14, C9, C5, C12, C13, C20, and C4 atoms are the most favorable sites for electron donation. Hence, heterocyclic rings are the most reactive sites for electron donor–acceptor interactions. The C20 and N29 atoms have a very positive dual descriptor (∆f) value, indicating a proclivity to donate electrons. According to the highly negativity dual descriptor (∆f), the C20 and C4 atoms are the most favorable sites for accepting electrons. Because it has the highest values of f− and f+ in both enol and keto forms, the C6 atom is a suitable site for both electrophilic and nucleophilic attacks.

For rotamer form, it is observed that C9, C5, C12, C13, C14, and C6 atoms have a higher f– value, which shows the possible site for electrophilic attack, and C6, C12, C4, C13, and C14 have a higher f+ value which indicates the possible site for nucleophilic attack. The dual descriptor (∆f) values of C12, C13, and C6 are highly positive; they have the tendency to donate electrons. On the contrary, N17 and O15 are susceptible to acquiring electrons.

The electrostatic potential (ESP) and charge distribution surface are broadly used to determine the reactivity of a given molecule and its expected interaction with other systems. The NBO charges of the investigated structures have been calculated in Table 6 in the gas phase at the B3LYP/6-31G (d, p) level. In agreement with the Fukui function results, the NBO charges calculations for 4-(methylsulfanyl)-3[(1Z)-1-(2 phenylhydrazinylidene) ethyl] quinoline-2(1H)-one show that higher negative charges are located on the N and O atoms and a higher positive charge is obtained by H. Three carbon atoms (C12, C13, and C14) have a positive charge due to the presence of highly electronegative nitrogen and oxygen atoms. Table 6 showing the enol form has the highest negative charges on the O and N atoms with the highest electronegativity. Thus, the enol form has the highest potential to act as a bidentate ligand.

ESP maps are described by the charged regions in the molecule. The different colours represent different values of the electrostatic potential. The potential grows in the following order: red, orange, yellow, green, and blue. The red colour in ESP maps represents the most negative electrostatic potential, while the blue colour reflects the most positive electrostatic potential regions. The ESP surfaces of investigated structures obtained using B3LYP/6-31G(d,p) are shown in Fig. 4. In the enol and keto forms, the ESP show the localization of a significant negative charge on the O atom, while the blue colour appears around the H7 and H8 atoms of the phenylhydrazinylidene ring. Therefore, the O atom has the highest electron donation ability toward metal ions. A significant blue colour on the H attached to O atoms in rotamers has been observed, and the red colour that exists in the region between the N and O atoms is due to the decline of H-bonds.

Molecular ESP surfaces of the investigated structures at B3LYP/6-31G(d,p).

The natural population analysis (NPA) carried out on the electronic structures of the selected transition states to describe the distribution and arrangement of different electrons in various sub-shells atomic orbits. The accumulation of charges on the atom and the accumulation of electrons in the core, valence and rydberg sub-shells are summarized in Supplementary Figs. S3–S5 and Supplementary Tables S2–S3 in SI file for the keto–enol conversion in the gas phase and ethanol (in parentheses). According to Supplementary Tables S2 and S3, the most electronegative center charge are accumulated on O15, N17, N20 and N27 atoms − 0.72488 (− 0.72666), − 0.25119 (− 0.24709), − 0.62136 (− 0.61801), and − 0.42898 (− 0.42863), respectively that indicate that these atoms are highly electronegative and tend to donate electrons, while the most electropositive center charge are C13 and S16 atoms with 0.63879 (0.63425), and 0.23081 (0.23448) which means these atoms tend to accept electrons.

Investigation of our compound shows it has two protons that are attached to either a nitrogen or oxygen atom. The acidity and basicity of any molecule are required to explain its structure, reactivity, and different chemical properties. Furthermore, learning the acidity constants (pKa) is important for estimating the equilibrium constant (K) of different reactions, specifically those involving proton transfers. When determining the pKa experimentally is difficult, Scheme 2 shows how the thermodynamic free energies cycle can be used to approximate computational methods. The protonated (cation) and deprotonated (anion) structures of the 4-(methylsulfanyl)-3[(1Z)-1-(2 phenylhydrazinylidene) ethyl] quinoline-2(1H)-one are depicted in Fig. 5. AH, is usually neutral, denoted \({AH}_{2}^{+}\), typically has a net charge of + 1, for the protonated form, while the corresponding enol/rotamer or keto. The hydroxyl group of the phenylhydrazinylidene ring has been rotated to give the rotamer structure, according to the optimized structure of the deprotonated form. Therefore, the energy of the rotamer will be studied throughout the acidity constant calculation from the protonated form. The deprotonated form denoted \({A}^{-}\), typically has a net charge of − 1.

Thermodynamic cycle of the protonated structures AH2+ and AH+ at different phases (gas (g) and aqueous (s)) for pKa calculation.

The optimized structure of protonated (\({AH}_{2}^{+},\) cation) and deprotonated (A-, anion) structures of 4-(methylsulfanyl)-3[(1Z)-1-(2 phenylhydrazinylidene) ethyl] quinoline-2(1H)-one at B3LYP/6-31G(d,p) level.

The equations used for calculating pKa values are given below (Eqs. (13–16)):

where Gi(g), ΔGi(solv), and Gi(aq) are the standard free energies of the species "i" in the gas phase, the solvation free energy of "i", and the free energy change in aqueous phase, respectively. The \({G}_{{H}^{+}}\left(g\right)\) and \({\Delta G}_{{H}^{+}}\left(solv\right)\) terms are − 6.28 kcal/mol69,70 and − 265.90 kcal/mol71,72, respectively.

To study the strength of the OH and NH bonds, Tables 5 and 6 presented different charges on the N and O atoms for enol, keto, and enol-rotamer structures. Table 7 collects the calculated acidity constant (pKa) for the protonated and deprotonated structures in ethanol at B3LYP and M06-2X levels. From the obtained results, the values show that the correlation between the computational and experimental acidity constants indicates that the B3LYP/6-311++G(2d,2p) level gives the nearest pKa value to the experimental values.

The experimental pKa of 4-(methylsulfanyl)-3[(1Z)-1-(2 phenylhydrazinylidene) ethyl] quinoline-2(1H)-one was 12.717, and the protonation structure estimated pKa is 17.48, 12.52, and 9.32 at B3LYP/6-31G (d, p), B3LYP/6-311++G (2d, 2p), M06 By comparing the pKa obtained by deprotonation, the enol form is less acidic than its keto and rotamer structures. This can be credited to some parameters such as the presence or absence of hydrogen bonds, the strength of the OH and NH bonds, and the stability of the resulting conjugate base upon deprotonation.

In the gas phase, the intrinsic basicity can be calculated from the proton affinity (PA), which is the negative of the protonation reaction of AH.

The higher energy of the keto and rotamer forms indicates that the enol form of 4-(methylsulfanyl)-3[(1Z)-1-(2 phenylhydrazinylidene) ethyl] quinoline-2(1H)-one is stabilized and dominant in nature.

The photo physics and photochemistry of any molecule can be determined by its serious ability to act as a dye or sensor. Table S4 collects the first λmax of the 4-(methylsulfanyl)-3[(1Z)-1-(2 phenylhydrazinylidene) ethyl] quinoline-2(1H)-one using different solvation models (PCM, CPCM, and SMD) and different DFT methods (B3LYP, CAM-B3LYP, PBE, ωB97X-D). The obtained theoretical results are compared against the experimental ones. The TDDFT-PBE/SMD model shows a good agreement with the 1st and 2nd maximum excitation peaks (Eex is the difference between the computational and experimental 1st maximum and 2nd excitation energies are 0.06/0.07 and 0.09 eV, respectively)53.

Table 8 presents the values of Eex and f and the transition configurations of the intense peaks for enol, keto, and rotamer, while Figs. 6 and 7 show the Ultraviolet–Vissible (UV/Vis) absorption spectra of different tautomeric structures of 4 (methylsulfanyl)-3[(1Z)-1-(2 phenylhydrazinylidene) ethyl] quinoline-2(1H)-one and the HOMO and LUMO plots for keto, enol and rotamer at TD-PBE/6-311+G(d,p) level, respectively.

The UV/Vis absorption spectra for different structures keto, enol and rotamer of 4 (methylsulfanyl)-3[(1Z)-1-(2 phenylhydrazinylidene) ethyl] quinoline-2(1H)-one at TD-PBE/6-311+G(d,p) level.

HOMO, LUMO plot of keto, enol and rotamer at TD-PBE/6-311+G(d,p) level.

The presence of HBs in chemical compounds can contribute to their chemical stability73,74 and is considered a strong reason for the strength, shifting, and broadening of the absorption peaks. Figure 6 depicts the keto structure with the strongest HB and a large blue shift (until 280 nm), followed by the enol structure (until 300 nm), and finally the lower shift in the rotamer (until 500 nm). The blue shift and hyperchromic effect for keto and enol forms can be returns to the presence of HBs17.

As shown in Table 4, the lower Eg of the keto, enol, and rotamer have maximums in the electronic absorption spectra of keto relative to enol and rotamer, which are bathochromically shifted by 17 and 66 nm, respectively. In the UV–Vis spectra, the peaks of keto, enol, and rotamer are extended over 348–438 nm, 272–500 nm, and 280–479 nm, respectively. The strong electronic absorption of keto is assigned to the HOMO−3 to LUMO and HOMO−2 to LUMO+1 transitions. The maximum absorption peak for enol, which appears at a lower wavelength than keto, is attributed to HOMO−2 to LUMO+1 and HOMO−1 to LUMO+1 transitions. The broad electronic absorption of rotamer can be attributed to the transition of HOMO−1 to LUMO+2.

To analyze the nature of UV–Vis absorption, Fig. 8 shows the NTOs of keto, enol, rotamar and transition states for the high-intensity excited states at the PBE/6-311+G(d,p) level with solvent effects of acetonitrile through SMD. In Fig. 8, the occupied NTOs are referred to "hole", while the unoccupied NTOs are the "particles" transition orbitals. NTOs can provide a simple description of the excited state rather than the canonical orbitals. For the investigated structures, the dominant transitions are expected to be π − π* and n − π* excitations which makes the analysis of excitations are difficult. However, from Fig. 8, the holes of NTOs can reproduce the bands given in Fig. 6 and Table 8 and the holes are seem to delocalize on the whole molecular structure while the particles NTOs are mainly delocalized on benzene rings which enhances the π–π* excitation.

Natural transition orbitals (NTOs) of keto, enol, rotamar and transition states for the excitation with significant oscillation strengths at the PBE/6-311+G(d,p) level with solvent effects of acetonitrile through SMD. The occupied (holes) and unoccupied (electrons) NTO pairs are only that have contribution more than 50% to each excited state (λ is eigenvalues of the pairs).

The density functional methods (DFT) were used to investigate the tautomer and related rotamers of 4 (methylsulfanyl)-3 [(1Z)-1-(2 phenylhydrazinylidene) ethyl] quinoline-2(1H)-one. B3LYP and M06-2× connected with 6-31G (d, p) and 6-311++G (2d, 2p) basis sets have been used for analysis of different structural properties, stability, and aromaticity. The obtained thermo-kinetic results show a relative stability for keto form compared to other forms and high barriers required for enol formation in gas phase and ethanol under the applied temperature range. Using Eckert tunneling correction indicates a great contribution in keto-enol conversion compered to enol-rotomar conversion. Different sites for the nucleophilic and electrophilic reactions were allocated using Fukui functions. UV absorption spectra in ethanol and gas phases were investigated using the time-dependent density functional theory (TD-DFT) methods B3LYP, PBE, PBE0, CAM-B3LYP, M06-2X, b97X-D, and CIS. The TDDFT-PBE/SMD approach shows good harmony with the 1st and 2nd maximum excitation peaks, and the transitions are π–π* transitions. The remarkable chemical shift of a proton at 14.76 and 19.40 ppm in the nuclear magnetic resonance spectrum has been attributed to the existence of low-barrier hydrogen bonds (LBHBs) for the enol and keto forms, respectively.

All data generated through this study are included in this manuscript and the Supporting Information file.

Sarveswari, S. V., Vijayakumar, R., Siva, R. & Priya, R. Synthesis of 4-hydroxy-2(1H)-quinolone derived chalcones, pyrazolines and their antimicrobial, in silico antimalarial evaluations. Appl. Biochem. Biotechnol. 175, 43–64 (2015).

Article CAS PubMed Google Scholar

El-Shennawy, A. M., Mohamed, A. H. & Abass, M. Studies on parasitologic and haematologic activities of an enaminone derivative of 4-Hydroxyquinolin-2(1H)-one against murine schistosomiasis mansoni. Medscape Gen. Med. 9, 15–33 (2007).

Google Scholar

Cheng, P. et al. Synthesis and in vitro anti-hepatitis B virus activities of 4-aryl-6-chloro-quinolin-2-one and 5-aryl-7-chloro-1,4-benzodiazepine derivatives. Bioorg. Med. Chem. Lett. 18, 3787 (2008).

Article CAS PubMed Google Scholar

Kreimeyer, A., Laube, B., Sturgess, M., Goeldner, M. & Foucaud, B. Evaluation and biological properties of reactive ligands for the mapping of the glycine site on the N-methyl-d-aspartate (NMDA) receptor. J. Med. Chem. 42, 4394–4404 (1999).

Article CAS PubMed Google Scholar

Aubry, A. et al. Novel Gyrase mutations in quinolone-resistant and -hypersusceptible clinical isolates of mycobacterium tuberculosis: Functional analysis of mutant enzymes. Antimicrob. Agents Chemother. 50, 104–114 (2006).

Article CAS PubMed PubMed Central Google Scholar

Ukrainets, V. et al. Khim. Geterotsikl. Soedin. 8, 1113 (1996) (Chem. Abstr. 126, 59840x (1997)).

Mohamed, E. A., Ismail, M. M., Gabr, Y. & Abass, M. Synthesis of some multiazaheterocycles as substituents to quinolone moiety of specific biological activity. Chem. Papers 48, 285–292 (1994).

CAS Google Scholar

Abass, M., Hassanin, H. M., Allimony, H. A. & Hassan, H. Substituted quinolinones 27. * Regioselective synthesis of pyrazolo-, oxazolo-, and triazepinoquinoline derivatives. Chem. Hetero. Compd. 51(11/12), 1023–1029 (2015).

Article CAS Google Scholar

Vila, J. A., Arnautova, Y. A., Vorobjev, Y. & Scheraga, H. A. Assessing the fractions of tautomeric forms of the imidazole ring of histidine in proteins as a function of PH. Proc. Natl. Acad. Sci. 108(14), 5602–5607 (2011).

Article ADS CAS PubMed PubMed Central Google Scholar

Minkin, V. I., Olekhnovich, L. P. & Zhdanov, Y. A. Molecular design of tautomeric compounds. Acc. Chem. Res. 3(14), 210–217 (1981).

Article Google Scholar

Dege, N. et al. Quantum computational, spectroscopic investigations on N-(2-((2-chloro-4,5-dicyanophenyl) amino) ethyl)-4-methylbenzenesulfonamide by DFT/TD-DFT with different solvents, molecular docking and drug-likeness research. Colloids Surf. A 638, 128311 (2022).

Article CAS Google Scholar

Abdulridha, A. A. et al. Corrosion inhibition of carbon steel in 1 M H2SO4 using new Azo Schiff compound: Electrochemical, gravimetric, adsorption, surface and DFT studies. J. Mol. Liq. 315, 113690 (2020).

Article CAS Google Scholar

Balakit, A. A. et al. Synthesis, spectrophotometric and DFT studies of new Triazole schiff bases as selective naked-eye sensors for acetate anion. Supramol. Chem. 32, 519–526 (2020).

Article CAS Google Scholar

Soliman, H. N. & Yahia, I. S. Synthesis and technical analysis of 6-butyl-3-[(4-chlorophenyl) diazenyl]-4-hydroxy-2H-pyrano[3,2-c] quinoline-2,5(6H)-dione as a new organic semiconductor: Structural, optical, and electronic properties. Dyes Pigm. 176, 108199 (2020).

Article Google Scholar

Buhl, M. & Wiillen, C. V. Computational evidence for a new C84 isomer. Chem. Phys. Lett. 247, 63–68 (1995).

Article ADS Google Scholar

Schleyer, P. et al. Nucleus-independent chemical shifts: A simple and efficient aromaticity probe. J. Am. Chem. Soc. 118, 6317–6318 (1996).

Article CAS PubMed Google Scholar

Oda, M. et al. Synthesis, molecular structure, and properties of 2-(2-Hydroxyphenyl)-1-Azaazulene. Eur. J. Org. Chem. 6, 2231–2236 (2012).

Article Google Scholar

Becke, A. D. A new mixing of Hartree-Fock and local-density-functional theories. J. Chem. Phys. 98(2), 1372–1377 (1993).

Article ADS CAS Google Scholar

Lee, C., Yang, W. & Parr, R. G. Development of the colle-salvetti correlation-energy formula into a functional of the electron density. Phys. Rev. B 37(2), 785–789 (1988).

Article ADS CAS Google Scholar

Stephens, P. J., Devlin, F. J., Chabalowski, C. & Frisch, M. J. Ab initio calculation of vibrational absorption and circular dichroism spectra using density functional force fields. J. Phys. Chem. 98(45), 11623–11627 (1994).

Article CAS Google Scholar

Zhao, Y. & Truhlar, D. G. The M06 suite of density functionals for main group thermochemistry, thermochemical kinetics, noncovalent interactions, excited states, and transition elements: Two new functionals and systematic testing of four M06-class functionals and 12 other function. Theor. Chem. Acc. 120(1–3), 215–241 (2008).

Article CAS Google Scholar

Zhao, Y., Schultz, N. E. & Truhlar, D. G. Design of density functionals by combining the method of constraint satisfaction with parametrization for thermochemistry, thermochemical kinetics, and noncovalent interactions. J. Chem. Theory Comput. 2(2), 364–382 (2006).

Article PubMed Google Scholar

Steinfeld, J. I., Francisco, J. S. & Hase, W. L. Chemical Kinetics and Dynamics 2nd edn. (Prentice-Hall, 1999).

Google Scholar

Carstensen, H. H., Dean, A. M. & Deutschmann, O. Rate constants for the H abstraction from alkanes (R–H) by R′ O2 radicals: A systematic study on the impact of R and R′. Proc. Combust. Inst. 31, 149–157 (2007).

Article Google Scholar

Fernandez, R. A., Ellingson, B. A., Garrett, B. C. & Truhlar, D. G. Variational transition state theory with multidimensional tunneling. J. Rev. Comput. Chem. 23, 125–232 (2007).

Article Google Scholar

Abdel-Rahman, M. A. et al. Mechanistic insights of degradation of O-Anisidine carcinogenic pollutant initiated by OH radial attack: Theoretical investigations. New J. Chem. 45, 5907–5924 (2021).

Article CAS Google Scholar

Shojaei, S. H. R., Shiroudi, A. & Abdel-Rahman, M. A. Computational studies on thermo-kinetics aspects of pyrolysis of isopropyl acetate and its methyl, bromide and hydroxyl derivatives. Heliyon 8, e11274 (2022).

Article CAS PubMed PubMed Central Google Scholar

El-Gogary, T. M., Heikal, L. A., Abdel-Rahman, M. A. & El-Nahas, A. M. First-principle kinetic studies of unimolecular pyrolysis of isopropyl esters as biodiesel surrogates. Theoret. Chem. Acc. 140, 1–15 (2021).

Article Google Scholar

Abdel-Rahman, M. A., Shibl, M. F. & Mahmoud, M. A. M. Pyrolytic elimination of ethylene from ethoxyquinolines and ethoxyisoquinolines: A computational study. Sci. Rep. 13, 6248 (2023).

Article ADS CAS PubMed PubMed Central Google Scholar

Abdel-Rahman, M. A., Shiroudi, A., Kaya, S. & El-Nahas, A. M. Theoretical investigations on the unimolecular decomposition mechanisms of isopropyl acetate. J. Mol. Struc. 1262, 133006 (2022).

Article CAS Google Scholar

Abdel-Rahman, M. A. et al. A W1 computational study on the kinetics of initial pyrolysis of a biodiesel model: Methyl propanoate. New J. Chem. 45, 19531–19541 (2021).

Article CAS Google Scholar

Al-Otaibi, J. S. et al. Ab initio-based kinetics of hydrogen atom abstraction from methyl propionate by H and CH3 radicals: A biodiesel model. Struct. Chem. (2021).

Eckart, C. The penetration of a potential barrier by electrons. Phys. Rev. 35, 1303–1309 (1930).

Article ADS CAS MATH Google Scholar

Reed, A. E., Weinstock, R. B. & Weinhold, F. Natural population analysis. J. Chem. Phys. 831(10), 1736–154104 (1985).

Article ADS Google Scholar

Glendening, E. D., Reed, A. E., Carpenter, J. E. & Weinhold, F. NBO Version 3.1. (Gaussian Inc., 2003). https://nbo6.chem.wisc.edu/.

Marenich, A. V., Cramer, C. J. & Truhlar, D. G. Universal solvation model based on solute electron density and on a continuum model of the solvent defined by the bulk dielectric constant and atomic surface tensions. J. Phys. Chem. B 113(18), 6378–6396 (2009).

Article CAS PubMed Google Scholar

Frisch, M. J., et al. Gaussian 09 Revision A.1 (Gaussian Inc., 2009).

Wolinski, K., Hinton, J. F. & Pulay, P. Efficient implementation of the gauge-independent atomic orbital method for NMR chemical shift calculations. J. Am. Chem. Soc 112, 8251–8260 (1990).

Article CAS Google Scholar

Wolff, S. K. & Ziegler, T. Calculation of DFT-GIAO NMR shifts with the inclusion of spin-orbit coupling. J. Chem. Phys. 109(3), 895–905 (1998).

Article ADS CAS Google Scholar

Cheeseman, J. R. et al. A comparison of models for calculating nuclear magnetic resonance shielding tensors. J. Chem. Phys. 104(14), 5497–5509 (1996).

Article ADS CAS Google Scholar

Lucila, L., Martinez, L., Mitnik, G. & Borunda, E. O. Theoretical calculation of UV–Vis, ir spectra and reactivity properties of tamoxifen drug: A methodology comparison. MOJ Biorg. Org. Chem. 1(3), 87–95 (2017).

Google Scholar

Flores, M. C., Márquez, E. A. & Mora, J. R. Molecular modeling studies of bromopyrrole alkaloids as potential antimalarial compounds: A DFT approach. Med. Chem. Res. 27, 844–856 (2018).

Article CAS Google Scholar

Mishra, V. R. & Sekar, N. Photostability of coumarin laser dyes—A mechanistic study using global and local reactivity descriptors. Fluoresc 27, 1101–1108 (2017).

Article CAS Google Scholar

Das, S., Shedge, S. V. & Pal, S. Critical study of the charge transfer parameter for the calculation of interaction energy using the local hard-soft acid-base principle. J. Phys. Chem. A 117(42), 10933–10943 (2013).

Article CAS PubMed Google Scholar

Morell, C. & Toro-labbe, A. New dual descriptor for chemical reactivity. J. Phys. Chem. A 109, 205–212 (2005).

Article CAS PubMed Google Scholar

Perdew, J. P., Burke, K. & Ernzerhof, M. Generalized gradient approximation made simple. Phys. Rev. Lett. 77, 3865–3868 (1996).

Article ADS CAS PubMed Google Scholar

Ernzerhof, M. & Perdew, J. P. Generalized gradient approximation to the angle- and system-averaged exchange hole. J. Chem. Phys. 109(9), 3313–3320 (1998).

Article ADS CAS Google Scholar

Adamo, C. & Barone, V. Toward reliable density functional nethods without adjustable parameters: The PBE0 model. J. Chem. Phys. 110(13), 6158–6170 (1999).

Article ADS CAS Google Scholar

Jacquemin, D., Perpe, E. A., Scuseria, G. E., Ciofini, I. & Adamo, C. TD-DFT performance for the visible absorption spectra of organic dyes: Conventional versus long-range hybrids. J. Chem. Theory Comput. 4(1), 123–135 (2008).

Article CAS PubMed Google Scholar

O’Boyle, N. M., Tenderholt, A. L. & Langner, K. M. Cclib: A library for package-independent computational chemistry algorithms. J. Comput. Chem. 29, 839–845 (2008).

Article PubMed Google Scholar

Martin, R. L. Natural transition orbitals. J. Chem. Phys. 118, 4775–4777 (2003).

Article ADS CAS Google Scholar

Dennington, R., Keith, T. & Millam, J. GaussView, Version 5 (Semichem Inc., 2009). https://gaussian.com/gaussview6/.

Michelson, A. Z., Petronico, A. & Lee, J. K. 2-Pyridone and derivatives: Gas-phase acidity, proton affinity, tautomer preference, and leaving group ability. J. Org. Chem. 77, 1623–1631 (2012).

Article CAS PubMed Google Scholar

Walker, M., Harvey, A. J. A., Sen, A. & Dessent, C. E. H. Performance of M06, M06–2X, and M06-HF density functionals for conformationally flexible anionic clusters: M06 functionals perform better than B3LYP for a model system with dispersion and ionic hydrogen-bonding interactions. J. Phys. Chem. A 117(47), 12590–12600 (2013).

Article CAS PubMed Google Scholar

Mckenzie, R. H. et al. Effect of quantum nuclear motion on hydrogen bonding. J. Chem. Phys. 140, 174508 (2014).

Article ADS PubMed Google Scholar

Kong, X. et al. Proton probability distribution in the O ··· H ··· O low-barrier hydrogen bond: A combined solid-state NMR and quantum chemical computational study of Dibenzoylmethane and Curcumin. J. Phys. Chem. B 120, 11692–11704 (2016).

Article CAS PubMed Google Scholar

Perrin, C. L. & Nielson, J. B. "Strong" Hydrogen bonds in chemistry and biology. Annu. Rev. Phys. Chem. 48(1), 511–544 (1997).

Article ADS CAS PubMed Google Scholar

Gilli, P., Bertolasi, V., Pretto, L.& Gilli, G. The nature of solid-state N-H...O/O-H...N tautomeric competition in resonant systems. intramolecular proton transfer in low-barrier hydrogen bonds formed by the ...O=C-C=N-NH... ...HO-C=C-N=N... Ketohydrazone-Azoenol system. A variable-temperature X-Ray Cr. J. Am. Chem. Soc. 124, 13554–13567 (2002).

Cleland, W. W., Frey, P. A. & Gerlt, J. A. The low barrier hydrogen bond in enzymatic catalysis. J. Biol. Chem. 273(40), 25529–25533 (1998).

Article CAS PubMed Google Scholar

Smallwood, C. J. & Mcallister, M. A. Characterization of low-barrier hydrogen bonds. 7. Relationship between strength and geometry of short-strong hydrogen bonds. The Formic acid—Formate anion model system. An Ab initio and DFT investigation. J. Am. Chem. Soc. 119, 11277–11281 (1997).

Article CAS Google Scholar

Sobczyk, L., Chudoba, D., Tolstoy, P. M. & Filarowski, A. Some brief notes on theoretical and experimental investigations of intramolecular hydrogen bonding. Molecules 21(12), 1657–1667 (2016).

Article CAS PubMed PubMed Central Google Scholar

Tolstoy, P. M. et al. Geometries and tautomerism of OHN hydrogen bonds in aprotic solution probed by H/D isotope effects on 13C NMR chemical shifts. J. Phys. Chem. A 114, 10775–10782 (2010).

Article CAS PubMed Google Scholar

Kumar, G. A. & Mcallister, M. A. Theoretical investigation of the relationship between proton NMR chemical shift and hydrogen bond strength. J. Org. Chem. 63, 6968–6972 (1998).

Article CAS PubMed Google Scholar

Garcia-viloca, M., Gelabert, R., Gonzalez-Lafont, A., Moreno, M. & Lluch, J. M. Is an extremely low-field proton signal in the NMR spectrum conclusive evidence for a low-barrier hydrogen bond?. J. Phys. Chem. A 101(46), 8727–8733 (1997).

Article CAS Google Scholar

Hibbert, F. & Emsley, J. Hydrogen bonding and chemical reactivity. J. Adv. Phys. Org. Chem. 26, 255–379 (1990).

CAS Google Scholar

Cai, T. et al. Structure and enhanced reactivity Rates of the D5h Sc3N@C80 andLu3N@C80 metallofullerene isomers: The importance of the pyracylene motif. J. Am. Chem. Soc. 128, 8581–8589 (2006).

Article CAS PubMed Google Scholar

Chermette, H. Chemical reactivity indexes in density functional theory. J. Comput. Chem. 20, 129–154 (1999).

3.0.CO;2-A" data-track-action="article reference" href="https://doi.org/10.1002%2F%28SICI%291096-987X%2819990115%2920%3A1%3C129%3A%3AAID-JCC13%3E3.0.CO%3B2-A" aria-label="Article reference 67" data-doi="10.1002/(SICI)1096-987X(19990115)20:13.0.CO;2-A">Article CAS Google Scholar

Sánchez-márquez, J. et al. Introducing, "UCA-FUKUI" software: Reactivity index calculations. J. Mol. Model. 20, 1–13 (2014).

Article Google Scholar

Topol, I. A., Tawa, G. J., Burt, S. K. & Rashin, A. A. On the structure and thermodynamics of solvated monoatomic ions using a hybrid solvation model. J. Chem. Phys. 111(24), 10998–11014 (1999).

Article ADS CAS Google Scholar

Jang, Y. H., Sowers, L. C., Çaǧin, T. & Goddard, W. A. First principles calculation of PKa values for 5-substituted uracils. J. Phys. Chem. A 2001(105), 274–280 (2001).

Article Google Scholar

Tissandier, M. D. et al. The proton's absolute aqueous enthalpy and gibbs free energy of solvation from cluster-ion solvation data. J. Phys. Chem. A 102(40), 7787–7794 (1998).

Article CAS Google Scholar

Camaioni, D. M. & Schwerdtfeger, C. A. Comment on "Accurate experimental values for the free energies of hydration of H+, OH-, and H3O+". J. Phys. Chem. A 109, 10795–10797 (2005).

Article CAS PubMed Google Scholar

Abdel-Rahman, M. A., Al-Hashimi, N., Shibl, M. F., Yoshizawa, K. & El-Nahas, A. M. Thermochemistry and kinetics of the thermal degradation of 2-methoxyethanol as possible biofuel additives. Sci. Rep. 9, 4535 (2019).

Article ADS PubMed PubMed Central Google Scholar

Abdel-Rahman, M. A. et al. Computational studies on the thermodynamic and kinetic parameters of oxidation of 2-methoxyethanol biofuel via H-atom abstraction by methyl radical. Sci. Rep. 9, 15361 (2019).

Article ADS PubMed PubMed Central Google Scholar

Download references

Open access funding provided by The Science, Technology & Innovation Funding Authority (STDF) in cooperation with The Egyptian Knowledge Bank (EKB).

Chemistry Department, Faculty of Education, Ain Shams University, Roxy, Cairo, 11711, Egypt

Shimaa Abdel Halim

Chemistry Department, Faculty of Science, Suez University, Suez, 43518, Egypt

Mohamed A. Abdel-Rahman

You can also search for this author in PubMed Google Scholar

You can also search for this author in PubMed Google Scholar

The authors contributed equally to the manuscript.

Correspondence to Shimaa Abdel Halim or Mohamed A. Abdel-Rahman.

The authors declare no competing interests.

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/.

Reprints and Permissions

Halim, S.A., Abdel-Rahman, M.A. First-principles density functional theoretical study on the structures, reactivity and spectroscopic properties of (NH) and (OH) Tautomer's of 4-(methylsulfanyl)-3[(1Z)-1-(2-phenylhydrazinylidene) ethyl] quinoline-2(1H)-one. Sci Rep 13, 8909 (2023). https://doi.org/10.1038/s41598-023-35933-8

Download citation

Received: 26 November 2022

Accepted: 25 May 2023

Published: 01 June 2023

DOI: https://doi.org/10.1038/s41598-023-35933-8

Anyone you share the following link with will be able to read this content:

Sorry, a shareable link is not currently available for this article.

Provided by the Springer Nature SharedIt content-sharing initiative

By submitting a comment you agree to abide by our Terms and Community Guidelines. If you find something abusive or that does not comply with our terms or guidelines please flag it as inappropriate.