Structure-Thermodynamics-Antioxidant Activity Relationships of Selected Natural Phenolic Acids and Derivatives: An Experimental and Theoretical Evaluation

Yuzhen Chen; Huizhi Xiao; Jie Zheng; Guizhao Liang

doi:10.1371/journal.pone.0121276

Abstract

Phenolic acids and derivatives have potential biological functions, however, little is known about the structure-activity relationships and the underlying action mechanisms of these phenolic acids to date. Herein we investigate the structure-thermodynamics-antioxidant relationships of 20 natural phenolic acids and derivatives using DPPH^• scavenging assay, density functional theory calculations at the B3LYP/6-311++G(d,p) levels of theory, and quantitative structure-activity relationship (QSAR) modeling. Three main working mechanisms (HAT, SETPT and SPLET) are explored in four micro-environments (gas-phase, benzene, water and ethanol). Computed thermodynamics parameters (BDE, IP, PDE, PA and ETE) are compared with the experimental radical scavenging activities against DPPH^•. Available theoretical and experimental investigations have demonstrated that the extended delocalization and intra-molecular hydrogen bonds are the two main contributions to the stability of the radicals. The C = O or C = C in COOH, COOR, C = CCOOH and C = CCOOR groups, and orthodiphenolic functionalities are shown to favorably stabilize the specific radical species to enhance the radical scavenging activities, while the presence of the single OH in the ortho position of the COOH group disfavors the activities. HAT is the thermodynamically preferred mechanism in the gas phase and benzene, whereas SPLET in water and ethanol. Furthermore, our QSAR models robustly represent the structure-activity relationships of these explored compounds in polar media.

Citation: Chen Y, Xiao H, Zheng J, Liang G (2015) Structure-Thermodynamics-Antioxidant Activity Relationships of Selected Natural Phenolic Acids and Derivatives: An Experimental and Theoretical Evaluation. PLoS ONE 10(3): e0121276. https://doi.org/10.1371/journal.pone.0121276

Academic Editor: Dennis Salahub, University of Calgary, CANADA

Received: September 10, 2014; Accepted: January 29, 2015; Published: March 24, 2015

This is an open access article, free of all copyright, and may be freely reproduced, distributed, transmitted, modified, built upon, or otherwise used by anyone for any lawful purpose. The work is made available under the Creative Commons CC0 public domain dedication

Data Availability: All relevant data are within the paper and its Supporting Information files.

Funding: G.L. is thankful for financial support from the National Natural Science Foundation of China (No.10901169), the Natural Science Foundation Project of Chongqing CSTC (No.cstc2012gg-gjhz10003), and the Fundamental Research Funds for the Central Universities (No. CQDXWL-2014-Z009). J.Z. thanks the National Science Foundation of USA (CAREER Award No. CBET-0952624 and CBET-1158447). The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.

Competing interests: Prof. Jie Zheng, as one of corresponding authors in this manuscript and a PLOS ONE Academic Editor, declares no conflict of interest. This does not alter the authors’ adherence to all the PLOS ONE policies on sharing data and materials.

Introduction

Phenolic acids and derivatives are a subclass of a larger category of metabolites (commonly named as “phenolics”), which widely spread throughout the plant kingdom [1]. Although the basic skeleton containing a carboxylic acid function and a phenolic ring remains the same [1,2], the numbers and positions of hydroxyl groups and other substituents on the aromatic ring create significant variations in both structures and functions. Recent interests in phenolic acids and derivatives mainly come from their potential biological functions, such as anti-inflammatory, antiallergic, antimicrobial, anticarcinogenic and antiviral activities [2]. However, little is known about the detailed action mechanisms of this class of compounds [1,3]. Thus, a thorough investigation of approximately 8000 naturally occurring phenolics [2] will help us to identify the compounds with desirable functionalities important for public health [4].

Phenolic compounds generally exert their protective activities by three different mechanisms [5–9]: hydrogen atom transfer (HAT) (Eq 1), electron transfer-proton transfer (SETPT) (Eq 2 and 3), and sequential proton loss-electron transfer (SPLET) (Eq 4, 5 and 6). All the mechanisms are believed to play important roles in determining radical scavenging activities of antioxidants in various environmental conditions [5,10]. It has been shown that radical scavenging activities of phenolic antioxidants are related to the phenolic O-H bond dissociation enthalpy (BDE), ionization potential (IP), proton dissociation enthalpy (PDE), proton affinity (PA) and electron transfer enthalpy (ETE) [6,10–14].

(1)

(2)

(3)

(4)

(5)

(6)

Radical scavenging activities of phenolic antioxidants are largely influenced by their structural and environmental features in vivo [5]. In vitro determination of radical scavenging activities has often relied on the ability of phenolic compounds to quench the color of stable radicals such as the 2,2-diphenyl-1-picrylhydrazyl (DPPH) radical [15,16]. This reaction with DPPH is considerably affected by experimental conditions including the absolute and relative concentrations of DPPH and antioxidants, solvents, hydrogen bonding strength, room temperature, time, and unspecified pH [16]. Accordingly, the “one-size-fits-all” radical scavenging mechanisms may or may not work, or a combination of several mechanisms may provide a better explanation for radical scavenging activities of phenolic antioxidants [16]. Besides, lack of standardization in sample preparation, reaction conditions, analytical protocols and expression of antioxidant action leads to difficulty or impossibility to compare different results obtained from different laboratories.

The quantum chemistry and computation methodologies allow obtaining atomic-level structures and energetic information of the systems with accuracy equivalent to or greater than those obtained from experiments. Therefore, theoretical calculations have been widely used as a cogent tool for rational design of novel potential drugs and for investigation of the underlying structure-activity relationships of these drugs [12]. There are several successful examples of rational interpretation of structure-activity relationships of some natural antioxidants [11–13,17,18] and design of novel antioxidants [6,19,20] using powerful and economical quantum chemical methods especially density functional theory (DFT). It has been reported that since most of the DFT methods underestimate thermodynamics parameters such as BDEs, these methods are more reliable and suitable for relative calculations than for absolute calculations [21].

In parallel, quantitative structure-activity relationship (QSAR) is one of powerful computational methods for prediction of activities, determination of action mechanisms, design of drugs, materials, catalysts, and proteins/peptides with desirable activities and functions [22–28]. The QSAR methods have been used to explore the structure-activity relationships and action mechanisms of phenolic compounds [29]. A predictive and interpretable QSAR model can help to further understand mechanisms of action of the explored molecules toward the target systems. Since chemical descriptors are the core of QSAR modeling, significant efforts and progress have been made to develop a wide variety of chemical descriptors to describe different levels of chemical, physical, and structural characteristics of the target molecules/systems [27]. The comprehensibly physiochemical descriptors, i.e., BDE, IP, PDE, PA and ETE, calculated by DFT calculations have been employed to elucidate the structure-activity relationships of the investigated phenolic antioxidants [6,10,13,30]. Correlation methods are the second important part of QSAR modeling. A linear model has some advantages such as simplicity and interpretability over a non-interpretably nonlinear one [27], ensuring the interpretability of QSAR models.

