Journal of Computational Methods in Physics

Volume 2014 (2014), Article ID 623235, 6 pages

http://dx.doi.org/10.1155/2014/623235

Research Article

Comparative DFT Study of Phytochemical Constituents of the Fruits of Cucumis trigonus Roxb. and Cucumis sativus Linn.

Subarayan Bothi Gopalakrishnan,¹ Thangaraj Kalaiarasi,² and Ramasamy Subramanian³

¹Department of Chemistry, Noorul Islam University, Kumaracoil, Kanyakumari, Tamil Nadu 629 180, India

²Department of Pharmaceutical Chemistry, Manonmaniam Sundaranar University, Tirunelveli, Tamil Nadu 627152, India

³Centre for Scientific and Applied Research, School of Basic Engineering and Sciences, PSN College of Engineering and Technology, Tirunelveli, Tamil Nadu 627012, India

Received 12 June 2014; Revised 17 July 2014; Accepted 31 July 2014; Published 12 August 2014

Academic Editor: Xavier Ferrieres

Copyright © 2014 Subarayan Bothi Gopalakrishnan et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Abstract.

The hepatoprotective active phytochemical constituents from the ethanolic extracts of the fruits of Cucumis trigonus Roxb. and Cucumis sativus Linn. were identified by GC-MS analysis. The density functional theory (DFT) of these molecules was calculated by density functional B3LYP methods using B3LYP/6-311++G(d,p) basis set. The optimized geometries of phytochemical constituents were evaluated. Physicochemical properties such as HOMO, LUMO, ionization potential, electron affinity, electronegativity, electrochemical potential, hardness, softness, electrophilicity, total energy, and dipole moment have also been recorded. These are very important parameters to understand the chemical reactivity and biological activity of the phytochemical constituents. Glycodeoxycholic acid and 2-(2-methylcyclohexylidene)-hydrazinecarboxamide were found to be effective drugs selected on the basis of their HOMO and LUMO energy gap and softness. The effective properties of these compounds may be due to the presence of amino, carbonyl, and alcohol as a functional group.

1. Introduction

DFT is the connection between the density and the wave function which was developed by Hohenberg and Kohn and stated in the Hohenberg-Kohn theorems [1, 2]. The first theorem shows that the Hamiltonian is fully determined by the ground state density, which means that all the physical properties of the system are completely determined by the ground state density. The second theorem demonstrates that the energy functional can be used to find the ground state density. Importantly, the excited states need to be determined in a different way. During the last decade, computer simulations based on a quantum-mechanical description of the interaction between electrons and atomic nuclei have strongly influenced the development of material science [3]. The major advantage from a computational point of view is that the electronic density depends only on the three variables of space, while the wave function for an N electron system depends on 3N variables.

DFT uses the B3LYP/6-311++G(d,p) hybrid functional in the gas phase. The presence of phytochemical constituents in the ethanolic extracts of the fruits of Cucumis trigonus and Cucumis sativus is responsible for curing liver disorder. Silymarin is used as a potential hepatoprotective agent. In the present work, a comparative DFT study of phytochemical constituents identified from the ethanolic extracts of the fruits of Cucumis trigonus and Cucumis sativus has been performed. Gaussian 09W program and Gauss-View molecular visualization program package on the personal computer have been used.

2. Material and Methods

2.1. Collection of Plant Materials

The fruits of Cucumis trigonus and Cucumis sativus were collected in the month of March and July, respectively, from Alangulam, Tirunelveli District, Tamil Nadu, and identified by Professor P. Jayaraman, Plant Anatomy Research Centre, West Tambaram, Chennai, Tamil Nadu, India (Registration Numbers of the Authentification Certificate: PARC/2013/2048 and PARC/2013/2047).

The voucher specimens (MSU/PHAR/HER-140 and MSU/PHAR/HER-141) have been preserved in the Herbarium of the Department of Pharmaceutical Chemistry, Manonmaniam Sundaranar University, Tirunelveli, Tamil Nadu, India.

