Next Article in Journal
Exploring miRNAs’ Based Modeling Approach for Predicting PIRA in Multiple Sclerosis: A Comprehensive Analysis
Previous Article in Journal
Structural Insights into the Mechanisms Underlying Polyaminopathies
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Proton-Coupled Electron Transfer and Hydrogen Tunneling in Olive Oil Phenol Reactions

by
Jelena Torić
1,
Ana Karković Marković
1,
Stipe Mustać
1,
Anamarija Pulitika
2,
Cvijeta Jakobušić Brala
1,* and
Viktor Pilepić
1,*
1
Faculty of Pharmacy and Biochemistry, University of Zagreb, 10000 Zagreb, Croatia
2
Faculty of Chemical Engineering and Technology, University of Zagreb, 10000 Zagreb, Croatia
*
Authors to whom correspondence should be addressed.
Int. J. Mol. Sci. 2024, 25(12), 6341; https://doi.org/10.3390/ijms25126341
Submission received: 29 April 2024 / Revised: 4 June 2024 / Accepted: 5 June 2024 / Published: 7 June 2024
(This article belongs to the Section Bioactives and Nutraceuticals)

Abstract

:
Olive oil phenols are recognized as molecules with numerous positive health effects, many of which rely on their antioxidative activity, i.e., the ability to transfer hydrogen to radicals. Proton-coupled electron transfer reactions and hydrogen tunneling are ubiquitous in biological systems. Reactions of olive oil phenols, hydroxytyrosol, tyrosol, oleuropein, oleacein, oleocanthal, homovanillyl alcohol, vanillin, and a few phenolic acids with a DPPH• (2,2-diphenyl-1-picrylhydrazyl) radical in a 1,4-dioxane:water = 95:5 or 99:1 v/v solvent mixture were studied through an experimental kinetic analysis and computational chemistry calculations. The highest rate constants corresponding to the highest antioxidative activity are obtained for the ortho-diphenols hydroxytyrosol, oleuropein, and oleacein. The experimentally determined kinetic isotope effects (KIEs) for hydroxytyrosol, homovanillyl alcohol, and caffeic acid reactions are 16.0, 15.4, and 16.7, respectively. Based on these KIEs, thermodynamic activation parameters, and an intrinsic bond orbital (IBO) analysis along the IRC path calculations, we propose a proton-coupled electron transfer mechanism. The average local ionization energy and electron donor Fukui function obtained for the phenolic compounds show that the most reactive electron-donating sites are associated with π electrons above and below the aromatic ring, in support of the IBO analysis and proposed PCET reaction mechanism. Large KIEs and isotopic values of Arrhenius pre-exponential factor AH/AD determined for the hydroxytyrosol, homovanillyl alcohol, and caffeic acid reactions of 0.6, 1.3, and 0.3, respectively, reveal the involvement of hydrogen tunneling in the process.

Graphical Abstract

1. Introduction

Oxidative stress, often induced by an overproduction of free radicals, presents a high risk to human health [1,2]. Free radicals, as reactive compounds, initiate chain reactions, leading to damage to biological molecules such as lipids, proteins, and DNA. Oxidative stress within a lipid environment is particularly important, as lipids react with oxidants faster than proteins and DNA. Primary antioxidants react directly with free radicals, leaving less reactive radicals unable to damage biological molecules. Natural products are an important source of antioxidants, containing a vast diversity of molecules. The study of the relative antioxidative activity of these compounds is important not only to understand their protective effects but also to design efficient strategies against oxidative stress.
Olive oil, as a cornerstone of the Mediterranean diet, and olive oil phenols have attracted considerable attention in recent decades due to their beneficial pharmacological, medicinal, and biochemical properties [3,4,5,6,7]. The molecules that stand out are hydroxytyrosol [8] and its more complex derivatives oleacein [9] and oleocanthal [10,11], which exhibit cardioprotective, neuroprotective, anticancer, immunomodulatory, and other effects, many of which rely on their ability to react as primary antioxidants. The antioxidant activity of olive oil phenols has been extensively studied using various chemical approaches. These include different tests of antioxidative activity (such as DPPH, ORAC, and ABTS assays), accelerated oxidation in a lipid model system (OSI, oxidative stability index), electrochemical methods (flow injection analysis, amperometry, and cyclic voltammetry), calculating the bond dissociation enthalpy and ionization potential of phenolic hydroxy groups, and measuring the effect on LDL oxidation [12,13]. Their positive effects on lipid oxidation are particularly recognized. The European Food Safety Authority has published a health claim regarding the role of olive oil phenols in protecting LDL from oxidation in vivo: “A daily intake of 20 g of olive oil, which contains at least 5 mg of hydroxytyrosol or its derivatives (e.g., oleuropein and tyrosol) provides beneficial effects” [14].
A kinetic approach is very useful for studying trends in antioxidative activity. The rate constant is a measure of reactivity; the faster the reaction, the higher the antioxidant activity. Kinetics also accounts for reaction mechanisms and other important features, such as the tunneling effect [15,16]. Additionally, computational chemistry methods are often used in the study of antioxidative activity [17].
The antioxidative activity of phenol (ArOH) is related to its ability to transfer a hydrogen to a radical (R•) (Equation (1)) much faster than the chain-propagating reaction of peroxidation of organic compounds (X) (Equation (2)) [15,18]:
R• + ArOH → RH + ArO•
R• + X → RH + X•
In reactions with R•, ArOH converts to a phenoxyl radical (ArO•) and gives an unreactive species (RH). The formed phenoxyl radicals are stabilized by resonance, conjugation, and delocalization effects [15,19]. Reactions of phenols with DPPH• (2,2-diphenyl-1-picrylhydrazyl radical) are often used to evaluate the antioxidant activity of phenols [15,18,20,21,22,23]. They serve as a prototype model for the reactions of phenols with peroxyl radicals, ROO•, which are extremely important in biological systems. However, the study of reactions with peroxyl radicals is rather complex because it is necessary to produce a precursor, and reactions are often rapid [24,25]. DPPH• is commercially available, stable, and possesses a strong absorption band in the visible region (λmax ≈ 520 nm), making it suitable for an experimental study.
Formal hydrogen atom transfer (electron and proton transfer) from phenol to radical (Equation (1)) may proceed through different mechanisms [15]. Electron transfer (ET) and proton transfer (PT) can occur together in a single chemical step, such as proton-coupled electron transfer (PCET) or hydrogen atom transfer (HAT). In PCET, an electron and a proton are transferred in a single step as two separated particles, while in HAT, an electron and a proton are transferred together as a single entity (a hydrogen atom). The main distinction between these mechanisms is that in HAT, the donor and acceptor sites are the same for both electron and proton, while in PCET, an electron and a proton are transferred from different sets of orbitals. The HAT and PCET mechanisms can be differentiated using computational chemistry methods [26,27]. Another possibility is a consecutive mechanism, such as sequential proton loss electron transfer (SPLET), ET/PT, or PT/ET. Generally, HAT/PCET mechanisms are proposed as predominant in nonpolar solvents, while SPLET is favored in polar solvents due to increased dissociation of phenolic OH groups and greater reactivity of phenolate anions [12,18]. The reactions of phenols with peroxyl radicals (ROO•) usually proceed by a concerted (PCET/HAT) mechanism; thus, reaction rate constants in nonpolar solvents, such as dioxane, are the only useful parameters for predicting the antioxidant ability of ArOH [22]. Although dioxane does not represent a physiological medium taking into account its polarity, it could represent a reaction in physiological conditions because the expected reaction mechanism in dioxane is HAT/PCET [12,18], which is the same mechanism that is proposed for ROO• radicals [22]. However, the reaction of phenols with DPPH•, widely used to evaluate the antioxidative properties of synthetic and natural phenols, is typically performed in methanol, where the proposed mechanism is SPLET. This raises questions about the validity of the conclusions drawn from these studies.
Hydrogen transfer reactions of phenols have been extensively studied as model reactions, mimicking reactions that involve biological oxidation of tyrosine [28]. Borden et al. used the self-exchange reaction between phenol and phenoxyl radicals as a prototype reaction to demonstrate the difference between PCET and HAT reactions [29]. The PCET mechanism is ubiquitous in biological systems [30,31,32], and it has been proposed for numerous reactions of phenols as well as other antioxidants such as vitamin C and vitamin E, enzyme reactions, the cell respiration process, and complex processes during photosynthesis. PCET is also a reaction mechanism in many chemical and electrochemical processes relevant to solar energy devices [33,34,35]. An interesting recent question is the role of the inverted region in biological pathways [36,37].
Another important feature of hydrogen transfer reactions is hydrogen tunneling [38,39,40], which has been observed in PCET reactions of natural antioxidants such as vitamin C and vitamin E [41,42,43]. Nakanishi et al. have recently observed hydrogen tunneling in the reaction between α-tocopherol and DPPH• [43].
In this study, we present an experimental and computational kinetic study of the reactions of several olive oil phenols with DPPH• in a nonpolar solvent, 1,4-dioxane. We examined the reactions of olive oil phenols, well-recognized as important bioactive compounds, including the simple phenols hydroxytyrosol [44,45], tyrosol [46], homovanillyl alcohol [47,48], and their more complex derivatives, oleuropein [49], oleacein [9], oleocanthal [50,51], and phenolic acids, including caffeic acid [52], p-coumaric acid [53], ferulic acid [54], vanillic acid [55], and vanillin. We have determined the reaction rate constants and proposed the PCET mechanism with the involvement of hydrogen tunneling. To the best of our knowledge, this is the first experimental kinetic study of the antioxidative reactions of the olive oil phenols hydroxytyrosol, tyrosol, oleuropein, oleacein, and oleocanthal.

