Abstract

Drug discovery for the management of neurological disorders is a challenging arena in medicinal chemistry. Vibrational spectral studies of (2E)-2-[3-(1H-imidazol-1-yl)-1-phenylpropylidene]-N-(4-methylphenyl)hydrazinecarboxamide ((2E)-IPPMP) have been recorded and analyzed to identify the functional groups and intermolecular/intramolecular interactions of the title molecule. The blue shift of the C-H stretching wavenumber reveals the presence of improper C-H⋯O hydrogen bonding. The equilibrium geometry, harmonic vibrational wavenumbers, Frontier orbital energy, and natural bond orbital analyses have been carried out using density functional theory with a B3LYP/6-311++G(d,p) level of the basis set. The vibrational modes have been unambiguously assigned using potential energy distribution analysis. The scaled wavenumbers are in good agreement with the experimental results. Natural bond orbital analysis has confirmed the intermolecular/intramolecular charge transfer interactions. HOMO-LUMO analysis was carried out to explore charge delocalization on the (2E)-IPPMP molecule. A molecular docking study has supported the anticonvulsant activity of the title molecule.

1. Introduction

Epilepsy is a rather neurobiological group of disorders. It has multiple origins and aspects depending upon the affected brain areas. It affects nearly 50 million people of the world’s population [13]. The hydrazinecarboxamide derivatives have a wide spectrum of biological activities. Among these activities are anticancer and antioxidant [4], antifertility [5], antimicrobial [6, 7], anticonvulsant [8, 9], and anti-inflammatory [10] activities. The title molecule (2E)-2-[3-(1H-imidazol-1-yl)-1-phenyl-propylidene]-N-(4-methylphenyl)hydrazinecarboxamide ((2E)-IPPMP) was synthesized in our laboratory and its crystal structure was previously reported [11]. (2E)-IPPMP exhibited anticonvulsant activity with 83% and 50% seizure protection at a dose level of 718 μmol/kg in subcutaneous pentylenetetrazole (scPTZ) and maximal electroshock seizure (MES) screens, respectively, without any neurotoxicity [9].

Literature screening indicated that computational studies on the (2E)-IPPMP molecule have not yet been reported. Therefore, detailed investigations of structural properties and vibrational spectral analysis of the (2E)-IPPMP molecule were carried out in the present study using density functional theory (DFT) computations. Moreover, the biological activity of the title molecule has been predicted by molecular docking analysis. It is expected that the current investigations might support the development of new potent anticonvulsant agents.

2. Experimental

2.1. General

Melting point was recorded on a Gallenkamp melting point instrument and it is uncorrected. The Fourier transform infrared spectrum of MMIMI was recorded on a Perkin Elmer RXL spectrometer (Waltham, Massachusetts, USA) in the region 4000–400 cm−1, with samples in the KBr pellet method. The resolution of the spectrum was 2 cm−1. The FT-Raman spectrum was measured in the range 3500–50 cm−1 using a Bruker RFS 100/S FT-Raman spectrophotometer (Ettlingen, Germany) with a 1064 nm Nd:YAG laser source of 100 mW power (Göttingen, Germany).

2.2. Synthesis

A solution containing N-(4-methylphenyl)hydrazinecarboxamide [12] (1.65 g, 10 mmol), 3-(1H-imidazol-1-yl)-1-phenylpropan-1-one (2.00 g, 10 mmol) [13], and a few drops of glacial acetic acid in ethanol (15 mL) was stirred at room temperature for 18 h. The reaction mixture was evaporated under reduced pressure and the residue was crystallized from ethanol to give 1.67 g (48%) of the title compound as colorless crystals (m.p. 476–478 K) which were suitable for single crystal X-ray analysis. 1H and 13C NMR as well as the mass spectral data of the title compound 2 are in accordance with the previously reported ones [9].

2.3. Theoretical Calculations

All DFT calculations of the (2E)-IPPMP molecule were performed using the Gaussian 09 program package [14] at the Becke3-Lee-Yang-Parr (B3LYP) level with a 6-311++G(d,p) basis set [1517]. The structural parameters were computed in the gas phase as well as in the liquid phase using a polarizable continuum model (PCM) method. In order to correct the overestimations arising from some negative factors such as basis set truncation effect, neglecting electron correlations, and anharmonicity characters of the vibrational modes, the calculated wavenumbers were scaled using a uniform scaling factor of 0.9673 [18, 19]. Theoretical NMR calculation was performed on the basis of the GIAO (gauge-independent atomic orbitals) theory method using a Gaussian program.

