Metrological aspects of a gas-phase DFT/B3LYP quantum-chemical approach to prioritize radical scavenging activity among a group of olive oil phenols

Abstract

Aim: A protocol relying on quantum chemical calculations to assist prioritization of phenolic compounds as antioxidants in terms of hydrogen atom donation efficiency is presented. The use of reference compounds, an important metrological issue for a future harmonization and standardization of computational approaches in research is also considered. Methods: A density functional theory (DFT) approach, namely B3LYP/6-311G++(2d,2p)//B3LYP/6-31G in the gas-phase was used for structure optimization, frequency calculation and single point energy (SPE) calculation to obtain the bond dissociation enthalpy (BDE) value of the most active O-H of olive oil phenols. For some of them used as a test set and for comparison, BDE values were calculated using three more approaches, M05-2X in the gas-phase, Becke three-parameter Lee-Yang-Parr (B3LYP) and implicit solvent effects (n-heptane to model bulk oils) with the integral equation formalism version of the polarizable continuum model (IEF-PCM), Minnesota 05 functional with double nonlocal exchange (M05-2X) at a single step using 6-31G+(d) basis set and solvation model density (SMD) as solvation model. Phenol and Trolox were used as reference compounds for ΔBDE calculation. Results: The proposed protocol was faster by 1.35-, 1.6-, and 8.3-fold respectively than the other three and provided almost the same activity trend application to other type of olive oil phenols indicated that prioritization based on ΔBDE values was in accordance with the limited existing experimental findings in bulk oils, and the order of activity was generally in agreement with the structure-antioxidant activity criteria. Conclusions: Present findings highlighted the usefulness of quantum chemical calculations as a tool to screen/prioritize molecules with an established structure saving experimental effort and waste production. The expression of results relatively to phenol and Trolox BDE values, may assist comparisons among research findings and facilitate standardization. Based on the findings hydroxytyrosol and related compounds should be efficient hydrogen atom donors compared to other potent virgin olive oil phenols.

Keywords

Antioxidant activity assessment / phenols / DFT / B3LYP / bond dissociation enthalpy / gas-phase / radical scavenging / hydroxytyrosol

References

[1] Nenadis N, Pyrka I, Tsimidou MZ. The Contribution of Theoretical Prediction Studies to the Antioxidant Activity Assessment of the Bioactive Secoiridoids Encountered in Olive Tree Products and By—Products. Molecules. 2023; 28: 2267.

[2] Silva AS, Nabavi SF, Saeedi M, Nabavi SM, editors. Recent advances in natural products analysis. Amsterdam: Elsevier; 2020.

[3] Gathungu RM, Kautz R, Kristal BS, Bird SS, Vouros P. The integration of LC—MS and NMR for the analysis of low molecular weight trace analytes in complex matrices. Mass Spectrom Rev. 2020; 39: 35-54.

[4] Nenadis N, Tsimidou MZ. 14—Assessing the activity of natural food antioxidants. In: Decker EA, editor. Oxidation in foods and beverages and antioxidant applications. Woodhead Publishing; 2010. pp. 332-67.

[5] Leopoldini M, Marino T, Russo N, Toscano M. Antioxidant properties of phenolic compounds: H—atom versus electron transfer mechanism. J Phys Chem A. 2004; 108: 4916-22.

[6] Tsao R, Li H. Antioxidant properties in vitro and in vivo: Realistic assessments of efficacy of plant extracts. CABI Reviews. 2012; 28: 1-9.

[7] Bursch M, Mewes JM, Hansen A, Grimme S. Best—Practice DFT Protocols for Basic Molecular Computational Chemistry. Angew Chem Int Ed Engl. 2022; 61: e202205735.

[8] Spiegel M. Current Trends in Computational Quantum Chemistry Studies on Antioxidant Radical Scavenging Activity. J Chem Inf Model. 2022; 62: 2639-58.

[9] Galano A, Alvarez—Idaboy JR. A computational methodology for accurate predictions of rate constants in solution: application to the assessment of primary antioxidant activity. J Comput Chem. 2013; 34: 2430—45.

[10] Sytar O, Brestic M, Rai M, Shao HB. Plant phenolic compounds for food, pharmaceutical and cosmetics production. J Med Plants Res. 2012; 6: 2526-39.

[11] Mahfuz S, Shang Q, Piao X. Phenolic compounds as natural feed additives in poultry and swine diets: a review. J Anim Sci Biotechnol. 2021; 12: 48.

[12] Zagoskina NV, Zubova MY, Nechaeva TL, Kazantseva VV, Goncharuk EA, Katanskaya VM, et al. Polyphenols in Plants: Structure, Biosynthesis, Abiotic Stress Regulation, and Practical Applications (Review). Int J Mol Sci. 2023; 24: 13874.

[13] Frankel EN. Lipid oxidation. 2nd ed. Oxford: Woodhead Publishing; 2005.

[14] Regulation (EU) No 432/2012 of 16 May 2012 establishing a list of permitted health claims made on foods, other than those referring to the reduction of disease risk and to children’s development and health. Off J Eur Union. 2012; L136.

[15] Servili M, Sordini B, Esposto S, Urbani S, Veneziani G, Di Maio I, et al. Biological Activities of Phenolic Compounds of Extra Virgin Olive Oil. Antioxidants (Basel). 2013; 3: 1-23.

[16] Frisch MJ, Trucks GW, Schlegel HB, Scuseria GE, Robb MA, Cheeseman JR, et al. Gaussian 16. Revision C.02[software]. Gaussian, Inc.; 2016 [cited 2023 Jun 20]. Available from: https://gaussian.com/gaussian16/

