Next Article in Journal
Low Levels of IgM and IgA Recognizing Acetylated C1-Inhibitor Peptides Are Associated with Systemic Lupus Erythematosus in Taiwanese Women
Previous Article in Journal
Photocatalytic Atom Transfer Radical Addition to Olefins Utilizing Novel Photocatalysts
 
 
Comment published on 21 November 2022, see Molecules 2022, 27(22), 8092.
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Comparative Study of Antioxidant Potential of Selected Dietary Vitamins; Computational Insights

by
Dinesh R. Pandithavidana
* and
Samith B. Jayawardana
Department of Chemistry, Faculty of Science, University of Kelaniya, Kelaniya 11600, Sri Lanka
*
Author to whom correspondence should be addressed.
Molecules 2019, 24(9), 1646; https://doi.org/10.3390/molecules24091646
Submission received: 29 March 2019 / Revised: 21 April 2019 / Accepted: 22 April 2019 / Published: 26 April 2019
(This article belongs to the Section Computational and Theoretical Chemistry)

Abstract

:
Density functional theory (DFT) was used to explore the antioxidant properties of some naturally occurring dietary vitamins, and the reaction enthalpies related to various mechanisms of primary antioxidant action, i.e., hydrogen atom transfer, single electron transfer–proton transfer, and sequential proton loss–electron transfer were discussed in detail. B3LYP, M05-2X, and M06-2X functionals were utilized in this work. For aqueous phase studies, the integral equation formalism polarized continuum model (IEF–PCM) was employed. From the outcomes, hydrogen atom transfer (HAT) was the most probable mechanism for the antioxidant action of this class of compounds. Comparison of found results with experimental data (available in literature), vitamin C possesses the lowest enthalpy values for both proton affinity (PA) and bond dissociation energy (BDE)in the aqueous phase, suggesting it as the most promising candidate as an antioxidant. Accordingly, these computational insights encourage the design of structurally novel, simple vitamins which will be more economical and beneficial in the pharmaceutical industry.

1. Introduction

The origin of the harmful process called oxidative stress lies in the extreme production of free radicals such as reactive oxygen species (ROS),reactive nitrogen species (RNS), and reactive sulfur species (RSS) with half-lives of only a few nanoseconds, whose effects can seriously alter cell structures (e.g., membranes) and damage biomolecules such as lipids, lipoproteins, proteins, and nucleic acids [1,2,3]. Among various types of antioxidants, natural vitamins represent a wide group of chemically distinct, water-soluble, and biologically active compounds which serve to inhibit or delay the oxidation of important macromolecules of cells by scavenging those free radicals [4,5].These antioxidants may also have potential to protect against a number of disease conditions such as aging, atherosclerosis, cancer, autoimmune conditions, asthma, and arthritis. Some vitamins (vitamin C and vitamin B6) have been often been reported as antioxidants which have the capability to limit the oxidative damage in humans and lower the risk of certain chronic diseases [6,7,8].
It is broadly recognized that the most vital structural characteristic which facilitates effective antioxidant activity is the presence of one or more conjugated OH groups or COOH groups, which boosts the ability of such a molecule to quench the free radicals [4,5]. There are three major proposed mechanisms which can be used to clarify how antioxidants release the atomic hydrogen from their OH group or COOH group and scavenge free radicals; (1) the hydrogen atom transfer (HAT) mechanism,(2) the single electron transfer–proton transfer (SET–PT) mechanism, and (3) the sequential proton loss–electron transfer (SPLET) mechanism [6,9].Though bond dissociation energy (BDE) is used as a key factor to determine antioxidant capacity of the HAT mechanism, adiabatic ionization potential (IP) and proton dissociation enthalpy (PDE) are used to determine antioxidant efficiency of the SET–PT mechanism, while proton affinity (PA) and electron transfer enthalpy (ETE) are used as computational parameters to the investigate the antioxidant efficiency of the SPLET mechanism [10,11,12] as illustrated in Scheme 1. It should be apparent that all these mechanisms may or may not co-exist, and they are conditional upon the radical characters and solvent properties. Despite this, the net upshot of all these mechanisms is almost the same.
Although some computational studies have been carried out to investigate the antioxidant activity of individual vitamins and their derivatives [6,7], a comparative density functional theory (DFT) study of a group of dietary vitamins (in Figure 1) was not computed. In this investigation, the radical scavenging activity of vitamin A (retinol), vitamin B1 (thiamine), vitamin B3 (nicotinic-d4 acid), vitamin B6 (pyridoxine), and vitamin C (l-ascorbic acid) have been interpreted in terms of some key thermochemical parameters such as BDE, IP, PDE, PA, and ETE. Revelation of the relationship between structure and reactivity is another goal of this study. Efforts have been made to find out such correlation from these thermochemical parameters. Thus, the current investigation may be helpful in shed more light on this front.