2. Results and Discussion

In this study, we investigated the antioxidant reactions of simple olive oil phenols, hydroxytyrosol (HOTyr), tyrosol, homovanillyl alcohol (HVA), their more complex derivatives, oleuropein, oleacein, oleocanthal, and phenolic acids, caffeic acid, p-coumaric acid, ferulic acid, vanillic acid, and vanillin. The structures of these compounds are presented in Figure 1 below. We analyzed reactions with DPPH• (2,2-diphenyl-1picrylhydrazyl radical) in a nonpolar solvent mixture, 1,4-dioxane:water = 99:1 and 95:5 v/v.

2.1. Reaction of Olive Oil Phenols with DPPH•

Rate constants for reactions of phenols with DPPH• were determined by following the disappearance of DPPH• radical at 520 nm (see Supplementary Material, Figures S11 and S12). Kinetic measurements were conducted under pseudo-first-order conditions, with phenol in excess at least ten times. Second-order reaction rate constants, kArOH, were obtained from the slopes of plots of the observed pseudo-first-order rate constants vs. the concentration of phenol (see Supplementary Materials, Tables S1–S7, and Figures S1–S7). For each analyzed phenol reaction, the dependence was strongly linear, in accordance with the rate law:
d D P P H d t = 1 2 k A r O H D P P H A r O H
The observed rate constant is as follows:
k o b s = 2 k A r O H A r O H
The same rate law has been determined previously for reactions of various phenols and flavonoids with DPPH• by others [15]. The overall stoichiometry of the reaction is ArOH:DPPH• = 1:2. The stoichiometry was determined by mixing a small quantity of phenol with an excess of DPPH•. This is related to factor 2 in the rate law. The obtained rate law is consistent with the proposed reaction mechanism [27].
In the first slow step (kArOH), phenol reduces the DPPH• radical, giving a phenoxyl radical (ArO•) and reduced H-DPPH (Equation (5)):
ArOH + DPPH• → ArO• + H-DPPH (kArOH)
The phenoxyl radical formed in the relatively slow first step either reduces another DPPH• or reacts with another molecule of ArO• in the subsequent fast reaction [56]:
ArO• + DPPH• → products
ArO• + ArO• → dimer or disproportionation products
The rate parameters were determined in 1,4-dioxane:water = 99:1 and 95:5 v/v solvent mixtures. Olive oil phenols are particularly recognized for their antioxidant activity in lipid protection, and such a nonpolar environment can be represented by dioxane as a nonpolar solvent. Solvent impurities can dramatically influence the experimental measurement of the rate parameters of these reactions, as previously observed and explained by the influence of traces of basic and acidic impurities in the solvent [57]. Dioxane appears to be the solvent in which these problems are present to the smallest extent.
Rate constants and kinetic isotope effects for the reactions of phenols with DPPH• in a 1,4-dioxane:water (99:1 or 95:5 v/v) solvent mixture are presented in Table 1 and activation parameters in Table 2.
The highest rate constants are observed for the ortho-diphenols HOTyr, oleuropein, and oleacein. HOTyr is a simple phenol, and oleuropein and oleacein are its more complex derivatives with the same aromatic ring structure (Figure 1). Tyrosol, which lacks one OH group compared to HOTyr, shows a rate constant approximately 300 times smaller than HOTyr. Oleuropein and oleacein, which are complex derivatives of HOTyr, show greater rate constants than oleocanthal, which is a derivative of tyrosol. HVA, with a metoxy group in place of one OH group, has a rate constant ~15 times smaller than that of HOTyr. Similar observations are found in the case of phenolic acids. Caffeic acid, a carboxylic acid with two OH groups in the ortho-position on the aromatic ring, has a rate constant ~200 times larger than that of p-coumaric acid, which lacks one OH group compared to caffeic acid. Ferulic acid, with a methoxy group in the ortho-position relative to the OH group, has a rate constant ~3 times smaller than that of caffeic acid.
The antioxidative activity of the phenols reacting with the hydrogen atom transfer mechanism can be predicted by considering the bond dissociation enthalpy (BDE) of the phenolic OH bond [15]. According to BDE, the most efficient antioxidants are ortho-diphenols, which are related to the stabilization of phenoxy radicals by internal H-bonds [19,58]. The presence of an OH group in the ortho-position decreases BDE by ~40 kJ/mol, while the presence of a methoxy group in the ortho-position slightly stabilizes the radical. The BDE values in the gas phase for HOTyr, tyrosol, and phenol are 307, 343, and 347 kJ/mol, respectively [25]. Tyrosol BDE indicates the stabilizing effect of a –CH2CH2OH chain in the para-position, and HOTyr BDE reflects the additional effect of the ortho-OH group. BDEs in water and in benzene solutions follow the same trend as in the gas phase [19]. Our experimentally determined rate constants follow the trend in antioxidant activities based on BDE [19].

2.2. Reaction Mechanism of Hydrogen Transfer