[17] Becke AD. Density—functional thermochemistry I. The effect of the exchange—only gradient correction. J Chem Phys. 1992; 96: 2155-60.

[18] Lee C, Yang W, Parr RG. Development of the Colle—Salvetti correlation—energy formula into a functional of the electron density. Phys Rev B Condens Matter. 1988; 37: 785-9.

[19] McLean AD, Chandler GS. Contracted Gaussian basis sets for molecular calculations. I. Second row atoms, Z=11—18. J Chem Phys. 1980; 72: 5639-48.

[20] Krishnan R, Binkley JS, Seeger R, Pople JA. Self—consistent molecular orbital methods. XX. A basis set for correlated wave functions. J Chem Phys. 1980; 72: 650-4.

[21] Zhao Y, Schultz NE, Truhlar DG. Exchange—correlation functional with broad accuracy for metallic and nonmetallic compounds, kinetics, and noncovalent interactions. J Chem Phys. 2005; 123: 161103.

[22] Scalmani G, Frisch MJ. Continuous surface charge polarizable continuum models of solvation. I. General formalism. J Chem Phys. 2010; 132: 114110.

[23] Marenich AV, Cramer CJ, Truhlar DG. Universal solvation model based on solute electron density and on a continuum model of the solvent defined by the bulk dielectric constant and atomic surface tensions. J Phys Chem B. 2009; 113: 6378—96.

[24] Dávalos JZ, Valderrama—Negrón AC, Barrios JR, Freitas VLS, Ribeiro da Silva MDMC. Energetic and Structural Properties of Two Phenolic Antioxidants: Tyrosol and Hydroxytyrosol. J Phys Chem A. 2018; 122: 4130—7.

[25] Wright JS, Johnson ER, DiLabio GA. Predicting the activity of phenolic antioxidants: theoretical method, analysis of substituent effects, and application to major families of antioxidants. J Am Chem Soc. 2001; 123: 1173-83.

[26] Nenadis N, Sigalas MP. A DFT study on the radical scavenging activity of maritimetin and related aurones. J Phys Chem A. 2008; 112: 12196-202.

[27] Papadopoulos G, Boskou D. Antioxidant effect of natural phenols on olive oil. J Amer Oil Chem Soc. 1991; 68: 669-71.

[28] Bountagkidou OG, Ordoudi SA, Tsimidou MZ. Structure—antioxidant activity relationship study of natural hydroxybenzaldehydes using in vitro assays. Food Res Inter. 2010; 43: 2014—9.

[29] Cuvelier ME, Richard H, Berset C. Comparison of the antioxidative activity of some acid—phenols: structure—activity relationship. Biosci Biotechnol Biochem. 1992; 56: 324—5.

[30] Gordon MH, Paiva—Martins F, Almeida M. Antioxidant activity of hydroxytyrosol acetate compared with that of other olive oil polyphenols. J Agric Food Chem. 2001; 49: 2480-5.

[31] Hrádková I, Šmidrkal J, Filip V, Merkl R, Kabrdová E. Antioxidant stability of phenolic acids and their esters. Czech J Food Sci. 2009; 27: S41-4.

[32] Nenadis N, Tsimidou MZ. Contribution of DFT computed molecular descriptors in the study of radical scavenging activity trend of natural hydroxybenzaldehydes and corresponding acids. Food Res Inter. 2012; 48: 538—43.

[33] Tsimidou MZ, Nenadis N, Zhang H. Structure—radical scavenging activity relationships of flavonoids and phenolic acids. In: Boskou D, Gerothanasis IP, Kefalas P, editors. Natural antioxidant phenols: sources, structure—activity relationship, current trends in analysis and characterisation. Trivandrum (India): Research Singpost; 2006. pp. 29-67.

[34] Angelis A, Antoniadi L, Stathopoulos P, Halabalaki M, Skaltsounis LA. Oleocanthalic and oleaceinic acids: New compounds from extra virgin olive oil (EVOO). Phytochem Lett. 2018; 26: 190—4.

[35] Boskou D. Sources of natural phenolic antioxidants. Trends Food Sci Technol. 2006; 9: 505-12.

[36] Bendini A, Cerretani L, Carrasco—Pancorbo A, Gómez—Caravaca AM, Segura—Carretero A, Fernández—Gutiérrez A, et al. Phenolic molecules in virgin olive oils: a survey of their sensory properties, health effects, antioxidant activity and analytical methods. An overview of the last decade. Molecules. 2007; 12: 1679-719.

[37] Amine Ajal E, Chaji S, Moussafir S, Nejjari R, Soulaymani A, Bajoub A. Virgin Olive Oil Phenolic Compounds: Insights on Their Occurrence, Health—Promoting Properties and Bioavailability. In: Akram M, editor. Olive Oil. Rijeka: IntecOpen; 2021.

[38] Ramachandran KI, Deepa G, Namboori DK. Computational chemistry and molecular modeling: Principles and applications. Heidelberg: Springer Berlin; 2008.

[39] AOAC SMPR 2011.011 Standard Method Performance Requirements for in vitro Determination of Total Antioxidant Activity in Foods, Beverages, Food Ingredients, and Dietary Supplements. J AOAC Inter. 2012; 95: 1577.

[40] Zhang HY, Sun YM, Chen DZ. O—H bond dissociation energies of phenolic compounds are determined by field/inductive effect or resonance effect? A DFT study and its implication. Mol Informatics. 2001; 20: 148-52.