CC  Vol.9 No.1 , January 2021
DFT Study of Dimerization Sites in Imidazo[1,2-a]pyridinyl-chalcone Series
Abstract: Quantum chemistry methods were performed in order to characterize the chemical reactivity on series of imidazo[1,2-a]pyridinyl-chalcone (IPC). In particular, the B3LYP/6-311G(d) theory level has been used to determine parameters which characterize the global and local reactivity on five molecules of the series. These compounds differ from one to another with the aryl groups. There are: 1-(2-methylimidazo[1,2-a]pyridin-3-yl)-3-phenylprop-2-en-1-one, 3-(4-fluorophenyl)-1-(2-methylimidazo [1,2-a]pyridin-3-yl)prop-2-en-1-one, 3-[4-(dimethylamino)phenyl]-1-(2-methylimidazo [1,2-a]pyridin- 3-yl)prop-2-en-1-one, 3-(2,4-dichlorophenyl)-1-(2-methylimidazo [1,2-a]pyridin-3-yl)prop-2-en-1-one, 3-(2,4-dichlorophenyl)-1-(2-methylimidazo [1,2-a]pyridin-3-yl)prop-2-en-1-one. All results lead to finding out that local nucleophilicity and electrophilicity of compounds are not substituent-dependant contrarily to their global nucleophilicity which prove to be more sensitive to the electron-donating character of the substituents. 3-[4-(Dimethylamino) phenyl]-1-(2-methylimidazo[1,2-a]pyridin-3-yl)prop-2-en-1-one was identified as the unique nucleophile compound by global reactivity. Respectively, the carbon atoms C5 and C14 are the prediction sites of electrophilic and nucleophilic attacks in the molecular skeleton of both molecules. Identification of interactions centres on IPC series is of great importance for organic synthesis and medicinal chemistry where the molecular hybridization strategy is very often used to improve biological activities of interesting therapeutic systems.

1. Introduction

Chalcones bring together a plethora of natural and synthetic compounds, within which a keto-ethylenic system links two aromatic nuclei (aryl and/or heteroaryl rings) [1] from the carbon atom in β-position and that of the carbonyl function. This family of compounds arouses great interest in various fields of applications. Their true importance is based on the wide spectrum of biological activities [2] [3] [4] and on the many synthetic perspectives offered by their skeleton. The biological activities displayed by these compounds are mainly linked to the functional group enone [5] as well as to the electronic nature of the nuclei A and B due to the substituents [6] [7] [8] [9]. Imidazopyridinyl-chalcones (IPC) [10] are synthetic chalcones obtained with imidazo[1,2-a]pyridine as ring A and an aryl as ring B. Imidazopyridine nucleus is the building block in the design of several drugs [11]. In the majority of cases, chalcones have biological activities similar to those of imidazopyridine and derivatives. Undoubtedly, the chalcones carrying imidazopyridine nucleus prove to be potential sources of more biologically active molecules due to the possible conjugation effects of the two entities involved. Imidazopiridinyl-chalcones showed high nematicidal activities against drug-resistant strains of nematodes [10]. Their remarkable therapeutic properties against cancer [12] and certain bacteria [13] are also reported. These biological potentials coupled with the possible molecular diversity define a pharmacological dynamic profile which, increasingly, motivates studies of physicochemical characterizations [14] [15] and structure-activity relationships [16] in IPC series. Also, like most chalcones, the compounds IPC are susceptible to addition reactions [17]. Molecular dimerization based on the IPC entities is therefore a path that deserves serious attention for the development of novel therapeutic agents which could be probably stronger and less toxic. However, any reaction to be carried out requires precise knowledge of the chemoselectivity, the stereoselectivity and the regioselectivity. In this area, the tools of quantum chemistry and quantum theories of chemical reactivity are of great help in predicting and justifying the different selectivity. Among the theories developed and validated, nucleophilicity and electrophilicity concepts have proven to be solid in the study of molecular reactivity [18] [19] [20] [21] [22]. The different parameters defined by these concepts are global and local. In this work, we are interested in five molecules from the IPC series; four of which are substituted from the IPC nucleus by varying the aryl substituents. In Table 1, these compounds are summarized with their respective nematicidal potential pLC100.

The nematicidal potential pLC100 was evaluated following the expression:

p L C 100 = log 10 ( L C 100 M 10 3 ) (i)

where M (in g/mol) is the molecular molar mass and LC100 (in µg/mol), the lowest concentration of tested compound that completely blocked the development of 100% of Haemonchus contortus larvae [12].

Depending on the chemistry family of substituent directly linked to carbon in the aryl, the substituted compounds could be divided into two groups: amino group (Chalc3) and halogen group (Chalc2, Chalc4, Chalc5). Moreover, the

Table 1. Molecular structures of the studied imidazo[1,2-a]pyridinyl-chalcones.

