ABSTRACT

A theoretical study of the reactivity of quinoline-4-one derivatives is undertaken in order to understand the involved mechanisms. The calculations were carried out in gas phase and in N, N-Dimethylformamide (DMF) solution. The Density Functional Theory (DFT) with B3LYP functional associated to 6-311G (d) and 6-311+G (d) bases is used to perform these calculations. The results of the thermodynamic parameters showed that there is an equilibrium relation between the different tautomers. This equilibrium can be used to explain the failure to obtain tetrahydroquinoline from 5,8-dimethoxy-quinolin-4-one. Reactivity analysis from Frontier Molecular Orbitals theory and Fukui function calculations revealed that ketone forms are less reactive than enol ones. The methoxyl substituent decreases the acidity of the nitrogen and oxygen atoms of quinolin-4-one while the bromine increases the acidity of the same sites. These results foresee that nitrogen deprotonation in the case of the brominated compound is easier than in the case of methoxylated ones.

Keywords:

Quinolein-4-One, Tautomers, Equilibrium Constants, Global Descriptors, Dual Descriptors

1. Introduction

Life expectancy, as well as pollution of various kinds, associated with new technologies and different risk factors such as alcohol or tobacco, increase in the number of pathologies. Cancer is found among these pathologies [1] . Thus, many researchers focused themselves on new treatments finding. These researches help them to discover, the very interesting role of quinone compounds. They are often found in biological systems such as photosynthetic or enzymatic systems, and in chemistry (oxidation, conductive polymers, etc.) [2] [3] [4] . These molecules are present in many structures with high antitumoral potential [5] [6] [7] [8] [9] . According to the usefulness of the biological interest of this nucleus, several ways of synthesis have been invented in order to obtain it [6] [7] [8] [10] [11] [12] [13] . Quinones have the same basic motif which is cyclohexadiene type. They possess units having exocyclic double bonds, giving them a large number of physicochemical properties. The para-dimethoxyl tetrahydroquinoline ring is an essential element for the synthesis of these quinones. Thus, 5,8-dimethoxy-4,4-dimethyltetrahydroquinoline is presented as a key intermediate which is necessary for the synthesis of these quinones. Although that of tetrahydroquinoline is described, many synthetic ways have been envisaged with the purpose of its elaboration [11] [14] [15] [16] . Moreover, it has been shown that obtaining this key nucleus from 5,8-dimethoxy-quinolin-4-one by using the conditions of Reetz does not give the desired product [17] . These experimental studies have also shown that the reduction of the double bonds’ number of quinolin-4-one is made in order to obtain dihydroquinolin-4-one but it has not been observed [18] . A theoretical study of quinolein-4-one derivatives is therefore undertaken in order to try to understand and predict the stability and reactivity of these compounds. The calculations are carried out in the gaseous phase and in N, N-dimethylformamide (DMF) solution with the functional B3LYP of the DFT method associated to the following bases 6-311G (d) and 6-311+G (d). The experimental difficulties evoked in obtaining 5,8-dimethoxy-4,4-dimethyltetrahydroquinoline from quinoline lead us to envisage the theoretical calculation of enthalpies of formation, enthalpies of reaction and free reaction enthalpies. The equilibrium constants are determined. Currently, the method of functional density theory (DFT) has been accepted as a popular approach for the calculation of the structural characteristics and energies of molecules by the community [19] and for the efficiency and accuracy of the evaluation of molecular properties [20] . This work aims to determine theoretically the preferred sites of electrophilic and nucleophilic attack on the aromatic carbon atoms and heteroatoms contained in the benzimidazolyl nucleus by different quantum chemistry methods. The molecular structures of the tautomers of quinoline derivatives studied are shown in Figure 1 and Figure 2. Figure 1 shows the two studied tautomeric forms; that is, ketone and enol forms. As for Figure 2, it shows the ketone tautomers (Q-4-one, DMQ and DBrQ) and enol (Q-4-hydroxy, DMHQ and DBrHQ).

2. Material and Methods

2.1. Level of Calculation Theory