In this work, we investigated the structure-thermodynamics-antioxidant activity relationships of 20 natural phenolic acids and derivatives using combined experimental and computational approaches. The antioxidant activities of the compounds were evaluated for their total antioxidant capacity using the DPPH^• scavenging assay. With high-precision DFT calculations, five sets of thermodynamics parameters (BDE, IP, PDE, PA and ETE) were identified and used to evaluate three working mechanisms (HAT, SETPT and SPLET) under different micro-environments (gas-phase, benzene, water and ethanol). Moreover, the spin density in free radicals, as well as the highest occupied molecular orbital (HOMO) distribution, was also computed to better describe the radical scavenging reactivity of the studied molecules. We also established the QSAR models to characterize the structure-activity relationships of the antioxidants. This work provides structural-based insights into the action mechanisms of phenolics, which may help expand their applications in pharmaceutical and food science.

Principles and Methods

Chemicals

DPPH was purchased from Sigma Chemical Co. A total of 20 naturally occurring phenolic acids and derivatives (Table A in S1 File) were purchased from Aladdin Industrial Inc. These compounds contain 16 natural phenolic acids, with two distinguishing constitutive carbon frameworks: COOH and C = CCOOH, and 4 compounds referred to as phenolic acid derivatives, with a COOR structure, including propyl gallate (01), methyl gallate (09), methyl vanillate (15), and ferulic acid ethyl ester (18).

DPPH radical scavenging assay

The DPPH assay has been widely used for the measurement of free radical scavenging capacity of various natural products [15,16]. The DPPH radical is a stable organic free radical with adsorption band at 515–528 nm. It loses this adsorption when accepting an electron or a free radical species, which results in a visually noticeable discoloration from purple to yellow [31]. The DPPH radical scavenging assay was measured using the method of Sun and Ho [15]. Briefly, 2 mL DPPH solution (0.2 mmol/L, in ethanol) was incubated with different concentrations of the sample. The reaction mixture was shaken and incubated in the dark for 30 min at room temperature. The absorbance was read at 517 nm against ethanol. The control containing ethanol instead of the sample and the blank containing ethanol instead of DPPH solution were also made. The test was run in triplicate and the inhibition of the DPPH radical of the sample was calculated according to the following formula:

The percentage of DPPH radical scavenging activity was plotted against the sample concentration to acquire the IC₅₀ value, defined as the concentration of sample necessary to cause 50% inhibition. The DPPH measurements were assayed using a spectrophotometer (UNICO7200, Unico(Shanghai) Instrument Co., Ltd.).

DFT calculations

All geometries were completely optimized in all internal degrees of freedom using DFT calculations with the B3LYP [32] functional and 6-311++G(d,p) [33] basis set. For all of the radical systems, the unrestricted B3LYP/6-311++G(d,p) method was used. Vibrational frequencies were further computed at the same level to ensure no imaginary frequency for the optimized structures. For open-shell species, accuracy of the energy evaluation is sensitive to spin contamination. Here, spin contaminations of radicals were found in the 0.76–0.78 range, and then spin contaminants dropped to a correct value of 0.75 after the annihilation of the first spin contaminant.

Here, BDE, IP, PDE, PA and ETE were determined in the gas phase, benzene, water and ethanol solvents at 298 K based on the following expressions (Eq 7, 8, 9, 10 and 11):

(7)

(8)

(9)

(10)

(11)

Therein, the BDE values were used to estimate the reactivity of an ArOH in HAT. The IP and PDE values from the ArOH^•+ radical cation were calculated to describe the SETPT mechanism. The PA values of the phenoxide anion, ArO⁻, were used to characterize the reaction enthalpy of the first step, and ETEs for the reaction enthalpy of the following step, in the SPLET mechanism.

The solvent effects were computed using an integral equation formalism polarized continuum model IEF-PCM method [34]. The calculated gas-phase enthalpy of a proton and an electron is 1.483 and 0.752 kcal/mol, respectively [35]. Proton and electron solvation enthalpies were taken from Rimarcik et al.’s report [36], and hydrogen atom solvation enthalpies from Parker’s [37] and Bizarro et al’s report [38]. All calculations were performed using Gaussian 09 [39].

QSAR modeling and validation

Simple or multiple linear regression analysis (SLR or MLR) was performed to derive the QSAR models. The BDE, IP, PDE, PA and ETE descriptors were regarded as inputs of the models, and the negative logarithmic IC₅₀ values (pIC₅₀) were treated as dependent variables. The BDE-antioxidant activity relationships were described by SLR-based QSAR models, while the IP-PDE-antioxidant activity and the PA-ETE-antioxidant activity relationships were described by MLR-based models. The leave-one-out [40] cross validation was employed to test the predictive ability of the models. The modeling performance was assessed based on a serial of statistical parameters, including multiple correlation coefficients (R² for regression modeling and Q² for leave-one-out cross validation), standard error, and Fisher’s criterion, etc. The linear regression operation and the leave-one-out cross validation was performed by our in-house applied program.

Results and Discussion

Effects of intra-molecular hydrogen bonds on the stability of the radicals

Considering that hydrogen bonds play an important role in the stability of the studied parent phenolic compounds and their corresponding hydrogen atom-abstracted radicals, we comparably calculated hydrogen bonds between the phenol-O^• and the meta OH of both the radicals and their parent molecules. The energies of the optimized 7 parent molecules (Fig. A in S1 File) with hydrogen bond(s) between the phenol-OH and the meta OH are reduced by 3.1–8.7 kcal/mol relative to those of 7 parent molecules without hydrogen bond(s) between the phenol-OH and the meta OH. The corresponding radicals derived from 7 parent compounds display that a 4-O^•···H-O-3 hydrogen bond is formed between the 4-O^• center and its neighboring OH group (Fig. 1), demonstrating an important contribution to the stability of the radicals. The energy difference (Fig. 2) of the radicals optimized at the B3LYP/6-311++G(d,p) level shows that among 4 different reaction conditions, the relative energies of the 4-O^• centers forming a hydrogen bond with the meta OH are less than those of the 4-O^• centers without any hydrogen bond. This supported hydrogen bond(s) could help to stabilize these radicals and be beneficial to the radical scavenging capacity of the parent molecules. Thus, in comparison to vanillic acid (02), an O-H···O hydrogen bond between the O^• center and the meta OH in protocatechuic acid (10) was formed, thereby to stabilize the corresponding phenoxy radical, which could partly explain why the former has a low radical scavenging activity (pIC₅₀ = 2.44) relative to the latter (pIC₅₀ = 4.25). More importantly, the largest energy differences between the radicals with hydrogen bonds and without hydrogen bonds derived from one same compound were found in the gas phase, whereas the smallest ones in water (Fig. 2). We speculated that water molecules are likely to form strong hydrogen bonds with the corresponding 4-O^•, thereby to conductively stabilize the corresponding radicals.

Download:

Fig 1. Optimized structures of 7 compounds with orthodiphenolic functionalities calculated at the B3LYP/6-311++G(d,p) level in ethanol.

(left: a hydrogen bond is formed between the 4-O^• center and the meta OH. right: no hydrogen bond is formed between the 4-O^• center and the meta OH.)

https://doi.org/10.1371/journal.pone.0121276.g001

Download:

Fig 2. The energy difference (ΔE in kcal/mol) caused by the hydrogen bond between the O^• center and the meta OH for 7 radicals calculated at the B3LYP/6-311++G(d,p) levels of theory in 4 reaction media.

https://doi.org/10.1371/journal.pone.0121276.g002

Structure-BDE-activity relationships