3. Results and Discussion

3.1. Synthesis

The target compound (2E)-IPPMP was obtained in a three-step reaction sequence as given in Scheme 1.

3.2. Structural Geometry Analysis

The structure of the (2E)-IPPMP molecule was optimized using the B3LYP method with a 6-311++G(d,p) basis set. The optimized molecular structure of the isolated molecule is shown in Figure 1. The optimized geometrical parameters of the isolated (2E)-IPPMP molecule in the gas and water phases are given in Table 1. The calculated values were compared with the X-ray diffraction results. In the (2E)-IPPMP molecule the two phenyl rings are bridged by a hydrazinecarboxamide skeleton containing an imidazole ring. The C-N bond lengths C17-N16 (1.4179 Å) and C17-N19 (1.3792 Å) are shorter than the normal single C-N bond length 1.480 Å [20]. This discrepancy is due to the conjugation of p-type electrons of the carbonyl group and nitrogen atom, allowing the electrons to smear out along the C-N bond. In the para-disubstituted phenyl ring, the calculated C25-H44 (1.0789 Å) bond length is shorter than that of the other C-C bonds; also, O18⋯H44 distance is 2.25 Å, which is significantly shorter than the van der Waals radius (2.72 Å) [21] between O and H atoms, which indicates the possibility of C-H⋯O hydrogen bonding. The elongation of the C20-C21 (1.4004 Å) bond is due to the transfer of a lone pair of electrons from the amide nitrogen to the carbon atom. In the other phenyl ring, the C5-H31 (1.0825 Å) bond is shorter than that of the other C-C bonds; also, N15⋯H31 distance is 2.56 Å, which is significantly shorter than the van der Waals radius (2.75 Å) between N and H atoms, which indicates the possibility of C-H⋯O hydrogen bonding. In addition, C5-C6 (1.4024 Å) and C4-C5 (1.4012 Å) bond distances on either side of the cyanide group are appreciably greater than the other C-C bonds due to the resonance effect between the cyanide group and phenyl ring. In the imidazole ring, C14-N10 (1.3676 Å) and C14-N13 (1.3139 Å) bonds are relatively shorter due to the lone pair interaction from the nitrogen atom. The deviation of bond lengths of N15-H39 (1.0123 Å) and N19-H40 (1.0092 Å) is due to the different environment of the nitrogen atom. Conjugation of the carbonyl group with the hydrazine moiety would favor planarity, but the van der Waals repulsion between H39 and H40 hinders the achievement of coplanarity. A small deviation was obtained within the calculated structural parameters in the gas and liquid phases, which were due to the solvent interactions over the solution phase calculations.

3.3. Natural Bond Orbital Analysis

Natural bond orbital (NBO) analysis was performed using the NBO 3.1 program [22] as implemented in the Gaussian 09 program package for the DFT method. The corresponding results are presented in Table 2. NBO analysis has proved to be an effective tool for chemical interpretation of hyperconjugative interaction and electron density transfer (EDT) from a filled lone pair to an unfilled antibonding orbital in the hydrogen bonding system [2325]. The intramolecular C-H⋯O hydrogen bonding is formed due to the orbital overlap between and which results in intramolecular charge transfer (ICT) causing stabilization of H-bonded systems. This interaction results in increased electron density (ED) of the C-H antibonding orbital, which strengthens the C-H bond. NBO analysis confirms C-H⋯O intramolecular hydrogen bonding formed by the orbital overlap between a lone pair and antibonding orbital with stabilization energy of 0.78 kcal/mol. The most important interaction and () energies of and are 46.48 and 26.84 kcal/mol, respectively. This larger value shows ICT interactions of the molecule.

3.4. Vibrational Spectral Analysis

Computed vibrational wavenumbers and atomic displacements corresponding to the different normal modes are used to identify vibrational modes. The vibrational modes are assigned on the basis of potential energy distribution analysis using the VEDA4 program [26]. The experimental infrared (IR) and Raman spectra are shown in Figures 2 and 3. The calculated vibrational wavenumbers, measured IR, and Raman band positions and their detailed assignments are given in Table 3.

