# Synthesis, crystal structure, DFT studies and biological activity of (Z)-3-(3-bromophenyl)-1-(1,5-dimethyl-1H-pyrazol-3-yl)-3-hydroxyprop-2-en-1-one

Said Tighadouini, Redouane Benabbes, Monique Tillard, Driss Eddike, Khadija Haboubi, Khalid Karrouchi, Smaail Radi

<figure><img src="image_1.png" /></figure>

## To cite this version:

Said Tighadouini, Redouane Benabbes, Monique Tillard, Driss Eddike, Khadija Haboubi, et al. Synthesis, crystal structure, DFT studies and biological activity of (Z)-3-(3-bromophenyl)-1-(1,5-dimethyl-1H-pyrazol-3-yl)-3-hydroxyprop-2-en-1-one. Chemistry Central Journal, 2018, 12 (1), pp.122. 10.1186/s13065-018-0492-4. http://hal-01934919

HAL Id: hal-01934919

https://hal.science/hal-01934919

Submitted on 31 May 2021

HAL is a multi-disciplinary open access archive for the deposit and dissemination of scientific research documents, whether they are published or not. The documents may come from teaching and research institutions in France or abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est destinée au dépôt et à la diffusion de documents scientifiques de niveau recherche, publiés ou non, émanant des établissements d'enseignement et de recherche français ou étrangers, des laboratoires publics ou privés.

<figure><img src="image_2.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
RESEARCH ARTICLE

Open Access

<figure><img src="image_3.png" /></figure>

CrossMark

# Synthesis, crystal structure, DFT studies and biological activity of (Z)-3-(3-bromophenyl)-1-(1,5-dimethyl-1H-pyrazol-3-yl)-3-hydroxyprop-2-en-1-one

Said Tighadouini¹, Redouane Benabbes², Monique Tillard<fnref n="3" />*, Driss Eddike⁴, Khadija Haboubi⁵, Khalid Karrouchi<fnref n="6" />* and Smaail Radi¹

**Abstract**

**Background:** Nowadays, is emerging a new generation of highly promising inhibitors bearing the β-ketoenol functionality. The present work relates to the first synthesis, the structure determination, the DFT studies and the use of a new biomolecule designed with a β-ketoenol group bounded to a pyrazolic moiety.

**Result:** A novel β-ketoenol-pyrazole has been synthesized, well characterized and its structure was confirmed by single crystal X-ray diffraction. The electron densities and the HOMO–LUMO gap have been calculated using the DFT method with BLYP, PW91, PWC functionals and 6-31G* basis set. An evaluation of the molecule stability is provided by a NBO analysis and the calculated Fukui and Parr functions have been used to locate the reactive electrophile and nucleophile centers in the molecule. The synthesized compound, screened for its in vitro antifungal behavior against the *Fusarium oxysporum* f.sp. albedinis FAO fungal strains, shows a moderate activity with an inhibition percentage of 46%. The product was also tested against three bacterial strains (*Escherichia coli*, *Bacillus subtilis* and *Micrococcus luteus*), but no significant effect was observed against these organisms.

**Conclusions:** Density functional calculations are used to evaluate the HOMO–LUMO energy gap, the molecular electrostatic potential and to provide a natural bond orbital analysis. The measured antimicrobial activities encourage us to continue searching for other structures, likely to be good antifungal candidates.

**Keywords:** β-Keto-enol-pyrazole, Single-crystal structure, NBO analysis, Reactivity indices, Fukui and Parr functions, Biological activity

**Introduction**

Pyrazoles represent a class of compounds endowed with a great interest in many domains. They have been widely described in the literature as chelating ligands [1–6] and several works have been gathered in reviews [7–10]. According to numerous literature reports, these derivatives are also well-known as important

heterocyclic biologically active compounds, acting as antitumor [11], antiviral [12], anti-inflammatory [13] anti-anxiety [14] or antimicrobial [15] agents.

On the other hand, β-ketoenols form an important class of compounds, with an interest both in medical and pharmaceutical fields, regarded as drugs against HIV [16–18], cancer [19–22] and influenza [23] but also as antioxidant [24] and anti-inflammatory [25] substances. The β-ketoenol derivatives play also an important role in the development of coordination chemistry, as they are able to easily form stable complexes with most transition metals involving different

<footnote>*Correspondence: mtillard@univ-montp2.fr; khalid.karrouchi@um5s.net. ma</footnote>

<footnote n="3">³ ICGM, CNRS, Université de Montpellier, ENSCM, Montpellier, France</footnote>

<footnote n="6">⁶ Laboratoire de Chimie Thérapeutique, Faculté de médecine et de Pharmacie, Université Mohamed V, Rabat, Morocco</footnote>

<footnote>Full list of author information is available at the end of the article</footnote>
modes of coordination and different functionalities [26, 27].

The pyrazoles associated with β-ketoenol groups lead to compounds with promising properties in both medicinal and coordination chemistry fields. In our recent works, some heterocycles containing the β-ketoenol functionality have been reported, that show significant biological activity [28] as well as interesting coordination properties [29–33].

The intention of this work was to develop of a novel pyrazole-based compound bearing a β-ketoenol functionality. Its crystal structure was solved from X-ray single crystal data and DFT studies were realized. The compound was also evaluated for its in vitro antifungal activity against *Fusarium oxysporum* f.sp. albedinis FAO fungal strains and against three bacterial strains (*Escherichia coli*, *Bacillus subtilis* and *Micrococcus luteus*).

# Results and discussion

## Chemistry

The target biomolecule based on β-ketoenol and pyrazole entities was prepared by a one-pot in situ condensation method which is similar to the procedures given in our previous works [28]. A solution of pyrazolic carboxylate was added to a suspension of sodium in toluene, then 1-(3-bromophenyl)ethanone was added at 0 °C (Scheme 1). After 2-days stirring at room temperature, the resulting precipitate has been treated and neutralized. The extracted organic layer was concentrated, dried and purified by silica gel column chromatography (see "Experimental section" part for details).

The β-keto-enol form was confirmed by the ¹H-NMR analysis of the compound whose spectrum (Additional file 1: Figure S1) shows a strong signal assigned to the =C–H group of the keto-enol form at 6.54 ppm, it represents 85% of the compound. The diketone form is also present in a maximal proportion of 15% and was detected by the weak signal at 4.54 ppm which was attributed to the CH₂ group of the diketone form. Traces of the keto form have also been detected in DEPTQ-135, which shows quaternary carbon atoms (C) and CH₂ group as

very small negative signals (Additional file 1: Figures S2 and S3). Good quality crystals of the major β-ketoenol structure were grown from methanolic solution by slow evaporation. The FT-IR spectrum confirms the formation of the ketoenol form with an enolic band at 1531 cm⁻¹ (Additional file 1: Figure S4). Also in good agreement, the mass spectrum shows a molecular peak at 320.97. (Additional file 1: Figure S5).

## X-ray crystal structure description

Single crystals of (Z)-3-(3-bromophenyl)-1-(1,5-dimethyl-1H-pyrazol-3-yl)-3-hydroxyprop-2-en-1-one (1) were analyzed by X-ray diffraction in order to determine the compound structure.

The main crystal data are given with principal refinement parameters in Table 1 and the atom position and displacement parameters are listed in Table 2. The full CIF file deposited at the Cambridge Crystallographic Data Center (CCDC 1817604) is available at http://www.ccdc.cam.ac.uk/conts/retrieving.html.