The reaction mechanism of the reaction between phenol and DPPH• (Equation (5)) involves electron (ET) and proton (PT) transfer. There are a few possible mechanisms for the transfer of electrons and protons. ET and PT could occur concertedly in a single reaction step (the HAT/PCET mechanism) or stepwise in two reaction steps such as PT/ET, ET/PT, or SPLET. The mechanism of ET and PT can be deduced by considering kinetic isotope effects (KIEs) and thermochemical analysis [59].
KIEs in the reactions of phenols with DPPH• were measured in reaction solvent mixtures 1,4-dioxane:H2O (kArOH) and 1,4-dioxane:D2O (kArOD) (v/v = 99:1; 95:5).
KIE = kArOH/kArOD
ArOD was obtained by dissolving phenol in a 1,4-dioxane-D2O solvent mixture due to the fast exchange of phenol OH hydrogens with deuterium. kArOD was determined in the same way as kArOH. The kinetic isotope effects are presented in Table 1.
In all the analyzed reactions, we have determined the KIEs. The observed KIEs clearly indicate the involvement of proton transfer in the rate-determining step. In the sequential ET/PT and SPLET mechanisms, the rate-determining step is ET, and so the sequential ET/PT and SPLET mechanisms can be excluded based on the determined KIE.
Burton et al. determined KIEs in reactions of substituted phenols and peroxyl radicals: 2,3,5,6-tetramethyl-4-methoxyphenol (10.6), 2,6-di-tert-butyl-4-methoxyphenol (9.4), 2,6-di-tert-butyl-4-methylphenol (6.8), and α-tocopherol (4.0) [24]. Small KIEs (23.3) were also determined in reactions of substituted phenols with DPPH• in cyclohexane [27] and α-tocopherol with DPPH• in water (6.6) [43].
Another possibility is a consecutive mechanism, where the first step is a slow PT, followed by a fast ET. This can be rationalized through thermochemical analysis. In the thermochemical analysis, the standard reaction Gibbs energy for initial PT in the sequential PT/ET is compared to the activation Gibbs energy determined experimentally for the corresponding reaction, with the assumption that the activation Gibbs energy is always at least equal to or greater than the standard reaction Gibbs energy [31].
We have determined experimentally activation Gibbs energy and other thermodynamic activation parameters from the temperature dependence of rate constants in a range from 10 to 45 °C (Supplementary Materials, Tables S8–S10 and Figures S8–S10).
The determined activation Gibbs energies for the reactions of HOTyr, HVA, and caffeic acid are 69.0, 82.5, and 74.7 kJ/mol, respectively. A consecutive PT/ET mechanism cannot be excluded as a possible mechanism if the corresponding ΔrGѳ is below the value of the activation Gibbs energy. ΔrGѳ can be calculated from the pKa value of the corresponding phenol as ΔrGѳ = −RTlnKa. pKa values of the phenolic OH group in water are about 10 or greater. In nonpolar solvents like dioxane, it is reasonable to assume that the pKa values are significantly higher than those in water. Rossini et al. studied the influence of nonpolar solvents on the pKa value of different phenols and predicted an increase in pKa from the value of 10 in water to 18 in methanol and 26 in acetonitrile [60]. It can be expected that pKa in 1,4-dioxane is even higher, as it is a more nonpolar solvent than acetonitrile. The highest values of ΔrGѳ for HOTyr, HVA, and caffeic acid reactions, which could be expected in the case of PT/ET, are 69.0, 82.5, and 74.7 kJ/mol and correspond to the pKa values of 12.1, 14.5, and 13.1, respectively. Considering the previously mentioned assumption for the pKa values in 1,4-dioxane, it can be predicted that ΔrGʅ is definitely above these values, i.e., ΔG is lower than ΔrGѳ. Thus, a sequential PT/ET mechanism can also be excluded.
The obtained Arrhenius activation parameters are in the same range as those previously reported by Foti et al. for reactions of substituted phenols with DPPH• [27].
Altogether, it can be concluded that the mechanism of PT and ET in reactions of analyzed phenols with DPPH• is the concerted PCET or HAT mechanism.

2.3. Hydrogen Tunneling

A large KIE, above the value predicted semi-classically, is a typical sign of the involvement of hydrogen tunneling in the reaction [61]. According to Bell, the upper limit for KIE in the case of an OH bond is 8, calculated considering stretching vibrations, and 13 with bending vibrations included. The appearance of a KIE value greater than the upper limit is clear evidence of hydrogen tunneling. In the reactions of HOTyr, HVA, and caffeic acid, we have determined KIEs of 16.0, 15.4, and 16.7, respectively.
It is interesting that in the case of the structurally similar phenols HOTyr (16.0) and tyrosol (3.3), we have determined significantly different KIEs. In addition, there is a great difference in the rate constants between HOTyr and Tyr, with the HOTyr rate constant being 300-fold greater. This significant difference in reactivity seems to be related to the involvement of hydrogen tunneling in the process [61].
KIEs and hydrogen tunneling have been observed earlier in antioxidative reactions, such as the reaction of vitamin C with TEMPO• (KIE~30) [42] or hexacyanoferrate (III) [62].
An additional sign of the involvement of hydrogen tunneling in the process is the deviation in isotopic values of thermodynamic activation parameters AH/AD and ΔEa from semi-classically predicted values [61,63]. The semi-classical limits of AH/AD in a hydrogen-transfer process are 0.7–1.2, and the upper limit for ΔEa value is 5.1 kJ/mol, calculated from the difference between zero-point energies E0(D) − E0(H) for the dissociation of the OH bond. The isotopic values of thermodynamic activation parameters are presented in Table 3.
The observed AH/AD in the reactions of HOTyr, HVA, and caffeic acid are 0.6, 1.3, and 0.3, respectively, all of which fall outside of the predicted semiclassical limits, supporting the proposed hydrogen tunneling. This is further corroborated by the observed isotopic differences in the activation energy ΔEa, 8.1, 5.9, and 9.5 for HOTyr, HVA, and caffeic acid, respectively, also beyond the semiclassical limit of 5.1 kJ/mol. Hydrogen tunneling appears to be an important feature of the reactions of phenols with prominent antioxidative activity.

2.4. Computational Analysis

The calculations of transition state structures (TS) for the proposed hydrogen transfer from HOTyr, HVA, and caffeic acid to DPPH• were performed and further explored with an intrinsic reaction coordinate (IRC) [64] and intrinsic bond orbital (IBO) analysis [65,66,67]. The obtained results are summarized in Figure 2, Figure 3 and Figure 4 and Table 4 and Table 5, and also in the Supplementary Materials, Figures S13–S17 and Table S11. Compared to related TS structures for hydrogen transfer from phenol or p-metoxyphenol to DPPH• obtained by Foti et al. [27], the obtained TS structures for HVA (Figure 2), HOTyr, and caffeic acid (Figure S13) have slightly shorter H⋯O and a bit longer N⋯H distances (Table 4). The calculated Gibbs activation energy (Table 5) agrees very well for HVA with those obtained from the experiment (Table 2) in support of the proposed TS structure. The reasonable agreement between the calculated and experimental Gibbs activation energies for HOTyr and caffeic acid can be explained by the presence of a second OH group. The calculated KIEs, when corrected for tunneling, are close to the experimentally obtained values (Table 5).
Intrinsic bond orbital (IBO) analysis proposed by Knizia [66] enables the visualization of molecular orbitals and electronic structure changes in an intuitive way along reaction paths [65,66,67]. The IBO analysis proved to be successful in the interpretation of reaction mechanisms and also in reactions of hydrogen transfer, where it was shown that it can unambiguously distinguish between HAT and PCET processes [67,69,70]. The representative intrinsic bond orbitals (IBOs) involved in hydrogen transfer from HVA, HOTyr, and caffeic acid to DPPH• plotted along the intrinsic reaction coordinate (IRC) obtained from the DFT calculations are shown in Figure 3, Figures S14 and S15, respectively. For IBO analysis along the entire IRC, see also the corresponding video files in the Supplementary Materials. In all three cases, the electron transferred to the DPPH• N atom comes from the aromatic ring orbital of the phenolic compound, leaving an unpaired electron on it, while the proton from the phenolic OH group is transferred separately to the orbitals located on another N atom of DPPH•. It can be concluded that in these cases, proton transfer (PT) and electron transfer (ET) are separate from different reaction sites and involve different orbitals, which clearly points to PCET reactions.
Average local ionization energy Ī(r) [71] and electron donor Fukui function f¯(r) [72] are local reactivity descriptors used to identify molecular reactivity sites for reactions with electrophiles or radicals in a series of studies [73,74,75,76,77,78]. The Ī(r) and f¯(r) obtained for HVA are shown in Figure 4 and for HOTyr and caffeic acid in Figures S16 and S17. For all of these phenolic compounds, the most reactive electron-donating sites (blue areas) are associated with π electrons above and below the aromatic rings, corresponding to the phenol orbital from which the electron is transferred to DPPH•, as obtained in the IBO analysis.

3. Materials and Methods

