By the functional B3LYP and M05-2X of DFT and in two bases set, more and more extended (6-311G and 6-311G (d, p)), theoretical study of antioxidant properties of four hydrazones was carried out. The calculations made concern the geometrical, spectroscopic and electronic parameters of the molecules. Analysis of the results relating to the geometrical parameters was carried out by calculating interatomic distances, relative errors between calculated values and those obtained experimentally by X-ray diffraction found in the literature. The 13C NMR spectra were calculated by GIAO (Gauge Including Atomic Orbitals) methods, and the results were subjected to statistical analysis by calculating Mean Absolute Deviation (MAD), Root Mean. Square (RMS) and the correlation coefficient (R 2), in comparison with experimental spectra. The analysis of the results of calculations of various electronic parameters (hardness ( η), softness (S), electronegativity ( χ), electrophile index ( ω), energy gap (HOMO-LUMO)) reveals that, overall, the methods M05-2X/6-311G (d, p) and B3LYP/6-311G (d, p) found that (R) - (?) - carvone salicylhydrazone (N 2) is the most antioxidant molecule of the four molecules and also classify them according to their stability. This confirms the results obtained on the antitrypanosomal activity, the toxicity, the cytotoxicity and the selectivity of the synthesized compounds.
? |
Nowadays, hydrazones are an important class of biologically active compounds [
Numerous biological activities of hydrazones are well known. Among them are antibacterial [
The main objective of this work is to study theoretically the biological activities of hydrazones. It will specifically study the antioxidant properties of these molecules and compare the theoretical results with those published in the literature.
Four hydrazones were studied in this study:(S)-(+)-carvone salicylhydrazone,(R)-(−)-carvone salicylhydrazone,(S)-(+)-carvone benzohydrazone,(R)-(−)-carvone benzohydrazone.
The four molecules were synthesized in the Laboratory of Physical and Synthetic Organic Chemistry (LACOPS) of the Department of Chemistry of the Faculty of Science and Technology (FAST) of the University of Abomey-Calavi (UAC).
The studied molecules constituting systems of more than four atoms, the implementation of the DFT with functional ones such as B3LYP and MO5-2X will make it possible to calculate the spectroscopic, electronic and thermodynamic quantities of our molecules with a good precision.
In the case of this work, the studied molecules being of large size, we chose to use the base of Pople [
performed using the Gaussian 09 program [
For molecule with several degrees of freedom, the potential energy surface is a hyper-surface with several minima and maxima. The maxima correspond to transition states and the minima to equilibrium geometries or conformations of the molecule. In a state of minimal energy, the vibration frequencies of the molecule are all positive. The lowest of minima called global minima, represents the most stable conformation of the studied system. For each of the molecules studied, by the B3LYP and M05-2X functions of the DFT, geometric parameters relating to a minimum energy state were calculated, in the base 6-311G (d, p).
Also the relative error associated with the calculation of each parameter was evaluated. The lower this relative error, the closer the calculated theoretical value is to the experimental data.
The experimental data considered to evaluate the relative errors related to the calculation of the geometrical parameters of the molecules are those of hydrazide (see
The 1H NMR and 13C NMR spectra of molecules were calculated with respect to DMSO (dimethylsulfoxide) taken as a reference molecule by the [
base 6-311G (d, p). Calculations were made using the GIAO method.
1) Mean Absolute Deviation (MAD)
MAD = 1 N ∑ i = 1 N | δ i e x p − δ i c a l |
2) Root Mean Square (RMS)
RMS = ∑ i = 1 N ( δ i e x p − δ i c a l ) 2 N
3) correlation coefficient R2
R 2 = ∑ i = 1 N ( δ i e x p − δ i e x p ¯ ) ⋅ ( δ i c a l − δ i c a l ¯ ) ∑ i = 1 N ( δ i e x p − δ i e x p ¯ ) 2 ⋅ ∑ i = 1 N ( δ i c a l − δ i c a l ¯ ) 2
These three statistical quantities must build us on suitability of the combination Method/Base judicious for the study carried out. More lowly the values of MAD and RMS are, closer the theoretical results are to the experimental data, and the more convenient are the functional and the base. And, a value of the correlation coefficient (R2) close to 1, indicates a good quality of the theoretical results.
For each of the molecules denoted ArOH, the calculated electronic parameters were:
・ The energy of the highest occupied molecular orbital (Ehomo).
・ The energy of the lowest molecular orbital vacant (Elumo).
・ The G a p ( HOMO-LUMO ) = E lumo − E homo , [
・ The hardness (η) which expresses the resistance of a molecule to the change of its number of electrons or to the charge transfer. The harder the hardness, the less reactive the molecule is [
η = ( I P − E A ) 2
・ Softness (S), defined as the inverse of hardness [
S = 1 2 η
・ Electronegativity (χ), which measures the tendency of a chemical species to attract electrons [
χ = ( I P + E A ) 2
・ The electrophilic index (ω), which represents the stabilizing energy of a molecule saturated by electrons coming from its surroundings [
ω = χ 2 2 η
The dipole moment (μ) of a molecule expresses the value of its great polarity. It constitutes a very important index of reactivity which makes it possible to define the biological properties related to the interaction with the active sites of the molecule. In fact, the polar molecules are defined by a non-zero dipole moment whereas the apolar ones are defined by a zero dipole moment [
The results from this section constitute a theoretical tool for geometric, spectroscopic and electronic characterization of four hydrazones by quantum chemistry methods. It consisted in calculating the geomeric, spectroscopic, and electronic parameters relating to transformations to which they can lend themselves.
Tables 1-4 show the experimental and calculated values of some geometric parameters (lengths) of the (R) - (−) - carvone benzohydrazone molecule, the (S) - (+) - carvone benzohydrazone, (R) - (−) - carvone salicyldrazone and (S) - (+) - carvone salicyldrazone, on the one hand, and related relative errors, on the other hand. Analysis of these tables will make a comparative study of the experimental and calculated values.
Analysis of the results in
The analysis of the results in
At the end of the results of
δ (ppm) | EXP (ppm) [ | B3LYP | M05-2X | ||
---|---|---|---|---|---|
X | Δx | x | Δx | ||
C10C11 C11C13 | 1.39 1.39 | 1.40 1.39 | 0.01 0.00 | 1.39 1.39 | 0.01 0.00 |
C13C16 | 1.39 | 1.39 | 0.00 | 1.39 | 0.00 |
C15C16 | 1.47 | 1.40 | 0.07 | 1.39 | 0.07 |
C15C12 | 1.39 | 1.39 | 0.00 | 1.39 | 0.00 |
C12C10 | 1.38 | 1.39 | 0.01 | 1.39 | 0.01 |
C10C8 | 1.50 | 1.51 | 0.01 | 1.50 | 0.01 |
C8O9 C8N6 | 1.22 1.36 | 1.21 1.39 | 0.01 0.03 | 1.21 1.38 | 0.01 0.03 |
N6N5 | 1.39 | 1.36 | 0.04 | 1.35 | 0.04 |
δ (ppm) | EXP (ppm) [ | B3LYP | M05-2X | ||
---|---|---|---|---|---|
X | Δx | x | Δx | ||
C10C11 C11C13 | 1.39 1.39 | 1.40 1.39 | 0.01 0.00 | 1.39 1.39 | 0.00 0.00 |
C13C16 | 1.39 | 1.39 | 0.00 | 1.39 | 0.00 |
C15C16 | 1.47 | 1.40 | 0.07 | 1.39 | 0.08 |
C15C12 | 1.39 | 1.39 | 0.00 | 1.39 | 0.00 |
C12C10 | 1.38 | 1.39 | 0.01 | 1.39 | 0.01 |
C10C8 | 1.50 | 1.51 | 0.01 | 1.50 | 0.01 |
C8O9 C8N6 | 1.22 1.36 | 1.21 1.39 | 0.01 0.03 | 1.21 1.38 | 0.01 0.02 |
N6N5 | 1.39 | 1.35 | 0.04 | 1.35 | 0.04 |
δ (ppm) | EXP (ppm) [ | B3LYP | M05-2X | ||
---|---|---|---|---|---|
X | Δx | x | Δx | ||
C10C11 C11C13 | 1.39 1.39 | 1.40 1.39 | 0.01 0.00 | 1.39 1.39 | 0.01 0.00 |
C13C16 | 1.39 | 1.39 | 0.00 | 1.39 | 0.00 |
C15C16 | 1.47 | 1.40 | 0.07 | 1.39 | 0.08 |
C15C12 | 1.39 | 1.39 | 0.00 | 1.39 | 0.00 |
C12C10 | 1.38 | 1.39 | 0.01 | 1.39 | 0.01 |
C12O30 | 1.50 | 1.51 | 0.01 | 1.50 | 0.00 |
C10C8 | 1.39 | 1.36 | 0.03 | 1.36 | 0.03 |
C8O9 C8N6 | 1.22 1.36 | 1.21 1.39 | 0.01 0.03 | 1.21 1.38 | 0.01 0.03 |
N6N5 | 1.39 | 1.35 | 0.04 | 1.35 | 0.04 |
functional with respect to (R) - (−) - carvone salicylhydrazone.
The analysis of results of
At this stage of the calculations, we do not find any significant difference between structural results given by the two functional ones. Since the theoretical calculations are carried out in the gaseous phase and taking into account the effects of solvents on the atoms, the comparison of certain theoretical results with those of the experimental ones will be made on the basis of the statistical parameters obtained.
δ (ppm) | EXP (ppm) [ | B3LYP | M05-2X | ||
---|---|---|---|---|---|
X | Δx | x | Δx | ||
C10C11 C11C13 | 1.39 1.39 | 1.40 1.39 | 0.10 0.00 | 1.39 1.39 | 0.00 0.00 |
C13C16 | 1.39 | 1.39 | 0.00 | 1.39 | 0.00 |
C15C16 | 1.47 | 1.40 | 0.07 | 1.39 | 0.07 |
C15C12 | 1.39 | 1.39 | 0.00 | 1.39 | 0.00 |
C12C10 | 1.38 | 1.39 | 0.01 | 1.39 | 0.01 |
C12O30 C10C8 | 1.50 1.39 | 1.51 1.36 | 0.01 0.03 | 1.50 1.36 | 0.00 0.03 |
C8O9 C8N6 | 1.22 1.36 | 1.21 1.39 | 0.01 0.03 | 1.21 1.38 | 0.01 0.03 |
N6N5 | 1.39 | 1.36 | 0.03 | 1.35 | 0.03 |
In
Since the theoretical calculations are carried out in the gaseous phase and do not take into account the effects of solvents on the atoms, the comparison of theoretical results of
The values of absolute mean deviation (MAD), standard deviation (RMS) and correlation coefficient (R2), relative to the 13C experimental and theoretical NMR spectra of each of the molecules, were calculated by functional B3LYP and M05-2X using the GIAO method. The results are shown in
These results show that all correlation coefficient values are close to unity (1). This means that approximations taken into account in the calculation methods used are good. However, it can be generally observed that the absolute average deviation values calculated by the functional B3LYP are lower than those calculated by the functional M05-2X. For the calculated standard deviation values, the same trend is observed. The standard deviation values calculated by the functional B3LYP are lower than those calculated by the functional M05-2X. Finally, the functional M05-2X and the base 6-311G (d, p) seem really suitable for calculating the 13C NMR spectrum of the (R) - (−) - carvone benzohydrazone, of (S) - (+) - carvone benzohydrazone, (R) - (−) - carvone salicylhydrazone and (S) - (+) - carvone salicylhydrazone.
(R) - (−) - carvone benzohydrazone | (S) - (+) - carvone benzohydrazone | ||||||
---|---|---|---|---|---|---|---|
δ (ppm) | EXP (ppm) [ | 6-311G (d, p) | δ (ppm) | EXP (ppm) [ | 6-311G (d, p) | ||
B3LYP | M05-2X | B3LYP | M05-2X | ||||
C20 | 157.68 | 161.66 | 168.62 | C18 | 20.44 | 20.43 | 12.41 |
C2 | 33.16 | 38.01 | 28.07 | C7 | 20.45 | 24.63 | 15.68 |
C3 | 29.47 | 33.98 | 23.68 | C3 | 31.73 | 32.62 | 22.51 |
C36 | 153.99 | 161.66 | 166.87 | C2 | 29.47 | 39.04 | 29.77 |
C27 | 151.80 | 154.64 | 160.50 | C20 | 110.32 | 42.52 | 41.43 |
C24 | 147.67 | 149.86 | 158.39 | C1 | 40.10 | 45.79 | 36.65 |
C30 | 133.87 | 136.62 | 145.83 | C22 | 116.77 | 47.41 | 36.37 |
C8 | 20.45 | 20.09 | 12.59 | C26 | 33.16 | 47.88 | 38.11 |
C33 | 132.28 | 68.42 | 58.49 | C5 | 40.10. | 49.73 | 146.44 |
C10 | 40.10 | 42.93 | 32.82 | C23 | 133.87 | 51.54 | 42.09 |
C11 | 130.85 | 55.39 | 168.62 | C8 | 132.28 | 70.40 | 59.85 |
C12 | 116.77 | 47.92 | 36.53 | C25 | 119.59 | 138.30 | 39.71 |
C13 | 119.59 | 51.56 | 47.05 | C34 | 148.67 | 150.97 | 159.68 |
C15 | 117.95 | 49.61 | 41.55 | C37 | 153.99 | 151.91 | 158.44 |
C16 | 110.32 | 47.38 | 35.74 | C11 | 157.68 | 158.25 | 166.29 |
C32 | 20.44 | 20.55 | 9.14 | C30 | 161.80 | 162.41 | 167.73 |
(S) - (−)- carvone salicylhydrazone | (R) - (+)- carvone salicylhydrazone | ||||||
---|---|---|---|---|---|---|---|
δ (ppm) | EXP (ppm) [ | 6-311G (d, p) | δ (ppm) | EXP (ppm) [ | 6-311G (d, p) | ||
B3LYP | M05-2X | B3LYP | M05-2X | ||||
C18 | 20.45 | 18.78 | 08.67 | C20 | 162.85 | 161.92 | 169.85 |
C22 | 148.00 | 19.87 | 09.23 | C27 | 152.78 | 154.54 | 159.99 |
C7 | 29.10 | 26.68 | 17.79 | C2 | 20.44 | 38.12 | 169.85 |
C3 | 148.52 | 132.40 | 166.87 | C3 | 20.45 | 34.54 | 24.41 |
C2 | 147.17 | 139.01 | 30.29 | C35 | 31.73 | 68.80 | 59.20 |
C1 | 107.20 | 47.00 | 158.39 | C32 | 29.47 | 136.89 | 145.56 |
C26 | 163.15 | 148.04 | 35.43 | C24 | 148.52 | 150.17 | 158.06 |
C21 | 36.16 | 49.83 | 38.58 | C7 | 40.10 | 166.00 | 159.99 |
C20 | 31.75 | 54.46 | 48.46 | C8 | 16.77 | 20.95 | 13.43 |
C23 | 40.10 | 60.16 | 53.01 | C9 | 33.16 | 161.60 | 39.73 |
C25 | 33.16 | 64.95 | 55.63 | C10 | 40.10 | 48.40 | 44.92 |
C8 | 40.70 | 70.11 | 58.41 | C11 | 33.87 | 50.10 | 40.28 |
- | - | - | - | C12 | 32.28 | 20.66 | 12.91 |
- | - | - | - | C13 | 119.59 | 60.46 | 51.87 |
- | - | - | - | C15 | 48.67 | 65.87 | 58.30 |
- | - | - | - | C16 | 53.99 | 48.22 | 37.19 |
- | - | - | - | C34 | 17.68 | 20.63 | 09.62 |
BASES | MAD | RMS | R2 | |||
---|---|---|---|---|---|---|
B3LYP | M05-2X | B3LYP | M05-2X | B3LYP | M05-2X | |
N1 | 6.52 | 8.13 | 9.15 | 11.46 | 0.936 | 0.945 |
N2 | 4.53 | 5.21 | 7.17 | 8.65 | 0.949 | 0.936 |
M1 | 4.31 | 4.43 | 6.12 | 7.96 | 0.908 | 0.951 |
M2 | 4.22 | 3.13 | 6.02 | 7.12 | 0.918 | 0.925 |
It is found that the four molecules do not have the same values of the electronic parameters.
Note that:
• hardnes of the four molecules are weak and almost equal, which would mean that it is easy for these molecules to exchange electrons with the environment.
• Gap (HOMO-LUMO) of the enantiomer (R) - (−) - carvone salicylhydrazone is greater than the Gap (HOMO-LUMO) relative to enantiomer (S) - (+) - carvone salicylhydrazone. With B3LYP/6-311G (d), p) method; therefore (R) - (−) - carvone salicylhydrazone is more stable than (S) - (+) - carvone salicylhydrazone. As for the enantiomers of the benzohydrazone, they have the same values of the Gap, so same stability with this method. Similarly, we find that (R) - (−) - carvone benzohydrazone is more stable than (S) - (+) - carvone salicylhydrazone and (R) - (−) - carvone salicylhydrazone is more stable than (S) - (+) - carvone benzohydrazone.
• However, (S) - (+) - carvone salicylhydrazone is less stable than the other three compounds.
• On the other hand, with the method M05-2X/6-311G (d, p), the enantiomers of the carvone salicylhydrazone have the same stability. It is the same for the enantiomers of carvone benzohydrazone. In addition, the enantiomers of carvone salicylhydrazone are more stable than those of carvone benzohydrazone.
On the basis of these series of results, the following observations can be made:
➢ With the method B3LYP/6-311G (d, p), Gap values (HOMO-LUMO) were obtained in ascending order: (S) - (+) - carvone salicylhydrazone, (R) - (−) -
Molecules | B3LYP/6-311G (d, p) | |||||||
---|---|---|---|---|---|---|---|---|
EHOMO | ELUMO | Gap | χ | η | S | ω | μ | |
M2 | −0.21 | −0.06 | 0.15 | 0.13 | 0.08 | 6.25 | 0.11 | 3.77 |
N2 | −0.22 | −0.05 | 0.17 | 0.13 | 0.09 | 5.56 | 0.09 | 4.81 |
N1 | −0.22 | −0.06 | 0.16 | 0.14 | 0.08 | 6.25 | 0.12 | 4.68 |
M1 | −0.22 | −0.06 | 0.16 | 0.14 | 0.08 | 6.25 | 0.12 | 4.94 |
M05-2X/6-311G (d, p) | ||||||||
Molecules | EHOMO | ELUMO | Gap | χ | η | S | ω | μ |
M2 | −0.28 | −0.02 | 0.26 | 0.15 | 0.13 | 3.85 | 0.09 | 4.12 |
N2 | −0.27 | −0.01 | 0.26 | 0.14 | 0.13 | 3.85 | 0.08 | 5.17 |
N1 | −0.27 | −0.02 | 0.25 | 0.15 | 0.13 | 3.85 | 0.09 | 4.96 |
M1 | −0.27 | −0.02 | 0.25 | 0.15 | 0.13 | 3.85 | 0.09 | 5.23 |
carvone benzohydrazone, (S) - (+) - carvone benzohydrazone, (R) - (−) - carvone salicylhydrazone. These results indicate that (S) - (+) - carvone salicylhydrazone and (R) - (−) - carvone benzohydrazone ((S) - (+) - carvone benzohydrazone) are more antioxidant than (R) - (−) -. -carvone salicylhydrazone. For the M05-2X/6-311G (d, p) method, however, (R) - (-) - carvone benzohydrazone and (S) - (+) - carvone benzohydrazone are more antioxidant than the (S) - (+) - carvone salicylhydrazone and (R) - (−) - carvone salicylhydrazone.
➢ The results obtained by the DFT methods B3LYP/6-311G (d, p) for the calculation of softness (S), found that the highest values were obtained for (S) - (+) - carvone salicylhydrazone, (R) - (−) - carvone benzohydrazone, and the (S) - (+) - carvone benzohydrazone thus displaying them as the best antioxidants of the four molecules. There really was not a difference in the results given by the method M05-2X/6-311G (d,p).
➢ With the B3LYP/6-311G (d, p) method, the lowest values of the electronegativity index (χ) were obtained for (S) - (+) - carvone salicylhydrazone and (R) - (−) - carvone salicylhydrazone. This makes (S) - (+) - carvone salicylhydrazone and (R) - (−) - carvone salicylhydrazone appear to be the most antioxidant of the four molecules. On the other hand, at the M05-2X/6-311G (d, p) approximation levels, the lowest index value was obtained for (R) - (−) - carvone salicylhydrazone. This compound is the most antioxidant by the method M05-2X/6-311G (d, p).
➢ Electrophilic index calculations (ω) carried out by the B3LYP/6-311G (d, p) method revealed that (S) - (+) - carvone salicylhydrazone and (R) - (−) - carvone salicylhydrazone are the most electrophilic and therefore the strongest antioxidants. In addition, (R) - (−) - carvone salicylhydrazone is more antioxidant than (S) - (+) - carvone salicylhydrazone.
➢ By calculations M05-2X/6-311G (d, p), it is rather the (R) - (−) - carvone salicylhydrazone which appears as the most antioxidant molecules.
➢ Values of the dipolar moment obtained made it possible to note that the (S) - (+) - carvone salicylhydrazone is the most polar of the four molecules whatever the method of calculation used.
The analysis of the results of calculations of the various electronic parameters, electronic affinities, ionization energy, hardness, softness, electronegativity, electrophile index, energy gap (HOMO-LUMO) reveals that, overall, the methods M05-2X/6-311G (d, p) and B3LYP/6-311G (d, p) showed that (S) - (+) - carvone salicylhydrazone and the (R) - (−) - carvone salicylhydrazone are the most antioxidant molecules of the four. In addition, (R) - (−) - carvone salicylhydrazone (N2) is more antioxidant than (S) - (+) - carvone salicylhydrazone (M2). The difference in the antioxidant powers observed in (R) - (−) - carvone salicylhydrazone (N2) and (S) - (+) - carvone salicylhydrazone (M2) is due to the descriptors of (R) and (S) configurations.
Theoretical study of the antioxidant properties of the four hydrazones (R) - (−) - carvone benzohydrazone, (S) - (+) - carvonebenzohydrazone, (R) - (−) - carvone salicylhydrazone and ((S) - (+) - carvone salicylhydrazone) was carried out, by the functional groups B3LYP and M05-2X of the DFT, in the bases 6-311G and 6-311G (d, p). For this, different parameters have been calculated. It’s about:
- Geometric parameters of the molecules. The comparison of the calculated values with experimental X-ray diffraction data from the literature did not tangibly retain one of the functionalities as the most suitable for the study.
- Spectroscopic parameters relating to the 13C NMR spectra of the three molecules. The statistical analysis of chemical shift values calculated by the GIAO method revealed that the M05-2X functional appears to be the most appropriate for calculating the 13C NMR spectrum of each molecule.
- Electronic parameters such as hardness (η), softness (S), electronegativity (χ), electrophilic index (ω), energy gap (HOMO-LUMO) of the molecules. The theoretical results obtained have shown that (R) - (−) - carvone salicylhydrazone (N2) is the most antioxidant molecule of the four and also classified according to their stability. The experimental results are confirmed by the theoretical results. The method M05-2X/6-311G (d, p) appears best indicated in the calculation of electronic parameters.
In perspective, we plan to:
- To enhance the biological activity of the compounds by formation of the metal complexes.
- To carry out quantum-chemical studies on the complexes formed.
Houngue, M.T.A.K., N’bouke, M., Atchade, B., Doco, R.C., Kuevi, U.A., Kpotin, G.A., Kpoviessi, S.D.S., Atohoun, Y.G.S., Badawi, M. and Mensah, J.-B. (2018) Quantum Chemical Studies of Some Hydrazone Derivatives. Computa- tional Chemistry, 6, 1-14. https://doi.org/10.4236/cc.2018.61001