3.4.1. Phenyl Ring Vibrations

Aromatic C-H stretching vibrations generally absorb in the region 3080–3010 cm−1 [27]. The observed weak IR band at 3119 cm−1 and Raman band at 3121 cm−1 correspond to aromatic C-H stretching mode. The blue shift of the C-H stretching wavenumber is due to weak intramolecular C-H⋯O hydrogen bonding. Aromatic C=C stretching vibrations occur in the region 1625–1430 cm−1 [28]. A medium IR band at 1539 cm−1 and a strong Raman band at 1591 cm−1 are observed, which correspond to aromatic C=C stretching mode. In-plane C-H deformation vibrations appear in the region 1290–1000 cm−1 [27]. Observed IR and Raman bands at 1288 and 1275 cm−1 are assigned to in-plane C-H deformation. In-plane ring deformation vibrations appear in the region 650–615 cm−1. Observed IR and Raman bands at 615 and 640 cm−1 are assigned to in-plane ring deformation.

3.4.2. Methylene Vibrations

Asymmetric and symmetric CH2 stretching vibrations normally appear strongly at about 2926 and 2855 cm−1 [29]. The Raman band at 2941 cm−1 is assigned to CH2 symmetric stretching vibration. Methylene scissoring vibrations normally appear in the region 1465–1445 cm−1 [29]. A medium band observed in IR at 1401 cm−1 is attributed to methylene scissoring mode. The twisting, wagging vibrations appear in the region 1422–719 cm−1 [30]. The observed strong IR band at 1141 cm−1 and weak Raman band at 1177 cm−1 are assigned to CH2 twisting modes for methylene. Wagging is observed in IR at 1352 cm−1 and Raman at 1361 cm−1.

3.4.3. Methyl Vibrations

The asymmetric C-H stretching mode of CH3 generally occurs at 2982–2962 cm−1 and CH3 symmetric stretching at 2882–2862 cm−1 [29]. A weak Raman band at 2977 cm−1 is attributed to CH3 symmetric stretching. Asymmetric bending vibrations of methyl groups occur in the region 1470–1450 cm−1 [28]. The observed weak Raman band at 1447 cm−1 is assigned to methyl scissoring vibration.

3.4.4. Secondary Amide Vibrations

Carbonyl stretching vibration in a secondary amide is expected in the region 1680–1630 cm−1 [28]. A very strong band at 1654 cm−1 is assigned to C=O stretching. N-H stretching vibrations generally appear in the region 3370–3170 cm−1. Observed IR and Raman bands at 3141 and 3140 cm−1 are attributed to N-H stretching. The in-plane N-H bending vibration usually appears from 1570 to 1515 cm−1 [28]. The observed very strong IR band at 1523 cm−1 and weak Raman band at 1527 cm−1 are assigned to in-plane N-H bending vibration.

3.4.5. Imidazole Vibrations

Imidazole C-H stretching vibrations are expected in the region 3145–3115 cm−1 [31, 32]. Observed Raman bands at 3121 and 3138 cm−1 and IR bands at 3119 and 3138 cm−1 are assigned to C-H stretching mode. The observed Raman band at 1338 cm−1 is attributed to C-C stretching mode [33].

3.4.6. Hydrazine Vibrations

In accordance with earlier reports and in agreement with the calculation, a weak band observed at 3207 cm−1 is assigned to hydrazine N-H stretching [27, 34]. Observed bands in IR and Raman at 1089 and 1088 cm−1 are assigned to N-N stretching vibration.

3.4.7. Skeletal Mode Vibrations

C-N and C-C stretching vibrations generally occur in the region 1150–850 cm−1 [34]. The weak IR and Raman bands observed at 1011 and 998 cm−1 are attributed to C-C stretching.

3.5. Frontier Molecular Orbital Energy Analysis

HOMO (highest occupied molecular orbital) and LUMO (lowest unoccupied molecular orbital) are important in defining the reactivity of a chemical species. The energy of HOMO indicates nucleophilicity and LUMO indicates electrophilicity [35]. HOMO is spread over the methyl phenyl, hydrazinecarboxamide fragment and LUMO is located on ph 2. This shows the charge transfer between the two rings through the hydrazinecarboxamide path. The HOMO (−6.61 eV) and LUMO (−1.61 eV) energies reflect the charge transfer within the molecule. The HOMO-LUMO energy gap is 4.46 eV. The frontier molecular orbital diagrams are shown in Figure 4.