2. Results and Discussion

To validate the computational procedure, a density functional theory (DFT) method with B3LYP [13], M05-2X [14], and M06-2X [15] functionals were employed to compute enthalpies of BDE, IP, PDE, PA, and ETE of vitamin A, vitaminB1, vitaminB3, vitaminB6, and vitamin C. When M05-2X and M06-2X functionals were used, the BDE, IP, PDE, PA, and ETE data were very close in comparison to the B3LYP functional as illustrated in Table 1 (Supplementary materials are also available online). Although thermodynamic properties computed from the B3LYP functional had deviated slightly, it was found that “growing orders” of these properties were identical for all three functionals. Thus, mainly results of M05-2X and M06-2Xfunctionalshave been discussed in the further text.
Since these vitamin molecules exhibit their antioxidant potential in an aqueous cellular environment in reality [5,16], optimized structures of vitamins in aqueous phase (water) have been used for the computational investigations as they are shown in Figure 2.
The comparison of bond dissociation enthalpies (BDE) of these dietary vitamins in aqueous phase is shown in Figure 3. The -OH bond dissociation enthalpies grew in this order: Vitamin C < vitamin A < vitamin B6 < vitamin B1 < vitamin-B3.
As BDE is inversely proportional to the antioxidant activity of the compound, having the lowest value for BDE, vitamin C exhibited the highest antioxidant potential through the HAT mechanism. The presence of a conjugated vinylic hydroxyl group in vitamin C favors the formation of an alkoxy radical (RO) which is stabilized through the resonance effect.
It could be assumed that the electron-withdrawing effect and dimer formation effect (it was computed that formation of strong intermolecular hydrogen bonding is thermodynamically more favorable) of the carboxylic (-COOH) group in vitaminB3 may have caused the highest BDE value (and showing the lowest antioxidant potential) when compared to other vitamins. Due to the presence of allylic or vinylic -OH groups, vitamin A and vitaminB6 demonstrated some moderate effect.
One of the numerical parameters associated with the SET-PT mechanism is the ionization potential (IP). Data in Figure 4 indicate how the IPs varies between these dietary vitamins in aqueous phase. It was revealed that vitamin A should be the best electron donor among the studied vitamins. Due to possessing the lowest IP value, vitamin A illustrates the best antioxidant activity according to the first step of the SET-PT mechanism.
Aqueous phase ionization potentials increased in the following order: Vitamin A < vitamin C < vitamin B6 < vitamin B3 < vitamin B1.
Again, the highest IPs were computed for vitaminB1and vitaminB3 due to the existence of a positive charge on vitaminB1and the presence of an electron-withdrawing COOH groupinvitaminB3. The stepwise mechanism of SET–PT is allied to the predisposition of antioxidants in giving electrons. In the case of this mechanism, an electron is transferred from the antioxidant to the free radical leading to the formation of the radical cation ROH+● which consequently deprotonates. Therefore, adiabatic ionization potential (IP) and proton dissociation enthalpy (PDE) are the most significant thermodynamic properties for analyzing the feasibility of this mechanism. The computed proton dissociation enthalpies (PDEs) corresponding to the second step of SET–PT mechanism followed a somewhat different order than the IPs (as shown in Figure 5): Vitamin B1< vitamin C< vitamin B6 < vitamin B3 < vitamin A.
The amount of energy required to form the radical cation is significantly lower than that needed to accomplish the second step of the SET-PT mechanism [11,12]. Accordingly, the deprotonation of the radical cation ROH+● is the most apparent step that limits the reaction rate of such a mechanism in aqueous medium (polar solvents).In general, IPs and PDEs often demonstrate nearly “mirror trends” as the lower the IP, the higher the PDE [17,18].Thus, computed data reveal that the IPs and PDEs of vitaminB1 and vitamin A are mirror trends of each other. Vitamins with low adiabatic ionization potential (IP) values are more susceptible to ionization and have a stronger antioxidant potential. Anyhow, the energies required to accomplish the whole SET–PT mechanism propose that vitaminB1 and vitaminB3 have relatively low antioxidant activity.
When the “combine effect” was considered (as the SET–PT mechanism = IP + PDE), vitamin C exhibited the most favorable antioxidant activity, and vitaminB3was thermodynamically unfavorable to follow the SET–PT mechanism (Table 1).
Deprotonation of the -OH group or -COOH group represents the first step of the SPLET mechanism. The acidity, the nature of hydrogen bonding, and the presence of electron-withdrawing groups are key influences which determine the enthalpy of proton affinity (PA). According to computed results, proton affinities of these selected vitamins shown in Figure 6 increase in this order: Vitamin C < vitamin B3 < vitamin B6 < vitamin B1 < vitamin A.
It could be assumed that the presence of the vinylic conjugated -OH group in vitamin C and the acidic -COOH group in vitaminB3 lowered proton affinities (PA) of these two vitamins in comparison to the rest. The development of the alkoxide anion in the first step of the SPLET mechanism (PA values of Figure 6) requires much more energy to occur than that required in the electron transfer from the alkoxide anion to the free radical (ETEs of Figure 7), showing that the first step must be the slowest one for this SPLET mechanism
To complete the second step of SPLET mechanism, electron transfer should take place. Electron transfer enthalpies (ETEs) of the vitamins were ascending in this order: Vitamin A < vitamin B6 < vitamin B1 < vitamin C < vitamin B3.
In the case of vitamin A, its higher PA is accompanied with considerably lower ETE in comparison to the rest of the studied vitamins. Hence, the deprotonated form of vitamin A should be the finest electron donor from a thermodynamic point of view. Due to the presence of electron-withdrawing carboxylate anions, vitaminB3 exhibited the highest ETE in aqueous phase. Since the SPLET mechanism is the net effect of PA and ETE, vitamin C was again confirmed to possess the highest antioxidant potential when compared to the other dietary vitamins studied here.

