Comparison Chosen DFT Methods with Six Basis Sets on Harmine Molecule

Journal of the Korean Chemical Society.
2013.
Feb,
57(1):
20-24

- Received : July 25, 2012
- Accepted : November 23, 2012
- Published : February 20, 2013

Download

PDF

e-PUB

PubReader

PPT

Export by style

Article

Metrics

Cited by

TagCloud

The effect of the DFT methods and basis sets on harmine molecule has been investigated. 26 DFT methods with 6 basis sets were used. Two main comparison chosen DFT methods in this work has been investigated. It is concluded that this contribution is very important and cannot be neglected. In the following analysis, changes in energy levels were investigated by two different methods. Considering a specific basis set, changes in energy levels were obtained using different DFT methods. A specific DFT method is chosen and changes in energy levels have been investigated by means of various basis sets. Effect of the choice of basis sets on geometrical parameters on harmine molecule has also been investigated.
^{1}
With this theory, the properties of a many-electron system can be determined by using functional, i.e. functions of another function, which in this case are the spatially dependent electron densities. Hence the name density functional theory comes from the use of functional of the electron density. DFT is among the most popular and versatile methods available in condensed-matter physics, computational physics, and computational chemistry. DFT has been very popular for calculations in solid-state physics since the 1970s. However, DFT was not considered accurate enough for calculations in quantum chemistry until the 1990s, when the approximations used in the theory were greatly refined to improve the exchange and correlation interactions modeling. In many cases the results of DFT calculations for solidstate systems agree quite satisfactorily with experimental data. Computational costs are relatively low when compared to traditional methods, such as Hartree–Fock theory and its descendants based on the complex many-electron wave function. Since there are a large variety of DFT methods, choosing appropriate method and basis set are of crucial importance especially when dealing with large biological molecules with hetero-aromatic rings. Previous studies show that select and use of appropriate methods and computational tasks and basis sets are related to the structure under investigation.
^{12,23}
The aim of this work is the comparison of a series of different DFT methods and basis sets, in order to find a more suitable method and basis set for similar investigation. For this purpose we chose a medium size molecule, harmine (
. 1
), (7-Methoxy-1-methyl-9H-pyrido [3,4-b] indole) that is comparable with biological monomers such as nucleotides and amino acids.
^{4}
Harmine originally isolated from the seeds of Peganum harmala, is a tricyclic compound belonging to the β-carboline alkaloids. β-carboline such as harmine and harmaline constitute over than 60% of the alkaloids present in the plant seeds.
Scheme of structure of the Harmine molecule.
These alkaloids possess a broad range of pharmacological activities, such as anxiolytic and behavioral effects.
^{5,6}
Some studies demonstrated that harmine possessed significant anti-tumor potential both in vitro and in vivo.
^{7−9}
Anti-bactericidal
^{10}
and anti-inflammation.
^{11,12}
Also harmine has effects on central nervous system,
^{13,14}
bone system
^{15}
and enzyme systems.
^{16,17}
^{18}
Major computations were carried out using different DFT methods in combination with different basis sets. Also a MP2 computations under 6−31+G* and 6−31G* basis sets was carried out for comparison. For all computations the starting geometry of harmine was the same and all computations were carried out in gas phase.
A brief introduction about the different DFT methods is given in the
1
. More comprehensive details can be found in.
^{19−23 }
According to
1
the DFT methods used in this study can be divided into four groups on the basis of their exchange functional. They are hybrid,
^{24 }
Beckes 88 exchange functional,
^{25}
the Xa exchange functional,
^{26}
and Slaters exchange functional,
^{27}
based DFT methods. These four groups have a few different DFT methods that are derived from the three correlation functional of Perdew 86 (P86),
^{28}
Perdew-Wang 91(P91),
^{29}
and Lee-Yang-Parr (LYP).
^{30}
They are also combined with Perdew Local (PL),
^{30}
and Vosko-Wilk-Nusair (VWN)
^{31}
correlation functional. Since the purpose of this calculation on harmine molecule is to compare DFT methods with different basis sets, 26 DFT methods that mentioned in
1
were chosen along with 6 different basis sets to calculate the exact level of optimized energy. We used different basis sets for all the atoms. A final refinement of the basis set implies the introduction of diffuse functions
^{32}
for heavy atoms which are supplemented with diffuse s- and p-type functions which leads to 6−31+G* and 6−31G** sets. A higher degree of splitting is possible with the 6−311G basis set producing the resulting 6−311G* and 6−311G** sets.
Optimized harmine energy (J) obtained from MP2 computations
Full list of the DFT methods available in Gaussian 03
Optimized Harmine energy (J) obtained from DFT methods and basis sets
Bond angels:
Bond angels of harmine molecule were investigated using 26 DFT methods and basis sets. Proposed angels were 10C−11N−12C and 8C−7N−4C (
. 1
). The bond angels show small dependency in the choice of the basis sets for smaller ring. Results for C−N−C bond angels are illustrated in Figures 9-34(result are presented in supplementary data file). It shows very small discrepancy between all the basis sets. Results obtained for 6−31G*, 6−31G**, 6−311G* and 6−311G** are very similar and the results for 6−31G and 6−311G are exactly equivalent.
Bond length:
Bond lengths of harmine molecule are investigated using 26 DFT methods and 6 basis sets. The bond length of 11N−12C, 10C−11N, 7N−4C and 8C−7N for basis sets comparison is proposed (
. 1
). Results obtained for all 28 DFT methods are very similar, but the results showed a very different behavior for PW91PL, PW91P86, PW91VWN, PW91LYP, XaPL and XaVWN. The minimal basis sets (6−31G and 6−311G) showed a maximum length of C−N bonds for both rings. Figures 35 to 63 (result are presented in supplementary data file).
harmine optimized energy (J) obtained from DFT methods with six basis sets.
harmine molecule optimized energy using by Beckes combination functional on basis sets.
harmine molecule optimized energy using by G96 combination functional on basis sets.
harmine molecule optimized energy using by MPW combination functional on basis sets.
harmine molecule optimized energy using by PW91 combination functional on basis sets.
harmine molecule optimized energy using by S exchange functional on basis sets.
harmine molecule optimized energy using by Xa exchange functional on basis sets.
MP2 Calculation:
The optimized energy of harmine molecule was calculated using MP2 method with 6-31G* and 6-31+G* basis sets. Results are shown in
2
.
Our MP2 results show that there is no significant difference between MP2 and difference DFT energies. In other words; different DFT methods have the same accuracy as MP2 in energy calculations. So that DFT methods are preferred, because of some other advantages such as computation time.
Kurita N.
,
Ivanov P. M.
2000
Correlated Ab initio Molecular Orbital (MP3, MP4) and Density Functional (PW91, MPW91) Studies on the Conformations
J. Mol. Struct.
554
184 -