1,4-dioxane, 2,2-diphenyl-1-picrylhydrazyl (DPPH•), homovanillyl alcohol, tyrosol, caffeic acid, p-coumaric acid, HCl (1 M, Titrisol), and DCl (99% D) were obtained from Merck, KGaA, (Darmstadt, Germany), and vanillic acid, ferulic acid, and vanillin were from Fluka Chemie GmbH (Buchs, Switzerland). Hydroxytyrosol and oleuropein were purchased from Extrasynthese (Genay, France), and oleacein and oleocanthal from PhytoLab GmbH & Co. KG (Vestenbergsgreuth, Germany). All chemicals used were of analytical grade. Water and heavy water (Aldrich, 99.9%) were twice distilled and carbon dioxide- and oxygen-free (bubbled with 99.999% N2).
Kinetic Measurements. Pseudo-first-order rate constants for the reaction of phenols with the DPPH• radical have been determined spectrophotometrically by monitoring the decrease in absorbance of DPPH• at 525 nm. An Avantes AvaSpec 2048L StarLine spectrometer (Avantes B.V., Apeldoorn, The Netherlands) provided with a Quantum Northwest QPOD Temperature-Controlled Sample Compartment for Fiber Optic Spectroscopy was used throughout to collect spectral and absorbance-time data. Kinetic measurements were performed under carefully maintained temperature conditions (within the limits of ±0.1 °C). In a typical kinetic run, at least 200 pairs of absorbance-time data were collected and fitted to the common least-squares algorithm. At least three to four observed pseudo-first rate constants were always used to calculate the corresponding rate parameters under the specified reaction conditions. The measurements were performed under pseudo-first-order conditions, taking phenol concentrations to be at least 15-fold in excess, except in the case of secoiridoids, which were 2.5- to 12-fold in excess at their lowest concentrations. Very good pseudo-first-order kinetics were obtained throughout. The stock solutions of reactants were usually prepared by dissolving phenol or DPPH• in 1,4-dioxane. The reaction solution was thermostated for 20 min before the kinetics were started. The reaction was initiated by the addition of an appropriate volume of DPPH• stock solution to a thermostated reaction solution. Measured kinetic isotope effects are corrected considering the content of H in the 1,4-dioxane-D2O solvent mixture.
Computational method details. The DFT calculations were carried out using the GAUSSIAN 16 software package [79] at the (U)B3LYP/6-311++G(2d,2p) level of theory, applying very tight convergence criteria. Non-specific solvent effects have been estimated by using the polarizable continuum model (PCM) of the self-consistent reaction field (SCRF) method [80] with 1,4-dioxane as a solvent. The stationary points obtained by geometry optimizations were confirmed either as minima or saddle points by vibrational analysis at the same theory level. The obtained frequencies were scaled by a factor of 0.9679 [68]. The transition state structures were further explored with an intrinsic reaction coordinate (IRC) analysis [64] at the same theory level. The obtained reaction path structures were examined with an intrinsic bond orbital (IBO) analysis utilizing PBE/def2-TZVPPD and the IBO localization procedure [65,66,67] with the IBOVIEW program [81]. Local reactivity descriptors, the average local ionization energy, and the electron donor Fukui function have been calculated from the wave functions in the checkpoint files obtained from the above B3LYP/6-311++G(2d,2p) calculations by the MULTIWFN program [82]. An estimation of the tunneling correction for the calculated reaction constants and KIE was performed from a vibration analysis using the Wigner method [83] with the TAMkin program [84].

4. Conclusions

Olive oil ortho-diphenols, hydroxytyrosol, oleuropein, and oleacein show the highest antioxidant activity, with the highest values of reaction rate constants. The electron and proton transfers in these reactions proceed via a proton-coupled electron transfer mechanism. Additionally, large kinetic isotope effects have been observed in the reactions of hydroxytyrosol, homovanillyl alcohol, and caffeic acid, indicating the involvement of hydrogen tunneling in the process. The intrinsic bond orbital (IBO) analysis along the IRC path calculations supports the proposed reaction mechanism. The average local ionization energy and electron donor Fukui function obtained for the phenolic compounds show that the most reactive electron-donating sites are associated with π electrons above and below the aromatic ring, in support of the IBO analysis and the proposed PCET reaction mechanism.

Supplementary Materials

The supporting information can be downloaded at https://www.mdpi.com/article/10.3390/ijms25126341/s1.

Author Contributions

Conceptualization, C.J.B. and V.P.; methodology, C.J.B. and V.P.; software, V.P.; validation, C.J.B., V.P. and A.K.M.; formal analysis, J.T., S.M. and A.P.; investigation, C.J.B., V.P. and J.T.; resources, C.J.B., V.P. and A.K.M.; data curation, J.T., V.P. and S.M.; writing—original draft preparation, C.J.B. and V.P.; writing—review and editing, C.J.B., V.P., J.T. and A.K.M.; visualization, S.M. and V.P.; supervision, C.J.B. and V.P.; project administration, C.J.B. and V.P.; funding acquisition, C.J.B. and V.P. All authors have read and agreed to the published version of the manuscript.

Funding

This work was supported by the project FarmInova (KK.01.1.1.02.0021), funded by the European Regional Development Fund.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

The original contributions presented in the study are included in the article/Supplementary Materials, further inquiries can be directed to the corresponding author/s.

Acknowledgments