We calculated a total of 31 O-H BDE values generated from 20 parent phenolic acids and derivatives with intra-molecular hydrogen bonds at the B3LYP/6-311++G(d,p) level of theory. As displayed in Table 1, the O-H BDE values exhibit an analogous size for each compound whatever the reaction medium. This indicated that from the thermodynamic point of view, the HAT mechanism has a close occurrence probability in 4 reaction media.

Download:

Table 1. BDEs of 20 investigated phenolic compounds and phenol calculated at the B3LYP/6-311++G(d,p) level ^{^a}.

https://doi.org/10.1371/journal.pone.0121276.t001

We investigated the conformational, electronic, and geometrical features of phenolic systems, which are of crucial importance for understanding their structure-antioxidant activity relationships [10]. For three compounds containing 3 adjacent phenolic OHs, such as propyl gallate (01), gallic acid (03) and methyl gallate (09), the 4-OH BDE values are generally smaller than the corresponding 3-OH and 5-OH BDE values (Table 1). Moreover, the energies of the 4-O^• radicals are smaller than those of the 3-O^• and 5-O^• radicals (Table 2), leading to the fact that for the same compound, the O^• centers in the para position to 3 groups (-COOCH₂CH₂CH₃, -COOH and COOCH₃) are more stable. For the spin density distribution of propyl gallate (01) displayed in Fig. 3, the spin density (0.325) in the 4-O^• center is lower than those of the 3-O^• (0.373) and 5-O^• (0.358) centers. Similarly, gallic acid (03) has a lower spin density (0.326) in the 4-O^• center than those in the 3-O^• and 5-O^• centers (0.360 and 0.380), while methyl gallate (09) has a lower spin density of 0.326 in the 4-O^• center than those in the 3-O^• and 5-O^• centers (0.375 and 0.357). As can be seen, in the same molecule, a low spin density of the O^• center indicates that the spin density experiences a more extended electronic delocalization, which is conducive to the stability of the radicals, thus leading to the smaller the corresponding BDE value.

Download:

Fig 3. Spin density values of phenoxy radicals of 20 investigated phenolic compounds and phenol calculated at the B3LYP/6-311++G(d,p) levels of theory in ethanol.

https://doi.org/10.1371/journal.pone.0121276.g003

Download:

Table 2. Absolute (E in kcal/mol) and relative energies (ΔE in kcal/mol) of various hydrogen atom-abstracted radicals in 4 reaction environments calculated at the B3LYP/6-311++G(d, p) levels of theory^{^a}.

https://doi.org/10.1371/journal.pone.0121276.t002

The distinct extended delocalization affects the spin density distribution of the corresponding radicals, thereby to bring out different BDEs and subsequent radical-scavenging activities for these investigated phenolic compounds. The extended delocalization in these radical systems is composed by the radical center, whole aromatic ring, relevant phenolic hydroxyl(s), and (C = C)COOH and COOR group(s). We found that the C = O spin density in 4-O^• distinctly differed from that of 3-O^• and 5-O^•. For example, propyl gallate (01) has a C = O spin density of (-0.038, 0.049) for 4-O^•, (0.034, 0.004) for 3-O^•, and (0.028, -0.010) for 5-O^•, respectively. Similarly, the C = O spin densities of 4-O^•, 3-O^• and 5-O^• are (-0.022, 0.053), (0.014, 0.006) and (0.022, 0.003) in Gallic acid (03) while (-0.036, 0.047), (0.025, -0.011) and (0.017, 0.006) in methyl gallate (09), respectively. Overall, the C = O in the carboxylic acid could greatly strengthen the para 4-O^• center to form a more electronic delocalization relative to the meta 3-O^• and 5-O^• center. The distinctly extended C = O delocalization can be seen in the 4-O^• centers of other 10 compounds (Fig. 3), including vanillic acid (02, 4-O^•), caffeic acid (04, 4-O^•), sinapic acid (05, 4-O^•), chlorogenic acid (06, 4-O^•), syringic acid (08, 4-O^•), protocatechuic acid (10, 4-O^•), ferulic acid (12, 4-O^•), p-coumaric acid (13, 4-O^•), methyl vanillate (15, 4-O^•) and ferulic acid ethyl ester (18, 4-O^•).

The CH = CH bridge in C = CCOOH or C = CCOOR takes part in the extended delocalization, particularly in its para O^• center. Relative to gallic acid (02), ferulic acid (12) possesses a conjugated CH = CH bridge between the phenyl and carbonyl group, which favors a resonance and conjugation effect on its para radical center. Therefore, the corresponding radical center has a lower spin density, leading to a lower 4-OH BDE value. Similarly, 4-O^• of caffeic acid (04) has a lower spin density (0.274) than 4-O^• (0.326) of protocatechuic acid (10), thus the OH BDE value (75.1 kcal/mol) of the former is the lower that of the latter (83.2 kcal/mol). The spin density distribution where the CH = CH bridge participates can be seen from other 4 radicals, including sinapic acid (05, 4-O^•), chlorogenic acid (06, 4-O^•), p-coumaric acid (13, 4-O^•) and ferulic acid ethyl ester (18, 4-O^•) as displayed in Fig. 3.

Orthodiphenolic functionalities also affected the structure-antioxidant activity relationships of 7 phenolic compounds, including propyl gallate (01), gallic acid (03), caffeic acid (04), chlorogenic acid (06), methyl gallate (09), protocatechuic acid (10) and ellagic acid (20). Orthodiphenolic functionalities could well secure a relatively complete spin density delocalization [13], which makes these compounds have higher antioxidant activities than the compounds with a single free hydroxy group on the ring. One of the most representative examples is that two orthodiphenolic functionalities of ellagic acid (20) mutually promote the extended spin density delocalization as shown in Fig. 3. This leads the 5-O^• radical of this compound to be more stable, and the compound has the highest antioxidant activity (pIC₅₀ = 5.00) among all the investigated compounds.

For 4 compounds with low antioxidant activities, including salicylic acid (07, pIC₅₀ = 1.10), 3-methylsalicylic acid (14, pIC₅₀ = 1.15), 3,5-dinitro salicylic acid (16, pIC₅₀ = 0.61), and 4-methylsalicylic acid (19, pIC₅₀ = 0.75), one common feature is that each compound has a single hydroxyl group in the ortho position next to the COOH group on the phenolic ring. These four compounds have a lower hydrogen donating ability than phenol under 4 different reaction environments due to their higher BDE values than phenol. The spin density (0.381) of phenol-O^• is relatively close to those of 2-O^• of these 4 parent compounds (0.375, 0.352, 0.354 and 0.365) (Fig. 3). Although there is a certain delocalization involving the O^• center, the C atoms of the ring, and the C = O constitution in COOH, this delocalization is not significant.

Structure-IP-PDE-activity relationships

Unlike BDEs, IPs appeared to be easily influenced by the solvent polarity since polar solvents may affect charge separation in a molecule [14]. As can be seen in Table 3, the IP values are sorted in the order of gas>benzene>ethanol>water for the same molecule, indicating the electron donating ability is more favorable in polar media. Cation radicals are sensitive to the polarity of different solvents, thus as expected, the IP values in water are lower than those in ethanol. As stated above, the IP is reaction enthalpies related to the first step of the SETPT mechanism, thus BDEs and IPs determine the thermodynamically preferred reaction pathway. The calculated IP values (Table 3) are significantly larger than the corresponding BDE values (Table 1), from the thermodynamic point of view, suggesting the occurrence probability of the HAT mechanism is greater than that of the SETPT mechanism in 4 reaction micro-environments. The PDEs under different media reveal different values, especially, significant average fall is found in water and ethanol. We also observed the PDEs in water are greater than those in ethanol, thus the proton dissociation ability of the studied compounds in ethanol is slightly stronger than that in water.

