Journal of the Mexican Chemical Society, vol. 63, no. 4, 2019
Sociedad Química de México A.C.
Mirjana S. Jankulovska * msjankulovska@gmail.com
Cyril and Methodius University, Faculty of Agricultural Sciences and Food in Skopje, Macedonia. e-mail: msjankulovska@gmail.com (Milena Jankulovska)., Macedonia
Ilinka Spirevska spirevska@yahoo.com
Cyril and Methodius University, Faculty of Natural Sciences and Mathematics in Skopje, Macedonia. e-mail: spirevska@yahoo.com (Ilinka Spirevska)., Macedonia
Vesna Dimova vdimova@tmf.ukim.edu
Cyril and Methodius University, Faculty of Technology and Metallurgy in Skopje, Macedonia. e-mail: vdimova@tmf.ukim.edu (Vesna Dimova)., Macedonia
Milena Jankulovska
University Ss. Kliment Ohridski, Chemical-Bromatological-Toxicological Analysis Laboratory, Higher Medical School in Bitola, Macedonia., Macedonia
Date received: 03 March 2019
Date accepted: 15 July 2019
Abstract.: The spectral behavior of some p-nitro-p-substituted benzoylhydrazones in the perchloric acid media was followed, applying the UV spectroscopy. The position of the absorption maximum in the spectra was defined in acidic media and the electronic transitions were discussed, as well (1<pH<7). The equilibrium between neutral and protonated form was investigated in the ethanol-water (V/V, 1:1) solutions. The observed changes in the UV spectra suggested that protonation process took place in one step. The pH region of protonation ranges between 1.4 and 2.9. Using the changes in the UV spectra which appear as a result of the protonation reaction the stoichiometric dissociation constants were determined numerically (pKBH+ = n·pH + logI) and graphically (intercept of the dependence of logI on pH). Thermodynamic dissociation constants were estimated as an intercept of dependence of pKBH+ on square root of the ionic strength. In order to achieve that, measurements were performed at different ionic strengths: 0.1, 0.25 and 0.5 mol/dm3, adjusted with sodium perchlorate. The obtained thermodynamic pKBH+ values ranged between 2.07 and 2.58. In order to predict proton transfer at a given pH, semiempirical methods AM1 and PM3 were applied. The influence of the substituents present in the p-position of the benzene ring on pKBH+ values of investigated hydrazones was discussed, too. Total energy, binding energy, enthalpy of formation, Gibbs energies of formation, atomic charge and proton affinity values were used to predict protonation site in hydrazone molecule. Furthermore, the stability and the proton affinity of the isomers (E and Z) in which hydrazones exist and their protonated forms were defined.
Keywords: p-nitro-p-substituted benzoylhydrazones, dissociation constants, protonation, UV spectroscopy, semiempirical AM1 and PM3 methods.
Resumen.: Se siguió el comportamiento espectral de algunas benzoilhidrazonas p-nitro-p sustituidas utilizando ácido perclórico como medio y aplicando espectroscopía UV. La absorción máxima en los espectros se definió en medios ácidos y también se discutieron las transiciones electrónicas (1<pH<7). Se investigó el equilibrio entre la forma neutra y protonada en soluciones de etanol-agua (V / V, 1: 1). Los cambios observados en los espectros UV sugirieron que el proceso de protonación tuvo lugar en un sólo paso. La región de protonación oscila en el intervalo de pH entre 1.4 y 2.9. Utilizando los cambios en el espectro de UV que surgen como resultado de la reacción de protonación, las constantes de disociación estequiométricas fueron determinadas numérica (pKBH+ = n · pH + logI) y gráficamente (intercepción de la dependencia de logI en el pH). Las constantes de disociación termodinámicas fueron estimadas como una intersección en la dependencia de pKBH+ con la raíz cuadrada de la fuerza iónica. Para conseguir lo anterior, se realizaron mediciones a diferentes fuerzas iónicas: 0.1, 0.25 y 0.5 mol / dm3, ajustadas con perclorato de sodio. Los valores termodinámicos de pKBH+ obtenidos oscilaron entre 2.07 y 2.58. Para predecir la transferencia de protones a un pH dado, se utilizaron los métodos semiempíricos AM1 y PM3. También, se discutió la influencia de los sustituyentes presentes en la posición p del anillo de benceno de los valores pKBH+ de las hidrazonas investigadas. Para predecir el sitio de protonación en la molécula de hidrazona, se utilizaron los valores de energía total, energía de unión, entalpía de formación, energía de Gibbs, carga atómica y afinidad protónica. Además, se definió la estabilidad y la afinidad protónica de isómeros (E y Z) en las cuales existen las hidrazonas y sus formas protonadas.
Palabras clave: Benzoilhidrazonas p-nitro-p sustituidas, constantes de disociación, protonación, espectroscopía UV, métodos semi empíricos AM1 y PM3.
Introduction
Aroylhydrazones and their derivatives constitute a versatile class of compounds in organic chemistry. They are a class of azomethines having the structure: Ar-CO-NH-N=CR1R2. Aroylhydrazones contain two interlinked nitrogen atoms of different nature and nucleophilic character, and a carbon atom which may act as either electrophile or nucleophile, according to the reaction environment [1,2]. In their structure is also presented a C=N double bond that is conjugated with a lone electron pair of the terminal nitrogen atom. These structural fragments are mainly responsible for the physical and the chemical properties of hydrazones. Both nitrogen atoms of the hydrazone group are nucleophilic, although the amino type nitrogen is more reactive. The carbon atom of hydrazone group has both electrophilic and nucleophilic character. Due to their structure hydrazones are important compounds for drug design and, as possible ligands for metal complexes, organocatalysis and for the syntheses of heterocyclic compounds [3]. There is a growing interest in the structural studies of Aroylhydrazones, as they show a broad spectrum of applications in pharmaceutical and industrial fields [4]. These observations increased the need for the development of new hydrazones that possess varied biological activities [5-7]. Literature studies revealed that hydrazones, and various substituted hydrazones are associated with a broad spectrum of biological activities such as: antioxidant, antibacterial, antiviral, analgesic, antiplatelet, antimicrobial, and anticancer activities etc [8,9].
In order to combat diseases with minimal toxicity and maximal effects, hydrazones are synthesized as drugs by many researchers. These predictions have provided therapeutic pathway for developing a new effective biologically active hydrazones [10]. Hydrazones are very effective organic compounds, but at the same time they are used as intermediates to prepare some biological active compounds by using the active hydrogen atom of CONHN=CH- azomethine group [6,8].
Hydrazones also have analytical applications, which attract so many researchers [11]. They are very important group of analytical regents for determination of various metal ions, by using different analytical techniques [12]. Hydrazone derivatives may act as multidentate ligands and their transition metals complexes were already used in the treatment of tuberculosis, in colorimetric or fluorimetric determinations, as well as, in applications involving catalytic processes [13,14]. Hydrazones are also important intermediates in the syntheses of nitrogen containing heterocyclic compounds [3]. Furthermore, hydrazones possess pesticide effects (herbicides, insecticides, nematocides, rodenticides) and plant growth regulators activity [15]. Many hydrazone derivatives have been reported to possess broad spectrum of insecticidal activities, and they can be used as active ingredients for controlling agricultural and horticultural pests [16,17].
The biological activity of hydrazones and their derivatives depends on the ionic form which exists in the solution. Hydrazones have acidic and basic functional groups and their ionization state is controlled by both, pH of the solution, and the value of dissociation constants. Dissociation constant is the most valuable parameter to understand chemical phenomenon such as: biological activity, absorption and extent of ionization of compound in different pH. So, it is the key parameter in drug development and optimization of the dose [18]. Dissociation constants are also of vital importance in understanding and optimization of analytical procedures like acid-base titration, solvent extraction and complex formation. The accurate determination of dissociation constant values is often required in many chemical and biochemical areas.
There are various techniques available for determination of dissociation constant such as: potentiometric titration, nuclear magnetic resonance (NMR), capillary electrophoresis (CE), chromatography (HPLC), hyper-Rayleigh scattering technique, UV spectroscopic methods etc. The ultraviolet-visible spectrophotometry is one of the most frequently employed techniques for determination of dissociation constant [18]. Knowing the pKBH+ values of a drug, and the pH of the absorption site, the relative amount of ionized and unionized drug in a solution at a particular pH, and the percent of drug ionized at this pH, can be determined [19,20]. Because of that, there is a considerable interest in using experimental methods for determination of the pKBH+ values of different types of organic compounds.
The aim of the present work was to investigate the acid-base behavior of some p-nitro-p-substituted benzoylhydrazones, in perchloric acid media (7<pH<1), as well as, to determine the stoichiometric dissociation constant of protonated form at different ionic strength, and thermodynamic dissociation constant values at constant ionic strength. In order to eliminate the influence of the solvent on the protonation reaction, method of characteristic vector analysis (CVA) was employed [21,22]. Protonation site of hydrazone molecule, stability of protonated form of E and Z isomers and influence of the substituent were determined using semiempirical data obtained applying semiempirical methods AM1 and PM3 [23,24].
Experimental
Chemicals and instrumentation
The investigated benzoylhydrazones were purified by recrystallization from ethanol or diluted ethanol (ethanol/water = 1/1, V/V). The purity was tested by measuring the melting point, as well as, by elemental analysis. The perchloric acid, sodium perchlorate and ethanol were of analytical grade p.a. (Merck), and they were used without further purification. A digital pH meter with glass electrode was used for measurements of the pH values (pH range from 1 to 14). The spectral measurements were carried out on a Varian Cary 50 spectrophotometer controlled by a computer and equipped with a 1 cm path length quartz cell, in the wavelength region from 190 to 400 nm. The maximum scan rate was 24 000 nm/min with the resolution of 1.5 nm. An Excel program was applied for calculation of the dissociation constants, while the UV spectra were obtained with computer program Grams Version 4.10.
Stock and test solutions
The stock solutions of investigated hydrazones were prepared by dissolving an accurate weight of the compound in the appropriate volume of 98 % ethanol, in a volumetric flask of 250 cm3. The concentration of the investigated hydrazones in these solutions was about 1·10-3 mol/dm3. The stock solutions were stable over a long period of time, under ordinary conditions, which was confirmed with the recorded UV spectra. The test solution was prepared from the stock solution with appropriate dilution in order to achieve concentration of investigated hydrazones in the test solution of about 3·10-5 mol/dm3. Because of the low solubility of investigated hydrazones in aqueous solutions the test solutions were prepared in redistilled water and ethanol with volume ratio 1:1 (V/V). The stability of these solutions was satisfactory for only 24 hours of time. A constant amount of sodium perchlorate (1 mol/dm3) was added in the test solutions, in order to maintain ionic strength (0.1, 0.25 and 0.5 mol/dm3). The pH values of the test solution were adjusted in the range between 1 and 7 with a suitable amount of perchloric acid (0.5 mol/dm3). The UV spectra were taken immediately after preparation of the test solutions, at room temperature. After that, the pH of each test solution was measured. At the same time, the blanks were prepared with the same composition as the test solutions. Three series of the test solutions were prepared, and all the measurements were performed in triplicate.
Structure of investigated hydrazones
Hydrazones possess an azomethine -NHN=CH group, thus they are considered as derivatives of aldehydes and ketones in which the oxygen atom has been replaced by the −NNH2 functional group [25]. Actually, hydrazones are formed by the reaction of hydrazine or hydrazide with aldehydes and ketones in ethanol [9,26]. In this investigation we used five p-nitro-p-substituted benzoylhydrazones previously synthesized in our laboratory [27]. The structural and molecular formula of investigated hydrazones as well as melting points are given in Table 1.
The investigated hydrazones were structurally characterized by UV spectroscopy, infrared spectroscopy (IR), nuclear magnetic resonance (1H NMR and 13C NMR), as well as, by elemental analysis [28].
pKBH+ calculations
The dissociation constants of investigated hydrazones were calculated using the Henderson - Hasselbach equation (Eq. 1) [20].
where, pKBH+ is the dissociation constant of the protonated form, I is the ionization ratio, and n is number of transferred protons.
The ionization ratio I (I = c(BH+)/c(B)) is actually the ratio between the concentration of the protonated (BH+) and neutral (B) form of the hydrazones. The calculations of ionization ratio were performed using the molar absorption coefficient values and the absorbance data obtained at four selected wavelengths. In accordance with Beer’s law (Eq. 2) an over determined system of four equations (absorbance data) with two unknown parameters (concentration of neutral and protonated form) was obtained.
Absorbance value which corresponds to the absorption maximum was used in calculations. Actually, the obtained results were more accurate when four absorbance values were used. The calculations of the pKBH+ values were performed from the absorbance data at four selected wavelengths around the absorption maximum that appears in neutral and acidic media. Following wavelengths were used: H1 (258, 266, 324 and 332 nm), H2 (242, 250, 324 and 332 nm), H3 (258, 266, 330 and 338 nm), H4 (242, 250, 332 and 324 nm) and H5 (258, 266, 332 and 340 nm). The obtained overdetermined system had this form (Eq. 3):
where, c(B) is the concentration of neutral form, c(BH+) is the concentration of protonated form, ε(B) is the molar absorption coefficient of neutral form, ελ4(BH+) is the molar absorption coefficient of protonated form. The obtained over determined system was calculated using Excel. The molar absorption coefficients were calculated in accordance with equation (4).
In order to determine the molar absorption coefficients, the absorbance values of the neutral (pH 4.5) and protonated (pH 1.4) form were used. For this purpose, the UV spectra were recorded at three different concentrations (2.40·10-5, 3.00·10-5 and 3.60·10-5 mol/dm3) of the investigated hydrazones in solution at ionic strengths 0.1 mol/dm3, 0.25 mol/dm3 and 0.5 mol/dm3.
The influence of the solvent can be eliminated by employing method of Characteristic Vector Analysis (CVA), which has been tested and proven to be applicable in the analysis of spectroscopic data, to investigate problems dealing with protonation [21,22]. The pK BH+ values of the investigated hydrazones were determined graphically, as an intercept of the dependence of logI on pH [29]. When c(BH+) = c(B), logI = 0, and the graphically determined pKBH+ value is equal to the pH value of the solution. The thermodynamic pKBH+ values were evaluated as an intercept with extrapolation of the curve to zero ionic strength [30].
Semiempirical calculations
There is a great interest in applying semiempirical methods to predict proton transfer at a given pH, for different types of organic molecules. Optimization of the geometry of investigated hydrazones and calculation of total energy, binding energy, enthalpy of formation, Gibbs energies of formation, atomic charge, proton affinity values was performed using AM1 (Austin Model 1) and PM3 (Parametric Method 3) semiempirical methods [23,31]. These parameters were further used for prediction of protonation site, influence of the substituents in the para position in the benzene ring of investigated hydrazones, as well as the stability of the isomers (E and Z) and their protonated forms. Proton affinity values were calculated according to following equations:
where, PA is the proton affinity, ΔHof (B)/ΔGf(B) is the heat/Gibbs energy of formation for the molecule, ΔHof(BH+)/ΔGf (BH+) is the heat/Gibbs energy of formation for the cation and ΔHof(H+)/ΔGf(H+) is the heat/Gibbs energy for the proton formation and has the value of 367.2 /362.57 kcal/mol, kcal/mol [32]. Theoretical calculations were carried out at the restricted Hartree-Fock level (RHF) in the HYPERCHEM program. All structures were optimized to a gradient norm of < 0.1.
Results and discussion
Spectral behavior of p-nitro-p-substituted benzoylhydrazones in acidic media
The experimental UV spectra of hydrazone N-p-nitrobenzaldehyde-p-etoxybenzoyl-hydrazone (H3) obtained in perchloric acid media in the pH region from 1.2 to 4.9 (water:ethanol = 1:1, V/V) maintaining the ionic strength of 0.1 mol/dm3 are shown in Fig. 1. Similar changes in the UV spectra of this hydrazone were noticed at 0.25 and 0.5 mol/dm3 ionic strength. The spectral behavior of the other investigated hydrazones (H1, H2, H4 and H5, see Table 1) at the same conditions was similar and their UV spectra are dropped out.
From the Fig. 1 it can be seen that in the weak acid media (pH 4.9), in the UV spectra of hydrazone H3 there are two absorption bands. The first band is located at 198 nm wavelength, while the second one has a maximum at 330 nm. According to the literature data the first absorption band is as a result of π→π* electron transitions, mainly in the aromatic structure of the molecules [33]. The second absorption band (330 nm) is probably due to the n→π* electron transitions [34]. Literature data revealed that when the nitrogen atom with lone electron pair is present as it is case with hydrazone nitrogen, than there is a huge probability of n→π electron transitions [35].
By increasing the acidity of the solution a hypochromic effect was noticed for the first absorption band (198 nm), while its position did not change. The first band has similar position for all investigated hydrazones suggesting that the substituents present in the investigated hydrazones have no important influence on this band. On the other hand the significant changes were observed in the intensity and the position of the second absorption band (330 nm). Namely, by increasing the acidity of the solutions the intensity of this band decreased and hypsochromic shift was observed. This absorption band disappeared at pH values less than 1.65. At the same time, in acidic media (pH < 2.77) a new absorption band at a wavelength of 260 nm appeared. The intensity of this band increased, by increasing the acidity of the solution (See Fig. 1). At pH values lower than 2 there were no further changes in the intensity and the position of this band. The described changes in the UV spectra of hydrazone H3 are due to the reaction of protonation in acidic media.
It can be noticed the existence of isosbestic points at 240 and 280 nm wavelength, in the UV spectra of the hydrazone H3 (Fig. 1), but they are not clearly defined. One reason for this can be the influence of the solvent taking into consideration that the protonation reaction is reversible i.e. the neutral and protonated form are in equilibrium in the solution. The influence of the solvent from the experimental spectra can be eliminated applying CVA method [21,22]. The reconstructed UV spectra of the hydrazone H3 at ionic strength 0.1 mol/dm3 are shown in the Fig. 2.
In the reconstructed spectra it can be observed well defined isosbestic point at 240 and 285 nm. The existence of isosbestic points additionally confirmed that in the solution there are two forms (neutral and protonated) of investigated hydrazones. The exact position of the absorption bands (neutral and protonated form) of investigated hydrazones (H1 -H5), as well as pH region of protonation is listed in Table 2.
The results presented in Table 2 demonstrated that by increasing the acidity of the solutions, the absorption band which appeared at 330 nm shifted from 66 nm (H1 and H5) to 84 nm (H2) i.e. towards shorter wavelengths. The pH region of protonation shows that there are no important differences in the acid-base behavior of investigated hydrazones in acidic media (Table 2). The dependence of the absorbance values on pH of the solution for investigated hydrazones (H1 - H5 ) at 0.1 mol/dm3 ionic strength are given in Fig. 3.
The dependence of the absorbance on pH of the solution has sigmoidal shape with one plateau indicating that the protonation reaction took place in one step. The upper plateau corresponds to the absorption of the neutral form which exists in the solution. As a result of protonation process the absorbance values decreased when the acidity of the solution increased. There were no noticed changes in the absorbance values, with the further changes in the acidity of the solution (pH ≤ 1.5). At this pH only the protonated form of the compounds is present in the solution. The substituents present in the molecule of the investigated hydrazones have not significant influence on the protonation process, presented with Scheme 1 [36]:
Semiempirical data
The influence of the nitro group present in the hydrazone molecule on the protonation site, was investigated by means of semiempirical calculation applying AM1 and PM3 methods. The obtained values of Gibbs energy of formation, enthalpy of formation and proton affinity values are listed in Table 3.
Comparing the proton affinity values calculated using AM1 and PM3 methods (see Table 3) it can be seen that the probable protonation site (higher proton affinity values) in the molecule of investigated hydrazone is sp2 hybridized nitrogen atom. The PA values are lower for hydrazone H4 in comparison with PA values of other investigated hydrazones. This indicated that presence of the chlorine atom in the hydrazone molecule (H4) causes a decrease in the affinity of sp2 nitrogen atom toward proton. The proton affinity values also suggested that the protonation reaction can be written according to Scheme 1.
In order to predict the stability of the Z and E isomers the values of total energy and proton affinity were compared. The calculation of proton affinity was made by the equation (6) given in the experimental part in accordance with AM1 semiempirical method. The obtained values are given in Table 4.
The total energy values demonstrated that the Z isomer has a little bit higher total energy values. This means that the Z isomer is slightly more stable in neutral media in comparison with the E isomer. The situation is different in acidic media where the cation is formed due to protonation of sp2 hybridized nitrogen atom. Namely, in acidic media there is no difference in the total energy values i.e. in the stability of the Z and E isomers. From the Table 4 it can be seen that the E isomer has higher PA values compared to Z isomer. This suggested that the E isomer can be easily protonated in acidic media.
Dissociation constants of protonated form, pKBH+
Stoichiometric pKBH+ values of N-p-nitrobenzaldehyde-p-substituted benzoylhydrazones (H1-H5) calculated from the absorbance values of the experimental and reconstructed spectra at ionic strength of 0.1, 0.25 and 0.5 mol/dm3 are listed in Tables 5 and 6, respectively. Statistical data (standard deviation (SD), relative standard deviation (RSD), coefficient of determination (R2)) and the range in which the obtained results are placed with confidence level of 0.05 (95%) are also given in Tables 5 and 6. The graphically obtained dissociation constants and the thermodynamic dissociation constants are presented in the Tables 5 and 6, as well. The dissociation constants values were obtained graphically, as an intercept from the dependence of logI on pH of the solution (Fig. 4).
The thermodynamic dissociation constant (table 3 and 4) were evaluated as an intercept from the dependence , with extrapolation to zero ionic strength [37, 38]. For illustration, the way of determination of thermodynamic dissociation constants for hydrazone H1 is given in Fig. 5.
The dependence can be used for water solution as well as for solution where the volume ratio of ethanol is 50 %.
There are no significant differences in the values of stoichiometric dissociation constants of all investigated hydrazones (H1 - H5). The calculated pKBH+ values are identical to those graphically obtained (See Tables 5 and 6). This means that dissociation constants can be successfully determined, in the both ways. This was also confirmed with the results of the t-test i.e. the value of t-parameter (-0.05) was below critical value (2.306) of this parameter. At the same time, there were no important differences in the pKBH+ values determined from the absorbance data of the experimental UV spectra with those determined from the absorbance data of the reconstructed spectra (See Tables 5 and 6). By increasing the value of the ionic strength, the pKBH+ values slightly increased which is in accordance with the theoretical expectations.
By comparison of thermodynamic dissociation constants it can be seen that the values ranged from 2.05 for hydrazone H4 to 2.55 for hydrazone H2. The following order of thermodynamic pKBH+ values was observed (Table 5):
The hydrazones H4 (−Cl) and H5 (−OH) have a little lower pKBH+ values i.e. theirs protonated forms are stronger acids compared with hydrazones H1 -H3. This is probably due to the influence of the present substituents in the molecule. Actually, the protonated form of hydrazone H4 which has chlorine atom in para position of benzene ring is the strongest acid compared with other investigated hydrazones. The weakest acid is hydrazone H2 (−CH3).
Conclusion
The spectral behavior of five p-nitro-p-substituted hydrazones was investigated in perchloric acid media, using the methods of UV spectroscopy. The changes in the UV spectra indicated that protonation reaction occurred in acidic media. The dependence of absorbance on pH indicated that the protonation process took place in one step. The influence of the solvent (50 % ethanol) was eliminated using the method of characteristic vector analysis (CVA). The changes in the UV spectra were used for calculation of dissociation constants. The pKBH+ values were determined numerically and graphically. In order to determine the thermodynamic dissociation constants the measurements were performed at three different ionic strengths (0.1, 0.25 and 0.5 mol/dm3). It should be emphasized that there are no big differences in the pKBH+ values of all investigated hydrazones suggesting that the substituents (−H, −CH3, −OCH3, −Cl and −OH) have no significant influence. The proton affinity values indicated that the probable protonation site in the molecule of investigated hydrazone is sp2 hybridized nitrogen atom, and the E isomer can be easily protonated. The Z isomer is slightly more stable in comparison with the E isomer in neutral media, while in acidic media there is no difference in stability of the Z and E isomers.
References
1. Brehme, R.; Enders, D.; Fernandez R.; Lassaletta, J. M. Eur. J. Org. Chem. 2007, 34, 5629-5660.
2. Belskaya, N. P.; Dehaen W.; Bakulev, V. A. ARKIVOC, 2010, (i), 275-332.
3. Banerjee, S.; Mondal, S.; Chakraborty, W.; Sen, S.; Gachhui, R.; Butcher, R. J.; Slawin, A. M. Z.; Mandal C.; Mitra, S. Polyhedron, 2009, 28(13), 2785-2793.
4. Mao, J.; Wang, Y.; Wan, B.; Kozikowski A. P.; Franzblau, S. G. Chem. Med. Chem. 2007, 2(11), 1624-1630.
5. Ali, Md. R.; Marella, A.; Alam, Md. T.; Naz, R.; Akhter, M.; Shaquiquzzaman, Md.; Saha R.; Tanwar, O.; Alam Md. M.; Hooda, J. Indonesian J. Pharm., 2012, 23, 193-202.
6. Kumar, P.; Rai, A.; Singh, M.; Kumar, D.; Sahdev A. K.; Raj V. EC Pharm. Sci. 2016, 23, 278 (2016).
7. Kumar, N.; Chauhan, L. S.; Dashora N.; Sharma, C. S. Sch. Acad. J. Pharm., 2014, 3(5), 366-373.
8. Padmini, K.; Jaya Preethi, P.; Divya, M.; Rohini, P.; Lohita, M.; Swetha, K.; Kaladar, P. Int. J Pharm. Res. & Rev. 2013, 2(8), 43-58.
9. Rollas S.; Küçükgüzel, S.G. Molecules, 2017, 12, 1910-1939.
10. Singh, R. B.; Jain P.; Singh, R. P. Talanta, 1982, 29, 77-84.
11. Suvarapu, L. N.; Seo1, Y. K.; Baek S. O.; Ammireddy, V. R. E-Journal of Chemistry, 2012, 9(3), 1288-1304.
12. Naskar, S.; Naskar, S.; Mondal, S.; Majhi, P. K.; Drew M. G. B.; Chattopadhyay, S. K. Inorg. Chim. Acta, 2011, 371, 100-106.
13. Singh M.; Raghav, N. Int. J. Pharm. Pharm. Sci. 2011, 3(4), 26-32.
14. Avila Terra, L. H. S.; da Cunha Areias, M. C.; Gaubeur, I.; Encarnacion, M.; Suarez-iha, V. Spectros. Lett. 1999, 32, 257-271.
15. Mohan, M.; Gupta, M. P.; Chandra L.; Jha, N. K. Inorg. Chim. Acta, 1988, 151, 61-68
16. Aggarwal, N.; Kumar, R.; Srivastva, C.; Dureja P.; Khurana, J. M. J. Agric. Food Chem. 2010, 58, 2056-3061.
17. Liu, M.; Wang, Y.; Wangyang, W. S. Z.; Liu, F.; Cui Y. L.; Duan, Y. S. J Agric Food Chem. 2010, 58, 6859-6863.
18. Pathare, B.; Tambe, V.; Dhole S.; Patil, V. Int. J. Pharm. 2014, 4(1), 278-285.
19. Beckett A. H.; Stenlake, J. B. Practical pharmaceutical chemistry. 4th Edition - part one, CBS publishers and distributors PVT LTD, New Delhi, 1997, 88-120.
20. Brahmankar D. M.; Jaiswal, S. B. Biopharmaceutics & pharmacokinetics. 2nd Ed., Vallabh prakashan, Delhi, 2009, 42-60.
21. Zalewski R. I.; Géribaldi, S. J. Chem. Soc., Perkin Trans. 1988, 2, 113-115.
22. Garcia, B.; Casado, R. M.; Castillo, J.; Ibeas, S.; Domingo, I.; Leal, J. M. J. Phys. Org. Chem., 1993, 6, 101-106.
23. Dewar, M. J. S.; Dieter, K. M. J. Am. Chem. Soc. 1986, 108, 8075-8086.
24. Steawart, J. J. P. J. Comput. Chem. 1989, 10, 209-216.
25. Ienascu, I. M. C.; Lupea, A. X.; Popescu, I. M.; Padure M. A.; Zamfir, A. D. J. Serb. Chem. Soc. 2009, 74, 847-855.
26. Rajput A. P.; Rajput, S. S. Int. J Pharm. Tech. Research, 2009, 1(4), 1605-1611.
27. Jankulovska, M.; Čolančeska-Rağenoviќ, K.; Dimova, V.; Spirevska I.; Makreski P. Org. Chem., An Ind. J. 2012, 8, 326-334.
28. Jankulovska, M.; Spirevska I.; Ragenoviќ, K. Č. Bull. Chem. Technol. Macedonia, 2006, 25, 29-37.
29. Davis C. T.; Geissman, T. A. J. Am. Chem. Soc. 1954, 76, 3507-3511.
30. E. J. King, Acid-Base Equilibria, Pergamon Press, Oxford, 1965, 1.
31. Dewar, M. J. S.; Zoebisch, E. G.; Healy E. F.; Stewart, J. J. P. J. Am. Chem. Soc. 1985, 107, 3902-3909.
32. Kireev, V. A. Methods of Practical Calculations in Thermodynamics of Chemical Reactions, Khimiya, Moscow, 1975, 1.
33. Kristallovich, E. L.; Eshimbetov, A. G.; Chuvylkin, V. D.; Belenkii L. I.; Shakhidoyatov, Kh. M. Chem. Nat. Compd. 2003, 39(5), 495-500.
34. Echevarria, A.; Nascimento, M.; Gerônimo, V.; Miller J.; Giesbrecht, A. J. Braz. Chem. Soc. 1999, 10(1), 60-64.
35. Amrallah, A. H.; Abdalla N. A.; Haty, E. El. J. Chin. Chem. Soc. 2006, 53, 697-706.
36. Jankulovska M. S.; Spirevska, I.; Maced. J. Chem. Chem. Eng. 2014, 33(1), 85-96.
37. Perisić-Janjić, N. U.; Lazarević, M.; Janjić, J.; Klisareva, Lj.; U. Scientist Phyl. Sciences. 1995, 7, 64-68.
38. Gaubeur, I.; Vincenza, M. R.; Iha, K.; Vázquez, M. E. S. I.; Eclética química. 2000, 25, 63-76.
Author notes
*Corresponding author: Mirjana S. Jankulovska e-mail: jankulovska_m@yahoo.com