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

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 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 \(\left( {1\bar{1}\bar{1}} \right)\) and separated by ~ 3.5 Å, each 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) 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 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 outs 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.

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 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 _k ⁺ and P _k ⁻ are powerful tools to study the molecular reactivity and they are well adapted to localize 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_HOMO + E_LUMO) = –4.37 eV) and the chemical hardness η (η = E_LUMO− E_HOMO = 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_HOMO − E_HOMO(TCE) where E_HOMO(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

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.

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-metoxy 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.

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 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 G03W 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. Diffracted 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^o 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.