Download:

Table 3. IPs and PDEs of 20 investigated phenolic compounds and phenol calculated at the B3LYP/6-311++G(d,p) level ^{^a}.

https://doi.org/10.1371/journal.pone.0121276.t003

The explored compounds have different reaction capability for the SETPT mechanism due to their different IP values (Table 3) in 4 different environments. It should be noted that 4 compounds, including salicylic acid (07), 3-methylsalicylic acid (14), 3,5-dinitro salicylic acid (16) and 4-methylsalicylic acid (19), have larger IPs than phenol. Together with the evidence of the 4 compounds with larger BDE values than phenol (Table 1), we concluded the 4 compounds, which all possess a single phenolic OH group in the ortho position to the COOH group, have a similar reactivity tendency for both the HAT and SETPT mechanisms. In particular, 3,5-dinitro salicylic acid (16) has both the largest BDE values (97.1 and 97.2 kcal/mol, Table 1) and the largest IP values (154.4 and 163.1 kcal/mol, Table 3) in water and ethanol, respectively, demonstrating this compound has a lowest reaction probability for both the HAT and SETPT mechanisms in polar solution. Two strong electron-drawing NO₂ groups in this compound may induce a reduction in the electron density of the phenolic ring, so that this compound has a lowest reaction probability of the two mechanisms among all the investigated compounds.

Structure-PA-ETE-activity relationships

According to the SPLET mechanism, we examined the deprotonation of phenolic OH group(s) and electron transfer tendency in the gas phase and solvents by calculating PAs and ETEs. Table 4 displays the PA values decrease significantly in the order of gas>benzene>water>ethanol for the same molecule, due to the high solvation enthalpies of proton. Therefore, the deprotonation process is more likely to occur in polar solvents, such as water and ethanol. A number of studies have shown that solvents induce significant changes in enthalpies of charged species, which dominantly affects SETPT and SPLET energetics [18,41]. In benzene, PAs (Table 4) are lower than IPs (Table 3), but still higher than the corresponding BDEs (Table 1). By combining with the analysis above for both the HAT and SETPT mechanisms, HAT appears to represent the most thermodynamically probable reaction pathway. However, in water and ethanol, PAs are considerably lower than BDEs and IPs, thus SPLET represents the thermodynamically preferred mechanism. In comparison to the IP values (Table 3) of the neutral form, the corresponding ETE values are significantly lower. Hence, the single electron transfer process from the anionic form is more preferable than that from the neutral form.

Download:

Table 4. PAs and ETEs of 20 investigated phenolic compounds and phenol calculated at the B3LYP/6-311++G(d,p) level ^{^a}.

https://doi.org/10.1371/journal.pone.0121276.t004

All the compounds have lower PA values than phenol (345.1, 102.5, 48.1 and 44.2 kcal/mol) in 4 different environments (Table 4), revealing the compounds possess the ability to scavenge DPPH radicals to some degrees as they take effects on the radical scavenging by the SPLET mechanism, if possible. Three compounds, including salicylic acid (07, pIC₅₀ = 1.10), 3-methylsalicylic acid (14, pIC₅₀ = 1.15) and 4-methylsalicylic acid (19, pIC₅₀ = 0.75), have relatively high PA values (42.5, 43.3 and 38.2 kcal/mol) in ethanol (Table 4), suggesting these compounds are not susceptible to the SPLET mechanism. As described above, the 3 compounds were unfavorable for both the HAT and SETPT mechanisms, because the 3 compounds with a common characteristic, i.e., a single OH group in the ortho position to the COOH group, have relatively low antioxidant activities in ethanol.

One of interesting findings is that 3,5-dinitro salicylic acid (16) has the lowest PA values (295.9, 62.2, 19.1 and 14.3 kcal/mol), suggesting this compound can release H⁺ from the 2-OH group much easier than other compounds. Two NO₂ groups, as strong electron-withdrawing groups contained in this compound, bring out a reduction in the electron density of the phenol ring, thereby making the strongest H⁺ donation ability. The first step is easy to occur, however, the ETE value of the second step is the highest among all the compounds, meaning that the second step is very difficult to react, as evidenced by a relatively low activity (pIC₅₀ = 0.75). In short, substituents with different electronegativities can alter the electron density of the phenolic ring, thereby to influence the abstraction of H^• and H⁺ release, and the occurrence probability of the individual mechanism.

HOMO orbitals

Within the framework of molecular orbital theory, important information on the working mechanisms of antioxidants can be derived from the frontier highest occupied molecular orbital (HOMO) energy. Generally, the lower the HOMO energy is, the weaker the molecule donating electron ability is. On the contrary, a higher HOMO energy implies that the molecule is a good electron donor [42]. Because the H abstraction reaction involves electron transfer, the HOMO composition of a phenolic compound can provide a qualitative data to identify its active site for the scavenging radical activity [11,43]. For the compounds with two or more phenol OH groups, Fig. 4 shows there are more atomic polar tensors [44] atomic charges in the more probable electron-donating sites; moreover, those easily attacked sites consistently respond to the stable O^• centers with small spin density values as displayed in Fig. 3.

Download:

Fig 4. HOMO of 20 investigated phenolic compounds calculated at the B3LYP/6-311++G(d,p) levels of theory in ethanol.

The numbers indicates atomic polar tensor charges.

https://doi.org/10.1371/journal.pone.0121276.g004

Interpretation of action mechanisms by QSAR models

On the basis of the calculated thermodynamics parameters in polar media (water and ethanol) by DFT calculations, we established QSAR models to explore the characterization ability of the parameters on the reaction characteristics, thereby to further help interpret the mechanisms of the studied compounds. Table 5 displays all the QSAR models exhibit qualified fitting and predictive abilities because regression modeling coefficients and cross-validated correlation coefficients fall in the range of Q²>0.500 and R²>0.600 [45], respectively. These demonstrated five sets of calculated thermodynamics parameters could characterize the antioxidant activities of the compounds, and our QSAR models favorably describe the thermodynamics-activity relationships of these natural phenolic compounds in polar media. Particularly, the PA-ETE-based models have larger R²_water of 0.739, R²_ethanol of 0.745, Q²_water of 0.674 and Q²_ethanol of 0.689) than those of the BDE-based and IP-PDE-based models, suggesting PAs and ETEs are more effective descriptors for explaining the radical scavenging activities of these compounds relative to BDE and IP-PDE descriptors in polar solvents [11].

Download:

Table 5. QSAR modeling results of the thermodynamics-activity relationship for 20 phenolic compounds in water and ethanol ^{^a}.

https://doi.org/10.1371/journal.pone.0121276.t005

As shown in Table 5, the absolute values of the t-test values for all the variables (BDE, IP, PDE, PA and ETE) are larger than 2.000, thus the variables significantly influence the dependent variables (pIC₅₀). Each variable in the QSAR models has a negative coefficient, indicating the variables have negative effects on the radical scavenging activities, consistent with the well known fact that the smaller the thermodynamic parameters, the larger the corresponding activities. We found the PDE coefficients (−0.292 and −0.325) are smaller than the IP coefficients (-0.246 and -0.250) for the IP-PDE-activity models in both water and ethanol, while the PA coefficients (-0.293 and -0.299) are smaller than the ETE coefficients (-0.238 and -0.244) for the PA-ETE-activity models in both water and ethanol. Therefore, the proton dissociation ability plays a significant effect on the activities of the explored phenolic compounds in polar media such as water and ethanol.