Theoretical calculations are made by using the software called Gaussian 03 [21] .

The different molecular structures have been optimized using the Functional Density Theory (DFT) [22] [23] [24] . The level of calculation used is B3LYP associated with bases 6-311G (d) and 6-311+G (d). The calculations are carried out gaseous phase and solution. N, N-Dimethylformamide (DMF) being one of the solvents used for the synthesis of these quinolines [18] , the calculations in solution are thus carried out in DMF. The Polarized Continuum Model (PCM) [25] [26] is taken as a solvation model. The optimization of the molecular geometries is followed by calculations of vibration frequencies at the same level of the theory in order to verify that the structures obtained are minimal.

2.2. Thermodynamic Parameters of Reaction

Knowledge of variations in energy contributions to internal energy which occurs between the products and the reagents helps to characterize energetically a chemical reaction. For a given energy parameter X, its variation is determined according to the following relation (1):

${\Delta }_{r}X={{\displaystyle \sum }}^{\text{}}X\left(\text{products}\right)-{{\displaystyle \sum }}^{\text{}}X\left(\text{reactants}\right)$ (1)

The thermodynamic determined quantities in this study are the reaction enthalpies ( ${\Delta }_{r}H$ ), the reaction free enthalpies ( ${\Delta }_{r}G$ ) and the equilibrium constants (K_eq). These quantities are calculated by using relations (2) and (3) below:

${\Delta }_{r}G={{\displaystyle \sum }}^{\text{}}{\Delta }_{f}G\left(\text{Products}\right)-{{\displaystyle \sum }}^{\text{}}{\Delta }_{f}G\left(\text{Reactants}\right)$ (2)

${\Delta }_{r}H={{\displaystyle \sum }}^{\text{}}{\Delta }_{f}H\left(\text{Products}\right)-{{\displaystyle \sum }}^{\text{}}{\Delta }_{f}H\left(\text{Reactants}\right)$ (3)

The equilibrium constant denoted K_eq characterizes the equilibrium state of a chemical reaction. The value of the equilibrium constant depends solely on the considered chemical reaction and the temperature. In this work, the equilibrium constants are determined at 298 K from the following relation:

$K_{eq}={\exp }^{\left(-{\Delta }_{r}G\left(T\right)/RT\right)}$

2.3. Reactivity Parameters Derived from Conceptual DFT

The Frontier Molecular Orbital (FMO) theory allows to understand the chemical reactivity. In particular, the prediction of the acidity of the different tautomers is made from the calculation of the deprotonation energies. The study of the reactivity of these molecules is carried out by taking into account the FMO [27] . This theory indicates the ease with which a molecule can be transformed or can react with other molecules. It gives a prediction of the excitation properties of a molecule. As a result, it constitutes quantum parameters for the determination of molecular reactivity [28] [29] . The smaller the energy gap between the HOMO and LUMO orbitals, the easier the molecule can be excited. Similarly a large gap HOMO-LUMO indicates a high stability of the studied molecule which implies its low reactivity during chemical reactions. The condensed Fukui functions are obtained using the Yang and Mortier procedure [30] [31] [32] based on the finite difference method. The different values of local descriptors are calculated from equations (4) and (5) below:

$f_k^{+}=q_k\left(N+1\right)-q_k\left(N\right)$ (4)

$f_k^{-}=q_k\left(N\right)-q_k\left(N-1\right)$ (5)

where: $q_k\left(N\right)$ : Electronic population of the atom k in the neutral form; $q_k\left(N+1\right)$ : Electronic population of the atom k in the anionic form; $q_k\left(N-1\right)$ : Electronic population of the atom k in the cationic form.

We can gather all the information contained in the two Fukui functions in a single function called dual descriptor. This function is represented by expression (6) which is below.

$\Delta f=f_k^{+}-f_k^{-}$ (6)

where Δf is a function which represent the reactivity of the molecule, it is positive in the electrophilic areas and negative in the nucleophilic ones [33] .