INTRODUCTION

Density functional theory (DFT) is a quantum mechanical modeling method used in physics and all fields of chemistry to investigate the electronic structure (principally the ground state) of many-body systems, in particular atoms, molecules, and the condensed phases.
PPT Slide

Lager Image

METHOD AND THEORY

Gaussian 03 program was used for all quantum chemical computations performed in this study.
Optimized harmine energy (J) obtained from MP2 computations

PPT Slide

Lager Image

Full list of the DFT methods available in Gaussian 03

PPT Slide

Lager Image

Optimized Harmine energy (J) obtained from DFT methods and basis sets

PPT Slide

Lager Image

RESULTS AND DISCUSSION

The purpose of this calculation on harmine molecule is to compare DFT methods with different basis sets. For this purpose, 26 DFT methods mentioned in
1
were chosen along in combination with 6 different basis sets to calculate the exact level of optimized energy. The results are shown in
3
. According to these results, the minimum energy is related to SLYP followed by SP86, SPL, SVWN5 and XaLYP respectively. On the other hand, the maximum energy is related to BVWN followed by MPWVWN, G96VWN, PW91VWN and G96PL, respectively. Chosen DFT methods listed in
2
were investigated in two different manners which yield two sets of results. DFT methods that show linear increase in rate of energy with a specific basis set are G96P86, G96LYP, PW91P86, PW91LYP and MPWLYP, which are illustrated in plot 2. On the other hand, there are DFT methods that was investigated with each basis set one by one. Based on plots 3 to 8 similar results were obtained for each method. In this manner the minimum energy level was related to basis set 6−31G, and the maximum energy level was related with higher degree of splitting such as 6−311G* and 6−311G**. Different behavior of basis set 6−31G* is noticeable in these plots.
- Effect of the Basis Sets on Geometrical Parameters

PPT Slide

Lager Image

PPT Slide

Lager Image

PPT Slide

Lager Image

PPT Slide

Lager Image

PPT Slide

Lager Image

PPT Slide

Lager Image

PPT Slide

Lager Image

CONCLUSION

Optimum energy level of harmine molecule is compared using 26 DFT methods and 6 basis sets (
2
). Results estimate that the minimum energy occurred while using Slater functional (SLYP, SP86 and SPL, Respectively) and the maximum energy occurred while using combination functional (PW91VWN, G96VWN, BVWN and MPWVWN, respectively).
1
and
2
and figures show the different rate of energy level changes for LYP, P86, PL and VWN, with exchange or combination functional. It is concluded that the increase in exchange functional energy level varies with correlation functional and it yields same results for every exchange functional. The maximum energy level was obtained from exchange functional with VWN and PL, P86 correlation functional, respectively. On the other hand, the minimum energy level is obtained with LYP correlation functional. It appears from these results that for compression chosen DFT methods in this work harmine molecule do not sufficient. For more accurate results in investigation of methods with different basis sets, we recommend using bigger bio molecules specially proteins and nucleotides acids.
Capelle K.
2006
A Bird’s-Eye View of Density-Functional Theory
Braz. J. Phy.
36
(4A)

Burke K.
2012
Perspective on Density Functional Theory
J. Chem. Phy.
136
150901 -
** DOI : 10.1063/1.4704546**

Nedelec J. M.
,
Hench L. L.
2000
Effect of Basis Set and of Electronic Correlation on Silica Rings
J. NonCryst. Solids
277
106 -
113
** DOI : 10.1016/S0022-3093(00)00306-9**

Murray T. D.
,
Berger A.
1997
Alcohol Withdrawal
Va. Med. Q
124
184 -
187

Brierley D. I.
,
Davidson C.
2012
Developments in Harmine Pharmacology Implications for Ayahuasca Use and Drug-dependence Treatment
Progress in Neuro-Psychopharmacology & Biological Psychiatry
39
263 -
272
** DOI : 10.1016/j.pnpbp.2012.06.001**

Almbladh C. O.
,
Barth U. V.
1985
Phys. Rev.
1331
3231 -

Perez Martin J. M.
,
Labrador V.
,
Fernandez Freire P.
,
Molero M. L.
,
Hazen M. J.
2004
Ultrastructural Changes Induced in HeLa Cells After Phototoxic Treatment with Harmine
J. Appl. Toxicol.
24
197 -
201
** DOI : 10.1002/jat.972**

Levy M.
1982
Phys. Rev.
A26
1200 -

Kohnand W.
,
Sham L.
1965
Phys. Rev.
140
A1133 -
** DOI : 10.1103/PhysRev.140.A1133**

Ashad N.
,
Zitterl-Eqlseer K.
,
Hasnain S.
,
Hess M.
2008
Effect of Peganum Harmala or Its Beta-carboline Alkaloids on Certain Antibiotic Resistant Strains of Bacteria and Protozoa from Poultry
Phytother. Res.
22
(11)
1533 -
1538
** DOI : 10.1002/ptr.2528**

Shahverdi A. R.
,
Monsef-Esfahani H. R.
,
Nickavar B.
,
Bitarafan L.
,
Khodaee S.
,
Khoshakhlagh N.
2005
Antimicrobial Activity and Main Chemical Composition of two Smoke Condensates from Peganum Harmala Seeds
Z. Naturforsch., C: J. Biosci.
60
(9-10)
707 -

Fortunato J. J. et al.
2009
A Cute Harmine Administration Induces Antidepressive-like Effects and Increases BDNF Levels in the Rat Hippocampus
Prog. Neuropsychopharmacol Biol. Psychiatry
33
(8)
1425 -
1430
** DOI : 10.1016/j.pnpbp.2009.07.021**