This research was performed using the advanced computing service provided by the University of Zagreb University Computing Centre (SRCE).

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Lushchak, V.I. Free Radicals, Reactive Oxygen Species, Oxidative Stress and Its Classification. Chem.-Biol. Interact. 2014, 224, 164–175. [Google Scholar] [CrossRef]
  2. Jomova, K.; Raptova, R.; Alomar, S.Y.; Alwasel, S.H.; Nepovimova, E.; Kuca, K.; Valko, M. Reactive Oxygen Species, Toxicity, Oxidative Stress, and Antioxidants: Chronic Diseases and Aging. Arch. Toxicol. 2023, 97, 2499–2574. [Google Scholar] [CrossRef]
  3. Finicelli, M.; Di Salle, A.; Galderisi, U.; Peluso, G. The Mediterranean Diet: An Update of the Clinical Trials. Nutrients 2022, 14, 2956. [Google Scholar] [CrossRef]
  4. Bucciantini, M.; Leri, M.; Nardiello, P.; Casamenti, F.; Stefani, M. Olive Polyphenols: Antioxidant and Anti-Inflammatory Properties. Antioxidants 2021, 10, 1044. [Google Scholar] [CrossRef]
  5. Gorzynik-Debicka, M.; Przychodzen, P.; Cappello, F.; Kuban-Jankowska, A.; Marino Gammazza, A.; Knap, N.; Wozniak, M.; Gorska-Ponikowska, M. Potential Health Benefits of Olive Oil and Plant Polyphenols. Int. J. Mol. Sci. 2018, 19, 686. [Google Scholar] [CrossRef]
  6. Almanza-Aguilera, E.; Cano, A.; Gil-Lespinard, M.; Burguera, N.; Zamora-Ros, R.; Agudo, A.; Farras, M. Mediterranean Diet and Olive Oil, Microbiota, and Obesity-Related Cancers. From Mechanisms to Prevention. Semin. Cancer Biol. 2023, 95, 103–119. [Google Scholar] [CrossRef]
  7. Losada-Barreiro, S.; Sezgin-Bayindir, Z.; Paiva-Martins, F.; Bravo-Díaz, C. Biochemistry of Antioxidants: Mechanisms and Pharmaceutical Applications. Biomedicines 2022, 10, 3051. [Google Scholar] [CrossRef]
  8. Gervasi, F.; Pojero, F. Use of Oleuropein and Hydroxytyrosol for Cancer Prevention and Treatment: Considerations about How Bioavailability and Metabolism Impact Their Adoption in Clinical Routine. Biomedicines 2024, 12, 502. [Google Scholar] [CrossRef]
  9. Huang, Y.; Guan, Q.; Zhang, Z.; Wang, P.; Li, C. Oleacein: A Comprehensive Review of Its Extraction, Purification, Absorption, Metabolism, and Health Effects. Food Chem. 2024, 433, 137334. [Google Scholar] [CrossRef]
  10. Beauchamp, G.K.; Keast, R.S.J.; Morel, D.; Lin, J.; Pika, J.; Han, Q.; Lee, C.-H.; Smith, A.B.; Breslin, P.A.S. Ibuprofen-like Activity in Extra-Virgin Olive Oil: Phytochemistry. Nature 2005, 437, 45–46. [Google Scholar] [CrossRef]
  11. González-Rodríguez, M.; Ait Edjoudi, D.; Cordero-Barreal, A.; Farrag, M.; Varela-García, M.; Torrijos-Pulpón, C.; Ruiz-Fernández, C.; Capuozzo, M.; Ottaiano, A.; Lago, F.; et al. Oleocanthal, an Antioxidant Phenolic Compound in Extra Virgin Olive Oil (EVOO): A Comprehensive Systematic Review of Its Potential in Inflammation and Cancer. Antioxidants 2023, 12, 2112. [Google Scholar] [CrossRef]
  12. Rietjens, S.J.; Bast, A. New Insights into Controversies on the Antioxidant Potential of the Olive Oil Antioxidant Hydroxytyrosol. J. Agric. Food Chem. 2007, 55, 7609–7614. [Google Scholar] [CrossRef]
  13. Bendini, A.; Cerretani, L.; Carrasco-Pancorbo, A.; Gómez Caravaca, A.M.; Segura-Carretero, A.; Fernández-Gutiérrez, A.; Lercker, A. 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–1719. [Google Scholar] [CrossRef]
  14. EFSA Panel on Dietetic Products, Nutrition and Allergies (NDA) Scientific Opinion on the Substantiation of Health Claims Related to Vitamin E and Protection of DNA, Proteins and Lipids from Oxidative Damage (ID 160, 162, 1947), Maintenance of the Normal Function of the Immune System (ID 161, 163), Maintenance of Norm: Vitamin E Related Health Claims. EFSA J. 2010, 8, 1816. [CrossRef]
  15. Ingold, K.U.; Pratt, D.A. Advances in Radical-Trapping Antioxidant Chemistry in the 21st Century: A Kinetics and Mechanisms Perspective. Chem. Rev. 2014, 114, 9022–9046. [Google Scholar] [CrossRef]
  16. Galano, A.; Raúl Alvarez-Idaboy, J. Computational Strategies for Predicting Free Radical Scavengers’ Protection against Oxidative Stress: Where Are We and What Might Follow? Int. J. Quantum Chem. 2019, 119, e25665. [Google Scholar] [CrossRef]
  17. Nenadis, N.; Pyrka, I.; Tsimidou, M.Z. 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. [Google Scholar] [CrossRef]
  18. Amorati, R.; Valgimigli, L. Modulation of the Antioxidant Activity of Phenols by Non-Covalent Interactions. Org. Biomol. Chem. 2012, 10, 4147–4158. [Google Scholar] [CrossRef]
  19. Leopoldini, M.; Russo, N.; Toscano, M. The Molecular Basis of Working Mechanism of Natural Polyphenolic Antioxidants. Food Chem. 2011, 125, 288–306. [Google Scholar] [CrossRef]
  20. Valgimigli, L.; Banks, J.T.; Ingold, K.U.; Lusztyk, J. Kinetic Solvent Effects on Hydroxylic Hydrogen Atom Abstractions Are Independent of the Nature of the Abstracting Radical. Two Extreme Tests Using Vitamin E and Phenol. J. Am. Chem. Soc. 1995, 117, 9966–9971. [Google Scholar] [CrossRef]
  21. Litwinienko, G.; Ingold, K.U. Solvent Effects on the Rates and Mechanisms of Reaction of Phenols with Free Radicals. Acc. Chem. Res. 2007, 40, 222–230. [Google Scholar] [CrossRef]
  22. Foti, M.C. Use and Abuse of the DPPH Radical. J. Agric. Food Chem. 2015, 63, 8765–8776. [Google Scholar] [CrossRef]
  23. Przybylski, P.; Konopko, A.; Łętowski, P.; Jodko-Piórecka, K.; Litwinienko, G. Concentration-Dependent HAT/ET Mechanism of the Reaction of Phenols with 2,2-Diphenyl-1-Picrylhydrazyl (Dpph˙) in Methanol. RSC Adv. 2022, 12, 8131–8136. [Google Scholar] [CrossRef]
  24. Burton, G.W.; Ingold, K.U. Autoxidation of Biological Molecules. 1. Antioxidant Activity of Vitamin E and Related Chain-Breaking Phenolic Antioxidants In Vitro. J. Am. Chem. Soc. 1981, 103, 6472–6477. [Google Scholar] [CrossRef]
  25. Burton, G.W.; Doba, T.; Gabe, E.; Hughes, L.; Lee, F.L.; Prasad, L.; Ingold, K.U. Autoxidation of Biological Molecules. 4. Maximizing the Antioxidant Activity of Phenols. J. Am. Chem. Soc. 1985, 107, 7053–7065. [Google Scholar] [CrossRef]
  26. Hammes-Schiffer, S.; Soudackov, A.V. Proton-Coupled Electron Transfer in Solution, Proteins, and Electrochemistry. J. Phys. Chem. B 2008, 112, 14108–14123. [Google Scholar] [CrossRef]
  27. Foti, M.C.; Daquino, C.; Mackie, I.D.; DiLabio, G.A.; Ingold, K.U. Reaction of Phenols with the 2,2-Diphenyl-1-Picrylhydrazyl Radical. Kinetics and DFT Calculations Applied To Determine ArO-H Bond Dissociation Enthalpies and Reaction Mechanism. J. Org. Chem. 2008, 73, 9270–9282. [Google Scholar] [CrossRef]
  28. Stubbe, J.; van der Donk, W.A. Protein Radicals in Enzyme Catalysis. Chem. Rev. 1998, 98, 705–762. [Google Scholar] [CrossRef]
  29. Mayer, J.M.; Hrovat, D.A.; Thomas, J.L.; Borden, W.T. Proton-Coupled Electron Transfer versus Hydrogen Atom Transfer in Benzyl/Toluene, Methoxyl/Methanol, and Phenoxyl/Phenol Self-Exchange Reactions. J. Am. Chem. Soc. 2002, 124, 11142–11147. [Google Scholar] [CrossRef]
  30. Weinberg, D.R.; Gagliardi, C.J.; Hull, J.F.; Murphy, C.F.; Kent, C.A.; Westlake, B.C.; Paul, A.; Ess, D.H.; McCafferty, D.G.; Meyer, T.J. Proton-Coupled Electron Transfer. Chem. Rev. 2012, 112, 4016–4093. [Google Scholar] [CrossRef]
  31. Darcy, J.W.; Koronkiewicz, B.; Parada, G.A.; Mayer, J.M. A Continuum of Proton-Coupled Electron Transfer Reactivity. Acc. Chem. Res. 2018, 51, 2391–2399. [Google Scholar] [CrossRef]
  32. Hammes-Schiffer, S. Proton-Coupled Electron Transfer: Moving Together and Charging Forward. J. Am. Chem. Soc. 2015, 137, 8860–8871. [Google Scholar] [CrossRef]
  33. Gentry, E.C.; Knowles, R.R. Synthetic Applications of Proton-Coupled Electron Transfer. Acc. Chem. Res. 2016, 49, 1546–1556. [Google Scholar] [CrossRef]
  34. Kaila, V.R.I.; Verkhovsky, M.I.; Wikström, M. Proton-Coupled Electron Transfer in Cytochrome Oxidase. Chem. Rev. 2010, 110, 7062–7081. [Google Scholar] [CrossRef]
  35. Markle, T.F.; Darcy, J.W.; Mayer, J.M. A New Strategy to Efficiently Cleave and Form C–H Bonds Using Proton-Coupled Electron Transfer. Sci. Adv. 2018, 4, eaat5776. [Google Scholar] [CrossRef]
  36. Parada, G.A.; Goldsmith, Z.K.; Kolmar, S.; Pettersson Rimgard, B.; Mercado, B.Q.; Hammarström, L.; Hammes-Schiffer, S.; Mayer, J.M. Concerted Proton-Electron Transfer Reactions in the Marcus Inverted Region. Science 2019, 364, 471–475. [Google Scholar] [CrossRef]
  37. Miller, J.L. Inverted Kinetics Seen in Concerted Charge Transfer. Phys. Today 2019, 72, 16–19. [Google Scholar] [CrossRef]
  38. Truhlar, D.G. Tunneling in Enzymatic and Nonenzymatic Hydrogen Transfer Reactions. J. Phys. Org. Chem. 2010, 23, 660–676. [Google Scholar] [CrossRef]
  39. Layfield, J.P.; Hammes-Schiffer, S. Hydrogen Tunneling in Enzymes and Biomimetic Models. Chem. Rev. 2014, 114, 3466–3494. [Google Scholar] [CrossRef]
  40. Klinman, J.P.; Offenbacher, A.R. Understanding Biological Hydrogen Transfer Through the Lens of Temperature Dependent Kinetic Isotope Effects. Acc. Chem. Res. 2018, 51, 1966–1974. [Google Scholar] [CrossRef]
  41. Karković Marković, A.; Jakobušić Brala, C.; Pilepić, V.; Uršić, S. Kinetic Isotope Effects and Hydrogen Tunnelling in PCET Oxidations of Ascorbate: New Insights into Aqueous Chemistry? Molecules 2020, 25, 1443. [Google Scholar] [CrossRef]
  42. Sajenko, I.; Pilepić, V.; Jakobušić Brala, C.; Uršić, S. Solvent Dependence of the Kinetic Isotope Effect in the Reaction of Ascorbate with the 2,2,6,6-Tetramethylpiperidine-1-Oxyl Radical: Tunnelling in a Small Molecule Reaction. J. Phys. Chem. A 2010, 114, 3423–3430. [Google Scholar] [CrossRef]
  43. Nakanishi, I.; Shoji, Y.; Ohkubo, K.; Fukuzumi, S. Tunneling in the Hydrogen-Transfer Reaction from a Vitamin E Analog to an Inclusion Complex of 2,2-Diphenyl-1-Picrylhydrazyl Radical with β-Cyclodextrin in an Aqueous Buffer Solution at Ambient Temperature. Antioxidants 2021, 10, 1966. [Google Scholar] [CrossRef]
  44. Robles-Almazan, M.; Pulido-Moran, M.; Moreno-Fernandez, J.; Ramirez-Tortosa, C.; Rodriguez-Garcia, C.; Quiles, J.L.; Ramirez-Tortosa, M. Hydroxytyrosol: Bioavailability, Toxicity, and Clinical Applications. Food Res. Int. 2018, 105, 654–667. [Google Scholar] [CrossRef]
  45. Wang, Q.; Wang, C.; Abdullah; Tian, W.; Qiu, Z.; Song, M.; Cao, Y.; Xiao, J. Hydroxytyrosol Alleviates Dextran Sulfate Sodium-Induced Colitis by Modulating Inflammatory Responses, Intestinal Barrier, and Microbiome. J. Agric. Food Chem. 2022, 70, 2241–2252. [Google Scholar] [CrossRef]
  46. Caruso, G.; Torrisi, S.A.; Mogavero, M.P.; Currenti, W.; Castellano, S.; Godos, J.; Ferri, R.; Galvano, F.; Leggio, G.M.; Grosso, G.; et al. Polyphenols and Neuroprotection: Therapeutic Implications for Cognitive Decline. Pharmacol. Therapeut. 2022, 232, 108013. [Google Scholar] [CrossRef]
  47. Paiva-Martins, F.; Silva, A.; Almeida, V.; Carvalheira, M.; Serra, C.; Rodrígues-Borges, J.; Fernandes, J.; Belo, L.; Santos-Silva, A. Protective Activity of Hydroxytyrosol Metabolites on Erythrocyte Oxidative-Induced Hemolysis. J. Agric. Food Chem. 2013, 61, 6636–6642. [Google Scholar] [CrossRef]
  48. Maiuolo, J.; Costanzo, P.; Masullo, M.; D’Errico, A.; Nasso, R.; Bonacci, S.; Mollace, V.; Oliverio, M.; Arcone, R. Hydroxytyrosol–Donepezil Hybrids Play a Protective Role in an In Vitro Induced Alzheimer’s Disease Model and in Neuronal Differentiated Human SH-SY5Y Neuroblastoma Cells. Int. J. Mol. Sci. 2023, 24, 13461. [Google Scholar] [CrossRef]
  49. Marcucci, G.; Domazetovic, V.; Nediani, C.; Ruzzolini, J.; Favre, C.; Brandi, M.L. Oxidative Stress and Natural Antioxidants in Osteoporosis: Novel Preventive and Therapeutic Approaches. Antioxidants 2023, 12, 373. [Google Scholar] [CrossRef]
  50. Francisco, V.; Ruiz-Fernandez, C.; Lahera, V.; Lago, F.; Pino, J.; Skaltsounis, L.; Angel Gonzalez-Gay, M.; Mobasheri, A.; Gomez, R.; Scotece, M.; et al. Natural Molecules for Healthy Lifestyles: Oleocanthal from Extra Virgin Olive Oil. J. Agric. Food Chem. 2019, 67, 3845–3853. [Google Scholar] [CrossRef]
  51. Filardo, S.; Roberto, M.; Di Risola, D.; Mosca, L.; Di Pietro, M.; Sessa, R. Olea Europaea L-Derived Secoiridoids: Beneficial Health Effects and Potential Therapeutic Approaches. Pharmacol. Therapeut. 2024, 254, 108595. [Google Scholar] [CrossRef]
  52. Ganguly, R.; Singh, S.V.; Jaiswal, K.; Kumar, R.; Pandey, A.K. Modulatory Effect of Caffeic Acid in Alleviating Diabetes and Associated Complications. World J. Diab. 2023, 14, 62–75. [Google Scholar] [CrossRef]
  53. Rahman, M.; Rahaman, S.; Islam, R.; Rahman, F.; Mithi, F.M.; Alqahtani, T.; Almikhlafi, M.A.; Alghamdi, S.Q.; Alruwaili, A.S.; Hossain, S.; et al. Role of Phenolic Compounds in Human Disease: Current Knowledge and Future Prospects. Molecules 2022, 27, 233. [Google Scholar] [CrossRef]
  54. Guo, C.; Huang, Q.; Wang, Y.; Yao, Y.; Li, J.; Chen, J.; Wu, M.; Zhang, Z.; Mingyao, E.; Qi, H.; et al. Therapeutic Application of Natural Products: NAD<SUP>+</SUP> Metabolism as Potential Target. Phytomedicine 2023, 114, 154768. [Google Scholar] [CrossRef]
  55. Kaur, J.; Gulati, M.; Singh, S.K.; Kuppusamy, G.; Kapoor, B.; Mishra, V.; Gupta, S.; Arshad, M.F.; Porwal, O.; Jha, N.K.; et al. Discovering Multifaceted Role of Vanillic Acid beyond Flavours: Nutraceutical and Therapeutic Potential. Trends Food Sci. Tech. 2022, 122, 187–200. [Google Scholar] [CrossRef]
  56. Foti, M.; Ingold, K.U.; Lusztyk, J. The Surprisingly High Reactivity of Phenoxyl Radicals. J. Am. Chem. Soc. 1994, 116, 9440–9447. [Google Scholar] [CrossRef]
  57. Foti, M.C.; Daquino, C.; Geraci, C. Electron-Transfer Reaction of Cinnamic Acids and Their Methyl Esters with the DPPH Radical in Alcoholic Solutions. J. Org. Chem. 2004, 69, 2309–2314. [Google Scholar] [CrossRef]
  58. Semidalas, C.; Semidalas, E.; Matsoukas, M.T.; Nixarlidis, C.; Zoumpoulakis, P. In Silico Studies Reveal the Mechanisms behind the Antioxidant and Anti-Inflammatory Activities of Hydroxytyrosol. Med. Chem. Res. 2016, 25, 2498–2511. [Google Scholar] [CrossRef]
  59. Mayer, J.M. Proton-Coupled Electron Transfer: A Reaction Chemist’s View. Ann. Rev. Phys. Chem. 2004, 55, 363–390. [Google Scholar] [CrossRef]
  60. Rossini, E.; Bochevarov, A.D.; Knapp, E.W. Empirical Conversion of pKa Values between Different Solvents and Interpretation of the Parameters: Application to Water, Acetonitrile, Dimethyl Sulfoxide, and Methanol. ACS Omega 2018, 3, 1653–1662. [Google Scholar] [CrossRef]
  61. Bell, R.P. The Tunnel Effect in Chemistry; Springer US: Boston, MA, USA, 1980; ISBN 978-0-412-21340-3. [Google Scholar]
  62. Jakobušić Brala, C.; Pilepić, V.; Sajenko, I.; Karković, A.; Uršić, S. Ions Can Move a Proton-Coupled Electron-Transfer Reaction into Tunneling Regime. Helv. Chim. Acta 2011, 94, 1718–1731. [Google Scholar] [CrossRef]
  63. Kohen, A. Kinetic Isotope Effects as Probes for Hydrogen Tunneling in Enzyme Catalysis. In Isotope Effects in Chemistry and Biology; Kohen, A., Limbach, H.H., Eds.; Taylor & Francis: London, UK; CRC Press: New York, NY, USA, 2006; pp. 744–764. [Google Scholar]
  64. Fukui, K. The Path of Chemical Reactions—The IRC Approach. Acc. Chem. Res. 1981, 14, 363–368. [Google Scholar] [CrossRef]
  65. Knizia, G.; Klein, J.E.M.N. Electron Flow in Reaction Mechanisms—Revealed from First Principles. Ang. Chem. Int. Ed. 2015, 54, 5518–5522. [Google Scholar] [CrossRef]
  66. Knizia, G. Intrinsic Atomic Orbitals: An Unbiased Bridge between Quantum Theory and Chemical Concepts. J. Chem. Theory Comput. 2013, 9, 4834–4843. [Google Scholar] [CrossRef]
  67. Klein, J.E.M.N.; Knizia, G. cPCET versus HAT: A Direct Theoretical Method for Distinguishing X–H Bond-Activation Mechanisms. Angew. Chem. Int. Ed. 2018, 57, 11913–11917. [Google Scholar] [CrossRef]
  68. Andersson, M.P.; Uvdal, P. New Scale Factors for Harmonic Vibrational Frequencies Using the B3LYP Density Functional Method with the Triple-ζ Basis Set 6-311+G(d,p). J. Phys. Chem. A 2005, 109, 2937–2941. [Google Scholar] [CrossRef]
  69. Cheng, Y.-Y.; Cui, C.-X. Theoretical Study on Hydrogen Transfer in the Dissociation of Dimethyl Disulfide Radical Cations. Phys. Chem. Chem. Phys. 2023, 25, 3780–3788. [Google Scholar] [CrossRef]
  70. Hess, K.M.; Leach, I.F.; Wijtenhorst, L.; Lee, H.; Klein, J.E.M.N. Valence Tautomerism Induced Proton Coupled Electron Transfer:X−H Bond Oxidation with a Dinuclear Au(II) Hydroxide Complex. Angew. Chem. Int. Ed. 2024, 63, e202318916. [Google Scholar] [CrossRef]
  71. Sjoberg, P.; Murray, J.S.; Brinck, T.; Politzer, P. Average Local Ionization Energies on the Molecular Surfaces of Aromatic Systems as Guides to Chemical Reactivity. Can. J. Chem. 1990, 68, 1440–1443. [Google Scholar] [CrossRef]
  72. Parr, R.; Yang, W. Density Functional-Approach to the Frontier-Electron Theory of Chemical-Reactivity. J. Am. Chem. Soc. 1984, 106, 4049–4050. [Google Scholar] [CrossRef]
  73. Politzer, P.; Murray, J.S.; Bulat, F.A. Average Local Ionization Energy: A Review. J. Mol. Model. 2010, 16, 1731–1742. [Google Scholar] [CrossRef]
  74. Fuentealba, P.; Cárdenas, C. Chapter 14—On the Analysis of the Fukui Function. In Chemical Reactivity; Kaya, S., von Szentpály, L., Serdaroğlu, G., Guo, L., Eds.; Elsevier: Amsterdam, The Netherlands, 2023; pp. 421–432. ISBN 978-0-323-90257-1. [Google Scholar]
  75. Chattaraj, P.K. Chemical Reactivity Theory: A Density Functional View; CRC Press/Taylor & Francis: Boca Raton, FL, USA, 2009; ISBN 978-1-4200-6543-5. [Google Scholar]
  76. Jakobušić Brala, C.; Fabijanić, I.; Marković, A.K.; Pilepić, V. The Average Local Ionization Energy and Fukui Function of L-Ascorbate, the Local Reactivity Descriptors of Antioxidant Reactivity. Comp. Theor. Chem. 2014, 1049, 1–6. [Google Scholar] [CrossRef]
  77. Fabijanić, I.; Brala, C.J.; Pilepić, V. The DFT Local Reactivity Descriptors of α-Tocopherol. J. Mol. Model. 2015, 21, 99. [Google Scholar] [CrossRef]
  78. Pilepić, V.; Uršić, S. Nucleophilic Reactivity of the Nitroso Group. Fukui Function DFT Calculations for Nitrosobenzene and 2-Methyl-2-Nitrosopropane. Theochem.-J. Mol. Struct. 2001, 538, 41–49. [Google Scholar] [CrossRef]
  79. Frisch, M.J.; Trucks, G.W.; Schlegel, H.B.; Scuseria, G.E.; Robb, M.A.; Cheeseman, J.R.; Scalmani, G.; Barone, V.; Petersson, G.A.; Nakatsuji, H.; et al. Gaussian 16, Revision C.01; Gaussian, Inc.: Wallingford, CT, USA, 2016. [Google Scholar]
  80. Tomasi, J.; Mennucci, B.; Cammi, R. Quantum Mechanical Continuum Solvation Models. Chem. Rev. 2005, 105, 2999–3094. [Google Scholar] [CrossRef]
  81. Knizia, G.; IboView. Available online: http://www.iboview.org (accessed on 25 March 2024).
  82. Lu, T.; Chen, F. Multiwfn: A Multifunctional Wavefunction Analyzer. J. Comput. Chem. 2012, 33, 580–592. [Google Scholar] [CrossRef]
  83. Wigner, E. On the Quantum Correction For Thermodynamic Equilibrium. Phys. Rev. 1932, 40, 749–759. [Google Scholar] [CrossRef]
  84. Ghysels, A.; Verstraelen, T.; Hemelsoet, K.; Waroquier, M.; Van Speybroeck, V. TAMkin: A Versatile Package for Vibrational Analysis and Chemical Kinetics. J. Chem. Inf. Model. 2010, 50, 1736–1750, TAMkin version 1.2.6. Available online: https://molmod.github.io/tamkin/index.html (accessed on 3 June 2024).