Our QSAR modeling results reveal the fitting and predictive abilities of the PA-ETE-activity models are distinctly higher those of BDE-activity models (Table 5), which is supported by the conclusion on the working mechanisms of the explored compounds by DFT calculations, specifically, from the thermodynamic point of view, the SPLET mechanism would play a dominant role in water and ethanol. Therefore, one could try to choose the PA-ETE-activity model to predict antioxidant activities of new phenolic acids and derivatives in polar solvents.

Conclusions

We employ a combination of computational and experimental approaches, including DPPH^• scavenging activity determination, DFT calculations and QSAR modeling, to investigate the structure-thermodynamic-antioxidant activity relationships of 20 phenolic acids and derivatives in different reaction environments. The conclusions are as follows: (1) The phenolic radicals can be stabilized by two main factors: the possibility to establish intra-molecular hydrogen bonds and the extended delocalization and conjugation of the electrons enhanced by resonance phenomena. (2) These calculated descriptors (BDE, IP, PDE, PA and ETE) can be used to effectively describe the mechanisms (HAT, SETPT and SPLET) of the studied compounds in individual micro-environments, which demonstrates HAT is the more favorable mechanism in the gas phase and benzene, whereas the SPLET mechanism prefers in water and ethanol. (3) Our QSAR models can characterize the structure-antioxidant activity relationships of the studied compounds in polar media.

Our future work will focus on the structure-activity relationships of more phenolic acids and derivatives using combined experimental and computational approaches. Investigations of the reaction kinetics of the phenolic compounds with DPPH^• are of significance not only for enriching the knowledge of chemical mechanism(s) of action, but also for stimulating the discovery of effective drugs, food additives, or other functional molecules.

Supporting Information

S1 File. Table A.

Structures and experimental radical scavenging activities of selected phenolic acids and derivatives. Fig. A. Optimized structures of seven compounds with hydrogen bond(s) between orthodiphenolic functionalities calculated at the B3LYP/6-311++G(d,p) levels of theory in ethanol.

https://doi.org/10.1371/journal.pone.0121276.s001

(DOC)

Author Contributions

Conceived and designed the experiments: JZ GZL. Performed the experiments: YZC HZX. Analyzed the data: YZC HZX JZ GZL. Contributed reagents/materials/analysis tools: JZ GZL. Wrote the paper: YZC HZX JZ GZL.

References