3.6. Natural Population Analysis

Natural population analysis provides an effective method to calculate atomic charges and electron distribution within a molecule [36]. The net atomic charges of the (2E)-IPPMP molecule obtained by natural population analysis are plotted in Figure 5. All hydrogen atoms have a net positive charge. The atoms H40 (0.3852e) and H39 (0.3788e) show more positive charge than other hydrogen atoms due to their attachment with a nitrogen atom. Among the hydrogen atoms (~0.2023e) of the phenyl ring, H44 (~0.2452e) shows the highest positive charge, being involved in C-H⋯O intramolecular hydrogen bonding. All carbon atoms are negatively charged except C7, C14, C17, and C20 due to their attachment with electronegative nitrogen or oxygen atoms. The atom C17 (0.7991e) shows more positive charge and N19 (−0.6181e) shows more negative charge, indicating charge delocalization in the molecule.

3.7. NMR Analysis

The scaled and experimental NMR (1H and 13C) chemical shift values for the (2E)-IPPMP molecule are presented in Table 4. The phenyl and imidazole ring carbon signals usually appear in the region 115–150 ppm. In this molecule, imidazole ring carbon signals are obtained at 119.4, 126.4, and 137.3 ppm, which are predicted at 119.36, 128.75, and 140.01 ppm, respectively.

Phenyl carbon signals were observed at 120.0, 128.4, 128.8, 128.9, 131.5, 136.3, 136.8, and 137.3 ppm, while their respective calculated ones were obtained at 115.12, 128.13, 128.53, 129.41, 137.03, 132.85, 138.17, and 140.01 ppm. A carbonyl carbon signal was seen at 144.9 ppm and its computed one was obtained at 150.09 ppm. This deviation may occur due to the presence of amide⋯amide interactions in the crystalline state. On the other hand, the computed 1H chemical shift values for the title molecule showed good agreement with the experimental ones (Table 4).

3.8. Molecular Docking Analysis

The title molecule (2E)-IPPMP was energy minimized based on the DFT method. Molecular docking was performed using AutoDock 4.2. A target protein (PDB ID: 1EOU) for antiepileptic agents was selected for the present docking analysis [37, 38]. The protein data bank file of the target protein was downloaded from the Research Collaboratory for Structural Bioinformatics (RCSB) database, with a resolution of 2.1 Å. Protein preparation was carried out by the following steps: (i) all water molecules were removed, (ii) hydrogen atoms were added to the crystal structure, (iii) Kollman charge was added, and (iv) a previously docked inhibitor (fructose-based sugar sulfamate RWJ-37497) was removed from the protein. Rigid protein and flexible ligand dockings were carried out using the AutoDock 4.2 program package [39] with the Lamarckian genetic algorithm, applying the following protocol: trials of 100 dockings, energy evaluations of 25000000, population size of 200, mutation rate of 0.02, crossover rate of 0.8, and elitism value of 1. The docking results were evaluated by sorting the docked conformations according to their predicted binding free energy. The protein-ligand interaction complex is given in Figure 6, displaying the conformer with the best predicted binding free energy (−7.94 kcal/mol). The amino acids ASN11, TYR7, and LYS169 in the active sites of the target protein bind with the (2E)-IPPMP ligand by N-H⋯O and N-H⋯N hydrogen bonds. These preliminary results support the exhibited anticonvulsant activity of the title molecule.

4. Conclusions

The geometry optimization and harmonic wavenumbers of the (2E)-IPPMP molecule have been performed at ground state calculations using the DFT method. FT-IR and FT-Raman measurements helped functional group identification of the title molecule. The fundamental wavenumbers are in good agreement with the theoretical results. Shifting of vibrational wavenumbers and hyperconjugative results confirmed the presence of intermolecular/intramolecular interactions in the molecule. The molecular docking results predicted the anticonvulsant activity of the (2E)-IPPMP molecule due to its ability to interact with a target protein (1EOU) for anticonvulsants. The results of the current study will support the development of new drug-like candidates in the anticonvulsant research area.

Competing Interests

The authors have declared that there are no competing interests.

Acknowledgments

The authors would like to extend their sincere appreciation to the Deanship of Scientific Research at King Saud University for its funding of this research through the Research Group Project no. RGP-196.