different substituents of the aryls (Table 1) are electron-withdrawing by inductive effect or electron-donating by mesomer effect. Furthermore, according to the difference in electronegativity between two fragments engaged in a chemical bond, the fragment rich in electrons because of its electron-withdrawing power is directed towards an electrons poor site while the other electrons deficient fragment points to a site in excess of electrons. They are, respectively, the nucleophilic fragment and the electrophilic fragment. Thus, a chemical bond is made by a nucleophile and its reaction partner (electrophile). Therefore, the very useful chemical concept of electrophilicity and nucleophilicity help the chemist in this quest to rationalize the electronic aspect of reactivity, selectivity and the effect of substituents in a reaction. Moreover, the different nematicidal potentials reflect the effect of substitution on activity in the IPC series. By implementing the methods of quantum chemistry, this work consists in determining the physicochemical properties that characterize the global and local chemical reactivity of five molecules of the IPC series as well as the influence of the aryl substituent. The aim is to identify, for all molecules, the preferential nucleophilic and electrophilic interaction favourite sites of the IPC nucleus not only to the formation of a covalent bond but also to the protection of the pharmacophore with a view to dimerization in IPC series.

2. Calculation Methods

Several theories have been developed for the study of chemical reactivity. In this work, we have used the indices and descriptors that characterize the reactivity of IPC molecules either isolated (global reactivity: static indices) either in interaction with others molecules (local reactivity: dynamic indices). The concepts of global and local reactivity make it possible to elucidate the main questions of relative reactivity, between molecules or between atoms, of the same molecule. All the reactivity parameters used come from the conceptual DFT [23].

2.1. Static Indices of Overall Reactivity

The calculations related to the global hardness (η) and global electrophilicity index (ω) in the Koopmans approximation using the energies of the frontier molecular orbital HOMO and LUMO [24] [25] according to the expressions below:

μ = E HOMO + E LUMO 2 = χ (ii)

η = E LUMO E HOMO 2 (iii)

ω = μ 2 2 η (iv)

These global reactivity descriptors make it possible to rationalize the possible movements of transfer and/ or acquisition of electrons by the isolated molecules. The nucleophilicity of the molecules was defined by calculating the relative nucleophilicity index N [26] [27]. The nucleophilicity scale used is that proposed by Domingo et al. [28], referenced in relation to the tetracyanoethylene (TCE) molecule and defined by the relation:

N = E HOMO ( IPC ) E HOMO ( TCE ) (v)

By taking into account their different classification scales, the use of global electrophilicity and nucleophilicity indices makes it possible to solve the problem of compounds which can be both good nucleophiles and good electrophiles.

2.2. Dynamic Indices of Reactivity

The determination of the reactivity sites of the different IPC molecules was based on the Fukui indices calculated according to the method of the approximation of finite differences:

f k + = q k ( N + 1 ) q k ( N ) (vi)

f k = q k ( N ) q k ( N 1 ) (vii)

q k ( N ) , q k ( N 1 ) , q k ( N + 1 ) are the atomic net charges of the site kth in the neutral, cationic and anionic forms, respectively, of each IPC molecule. f k + and f k are used to identify the most favourable sites for nucleophilic attack and electrophilic attack, respectively. The electrophilic or nucleophilic character of a reactive site is great when the value of the Fukui function is high. However, a site can be both a good nucleophile and a good electrophile. In this case, the ambiguity is removed by determining the canonical dual descriptor, Δ f k , on the one hand. This canonical reactivity descriptor results from the difference between f k + and f k such as Δ f k = f k + f k [29]. The sign of Δ f k defines the reactivity character of the site. A site of the molecule whose electrophilic power is predominant will give a positive value; in the case of a predominant nucleophilic power, Δ f k will be negative. On the other hand, the local reactivity difference index Rk proposed by Chattaraj et al. [30] is also used. The values of this index make it possible to characterize the nature of the electrophilic or nucleophilic or ambiphilic power of a site within an organic molecule. The sign + or − or ± which precedes the Rk values indicates, respectively, an electrophilic or nucleophilic or ambiphilic character of the considered site k. The character is all the stronger the higher the absolute value of Rk. All these physicochemical parameters constitute, on the whole, solid bases for predicting chemical reactivity in IPC series, both qualitatively and quantitatively. Their determinations can only be made in the fundamental state of each of the molecular species studied, respectively. In this context, all geometry optimization calculations are followed by those of frequencies. The calculations of the electronic structures are carried out with the help of the Gaussian09 software [31] using the DFT method by its functional B3LYP and the functions basis 6-311G(d) implemented in the program. The global reactivity parameters are determined only for the IPC molecules when the local reactivity parameters have concerned both the molecules and their respective ionic states.

3. Results and Discussion

The frequencies calculation which followed the geometric optimization of each molecule IPC made it possible to note the absence of imaginary frequency in all the corresponding Hessian matrices. This confirmed that each geometry of structure obtained, at B3LYP/6-311G (d) level, characterizes the fundamental state of the molecule IPC. All the chemical reactivity parameters, for each molecule, have been determined in this state.

3.1. Study of Global Reactivity

Table 2 presents the energies of the frontier orbital (FO), HOMO and LUMO, as well as the energy gaps Δ E such that Δ E = E LUMO E HOMO of the molecules IPC studied are highlighted.

We observe a particularly demarcation between FOs levels of the Chalc3 compound and those of the other compounds. With reference to the unsubstituted molecule Chalc1, the FOs of the amino-aryl derivative (Chalc3) are all

Table 2. Energy values (in eV) of frontier molecular orbital of the studied imidazo[1,2-a]pyridinyl-chalcones.

higher in energy; those of halo-aryl derivatives (Chalc2, Chalc4 and Chalc5) are all lower. Generally, electron-donor substituent increases in the energy level of HOMO or LUMO unlike electron-withdrawing substituent which decreases in the energy level of these frontier orbitals. Substituents of the different aryl groups all have a double electronic effect, i.e. attractor by inductive effect and donor by mesomer effect. Consequently, in Chalc3, the electron-donor nature of dimethylamine is predominant over its electron-withdrawing nature. This observation is in accordance with the rule in organic chemistry which retains the primacy of the mesomer effect over the inductive effect in the event of a double electronic effect of a substituent. However, in the compounds Chalc2, Chalc4 and Chalc5, this rule is taken in default: the inductive effects of halogens appear to be predominant over their mesomer effects. Taking into account these observations, the substituted molecules can be divided into two groups according to the nature of the electronic effects exerted by the substituents. There is the set of compounds Chalc2, Chalc4 and Chalc5 whose substituents affected compound Chalc1 by an electron-withdrawing effect and the second group consisting of compound Chalc3 which is marked by an electron-donor effect. Moreover, HOMO or LUMO energies constitute precious indicators for comparing electron-donor or electro-acceptor powers, in a series of molecules. Chalc3 displays the higher HOMO energy value (−5.334 eV); it is the best electron-donor. Chalc4, which has the lowest LUMO energy (−2.525 eV), is the best electron-acceptor in the series. In the case of Chalc4, there are two chlorine atoms attached to the aryl group. This suggests that the higher the number of halogens on nucleus B, the stronger the overall electro-acceptor character. In addition, the electronic configurations of the atoms which constitute IPC molecules indicate that they are Lewis bases. They are therefore structurally electron-donor. The influence of the aryl substituent on this naturally electron-donor character was highlighted from the deviations of the HOMO energies. Figure 1 shows the diagram established according to the values of the HOMO energy variations in relation to the unsubstituted compound Chalc1.

The diagram illustrates well the electronic effect of each substituent in the aryl group. On the electron-donating scale, ternary amines are far more electron-donor by resonance than halogens, which are very much less electron-donors. The electron-donating power of the studied molecules is then strong when the aryl

Figure 1. Evolution diagram of HOMO energy of IPC.

substituent is a good electron-donor but very weak when the substituent is a strong electron-withdrawing agent. For example, in the case of halo-aryl compounds (Chalc2, Chalc4, Chalc5), the negative values Δ E HOMO clearly reflect the strong attenuation of the electron-donating character due to the presence of halogens and the tendency of the molecules rather towards the global electro-accepting character. Under these conditions, the compounds Chalc2, Chalc4 and Chalc5 become more susceptible to attack a nucleophile. The compound Chalc3 could play the role of electron-donor in IPC intermolecular interactions. Therefore, if we calculate the differences in energies according to the relationship:

Δ E C h a l c 3 C h a l c i = E LUMO ( C h a l c i ) E HOMO ( C h a l c 3 ) , i = 2 , 4 , 5 (viii)

the sequence of evolution in increasing order of the strength of intermolecular interactions coded IPC 3 i , between Chalc3 and another IPC (Chalci), is presented as follows:

Interaction : IPC 3 1 < IPC 3 5 < IPC 3 2 < IPC 3 4 (ix)

The strongest predicted interaction IPC 3 4 , is associated with an energy difference of 2.809 eV with the substituted molecule Chalc4 and that of lesser force, IPC 3 1 , with an energy difference of 3.069 eV implies the unsubstituted molecule Chalc1. In fact, polar interaction is all the stronger when the energy difference between nucleophile and electrophile is small. In the case of molecules in the IPC series, any substitution carried out at the level of the aryl group is always favourable to the strengthening of the nucleophilic-electrophilic interactions. The nucleophilic and electrophilic characters of the substituted compounds are then rationalized using the global nucleophilicity and electrophilicity indices. The values of electronegativity χ, of hardness η and of the global electrophilicity ω and nucleophilicity indices N are reported in Table 3.

According to the values in Table 3, two evolution trends can be defined by following the substituted compounds with reference to the unsubstituted compound Chalc1. Either the characterized chemical property is improved or it is degenerated under the substitution effect. Thus, the lowest values displayed by Chalc3 in the cases of electronegativity and global hardness show that the substituent dimethylamine not only decreased in the ability of the molecule to attract electrons, but also made it more flexible to the deformation of its electronic cloud in an electric field. This is all the more true because the compound Chalc3, identified as the best electron-donor, has the highest global nucleophilicity index (N = 3.879 eV) and the lowest global electrophilicity index (ω = 3.792 eV). In the

Table 3. Values of the global reactivity indices of the studied imidazo[1,2-a]pyridinyl-chalcone.

same way, the compound Chalc4, best electro-acceptor, displays the higher values of the global electrophilicity (ω = 5.178 eV) and electronegativity (χ = 4.363) indicating a great susceptibility of the latter to attract surrounding electrons. According to the general analysis, the halo-aryl derivatives have their values of electronegativity and of global electrophilicity higher than those of the unsubstituted compound Chalc1. Furthermore, their global nucleophilicity values all lower than those of Chalc1, confirm the electron-withdrawing effect exerted by the halogens. The sequences of evolution of the indices in the ascending order of the electrophilic and nucleophilic powers of the compounds are established, respectively, as follows:

ω : C h a l c 3 < C h a l c 1 < C h a l c 5 < C h a l c 2 < C h a l c 4 (x)

N : C h a l c 4 < C h a l c 5 < C h a l c 2 < C h a l c 1 < C h a l c 3 (xi)

Analysis of these evolution sequences shows that they are related to the LUMO and HOMO energies evolution orders, respectively. In fact, the electrophilic character becomes stronger and stronger when the LUMO energy is weaker and weaker, unlike the strength of the nucleophilic character which increases with the HOMO energy. This result is very interesting, especially since the energies of the frontier orbital are regularly associated with the electron-donating and electro-accepting powers of the compounds. Thus, for the IPC molecules, the energies of the LUMO and HOMO are presented as good descriptors of their electrophilic and nucleophilic powers, respectively. The compounds IPC are therefore ambiphilic. This character indicates that they aren’t only nucleophilic by nature but they have a non-negligible electrophilic character also, especially under the effect of substitution with halogens. In fact, apart from the compound Chalc3, the electrophilicity indices of the compounds Chalc2, Chalc4 and Chalc5 are all higher than their respective nucleophilicity indices. For these derivatives, the difference between these indices, which ranges from 1.507 eV (Chalc5) to 2.166 eV (Chalc4), is very significant whereas it is only 0.077 eV with Chalc3. The electrophilic character of compounds Chalc2, Chalc4 and Chalc5 is predominant over their nucleophilic character. The study of reactivity by the global approach of IPCs has made it possible to highlight the relative reactivity of molecules. However, although the global reactivity criteria provide information on the behaviour of a chemical species in its entity, the local reactivity criteria remain very important for the study of the reactivity and selectivity trends of specific sites within a molecule. The local reactivity approach therefore makes it possible to clearly identify the privileged atoms by which a molecule will establish the bonds during molecular interactions. In the case of this study, the isodensity maps relating to the frontier orbitals (HOMO, LUMO) of the different molecules clearly indicate that the HOMO is located mainly on the imidazopyridine nucleus when the LUMO is largely around of the keto-ethylenic system. This therefore supposes that the probable transfer of the electron density takes place from the atoms of the imidazopyridine nucleus to the atoms which constitute the keto-ethylenic system.

3.2. Study of Local Reactivity

In this study, the interest was focused on the heavy atoms of the basic skeleton of imidazo[1,2-a]pyridinyl-chalcone derivatives. Figure 2 shows the numbering of the atoms in the molecules considered.

The analysis of the isodensity maps of HOMO and LUMO, respectively, allowed identifying the potential interaction sites in the compounds studied. From one molecule to another, the atoms which contribute highly to frontier orbital formation stay the same. In fact, the largest lobes characterize the highest electron density in the case of HOMO and the lowest electron density in the case of LUMO, around the atom implicated. On the maps of HOMOs, the four largest lobes contain the carbon atoms C2, C6, C7, C9, C5 and the nitrogen N3. All of these atoms are the likely nucleophile sites in IPCs. On the maps of LUMOs, the carbon atoms C11, C13 and C14 which belong to the two largest lobes are the likely electrophile sites. Figure 3, related to the compound Chalc1, allows illustrating these observations.

Figure 2. Numbering of the atoms in the basic structural skeleton of imidazo[1,2-a]pyridinyl-chalcones.

Figure 3. Isodensity maps of HOMO and LUMO of the compound Chalc1.

Regarding the site C11, it belongs to the carbonyl bond, the presence of which is very useful for the manifestation or even the preservation of the biological properties of chalcones. Thus, the different atoms whose reactivity has been demonstrated are the atoms C2, C7, C8, C9, C5, C13, C14 and N3 of the imidazo[1,2-a]pyridinyl-chalcone nucleus.

The local reactivity indices of all of the interaction sites are determined according to Fukui theory. Their tendency to behave like an electrophile or a nucleophile will be estimated by following the nature of electrons movement caused by the aryl substituents. In Table 4, the different values of the dual index Δ f k ,

Table 4. Local reactivity indices of the interaction sites within the molecules of imidazo[1,2-a]pyridinyl-chalcones calculated at B3LYP/6-311G (d) level.

of the Fukui functions f k + and f k and of the local reactivity difference index Rk obtained for each of these atoms within the molecules are reported. The listed atoms are ordered (from left to right) according to the decrease in the values of the dual index Δ f k so considered either in a decreasing order of their nucleophilic reactivity or in an increasing order of their electrophilic reactivity, within each molecule considered.

The dual reactivity index Δ f k = f k + f k is one of the very precious and useful parameters for rationalizing the sensitivity of a site within a given molecule to an electrophilic or nucleophilic attack, especially when this site displays a double tendency of reactivity. Also, based on the indices analysis, Δ f k and Rk, the sign provides information on the charge inequalities between nucleus and electronic cloud and, consequently, on the predominant reactivity of site k within the molecule. According to Table 4, and considering all five molecules of the series, the higher Δ f k values are observed for the carbon atom C14 followed closely by the nitrogen atom N3 when the lower values are found with the atoms of carbon C2 and C5, within each molecule. Concerning the atom C14, the Δ f k values are all positive, thus showing that this atom has a great susceptibility to increase its electronic density in the event of electrons gained by the molecule. C14 is then an electron deficient site and its reactivity is therefore very strong towards a nucleophile. Moreover, the Fukui function f k + which makes it possible to highlight the electrophilic reactivity of a site relatively to a nucleophilic attack displays its highest values with the atom C14 whatever the molecule. Regarding the atom C5, the Δ f k negative values reflect its susceptibility to have a high electron density even if the molecule loses electrons; it is then rich in electrons and therefore the most reactive centre towards an electrophile. Also, the higher values of the Fukui function f k recorded at this site, whatever the substitution, corroborate this observation. Consequently, the carbon atoms C14 and C5 constitute, in the order, the most reactive electrophilic and nucleophilic centres, in the IPC series. All the sites k with positive values of Δ f k and Rk indicate the potentially electrophilic centres, while the negative values characterize the potentially nucleophilic centres, within the molecules. Thus, analysing Δ f k and Rk values, in the decreasing order of electrophilic reactivity, the site C2 directly follows the site C14 when the site N3 occupies the second place of nucleophilic reactivity after the site C5, in the decreasing order, on all the molecules. Indeed, and particularly within the compound Chalc3, although | Δ f C 13 | = 0.091 > | Δ f C 7 | = 0.048 indicates a higher nucleophilic reactivity of C13 relatively to C7, the value R k = ± 0.249 eV for atom C13 rather reflects a no less stable nucleophilic behaviour. In fact, the value R k = ± 0.249 eV is characteristic of an ambiphilic site, i.e. which behaves as both nucleophile and electrophile. This result agrees with the local chemical reactivity in α, β-unsaturated keto-ethylenic systems where the localization of π electrons due to conjugation keep the α-position carbon C13 in the same ethylenic environment; this carbon atom doesn’t lose or gain electrons during the movement of electrons. Another notable remark concerns the Rk values for site C8. Another remarkable observation is the ± sign which precedes all the Rk values for C8 site, on all of the studied molecules. The presence of this sign reveals that the carbon C8 is an ambiphilic site within the compounds IPC. In addition, the inductive and resonant effects exerted by the substituents attached to the aryl nuclei can possibly influence the sensitivity to attack of the predicted privileged sites C5 and C14. The sense of influence on the reactivity of the sites is then highlighted for the best prediction of intermolecular interactions in the IPCs series. The comparative study was carried out based on the local reactivity difference index Rk. The higher the absolute value of the index Rk of site k, the greater the activation of the site nucleophilicity. At the C5 site, the evolution sequence observed in ascending order of the absolute values of this index is presented below:

| R C 5 | : C h a l c 3 < C h a l c 2 < C h a l c 5 < C h a l c 4 (xii)

Absolute indices of the halo-aryl compounds IPC are all higher than that of the amino-aryl compound IPC. Thereby, the nucleophilic reactivity of the site C5 within the molecules IPC is higher in the presence of halogens. This sequence of evolution is opposite to that of global nucleophilicity for these same compounds. The nucleophilic character of the site C5 is more and more strengthened when the compound IPC becomes less and less nucleophilic. The high nucleophilicity of the site C5 within a molecule is therefore not directly linked to the high nucleophilicity of the molecule. The best electron-donor compound in the series, Chalc3, has the weak electron-donor nucleophilic site C5 while the best electro-acceptor compound in the series, Chalc4, has the best electron-donor nucleophilic site C5. Consequently, substitution with a group having a strong electron-donating character proves unfavourable to the nucleophilic reactivity of site C5. For example, in the case of Chalc3, the substituent group of dimethylamine (strongly electron-donor) caused the nucleophilicity of the site to regress by around 43.48% compared to 1.80% in Chalc2 with fluorine (weakly electron-donor).

As regards the electrophilic site C14, considering the index Rk values, the increasing order of electrophilic reactivity within the different molecules is:

| R C 14 | : C h a l c 5 < C h a l c 4 < C h a l c 2 < C h a l c 3 (xiii)

This sequence of evolution for the electrophilic reactivity of the site C14 within the compounds IPC is not in agreement with the evolution of the global electrophilic behaviour of the studied molecules. However, it seems to correspond to the nucleophilic behaviour evolution of these compounds. The highest electrophilic reactivity of the site C14 is recorded with the best nucleophilic compound Chalc3 in the series when the weakest is found with Chalc5, a less nucleophilic compound. Consequently, substitution on the aryl nucleus with a strongly electron-donating fragment is favourable to strengthening the ability of the electrophilic site C14 to bind to a nucleophilic site. For example, the improving effect of the electrophilic reactivity of C14 is around 39.49% with dimethylamine compared to 9.58% with fluorine, by substitution.

Furthermore, a structure-activity study was undertaken according to the local reactivity descriptors of the identified preferential sites C5 and C14 and the nematicidal potentials of the molecules. In general, no clearly established relationship could be identified with the series of five molecules. However, with the subgroup of halo-aryl derivatives, (Chalc2, Chalc4 and Chalc5), structure-activity relationships are revealed. The sites C5 and C14 evolution sequences were therefore analysed relatively to that of nematicidal potentials of the compounds. Higher nematicidal potential corresponds to the more biologically active molecule against Haemonchus contortus larvae. The nematicidal activity sequence of evolution following the increasing values of nematicidal potentials pLC100 can be defined as below:

p L C 100 : C h a l c 2 < C h a l c 5 < C h a l c 4 (xiv)

The above sequence is opposite to that of the global nucleophilicity N in the subgroup of halogenated derivatives of IPC. In other words, in the series of halogenated IPCs, a weaker global nucleophilic character leads to better nematicidal molecule. Regarding the sites reactivity within the halogenated derivatives, the sequences of evolution according to the increasing values of local difference reactivity index Rk and Fukui positive function f k + are presented, respectively:

| R C 5 | : C h a l c 2 < C h a l c 5 < C h a l c 4 (xv)

f C 14 + : C h a l c 4 < C h a l c 5 < C h a l c 2 (xvi)

The evolution sequences (xv) and (xvi) are opposite to each other. However, the first is similar to that of nematicidal activity. In this case, the nematicidal activity increases with the sensitivity of the site C5 towards an electrophile site. In the second case, more the halogen strengthens the electrophilicity of the site C14, less the nematicidal activity of the molecule is strong. In other words, the nematicidal activity of IPCs becomes strong as the sensitivity of the site C14 to a nucleophile becomes weaker and weaker. Therefore, under the substitution effect, the strength of the electrophilic nature of carbon atom C14 appears as a modulating parameter of the nematicidal activity in the IPCs series.

Moreover, the substitution effects on the nucleophilicity of nitrogen atom N3 have been demonstrated. Previous studies reported that atom N3 was the preferential site for protonation as well as for hydrogen bonding in the IPC series [14] [15]. In the present study, this nitrogen atom appears as the only heteroatom among the potential sites of reactivity identified. And, it is the second best nucleophilic centre of reactivity within the studied substituted IPCs. Its nucleophilic character is therefore of particular importance for the intermolecular interactions involving IPC molecules. In the increasing order of Fukui negative function f k values, the sequence of evolution for the nucleophilic reactivity of nitrogen N3 is the following:

f N 3 : C h a l c 3 < C h a l c 2 < C h a l c 4 < C h a l c 5 (xvii)

Chalc3 is the IPC derivative with dimethylamino-aryl while the compounds Chalc2, Chalc4 and Chalc5 belong to the halo-aryl group. Thereby, the above sequence indicates that halogens better promote the nucleophilic character of nitrogen in the IPCs series by the substitution. Nevertheless, the sequence (xvii) of evolution for the nucleophilic character of the site N3 is in the opposite direction to that of the electrophilic character of the site C14 (xiii) defined following the index Rk values. This fact means that any substitution which reinforces the nucleophilic character of nitrogen N3 disadvantages the electrophilic character of the carbon C14. Such a phenomenon is therefore beneficial to the improvement of the nematicidal activity in the IPC series.

All of these observations are very interesting. They could help in the implementation of strategies for improving nematicidal properties in the IPC series.

4. Conclusion

The chemical reactivity of the compounds in imidazo[1,2-a]pyridinyl-chalcone series is predictable by quantum chemistry methods using conceptual DFT reactivity descriptors. In this series, the effects of substitution on the determined parameters of chemical reactivity are notable especially at the molecular level. On the one hand, the determination of the global descriptors made it possible to show that the presence of electron-donor substituent reinforces the nucleophilic nature of the imidazo[1,2-a]pyridinyl-chalcone compounds. The molecule 3-[4-(dimethylamino)phenyl]-1-(2-methylimidazo[1,2-a]pyridin-3-yl)prop-2-en-1-one has been identified as the most reactive towards electrophilic molecule among the five molecules in the series. On the other hand, the presence of electron-withdrawing substituents such as halogens makes the imidazo[1,2-a] pyridinyl-chalcone compounds more sensitive to nucleophilic molecule attacks. The molecule 3-(2,4-dichlorophenyl)-1-(2-methylimidazo[1,2-a]pyridin-3-yl)prop-2-en-1-one was identified as the most reactive towards any nucleophilic molecule. However, based upon nucleophilicity and electrophilicity index scales, calculations predicted an ambiphilic behaviour of molecules in imidazo[1,2-a]pyridinyl-chalcone series. Thus, these compounds can participate in polar reactions as well as electrophile and nucleophile. From a local point of view, the Fukui indices have made it possible to identify the two preferential nucleophilic and electrophilic centres of interaction for these molecules. The carbon C5 of the imidazo[1,2-a]pyridine nucleus which is directly linked to the carbonyl carbon has been identified as the more potential nucleophilic centre. The carbon C14, in β-position in the keto-ethylenic system was designated as the more potential electrophilic centre. From one molecule IPC to another, the electrophilic or nucleophilic natures as well as the order position of these two preferential sites were independent of the substituents nature and electronic effects. Therefore, in attractive intermolecular interactions in the IPC series, the most likely bond will be established between carbon atoms C14 and C5, regardless of the molecules. This observation was of capital importance in the sense that the formation of dimers based on the IPCs obeyed the same mode of connection whatever the substitution on the aryl group. Finally, the structure-activity analysis revealed that the inactivation of the electrophilic character of the site C14 by substitution on the B nucleus or by formation of a bond with this site would be favourable for improving the IPC series. The results of this study provide an interesting basis for the study of dimerization in the IPC series.

Cite this paper: Konate, B. , Affi, S. and Ziao, N. (2021) DFT Study of Dimerization Sites in Imidazo[1,2-a]pyridinyl-chalcone Series. Computational Chemistry, 9, 1-17. doi: 10.4236/cc.2021.91001.

[1]   Chandrabose, K., Narayana, S.H., Narayana, M., Sakthivel, R., Uma Vanam, Elangovan, M., Devarajan, K. and Piyush, T. (2015) Advances in Chalcones with Anticancer Activities. Recent Patents on Anti-Cancer Drug Discovery, 10, 97-115.

[2]   Dimmock, J.R., Elias, D.W., Beazely, M.A. and Kandepu, N.M. (1999) Bioactivities of Chalcones. Current Medicinal Chemistry, 6, 1125-1149.

[3]   Batovska, D.I. and Todorova, I.T. (2010) Trends in Utilization of the Pharmacological Potential of Chalcones. Current Clinical Pharmacology, 5, 1-29.

[4]   Sahu, N.K., Balbhadra, S.S., Choudhary, J. and Kohli, D.V. (2012) Exploring Pharmacological Significance of Chalcone Scaffold: A Review. Current Medicinal Chemistry, 19, 209-225.

[5]   Dhar, D.N. (1981) The Chemistry of Chalcones and Related Compounds. Wiley, New York.

[6]   Das, B.P., Choudhury, T.R., Das, G.K., Choudhury, B. and Chowdhury, D.N. (1995) Comparative Larvicidal Activity of Some Chalcones on the Larvae of Culex Quinquefasciatus. Environment & Ecology, 13, 694-697.

[7]   Lahtchev, K.L., Batovska, D.I., Parushev, S.P., Ubiyvovk, V.M. and Sibirny, A.A. (2008) Antifungal Activity of Chalcones: A Mechanistic Study Using Various Yeast Strains. European Journal of Medicinal Chemistry, 43, 2220-2228.

[8]   Bag, S., Ramar, S. and Degani, M.S. (2009) Synthesis and Biological Evaluation of α, β-Unsaturated Ketone as Potential Antifungal Agents. Medicinal Chemistry Research, 18, 309-316.

[9]   Awasthi, S.K., Mishra, N., Dixit, S.K., Singh, A., Yadav, M., Yadav, S.S. and Rathaur, S. (2009) Antifilarial Activity of 1,3-Diarylpropen-1-one: Effect on Glutathione-S- Transferase, A Phase-II Detoxification Enzyme. The American Journal of Tropical Medicine and Hygiene, 80, 764-768.

[10]   Sissouma, D., Ouattara, M., Koné, M.W., Menan, H.E., Adjou, A. and Ouattara, L. (2011) Synthesis and in Vitro Nematicidal Activity of New Chalcones Vectorised by Imidazopyridine. Africa Journal of Pharmacy and Pharmacology, 5, 2086-2093.

[11]   Bagdi, A.K., Santra, S., Monir, K. and Hajra, A. (2015) Synthesis of Imidazo[1,2-a] pyridines: A Decade Update. Chemical Communications, 51, 1555-1575.

[12]   Kuthyala, S., Nagaraja, G.K., Sheik, S., Hanumanthappa, M. and Kumar, M.S. (2019) Synthesis of Imidazo[1, 2-a]pyridine-chalcones as Potent Inhibitors against A549 Cell Line and Their Crystal Studies, Journal of Molecular Structure, 1177, 381-390.

[13]   Rao, N.S., Kistareddy, C., Balram, B. and Ram, B. (2012) Synthesis and Antibacterial Activity of Novel Imidazo[1,2-a]pyrimidine and Imidazo[1,2-a]pyridine Chalcones Derivatives. Der Pharma Chemica, 4, 2408-2415.

[14]   Thomas, A.S., Nahossé, Z. and Kafoumba, B. (2015) Determination, par des méthodes ab initio et DFT, des sites et énergies de protonation d’une série de molécules d’imidazopyridinyl-chalcones substituées. European Scientific Journal, 11, 148-158.

[15]   Affi, T.S., Ziao, N., Ouattara, M., Sissouma, D. and Yapo, K. (2014) Caractérisation théorique des sites d’interaction par liaison hydrogène de 3-(4-isopropylphenyl)-1- (2-méthylimidazopyridin-3-yl)prop-2-en-1-one et de 3-(2-méthoxyphenyl)-1-(2- méthylimidazopyridin-3-yl)prop-2-en-1-one. European Journal of Scientific Research, 123, 340-347.