Farzin D.
,
Salimi E.
2009
Effects of Harmane, Norharmane and Harmine on Apomorphine-induced Pecking Behavior in Chick
Mazandaran Uni. Med. Sciences
19
(70)
747 -

Yonezawa T.
2011
αβ-carboline Alkaloid, Inhibits Osteoclast Differentiation and Bone Resorption in Vitro and in Vivo
Eur J. Pharmacol
650
(2-3)
511 -
518
** DOI : 10.1016/j.ejphar.2010.10.048**

Adayev T.
,
Wegiel J.
,
Hwang Y. W.
2011
Harmine is an ATPCompetitive Inhibitor for Dual-Specificity Tyrosine Phosphorylation- Regulated Kinase 1A [Dyrk1A]
Arch. Biochem. Biophys.
507
(2)
212 -
218
** DOI : 10.1016/j.abb.2010.12.024**

Heshmati E.
,
Abdolmaleki A.
,
Mozdarani H.
,
Sarvestani A. S.
2009
Effects of Halogen Substitution on Watson-Crick Base Pairing: a Possible Mechanism for Radiosensitivity
Bioorg. Med. Chem. Lett.
19
(17)
5256 -
5260
** DOI : 10.1016/j.bmcl.2009.06.105**

Göckler N.
,
Jofre G.
,
Papadopoulos C.
,
Soppa U.
,
Tejedor F. J.
,
Becker W.
2009
Harmine Specifically Inhibits Protein Kinase DYRK1A and Interferes with Neurite Formation
FEBS J
276
(21)
6324 -
6237
** DOI : 10.1111/j.1742-4658.2009.07346.x**

1998
Ab initio and Density Functional Theory Studies of the Potential Energy Sarface of the Linc→Linc Isomerization
J. Mol. Struct.
428
42 -

Jursic B. S.
1996
Computing Transitoin State Structures whit Density Functionl Theory Methods
Theoretical and Computational Chemsitry
4
709 -
** DOI : 10.1016/S1380-7323(96)80101-4**

Jursic B. S.
1999
Density Functional Theory and Complete Basis Set Ab initio Evaluation of Proton Affinity for Some Selected Chemical System
J. Mol. Struct.
487
194 -

Kurita N.
,
Ivanov P. M.
2000
Correlated Ab initio Molecular Orbital (MP3, MP4) and Density Functional (PW91, MPW91) Studies on the Conformations
J. Mol. Struct.
554
184 -

Ozimiski W.
,
Garnuszek P.
,
Bednarek E.
,
Dobrowol C. S. J.
2007
The Platinum Complexes whit Histamin: Pt(II) (Iodo-Hist) cl2and Pt(IV) (Hist)2cl2
Inorg. Chim. Acta
360
19047 -

Becke A. D.
1993
J. Chem. Phys.
9X
5648 -

Becke A. D.
1988
Phys. Rev. A.
38
3098 -
** DOI : 10.1103/PhysRevA.38.3098**

Perdew J. P.
1986
Phys. Rev. B
33
8822 -
** DOI : 10.1103/PhysRevB.33.8822**

Perdew J. P.
,
Wang Y.
1992
Phys. Rev. B
45
13244 -
** DOI : 10.1103/PhysRevB.45.13244**

Perdew J. P.
,
Zunger A.
19
Phys. Rev. B.
23
5048 -

Vosko S. H.
,
Wilk L.
,
Nusair M.
1980
Can. J. Phys.
58
1200 -
** DOI : 10.1139/p80-159**

Nedelec J. M.
,
Hench L. L.
2000
Journal of Non- Crystaline Solids
277
106 -
** DOI : 10.1016/S0022-3093(00)00306-9**

Citing 'Comparison Chosen DFT Methods with Six Basis Sets on Harmine Molecule
'

@article{ JCGMDC_2013_v57n1_20}
,title={Comparison Chosen DFT Methods with Six Basis Sets on Harmine Molecule}
,volume={1}
, url={http://dx.doi.org/10.5012/jkcs.2013.57.1.20}, DOI={10.5012/jkcs.2013.57.1.20}
, number= {1}
, journal={Journal of the Korean Chemical Society}
, publisher={Korean Chemical Society}
, author={Malekzadeh, Maryam
and
Heshmati, Emran}
, year={2013}
, month={Feb}