1. Robbins RJ. Phenolic acids in foods: An overview of analytical methodology. J Agric Food Chem 2003;51: 2866–2887. pmid:12720366
- View Article
- PubMed/NCBI
- Google Scholar
2. Jayasena T, Poljak A, Smythe G, Braidy N, Munch G, Sachdev P. The role of polyphenols in the modulation of sirtuins and other pathways involved in Alzheimer's disease. Ageing Res Rev 2013;12: 867–883. pmid:23831960
- View Article
- PubMed/NCBI
- Google Scholar
3. Cherblanc FL, Davidson RWM, Di Fruscia P, Srimongkolpithak N, Fuchter MJ. Perspectives on natural product epigenetic modulators in chemical biology and medicine. Nat Prod Rep 2013;30: 605–624. pmid:23396528
- View Article
- PubMed/NCBI
- Google Scholar
4. Manach C, Scalbert A, Morand C, Remesy C, Jimenez L. Polyphenols: food sources and bioavailability. Am J Clin Nutr 2004;79: 727–747. pmid:15113710
- View Article
- PubMed/NCBI
- Google Scholar
5. Ajitha MJ, Mohanlal S, Suresh CH, Jayalekshmy A. DPPH radical scavenging activity of Tricin and its conjugates isolated from "javara" Rice Bran: A density functional theory study. J Agric Food Chem 2012;60: 3693–3699. pmid:22397636
- View Article
- PubMed/NCBI
- Google Scholar
6. 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–1183. pmid:11456671
- View Article
- PubMed/NCBI
- Google Scholar
7. Fiorucci S, Golebiowski J, Cabrol-Bass D, Antonczak S. DFT study of quercetin activated forms involved in antiradical, antioxidant, and prooxidant biological processes. J Agric Food Chem 2007;55: 903–911. pmid:17263492
- View Article
- PubMed/NCBI
- Google Scholar
8. Litwinienko G, Ingold KU. Solvent effects on the rates and mechanisms of reaction of phenols with free radicals. Acc Chem Res 2007;40: 222–230. pmid:17370994
- View Article
- PubMed/NCBI
- Google Scholar
9. Tishchenko O, Truhlar DG, Ceulemans A, Nguyen MT. A unified perspective on the hydrogen atom transfer and proton-coupled electron transfer mechanisms in terms of topographic features of the ground and excited potential energy surfaces as exemplified by the reaction between phenol and radicals. J Am Chem Soc 2008;130: 7000–7010. pmid:18465862
- View Article
- PubMed/NCBI
- Google Scholar
10. Leopoldini M, Russo N, Toscano M. The molecular basis of working mechanism of natural polyphenolic antioxidants. Food Chem 2011;125: 288–306.
- View Article
- Google Scholar
11. Xue Y, Zheng Y, An L, Dou Y, Liu Y. Density functional theory study of the structure-antioxidant activity of polyphenolic deoxybenzoins. Food Chem 2014;151: 198–206. pmid:24423521
- View Article
- PubMed/NCBI
- Google Scholar
12. Javan AJ, Javan MJ, Tehrani ZA. Theoretical investigation on antioxidant activity of bromophenols from the marine red alga Rhodomela confervoides: H-atom vs electron transfer mechanism. J Agric Food Chem 2013;61: 1534–1541. pmid:23347257
- View Article
- PubMed/NCBI
- Google Scholar
13. 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. pmid:23870923
- View Article
- PubMed/NCBI
- Google Scholar
14. 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 Int 2012;48: 538–543.
- View Article
- Google Scholar
15. Sun T, Ho C-T. Antioxidant activities of buckwheat extracts. Food Chem 2005;90: 743–749.
- View Article
- Google Scholar
16. Xie J, Schaich KM. Re-evaluation of the 2,2-Diphenyl-1-picrylhydrazyl free radical (DPPH) assay for antioxidant activity. J Agric Food Chem 2014;62: 4251–4260. pmid:24738928
- View Article
- PubMed/NCBI
- Google Scholar
17. Nenadis N, Zhang HY, Tsimidou MZ. Structure-antioxidant activity relationship of ferulic acid derivatives: effect of carbon side chain characteristic groups. J Agric Food Chem 2003;51: 1874–1879. pmid:12643644
- View Article
- PubMed/NCBI
- Google Scholar
18. Lengyel J, Rimarcik J, Vaganek A, Klein E. On the radical scavenging activity of isoflavones: thermodynamics of O-H bond cleavage. Phys Chem Chem Phys 2013;15: 10895–10903. pmid:23698223
- View Article
- PubMed/NCBI
- Google Scholar
19. Pratt DA, DiLabio GA, Brigati G, Pedulli GF, Valgimigli L. 5-Pyrimidinols: novel chain-breaking antioxidants more effective than phenols. J Am Chem Soc 2001;123: 4625–4626. pmid:11457259
- View Article
- PubMed/NCBI
- Google Scholar
20. Wijtmans M, Pratt DA, Valgimigli L, DiLabio GA, Pedulli GF, Porter NA. 6-Amino-3-pyridinols: towards diffusion-controlled chain-breaking antioxidants. Angew Chem Int Ed Engl 2003;42: 4370–4373. pmid:14502714
- View Article
- PubMed/NCBI
- Google Scholar
21. Feng Y, Liu L, Wang J-T, Huang H, Guo Q-X. Assessment of experimental bond dissociation energies using composite ab initio methods and evaluation of the performances of density functional methods in the calculation of bond dissociation energies. J Chem Inf Comput Sci 2003;43: 2005–2013. pmid:14632451
- View Article
- PubMed/NCBI
- Google Scholar
22. Hansch C, Maloney PP, Fujita T, Muir RM. Correlation of biological activity of phenoxyacetic acids with hammett substituent constants and partition coefficients. Nature 1962;194: 178–180.
- View Article
- Google Scholar
23. Liang G, Yang L, Chen Z, Mei H, Shu M, Li Z. A set of new amino acid descriptors applied in prediction of MHC class I binding peptides. Eur J Med Chem 2009;44: 1144–1154. pmid:18662841
- View Article
- PubMed/NCBI
- Google Scholar
24. Puzyn T, Leszczynska D, Leszczynski J. Toward the development of "nano-QSARs": advances and challenges. Small 2009;5: 2494–2509. pmid:19787675
- View Article
- PubMed/NCBI
- Google Scholar
25. Puzyn T, Rasulev B, Gajewicz A, Hu X, Dasari TP, Michalkova A, et al. Using nano-QSAR to predict the cytotoxicity of metal oxide nanoparticles. Nat Nanotechnol 2011;6: 175–178. pmid:21317892
- View Article
- PubMed/NCBI
- Google Scholar
26. Ianni JC, Annamalai V, Phuan P-W, Panda M, Kozlowski MC. A priori theoretical prediction of selectivity in asymmetric catalysis: design of chiral catalysts by using quantum molecular interaction fields. Angew Chem Int Ed 2006;45: 5502–5505. pmid:16906498
- View Article
- PubMed/NCBI
- Google Scholar
27. Cherkasov A, Muratov EN, Fourches D, Varnek A, Baskin II, Cronin M, et al. QSAR modeling: Where have you been? Where are you going to? J Med Chem 2014;57: 4977–5010. pmid:24351051
- View Article
- PubMed/NCBI
- Google Scholar
28. Cumming JG, Davis AM, Muresan S, Haeberlein M, Chen H. Chemical predictive modelling to improve compound quality. Nat Rev Drug Discov 2013;12: 948–962. pmid:24287782
- View Article
- PubMed/NCBI
- Google Scholar
29. Amić A, Marković Z, Dimitrić Marković JM, Stepanić V, Lučić B, Amić D. Towards an improved prediction of the free radical scavenging potency of flavonoids: The significance of double PCET mechanisms. Food Chem 2014;152: 578–585. pmid:24444978
- View Article
- PubMed/NCBI
- Google Scholar
30. Stepanic V, Troselj KG, Lucic B, Markovic Z, Amic D. Bond dissociation free energy as a general parameter for flavonoid radical scavenging activity. Food Chem 2013;141: 1562–1570. pmid:23790952
- View Article
- PubMed/NCBI
- Google Scholar
31. Sánchez-Moreno C. Review: Methods used to evaluate the free radical scavenging activity in foods and biological systems. Food Sci Technol Int 2002;8: 121–137.
- View Article
- Google Scholar
32. Becke AD. Density‐functional thermochemistry. III. The role of exact exchange. J Chem Phys 1993;98: 5648–5652.
- View Article
- Google Scholar
33. Binkley JS, Pople JA, Hehre WJ. Self-consistent molecular orbital methods. 21. Small split-valence basis sets for first-row elements. J Am Chem Soc 1980;102: 939–947.
- View Article
- Google Scholar
34. Cancès E, Mennucci B, Tomasi J. A new integral equation formalism for the polarizable continuum model: Theoretical background and applications to isotropic and anisotropic dielectrics. J Chem Phys 1997;107: 3032–3041.
- View Article
- Google Scholar
35. Bartmess JE. Thermodynamics of the electron and the proton. J Phys Chem 1994;98: 6420–6424.
- View Article
- Google Scholar
36. Rimarcik J, Lukeš V, Klein E, Ilčin M. Study of the solvent effect on the enthalpies of homolytic and heterolytic N-H bond cleavage in p-phenylenediamine and tetracyano-p-phenylenediamine. J Mol Struct THEOCHEM 2010;952: 25–30.
- View Article
- Google Scholar
37. Parker VD. Homolytic bond (H-A) dissociation free energies in solution. Applications of the standard potential of the (H+/H.bul.) couple. J Am Chem Soc 1992;114: 7458–7462.
- View Article
- Google Scholar
38. Bizarro MM, Costa Cabral BJ, Borges dos Santos RM, Martinho Simões JA. Substituent effects on the O–H bond dissociation enthalpies in phenolic compounds: Agreements and controversies. Pure Appl Chem 1999;71: 1249–1256.
- View Article
- Google Scholar
39. Frisch MJ, Trucks GW, Schlegel HB, Scuseria GE, Robb MA, Cheeseman JR, et al. Gaussian 09 2009: Gaussian, Inc.: Wallingford CT.
40. Chou KC, Shen HB. Recent progress in protein subcellular location prediction. Anal Biochem 2007;370: 1–16. pmid:17698024
- View Article
- PubMed/NCBI
- Google Scholar
41. Leopoldini M, Pitarch IP, Russo N, Toscano M. Structure, conformation, and electronic properties of Apigenin, Luteolin, and Taxifolin antioxidants. A first principle theoretical study. J Phys Chem A 2003;108: 92–96.
- View Article
- Google Scholar
42. Karelson M, Lobanov VS, Katritzky AR. Quantum-chemical descriptors in QSAR/QSPR studies. Chem Rev 1996;96: 1027–1044. pmid:11848779
- View Article
- PubMed/NCBI
- Google Scholar
43. Mazzone G, Toscano M, Russo N. Density functional predictions of antioxidant activity and UV spectral features of Nasutin A, Isonasutin, Ellagic Acid, and one of its possible derivatives. J Agric Food Chem 2013;61: 9650–9657. pmid:24024615
- View Article
- PubMed/NCBI
- Google Scholar
44. Cioslowski J. A new population analysis based on atomic polar tensors. J Am Chem Soc 1989;111: 8333–8336.
- View Article
- Google Scholar
45. Golbraikh A, Tropsha A. Beware of q2! J Mol Graph Model 2002;20: 269–276. pmid:11858635
- View Article
- PubMed/NCBI
- Google Scholar

[ref1] 1. Robbins RJ. Phenolic acids in foods: An overview of analytical methodology. J Agric Food Chem 2003;51: 2866–2887. pmid:12720366
View Article
PubMed/NCBI
Google Scholar

[2] View Article

[3] PubMed/NCBI

[4] Google Scholar

[ref2] 2. Jayasena T, Poljak A, Smythe G, Braidy N, Munch G, Sachdev P. The role of polyphenols in the modulation of sirtuins and other pathways involved in Alzheimer's disease. Ageing Res Rev 2013;12: 867–883. pmid:23831960
View Article
PubMed/NCBI
Google Scholar