The asymmetric unit contains two independent molecules. Each molecule is built with two rings, a bromophenyl ring bonded (at meta position) to a dimethyl pyrazole ring through a central core unit –CO–C–COH– (Fig. 1). According to the root mean square deviations of the fitted atoms in each group, ranging from 0.003 to 0.008, these three units are planar. The dihedral angles between the central core and the two rings, between 4.3 and 8.7°, indicate a slight deviation to flatness within each of the independent molecules. The two independent molecular units are almost coplanar, as shown by the angle of only 1.96(3)° measured between their mean planes (calculated with all the non-H atoms). The bonds lengths and angles measured in the two molecules are very close and in range of values found in the literature for similar compounds [34–38].

The molecular arrangement in the solid is such as one unit cell contains four molecules, two-by-two symmetry-related. As a consequence of the co-planarity of the independent molecular units, the packing in the crystal results in a layered arrangement shown in Fig. 2. Within planes parallel to (111) and separated by ~3.5 Å, each

<figure><img src="image_4.png" /><figcaption>**Scheme 1** Synthesis of the target compound 1</figcaption></figure>
molecule is surrounded by six homologous units, this molecular organization is studied later in this work for evaluation of in-plane interactions. On the other hand, no π-stacking interaction must be expected because of relative position of the successive planes without ring overlapping.

### DFT calculations

To investigate the molecular geometry and the electron distribution in the solid, density functional theory (DFT)

<figure><img src="image_5.png" /><figcaption>**Fig. 1** Ortep molecular representation of C₁₄ H₁₃ Br N₂ O₂ (30% probability ellipsoids)</figcaption></figure>

<table><caption>**Table 1** Crystal data and structure refinement for C₁₄ H₁₃ Br N₂ O₂</caption><tr><td>CCDC deposit number</td><td>1817604</td></tr><tr><td>Formula, M, Z</td><td>C₁₄ H₁₃ Br N₂ O₂, 321.17, 4</td></tr><tr><td>Space group</td><td>Triclinic, P<span style="text-decoration: overline">1</span></td></tr><tr><td>Lattice</td><td>a = 11.1458(7), b = 11.6337(3), c = 12.7221(9) Å, α = 112.075(2), β = 105.637(2), γ = 103.793(2)°</td></tr><tr><td>θ range</td><td>1.88 to 29.09°</td></tr><tr><td>Reflections</td><td>18,366 collected/6297 unique [R(int) = 0.0414]</td></tr><tr><td>Crystal</td><td>colorless, 0.20 × 0.13 × 0.12 mm</td></tr><tr><td>Data/parameters</td><td>6297/357</td></tr><tr><td>R indices [I > 2σ(I)]</td><td>R1 = 0.0451, wR2 = 0.0968</td></tr><tr><td>R indices (all data)</td><td>R1 = 0.1012, wR2 = 0.1068</td></tr><tr><td>Δρ Fourier residuals</td><td>0.55/−0.33 e.Å<sup>−3</sup></td></tr></table>

calculations were carried out using the program Dmol³ at three DFT levels: with PW91 or BLYP functional within the GGA (generalized gradient approximation) and with PWC functional within the LDA (local density approximation) [39–41]. Double numerical plus polarization DNP basis sets were taken in all calculations.

Full geometry optimization (LDA-PWC) by minimization of the total energy was first carried out, both starting from experimental geometry and from three planar geometries built by rotation of the rings around the C1–C9 and C3–C4 bonds (Fig. 3).

The total energies calculated indicate a lower stability, by 0.16 and 0.17 eV, for the M2 and M3 models compared to experimental geometry. This potentially results from

<table><caption>**Table 2** Atomic coordinates (× 10⁴) and equivalent isotropic displacement parameters (Å² × 10³) for C₁₄ H₁₃ Br N₂ O₂. Ueq is defined as 1/3 of the trace of the orthogonalized Uij tensor</caption><thead><tr><th rowspan="2"></th><th colspan="4">Molecule 1</th><th colspan="4">Molecule 2</th></tr><tr><th>x</th><th>y</th><th>z</th><th>Ueq</th><th>x</th><th>y</th><th>z</th><th>Ueq</th></tr></thead><tbody><tr><td>Br</td><td>56(1)</td><td>6869(1)</td><td>227(1)</td><td>66(1)</td><td>1104(1)</td><td>12,097(1)</td><td>−3397(1)</td><td>77(1)</td></tr><tr><td>O1</td><td>−4056(2)</td><td>6644(2)</td><td>−3183(2)</td><td>56(1)</td><td>4307(2)</td><td>11,637(2)</td><td>182(2)</td><td>55(1)</td></tr><tr><td>O2</td><td>−6166(2)</td><td>6129(2)</td><td>−4913(2)</td><td>55(1)</td><td>5549(2)</td><td>11,122(2)</td><td>1789(2)</td><td>57(1)</td></tr><tr><td>N1</td><td>−10,238(2)</td><td>3242(2)</td><td>−5973(2)</td><td>38(1)</td><td>3715(3)</td><td>8399(2)</td><td>3102(2)</td><td>44(1)</td></tr><tr><td>N2</td><td>−8973(2)</td><td>3673(3)</td><td>−5150(2)</td><td>40(1)</td><td>3306(3)</td><td>8867(3)</td><td>2305(2)</td><td>44(1)</td></tr><tr><td>C1</td><td>−4761(3)</td><td>5719(3)</td><td>−3007(3)</td><td>40(1)</td><td>3195(3)</td><td>10,766(3)</td><td>28(3)</td><td>42(1)</td></tr><tr><td>C2</td><td>−6118(3)</td><td>5006(3)</td><td>−3717(3)</td><td>42(1)</td><td>3209(4)</td><td>10,106(3)</td><td>720(3)</td><td>44(1)</td></tr><tr><td>C3</td><td>−6786(3)</td><td>5265(3)</td><td>−4677(3)</td><td>42(1)</td><td>4436(3)</td><td>10,317(3)</td><td>1601(3)</td><td>42(1)</td></tr><tr><td>C4</td><td>−8234(3)</td><td>4514(3)</td><td>−5421(3)</td><td>37(1)</td><td>4448(3)</td><td>9595(3)</td><td>2341(3)</td><td>40(1)</td></tr><tr><td>C5</td><td>−9038(3)</td><td>4600(3)</td><td>−6420(3)</td><td>43(1)</td><td>5559(3)</td><td>9567(3)</td><td>3141(3)</td><td>46(1)</td></tr><tr><td>C6</td><td>−10,317(3)</td><td>3787(3)</td><td>−6759(3)</td><td>39(1)</td><td>5077(3)</td><td>8798(3)</td><td>3630(3)</td><td>39(1)</td></tr><tr><td>C7</td><td>−11,613(3)</td><td>3477(4)</td><td>−7736(3)</td><td>53(1)</td><td>5792(4)</td><td>8413(3)</td><td>4543(3)</td><td>55(1)</td></tr><tr><td>C8</td><td>−11,320(3)</td><td>2284(3)</td><td>−5949(3)</td><td>50(1)</td><td>2700(3)</td><td>7551(4)</td><td>3293(3)</td><td>62(1)</td></tr><tr><td>C9</td><td>−3978(3)</td><td>5493(3)</td><td>−2014(3)</td><td>41(1)</td><td>1964(3)</td><td>10,554(3)</td><td>−957(3)</td><td>42(1)</td></tr><tr><td>C10</td><td>−4596(3)</td><td>4618(3)</td><td>−1630(3)</td><td>53(1)</td><td>700(3)</td><td>9632(3)</td><td>−1272(3)</td><td>51(1)</td></tr><tr><td>C11</td><td>−3823(4)</td><td>4409(4)</td><td>−715(3)</td><td>61(1)</td><td>−430(4)</td><td>9417(4)</td><td>−2231(4)</td><td>61(1)</td></tr><tr><td>C12</td><td>−2436(4)</td><td>5076(4)</td><td>−175(3)</td><td>59(1)</td><td>−292(4)</td><td>10,156(4)</td><td>−2853(3)</td><td>65(1)</td></tr><tr><td>C13</td><td>−1844(3)</td><td>5938(3)</td><td>−548(3)</td><td>47(1)</td><td>950(4)</td><td>11,076(3)</td><td>−2530(3)</td><td>51(1)</td></tr><tr><td>C14</td><td>−2599(3)</td><td>6146(3)</td><td>−1451(3)</td><td>43(1)</td><td>2073(3)</td><td>11,293(3)</td><td>−1593(3)</td><td>46(1)</td></tr></tbody></table>
the occurrence of some N···O repulsive interactions in these configurations. Instead, the M1 model is almost as stable (only differing by 0.001 eV) as the experimental molecule in which the central core and the rings are rigorously coplanar (dihedral angles lower than 2°).

Trying to evaluate intermolecular interactions within the solid, geometry was optimized for a large molecular fragment consisting of a molecule and its six surrounding neighbors. The three dimensional contour of the total density, drawn at the 0.04e⁻/Å³ isolevel and mapped with the electron deformation density, has been represented in Fig. 4. The deformation density, computed as the total electron density with the density of isolated atoms subtracted, points out the electron localization as positive regions and the electron losses as negative regions. Looking at its representation, one can conclude that no strong intermolecular interaction exists in this compound. Nevertheless, the positive domains (yellow zones in Fig. 4) indicate the occurrence of an interaction between the bromine atom of a molecule and the nitrogen N2 atom of a neighboring molecule.

As expected, the periodic calculations in the solid confirm the absence of any bonding density between the molecular planes separated by 3.5 Å. Selected bonds and angles given in Table 3 illustrate the structural, and rather weak, packing constraints onto the molecular geometry. It is interesting to remark that the dihedral angle between mean planes of the two independent molecules in the unit cell is about 4° in the optimized structures, a

value which is nearly twice the angle in the experimental crystal. On the other hand, the experimental torsion angles C4–C3–C2–C1 and C3–C2–C1–C9, which reflect the molecule flatness, range from 176.8 to 179.6° and do not differ much from the values in optimized molecular (177.8–179.9°) and 3D periodic (175.6–179.8°) models.

Within the overall context of the FMO theory, the energetic level and the form of the frontier orbitals are relevant parameters for analysis of the molecular reactivity [42]. Regardless of the DFT level of theory, an HOMO–LUMO energy separation of about 2.3 eV has been calculated with DMol³ for the molecule (2.7 eV in the crystal), this value should be compared with the gap of 3.46 eV measured experimentally from UV experiments (absorption peak at 358 nm) (Additional file 1: Figure S6). Both orbitals display a π-type character mainly localized on the central part of the molecule, which is bonding at the HOMO and antibonding at the LUMO levels.

## Fukui and Parr functional analysis

Calculated electron deformation density is strongly related to molecular electrostatic potential and for this reason may be equally used to discuss the reactivity [43]. The latter is also seen as a useful tool in evaluation of the regiochemistry, especially in reactions that are dominated by electrostatic effects. The electrophilic f(−) and nucleophilic f(+) Fukui functions, whose extreme values reflect the ability for an electrophilic or a nucleophilic attack, are defined as electron density derivatives with

<figure><img src="image_6.png" /><figcaption>**Fig. 2** Molecular packing in the triclinic lattice enhancing the peculiar layered arrangement of molecules in planes parallel to (111)</figcaption></figure>
<figure><img src="image_7.png" /></figure>

Exp

<figure><img src="image_8.png" /></figure>

M1

<figure><img src="image_9.png" /></figure>

M3

<figure><img src="image_10.png" /></figure>

M2

<figure><img src="image_11.png" /><figcaption>**Fig. 4** Representation of the 3D isosurface electron density (volumic contour) mapped with the deformation density. Positive value (yellow domains) indicate electron localization</figcaption></figure>

respect to the number of electrons at a constant potential. They can give a measure of the local reactivity and they have been considered suitable to rationalize the regioselectivity [44–47].

Though, it has been stated that regioselectivity in polar reactions should be predictable alternately using the local electrophilicity [48]. Then electrophilic and nucleophilic Parr functions P<sub>k</sub><sup>+</sup> and P<sub>k</sub><sup>-</sup> are powerful tools to study the molecular reactivity and they are well adapted to localize
<table><caption>Table 3 Selected bond lengths [Å] and angles [°] for C₁₄ H₁₃ Br N₂ O₂</caption><thead><tr><th></th><th>Mole 1/mole 2</th><th>PW91</th><th>BLYP</th><th>PWC</th><th>Solid PWC</th><th>Solid BLYP</th><th>Solid PW91</th></tr></thead><tbody><tr><td>Br–C13</td><td><b>1.900(3)/1.910(3)</b></td><td>1.928</td><td>1.955</td><td>1.899</td><td>1.913–1.916</td><td>1.950–1.957</td><td>1.928–1.935</td></tr><tr><td>O2–C3</td><td><b>1.260(4)/1.262(3)</b></td><td>1.275</td><td>1.276</td><td>1.280</td><td>1.282–1.283</td><td>1.281–1.280</td><td>1.280–1.277</td></tr><tr><td>N1–N2</td><td><b>1.345(3)/1.348(3)</b></td><td>1.344</td><td>1.358</td><td>1.326</td><td>1.328–1.330</td><td>1.353–1.357</td><td>1.341–1.345</td></tr><tr><td>N1–C6</td><td><b>1.365(4)/1.361(4)</b></td><td>1.378</td><td>1.386</td><td>1.365</td><td>1.364–1.365</td><td>1.379–1.381</td><td>1.373–1.375</td></tr><tr><td>N1–C8</td><td><b>1.450(4)/1.452(4)</b></td><td>1.451</td><td>1.464</td><td>1.431</td><td>1.434–1.435</td><td>1.461–1.463</td><td>1.449–1.451</td></tr><tr><td>N2–C4</td><td><b>1.333(4)/1.329(4)</b></td><td>1.353</td><td>1.336</td><td>1.338</td><td>1.342–1.341</td><td>1.356–1.354</td><td>1.351–1.350</td></tr><tr><td>O1–C1</td><td><b>1.304(4)/1.312(4)</b></td><td>1.328</td><td>1.343</td><td>1.297</td><td>1.303–1.301</td><td>1.342–1.342</td><td>1.328–1.328</td></tr><tr><td>N2–N1–C6</td><td><b>112.8(2)/113.5(2)</b></td><td>112.3</td><td>112.7</td><td>111.3</td><td>113.2–113.1</td><td>112.9–112.7</td><td>113.1–113.2</td></tr><tr><td>C4–N2–N1</td><td><b>104.3(2)/103.8(2)</b></td><td>104.7</td><td>104.7</td><td>104.8</td><td>104.9–105.0</td><td>104.6–104.8</td><td>104.8–104.7</td></tr><tr><td>mol1–mol2 angle</td><td><b>1.96</b></td><td></td><td></td><td></td><td>4.4</td><td>3.5</td><td>4.2</td></tr></tbody></table>