2.2. Identification of Photochemical Constituents

Nine phytochemical constituents, namely, demeclocycline, glycodeoxycholic acid, 3α,7α,12α-trihydroxycoprostanic acid, chlortetracycline, azafrin methyl ester, giganteumgenin N, phorbol 12,13-dihexanoate, astaxanthin, and tetrahydrospirilloxanthin, were identified from the ethanolic extract of the fruits of Cucumis trigonus. Six phytochemical constituents, 2,3-dihydro-3,5-dihydroxy-6-methyl-4H-pyran-4-one, 5-(hydroxymethyl)-2-furancarboxylic aldehyde, 4-hydroxy-3-methyl-2-butenyl-acetate, 2-(2-methylcyclohexylidene)-hydrazide carboxamide, n-hexadecanoic acid, and 1,2-benzenedicarboxylic acid-diisooctyl ester, were identified from the ethanolic extract of the fruits of Cucumis sativus using GC-MS [4, 5] analysis.

2.3. Computational Methods

All calculations were performed using Gaussian 09 software package [6], at the B3LYP/6-311++G(d,p) level of theory. Transitions to the lowest excited singlet electronic states of phytochemical constituents from the ethanolic extracts of the fruits of Cucumis trigonus and Cucumis sativus were computed by using the gradient corrected DFT with the three-parameter hybrid functional Becke3 (B3) for the exchange part and the Lee-Yang-Parr (LYP) correlation function, and HOMO-LUMO energy level calculations have been carried out in the present investigation, using 6-311++G(d,p) Basis set with Gaussian 09W program package [7–9].

The chemical reactivity descriptors were calculated using DFT. These are very important physical parameters to understand chemical and biological activities of the phytochemical constituents. The calculated HOMO-LUMO orbital energies can be used to estimate the ionization energy [10], electron affinity [10], electronegativity [11], electronic chemical potential [11], molecular hardness [12], molecular softness [12], and electrophilicity index [13] using the following equations:

 ionization energy (IE) = , electron affinity (EA) = , electronegativity () = (IE + EA)/2, electronic chemical potential , chemical hardness () = (IE− EA)/2, chemical softness , electrophilicity index .

3. Results and Discussion

DFT has long been used to predict certain functions such as ionization potentials [14], electron affinities [14–16], orbital energies [17], and molecular structures [18]. The frontier orbitals (HOMO and LUMO) of a chemical species are very important in defining its reactivity. Higher value of HOMO of a molecule has a tendency to donate electrons to appropriate acceptor molecule with low energy, empty molecular orbitals. , on the other hand, indicates the ability of the molecule to accept electrons. The binding ability of the molecule increases with increasing HOMO and decreasing LUMO energy values. Thus, the lower the value of is, the most probable the molecule would accept electrons. When the HOMO-LUMO energy gap () decreases, the reactivity of the molecule increases leading to a decrease in the stability of the molecule. HOMO-LUMO gap is generally the lowest energy electronic excitation that is possible in a molecule.

Chemical hardness is associated with the stability and reactivity of a chemical system. The larger the HOMO-LUMO energy gap, the harder and more stable and less reactive the molecule. The smaller the HOMO-LUMO energy gap, the softer and less stable and more reactive the molecule. Electronic chemical potential is the negative value of electronegativity of a molecule. The greater the electronic chemical potential, the less stable or more reactive, the molecule will be. Electrophilicity index () [13], measures the propensity or capacity of a species to accept electrons [19, 20]. It is a measure of the stabilization in energy after a system accepts additional amount of electronic charge from the environment [21]. More stable molecule corresponds to lesser electrophilicity.