[6] View Article

[7] PubMed/NCBI

[8] Google Scholar

[ref3] 3. Cherblanc FL, Davidson RWM, Di Fruscia P, Srimongkolpithak N, Fuchter MJ. Perspectives on natural product epigenetic modulators in chemical biology and medicine. Nat Prod Rep 2013;30: 605–624. pmid:23396528
View Article
PubMed/NCBI
Google Scholar

[10] View Article

[11] PubMed/NCBI

[12] Google Scholar

[ref4] 4. Manach C, Scalbert A, Morand C, Remesy C, Jimenez L. Polyphenols: food sources and bioavailability. Am J Clin Nutr 2004;79: 727–747. pmid:15113710
View Article
PubMed/NCBI
Google Scholar

[14] View Article

[15] PubMed/NCBI

[16] Google Scholar

[ref5] 5. Ajitha MJ, Mohanlal S, Suresh CH, Jayalekshmy A. DPPH radical scavenging activity of Tricin and its conjugates isolated from "javara" Rice Bran: A density functional theory study. J Agric Food Chem 2012;60: 3693–3699. pmid:22397636
View Article
PubMed/NCBI
Google Scholar

[18] View Article

[19] PubMed/NCBI

[20] Google Scholar

[ref6] 6. 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–1183. pmid:11456671
View Article
PubMed/NCBI
Google Scholar

[22] View Article

[23] PubMed/NCBI

[24] Google Scholar

[ref7] 7. Fiorucci S, Golebiowski J, Cabrol-Bass D, Antonczak S. DFT study of quercetin activated forms involved in antiradical, antioxidant, and prooxidant biological processes. J Agric Food Chem 2007;55: 903–911. pmid:17263492
View Article
PubMed/NCBI
Google Scholar

[26] View Article

[27] PubMed/NCBI

[28] Google Scholar

[ref8] 8. Litwinienko G, Ingold KU. Solvent effects on the rates and mechanisms of reaction of phenols with free radicals. Acc Chem Res 2007;40: 222–230. pmid:17370994
View Article
PubMed/NCBI
Google Scholar

[30] View Article

[31] PubMed/NCBI

[32] Google Scholar

[ref9] 9. Tishchenko O, Truhlar DG, Ceulemans A, Nguyen MT. A unified perspective on the hydrogen atom transfer and proton-coupled electron transfer mechanisms in terms of topographic features of the ground and excited potential energy surfaces as exemplified by the reaction between phenol and radicals. J Am Chem Soc 2008;130: 7000–7010. pmid:18465862
View Article
PubMed/NCBI
Google Scholar

[34] View Article

[35] PubMed/NCBI

[36] Google Scholar

[ref10] 10. Leopoldini M, Russo N, Toscano M. The molecular basis of working mechanism of natural polyphenolic antioxidants. Food Chem 2011;125: 288–306.
View Article
Google Scholar

[38] View Article

[39] Google Scholar

[ref11] 11. Xue Y, Zheng Y, An L, Dou Y, Liu Y. Density functional theory study of the structure-antioxidant activity of polyphenolic deoxybenzoins. Food Chem 2014;151: 198–206. pmid:24423521
View Article
PubMed/NCBI
Google Scholar

[41] View Article

[42] PubMed/NCBI

[43] Google Scholar

[ref12] 12. Javan AJ, Javan MJ, Tehrani ZA. Theoretical investigation on antioxidant activity of bromophenols from the marine red alga Rhodomela confervoides: H-atom vs electron transfer mechanism. J Agric Food Chem 2013;61: 1534–1541. pmid:23347257
View Article
PubMed/NCBI
Google Scholar

[45] View Article

[46] PubMed/NCBI

[47] Google Scholar

[ref13] 13. 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. pmid:23870923
View Article
PubMed/NCBI
Google Scholar

[49] View Article

[50] PubMed/NCBI

[51] Google Scholar

[ref14] 14. 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 Int 2012;48: 538–543.
View Article
Google Scholar

[53] View Article

[54] Google Scholar

[ref15] 15. Sun T, Ho C-T. Antioxidant activities of buckwheat extracts. Food Chem 2005;90: 743–749.
View Article
Google Scholar

[56] View Article

[57] Google Scholar

[ref16] 16. Xie J, Schaich KM. Re-evaluation of the 2,2-Diphenyl-1-picrylhydrazyl free radical (DPPH) assay for antioxidant activity. J Agric Food Chem 2014;62: 4251–4260. pmid:24738928
View Article
PubMed/NCBI
Google Scholar

[59] View Article

[60] PubMed/NCBI

[61] Google Scholar

[ref17] 17. Nenadis N, Zhang HY, Tsimidou MZ. Structure-antioxidant activity relationship of ferulic acid derivatives: effect of carbon side chain characteristic groups. J Agric Food Chem 2003;51: 1874–1879. pmid:12643644
View Article
PubMed/NCBI
Google Scholar

[63] View Article

[64] PubMed/NCBI

[65] Google Scholar

[ref18] 18. Lengyel J, Rimarcik J, Vaganek A, Klein E. On the radical scavenging activity of isoflavones: thermodynamics of O-H bond cleavage. Phys Chem Chem Phys 2013;15: 10895–10903. pmid:23698223
View Article
PubMed/NCBI
Google Scholar

[67] View Article

[68] PubMed/NCBI

[69] Google Scholar

[ref19] 19. Pratt DA, DiLabio GA, Brigati G, Pedulli GF, Valgimigli L. 5-Pyrimidinols: novel chain-breaking antioxidants more effective than phenols. J Am Chem Soc 2001;123: 4625–4626. pmid:11457259
View Article
PubMed/NCBI
Google Scholar

[71] View Article

[72] PubMed/NCBI

[73] Google Scholar

[ref20] 20. Wijtmans M, Pratt DA, Valgimigli L, DiLabio GA, Pedulli GF, Porter NA. 6-Amino-3-pyridinols: towards diffusion-controlled chain-breaking antioxidants. Angew Chem Int Ed Engl 2003;42: 4370–4373. pmid:14502714
View Article
PubMed/NCBI
Google Scholar

[75] View Article

[76] PubMed/NCBI

[77] Google Scholar

[ref21] 21. Feng Y, Liu L, Wang J-T, Huang H, Guo Q-X. Assessment of experimental bond dissociation energies using composite ab initio methods and evaluation of the performances of density functional methods in the calculation of bond dissociation energies. J Chem Inf Comput Sci 2003;43: 2005–2013. pmid:14632451
View Article
PubMed/NCBI
Google Scholar

[79] View Article

[80] PubMed/NCBI

[81] Google Scholar

[ref22] 22. Hansch C, Maloney PP, Fujita T, Muir RM. Correlation of biological activity of phenoxyacetic acids with hammett substituent constants and partition coefficients. Nature 1962;194: 178–180.
View Article
Google Scholar

[83] View Article

[84] Google Scholar

[ref23] 23. Liang G, Yang L, Chen Z, Mei H, Shu M, Li Z. A set of new amino acid descriptors applied in prediction of MHC class I binding peptides. Eur J Med Chem 2009;44: 1144–1154. pmid:18662841
View Article
PubMed/NCBI
Google Scholar

[86] View Article

[87] PubMed/NCBI

[88] Google Scholar

[ref24] 24. Puzyn T, Leszczynska D, Leszczynski J. Toward the development of "nano-QSARs": advances and challenges. Small 2009;5: 2494–2509. pmid:19787675
View Article
PubMed/NCBI
Google Scholar

[90] View Article

[91] PubMed/NCBI

[92] Google Scholar

[ref25] 25. Puzyn T, Rasulev B, Gajewicz A, Hu X, Dasari TP, Michalkova A, et al. Using nano-QSAR to predict the cytotoxicity of metal oxide nanoparticles. Nat Nanotechnol 2011;6: 175–178. pmid:21317892
View Article
PubMed/NCBI
Google Scholar

[94] View Article

[95] PubMed/NCBI

[96] Google Scholar

[ref26] 26. Ianni JC, Annamalai V, Phuan P-W, Panda M, Kozlowski MC. A priori theoretical prediction of selectivity in asymmetric catalysis: design of chiral catalysts by using quantum molecular interaction fields. Angew Chem Int Ed 2006;45: 5502–5505. pmid:16906498
View Article
PubMed/NCBI
Google Scholar

[98] View Article

[99] PubMed/NCBI

[100] Google Scholar

[ref27] 27. Cherkasov A, Muratov EN, Fourches D, Varnek A, Baskin II, Cronin M, et al. QSAR modeling: Where have you been? Where are you going to? J Med Chem 2014;57: 4977–5010. pmid:24351051
View Article
PubMed/NCBI
Google Scholar

[102] View Article

[103] PubMed/NCBI

[104] Google Scholar

[ref28] 28. Cumming JG, Davis AM, Muresan S, Haeberlein M, Chen H. Chemical predictive modelling to improve compound quality. Nat Rev Drug Discov 2013;12: 948–962. pmid:24287782
View Article
PubMed/NCBI
Google Scholar

[106] View Article

[107] PubMed/NCBI

[108] Google Scholar

[ref29] 29. Amić A, Marković Z, Dimitrić Marković JM, Stepanić V, Lučić B, Amić D. Towards an improved prediction of the free radical scavenging potency of flavonoids: The significance of double PCET mechanisms. Food Chem 2014;152: 578–585. pmid:24444978
View Article
PubMed/NCBI
Google Scholar

[110] View Article

[111] PubMed/NCBI

[112] Google Scholar

[ref30] 30. Stepanic V, Troselj KG, Lucic B, Markovic Z, Amic D. Bond dissociation free energy as a general parameter for flavonoid radical scavenging activity. Food Chem 2013;141: 1562–1570. pmid:23790952
View Article
PubMed/NCBI
Google Scholar

[114] View Article

[115] PubMed/NCBI

[116] Google Scholar

[ref31] 31. Sánchez-Moreno C. Review: Methods used to evaluate the free radical scavenging activity in foods and biological systems. Food Sci Technol Int 2002;8: 121–137.
View Article
Google Scholar

[118] View Article

[119] Google Scholar

[ref32] 32. Becke AD. Density‐functional thermochemistry. III. The role of exact exchange. J Chem Phys 1993;98: 5648–5652.
View Article
Google Scholar

[121] View Article

[122] Google Scholar

[ref33] 33. Binkley JS, Pople JA, Hehre WJ. Self-consistent molecular orbital methods. 21. Small split-valence basis sets for first-row elements. J Am Chem Soc 1980;102: 939–947.
View Article
Google Scholar

[124] View Article

[125] Google Scholar

[ref34] 34. Cancès E, Mennucci B, Tomasi J. A new integral equation formalism for the polarizable continuum model: Theoretical background and applications to isotropic and anisotropic dielectrics. J Chem Phys 1997;107: 3032–3041.
View Article
Google Scholar

[127] View Article

[128] Google Scholar

[ref35] 35. Bartmess JE. Thermodynamics of the electron and the proton. J Phys Chem 1994;98: 6420–6424.
View Article
Google Scholar

[130] View Article

[131] Google Scholar

[ref36] 36. Rimarcik J, Lukeš V, Klein E, Ilčin M. Study of the solvent effect on the enthalpies of homolytic and heterolytic N-H bond cleavage in p-phenylenediamine and tetracyano-p-phenylenediamine. J Mol Struct THEOCHEM 2010;952: 25–30.
View Article
Google Scholar

[133] View Article

[134] Google Scholar

[ref37] 37. Parker VD. Homolytic bond (H-A) dissociation free energies in solution. Applications of the standard potential of the (H+/H.bul.) couple. J Am Chem Soc 1992;114: 7458–7462.
View Article
Google Scholar

[136] View Article

[137] Google Scholar

[ref38] 38. Bizarro MM, Costa Cabral BJ, Borges dos Santos RM, Martinho Simões JA. Substituent effects on the O–H bond dissociation enthalpies in phenolic compounds: Agreements and controversies. Pure Appl Chem 1999;71: 1249–1256.
View Article
Google Scholar

[139] View Article

[140] Google Scholar

[ref39] 39. Frisch MJ, Trucks GW, Schlegel HB, Scuseria GE, Robb MA, Cheeseman JR, et al. Gaussian 09 2009: Gaussian, Inc.: Wallingford CT.

[ref40] 40. Chou KC, Shen HB. Recent progress in protein subcellular location prediction. Anal Biochem 2007;370: 1–16. pmid:17698024
View Article
PubMed/NCBI
Google Scholar

[143] View Article

[144] PubMed/NCBI

[145] Google Scholar

[ref41] 41. Leopoldini M, Pitarch IP, Russo N, Toscano M. Structure, conformation, and electronic properties of Apigenin, Luteolin, and Taxifolin antioxidants. A first principle theoretical study. J Phys Chem A 2003;108: 92–96.
View Article
Google Scholar

[147] View Article

[148] Google Scholar

[ref42] 42. Karelson M, Lobanov VS, Katritzky AR. Quantum-chemical descriptors in QSAR/QSPR studies. Chem Rev 1996;96: 1027–1044. pmid:11848779
View Article
PubMed/NCBI
Google Scholar

[150] View Article

[151] PubMed/NCBI

[152] Google Scholar

[ref43] 43. Mazzone G, Toscano M, Russo N. Density functional predictions of antioxidant activity and UV spectral features of Nasutin A, Isonasutin, Ellagic Acid, and one of its possible derivatives. J Agric Food Chem 2013;61: 9650–9657. pmid:24024615
View Article
PubMed/NCBI
Google Scholar

[154] View Article

[155] PubMed/NCBI

[156] Google Scholar

[ref44] 44. Cioslowski J. A new population analysis based on atomic polar tensors. J Am Chem Soc 1989;111: 8333–8336.
View Article
Google Scholar

[158] View Article

[159] Google Scholar

[ref45] 45. Golbraikh A, Tropsha A. Beware of q2! J Mol Graph Model 2002;20: 269–276. pmid:11858635
View Article
PubMed/NCBI
Google Scholar

[161] View Article

[162] PubMed/NCBI

[163] Google Scholar

Figures

Abstract

Introduction

Principles and Methods

Chemicals

DPPH radical scavenging assay

DFT calculations

QSAR modeling and validation

Results and Discussion

Effects of intra-molecular hydrogen bonds on the stability of the radicals

Structure-BDE-activity relationships

Structure-IP-PDE-activity relationships

Structure-PA-ETE-activity relationships

HOMO orbitals

Interpretation of action mechanisms by QSAR models

Conclusions

Supporting Information

S1 File. Table A.

Author Contributions

References