Experimental (bolditalic) values in the crystal are compared with values calculated (italic) in molecular or 3D models

the electrophile and nucleophile centers in an organic molecule [49, 50]. They have been computed from the atomic spin density (difference between α and β electron densities) for the radical anion and the radical cation having the geometry optimized for the neutral molecule.

The spatial distribution of the atomic spin density mapped on the electrostatic potential provides a graphical view of the localization of electrophilic and nucleophilic centers (Fig. 5). In the present case, the Parr functions quite well validate the results predicted by Fukui functions with an electrophile center at C2 atom (and at a less extent at N2 and C10 atoms) while the nucleophile center is situated at the central core of the molecule mainly at O2 atom, most likely to undergo a nucleophilic attack.

The global electrophilicity index ω, expressed as μ²/2η [51], takes here the value of 2.28. It has been obtained from the electronic chemical potential μ (μ=1/2 × (E<sub>HOMO</sub>+E<sub>LUMO</sub>) = -4.37 eV) and the chemical hardness η (η = E<sub>LUMO</sub> - E<sub>HOMO</sub> = 4.17 eV), both calculated with G03W code from the one-electron energies of the frontier orbitals. The maximal charge transfer (−μ/η) is found close to unity and the nucleophilicity index N close to 3. The latter is calculated comparatively to TCE taken as a reference from N = E<sub>HOMO</sub> - E<sub>HOMO</sub>(TCE) where

E<sub>HOMO</sub>(TCE) = 9.4083 eV is the energy calculated in the same conditions for tetracyanoethylene. Based on these calculated theoretical reactivity indices, the molecule is characterized with a moderate nucleophile and rather good electrophile character.

## Natural bond orbital (NBO) analysis

The NBO analysis is a helpful way to study the interactions among bonds and to examine the charge transfer resulting of conjugative interactions in a molecular system [52, 53]. The hyperconjugative interaction energy is deduced from a second-order perturbation approach [53]. Considering a donor (i) and an acceptor (j), the stabilization energy E(2) associated with delocalization is estimated as $q_i \times F(i, j)^2 / (\epsilon_i - \epsilon_j)$, where $q_i$ is the donor orbital occupancy, $\epsilon_i$ and $\epsilon_j$ are diagonal elements (orbital energies) and $F(i, j)$ is the off-diagonal NBO Fock matrix element. The larger the stabilization energy, the stronger the donor-to-acceptor interaction, i.e. more important is the electron-donor trend towards acceptor and greater is the extent of conjugation on the whole system. Table 4 summarizes the highest interactions between bonding and antibonding (Lewis/non-Lewis) natural orbitals as for example between the π C1–C2 donor and the π* O2–C3 acceptor with a stabilizing energy of 32.54 kcal.

<figure><img src="image_12.png" /><figcaption>**Fig. 5** Atomic spin density spatial distribution mapped onto the electrostatic potential showing the electrophilic (left image) and nucleophilic (right image) centers as red zones</figcaption></figure>
mol⁻¹. The lone pair (LP) orbitals are also seen to have important stabilizing contributions as illustrated with interaction of the lone pair at O1 atom with the π* C1–C2 acceptor (44.40 kcal.mol⁻¹).

## Biological activity

The in vitro antibacterial and antifungal activities of compound **1** were tested by the agar diffusion technique [54–56] using fungal strains (*Fusarium oxysporum* f.sp. albedinis *FAO*) and bacterial strains (*Escherichia coli*, *Bacillus subtilis*, and *Micrococcus luteus*). Tests were also performed for comparison on several compounds already prepared in our former works (Fig. 6).

Results of the tests carried out towards bacterial strains for the new compound **1** and also for the other products **2–6** revealed no significant effect against these organisms.

On the contrary, the structure **1** led to a moderate inhibiting antifungal activity of 46% occurring at 500 µL of sample. This antifungal effect remains rather modest with regard to the benomyl fungicide (94%) as reported in Table 5.

Generally, the results obtained for this new structure are in agreement with activities of the similar compounds **2–4** [28]. The inhibiting percentage of 46% found for **1** is in the range of values obtained for compounds **2–4**, between 36 and 54%. Instead, the antifungal

performances found for structures **5** and **6**, also belonging to the same family, reach values very close to the benomyl fungicide taken as reference. Such differences in the biological activity are obviously dependent on the radical group attached to the pyrazol ketoenol fragment. The best inhibiting activities are obtained with a methyl or a phenyl group while the presence of substituted phenyl groups (*m*-bromophenyl, *p*-methyl phenyl or *p*-methoxy phenyl) provide lower antifungal properties. These findings push us to dig deeper to find much more interesting molecules. Of course, various structural modifications to bring to these compounds as antifungal candidates are currently in progress.

## Experimental section

### General information

All solvents and other chemicals (purity > 99.5%, Aldrich, Saint-Louis, MO, USA) of analytical grade were used without further purification. An Xcalibur four circle CCD diffractometer (Oxford Diffraction, Abingdon, Oxfordshire, England) was used to collect the X-ray intensities diffracted by a parallelepiped colorless selected crystal (CNRS, ICGM, France). Elemental analyses were performed by the Microanalysis Centre Service (CNRST, Rabat, Morocco). Melting points were measured using a Büchi 510 m.p. apparatus (LCAE, Oujda, Morocco). ¹H and ¹³C-NMR spectra were recorded using an AC