The chemical potential μ measures the tendency of the electron cloud to escape from the molecule. It is a global property of the molecular system which is equal to the slope of the total electronic energy E. This energy is a function of the number of electrons N with an external potential V (r) constant. The studied parameter corresponds to the opposite of the electronegativity χ as defined by Pauling and Mulliken [34] [35] . This parameter is determined by the relation below (7)

$\mu =\left(\frac{\partial E}{\partial N}\right)=-\chi$ (7)

μ can be expressed as a function of the ionization potential IP and the electronic affinity EA and it is calculated by using the following relation (8).

$\mu =-\frac{\text{IP}+\text{EA}}{2}=-\chi$ (8)

The hardness ɳ $\eta$ is determined from the first derivative of the chemical potential. And it must be noted that the softness S is the inverse of the hardness [36] [37] . They are given by the following relation (9)

$\eta =\left(\frac{\partial \mu }{\partial N}\right)=\left(\frac{{\partial }^{2}E}{\partial E^2}\right)=\frac{1}{S}$ (9)

According to the theory of acids and bases, developed by Pearson [38] , these quantities can be expressed as a function of both the Ionization Potentials (IP) and Electronic Affinity (EA). It permits to determine the relation (10)

$\eta =\frac{1}{S}=\left(\text{IP}-\text{EA}\right)/2$ (10)

In order to characterize other electronic behaviors of the systems, we have interested ourselves to other descriptors such as the electrophile index ω which evaluate the ability of a molecule to generate an electron transfer. It is calculated according to the following formula [39] .

$\omega =\frac{{\mu }^2}{2\eta }=\frac{{\chi }^2}{2\eta }$ (11)

These descriptors can be calculated in a very simple way in the framework of Koopmans’ approximation [40] according to which:

$\text{IP}=-E_{\text{HOMO}}\text{andEA}=-E_{\text{LUMO}}$

$E_{\text{HOMO}}$ and $E_{\text{LUMO}}$ are the energies of the HOMO and LUMO orbitals.

2.4. Parameters of Tautomers Acidity Prediction

De protonation energies are widely used for the prediction of Brönsted acidity in gaseous and aqueous phases [41] [42] [43] [44] . The deprotonation reaction is defined by the following relationship (12).

$AH\to A^{-}+H^{+}$ (12)

The Gibbs energy $\Delta H_0^{298}$ of the protonic abstraction reaction is as follows (13):

$\Delta G_0^{298}=\Delta H_0^{298}-T\Delta S_0^{298}$ (13)

The contribution of entropy is given by the expression (14):

$-T\Delta S_0^{298}=-T\left[S\left(A^{-}\right)+S\left(H^{+}\right)-S\left(AH\right)\right]$ (14)

For a temperature of T = 298 K at standard pressure, the second term $TS\left(H^{+}\right)=32.5\text{Kj}{\text{.mol}}^{\text{-1}}$ gives the expression (15) [42] [43] .

$\Delta G_0^{298}=\Delta H_0^{298}-T\left[S\left(A^{-}\right)-S\left(AH\right)\right]-32.5$ (15)

which leads to the following relation (16):

$\Delta G_0^{298}=G\left(A^{-}\right)-G\left(AH\right)-32.5$ (16)

Previous theoretical work on the prediction of acidity has shown that gradient-corrected functional and hybrid functionals such as B3LYP associated with bases 6-311+G(d) and 6-311+G(3df, 3pd) lead to results in good agreement with the experimental data [45] [46] [47] [48] . This motivates our choice of the level of calculation implemented in this work.

3. Results and Discussion

The ketone-enol tautomeric equilibrium in the case of quinolin-4-one and the tautomeric structures which were the subject of this study are shown respectively in Figure 1 and Figure 2.

The performed calculations highlight the following parameters:

1) Thermodynamic parameters (Table 1);

2) The energies of boundaries molecular orbital and energetic gaps (Table 2);

3) Reactivity indices (Table 3);

4) The functions of Fukui (Table 4);

5) The deprotonation energies (Table 5).

3.1. Thermodynamic Parameters and Relative Stability