Figure 1. Structures of olive oil phenols.
Figure 1. Structures of olive oil phenols.
Ijms 25 06341 g001
Figure 2. Transition structure for the proposed hydrogen transfer from HVA to DPPH• obtained from the DFT calculations. The normal mode displacement vector shown (blue arrow) for the unique imaginary frequency of 1608.7i cm−1 is associated with a motion of the H atom of the O-H moiety of HVA to the N atom of DPPH•.
Figure 2. Transition structure for the proposed hydrogen transfer from HVA to DPPH• obtained from the DFT calculations. The normal mode displacement vector shown (blue arrow) for the unique imaginary frequency of 1608.7i cm−1 is associated with a motion of the H atom of the O-H moiety of HVA to the N atom of DPPH•.
Ijms 25 06341 g002
Figure 3. Changes in energy plotted along the IRC path (inset), α (green), and β (purple) intrinsic bond orbitals (IBO) involved in hydrogen transfer from HVA to DPPH• plotted along the intrinsic reaction coordinate (IRC) obtained from the DFT calculations. For IBO analysis along the entire IRC, see also the corresponding video files in the Supplementary Materials.
Figure 3. Changes in energy plotted along the IRC path (inset), α (green), and β (purple) intrinsic bond orbitals (IBO) involved in hydrogen transfer from HVA to DPPH• plotted along the intrinsic reaction coordinate (IRC) obtained from the DFT calculations. For IBO analysis along the entire IRC, see also the corresponding video files in the Supplementary Materials.
Ijms 25 06341 g003aIjms 25 06341 g003b
Figure 4. The average local ionization energy Ī(r) (a) and the electron donor Fukui function f ¯(r) (b) plotted on 0.001 a.u. electron density isosurface of HVA. The color scales correspond to the range of values 13.6 eV (red) to 4.5 eV (blue) for Ī(r) and 0.0 a.u. (red) to 0.0003 a.u. (blue) for f ¯(r). The lowest values of Ī(r) (blue areas) and the highest values of f ¯(r) (blue areas) indicate the most preferred sites for donating electrons. The lowest average local ionization energy ĪS,min has a value of 8.64 eV.
Figure 4. The average local ionization energy Ī(r) (a) and the electron donor Fukui function f ¯(r) (b) plotted on 0.001 a.u. electron density isosurface of HVA. The color scales correspond to the range of values 13.6 eV (red) to 4.5 eV (blue) for Ī(r) and 0.0 a.u. (red) to 0.0003 a.u. (blue) for f ¯(r). The lowest values of Ī(r) (blue areas) and the highest values of f ¯(r) (blue areas) indicate the most preferred sites for donating electrons. The lowest average local ionization energy ĪS,min has a value of 8.64 eV.
Ijms 25 06341 g004
Table 1. Rate constants (kArOH) and kinetic isotope effects (KIEs) for the reaction of phenol and DPPH• in 1,4-dioxane:water solvent mixtures (0.99:0.01 v/v or 0.95:0.05 v/v) at 25 °C.
Table 1. Rate constants (kArOH) and kinetic isotope effects (KIEs) for the reaction of phenol and DPPH• in 1,4-dioxane:water solvent mixtures (0.99:0.01 v/v or 0.95:0.05 v/v) at 25 °C.
PhenolkArOH/M−1 s−1KIE
Simple phenols
Hydroxytyrosol2.56 (0.07)16.0 (1.0)
Tyrosol0.0083 (0.0007)3.3 (0.4)
Homovanillyl alcohol0.150 (0.002)15.4 (0.3)
Secoiridoids
Oleuropein3.03 (0.10)3.2 (0.4)
Oleacein2.83 (0.36)1.7 (0.8)
Oleocanthal0.88 (0.05)1.8 (0.2)
Phenolic acids
Caffeic acid0.537 (0.009)16.7 (0.5)
p-coumaric acid 10.0025 (0.0001)3.4 (0.02)
Ferulic acid 20.161 (0.029)
Vanillic acid 20.0058 (0.0001)4.1 (0.1)
Vanillin0.00291 (0.00001)
1 0.2 M p-coumaric acid; 2 0.1 M phenol.
Table 2. Thermodynamic activation parameters for the reaction of phenol and DPPH• in 1,4-dioxane-H2O solvent mixtures (0.99:0.01 v/v or 0.95:0.05 (v/v)) at 25 °C.
Table 2. Thermodynamic activation parameters for the reaction of phenol and DPPH• in 1,4-dioxane-H2O solvent mixtures (0.99:0.01 v/v or 0.95:0.05 (v/v)) at 25 °C.
PhenolΔG/kJ mol−1ΔH/kJ mol−1ΔS/J K−1 mol−1Ea/kJ mol−1ln(A/M−1 s−1)
HOTyr69.0 (0.8)40.0 (0.7)−146.0 (2.2)45.2 (0.3)16.0 (0.1)
HVA82.5 (0.9)33.1 (1.2)−143.8 (4.2)35.6 (1.2)13.1 (0.5)
caffeic acid74.7 (1.2)40.8 (0.9)−113.4 (2.9)43.4 (0.9)16.8 (0.3)
Table 3. Isotopic difference in activation parameters for the reaction of phenol and DPPH• in 1,4-dioxane-H2O solvent mixtures (0.99:0.01 v/v or 0.95:0.05 v/v) at 25 °C.
Table 3. Isotopic difference in activation parameters for the reaction of phenol and DPPH• in 1,4-dioxane-H2O solvent mixtures (0.99:0.01 v/v or 0.95:0.05 v/v) at 25 °C.
PhenolΔΔG/kJ mol−1ΔEa (D, H)/kJ mol−1AH/AD
HOTyr6.9 (3.1)8.1 (3.0)0.6
HVA6.5 (1.8)5.9 (1.4)1.3
caffeic acid6.8 (1.8)9.5 (1.3)0.3
Table 4. Selected properties of the transition state structures obtained for hydrogen transfer from the 4-OH moiety of the phenols to the DPPH• N atom obtained from the DFT calculations.
Table 4. Selected properties of the transition state structures obtained for hydrogen transfer from the 4-OH moiety of the phenols to the DPPH• N atom obtained from the DFT calculations.
PhenolN⋯H/Å aH⋯O/Å bN⋯H⋯O/° cvH/cm−1 dvD/cm−1 d
HOTyr1.3021.177163.31646.6i1218.8i
HVA1.3131.167162.51608.7i1192.4i
caffeic acid1.2751.202162.51689.0i1246.9i
a Distance between the N atom of DPPH• and the H atom. b Distance between the O atom of the phenol and the H atom. c N⋯H⋯O angle. d Normal mode imaginary frequency associated with the proton or deuterium transfer scaled by a factor of 0.9679 [68].
Table 5. Activation and reaction Gibbs energy and KIE for hydrogen transfer from the 4-OH moiety of the phenols to the DPPH• N atom obtained from the DFT calculations.
Table 5. Activation and reaction Gibbs energy and KIE for hydrogen transfer from the 4-OH moiety of the phenols to the DPPH• N atom obtained from the DFT calculations.
PhenolΔGH/kJ mol−1 ΔGD/kJ mol−1 ΔrG/kJ mol−1KIESC aKIETUNN b
HOTyr+79.0+83.4−1.36.018.9
HVA+81.5+85.610.35.219.2
caffeic acid+77.5+83.4−2.64.919.6
a Kinetic isotope effects obtained from the calculated ΔG for H and D substituted reactants. b Kinetic isotope effects corrected for H atom tunneling using the Wigner method.
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Torić, J.; Karković Marković, A.; Mustać, S.; Pulitika, A.; Jakobušić Brala, C.; Pilepić, V. Proton-Coupled Electron Transfer and Hydrogen Tunneling in Olive Oil Phenol Reactions. Int. J. Mol. Sci. 2024, 25, 6341. https://doi.org/10.3390/ijms25126341

AMA Style

Torić J, Karković Marković A, Mustać S, Pulitika A, Jakobušić Brala C, Pilepić V. Proton-Coupled Electron Transfer and Hydrogen Tunneling in Olive Oil Phenol Reactions. International Journal of Molecular Sciences. 2024; 25(12):6341. https://doi.org/10.3390/ijms25126341

Chicago/Turabian Style

Torić, Jelena, Ana Karković Marković, Stipe Mustać, Anamarija Pulitika, Cvijeta Jakobušić Brala, and Viktor Pilepić. 2024. "Proton-Coupled Electron Transfer and Hydrogen Tunneling in Olive Oil Phenol Reactions" International Journal of Molecular Sciences 25, no. 12: 6341. https://doi.org/10.3390/ijms25126341

APA Style

Torić, J., Karković Marković, A., Mustać, S., Pulitika, A., Jakobušić Brala, C., & Pilepić, V. (2024). Proton-Coupled Electron Transfer and Hydrogen Tunneling in Olive Oil Phenol Reactions. International Journal of Molecular Sciences, 25(12), 6341. https://doi.org/10.3390/ijms25126341

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop