THE THEORETICAL STUDIES OF MOLECULAR ORBITALS AND THERMODYNAMIC PROPERTIES OF ANTIFUNGAL DRUG
HTML Full TextTHE THEORETICAL STUDIES OF MOLECULAR ORBITALS AND THERMODYNAMIC PROPERTIES OF ANTIFUNGAL DRUG
Manju Verma, Anjali Pandey, Ritu, Manoj Kumar Verma and Joy Sarkar *
Department of Chemistry, University of Lucknow, Lucknow - 226007, Uttar Pradesh, India.
ABSTRACT: In this work, the molecular structure optimization, electronic properties like local chemical reactivity descriptors (Fukui function analysis for reactive site prediction), and thermodynamic properties (Heat capacity, entropy, zero points vibrational energy, thermal energy, translational, rotational and vibrational energy) of antifungal drugs were theoretically computed. In this work also, electronic structure, frontier molecular orbital energies (HOMO-LUMO gap for the chemical activity of drug), molecular electrostatic potential (useful for recognizing electrophilic and nucleophilic attack) and mulliken charge has been carried out. Atomic charges dipole moment, molecular polarizability, electronic parameter, and refractivity of molecular systems were calculated by the DFT method. The reactive site prediction and hydrogen bond interaction analysis was determined by using the Fukui function analysis. In this study, DFT/B3LYP/6-31+G and DFT/CAM-B3LYP/6-31+G (d, p) methods were used.
Keywords: |
Antifungal drugs (Miconazole) Thermodynamic properties, FMO energies, Fukui function
INTRODUCTION: Miconazole is an antifungal drug, in the market, it is known as Monistat; it is an antifungal medicine used for ringworm and yeast infections of the skin 1. This drug is used for ringworm of the body, feet (athlete's foot), and groin (jock itch) 1. The side effects of this drug include irritation or itchiness of the area in which it was applied. Use in pregnancy is believed to be safe for the baby 2. Miconazole is the imidazole of the medication family. The function of this drug to decrease the ability of fungi to make ergosterol 1. Miconazole was patented in 1968 and approved for medical use in 1971 3. It is on the World Health Organization's List of Essential Medicines, the most effective and safe medicines needed in a health system 4.
Computational Method: The Computational study was obtained by the DFT method with the hybrid functional B3LYP (A.D. Becke 1993) by using the 6-31G basis set. First, the drug is optimized then by optimized structure geometrical parameters are used for the calculation of energy (HOMO and LUMO), energy gap, molecular electrostatic potential (MESP) to characterize the stability of the molecule and reactive sites of the molecule.
RESULTS AND DISCUSSION:
Molecular Geometry: Fig. 1 shows the structure and atom numbering of Miconazole.
FIG. 1: OPTIMIZED STRUCTURE OF MICONAZOLE
Molecular Orbital Analysis: HOMO (highest occupied molecular orbital) is the highest energy orbital. The outermost electrons of orbital HOMO can act as an electron donor. LUMO is the innermost orbital. It is the lowest energy orbital. LUMO can act as electron acceptor 5. The transition of electrons of reacting species is due to an interaction between HOMO and LUMO; it is defined by FMO 6. More value of the energy of HOMO is showing the tendency of the molecule to donate electrons to an appropriate acceptor molecule of lower molecular orbital energy. The less values of energy of LUMO indicate more acceptances of electrons. The difference between the energy of HOMO and LUMO shows an important stability index of electronic transition. Low gap value indicates the higher electronic transition and vice versa. The energy of HOMO and LUMO are - 6.43227eV and -1.33826 eV respectively, and the energy gap (ΔE) is 5.09401 eV. The gap of energy shows the chemical activity of the compound. Analysis of Fig. 2 shows that the distribution of two energies HOMO and LUMO. We can see from Fig. 2 that the electron density of HOMO is mostly distributed over the -imidazole while in LUMO the charge density is mainly accumulated on dichlorophenyl ethyl.
FIG. 2: THE ATOMIC ORBITAL COMPOSITIONS OF THE FRONTIER MOLECULAR ORBITAL (HOMO AND LUMO) FOR MICONAZOLE
Molecular Electrostatic Potential: MESP refers to the electronic density, and it is very useful for recognizing the electrophilic and nucleophilic attack as well as hydrogen-bonding interactions 7-8. To know the active sites for electrophilic and nucleophilic attack for Miconazole, MESP is calculated by applying the DFT/B3LYP method and 6-31G (d) basis set for the optimized geometry. The negative (red) part of molecular electrostatic potential (MESP) represents electrophilic and the positive (blue) part for nucleophilic reactivity (as shown in Fig. 3. It provides a pictorial method to recognize the relative polarity of a molecule and help as useful information to interpret hydrogen bonding, reactivity, and structure-activity relation-ship of molecules. It is the potential energy of a proton at a definite area nearby a molecule. Various values of the electrostatic potential at the surface of a molecule referred to the various colors. In brief, the attractive potential (or negative) arises in red-colored part, and those of repulsive (or positive) potential appear in blue. As can be seen in Fig. 3, there are two possible sites on the compound for an electrophilic attack. The negative part majorly over the Nitrogen atom of Imidazole ring. The positive part is located on the ethyl group attached Imidazole ring in a very small amount (site for nucleophilic attack). The maximum values of negative and positive parts are -6.81 × 10-2 a.u and 6.81 × 10-2 a.u, respectively. These results refer that the combined study of the structural and electronic features of a molecule is useful for recognizing the structure- performance relationship, thus allowing for the design of new and better heterocyclic drugs with applications.
FIG. 3: 3D PLOT OF THE MOLECULAR ELECTRO-STATIC POTENTIAL OF MICONAZOLE
Thermodynamical Analysis: Statistical thermo-dynamic functions mainly heat capacity and entropy were calculated for the molecule at varying temperatures (100 to 500 K) and summarized in Table 1. The correlation graph between these thermodynamic measurements and temperatures (T) are shown in Figs. (a) and (b). The zero-point vibrational energy (ZPVEs), thermal energy, rotational constant, molar heat capacity, entropy, and enthalpy at room temperature for Miconazole were obtained and indexed in Table 2. It is obvious from our observations that the calculated ZPVE energy is lower in the B3LYP (177.60 kcal mol-1) than CAM-B3LYP 181.38 (kcal mol-1) method. However, the calculated molar heat and entropy were 85.548, 163.936, and 82.344, 161.192 cal mol-1 k-1 respectively at B3LYP and CAM-B3LYP/6-31G (d, p) hybrid functional.
FIG. 4: (A) & (B) CORRELATION GRAPHS OF HEAT CAPACITY AND ENTROPY CALCULATED AT VARIOUS TEMPERATURES USING B3LYPAND CAM-B3LYP/6-31G (D, P) OF MICONAZOLE
TABLE 1: THERMODYNAMIC FUNCTIONS AT DIFFERENT TEMPERATURE OF MICONAZOLE DRUG EMPLOYING B3LYP AND CAM-B3LYP/6-31 G (D, P) METHOD
Temperature(T) (K) | B3LYP/6-31-G (d,p) | CAM-B3LYP/6-31-G (d,p) | ||
Temperature | Heat capacity (CV) (Cal/mol K) | Entropy (S)
Cal/Mol K) |
Heat capacity (CV) (Cal/mol K) | Entropy (S)
Cal/Mol K) |
100 | 35.869 | 100.159 | 35 .869 | 99.074 |
200 | 61.025 | 134.191 | 59.037 | 132.509 |
298 | 85.548 | 163.936 | 82.344 | 161.192 |
300 | 86.001 | 164.479 | 86.001 | 164.479 |
400 | 108.940 | 193.014 | 108.940 | 193.014 |
500 | 127.944 | 219.886 | 124.914 | 215.441 |
TABLE 2: CALCULATED THERMODYNAMIC PARAMETERS OF MICONAZOLE EMPLOYING B3LYP AND CAM-B3LYP/6-31 G (d, p) METHODS-(AT ROOM TEMPERATURE)-
Parameters | B3LYP/6-31 G (d,p) | CAM-B3LYP/6-31G (d,p) | ||
Zero-point vibrational energy (Kcal/mol) | 177.60709 | 181.38253 | ||
Rotational temperatures (K) | 0.01073 | 0.01073 | ||
0.00384 | 0.00384 | |||
0.00329 | 0.00329 | |||
Rotational constants (GHZ) | ||||
X | 0.22361 | 0.22361 | ||
Y | 0.08003 | 0.08003 | ||
Z | 0.06847 | 0.06847 | ||
Thermal energy (Kcal/mol) Total | 191.773 | 195.187 | ||
Translational | 0.889 | 0.889 | ||
Rotational | 0.889 | 0.889 | ||
Vibrational | 189.996 | 193.410 | ||
Molar capacity at constant volume (cal mol/K) | ||||
Total | 85.548 | 82.344 | ||
Translational | 2.981 | 2.981 | ||
Rotational | 2.981 | 2.981 | ||
Vibrational | 79.586 | 76.383 | ||
Entropy (cal mol/K) | ||||
Total | 163.936 | 161.192 | ||
Translational | 43.953 | 43.953 | ||
Rotational | 36.815 | 36.815 | ||
Vibrational | 83.168 | 80.424 | ||
Local Reactivity Descriptors: Fukui function (FF) is a function that is related to electron density; it gives information regarding the local site reactivity of the molecule. Fukui function is related to qualitative descriptors of reactivity of different atoms in the molecule.
A study by Ayers and Parr 9-10 has shown that the greater value of FF gives information that attack is by soft reagents, and lesser value of FF gives information that attacked by hard reagents. It is calculated by the DFT method by using the B3LYP/6-31 G (d, p) basis set. With the use of Mulliken atomic charges of cationic and anionic states, local Fukui functions, local softness value, and local electrophilicity indices have been calculated using the following equation:
f k+ = [q (N+1) – q (N
(For nucleophilic attack)
fk- = [q (N) – q (N-1)]
(For electrophilic attack)
fk0 = [q (N+1) – q (N-1)]/2
Fukui functions for electrophilic and nucleophilic attacks have been calculated by (for free radical attack) Local softness values, and electrophilicity indices are calculated by using
sk+ =Sfk + sk- =Sfk-
ωk+= ωfk + ωk -= ωfk-
Where + and − signs show nucleophilic and electrophilic attack, respectively 11. Fukui functions, local softness values, and local electrophilicity for selected atomic sites in Miconazole have been listed in Table 3.
TABLE 3: FUKUI FUNCTIONS ( fk+, fk-), LOCAL SOFTNESS (sk+, sk-), AND LOCAL ELECTROPHILICITY INDICES (ωk+, ωk-) FOR MICONAZOLE, USING MULLIKEN POPULATION ANALYSIS AT B3LYP/6-311 G (d, p) LEVEL
1 C | 0.077958 | 0.124765 | -0.003145 | 0.046807 | 0.081103 | 0.012066 | 0.0209 | 0.171346 | 0.296894 |
2 C | 0.027475 | 0.103045 | 0.050806 | 0.07557 | -0.02333 | 0.01948 | -0.00601 | 0.276639 | -0.08541 |
3 C | -0.040165 | 0.014959 | -0.012937 | 0.055124 | -0.02723 | 0.01421 | -0.00702 | 0.201792 | -0.09967 |
4 C | -0.358275 | -0.276161 | -0.295442 | 0.082114 | -0.06283 | 0.021167 | -0.01619 | 0.300595 | -0.23001 |
5 C | 0.135254 | 0.171129 | 0.076714 | 0.035875 | 0.05854 | 0.009248 | 0.015086 | 0.131328 | 0.214297 |
6 C | -0.3382 | -0.26388 | -0.271881 | 0.07432 | -0.06632 | 0.019158 | -0.01709 | 0.272063 | -0.24277 |
10 C | 0.290993 | 0.35731 | 0.291385 | 0.066317 | -0.00039 | 0.017095 | -0.0001 | 0.242767 | -0.00143 |
11 O | -0.469877 | -0.48188 | -0.500405 | -0.012 | 0.030528 | -0.00309 | 0.007867 | -0.04394 | 0.111754 |
12 C | 0.275691 | 0.215638 | 0.179239 | -0.06005 | 0.096452 | -0.01548 | 0.024856 | -0.21984 | 0.353082 |
13 C | 0.447507 | 0.33353 | 0.250611 | -0.11398 | 0.196896 | -0.02938 | 0.05074 | -0.41724 | 0.720777 |
15 C | -0.33714 | 0.022479 | -0.004754 | 0.359619 | -0.33239 | 0.092703 | -0.08566 | 1.316457 | -1.21677 |
16 C | 0.107204 | 0.063483 | -0.011429 | -0.04372 | 0.118633 | -0.01127 | 0.030572 | -0.16005 | 0.43428 |
17 C | 0.073343 | 0.112247 | -0.000471 | 0.038904 | 0.073814 | 0.010029 | 0.019022 | 0.142416 | 0.270211 |
18 C | -0.334386 | -0.271268 | -0.297754 | 0.063118 | -0.03663 | 0.016271 | -0.00944 | 0.231056 | -0.1341 |
19 C | 0.142363 | 0.16912 | 0.065381 | 0.026757 | 0.076982 | 0.006897 | 0.019838 | 0.097949 | 0.281808 |
20 C | -0.280295 | -0.291753 | -0.294873 | -0.01146 | 0.014578 | -0.00295 | 0.003757 | -0.04194 | 0.053366 |
24 C | 0.4305 | 0.629336 | 0.487154 | 0.198836 | -0.05665 | 0.051256 | -0.0146 | 0.727879 | -0.20739 |
25 C | 0.120872 | 0.250926 | 0.089197 | 0.130054 | 0.031675 | 0.033525 | 0.008163 | 0.476089 | 0.115953 |
26 C | 0.176179 | 0.387165 | 0.290773 | 0.210986 | -0.11459 | 0.054388 | -0.02953 | 0.772356 | -0.41949 |
27 N | -0.515518 | -0.667325 | -0.658416 | -0.15181 | 0.142898 | -0.03913 | 0.036825 | -0.55572 | 0.523107 |
31 N | -0.413434 | -0.464999 | -0.529115 | -0.05157 | 0.115681 | -0.01329 | 0.029811 | -0.18876 | 0.423473 |
36 Cl | 0.173539 | 0.226719 | 0.024944 | 0.05318 | 0.148595 | 0.013709 | 0.038293 | 0.194676 | 0.543962 |
37 Cl | 0.210206 | 0.18331 | 0.044313 | -0.0269 | 0.165893 | -0.00693 | 0.042751 | -0.09846 | 0.607285 |
38 Cl | 0.215865 | 0.156889 | 0.031927 | -0.05898 | 0.183938 | -0.0152 | 0.047401 | -0.21589 | 0.673342 |
39 Cl | 0.182342 | 0.195217 | -0.001824 | 0.012875 | 0.184166 | 0.003319 | 0.04746 | 0.047132 | 0.674176 |
Mulliken Charges: Mulliken charges originate from the Mulliken population analysis. Effective atomic charge calculation plays a main role in the quantum mechanical calculations of molecular systems. Due to the change of atomic charges dipole moment, molecular polarizability, electronic parameter and refractivity of molecular systems can be calculated. Mulliken charges calculated by DFT at B3LYP /6-31G (d, p) and CAM B3LYP/6-31 basis set 12-13. It is calculated to describe the electron population of each atom and to see charge distributions over the molecule as listed in Table 4 and plotted in Fig. 5.
TABLE 4: THE MULLIKEN CHARGE DISTRIBUTION CALCULATED AT B3LYP AND CAM- B3LYP/6-31G (d, p) METHODS OF MICONAZOLE
Atom number | Atomic charges (Mulliken) | |
B3LYP/6-31
G (d,p) |
CAM-B3LYP/6-31 G (d,p) | |
1 C | 0.067417 | 0.053462 |
2 C | 0.069334 | 0.044961 |
3 C | 0.092061 | 0.064328 |
4 C | -0.286273 | -0.137781 |
5 C | 0.112025 | 0.093869 |
6 C | -0.231068 | -0.089708 |
10 C | 0.278123 | 0.313927 |
11 O | -0.537491 | -0.533968 |
12 C | 0.192637 | 0.241123 |
13 C | 0.259908 | 0.211375 |
15 C | 0.142166 | 0.111708 |
16 C | 0.012239 | 0.010598 |
17 C | 0.090934 | 0.071928 |
18 C | -0.235925 | -0.092115 |
19 C | 0.125731 | 0.105845 |
20 C | -0.287242 | -0.144263 |
24 C | 0.390729 | 0.365357 |
25 C | 0.064145 | 0.103537 |
26 C | 0.284624 | 0.184785 |
27 N | -0.595935 | -0.424533 |
31 N | -0.40215 | -0.450891 |
36 Cl | 0.083688 | -0.034883 |
37 Cl | 0.094883 | -0.031801 |
38 Cl | 0.114636 | -0.015504 |
39 Cl | 0.100804 | -0.021354 |
FIG. 5: MULLIKEN CHARGE DISTRIBUTION IN MICONAZOLE
CONCLUSION: In this theoretical study analysis of geometries and electronic properties of Miconazole is performed using DFT/B3LYP method at a 6-31 G (d,p) basis set. The effect of the substituted group on the structure and electronic properties are studied. The stabilization energy and the calculated HOMO and LUMO energies indicated charge transfer in the molecule, which in turn indicated its bioactive properties. Thus, the theoretical study indicates that these molecules are polar and active molecule, and they may interact with its environment strongly. The indications of the theoretical study reveal useful information about the reactivity of such molecules and give good information about the active sites in the molecules and clarify the sites of molecules that undergo nucleophilic substitution or electrophilic substitution reactions.
ACKNOWLEDGEMENT: Authors are thankful to Head of Department of chemistry, Lucknow University, Lucknow for providing laboratory, Spectral, and computational Facilities.
CONFLICTS OF INTEREST: The authors confirm that this article content has no conflicts of interest.
REFERENCES:
- Mapelli ET, Cerri A, Bombonato C and Menni S: Tinea capitis in the paediatric population in Milan, Italy: the emergence of Trichophyton violaceum. Mycopathologia 2013; 176: 243.
- Stephansson O, Granath F, Svensson T, Haglund B, Ekbom A and Kieler H: Drug use during pregnancy in Sweden - assessed by the prescribed drug register and the medical birth register. Clin Epidemiol 2011; 3: 43-50.
- Janos F, Robin GC, John Wiley & Sons, Analogue- based Drug Discovery 2006; 502.
- Wirtz VJ, Hogerzeil HV, Gray AL, Bigdeli M, de Joncheere CP, Ewen MA, et al. Essential medicines for universal health coverage. Lancet 2017; 389: 403-76.
- RM Issa, MK Awad and Atlam FM: “Quantum Chemical Studies on the Inhibition of Corrosion of Copper Surface by Substituted Uracils,” Applied Surface Science 2008; 255(5): 2433-41.
- Ashry ESHEl, El Nemr A, Esawy SA and Ragab S: “Corrosion Inhibitors: Part II: Quantum Chemical Studies on the Corrosion Inhibitions of Steel in Acidic Medium by Some Triazole, Oxadiazole and Thiadiazole Derivatives,” Electrochimica Acta 2006; 51(19): 3957-68.
- Murray PP and Clark JST: Halogen bonding and other σ-hole interactions: A Perspective Phys Chem Chem Phys 2013; 15: 11178-89.
- larcón-Espósito J, Contreras R and Campodónico PR: Iso-solvation effects in mixtures of ionic liquids on the kinetics of a model SNAr reaction. New J Chem 2017; 41: 13435-41.
- Ayers PW and Parr RG: “Variational principles for describing chemical reactions: the Fukui function and chemical hardness revisited,” Journal of the American Chemical Society 2000; 122(9): 2010-18.
- Gillis EP, Eastman KJ, Hill MD, Donnelly DJ and Meanwell NA: Applications of fluorine in medicinal chemistry. J Med Chem 2015; 58: 8315-59.
- Chattaraj PK, Maiti B and Sarkar UJ: Phys Chem A 2003; 107; 4973.
- Balachandran V, Karpagam V, Santhi G, Revathi B, Ilango G and Kavimani M: Spectro Acta Part A: Mol Biomol Spectrosc 2015; 137: 165.
- Prabhaharan M, Prabakaran AR, Gunasekaran S and Srinivasan S: Spectrochim. Acta Part A: Mol Biomol Spectrosc 2015; 136: 494.
How to cite this article:
Verma M, Pandey A, Ritu, Verma MK and Sarkar J: The theoretical studies of molecular orbitals and thermodynamic properties of antifungal drug. Int J Pharm Sci & Res 2020; 11(12): 6488-93. doi: 10.13040/IJPSR.0975-8232.11(12).6488-93.
All © 2013 are reserved by the International Journal of Pharmaceutical Sciences and Research. This Journal licensed under a Creative Commons Attribution-NonCommercial-ShareAlike 3.0 Unported License.
Article Information
65
6488-6493
747
517
English
IJPSR
M. Verma, A. Pandey, Ritu, M. K. Verma and J. Sarkar *
Department of Chemistry, University of Lucknow, Lucknow, Uttar Pradesh, India.
dr.jsarkar4@gmail.com
04 November 2019
06 May 2020
05 November 2020
10.13040/IJPSR.0975-8232.11(12).6488-93
01 December 2020