3. Computational Details

All DFT calculations have been performed using the Gaussian 09 program package (University of Colombo, Department of Chemistry, Colombo-03, Sri Lanka). The geometry of all vitamins, including their radicals, radical cations, and anions, have been fully optimized by employing the hybrid functional of B3LYP, M05-2X, and M06-2X with the basis set of 6-311++G(d, p) (without any constraints) [19]. Optimized structures were confirmed to be real minima by frequency analyses. Aqueous phase data were obtained for geometries of studied optimized species, employing integral equation formalism polarized continuum model (IEF–PCM) [20]. The following thermodynamic properties were computed at T = 298.15 K with non-scaled zero-point energies (ZPE).BDE, IP, PDE, PA, and ETE have been computed in water at 298.15 K following the formulae already applied on similar molecules [20,21], according to the following expressions:
BDE =H(RO) + H(H) − H(ROH)
IP = H(ROH●+) + H(e) −H(ROH)
PDE= H (RO) + H(H+) −H(ROH●+)
PA = H(RO) + H(H+) −H(ROH)
ETE = H(RO) + H(e) −H(RO)
H(ROH), H(RO), and H(ROH●+) represent the protonated form, the radical form, and the radical cation form of the vitamin molecule, respectively. The calculated H(H+) and H(e) enthalpies were obtained as 6.197 kJmol−1 and 3.145 kJmol−1, respectively from the literature [21] available.

4. Conclusions

In this computational investigation, antioxidant potentials of the five natural dietary vitamins were explored from the thermodynamic point of view. HAT, SET–PT, and SPLET mechanisms were used to illustrate their radical scavenging activities in the aqueous phase. Employed functionals provided aqueous phase reaction enthalpies in very good mutual agreement. However, results demonstrated that the enthalpies from the three considered mechanisms do have some deviation, which may provide a primary clue to which mechanism is favored over the other.
In water (aqueous phase), deprotonation of the OH or COOH group represented the preferred process. Employing the IEF–PCM model as an implicit solvent model, vitaminB3 and vitaminB1 displayed lower radical scavenging/antioxidant capabilities according to all three mechanisms employed. This was more pronounced in the mechanisms of HAT and SPLET, having relatively low energy pathways when compared to the SET–PT mechanism. However, in some cases, the HAT mechanism (described by BDE) was more applicable than both SET–PT and SPLET mechanisms. This is because smaller values of BDE were more pronounced in the HAT mechanism when compared to relatively higher IP and PA values. Thus, among the selected dietary vitamin molecules, vitamin C has the lowest values for both PA and BDE in the aqueous phase. In general, results of computational investigation positively agreed with the available experimental work, suggesting that vitamin C possesses the highest antioxidant effect in the aqueous environment. It can be concluded that modification of the chemical structure of vitamin C provides insight into the design of structurally novel, simple vitamins which will be more economical and beneficial in pharmaceutical industry.

Supplementary Materials

The following is available online at https://www.mdpi.com/1420-3049/24/9/1646/s1, Table S1: Reaction enthalpies in kJ mol−1 for the selected dietary vitamins, employing functional of B3LYP, M05-2X and M06-2X in aqueous medium (water).

Author Contributions

Conceptualization, D.R.P.; methodology, S.B.J.; software, S.B.J.; validation, D.R.P. and S.B.J.; formal analysis, S.B.J.; investigation, D.R.P.; resources, D.R.P.; data curation, D.R.P and S.B.J.; writing—original draft preparation, D.R.P.; writing—review and editing, D.R.P.; visualization, D.R.P.; supervision, D.R.P.; project administration, D.R.P.

Funding

This research received no external funding.

Acknowledgments

Authors gratefully acknowledge the support from University of Colombo, Colombo-03, Sri Lanka for providing opportunity to do the DFT calculations (using Gaussian 09 program package) at the Computer Laboratory of the Department of Chemistry, University of Colombo.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Nimse, S.B.; Pal, D. Free radicals, natural antioxidants, and their reaction mechanisms. RSC Adv. 2015, 5, 27986–28006. [Google Scholar] [CrossRef] [Green Version]
  2. Lobo, V.; Patil, A.; Phatak, A.; Chandra, N. Free radicals, antioxidants and functional foods: Impact on human health. Pharm. Rev. 2010, 4, 118–126. [Google Scholar] [Green Version]
  3. Dizdaroglu, M. Oxidative damage to DNA in mammalian chromatin. Mutat Res. DNAging. 1992, 275, 331–342. [Google Scholar] [CrossRef]
  4. Mohajeri, A.; Asemani, S.S. Theoretical investigation on antioxidant activity of vitamins and phenolic acids for designing a novel antioxidant. J. Mol. Struct. 2009, 930, 15–20. [Google Scholar] [CrossRef]
  5. Brewer, M.S. Natural Antioxidants: Sources, Compounds, Mechanisms of Action, and Potential Applications. Compr. Rev. Food Sci. Food Saf. 2011, 10, 221–247. [Google Scholar] [CrossRef] [Green Version]
  6. Skorna, P.; Rimarcík, J.; Poliak, P.; Lukes, V.; Klein, E. Thermodynamic study of vitamin B6 antioxidant potential. Comp. Theor. Chem. 2016, 1077, 32–38. [Google Scholar] [CrossRef]
  7. Bendich, A.; Machlin, L.J.; Scandurra, O.; Burton, G.W.; Wayner, D.D.M. The antioxidant role of vitamin C. Adv. Free Radical Biol. Med. 1986, 2, 419–444. [Google Scholar]
  8. McCall, M.R.; Frei, B. Can antioxidant vitamins materially reduce oxidative damage in humans? Free Radic. Biol. Med. 1999, 26, 1034–1053. [Google Scholar] [CrossRef]
  9. Mazzone, G.; Russo, N.; Toscano, M. Antioxidant properties comparative study of natural hydroxycinnamic acids and structurally modified derivatives: Computational insights. Comp. Theor. Chem. 2016, 1077, 39–47. [Google Scholar] [CrossRef]
  10. Borgohain, R.; Guha, A.K.; Pratihar, S.; Handique, J.G. Antioxidant activity of some phenolic aldehydes and their diamine derivatives: A DFT study. Comp. Theor. Chem. 2015, 1060, 17–23. [Google Scholar] [CrossRef]
  11. Szeląg, M.; Mikulski, D.; Molski, M. Quantum-chemical investigation of the structure and the antioxidant properties of α-lipoic acid and its metabolites. J. Mol. Model. 2012, 18, 2907–2916. [Google Scholar] [CrossRef]
  12. Litwinienko, G.; Ingold, K.U. Abnormal solvent effects on hydrogen atom abstraction.2.Resolution of the curcumin antioxidant controversy. The role of sequential proton loss electron transfer. J. Org. Chem. 2004, 69, 5888–5896. [Google Scholar] [CrossRef] [PubMed]
  13. Becke, A. Density-functional thermochemistry. III. The role of exact exchange. J. Chem. Phys. 1993, 98, 5648–5652. [Google Scholar] [CrossRef]
  14. Zhao, Y.; Schultz, N.E.; Truhlar, D.G. Design of density functionals by combining the method of constraint satisfaction with parametrization for thermochemistry, thermochemical kinetics, and noncovalent interactions. J. Chem. Theory Comput. 2006, 2, 364–382. [Google Scholar] [CrossRef] [PubMed]
  15. Zhao, Y.; Truhlar, D.G. The M06 suite of density functionals for main group thermochemistry, thermochemical kinetics, noncovalent interactions, excited states, and transition elements: Two new functionals and systematic testing off our M06-class functionals and 12 other functionals. Theor. Chem. Acc. 2008, 120, 215–241. [Google Scholar]
  16. Kamat, J.P.; Devasagayam, T.P.A. Nicotinamide (vitamin B3) as an effective antioxidant against oxidative damage in rat brain mitochondria. Redox Rep. 1999, 4, 179–184. [Google Scholar] [CrossRef] [PubMed]
  17. Klein, E.; Lukeš, V. DFT/B3LYP study of the substituent effect on the reaction enthalpies of the individual steps of single electron transfer–proton transfer and sequential proton loss electron transfer mechanisms of phenols antioxidant action. J. Phys. Chem. A 2006, 110, 12312–12320. [Google Scholar] [CrossRef]
  18. Klein, E.; Lukeš, V.; Ilcin, M. DFT/B3LYP study of tocopherols and chromans antioxidant action energetics. Chem. Phys. 2007, 336, 51–57. [Google Scholar] [CrossRef]
  19. Becke, A.D. Density-functional exchange-energy approximation with correct asymptotic behavior. Phys. Rev. A 1988, 38, 3098–3100. [Google Scholar] [CrossRef]
  20. Mazzone, G.; Malaj, N.; Russo, N.; Toscano, M. Density functional study of the antioxidant activity of some recently synthesized resveratrol analogues. Food Chem. 2013, 141, 2017–2024. [Google Scholar] [CrossRef]
  21. Mazzone, G.; Malaj, N.; Galano, A.; Russo, N.; Toscano, M. Antioxidant properties of several coumarin–chalcone hybrids from theoretical insights. RSC Adv. 2015, 5, 565–575. [Google Scholar] [CrossRef]