The results of density functional theory of phytochemical constituents identified from the ethanolic extracts of the fruits of Cucumis trigonus and Cucumis sativus are presented in Tables 1 and 2, respectively. By DFT studies, the presence of maximum number of phytochemical constituents in the ethanolic extract of the fruits of Cucumis trigonus showed more stability and biological activity than the ethanolic extract of the fruits of Cucumis sativus. The values of LUMO-HOMO energy gap reflect the chemical activity of the molecule. The larger the HOMO-LUMO energy gap, the harder and more stable and less reactive the molecule [22, 23]. The decrease in the HOMO and LUMO energy gap explains the eventual charge transfer interaction taking place within the molecule due to the strong electron-accepting ability of the electron acceptor group. The strong charge transfer interaction is responsible for the bioactivity of the molecule.

Table 1: Global chemical reactivity indices of phytochemical constituents of the ethanolic extract of the fruits of Cucumis  trigonus Roxb.

Serial number
Compounds
HOMO

(eV)
LUMO

(eV)
Energy

Gap
Ionization

potential (IP)

(eV)
Electron affinity (EA)

(eV)
Electronegativity ()  

(eV)
Electrochemical potential ()  

(eV)
Hardness ()  

(eV)
Softness ()  

(eV)
Electrophilicity ()  

(eV)
Total energy

(A.U)
Dipole moment

(Debye)

1
Demeclocycline
−0.2016
−0.0763
0.1253
0.2016
0.0763
0.1389
−0.1389
0.0626
15.9630
0.1540
−1973.77
3.4142

2
Glycodeoxycholic acid
−0.0827
−0.0578
0.0249
0.0827
0.0578
0.0703
−0.0703
0.0125
80.3213
0.1982
−2054.06
9.7149

3
3,7,12-Trihydroxycoprostanic

acid
−0.2252
−0.0118
0.2134
0.2252
0.0118
0.1185
−0.1185
0.1067
9.3721
0.0658
−1424.71
2.8593

4
Chlortetracycline
−0.1069
−0.0681
0.0388
0.1069
0.0681
0.0875
−0.0875
0.0194
51.5464
0.1973
−1796.74
4.4521

5
Azafrin

methyl ester
−0.2032
−0.0606
0.1426
0.2032
0.0606
0.1319
−0.1319
0.0713
14.0252
0.1220
−1384.08
3.3121

6
Giganteumgenin N
−0.2324
−0.0047
0.2277
0.2324
0.0047
0.1186
−0.1186
0.1139
8.7835
0.0617
−1439.54
3.6557

7
Phorbol

12,13-dihexanoate
−0.0718
−0.1406
−0.0688
0.0718
0.1406
0.1062
−0.1062
−0.0344
−29.0698
−0.1639
−1837.52
4.6762

8
Astaxanthin
−0.1916
−0.0667
0.1249
0.1916
0.0667
0.1292
−0.1292
0.0625
16.0128
0.1335
−1845.98
3.1631

9
Tetrahydrospirilloxanthin
−0.1801
−0.0509
0.1292
0.1801
0.0509
0.1155
−0.1155
0.0646
15.4799
0.1033
−1779.20
2.5202

10
Silymarin (std.)
−4.9624
−1.6207
3.3417
4.9624
1.6207
3.2915
−3.2916
1.6708
0.5985
3.2422
−1717.58
8.4093

Table 2: Global chemical reactivity indices of phytochemical constituents of the ethanolic extract of the fruits of Cucumis  sativus Linn.

Serial number
Compounds
HOMO

(eV)
LUMO

(eV)
Energy

gap
Ionization

potential (IP)

(eV)
Electron

affinity (EA)

(eV)
Electronegativity () 

(eV)
Electrochemical

potential () 

(eV)
Hardness ()  

(eV)
Softness () 

(eV)
Electrophilicity () 

(eV)
Total

energy (A.U)
Dipole

moment

(Debye)

