Theoretical Prediction Strategies to the Antioxidant Activity Assessment: Comparison
Please note this is a comparison between Version 1 by Nikolaos Nenadis and Version 2 by Camila Xu.

Theoretical examination of the antioxidant activity of the individual or selected compounds identified in a matrix can be a green mid-tool for their prioritization regarding antioxidant behavior before seeking experimental in vitro or in vivo data.

  • antioxidant activity
  • theoretical predictions
  • phenolic compounds
  • olive tree products
  • olive tree by-products
  • oleuropein
  • hydroxytyrosol
  • tyrosol
  • secoiridoids
  • density functional theory

1. Introduction

The interest in natural compounds that exhibit antioxidant activity in food matrices [1] is extended also to the benefits such compounds can exert on humans or animal species after consumption [2][3][2,3]. For example, the high oxidative stability of virgin olive oil (VOO) at first was attributed to its rich in monounsaturated fatty acids character and then to the critical role of certain groups of intrinsic antioxidants such as the phenolic compounds present in the polar fraction of the oil. For more than two decades, the interest has been focused on the in vivo properties of the secoiridoids oleuropein, ligstroside, and related derivatives, which can be found not only in the oil but also in other edible products of the olive tree or can be isolated from by-products of the olive industry [4][5][4,5].
Antioxidants, in general, exhibit their activity even when present in minute concentrations. The phenolic-type antioxidants are considered primary antioxidants because they react directly with reactive oxygen and nitrogen species via hydrogen atom and/or electron donation and have gained the interest of food, medical, pharmaceutical, and cosmeceutical researchers all over the world [6].

2. Theoretical Prediction Strategies

As a result of advances in computer science and software engineering, it is feasible to accomplish theoretical calculations for compounds of medium or even larger size in a reasonable amount of time, making them suitable even for routine studies. The overall number of publications regarding antioxidants that involve theoretical calculations is steadily increasing, but the latter have not found their place in the agriculture and food sector yet. This is indicated by the bibliometrics shown in Figure 1.
Figure 1. (a) Publications categorized by subject area and (b) total publications per year appearing in Scopus database using “density functional theory” and “antioxidant” as keywords (accessed 21 November 2022, total documents, 995).
In most of the relevant publications, emphasis is given to the study of the radical scavenging process that coincides with the major mechanisms of action of primary antioxidants, such as the phenolic compounds [7][8][9][10][11][12][13][14][7,15,16,17,18,19,20,21]. Different strategies can be applied to theoretical studies so that the derived data are reported differently, for example, in terms of intrinsic reactivity indices, or using thermochemical or kinetic parameters. In all cases, knowledge of the possible reactions involved and the factors that may influence the reactivity (e.g., environment) is a prerequisite [15][22]. The theoretical approach usually begins with the conformer analysis to locate the most stable compound (global minimum). This requires the examination of all possible structures, often including the determination of enthalpy barriers upon, e.g., rotation of one or more dihedral angles. More elaborate is the potential energy surface (PES) examination, which results in a plot of the energy versus geometry with concomitant rotation of dihedrals. Such an approach can be useful for equilibrium geometries determination and transition states [16][17][23,24]. An example of targeted PES examination with a concomitant rotation of two dihedrals is given in Figure 2 for the ligstroside derivative tyrosol. The contour map indicates the place of global minimum, local minima, and destabilization.
Figure 2. (a) The varying dihedrals of tyrosol (SC1: C5-C6-C8-C9, SC2: C8-C9-O10-H20) selected for the potential energy surface (PES) scan at B3LYP/6-31G in gas phase and (b) 3D plot of PES [18][25].
Certainly, it is possible to locate more than one conformer, in case the difference in enthalpy from the global minimum is small, suggesting that interconversion can be feasible. Except for the rotation of groups and side chains, tautomerism and shifting from E to Z configuration or bending of dihedral angles are features to be also explored, as highlighted by Spiegel [7]. All of them can affect hydrogen bonding, electron delocalization, polarizability, dipole moment value, as well as steric hindrance phenomena. Considering that an optimum structure has been obtained, then a first-cut approach resides in the frontier orbital examination. The highest molecular orbital (HOMO) energy is related to electron donation, and the lowest unoccupied molecular orbital (LUMO) to electron acceptance [7][19][7,26]. Thus, the energy gap between the two orbitals, electronegativity (χ), chemical potential (μ), global hardness (η), global softness (S), and electrophilicity (ω) can be calculated to have an estimate of molecular reactivity. However, electron correlation (a measure of the influence in an electron movement due to the presence of all other electrons in the system), which is critical for accurate and quantitative evaluation of molecular energies, is not considered and may affect the validity of conclusions [19][26]. In addition, visualization of HOMO may serve in detecting the structural site of the attack, even though the electrophilic/nucleophilic sites are better revealed via solving the Fukui functionals [20][27] and more accurately by applying the dual descriptor proposed by Martínez-Araya [21][28]. Another approach is the examination of the so-called intrinsic reactivity of the antioxidants, which is most frequently employed by different research groups [15][22]. Such an approach requires the calculation of enthalpy of derived species, namely phenoxy radical, cation radical, anion, hydrogen atom, and proton, even of further derivatives depending on whether the structure favors more than one proton or electron donation in a stepwise manner (e.g., double-HAT, double-SPLET, SPLHAT) [14][21], which in turn are used to calculate the indices given in Table 1.
Table 1. Common molecular indices related to the radical scavenging potential of compounds and the formulae for their calculation.