The calculated thermodynamic parameters are free enthalpies and enthalpies of formation. From these thermodynamic quantities, the reaction enthalpies and the equilibrium constants are determined. The order of relative stability of the tautomers is established from the enthalpy energies of formation of the different compounds. The values of the thermodynamic quantities calculated in gaseous phase and in DMF solution are gathered in Table 1.

Keq1 and Keq2 are respectively equilibrium constants in directions 1 and 2

(Figure 1). The values of the equilibrium constants mentioned in Table 1 are approximately equal to 1 which indicates that the reaction is in its equilibrium domain. Experimental data have shown that the probable existence of a tautomer relationship could explain the failure of tetrahydroquinoline to be obtained from 5,8-dimethoxy-quinolin-4-one by Reetz conditions [17] . Equilibrium constancy calculations confirm this hypothesis. They highlight a tautomeric equilibrium relationship between different tautomers. For two given tautomers, the results show that the ketone forms (Q-4-one, DMQ, DBrQ) have the lowest formation enthalpy values (Table 1). These are therefore more stable than the enol forms. In the case of 5,8-dimethoxy-quinolin-4-one, the enthalpies of formation of tautomers have very similar values. Both ketone and enol forms have a very similar stability. The reaction enthalpies are all positive which means that the equilibrium reaction is an endothermic one.

3.2. Chemical Reactivity

3.2.1. Theory of Frontier Molecular Orbitals

Table 2 contains the energies of the HOMO LUMO boundary orbitals and the HOMO-LUMO energy gaps.

The energy gap values decrease in the DMF solvent (Table 2) compared with those obtained in the gas phase. The molecules are therefore more reactive in solution. The obtained values of the ketone tautomers are greater than those obtained with enols in the two areas. The ketone forms are less polarizable than tautomers ones. And they have low chemical reactivity with high kinetic stability compared to enol tautomers. The substituents on quinoline cause a decrease in the energy gap. The smallest energy gap is obtained with bromine. This substituent therefore increases the reactivity of quinolin-4-one. At the enol form level, the quinoline-4-one methoxylated has the lowest value of ΔE which means that the methoxyl group leads to the most reactive enol form. We note however that among two tautomers, the different energy gaps are close to each other. This denotes that the tautomers reactivity is very closed to each other. This can explain the equivalence of N- and O-benzylation yielding rates in the benzylation reaction of quinolin-4-one [17] .

3.2.2. Global Indices of Reactivity

The study of the global reactivity of tautomers is based on the computation of global indices deduced from electronic properties. Global indices of the reactivity of studied tautomers are shown in Table 3.

According to these parameters, the chemical reactivity varies with the structure of the tautomers. Global softness (S) values of ketonic tautomers are lower than enol ones, indicating that ketone tautomers are soft. Also, we note that ketones have their electronegativity values (χ) lower than those of enols; except (Q-4-one) in DMF (3.823 eV). They are therefore good donors of electrons. In terms of ketone forms, the brominated derivative has the greatest value of electronegativity. It is therefore the best electron acceptor compared to the methoxylated derivatives. These data show that the bromine atom contributes to increase this index of reactivity. Concerning the methoxyl substituent it causes a decrease in the value of the same parameter. In addition, electrophile index

values of ketonic tautomers are lower than those of enol indicating that ketonic tautomers are less electrophilic.

3.2.3. Fukui Functions

The calculation of Fukui’s functions was performed using Hirschfeld populations [49] . The results are summarized in Table 4. This calculation relates to all the atoms of the quinolin-4-one ring (Figure 1).

The calculated Fukui functions indicate positive values of Δf at the C2, C3, C4, C5, C7 and C8 carbon atoms in all molecules (Table 4). These atoms are electrophilic sites. A nucleophilic attack in these molecules will happen on C4 and C5 carbons which have the highest values of Δf. Failure of C4 carbon dimethylation and reduction of the C2 = C3 double bond of the quinoline ring may be due to the fact that these atoms are nucleophilic sites. They are therefore favorable to an electrophilic attack. The N1, C8, C9 and O11 atoms have a negative value of Δf. These atoms are nucleophilic sites. They are susceptible to electrophilic attack. The nitrogen atom N1 has the highest negative value of Δf which means that an electrophilic attack will happen preferentially on this atom. The probable sites of nucleophilic and electrophilic attack are the same in all the molecules. The substitution does not modify the reactivity sites of these molecules.

