Antioxidant activity of some phenolic aldehydes and their diimine derivatives: A DFT study
Accepted Manuscript Antioxidant Activity of Some Phenolic Aldehydes and Their Diimine Derivatives: A DFT Study Romesh Borgohain, Ankur Kanti Guha, San...
Accepted Manuscript Antioxidant Activity of Some Phenolic Aldehydes and Their Diimine Derivatives: A DFT Study Romesh Borgohain, Ankur Kanti Guha, Sanjay Pratihar, Jyotirekha G. Handique PII: DOI: Reference:
Please cite this article as: R. Borgohain, A.K. Guha, S. Pratihar, J.G. Handique, Antioxidant Activity of Some Phenolic Aldehydes and Their Diimine Derivatives: A DFT Study, Computational & Theoretical Chemistry (2015), doi: http://dx.doi.org/10.1016/j.comptc.2015.02.014
This is a PDF file of an unedited manuscript that has been accepted for publication. As a service to our customers we are providing this early version of the manuscript. The manuscript will undergo copyediting, typesetting, and review of the resulting proof before it is published in its final form. Please note that during the production process errors may be discovered which could affect the content, and all legal disclaimers that apply to the journal pertain.
Antioxidant Activity of Some Phenolic Aldehydes and Their Diimine Derivatives: A DFT Study Romesh Borgohaina, Ankur Kanti Guhab, Sanjay Pratiharc and Jyotirekha G. Handiquea* a
Department of Chemistry, Dibrugarh University, Dibrugarh-786004, Assam, India
b
Department of Chemistry, Indian Institute of Technology, Kanpur, Uttar Pradesh, India-
208016 c
Department of Chemical Sciences, Tezpur University, Tezpur, Assam, India-784028
1. INTRODUCTION The antioxidant activity of a compound indicates its ability to inhibit or prevent the oxidation of oxidizable material by scavenging free radicals like ROS (Reactive Oxygen Species), RNS (Reactive Nitrogen Species) and thereby protects the body from oxidative damage caused by the free radicals [1]. There are mainly two pathways through which antioxidants can scavenge the free radicals viz., (1) hydrogen atom transfer (HAT) and (2) single electron transfer (SET) mechanism. Bond dissociation energy (BDE) and ionization potential (IP) are the key factors that determine the mechanism and efficiency of antioxidants [2]. However, two other mechanisms of radical scavenging known as single electron transfer followed by proton transfer (SET-PT) and sequential proton loss electron transfer (SPLET) may also occur [3]. HAT-based mechanism involves the measurement of the classical ability of an antioxidant to scavenge free radicals by hydrogen donation to form stable compounds. This mechanism is more relevant to the radical chain breaking antioxidant capacity.
AH + R RH + A (Where AH = any hydrogen donor and R is a radical. And A is more stable than R) Relative reactivity of HAT mechanism is determined by the bond dissociation energy of the H-donating group in the potential antioxidant. HAT reactions are solvent and pH dependent and are generally quite rapid, typically completed in seconds to minutes. SET-based mechanism deals with the ability of a potential antioxidant to transfer one electron to reduce any compounds, including metals, carbonyls and radicals; M(III) + AH AH + M(II) Deprotonation and ionization potential of the reactive functional group determines the relative reactivity in SET-based mechanism. SET reactions are pH dependent showing a decrease with increasing pH, reflecting increased electron donating capacity with deprotonation. In comparison to HAT, the SET mechanism is strongly solvent dependent due to solvent stabilization of the charged species [4,5]. In SET-PT mechanism, the first step involves the transfer of an electron from the antioxidant to form a radical cation and in the next step, deprotonation from the radical cation occurs. AH AH+ + eAH+ A + H+ Again, in SPLET mechanism loss of proton from the antioxidant occurs followed by transfer of an electron. AH A- + H+ A- A + eThus, compounds having antioxidant property can scavenge free radicals by releasing their phenolic hydrogen atoms following these four mechanisms. It is to be noted that all these mechanisms may or may not co-exist and they are dependent upon the
solvent properties and radical characters. However, the net result of all these mechanisms is same [3]. Phenolic compounds have great importance as antioxidants to inhibit the oxidation of materials having commercial and biological importance [6]. They possess different antioxidant activity potentials because of their phenolic hydroxyl groups which can act as a hydrogen or electron donor [7]. Thus they are used as universal remedy for a wide range of lifestyle-related diseases such as early aging, cancer, diabetes, cardiovascular and other degenerative diseases [8]. With the development of quantum chemistry and computational methodologies, it becomes easier to calculate BDEs and IPs with accuracy equivalent to or greater than those obtained from experiments. Therefore, theoretical calculations could be used as a useful tool for predicting antioxidant activity of compounds as well as to design novel potential antioxidants. Wright and his co-workers [2] studied the antioxidant activity of various phenolic compounds using density functional theory and made a comparison between HAT and SET mechanisms. They concluded that HAT mechanism predominates in most cases [5]. Reis and his co-workers also theoretically evaluated the antioxidant activity of various phenolic compounds using different levels of calculations [9]. Again, Javan et. al. investigated computationally the H-atom vs electron transfer mechanism for antioxidant activity of some bromophenols and compared their observation to those for which the antioxidant activity was previously evaluated [10]. Analysis of antioxidant activity of various gallic acid derivatives as well as the values of some density–based reactivity descriptors was carried out by Kar et. al. and they found that phenolic compounds are more stable in gas phase than in solvent medium [11]. The antioxidant radical scavenging capacity of pyranoanthocyanins was investigated by Russo et. al. [12] They reported that compounds forming radicals, which are stabilized by intramolecular H-
bonds, generally favors HAT mechanism while compounds having a higher degree of electron delocalization favors SET mechanism. Similarly, theoretical calculations were performed by Mohajeri et. al. to predict the antioxidant activity of some compounds like phenolic acids, vitamins etc. and they found that in case of HAT mechanism, antioxidants with lower BDE undergo hydrogen abstraction with low barrier and considerable exothermicity [13]. Alberto et. al studied the antioxidant character of trolox in different solvents using DFT. They considered several reaction mechanisms and different chemical nature of the free radicals and found good agreement of their results with experimentally available data [14]. Again, Leopoldini et. al investigated the effectiveness of caffeic acid in scavenging hydroxyl radical using hybrid density functional theory. They also analyzed different pathways of radical scavenging and concluded that caffeic acid is a very good radical scavenger[15]. Thus, study of the antioxidant capacity of phenolic compounds and their derivatives is continuing and expected to be continued in future. Several theoretical and experimental works have been performed to evaluate the antioxidant activity of the phenolic compounds [16-20]. In our present study, we have chosen three plant-based phenolic aldehydes [21-23] viz., 3-hydroxybenzaldehyde (I), 3,4-dihydroxybenzaldehyde (II) and vanillin (III), having antioxidant property and designed some diimines [24,25] using three ,-diamines of different spacer lengths viz., 1,2-diaminoethane (IV), 1,4-diaminobutane (V) and 1,6diaminohexane (VI) (Scheme 1). We have calculated the O-H bond dissociation energy (BDE) and vertical ionization potentials (IP V) of these derivatives (Scheme 1). Since, antioxidants generally work in the physiological liquids, we therefore, included water alongwith benzene (having dielectric constants of 80 and 2.25 respectively) in our calculations. Compounds are named as shown in Scheme 1. For example, diimines of 3hydroxybenzaldehyde (I) with n = 1, 2 and 3 will be referred to as Ia, Ib and Ic
respectively and the corresponding radical of (I) will be referred to as I-rad. Similarly, the diimine of 3,4-dihydroxybenzaldehyde (II) with n = 1 will be referred to as IIa and its corresponding radical will be referred to as IIa-rad. 2. COMPUTATIONAL DETAILS All the structures were fully optimized at gradient corrected BP86 level of theory [26]. We used both 6-31+G(d) and 6-311++G(d,p) basis set for all the atoms. Frequency calculations were performed at the same level of theory to characterize the nature of the stationary points. All the structures were found to be local minimum with real values of frequencies. Vertical ionization potential IP v, has been calculated by running single point calculations by taking the optimized geometry of the neutral and the cation at the same geometry of the neutral. All the calculations were performed using the Gaussian03 suite of program [27]. We have also calculated single point energy calculation at B3LYP/6311++G(d,p) on the optimized geometries to calculate BDE. Natural bond orbital (NBO) calculations have been performed at BP86/6-31+G(d) level of theory to characterize the electronic structures of these molecules [28]. Hyperconjugation increases the stability of these compounds and hence, stability arising from hyperconjugation has been calculated using the NBO method. Hyperconjugation arising from the overlap between the occupied and unoccupied orbital leads to a stability of the molecules when these orbitals are properly oriented. This interaction can be quantitatively expressed by the second-order perturbation interaction energy (E2) [29]. It can be deduced from the second-order perturbation approach:
where qi is the donor orbital occupancy, εi, εj are diagonal elements (orbital energies), and F(i, j) is the off-diagonal NBO Fock matrix element. Spin densities of the radicals were
calculated at BP86/6-31+G(d) level of theory and plotted to examine the extent of delocalization of the unpaired electron throughout the molecular systems. Solvent effects were computed within the framework of polarized continuum model (PCM) as implemented in Gaussian 03 [30] in BP86/6-31+G(d) level of theory. Owing to the nature of phenolics studied with different spacer lengths (Scheme 1), conformational flexibility is expected. Hence, we also included dispersion correction of Grimme et al using the ωB97XD [31,32] to check the conformational flexibility of some molecules (Ia-c, IIa and IIIb, Scheme 1). 3. RESULTS AND DISCUSSION 3.1 Bond Dissociation Energies (BDE) and Ionization Potentials (IP V): To calculate the BDE and IP V values, minimum energy conformations calculated at BP86/6-31+G(d) are used. Further, dispersion corrected functional such as ωB97XD has been used to reoptimize some of the structures (Ia-c, IIa and IIIa in Scheme 1). It is interesting to note that the minimum energy conformations of these molecules are almost similar. Figure 1 represents the minimum energy conformations calculated at BP86 and ωB97XD level of theory using the 6-31+G(d) basis set. It is evident from Figure 1 that there is very minor change in the dihedral angles in the minimum energy structures after addition of dispersion correction. Thus, it is expected that the geometries of these molecules will not show major conformational changes after addition of dispersion correction and hence, discussion of the rest of the text will be based on BP86 results.
In the HAT mechanism, the bond dissociation energy (BDE) of the O-H bond is an important parameter in evaluating the antioxidant property because the strength of the O-H bonds determines the effectiveness of H abstraction of a particular phenolic compound. This H atom can be transferred to the reactive radical intermediates such as hydroxyl,
alkoxyl, peroxyl, and hydroperoxyl radicals formed during degradation reactions [33]. Table 1 contains BDE values calculated at the BP86/6-31+G(d) level of theory in gas phase and in solvent (water and benzene) medium. It should be noted that the calculated BDE value for phenol (ca. 87.53 kcal/mol in gas phase) deviates from the experimental value by ~1.0 kcal/mol (experimental BDE of phenol is 88.3 ± 0.8) [34]. This might be due to the effect of basis set or functional used. Hence, we also calculated the BDE of phenol at B3LYP [35] using a larger basis set 6-311++G(d,p). The calculated value of BDE at B3LYP/6-311++G(d,p) level of theory is found to be 81.7 kcal/mol which also deviates from the experimental value by ~7.0 kcal/mol. Thus, it appears that BP86/631+G(d) produces better results and hence, in this present work, BP86/6-31+G(d) (in gas phase and in solvents like water, benzene) and BP86/6-311++ G(d,p) (only in gas phase) values are reported owing to the large size of the phenolics considered in this study. Calculations with 6-311++G(d,p) level produces slightly higher BDE values, however, a similar trend with 6-31+G(d) calculated values is maintained. From compound I to II, there is a decrease in BDE with increasing number of –OH group. However, there is again an increase in BDE in case of compound III due to blocking of one –OH group by methyl group. The calculated BDE values (Table 1) are found to be highest in water and lowest in benzene although the differences are not so significant. It is interesting to note that increasing the spacer length via amino substitution decreases the BDE of the phenolics with an exception of (IIIc) where a slight increase in BDE value is observed compared to the parent vanillin. Thus, it is clear that increasing the spacer length results in more favorable H-atom abstraction of the parent phenolic compounds. The calculation of ionization potentials (IP V) of these antioxidants is another important parameter which provides hint towards the possible single electron transfer
(SET) mechanism of these compounds. Lower the value of ionization potential, higher is the probability of a particular compound to undergo single electron transfer resulting in the formation of ArOH•+. The calculated vertical ionization potential (IP V) values of all these phenolics are found to be highest in the gas phase and lowest in water medium. Similar IPV values are obtained with both 6-31+G(d) and 6-311++G(d,p) basis sets. The effect of high solvent polarity on lowering the ionization potential values were previously observed by Belcastro et. Al [36]. This implies that solvent polarity plays a crucial role in dictating the single electron transfer mechanism of an antioxidant. The polarity of the solvent has a significant effect on the IP V values of these compounds unlike BDE values where no such significant variation was observed. Like BDE values, increase in spacer length via amino substitution leads to decrease in IPV values of the phenolics, thereby, favoring single electron transfer.
Figure 2 represents the spin density plots of the hydrogen atom-abstracted radicals. It is evident from Figure 1 that the unpaired electron in all these radicals is delocalized to a great extent to the neighboring aromatic rings. This delocalization stabilizes the radicals and makes the phenolic compounds suitable towards hydrogen release to other potent radicals formed during degradation reactions [7] and thus, these phenolic compounds may act as potent antioxidants. The spin density plots of the cationic radicals formed during single electron transfer are shown in supporting information (Figure S1), which also reveals such delocalization. Moreover, Table 2 collects all the energetics arising from the hyperconjugative interaction between the non bonding orbital of oxygen atom to the neighboring σ* or π* orbitals of the C-C bonds of the aromatic rings. This hyperconjugative interaction adds to the stability of these radicals.
It is evident from Table 2 that all these radicals enjoy the energetic stabilization arising from the hyperconjugative interaction between the non bonding orbital of oxygen atom to the neighboring σ* or π* orbitals of the C-C bonds of the aromatic rings. It is also evident that the interaction energy arising from this hyperconjugation increases with the introduction of amino spacer in the molecules. Hence, with increase in amino spacer length, the stability of the radicals increases thereby, rendering them to show superior antioxidant ability. The calculated dipole moments also increase with the increase in amino spacer lengths, indicating higher charge rearrangement and thus, contributing to higher stability of the radicals. Moreover, the calculated dipole moments of all the radicals increase with increase in dielectric constant of the medium. The calculated values of the dipole moments in high dielectric water solvent (ε = 80.0) are found to be highest for all the radicals indicating that in high dielectric solvents like water, significant reorganization of charge distribution takes place which render higher stability to these radicals. 3.2 Frontier Molecular Orbitals: According to the frontier orbital theory in DFT, the higher the energy of the HOMO (highest occupied molecular orbital), the easier they lose electrons and the reaction becomes faster [37]. Thus, the energy of the HOMO provides a quantitative measure of the electron donating ability of the antioxidants and thereby, providing a measure of these species to undergo SET mechanisms. On the other hand, the lower the energy of the LUMO (lowest unoccupied molecular orbital), the higher will be the ability of the species to accept electrons. In general, a higher HOMO of the phenolic compounds (ArOH) under study and higher LUMO of the radicals (ArO•) makes the forward reaction of Eq. 1 highly favorable and
the phenolic compounds must show antioxidant activity. Table 3 collects the energies of the HOMO of the phenolics and LUMO of the corresponding radicals. R• + ArOH → RH + ArO•
(1)
where R• is any reactive radical intermediate (here DPPH•) formed during degradation reactions [28] and ArOH is the phenolic compounds under consideration.
It is evident from Table 3 that all the phenolic compounds under consideration have higher HOMO than the RH and higher LUMO than R•. This indicates that all these phenolic compounds under consideration can neutralize the free radical according to Eq. 1 and thus, they may act as potent antioxidants. The calculated energies of the HOMO and LUMO of the phenolic compounds and their radicals, respectively, are found to increase with the increase in polarity of the medium, i.e., from gas to water. This is also in accordance with the calculated IPV values (Table 1) which were found to be lowest in water medium. Thus, both the frontier orbital approach and calculated IPV values indicate that these phenolic compounds have higher antioxidant activity in polar solvents. 4. CONCLUSION Quantum chemical calculations have been carried out to investigate the antioxidant capacity of some phenolic compounds and some of their designed diimine derivatives with different spacer lengths. The calculated BDE suggests that with increase in spacer length, the H-atom transfer becomes more prominent. Moreover, the polarity of the medium also lowers the O-H BDE values suggesting their superior H-atom transfer capacity in polar solvents. The calculated vertical ionization potential (IP V) values suggest that with increase in spacer length, single electron transfer from these phenolics becomes easy and thereby making them superior antioxidants. Like BDE, the increase in polarity of the medium decreases the IPV values and makes the single electron transfer more favorable.
Hyperconjugative interactions involving the non-bonding orbital of the phenolic oxygen to the σ* or π* orbitals of the C-C bonds of the aromatic ring are found to stabilize the radicals. The frontier orbital approach reveals that all these phenolics have higher lying HOMO and higher lying LUMO compared to DPPH free radical and thereby, they are expected to show superior antioxidant capacity. The syntheses of these designed diimine derivatives are going on in our laboratory and will be reported in due time. ACKNOWLEDGEMENT The authors are thankful to Department of Science and Technology, Govt. of India for financial assistance (Project No. SR/S1/OC-11/2008, dated 29-10-2008). SUPPORTING INFORMATION Cartesian coordinates of the phenolic compounds in gas phase and spin density plot of the cationic radicals (Figure S1). REFERENCES (1) J. Dai, R. J. Mumper, Plant Phenolics: Extraction, Analysis and Their Antioxidant and Anticancer Properties, Molecules 15(2010) 7313–7352. (2) J. S. Wright, E. R. Johnson, G. A. Di Labio, Predicting the Activity of Phenolic Antioxidants: Theoretical Method, Analysis of Substituent Effects, and Application to Major Families of Antioxidants. J. Am. Chem. Soc. 123 (2001) 1173-1183. (3) Z. Markovic, D. Milenkovic, J. Dorovic, J. M. D. Markovic, V. Stepanic, B. Lucic, D. Amic; PM6 and DFT study of free radical scavenging activity of morin; Food Chemistry 134 (2012) 1754–1760 (4) B. Ou, D. Huang, H. M. Woodill, J. A. Flanagan, E. K. Deemer, Analysis of Antioxidant Activities of Common Vegetables Employing Oxygen Radical Absorbance Capacity (ORAC) and Ferric Reducing Antioxidant Power (FRAP) Assays: A Comparative Study. J. Agric. Food Chem. 50 (2002) 3122-3128.
(5) R. L. Prior, X. Wu, K. Schaich, Identification and Characterization of Anthocyanins by High-Performance Liquid Chromatography−Electrospray Ionization−Tandem Mass Spectrometry in Common Foods in the United States: Vegetables, Nuts, and Grains. J. Agric. Food Chem. 53 (2005) 3101-3113. (6) B. Halliwell, R. Asechbach, J. Loliger, O. I. Aruoma, The characterization of antioxidants. Food Chem. Toxicol. 33(1995) 601-617 (7) B. H. Havsteen, The biochemistry and medicinal significance of the flavonoids. Pharmacol. Ther. 96 (2002) 67–202. (8) H. A. Jung, J. C. Park, H. Y. Chung, J. Kim, J. S. Choi, Antioxidant flavonoids and chlorogenic acid from the leaves of Eriobotrya japonica. Arch Pharm Res 22 (1999) 213– 218 (9) M. Reis, B. Lobato, J. Lameira, A. S. Santos, C. N. Alves, A theoretical study of phenolic compounds with antioxidant properties. Eur. J. Med. Chem. 42 (2007) 440-446. (10) A. J. Javan, M. J. Javan, J. A. Tehrani , Theoretical Investigation on Antioxidant Activity of Bromophenols from the Marine Red Alga Rhodomela confervoides: H-Atom vs Electron Transfer Mechanism. J. Agric. Food Chem. 61(2013) 1534−1541. (11) D. Kalita, R. Kar, J. G. Handique, A Theoretical Study On The Antioxidant Property Of Gallic Acid And Its Derivatives. J. Theor. Comput. Chem. 11 (2012) 391-402. (12) M. Leopoldini, F. Rondinelli, N. Russo, A. Toscano, Pyranoanthocyanins: A Theoretical Investigation on Their Antioxidant Activity. J. Agric. Food Chem. 58 (2010) 8862–8871. (13) A. Mohajeri, S. S. Asemani, Theoretical investigation on antioxidant activity of vitamins and phenolic acids for designing a novel antioxidant. J. Mol. Struct. 930 (2009) 15–20.
(14) M. E. Alberto, N. Russo, A. Grand, A. Galano, A physicochemical examination of the free radical scavenging activity of Trolox: mechanism, kinetics andinfluence of the environment. Phys. Chem. Chem. Phys., 15(2013) 4642-4650 (15) M. Leopoldini, S. G. Chiodo, N. Russo, M. Toscano, Detailed Investigation of the OH Radical Quenching by Natural Antioxidant Caffeic Acid Studied by Quantum Mechanical Models. J. Chem. Theory Comput. 7(2011) 4218-4233 (16) C. Iuga, J. R. A. Idaboy, N. Russo, Antioxidant Activity of trans-Resveratrol toward Hydroxyl and Hydroperoxyl Radicals: A Quantum Chemical and Computational Kinetics Study. J. Org. Chem. 77(2012) 3868-3877 (17) M. Leopoldini, N. Russo, M. Toscano, The molecular basis of working mechanism of natural polyphenolic antioxidants. Food Chemistry 125 (2011) 288-306 (18) (a) A. N. Queiroz, A. Q. Bruno, G. W. M. Jr. Moraes, R. S. Borges, A theoretical antioxidant pharmacophore for resveratrol. Eur. J. Med. Chem. 44 (2009) 1644–1649. (b) E. Anouar, P. Kosˇinová, D. Kozlowski, R. Mokrini, J. L. Duroux, P. Trouillas, New aspects of the antioxidant properties of phenolic acids: A combined theoretical and experimental approach. Phys. Chem. Chem. Phys. 11(2009) 7659–7668 (c) S. Dudonné, X. Vitrac, P. Coutiére, M. Woillez, J. –M. Mérillon, Comparative Study of Antioxidant Properties and Total Phenolic Content of 30 Plant Extracts of Industrial Interest Using DPPH, ABTS, FRAP, SOD, and ORAC Assays. J. Agric. Food Chem. 57(2009) 1768– 1774 (d) B. Rojano, J. Saez, G. Schinella, J. Quijano, E. Vélez, A. Gil, R. Notario, Experimental and theoretical determination of the antioxidant properties of isoespintanol (2-isopropyl-3,6-dimethoxy-5-methylphenol). J. Mol. Struct. 877 (2008) 1–6. (19) (a) P. Trouillas, P. Marsal, A. Svobodová, J. Vostálová, R. Gazˇák, J. Hrbá, P. Sedmera, V. K. R. Lazzaroni, J. –L. Duroux, D. Walterova, Mechanism of the Antioxidant Action of Silybin and 2,3-Dehydrosilybin Flavonolignans: A Joint Experimental and
Theoretical Study. J. Phys. Chem. A 112(2008) 1054-1063 (b) I. Tejero, N. G. García, A. G. Lafont, J. M. Lluch, Tunneling in Green Tea: Understanding the Antioxidant Activity of Catechol-Containing Compounds. A Variational Transition-State Theory Study. J. Am. Chem. Soc. 129 (2007) 5846-5854. (20) (a) E. G. Bakalbassis, A. T. Lithoxoidou, A. P. Vafiadis, Theoretical Calculation of Accurate Absolute and Relative Gas- and Liquid-Phase O-H Bond Dissociation Enthalpies of 2-Mono- and 2,6-Disubstituted Phenols, Using DFT/B3LYP. J. Phys. Chem. A 107(2003) 8594-8606 (b) M. Leopoldini, I. P. Pitarch, N. Russo, M. Toscano, Structure,
Conformation, and Electronic Properties of Apigenin, Luteolin, and Taxifolin Antioxidants. A First Principle Theoretical Study. J. Phys. Chem. A 108(2004), 92-96 (c) J. S. Wright, Predicting the antioxidant activity of curcumin and curcuminoids; J. Mol. Struct. (THEOCHEM) 591(2002) 207–217 (d) Z. Cheng, J. Ren, Y. Li, W. Chang, Z. Chen, Study on the Multiple Mechanisms Underlying the Reaction Between Hydroxyl Radical and Phenolic Compounds by Qualitative Structure and Activity Relationship. Bioor. Med. Chem. 10(2002) 4067–4073 (e) H. –G. Korth, M. I. Heer, P. Mulder, A DFT Study on Intramolecular Hydrogen Bonding in 2-Substituted Phenols: Conformations, Enthalpies, and Correlation with Solute Parameters. J. Phys. Chem. A 106(2002) 87798789. (21) T. Herraiz, J. Galisteo, C. Chamorro, L-Tryptophan Reacts with Naturally Occurring and Food-Occurring Phenolic Aldehydes To Give Phenolic Tetrahydro-α-carboline Alkaloids: Activity as Antioxidants and Free Radical Scavengers. J. Agric. Food Chem. 51(2003) 2168-2173. (22) D. J. Fitzgerald, M. Stratford, M. J. Gasson, A. Narbad, Structure-Function Analysis of the Vanillin Molecule and Its Antifungal Properties. J. Agric. Food Chem. 53(2005) 1769-1775.
(23) M. Watanabe, Y. Ohshita, T. Tsushida, Antioxidant Compounds from Buckwheat (Fagopyrum esculentum Möench) Hulls. J. Agric. Food Chem. 45(1997) 1039-1044. (24) C. Y. Lee, A. Sharma, J. E. Cheong, J. L. Nelson, Synthesis and antioxidant properties of dendritic polyphenols. Bioorg. Med. Chem. Lett. 19 (2009) 6326–6330. (25) C. Y. Lee, A. Sharma, R. L. Uzarski, J. E. Cheong, H. Xu, R. A. Held, S. K. Upadhaya, J. L. Nelson, Potent antioxidant dendrimers lacking pro-oxidant activity. Free Radical Biol. Med. 50 (2011) 918–925. (26) (a) A. D. Becke, Density-functional exchange-energy approximation with correct asymptotic behavior. Phys. Rev. A 38 (1988) 3098-3100. (b) J. P. Perdew, Density-functional approximation for the correlation energy of the inhomogeneous electron gas. Phys. Rev. B 33 (1986) 8822-8824. (27) M. J. Frisch, G. W. Trucks, H. B. Schlegel, G. E. Scuseria, M. A. Robb, J. R. Cheeseman, J. A. Montgomery, T. Vreven, K. N. Kudin, J. C. Burant, J. M. Millam, S. S. Iyengar, J. Tomasi, V. Barone, B. Mennucci, M. Cossi, G. Scalmani, N. Rega, G. A. Petersson, H. Nakatsuji, M. Hada, M. Ehara, K. Toyota, R. Fukuda, J. Hasegawa, M. Ishida, T. Nakajima, Y. Honda, O. Kitao, H. Nakai, M. Klene, X. Li, J. E. Knox, H. P. Hratchian, J. B. Cross, V. Bakken, C. Adamo, J. Jaramillo, R. Gomperts, R. E. Stratmann, O. Yazyev, A. J. Austin, R. Cammi, C. Pomelli, J. W. Ochterski, P. Y. Ayala, K. Morokuma, G. A. Voth, P. J. Salvador, J. Dannenberg, V. G. Zakrzewski, S. Dapprich, A. D. Daniels, M. C. Strain, O. Farkas, D. K. Malick, A. D. Rabuck, K. Raghavachari, J. B. Foresman, J. V. Ortiz, Q. Cui, A. G. Baboul, S. Clifford, J. Cioslowski, B. B. Stefanov, G. Liu, A. Liashenko, P. Piskorz, I. Komaromi, R. L. Martin, D. J. Fox, T. Keith, M. A. AlLaham, C. Y. Peng, A. Nanayakkara, M. Challacombe, P. M. W. Gill, B. Johnson, W. Chen, M. W. Wong, C. Gonzalez, J. A. Pople, Gaussian 03, revision D.01, Gaussian, Inc.: Pittsburgh, PA, 2004.
(28) (a) E. D.Glendening, A. E.Reed, J. E.Carpenter, F.Weinhold, NBO Program 3.1: University of Wisconsin: W. T.Madison,1988. (b) A. E.Reed, F.Weinhold, L. A.Curtiss, Intermolecular interactions from a natural bond orbital, donor-acceptor viewpoint. Chem. Rev. 88 (1988) 899–926. (29) (a) A. E.Reed, R. B.Weinstock, F. Weinhold, Natural Atomic Orbitals and Natural Population Analysis. J. Chem. Phys. 83(1985) 735−747 (b) A. E. Reed, F. Weinhold, Natural localized molecular orbitals. J. Chem. Phys. 83(1985) 1736−1741 (c) A. E. Reed, L. A. Curtiss, F. Weinhold, Intermolecular interactions from a natural bond orbital, donor−acceptor viewpoint. Chem. Rev. 88 (1988) 899−926 (d) A. E. Reed, P. R. Schleyer, Chemical bonding in hypervalent molecules. The dominance of ionic bonding and negative hyperconjugation over d-orbital participation. J. Am. Chem. Soc. 112(1990) 1434−1445. (30) J. Tomasi, M. Persico, Molecular Interactions in Solution: An Overview of Methods Based on Continuous Distributions of the Solvent. Chem. Rev. 94(1994) 2027–2094. (31) S. Grimme, Semiempirical GGA-type density functional constructed with a longrange dispersion correction, J. Comp. Chem 27(2006) 1787-1799. (32) S. Grimme, Semiempirical hybrid density functional with perturbative second-order correlation, J. Chem. Phys. 124(2006) 034108. (33) Q. Zhu, X. M. Zhang, A. Fry, Bond dissociation energies of antioxidants. Polym. Degrad. Stab. 57 (1997) 43−50. (34) M. Lucarini, G. F. Pudulli, M. Cipollone, Bond Dissociation Enthalpy of alphaTocopherol and Other Phenolic Antioxidants. J. Org. Chem. 59 (1994) 5063-5070. (35) (a) A. D. Becke, Density-functional thermochemistry. III. The role of exact exchange. J. Chem. Phys. 98(1993) 5648–5652. (b) C. Lee, W. Yang, R. G. Parr, Development of
the Colle-Salvetti correlation-energy formula into a functional of the electron density. Phys. Rev. B 37(1988) 785–789. (36) M. Belcastro, T. Marino, N. Ruusso, M. Toscano, Structural and electronic characterization of antioxidants from marine organisms Theor. Chem. Acc. 115(2006) 361-369. (37) J. F. Janak, Proof that ∂E∂ni=ε in density-functional theory. Phys Rev B 18 (1978) 7165–7168.
Captions of figures
Scheme 1. Structures of the phenolic compounds considered in the study Figure 1. Selected dihedral angles (in degrees) of the minimum energy structures of Ia-c, IIa and IIIa at BP86/6-31+G(d) and ωB97XD/6-31+G(d) level of theory. The values within parenthesis refers to ωB97XD/6-31+G(d) level of theory. Figure 2. Plot of spin density of hydrogen atom-abstracted radicals with contour value of 0.05
Figures
Scheme 1. Structures of the phenolic compounds considered in the study
Figure 1. Selected dihedral angles (in degrees) of the minimum energy structures of Ia-c, IIa and IIIa at BP86/6-31+G(d) and ωB97XD/6-31+G(d) level of theory. The values within parenthesis refers to ωB97XD/6-31+G(d) level of theory.
Figure 2. Plot of spin density of hydrogen atom-abstracted radicals with contour value of 0.05.
Table 1. BP86/6-31+G(d) and BP86/6-311++G(d,p) calculated bond dissociation energies (BDE in kcal/mol) and vertical ionization potential (IP V, in kcal/mol) at BP86/6-31+G(d) of the studied phenolics. Compounds
BDE Water Benzene
Gas 6-31 6-311 +G(d) ++G(d,p) I
97.6
102.4
100.1
99.2
Gas 6-31 6-311 +G(d) ++ G(d,p) 276.7 276.9
II
83.3
91.7
89.5
85.5
265.1
III
93.3
97.7
98.3
96.4
Ia
88.8
93.3
91.4
Ib
92.3
96.4
Ic
95.2
IIa
IPV Water Benzene
228.2
272.1
265.4
230.6
260.5
191.4
192.7
145.2
186.7
89.5
182.1
187.6
136.0
177.5
94.6
92.9
159.1
161.3
126.8
147.5
177.2
113.9
97.2
242.1
248.9
147.5
235.2
75.5
84.9
81.2
77.3
186.7
181.5
80.7
182.1
IIb
85.3
89.4
86.5
86.1
184.4
187.2
96.8
179.8
IIc
79.6
90.3
102.7
90.2
223.6
225.9
117.6
212.1
IIIa
77.5
81.7
82.8
79.3
175.2
177.7
131.4
166.0
IIIb
85.4
89.2
89.5
86.7
172.9
180.3
140.6
168.3
IIIc
106.5
144.9
105.5
104.8
131.4
136.0
96.8
124.5
Table 2. Interaction energies (kcal/mol) arising from the hyperconjugative interaction between the non bonding orbital of oxygen atom to the σ* or π* orbitals of the C-C bonds of the aromatic rings in the radicals in gas phase and the calculated dipole moments (in Debye) in different mediums. Radicals
E2 (kcal/mol)
μ Gas
Benzene
Water
I-rad
7.75
1.63
1.99
2.55
II-rad
8.69
1.20
1.51
2.13
III-rad
9.35
2.32
2.61
3.23
Ia-rad
14.91
3.88
4.79
6.51
Ib-rad
14.00
4.42
5.37
7.33
Ic-rad
14.32
7.22
7.89
9.36
IIa-rad
14.45
3.83
4.41
5.31
IIb-rad
14.22
5.96
9.32
11.72
IIc-rad
13.85
7.36
8.03
10.68
IIIa-rad
16.07
5.25
6.44
8.63
IIIb-rad
13.96
7.40
8.75
11.18
IIIc-rad
8.82
7.12
8.32
8.89
Table 3. Energies of HOMO and LUMO of the phenolics under study and their radicals in kcal/mol. Phenolics
HOMO
Radicals
Gas
Benzene
Water
I
-148.1
-147.3
-146.3
II
-146.5
-144.1
III
-146.7
Ia
LUMO Gas
Benzene
Water
I-rad
-69.6
-69.5
-69.1
-143.4
II-rad
-97.0
-96.8
-96.8
-145.8
-145.1
III-rad
-62.1
-61.7
-60.8
-129.9
-129.2
-128.5
Ia-rad
-53.3
-52.8
-52.3
Ib
-123.0
-122.3
-121.7
Ib-rad
-53.4
-52.6
-51.9
Ic
-114.9
-114.3
-113.5
Ic-rad
-62.1
-61.9
-61.2
IIa
-121.3
-120.9
-120.2
IIa-rad
-51.6
-51.1
-50.6
IIb
-126.5
-125.8
-125.1
IIb-rad
-51.5
-51.2
-50.2
IIc
-106.2
-105.7
-105.2
IIc-rad
-62.6
-62.1
-61.5
IIIa
-115.7
-115.3
-114.7
IIIa-rad
-48.7
-48.1
-47.6
IIIb
-126.1
-126.1
-125.4
IIIb-rad
-48.5
-48.3
-47.9
IIIc
-100.1
-99.6
-98.3
IIIc-rad
-60.6
-60.1
-59.8
*RH in gas: HOMO = -176.3; in benzene: HOMO = -175.4; in water: HOMO = -175.1 R• in gas: LUMO = -93.9; in benzene: LUMO = -93.2; in water: LUMO = -91.7
Graphical abstract
Highlights : Both HAT and SET mechanism for antioxidant property of some diimine derivatives of phenolic aldehydes were studied using DFT. Increasing spacer length of the diimine derivatives makes them good antioxidants. Increase in polarity of the solvent mediums increases the antioxidant property of the compounds.