1
2,3-Dihydro-3,5-dihydroxy-6-methyl-4H-pyran-4-one
−5.9894
−1.5130
4.4764
5.9894
1.5130
3.7512
−3.7512
2.2382
0.4468
3.1435
−534.14
1.8267

2
5-(Hydroxymethyl)-2-furancarboxaldehyde
−6.3997
−1.4808
4.9189
6.3997
1.4808
3.9402
−3.9403
2.4594
0.4066
3.1563
−457.72
5.3055

3
4-Hydroxy-3-methyl-2-butenyl-acetate
−6.6793
−0.6219
6.0574
6.6793
0.6219
3.6506
−3.6506
3.0287
0.3302
2.2001
−498.80
3.7023

4
2-(2-Methylcyclohexylidene)-hydrazinecarboxamide
−4.1393
−1.6646
2.4747
4.1393
1.6646
2.9019
−2.9020
1.2373
0.8082
3.4030
−1937.65
9.5533

5
n-Hexadecanoic acid
−6.9537
−0.5192
6.4345
6.9537
0.5192
3.7364
−3.7365
3.2172
0.3108
2.1697
−779.25
4.7330

6
1,2-Benzenedicarboxylic acid-diisooctyl ester
−5.2598
−0.8407
4.4191
5.2598
0.8407
3.0502
−3.0503
2.2095
0.4526
2.1054
−553.19
6.0504

7
Silymarin (std.)
−4.9624
−1.6207
3.3417
4.9624
1.6207
3.2915
−3.2916
1.6708
0.5985
3.2422
−1717.58
8.4093

Among the phytochemical constituents isolated and identified from the ethanolic extract of the fruits of Cucumis trigonus, glycodeoxycholic acid (A) showed more stability and biological activity as it shows less energy gap, low hardness, and more softness when compared with the standard silymarin (C). The energy gap, hardness, and softness of glycodeoxycholic acid are 0.0249, 0.0125, and 80.3213, respectively.

In the case of the phytochemical constituents isolated and identified from the ethanolic extract of the fruits of Cucumis sativus, 2-(2-methylcyclohexylidene)-hydrazinecarboxamide (B) shows less energy gap (2.4747), low hardness (1.2373), and more softness (0.8082) when compared with those of standard silymarin (C). DFT optimized geometries and frontier molecular orbitals of phytochemical constituents using B3LYP/6-311++G(d,p) basis set are presented in Figure 1.

Figure 1: DFT optimized geometries and frontier molecular orbitals of compounds (A), (B), and (C).

From the theoretical study, it has been found that the biological activities of the phytochemical constituents isolated and identified from the ethanolic extract of the fruits of Cucumis trigonus were found to be more than that of the phytochemical constituents isolated and identified from the ethanolic extract of the fruits of Cucumis sativus. The activities were found to be more than that of the standard drug, silymarin. This is supported by the experimental studies carried out by our group earlier [24–26].

Administration of paracetamol to rats by oral route caused liver damage as indicated by a significant increase in serum enzymes. Treatment with the ethanolic extracts of the fruits of Cucumis trigonus and Cucumis sativus decreased the levels of elevated serum enzymes like serum glutamic oxalacetic transaminase (SGOT), serum glutamic pyruvate transaminase (SGPT), serum alkaline phosphatase (SALP), γ-glutamate transpeptidase (GGTP), total bilirubin, conjugated bilirubin, and unconjugated bilirubin and increased the levels of total protein, albumin, and globulin. In 500 mg/kg body weight dose of the ethanolic extract of the fruits of Cucumis trigonus treated groups, the levels of SGOT, SGPT, SALP, GGTP, total bilirubin, conjugated bilirubin, unconjugated bilirubin, total protein, albumin, and globulin were considerably reduced even lower than that of the Cucumis sativus treated groups and normal rats. The lowering of enzyme level is a definite indication of hepatoprotective action of the drug. The results indicate that the 500 mg/kg body weight of the ethanolic extract of the fruits of Cucumis trigonus possesses more hepatoprotective activity than the standard drug, silymarin.

The antimicrobial effects of Cucumis trigonus and Cucumis sativus fruit extracts were tested against the different strains. In Cucumis trigonus, the benzene extract was effective in the gram-negative bacteria and fungi. The ethanol extract showed the maximum activity against all the tested strains followed by water extract. In Cucumis sativus, there was no significant antibacterial activity of the pet ether, benzene, and chloroform extracts against the test organisms. The ethanol extract and water extract of the fruits of Cucumis sativus only possess significant inhibitory activities against all the gram-positive, gram-negative bacteria and fungi. The difference in the activity may be due to the different secondary metabolites present in the various extracts. From these results it may be concluded that the fruits of Cucumis trigonus showed better activity than the Cucumis sativus.

The biological activities such as hepatoprotective activity [24, 25] and antimicrobial activity [26] of phytochemical constituents identified from the ethanolic extract of the fruits of Cucumis trigonus were more than that of the phytochemical constituents identified from the ethanolic extract of the fruits of Cucumis sativus.

4. Conclusion

Chemical reactivity descriptors such as HOMO-LUMO energy cap, softness, total energy, chemical hardness, electronic chemical potential, and electrophilicity of bioactive compounds have been examined theoretically. From the results it may be concluded that phytochemical constituents of the fruits of Cucumis trigonus and Cucumis sativus exhibited more biological activity when compared to that of the standard drug, silymarin. Based on the analysis of geometries, electronic structures, the biological activity of phytochemical constituents of the ethanolic extract of the fruits of Cucumis trigonus was found to be higher than that of phytochemical constituents of the ethanolic extract of the fruits of Cucumis sativus. The results indicate that the density functional theory is a valuable tool for predicting biological activity of phytochemical constituents. The HOMO-LUMO energy gap reveals the stability index of compounds.

Conflict of Interests

The authors declare that there is no conflict of interests regarding the publication of this paper.

References

K. Fukui, T. Yonezawa, and H. Shingu, “A molecular orbital theory of reactivity in aromatic hydrocarbons,” Journal of Chemical Physics, vol. 20, no. 4, pp. 722–725, 1952.

L. H. Mendoza-Huizar and C. H. Rios-Reyes, “Chemical reactivity of Atrazine employing the Fukui function,” Journal of the Mexican Chemical Society, vol. 55, no. 3, pp. 142–147, 2011.

J. Hafner, “Materials simulations using VASP—a quantum perspective to materials science,” Computer Physics Communications, vol. 177, pp. 6–13, 2007.

S. Gopalakrishnan and T. Kalaiarasi, “Identification of chemical compounds from the fruits of Cucumis trigonus Roxb. by GC-MS analysis,” International Journal of Phytopharmacy, vol. 2, pp. 122–128, 2012.

S. Gopalakrishnan and T. Kalaiarasi, “Hepatoprotective activity of the fruits of Cucumis sativus Linn,” International Journal of Pharmaceutical Sciences Review and Research, vol. 20, no. 2, pp. 229–234, 2013.

M. J. Frisch, G. W. Trucks, H. B. Schlegel, et al., Gaussian 09, Revision A.01, Gaussian, Wallingford, Conn, USA, 2009.

A. D. Becke, “Density-functional thermochemistry. II. The effect of the Perdew-Wang generalized-gradient correlation correction,” The Journal of Chemical Physics, vol. 97, no. 12, pp. 9173–9177, 1992.

A. D. Becke, “Density-functional thermochemistry. III. The role of exact exchange,” The Journal of Chemical Physics, vol. 98, no. 7, pp. 5648–5652, 1993.

C. Lee, W. Yang, and R. G. Parr, “Development of the Colle-Salvetti correlation-energy formula into a functional of the electron density,” Physical Review B: Condensed Matter, vol. 37, no. 2, pp. 785–789, 1988.

J. B. Foresman and A. Frisch, Exploring Chemistry with Electronic Structure Methods, Gaussian, Pittsburg, Pa, USA, 1995.

L. Pauling, The Nature of the Chemical Bond, Cornell University Press, Ithaca, NY, USA, 1960.

P. Senet, “Chemical hardnesses of atoms and molecules from frontier orbitals,” Chemical Physics Letters, vol. 275, no. 5-6, pp. 527–532, 1997.

R. G. Parr and R. G. Pearson, “Absolute hardness: companion parameter to absolute electronegativity,” Journal of the American Chemical Society, vol. 105, no. 26, pp. 7512–7516, 1983.

C. Zhan, J. A. Nichols, and D. A. Dixon, “Ionization potential, electron affinity, electronegativity, hardness, and electron excitation energy: Molecular properties from density functional theory orbital energies,” The Journal of Physical Chemistry A, vol. 107, no. 20, pp. 4184–4195, 2003.

J. C. Rienstra-Kiracofe, G. S. Tschumper, H. F. Schaefer III, S. Nandi, and G. B. Ellison, “Atomic and molecular electron affinities: photoelectron experiments and theoretical computations,” Chemical Reviews, vol. 102, no. 1, pp. 231–282, 2002.

J. C. Rienstra-Kiracofe, C. J. Barden, S. T. Brown, and H. F. Schaefer III, “Electron affinities of polycyclic aromatic hydrocarbons,” The Journal of Physical Chemistry A, vol. 105, no. 3, pp. 524–528, 2001.

J. Muscat, A. Wander, and N. M. Harrison, “On the prediction of band gaps from hybrid functional theory,” Chemical Physics Letters, vol. 342, no. 3, pp. 397–401, 2001.

L. A. Curtiss, P. C. Redfern, K. Raghavachari, and J. A. Pople, “Gaussian-3X (G3X) theory: use of improved geometries, zero-point energies, and Hartree-Fock basis sets,” Journal of Chemical Physics, vol. 114, no. 1, pp. 108–117, 2001.

A. Vektariene, G. Vektaris, and J. Svoboda, “A theoretical approach to the nucleophilic behavior of benzofused thieno[3,2-b]furans using DFT and HF based reactivity descriptors,” Arkivoc, vol. 2009, no. 7, pp. 311–329, 2009.

R. G. Parr, L. V. Szentpály, and S. Liu, “Electrophilicity index,” Journal of the American Chemical Society, vol. 121, no. 9, pp. 1922–1924, 1999.

S. Liu, “Dynamic behavior of chemical reactivity indices in density functional theory: a Bohn-Oppenheimer quantum molecular dynamics study,” Journal of Chemical Sciences, vol. 117, no. 5, pp. 477–483, 2005.

P. K. Chattaraj and B. Maiti, “HSAB principle applied to the time evolution of chemical reactions,” Journal of the American Chemical Society, vol. 125, no. 9, pp. 2705–2710, 2003.

K. D. Sen and D. M. P. Mingos, Chemical Hardness: Structure and Bonding, Springer, Berlin, Germany, 1993.

S. Gopalakrishnan and T. Kalaiarasi, “Hepatoprotective activity of the ethanolic extract of the fruits of Cucumis trigonus Roxb,” International Journal of Pharmacy and Pharmaceutical Sciences, vol. 5, no. 2, pp. 268–272, 2013.

S. Gopalakrishnan and T. Kalaiarasi, “Hepatoprotective activity of the fruits of Cucumis sativus (L.),” International Journal of Pharmaceutical Sciences Review and Research, vol. 20, no. 2, pp. 229–234, 2013.

S. Gopalakrishnan, T. Kalaiarasi, and R. Rajameena, “Evaluation of antimicrobial activity of the fruits of Cucumis trigonus Roxb,” Research Journal of Pharmacy, vol. 3, no. 5, pp. 256–258, 2012.