3.3. Prediction of Acidity of Tautomers of Quinoline-4-one Derivatives

Gibbs free energies of deprotonation ( $\Delta G_0^{298}$ ) are calculated in the gas phase at B3LYP/6-311+G (d). The values of this thermodynamic quantity are grouped together in Table 5.

The grouped data in Table 5 show an increase in the free enthalpy of the deportation reaction with the methoxyl group. This energy decreases with bromine.

According to these data, we can establish the increasing order of the studied molecules by taking into account the free energy of deprotonation ( $\Delta G_0^{298}$ ) at 298 K: This order is: ΔG[DBrQ (O¹¹-H)] < ΔG[DBrQ (N¹-H)] < ΔG[Q-4-one (O¹¹-H)] < ΔG[Q-4-one (N¹-H)] < ΔG[DMQ (O¹¹-H)] < ΔG[DMQ (N¹-H)]. Bromine contributes to reducing the free energies of deprotonation of the nitrogen and oxygen atoms while the methoxyl group increases the same thermodynamic parameter on the same atoms. These results indicate that a deprotonation of the N¹-H nitrogen in the case of brominated quinolin-4-one will be easier than that of the nitrogen (N¹-H) of the methoxylated derivative. Work has shown that there is an inverse relationship between the amplitude of ΔG and the strength of the acid. The lower the value of ΔG, the more acidic the site is [41] [42] [43] [44] . In this way, the nitrogen atoms N¹-H and oxygen O¹¹-H in the brominated compounds are more acidic. These sites are less acidic in the methoxylated derivatives. Through the analysis of these free energies of deprotonation, we can establish the following decreasing order of studied molecules according to acidity: DBrQ (O¹¹-H) > DBrQ (N¹-H) > Q-4-one (O¹¹-H)> Q-4-one (N¹-H)] > DMQ (O¹¹-H)> DMQ (N¹-H). For a given quinoline, the deprotonation energy of the hydroxyl (O¹¹-H) of the enol form is lower than that of the nitrogen (N¹-H) of the ketone form.

4. Conclusion

The results of this theoretical study on the quinoline-4-one derivatives have shown that the ketone forms are therefore more stable than the enol ones. Experimental data have shown that the probable existence of a tautomerism can explain the failure to obtain tetrahydroquinoline from 5,8-dimethoxy-quinolin-4-one. Equilibrium constancy calculations confirmed this hypothesis. They highlight equilibrium relationship between different tautomers which is an endothermic reaction. Reactivity analysis from boundary molecular orbital theory has also shown that bromine (Br) increases the reactivity of quinoline. On the other hand, brominated quinolin-4-one is more reactive than its methoxylated analogues. Calculation of Fukui functions showed that C2 and C3 carbon atoms are electrophilic sites. This fact is a possible explanation for the failure of the reduction of the double bond in quinolin-4-one. These results indicate that deprotonation of N¹-H nitrogen in the case of brominated quinoline would be easier than that of nitrogen (N¹-H) of the methoxylated derivatives. The results predict a higher acidity of the N¹-H and O¹¹-H sites in the case of the brominated compounds. A theoretical study of the reaction mechanism of the N¹-H nitrogen deprotonation reaction of these molecules is envisaged in order to bring more precision to this first conclusion.

Acknowledgements

The authors would like to thank Lamoussa Ouattara for his contribution to the realization of this work.

Cite this paper

Bédé, A.L., Assoma, A.B., Yapo, K.D., Koné, M.G.-R., Koné, S., Koné, M., N’Guessan, B.R. and Bamba, E.-H.S. (2018) Theoretical Study by Density Functional Theory Method (DFT) of Stability, Tautomerism, Reactivity and Prediction of Acidity of Quinolein-4-One Derivatives. Computational Chemistry, 6, 57-70. https://doi.org/10.4236/cc.2018.63005

References