<table><caption>Table 4 Second-order perturbation analysis interactions in C₁₄ H₁₃ Br N₂ O₂</caption><thead><tr><th colspan="2">Donor (i)</th><th colspan="2">Acceptor (j)</th><th>E(2)</th><th>ε<sub>i</sub> - ε<sub>j</sub></th><th>F(i, j)</th></tr><tr><th colspan="2">NBO type occupation</th><th colspan="2">NBO type occupation</th><th>kcal.mol<sup>-1</sup></th><th>a.u.</th><th>a.u.</th></tr></thead><tbody><tr><td>O1</td><td>LP</td><td>C1–C2</td><td>π*</td><td>0.25875</td><td>44.40</td><td>0.36</td><td>0.114</td></tr><tr><td>N1</td><td>LP</td><td>C5–C6</td><td>π*</td><td>0.32998</td><td>38.56</td><td>0.31</td><td>0.101</td></tr><tr><td>C1–C2</td><td>π</td><td>O2–C3</td><td>π*</td><td>0.32146</td><td>32.54</td><td>0.28</td><td>0.087</td></tr><tr><td>N1</td><td>LP</td><td>N2–C4</td><td>π*</td><td>0.46366</td><td>31.65</td><td>0.28</td><td>0.084</td></tr><tr><td>C5–C6</td><td>π</td><td>N2–C4</td><td>π*</td><td>0.46366</td><td>29.85</td><td>0.27</td><td>0.085</td></tr><tr><td>C2–H2</td><td>σ</td><td>C13–C14</td><td>π*</td><td>0.35711</td><td>22.20</td><td>4.36</td><td>0.304</td></tr><tr><td>C11–C12</td><td>π</td><td>C9–C10</td><td>π*</td><td>0.37373</td><td>21.06</td><td>0.29</td><td>0.070</td></tr><tr><td>C9–C10</td><td>π</td><td>C11–C12</td><td>π*</td><td>0.32589</td><td>20.12</td><td>0.28</td><td>0.068</td></tr><tr><td>C13–C14</td><td>σ</td><td>C11–C12</td><td>π*</td><td>0.32589</td><td>20.03</td><td>0.30</td><td>0.070</td></tr><tr><td>C13–C14</td><td>σ</td><td>C9–C10</td><td>π*</td><td>0.37373</td><td>17.75</td><td>0.30</td><td>0.067</td></tr><tr><td>O2</td><td>LP</td><td>C3–C4</td><td>σ*</td><td>0.05429</td><td>16.30</td><td>0.75</td><td>0.100</td></tr><tr><td>N2–C4</td><td>π</td><td>O2–C3</td><td>π*</td><td>0.32146</td><td>15.55</td><td>0.31</td><td>0.064</td></tr><tr><td>C9–C10</td><td>π</td><td>C1–C2</td><td>π*</td><td>0.25875</td><td>14.87</td><td>0.29</td><td>0.060</td></tr><tr><td>N2–C4</td><td>π</td><td>C5–C6</td><td>π*</td><td>0.32998</td><td>13.19</td><td>0.34</td><td>0.062</td></tr><tr><td>O2</td><td>LP</td><td>O1–HO1</td><td>σ*</td><td>0.03204</td><td>12.68</td><td>1.12</td><td>0.108</td></tr><tr><td>O2</td><td>LP</td><td>C2–C3</td><td>σ*</td><td>0.04207</td><td>12.25</td><td>0.80</td><td>0.090</td></tr><tr><td>C1–C2</td><td>π</td><td>C9–C10</td><td>π*</td><td>0.37373</td><td>9.46</td><td>0.29</td><td>0.049</td></tr></tbody></table>

E(2) is the hyper conjugative interaction energy, εᵢ - εⱼ the energy difference between i and j NBO orbitals and F(i, j) the Fock matrix element between i and j NBO orbital

LP lone pair
<figure><img src="image_13.png" /><figcaption>**Fig. 6** Molecular structures of compounds 1–6</figcaption></figure>

<table><caption>**Table 5** Values of antifungal activity of the pyrazole-keto enol compounds against *Fusarium oxysporum* f.sp.</caption><thead><tr><th>Compounds</th><th>Volume withdrawn (µL)</th><th>Diameter of the strain in the presence of the drug (cm)</th><th>Inhibition (%)</th></tr></thead><tbody><tr><td rowspan="3">1</td><td>50</td><td>5</td><td>0</td></tr><tr><td>200</td><td>3.8</td><td>24</td></tr><tr><td>500</td><td>2.7</td><td>46</td></tr><tr><td rowspan="3">2</td><td>50</td><td>5</td><td>0</td></tr><tr><td>200</td><td>3.5</td><td>30</td></tr><tr><td>500</td><td>2.3</td><td>54</td></tr><tr><td rowspan="3">3</td><td>50</td><td>5</td><td>0</td></tr><tr><td>200</td><td>3.6</td><td>28</td></tr><tr><td>500</td><td>2.5</td><td>50</td></tr><tr><td rowspan="3">4</td><td>50</td><td>5</td><td>0</td></tr><tr><td>200</td><td>3.8</td><td>24</td></tr><tr><td>500</td><td>3.2</td><td>36</td></tr><tr><td rowspan="3">5</td><td>50</td><td>1.2</td><td>76</td></tr><tr><td>200</td><td>0.9</td><td>82</td></tr><tr><td>500</td><td>0.5</td><td>90</td></tr><tr><td rowspan="3">6</td><td>50</td><td>2.0</td><td>60</td></tr><tr><td>200</td><td>1.3</td><td>74</td></tr><tr><td>500</td><td>0.2</td><td>96</td></tr><tr><td rowspan="3">Benomyl</td><td>50</td><td>2.3</td><td>54</td></tr><tr><td>200</td><td>1.1</td><td>78</td></tr><tr><td>500</td><td>0.3</td><td>94</td></tr></tbody></table>

300 spectrometer (CNRST) (Bruker, LLN, Belgium) (300 MHz for ¹H and 75.47 MHz for ¹³C spectra). A JMS DX-300 mass spectrometer (JEOL, Rabat, Morocco) was used for the determination of molecular weights. Infrared (IR) spectra were recorded on a Shimadzu infrared spectrophotometer (LCAE, Oujda, Morocco) using the KBr disc technique. Geometry optimizations and DFT calculations were carried out using Dmol³ and Gaussian 093W programs (CNRS, ICGM, France).

### Procedure for the synthesis of pyrazole-ketoenol

To a solution of toluene (20 mL) containing metallic sodium (15.21 mmol) was added the pyrazolic carboxylate (12.01 mmol) solubilized in toluene (20 mL); then 3-bromophenyl methyl ketone (12.01 mmol) in toluene (10 mL) was added at 0 °C. The resulting mixture was stirred at room temperature for 2 days. The resulting precipitate was filtered, washed, dissolved in water, and neutralized with acetic acid to pH 5. The CH₂Cl₂ extracted fraction was dried over anhydrous sodium sulfate and concentrated to dry. The final product, as a white solid, was obtained after purification through silica gel column chromatography using CH₂Cl₂/MeOH in 28% yield. The β-ketoenol form was recrystallized from methanol (95%) to obtain the (Z)-3-(3-bromophenyl)-1-(1,5-dimethyl-1H-pyrazol-3-yl)-3-hydroxyprop-2-en-1-one compound which was confirmed by FT-IR, ¹H-NMR, ¹³C NMR, and mass spectroscopy. yield: 28%; m.p. 124 °C; Rf=0.52 (CH₂Cl₂/MeOH 9/1)/silica. IR (KBr, cm⁻¹): ν(OH)=3431; ν(C=O)=1676; ν (enolic C=C)=1531; ¹H NMR [CDCl₃, δ(ppm)]: 2.24 (s, 3H, Pz-CH₃); 3.78 (s, 3H, CH₃-N); 4.54 (s, 0.1H, keto, CH₂); 6.54 (s, 0.9H, enol, C-H); 7.39 (m, 3H, Ar-2H, Pz-H); 7.91(m, 2H, Ar-H). ¹³C NMR [CDCl₃, δ(ppm)]: 11.29 (1C, Pz-CH₃); 36.88 (1C, CH₃-N); 49.39 (1C, keto CH₂); 93.24 (1C, enol C-H); 106.14 (1C, =CH, Pz); 127.01 (1C, Ar-C3); 128.53 (1C, Ar-C6); 132.04 (1C, Ar-C5); 133.31(1C, Ar-C2); 134.97 (1C, Ar-C1); 140.37 (1C, Ar-C4); 147.72 (2C, PzN=C, PzN-C=); 181.99 (1C, C-OH); 183.19 (1C, C=O). Anal. Calcd. for C₁₄H₁₃BrN₂O₂: C 52.36, H 4.08, N 8.72. Found: C 52.25, H 4.12, N 8.65. *m/z*: 320.97.

### X-ray data collection and treatment

Fairly regularly shaped crystals were selected using a stereomicroscope equipped with a polarizing filter. Diffraction intensities were collected at room temperature within the complete diffraction sphere on the four-circle diffractometer (Mo-Kα radiation, λ = 0.71073 Å) and data
reduction was carried out using CrysAlis software [57]. The lattice dimensions and corresponding standard deviations were determined by least-squares method from the entire data set of reflections. Full-matrix least-squares refinements on F² used the complete data set of 18,366 collected reflections (including symmetry equivalent and redundant) of which 6297 are unique and 3027 observed according to the criterion I > 2σ(I). The diffracted intensities were corrected for Lorentz and polarization effects. The structure solution and subsequent refinements were performed using SHELX-2013 program packages [58]. Atoms positions and anisotropic displacement parameters were refined for all non-hydrogen atoms. The hydrogen atoms at the keto enol OH group were detected in the final Fourier difference and were treated as riding, following the HFIX/AFIX instructions, in the final refinement (even if they could have been freely refined). The hydrogen atoms have been considered with an isotropic displacement parameter equal to 1.2 times (1.5 for terminal –CH₃) the Ueq of the parent atom. Molecular pictures are drawn with ORTEP-3 for windows [59].

Full CIF file can be obtained free of charge via http://www.ccdc.cam.ac.uk/conts/retrieving.html (or from the CCDC, 12 Union Road, Cambridge CB2 1EZ, UK; Fax: +44 1223 336033; E-mail: deposit@ccdc.cam.ac.uk).

## Computation details

Calculations were performed within the framework of the density functional theory DFT using DMol³ module [60, 61] provided in Materials Studio software. A DNP basis set (Double Numerical with extra Polarization function on all atoms) was used in these calculations, its size is equivalent to the Gaussian 6-31G*. All-electrons geometry optimizations were performed at fine quality level, for both the molecule and the periodic crystal packing, with LDA-PWC, GGA-PW91 and GGA-BLYP functionals. The Fukui and Parr functions have been computed to give a description of the global reactivity. Gaussian 03 W optimizations using the Berny analytical gradient method with B3LYP functional and 6.31G+(d,p) basis set have been done prior to calculate the reactivity indices and to perform the NBO analysis [53].

## Anti-fungal tests

In vitro antibacterial and antifungal activities were tested by the agar diffusion technique (ADT) [54, 55]. ADT has been investigated using susceptibility test of NCCLS (National Committee for Clinical Laboratory Standards) recommended by the WHO and the French standard NF-U-47-107 AFNOR 2004. The agar media were inoculated with test organisms and a solution of the tested compound in DMSO/EtOH (50/50) was added to different concentrations in the culture media. The growth was

followed by counting the bacteria or the yeast colonies and by measuring the mycelium diameter. The inhibition percentage of a molecule is equal to the ratio of the colonies number (or the mycelium diameter of the culture) in the presence of a dose of the tested compound over the colonies number (or the mycelium diameter) of the reference culture, multiplied by 100. The minimum inhibition concentration (MIC) is the least dose of the compound which causes an inhibition of the micro-organism growth. Calculation of the concentration IC₅₀ was done using the same bacterial inocula mentioned above with decreasing concentration of the tested products. The D° value was measured for each culture at 625 nm. The inhibition percentage (%) is expressed as (D° – Dx)/D° × 100 where D° is the diameter of the mycelial growth of the culture witness and Dx the diameter of the mycelial growth in the presence of the product to be tested.

## Conclusion

A novel pyrazole-β-ketoenol compound has been synthesized, it has been characterized by NMR and IR techniques and its XRD single crystal structure was determined. Density functional calculations are used to evaluate the HOMO–LUMO energy gap, the molecular electrostatic potential (MEP) and to provide a natural bond orbital (NBO) analysis. From reactivity indices, the present molecule displays a moderate electrophile character. Computed Parr functions quite well agree with Fukui functions, indicating the position of the nucleophile and electrophile centers in the molecule. The title compound has been tested against *Fusarium oxysporum* f.sp. albedinis FAO fungal strains and three bacterial strains (*Escherichia coli*, *Bacillus subtilis*, and *Micrococcus luteus*). The measured activities encourage us to continue searching for other structures, likely to be good antifungal candidates.

## Additional file

Additional file 1: Figure S1. ¹H NMR spectrum of 1. Figure S2. ¹³C NMR spectrum of 1. Figure S3. ¹³C NMR-DEPTQ-135 spectrum of 1. Figure S4. Mass spectrum of 1. Figure S5. FT-IR spectrum of 1. Figure S6. UV-Vis spectrum of 1.

## Abbreviations

DFT: density functional theory; NBO: natural bond orbital; MEP: molecular electrostatic potential; HOMO: highest occupied molecular orbital; LUMO: lowest unoccupied molecular orbital; HIV: human immunodeficiency virus; LDA: local density approximation; DNP: double numerical plus polarization.

## Authors' contributions

ST, SR, KH and KK carried out of the experimental work and cooperated in the preparation of the manuscript. MT performed the density functional theory calculations. MT and DE determined the X-ray crystal structure. RB carried out the Biological activity. All authors read and approved the final manuscript.
Author details

¹ Laboratoire de Chimie Appliquée et Environnement (LCAE), Faculté des Sciences, Université Mohamed I, 60000 Oujda, Morocco. ² Département de Biologie, Faculté des Sciences, Université Mohamed I, 60000 Oujda, Morocco.

³ ICGM, CNRS, Université de Montpellier, ENSCM, Montpellier, France.

⁴ Laboratoire de Chimie du Solide Minéral et Analytique, Faculté des Sciences, Université Mohamed I, 60000 Oujda, Morocco. ⁵ ENSA-AL Hociema, Université Mohamed I, 60000 Oujda, Morocco. ⁶ Laboratoire de Chimie Thérapeutique, Faculté de médecine et de Pharmacie, Université Mohamed V, Rabat, Morocco.

Acknowledgements

The authors extend their appreciation to the PPR2-MESRSFC-CNRS-T10 project (Morocco) for its supporting this work.

Competing interests

The authors declare that they have no competing interests.

Availability of data and materials

All samples of the synthesized compounds as well as their data are available from the authors.

Consent for publication

Not applicable.

Ethics approval and consent to participate

Not applicable.

Funding

Not applicable.

Publisher's Note

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

Received: 20 April 2018 Accepted: 15 November 2018

Published online: 26 November 2018

References

1. Dowling C, Murphy VJ, Parkin G (1996) Bis(pyrazolelyethyl) ether ligation to zinc and cobalt: meridional vs facial coordination and the suitability of such ligands in providing a NNO donor set for modeling bioinorganic aspects of zinc chemistry. *Inorg Chem* 35:2415–2420

2. Haanstra WJ, Driessen WL, Van Roon M, Stoffels ALE, Reedijk J (1992) Co-ordination compounds with the N 2 S-donor ligand 1, 5-bis (3, 5-dimethyl-ylpyrazol-1-yl)-3-thiapentane. *J Chem Soc Dalton Trans* 3:481–486

3. Driessen WL, Wiesmeijer WGR, Schipper-Zablotskaja M, De Graaff RAG, Reedijk J (1989) Transition metal coordination compounds of two pyrazole-substituted ammonia ligands. X-ray structure of [Biss(1-pyrazolylmethyl)aminecobalt(II)] bis(nitrate). *Inorg Chim Acta* 162:233–238