Sample Availability: Samples of the compounds (optimized all dietary vitamins used here) are available from the authors.
Scheme 1. Schematic representation of studied reaction mechanisms and denotations used for corresponding reaction enthalpies.
Scheme 1. Schematic representation of studied reaction mechanisms and denotations used for corresponding reaction enthalpies.
Molecules 24 01646 sch001
Figure 1. Dietary vitamins; (a) vitamin A, (b) vitamin B1, (c) vitamin B3, (d) vitamin B6, and (e) vitamin C.
Figure 1. Dietary vitamins; (a) vitamin A, (b) vitamin B1, (c) vitamin B3, (d) vitamin B6, and (e) vitamin C.
Molecules 24 01646 g001
Figure 2. Geometries for (a) vitamin C, (b) vitamin B1, (c) vitamin B3, (d) vitamin B6, and (e) vitamin A in aqueous phase. Color atoms: for O, for C, for N, for S, for H.
Figure 2. Geometries for (a) vitamin C, (b) vitamin B1, (c) vitamin B3, (d) vitamin B6, and (e) vitamin A in aqueous phase. Color atoms: for O, for C, for N, for S, for H.
Molecules 24 01646 g002
Figure 3. Comparison of bond dissociation enthalpies (BDE) of the dietary vitamins.
Figure 3. Comparison of bond dissociation enthalpies (BDE) of the dietary vitamins.
Molecules 24 01646 g003
Figure 4. Comparison of ionization potentials (IP) of the dietary vitamins.
Figure 4. Comparison of ionization potentials (IP) of the dietary vitamins.
Molecules 24 01646 g004
Figure 5. Comparison of proton dissociation enthalpies (PDEs) of the dietary vitamins.
Figure 5. Comparison of proton dissociation enthalpies (PDEs) of the dietary vitamins.
Molecules 24 01646 g005
Figure 6. Comparison of proton affinities (PA) of the dietary vitamins.
Figure 6. Comparison of proton affinities (PA) of the dietary vitamins.
Molecules 24 01646 g006
Figure 7. Comparison of electron transfer enthalpies (ETE) of the dietary vitamins.
Figure 7. Comparison of electron transfer enthalpies (ETE) of the dietary vitamins.
Molecules 24 01646 g007
Table 1. Reaction enthalpies in kJ mol−1 in water for the selected dietary vitamins, employing mean functional values of M05-2X and M06-2X.
Table 1. Reaction enthalpies in kJ mol−1 in water for the selected dietary vitamins, employing mean functional values of M05-2X and M06-2X.
VitaminsHAT Mechanism (BDE)SET-PT Mechanism(IP + PDE)SPLET Mechanism(PA + ETE)
Vitamin A356.0 (±1.0)1742.0 (±2.0)532.0 (±1.0)
Vitamin B1435.0 (±0.0)1750.5 (±0.5)611.5 (±1.5)
Vitamin B3481.5 (±0.5)1790.5 (±1.5)658.0 (±2.0)
Vitamin B6360.5 (±1.5)1655.0 (±1.0)537.0 (±0.0)
Vitamin C318.5 (±0.5)1611.5 (±1.5)495.0 (±2.0)

Share and Cite

MDPI and ACS Style

Pandithavidana, D.R.; Jayawardana, S.B. Comparative Study of Antioxidant Potential of Selected Dietary Vitamins; Computational Insights. Molecules 2019, 24, 1646. https://doi.org/10.3390/molecules24091646

AMA Style

Pandithavidana DR, Jayawardana SB. Comparative Study of Antioxidant Potential of Selected Dietary Vitamins; Computational Insights. Molecules. 2019; 24(9):1646. https://doi.org/10.3390/molecules24091646

Chicago/Turabian Style

Pandithavidana, Dinesh R., and Samith B. Jayawardana. 2019. "Comparative Study of Antioxidant Potential of Selected Dietary Vitamins; Computational Insights" Molecules 24, no. 9: 1646. https://doi.org/10.3390/molecules24091646

APA Style

Pandithavidana, D. R., & Jayawardana, S. B. (2019). Comparative Study of Antioxidant Potential of Selected Dietary Vitamins; Computational Insights. Molecules, 24(9), 1646. https://doi.org/10.3390/molecules24091646

Article Metrics

Back to TopTop