[16]   Gadkari, S., Choudhari, P., Bhatia, M., Khetmar, S. and Jadhav, S. (2012) 3D QSAR, Pharmacophore Identification Studies on Series of Imidazopyridine Analogs as Nematicidal Activity. Pharmacophore, 3, 199-208.

[17]   Aksöz, B.E. and Ertan, R. (2011) Chemical and Structural Properties of Chalcones I. FABAD Journal of Pharmaceutical Sciences, 36, 223-242.

[18]   Parr, R.G., Szentpály, L.V. and Liu, S. (1999) Electrophilicity Index. Journal of the American Chemical Society, 121, 1922-1924.

[19]   Domingo, L.R., Aurell, M.J., Pérez, P. and Contreras, R. (2002) Quantitative Characterization of the Global Electrophilicity Power of Common Diene/Dienophile Pairs in Diels-Alder Reactions. Tetrahedron, 58, 4417-4423.

[20]   Eyring, H. and Polanyi, M. (2013) On Simple Gas Reactions. Zeitschrift für Physikalische Chemie, 227, 1221-1245.

[21]   Eyring, H. (1935) The Activated Complex in Chemical Reactions. The Journal of Chemical Physics, 3, 107-115.

[22]   Geerlings, P., De Proft, F. and Langenaeker, W. (2003) Conceptual Density Functional Theory. Chemical Reviews, 103, 1793-1873.

[23]   Parr, R.G. and Yang, W. (1989) Density Functional Theory of Atoms and Molecules. Oxford University Press, New York.

[24]   Parr, R.G. and Pearson, R.G. (1983) Absolute Hardness: Companion Parameter to Absolute Electronegativity. Journal of the American Chemical Society, 105, 7512-7516.

[25]   Domingo, L.R., Chamorro, E. and Pérez, P. (2008) Understanding the Reactivity of Captodative Ethylenes in Polar Cycloaddition Reactions; a Theoretical Study. The Journal of Organic Chemistry, 73, 4615-4624.

[26]   Pérez, P., Domingo, L.R., Duque-Noreña, M. and Chamorro, E. (2009) A Condensed-to-Atom Nucleophilicity Index. An Application to the Director Effects on the Electrophilic Aromatic Substitutions. Journal of Molecular Structure (Theochem), 895, 86-91.

[27]   Jaramillo, P., Domingo, L.R., Chamorro, E. and Pérez, P. (2008) A Further Exploration of a Nucleophilicity Index Based on the Gas-Phase Ionization Potentials. Journal of Molecular Structure (Theochem), 865, 68-72.

[28]   Domingo, L.R., Aurell, M.J., Pérez, P. and Contreras, R. (2002) Quantitative Characterization of the Local Electrophilicity of Organic Molecules. Understanding the Regioselectivity on Diels-Alder Reactions. The Journal of Physical Chemistry A, 106, 6871-6875.

[29]   Morell, C., Grand, A. and Toro-Labbé, A. (2005) New Dual Descriptor for Chemical Reactivity. The Journal of Physical Chemistry A, 109, 205-212.

[30]   Chattaraj, P.K., Duleya, S. and Domingo, L.R. (2012) Understanding Local Electrophilicity/Nucleophilicity Activation through a Single Reactivity Difference Index. Organic and Biomolecular Chemistry, 10, 2855-2861.

[31]   Frisch, M.J., Trucks, G.W., Schlegel, H.B., Scuseria, G.E., Robb, M.A., Cheeseman, J.R., Scalmani, G., Barone, V., Mennucci, B., Petersson, G.A., Nakatsuji, H., Caricato, M., Li, X., Hratchian, H.P., Izmaylov, A.F., Bloino, J., Zheng, G., Sonnenberg, J.L., Hada, M., Ehara, M., Toyota, K., Fukuda, R., Hasegawa, J., Ishida, M., Nakajima, T., Honda, Y., Kitao, O., Nakai, H., Vreven, T., Montgomery Jr., J.A., Peralta, J.E., Ogliaro, F., Bearpark, M., Heyd, J.J., Brothers, E., Kudin, K.N., Staroverov, V.N., Kobayashi, R., Normand, J., Raghavachari, K., Rendell, A., Burant, J.C., Iyengar, S.S., Tomasi, J., Cossi, M., Rega, N., Millam, J.M., Klene, M., Knox, J.E., Cross, J.B., Bakken, V., Adamo, C., Jaramillo, J., Gomperts, R., Stratmann, R.E., Yazyev, O., Austin, A.J., Cammi, R., Pomelli, C., Ochterski, J.W., Martin, R.L., Morokuma, K., Zakrzewski, V.G., Voth, G.A., Salvador, P., Dannenberg, J.J., Dapprich, S., Daniels, A.D., Farkas, O., Foresman, J.B., Ortiz, J.V., Cioslowski, J. and Fox, D.J. (2009) Gaussian 09, Revision A.1. Gaussian Inc., Wallingford.