4. Pennings YCM, Driessen WL, Reedijk J (1988) Copper(I) coordination compounds with two bidentate chelating pyrazole ligands. X-ray crystal structures of DI-µ-chloro-bis[[N, N-bis(1-pyrazolylmethyl)aminoethane] copper(I)] and [bis[N, N-bis(1-pyrazolylmethyl)aminoethane]copper(I) triflate. *Polyhedron* 7:2583–2589

5. Veldhuis JB, Driessen WL, Reedijk J (1986) A pyrazole derivative of aminoethane as a tridentate chelating ligand towards transition metals. The X-ray structure of [N, N-bis(pyrazol-1-ylmethyl)aminoethane] dibromocopper(II). *Chem Soc Dalton Trans* 3:537–541

6. Blonk HJ, Driessen WL, Reedijk J (1985) Transition-metal co-ordination compounds of a novel aniline-based pyrazole derivative. X-Ray crystal structures of [NN-bis(3,5-dimethylpyrazol-1-ylmethyl) aminobenzene]-dichlorocobalt(II) and -dibromocopper(II). *Chem Soc Dalton Trans* 8:1699–1705

7. Mukherjee R (2000) Coordination chemistry with pyrazole-based chelating ligands: molecular structural aspects. Coord Chem Rev 203:151–218

8. Trofimenko S (1993) Recent advances in poly (pyrazolyl) borate (scorpionate) chemistry. *Chem Rev* 93:943–980

9. Trofimenko S (1986) The coordination chemistry of pyrazole-derived ligands. *Prog Inorg Chem* 34:115–210

10. Radi S (2016) Acyclic and macrocyclic polypyrazoles: Coordination chemistry. LAP Lambert Academic Publishing, Saarbrücken

11. Bouabdallah I, M'barek LA, Zyad A, Ramadan A, Zidane I, Melhaoui A (2006) Anticancer effect of three pyrazole derivatives. *Nat Prod Res* 20:1024–1030

12. Janus SL, Magdif AZ, Erik BP, Claus NN (1999) Synthesis of triazeno-pyrazole derivativesas potential inhibitors of HIV-1. *Monatsh Chem* 130:1167–1174

13. Tewari AK, Mishra A (2001) Synthesis and anti-inflammatory activities of N4, N5-disubstituted-3-methyl-1H-pyrazolo [3, 4-c] pyridazines. *Bioorg Med Chem* 9:715–718

14. Wustrow DJ, Capiris T, Rubin R, Knobelsdorf JA, Akunne H, Davis MD, MacKenzie R, Pugsley TA, Zoski KT, Heffner TG (1998) Pyrazolo [1, 5-a] pyrimidine CRF-1 receptor antagonists. *Bioorg Med Chem Lett* 8:2067–2070

15. Pimerova EV, Voronina EV (2001) Antimicrobial activity of pyrazoles and pyridazines obtained by intraction of 4-aryl-3-arylhydrazono-2, 4-dioxobutanoic acids and their esters with hydrazines. *Pharm Chem J* 35:602–604

16. Goldgur Y, Craigir R, Cohen GH, Fujiwara T, Yoshinaga T, Fujishita T, Sugimoto H, Endo T, Murai H, Davies DR (1999) Structure of the HIV-1 integrase catalytic domain complexed with an inhibitor: a platform for antiviral drug design. *Proc Natl Acad Sci USA* 96:13040–13043

17. Hazuda DJ, Felock P, Witmer M, Wolfe A, Stillmock K, Grobler JA, Espeseth A, Gabryelski L, Schleif W, Blau C (2000) Inhibitors of strand transfer that prevent integration and inhibit HIV-1 replication in cells. *Science* 287:646–650

18. Minassi A, Sanchez-Duffhues G, Collado JA, Munoz E (2013) Appendino G. Dissecting the pharmacophore of curcumin. Which structural element is critical for which action? *J Nat Prod* 76:1105–1112

19. Pluymers W, Pais G, Maele BV, Pannecouque C, Fikkert V, Burke JTR, de Clercq E, Witvrouw M, Neamati N, Debyser Z (2002) Inhibition of human immunodeficiency virus type 1 integration by diketo derivatives. *Antimicrob Agents Chemother* 46:3292–3297

20. Adams BK, Cai J, Armstrong J, Herold M, Lu YJ, Sun A, Snyder JP, Liotta DC, Jones DP (2005) EF24, a novel synthetic curcumin analog, induces apoptosis in cancer cells via a redox-dependent mechanism. *Anti-Cancer Drug.* 16:263–275

21. Tan KL, Ali A, Du Y, Fu H, Jin HX, Chin TM, Khan M, Go ML (2014) Synthesis and evaluation of bisbenzyldienedioxotetrahydrothiophranones as activators of endoplasmic reticulum (ER) stress signaling pathways and apoptotic cell death in acute promyelocytic leukemic cells. *J Med Chem* 57:5904–5918

22. Liang G, Shao L, Wang Y, Zhao C, Chu Y, Xiao J, Zhao Y, Li X, Yang S (2009) Exploration and synthesis of curcumin analogues with improved structural stability both in vitro and in vivo as cytotoxic agents. *Bioorg Med Chem* 17:2623–2631

23. Ishikawa Y, Fujii S (2011) Binding mode prediction and inhibitor design of anti-influenza virus diketo acids targeting metalloenzyme RNA polymerase by molecular docking. *Bioinformation.* 6:221–225

24. Anand P, Thomas SG, Kunnumakkara AB, Sundaram C, Harikumar KB, Sung B, Tharakan ST, Misra K, Priyadarsini IK, Rajasekharan KN (2008) Biological activities of curcumin and its analogues (Congeners) made by man and Mother Nature. *Biochem Pharmacol* 76:1590–1611

25. Liang G, Li X, Chen L, Yang S, Wu X, Studer E, Gurley E, Hylemon PB, Ye F, Li Y (2008) Synthesis and anti-inflammatory activities of mono-carbonyl analogues of curcumin. *Bioorg Med Chem Lett* 18:1525–1529

26. Hui YY, Shu HM, Hu HM, Song J, Yao HL, Yang XL, Wu QR, Yang ML, Xue GL (2010) Syntheses, structures and magnetic properties of tetranuclear and trinuclear nickel (II) complexes with β-diketone-functionalized pyridine-carboxylate ligand. *Inorg Chim Acta* 363:3238–3243

27. Verna PN, Juneja HD (2012) Synthesis and characterization of 1-(2', 4'-dihydroxy-5-nitrophenyl)-3-(pyridin-3-yl)- propane-1, 3-dione and its metal complexes. *J Chem Tech Res.* 4:1000–1006

28. Radi S, Tighadouini S, Feron O, Riant O, Bouakka M, Benabbes R, Mabkhot YN (2015) Synthesis of novel β-Keto-enol derivatives tethered pyrazole, pyridine and furan as new potential antifungal and anti-breast cancer agents. *Molecules* 20:20189–20194

29. Tighadouini S, Radi S, Bacquet M, Degoutin S, Zaghrïoui M, Jodeh S, Warad I (2017) Removal efficiency of Pb(II), Zn (II), Cd (II) and Cu (II) from
aqueous solution and natural water by ketoenol-pyrazole receptor functionalized silica hybrid adsorbent. Sep Sci Technol 52:608–621

30. Radi S, Tighadouini S, Bacquet M, Degoutin S, Garcia Y (2016) New hybrid material based on a silica-immobilised conjugated β-ketoenol-bipyridine receptor and its excellent Cu (II) adsorption capacity. Anal Methods 8:6923–6931

31. Radi S, Tighadouini S, Bacquet M, Degoutin S, Dacquin JP, Eddike D, Tillard M, Mabkhot YN (2016) β-keto-enol tethered pyridine and thiophene: synthesis, crystal structure determination and its organic immobilization on silica for efficient solid-liquid extraction of heavy metals. Molecules 7:888–900

32. Radi S, Tighadouini S, El-Massaoudi M, Bacquet M, Degoutin S, Revel B, Mabkhot YN (2015) Thermodynamics and kinetics of heavy metals adsorption on silica particles chemically modified by conjugated β-ketoenol furan. J Chem Eng Data 60:2915–2925

33. Radi S, Tighadouini S, Bacquet M, Degoutin S, Revel B, Zaghrioui M (2015) Quantitative removal of Zn (II) from aqueous solution and natural water using new silica-immobilized ketoenol-pyridine receptor. J Env Chem Eng. 3:1769–1778

34. Radi S, Tighadouini S, Eddike D, Tillard M, Mabkhot YN (2017) Crystal structure of (Z)-1-(1, 5-dimethyl-1H-pyrazol-3-yl)-3-hydroxy-3-(4-methoxyphenyl) prop-2-en-1-one, C15H16N2O3. Z Kristallogr NCS. 232:199–200

35. Radi S, Tighadouini S, Eddike D, Tillard M, Mabkhot YN (2017) Crystal structure of (Z)-1-(1, 5-dimethyl-1H-pyrazol-3-yl)-3-hydroxy-3-phenyl-prop-2-en-1-one, C14H14N2O2. Z Kristallogr NCS. 232:201–202

36. Radi S, Tighadouini S, Eddike D, Tillard M, Mabkhot YN (2017) Crystal structure of (Z)-1-(1, 5-dimethyl-1H-pyrazol-3-yl)-3-(4-ethoxyphenyl)-3-hydroxyprop-2-en-1-one, C16H18N2O3. Z Kristallogr NCS. 232:207–208

37. Radi S, Tighadouini S, Eddike D, Tillard M, Mabkhot YN (2017) Crystal structure of (Z)-3-hydroxy-3-(4-methoxyphenyl)-1-(pyridin-2-yl)prop-2-en-1-one, C15H13NO3. Z Kristallogr NCS. 232:209–210

38. Radi S, Tighadouini S, Eddike D, Tillard M, Mabkhot YN (2017) Crystal structure of (Z)-3-hydroxy-3-(4-methoxyphenyl)-1-(pyridin-2-yl)prop-2-en-1-one, C15H13NO3. Z Kristallogr NCS. 232:235–236

39. Becke AD (1988) A multicenter numerical integration scheme for polyatomic molecules. J Chem Phys. 88:2547–2553

40. Lee C, Yang W, Parr RG (1988) Development of the Colle-Salvetti correlation-energy formula into a functional of the electron density. Phys Rev B. 37:785–789

41. Perdew JP, Wang Y (1992) Accurate and simple analytic representation of the electron-gas correlation energy. Phys Rev B. 45:13244–13249

42. Yang W, Parr RG (1985) Hardness, softness, and the fukui function in the electronic theory of metals and catalysis. Proc Natl Acad Sci USA 82:6723–6726

43. Weber J, Roch M, Williams AF (1986) Molecular electrostatic potentials and electron deformation densities of chromium(V), molybdenum(VI) and niobium(V) tetraperoxo complexes. Chem Phys Let. 123:246–253

44. Berger B (2013) Using conceptual density functional theory to rationalize regioselectivity: a case study on the nucleophilic ring-opening of activated aziridines. Comput Theor Chem. 1010:11–18

45. Khaled KF (2010) Experimental, density function theory calculations and molecular dynamics simulations to investigate the adsorption of some thiourea derivatives on iron surface in nitric acid solutions. Appl Surf Sci 256:6753–6763

46. Koubsky T, Kalvoda L (2015) Application of ab initio molecular electronic structure calculations of radiolytic and hydrolytic stabilities of prospective extractants. J Radioanal Nucl Chem 304:227–235

47. Lakbaibi Z, Abou El Makarim H, Tabyaoui M, El Hajbi A (2014) Theoretical study of the formation of. α.-chloroglycidic esters in aliphatic series using the quantum DFT method with B3LYP/6-311G (d, p). Intern J Innov Appl Stud. 7:602–616

48. Cedillo A, Contreras R (2012) Stability and bonding in the borane-H2 complexes. J Mex Chem Soc 56:257–260

49. Domingo LR, Pérez P, Saez JA (2013) Understanding the local reactivity in polar organic reactions through electrophilic and nucleophilic Parr functions. RSC Adv. 3:1486–1494

50. Domingo LR, Ríos-Gutiérrez M, Pérez P (2016) Applications of the conceptual density functional theory indices to organic chemistry reactivity. Molecules 21:748–769

51. Thomson HW, Torkington P (1945) The vibrational spectra of esters and ketones. J Chem Soc 171:640–645

52. Snehalatha M, Ravikumar C, Joe IH, Sekar N, Jayakumar VS (2009) Spectroscopic analysis and DFT calculations of a food additive Carmoisine. Spectrochim Acta A: 72:654–662

53. Frisch MJ, Trucks GW, Schlegel HB, Scuseria GE, Robb MA, Cheeseman JR et al (2003) Gaussian 03, Revision E.01. Gaussian Inc., Pittsburgh

54. Carrod LP, Grady ED (1972) Antibiotics and chemotherapy, 3rd edn. Churchill Livingstone, Edinburgh, p 477

55. Neri F, Mari M, Brigati S (2006) Control of Penicillium expansum by plant volatile compounds. Plant Pathol 55:100–105

56. Sudha G, Priya MS, Shree RI, Vadivukkarasi S (2011) In vitro free radical scavenging activity of raw pepino fruit (Solanum muricatum aiton). Int J Curr Pharm Res. 3:137–140

57. Rigaku Oxford Diffraction (2015) CrysAlisPro software system, version 1.171.38.41. Rigaku Corporation, Oxford

58. Sheldrick GM (2013) SHELXS: A program for crystal structures solution; SHELXL: A program for refining crystal structures. University of Göttingen, Göttingen

59. Farrugia LJ (2012) ORTEP-3 for windows, an update. J Appl Crystallogr 45:849

60. Delley B (1990) An all-electron numerical method for solving the local density functional for polyatomic molecules. J Chem Phys. 92:508

61. Delley B (2000) From molecules to solids with the DMol³ approach. J Chem Phys. 113:7756

Ready to submit your research? Choose BMC and benefit from:

• fast, convenient online submission

• thorough peer review by experienced researchers in your field

• rapid publication on acceptance

• support for research data, including large and complex data types

• gold Open Access which fosters wider collaboration and increased citations

• maximum visibility for your research: over 100M website views per year

At BMC, research is always in progress.

Learn more biomedcentral.com/submissions

<figure><img src="image_14.png" /></figure>