Next Article in Journal
Instability of Alien Chromosome Introgressions in Wheat Associated with Improper Positioning in the Nucleus
Next Article in Special Issue
Computational Studies on Water-Catalyzed Mechanisms for Stereoinversion of Glutarimide Intermediates Formed from Glutamic Acid Residues in Aqueous Phase
Previous Article in Journal
Use of Multifactorial Treatments to Address the Challenge of Translating Experimental Myocardial Infarct Reduction Strategies
Previous Article in Special Issue
Substituent Effects on the Radical Scavenging Activity of Isoflavonoid
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

DFT Studies on the Antioxidant Activity of Naringenin and Its Derivatives: Effects of the Substituents at C3

1
College of Bee Science, Fujian Agriculture and Forestry University, Fuzhou 350002, China
2
Key Laboratory of Bioorganic Phosphorous Chemistry and Chemical Biology (Ministry of Education), Department of Chemistry, Tsinghua University, Beijing 100084, China
*
Author to whom correspondence should be addressed.
Int. J. Mol. Sci. 2019, 20(6), 1450; https://doi.org/10.3390/ijms20061450
Submission received: 18 February 2019 / Revised: 11 March 2019 / Accepted: 12 March 2019 / Published: 22 March 2019
(This article belongs to the Special Issue DFT Applications to Biomolecules and Complex Reactions)

Abstract

:
The radical scavenging activity of a flavonoid is largely influenced by its structure. The effects of the substituents at C3 position on the antioxidant activity of naringenin were carried out using the density functional theory (DFT) method. The reaction enthalpies related with the three well-established mechanisms were analyzed. Excellent correlations were found between the reaction enthalpies and Hammett sigma constants. Equations obtained from the linear regression can be helpful in the selection of suitable candidates for the synthesis of novel naringenin derivatives with enhanced antioxidant properties. In the gas and benzene phases, the antioxidant activity of naringenin was enhanced by the electron-donating substituents via weakening the bond dissociation enthalpy (BDE). In the water phase, it was strengthened by electron-withdrawing groups—via lowering the proton affinity (PA). The electronic effect of the substituent on the BDE of naringenin is mainly governed by the resonance effect, while that on the ionization potential (IP) and PA of naringenin is mainly controlled by the field/inductive effect.

1. Introduction

Excess formation of reactive free radicals by various enzymatic and non-enzymatic processes in the body has been found associated with the oxidation of biomolecules such as nucleic acids, proteins and lipids [1,2]. The oxidations of biomolecules are closely related with the initiation and progression of various diseases such as atherosclerosis, cardiovascular disease, neurodegenerative disease, cancer, and aging [1,3].
Antioxidants are chemical compounds that can trap free radical intermediates formed during the oxidative reactions and inhibit the oxidative reactions. The natural antioxidants have recently received notable attention because they can avoid or at least significantly reduce the oxidation of biomolecules by free radicals [4,5]. Among the multiple natural bioactive compounds, flavonoids as free radical chain-breaking antioxidants are well known. Flavonoids are a class of naturally occurring polyphenolic compounds that exist in a wide variety of foods and thus are important constituents of the human diet [6]. They are widely found in natural foods such as fruits, vegetables, cereals, teas, wines, honey, and bee pollen [6].
The flavonoid basic structure consists of two benzene rings (A and B) connected by an oxygen-containing pyrene ring (C) (Figure 1). To date, more than 6000 flavonoids have been identified. Flavonoids possess a broad range of biological activities including antibacterial, antiviral, anti-inflammatory and anti-ischemic activities [7,8,9,10], which are closely related with their free radical scavenging activities. It is reported that flavonoids play their anti-inflammatory action mainly via the ability to modulate production of free radical by phagocytic leukocytes [11]. Naringenin (4′,5,7-trihydroxyflavanone, Figure 1) is a dietary flavonoid abundant in foods such as citrus fruits, honey, and bee pollen. Naringenin has attracted increasing attention due to its positive health benefits in the human body. It has been reported to have several biological effects, such as monoamine oxidase inhibitory [12] and neuroprotective activities [13]. Monoamine oxidase inhibitory and neuroprotective activities of naringenin may provide protection against oxidative neurodegeneration and alleviate central nervous system disorders such as depression [14].
The antioxidant activity of phenolic compounds depends on their capacity to resist the detrimental effect of free radicals. There is a close relationship between antioxidant activity and structural properties. Among multiple structural features, the substituent effect is one of the most important factors that influence the antioxidant properties of phenolic compounds [15,16,17,18]. The studies on the influence of the substituent effect can be used to synthesize and select novel compounds with better antioxidant activity. Until now, there have been multiple studies that focus on the number and position of hydroxyl groups in relation to the biological activity of the flavonoids [19,20,21]. On the other hand, the substituent effect on the antioxidant activity of flavonoids is seldom known.
There are several mechanisms for the free radical scavenging action of phenolic compounds and the net result of all mechanisms is the same: transferring a hydrogen atom to the free radicals. The following mechanisms are widely accepted for phenolic compounds exerting their antioxidant activity [22,23,24]: (1) the hydrogen atom transfer (HAT) in which the radical abstracts the hydrogen by homolytic cleavage of the hydroxyl group (Equation (1)); (2) the single-electron transfer followed by proton transfer (SET-PT) takes place through two consecutive steps: an electron donates from the phenolic compounds (the first step, Equation (2)) following by a proton transferring from the formed cation radical (the second step, Equation (3)); (3) the sequential proton loss electron transfer (SPLET), which takes place in two steps: a proton transfers from the antioxidant (the first step, Equation (4)) following by an electron donating from the anion created in the first step (the second step, Equation (5)).
R• + ArOH → RH + ArO•
R• + ArOH → R + ArOH+
R + ArOH+• → RH + ArO•
ArOH → ArO + H+
ArO + R• + H+ → ArO• + RH
Understanding the role of substituent is of great importance for the preparation of novel phenolic compounds with enhanced antioxidant property. The computational method is a powerful tool that can provide valuable information about the antioxidant character of the phenolic compound. The work that is difficult to carry out by the experiment method can be performed by the computational method. Computational methods, especially the density functional theory (DFT), have been successfully used to evaluate the chemical properties, which are closely related with the antioxidant activity [22,23,24,25,26,27,28,29,30,31,32,33,34,35,36]. In this paper, the effects of various substituents on the antioxidant activity of naringenin were investigated by the DFT method. Reaction enthalpies related to HAT, SET-PT, and SPLET mechanisms were calculated. The derivatives of flavonoids with substituents at C3 positions are widely distributed in natural foods such as honey and bee pollens. In this work, various substituents (NH2, OMe, Me, OH, F, Cl, CHO, CF3, CN and NO2) covering the electron-withdrawing groups and electron-donating groups were placed at the C3 position of naringenin. The Hammett sigma constant has been one of the most widely used means to study and interpret the organic reactions and their mechanisms [18]. To find a relationship between substituents and the antioxidant activity, correlations of calculated reaction enthalpies with Hammett sigma constants of the substituents were investigated.

2. Results and Discussion

The atom numberings for the carbon atoms in the basic structure of a flavonoid are shown in Figure 1. For simplicity, the numberings for oxygen and hydrogen atoms are used to correspond to the number of attached carbon atoms. For example, atoms of OH group at C5′ are O5′ and H5′, at C3 are O3 and H3. The naringenin monomer is used as the same atom numbering as the flavonoid basic structure.
Reaction enthalpy characterizing the first step of the reaction mechanism is of great importance in evaluating the antioxidant action. In the three mechanisms, HAT is a one-step mechanism, while SET-PT and SPLET mechanisms take place via two steps. Hence, the following analyses mainly focus on the bond dissociation enthalpy (BDE), ionization potential (IP) and ionization potential (PA) characterizing the first steps of the three mechanisms.

2.1. HAT Mechanism

2.1.1. Calculated BDEs of Naringenin and its C3-Substituted Derivatives

In a HAT mechanism, a hydrogen atom abstracts from the phenolic hydroxyl group to the free radical via homolytic cleavage. In this mechanism, BDE is the parameter evaluating the antioxidant action of a compound. The lower the BDE, the stronger the antioxidant activity and the more important the role of the corresponding O−H group in the antioxidant action. The calculated O−H BDEs of naringenin and its C3-substituted derivatives are given in Table 1. The antioxidant activities of naringenin and its C3-OH substituted derivatives kaempferol (4′,3,5,7-tetrahydroxyflavanone) have been reported in previous work [37]. It is found that in the aqueous phase, the antioxidant activity of kaempferol is stronger than naringenin [37]. In this work, the BDEs of naringenin in different hydroxyl group are stronger than kaempferol in the water phase. It is consistent with the experimental results, thus, confirming the accuracy of our work.
For naringenin and its derivatives, they have a common feature: the lowest BDEs are at 4′−OH irrespective of the studied media. Therefore, the antioxidant capacities of naringenin and its derivatives are preferred to present by 4′−OH in the HAT mechanism. This observation may be attributed to that of the flavonoiud−O4′, where the radical is more stable in the studied phases.
ΔBDE representing the difference between BDE of the substituted naringenin and naringenin itself, reflects the effect of the substituents on O–H BDE. The calculated ΔBDEs are summarized in Table S1 (Supplementary Material). In the gas phase, the lowest effect of the substituent on the O-H BDE is found for Me group at C3 position, where BDEs vary only by 0.5, −0.8 and −0.2 kJ/mol for 4′-OH, 5-OH and 7-OH, respectively. The biggest effect of the substituent on the O-H BDE is found for NO2 group at C3 position, where BDEs vary by 7.2, 4.4, and 6.2 kJ/mol for 4′-OH, 5-OH, and 7-OH respectively. The similar trends are also observable in the benzene and water phases. Thus, the Me group at C3 position influences the antioxidant activity of naringenin weaker compared with other groups, while the NO2 group has the strongest impact on the antioxidant action.

2.1.2. Dependence of BDEs on Hammett Sigma Constants

The Hammett sigma constant has been one of the most widely used means to study and interpret the organic reactions and their mechanisms [18]. In this work, the Hammett sigma constants: σm and σp of the substituents were used. Table S2 (Supplementary Material) lists the σm and σp of the studied substituents collected by Hansch et al. [18]. They were obtained from the ionization of organic acids in solutions and have been successfully used to predict the equilibrium and rate constants of a family variety of reactions. To see the correlation of the Hammett sigma constants with the BDE of the substituted naringenin, the Pearson correlation was carried out. The P(BDEs, σ) are drawn as a histogram in Figure 2. In Figure 2, all of the P(BDEs, σ) are positive and larger than 0.8. Hence, the BDEs are positively and highly relevant with the Hammett sigma constants. The positive coefficients also character that electron-withdrawing groups at C3 position raise the BDE of naringenin, while electron-donating groups play the contrary role. The obtained results can be interpreted that electron-donating groups at the C3 position can stabilize the radical, while destabilizing the parent molecule. Hence, they decrease the O−H BDE. However, electron-withdrawing groups at the C3 position have the opposite effect, and their presence would lead to an increase in the O−H BDE. On the other hand, the electron-donating groups placed at the C3 position would strengthen the antioxidant activity of naringenin, while the electron-withdrawing groups have an opposite effect in a HAT mechanism.
In Figure 2, the P(BDEs, σp) are larger than P(BDEs, σm) and exceed 0.9. Therefore, the correlation between the BDE of naringenin and σp is better than that between the BDE and σm and the BDE of naringenin can be predicted better by σp. The equations obtained from the linear regressions between σp and BDE(4′−OH) are placed in the Figure 3A. They can be used to predict the BDE(4′−OH) by the Hammett sigma constants. As the 4′−OH is the strongest antioxidant hydroxyl group, the antioxidant activity of the substituted naringenin can also be predicted by the obtained equations.

2.1.3. The Electronic Effects of the Substituents on the BDEs

In this work, the electronic effects of the substituents on the BDEs have been considered. The electronic effect of the substituent is mainly composed of two parts: field/inductive effect represented by parameter F and resonance effect characterized by parameter R [18,38,39]. The Pearson correlation coefficients between the F/R and the BDEs are drawn as a histogram in Figure 2. In Figure 2, it can be seen that the correlation between BDE of naringenin and R is much better than that between BDE and F. The P(BDEs, R) exceeds 0.8. Therefore, the electronic effect of substituent at the C3 position on the BDE of naringenin is mainly governed by the resonance effect.

2.2. SET-PT Mechanism

2.2.1. Calculated IPs of Naringenin and its C3-Substituted Derivatives

Another important mechanism of the antioxidant process is SET-PT. The first step of this mechanism is transferring a single electron from a natural flavonoid. This step is significantly important for the SET-PT mechanism. It is characterized by IP. The lower the IP, the stronger the antioxidant activity. The calculated IPs for naringenin and its derivatives are given in Table 2. It is well known that the charged species are especially sensitive to the polarity change of the solvents. Thus, the IPs of naringenin and its derivatives change drastically with the increasing polarity of the environments, as can be seen in Table 2. The cause may be that the cation radicals of naringenin and its derivatives are more stable, and the delocalization and conjugation of the π-electrons are more delocalized in the polar environments.
ΔIP representing the difference between IP of the substituted naringenin and the naringenin itself reflects the effect of substituent on IP. The calculated ΔIPs are summarized in Table S3 (Supplementary Material). In the gas phase, the lowest effect of the substituent on the IP is found for Me group at C3 position, where IP varies only by 3.3 kJ/mol. The biggest effect of the substituent on the IP is found for NO2 group at C3 position, where IP changes by 41.4 kJ/mol. Similar trends are also observable in the benzene and water phases. Thus, the Me group at C3 position influences the antioxidant activity weaker compared with other groups, while the NO2 group has the strongest impact on the antioxidant action of naringenin.

2.2.2. Dependence of IPs on Hammett Sigma Constants

To see the correlation of the Hammett sigma constants with the IP of substituted naringenin, the Pearson correlation were carried out. The Pearson correlation coefficients between the σm/σp and the IPs are drawn as a histogram in Figure 2. In Figure 2, all of the P(IPs, σ) are positive and larger than 0.8. Hence, the IPs are positively and highly relevant with the Hammett sigma constants. The positive coefficient also characterizes that electron-donating groups at the C3 position reduce the IP of naringenin, while electron-withdrawing groups have the opposite effect. The obtained results can be interpreted that electron-donating groups at the C3 position can stabilize the cation radical and destabilize the parent molecule. Hence, they reduce the IPs. However, electron-donating groups at the C3 position have the opposite effect, and their presence would lead to an increase in the IP. Therefore, the electron-donating groups placed at the C3 position would enhance the antioxidant activity of naringenin, while the electron-withdrawing groups have an opposite effect in SET-PT mechanism.
In Figure 2, the P(IPs, σm) are larger than P(IPs, σp) and exceed 0.9. Therefore, the correlation between the IP in naringenin and σm is better than that between the IP and σp, and the IP of naringenin derivatives with the substituent at C3 position can be predicted well by σm. The equations obtained from the linear regressions between σm and IP are placed in Figure 3B. They can be used to predict the IP and antioxidant activity of the substituted naringenin by the σm. This can be useful in the selection of suitable candidates for the synthesis of novel naringenin derivatives with enhanced antioxidant properties.

2.2.3. The Electronic Effects of the Substituents on the IPs

In this work, the electronic effects of the substituents on the IPs have been considered. The Pearson correlation coefficients between the F/R and the IPs are drawn as a histogram in Figure 2. In Figure 2, it can be seen that the correlation between IP of naringenin and F is much better than that between IP and R. The P(IPs, F) exceed 0.9. Therefore, the electronic effect of substituent at the C3 position on the IP of naringenin is mainly governed by the field/inductive effect.

2.3. SPLET mechanism

2.3.1. Calculated PAs of Naringenin and its C3-Substituted Derivatives

In polar phases, SPLET mechanism plays an important role in scavenging free radicals of the phenolic compounds. The heterolytic O‒H bond dissociation is the first step of this mechanism. This step is significantly important for SPLET mechanism and is characterized by PA. The lower the PA, the stronger the antioxidant activity and the more important the role of the corresponding O−H bond in the antioxidant action. The calculated O−H PAs of naringenin and its derivatives are given in Table 3. The antioxidant activities of naringenin and its C3-OH substituted derivatives kaempferol (4′,3,5,7-tetrahydroxyflavanone) have been reported in previous work [37]. It is found that in the aqueous phase, the antioxidant activity of kaempferol is stronger than naringenin [37]. In this work, the PAs of naringenin in a different hydroxyl group are stronger than those kaempferol in the water phase. It is consistent with the experimental results, thus, confirming the accuracy of our work.
For naringenin and its derivatives, they have common features: the lowest PAs are at 7−OH irrespective of the studied media. Therefore, the antioxidant capacities of naringenin and its derivatives are preferred to present by 7−OH in the SPLET mechanism. This observation may be attributed to the fact that the acid strength of 7−OH is stronger and the flavonoiud−O7 is more stable in the studied phases.
The charged species are especially sensitive to the polarity change of the solvents and the solvation enthalpies of the proton and anion are relatively high. Thus, from the gas phase to the solvent phases, the PAs decrease drastically. For example, the average deviation of PA(4′−OH) between the gas phase and the water phase reaches 1244 kJ/mol. This result indicates that the first step of the SPLET mechanism for naringenin and its derivatives are more favored in the polar phases.
ΔPA representing the difference between PA of the substituted naringenin and naringenin itself, reflects the effect of substituent on O–H PA. The calculated ΔPAs are summarized in Table S4 (Supplementary Material). In the gas phase, the lowest effect of the substituent on the O-H PA is found for Me group at C3 position, where PAs vary only by 0.6, 2.0 and 2.9 kJ/mol for 4′-OH, 5-OH and 7-OH, respectively. The biggest effect of the substituent on the O-H PA is found for NO2 group at C3 position, where PAs vary by −8.5, −45.7 and −38.4 kJ/mol for 4′-OH, 5-OH and 7-OH respectively. The similar trends are also observable in the benzene and water phases. Thus, the Me group at C3 position influences the antioxidant activity of naringenin weaker compared with other groups, while the NO2 group has the strongest impact on the antioxidant action.

2.3.2. Dependence of PAs on Hammett Sigma Constants

To see the correlation of the Hammett sigma constants with the PAs of substituents naringenin, the Pearson correlation were carried out. The Pearson correlation coefficients between the σm/σp and the PAs are drawn as a histogram in Figure 2. In Figure 2, all of the P(PAs, σ) are negative and the absolute value of P(PAs, σ) are larger than 0.8. Hence, the PAs are negatively and highly relevant with the Hammett sigma constants. The negative value of the coefficient also characterizes that electron-donating groups at C3 position of naringenin increase the PA, while electron-withdrawing groups have an opposite effect. The obtained results can be interpreted that electron-donating groups at the C3 position of naringenin can stabilize the neutral molecule and destabilize the anion. Hence, they increase the O−H PA. However, electron-withdrawing groups in the C3 position have the opposite effect, and their presence would lead to a decrease in the O−H PA. Therefore, the electron-donating groups placed at the C3 position would lower the antioxidant activity of naringenin while the electron-donating groups have an opposite effect in SPLET mechanism.
In Figure 2, the P(PAs, σm) are larger than P(PAs, σp) and exceed 0.9. Therefore, the correlation between the PA of naringenin and σm is better than that between the PA and σp and the PA of O−H can be predicted better by σm. The equations obtained from the linear regressions between σm and PA(7−OH) are placed in the Figure 3C. They can be used to predict the PA(7−OH) by the Hammett sigma constants. As the 7−OH is the strongest antioxidant hydroxyl group in SPLET mechanism, the antioxidant activity of the substituted naringenin can also be predicted by the equations.

2.3.3. The Electronic Effects of the Substituents on the PAs

The Pearson correlation coefficients between the F/R and the PAs are drawn as a histogram in Figure 2. In Figure 2, it can observe that the correlation between PA and F is much better than that between PA and R. The P(PAs, F) exceeds 0.8. Therefore, the electronic effect of the substituent at the C3 position on the PA of naringenin is mainly governed by the field/inductive effect.

2.4. The antioxidant Activity Strength Influenced by the Substituted Groups

In general, free energy (∆rG = ∆rHTrS) is the criterion to distinguish the thermodynamically favored process. However, in the case of the studied mechanisms, the absolute values of the entropic term (−TrS) reach only few units of kJ/mol and the free energies are only shifted in comparison to the corresponding enthalpies (∆rH).
Reaction enthalpy characterizing the first step of reaction mechanism is of great importance in evaluating the antioxidant action. Thus, the BDE, IP and PA can be used to determine the thermodynamically preferred reaction pathway for the antioxidant reaction of naringenin and its C3-substituted derivatives. Analyzed the data in Table 1, Table 2 and Table 3, it can be deduced that the BDEs of 4′−OH (the strongest antioxidant hydroxyl group following the HAT mechanism) are significantly lower than the IPs and the PAs of 7−OH (the strongest antioxidant hydroxyl group following the SPLET mechanism) for naringenin and its derivatives in the gas and benzene phases. In the water phase, the PAs of 7−OH are significantly smaller than the lowest BDEs of 4′−OH and IPs. Hence, from the thermodynamic point of view, HAT is the most favored mechanism in the gas and benzene phases, while SPLET is more preferred than HAT and SET-PT in the water phase.
In the studied phases, the calculated IPs for naringenin and its derivatives are significantly higher than the lowest BDEs and the lowest PAs. Therefore, from the thermodynamic point, SET-PT is the least favored mechanism in the studied phases.
In HAT and SPLET mechanisms, the antioxidant activities of flavonoids are characterized by the lowest BDE and lowest PA, respectively. In the gas and benzene phases, the electron-donating substituents placed at the C3 position would enhance the antioxidant activity of naringenin via lowering the BDE of 4′−OH, while the electron-withdrawing groups have the opposite effect. In water phase, the electron-donating substituents placed at the C3 position would weaken the antioxidant activity of naringenin via enlarging the PA of 7−OH, while the electron-withdrawing groups have an opposite effect. Among the studied compounds, the NH2 substituted derivative is the strongest antioxidant in the gas and benzene phases, while the NO2 substituted derivative is the strongest one in water phase.

3. Materials and Methods

3.1. Computational Details

All calculations were carried out using Gaussian 09 program package [40]. The gas-phase geometries of neutral compounds and the corresponding radical or ionic structures were optimized using the M062X functional without any constraints. The 6-311++G ** basis set including diffuse and polarization functions was used. It allows balanced treatment of neutral compounds and corresponding radicals or ions. The harmonic vibrational frequencies are calculated to confirm that the optimized geometry correctly corresponds to a local minimum, which has only real frequencies. As the solvent would alter the thermodynamically favored mechanism, benzene and water phases were also considered in this work. The solvent effects were conducted by employing the SMD continuum solvent model based on the optimized gas-phase geometries. Single-point electronic energies were then performed at the M062X/6-311+G ** level of theory in different solvent environments based on the optimized geometries.
The reaction enthalpies related to HAT, SET-PT and SPLET mechanisms are usually denoted as follows:
BDE: bond dissociation enthalpy related to HAT. The calculated equation for BDE is:
BDE = H(ArO•) + H(H•) – H(ArOH)
IP: ionization potential related to the first step of SET-PT mechanism. The calculated equation for IP is:
IP = H(ArOH+•) + H(e) – H(ArOH)
PDE: proton dissociation enthalpy related to the second step of SET-PT mechanism. The calculated equation for PDE is:
PDE = H(ArO•) + H(H+) – H(ArOH+•)
PA: proton affinity of phenoxide anion related to the first step of SPLET. The calculated equation for PA is:
PA = H(ArO) + H(H+) – H(ArOH)
ETE: electron transfer enthalpy related to the second step of SPLET. The calculated equation for ETE is:
ETE = H(ArO•) + H(e) – H(ArO)
The calculated gas and solvent phases enthalpies for H+, e and H• were obtained from the references [41,42,43]. The calculated enthalpies for H+ in the gas, water and benzene phases are 6.2 kJ/mol, −1050.0 kJ/mol and −870.6 kJ/mol respectively. The calculated enthalpies for e in the gas, water and benzene phases are 3.1 kJ/mol, −15.1 kJ/mol and −6.6 kJ/mol respectively. The calculated enthalpies for H• in the gas, water and benzene phases are −1306.6 kJ/mol, −1310.6 kJ/mol and −1301.6 kJ/mol respectively.

3.2. Statistics Analysis

Pearson correlation is an effective method to study a broad class of relationships among variables [44]. It is based on the covariance matrix of the data and Pearson correlation coefficient is the parameter to evaluate the strength of relationship between two vectors. Normally, Pearson correlation coefficient (Px,y) between two vectors x and y is calculated using the following equation:
P x , y   =   x y     x y N ( x 2     ( x 2 ) N ) ( y 2     ( y 2 ) N )
where N refers to the size of the signature array and is 11 in this work. The Pearson correlation coefficient is symmetric: Px,y = Py,x. The Px,y is among 1 to −1. The more close to 1/−1 the Px,y, the more related the x and y. The Px,y equal to 1 or −1 corresponds to data points lying exactly on a line, or to a bivariate distribution entirely supported on a line.

4. Conclusions

Radical scavenging activity of the flavonoid is largely influenced by the structure of the compound. Understanding the role of the substituents is of great importance for the preparation of novel phenolic compounds with an enhanced antioxidant property. In this study, systematic calculations about the influence of the substituents at C3 position on the antioxidant reactivity of naringenin were carried out using DFT calculations.
The most stable geometries of the neutral compounds and the corresponding radical or ionic structures were optimized using the M062X/6-311++G ** method. Based on the optimized geometries, reaction enthalpies (BDE, IP and PA) related to HAT, SET-PT, and SPLET mechanisms have been discussed in detail. Excellent correlations were found between the BDE/IP/PA and Hammett sigma constants. Equations obtained from the linear regression can be useful in the selection of suitable candidates for the synthesis of novel flavonoids derivatives with enhanced antioxidant properties. Besides, it is also found that:
(1) Electron-withdrawing groups enhance BDE and IP, while the electron-donating groups reduce the BDE and IP. The effect of the substituents on PA is on the contrary.
(2) The electronic effect of the substituents on the BDEs of naringenin is mainly governed by the resonance effect, while that on the IPs and PAs of naringenin is mainly controlled by field/inductive effect.
(3) In the gas and benzene phases, HAT is the most favorable mechanism and the antioxidant capacities of naringenin and its derivatives are preferred to present by 4′−OH. In water phase, SPLET is more favored than HAT and SET-PT and the antioxidant capacities are preferred to present by 7−OH.
(4) In the gas and benzene phases, the electron-donating substituents placed at the C3 position would enhance the antioxidant activity of naringenin via weakening the BDE of 4′−OH. In the water phase, they will reduce the antioxidant activity by strengthening the PA of 7−OH. Contrary results occur for the electron-withdrawing group. Besides, the NH2 substituted derivative is the strongest antioxidant in the gas and benzene phases, while the NO2 substituted derivative is the strongest one in the water phase.

Supplementary Materials

Supplementary materials can be found at https://www.mdpi.com/1422-0067/20/6/1450/s1. Table S1: ΔBDE in kJ/mol obtained by the difference between BDE of the substituted naringenin and naringenin itself, Table S2: The Hammett sigma constants (σm and σp), the field/inductive parameter (F) and resonance parameter (R) of the substituents, Table S3: ΔIP in kJ/mol obtained by the difference between IP of the substituted naringenin and naringenin itself, Table S4: ΔPA in kJ/mol obtained by the difference between PA of the substituted naringenin and naringenin itself.

Author Contributions

Conceptualization, Y.-Z.Z. and G.D.; Methodology, Y.-Z.Z.; Software, G.D.; Validation, Y.-Z.Z., G.D. and R.G.; Formal Analysis, Y.-Z.Z.; Investigation, G.D.; Resources, G.D. and R.G.; Data Curation, R.G.; Writing-Original Draft Preparation, Y.-Z.Z. and G.D.; Writing-Review & Editing, G.D., R.G., D.-F.C. and Z.-M.F.; Visualization, Y.-Z.Z.; Supervision, D.-F.C.; Project Administration, Y.-Z.Z.; Funding Acquisition, Y.Z.Z. and D.-F.C.

Funding

This work was funded by the Earmarked Fund for China Agriculture Research System (CARS-44-KXJ7) and the Fujian Agriculture and Forestry University Foundation for excellent youth teachers (xjq201715).

Conflicts of Interest

The authors declare no conflict of interest.

Abbreviations

DFTDensity Functional Theory
HATHydrogen Atom Transfer
SPLETSequential Proton Loss Electron Transfer
SET-PTSingle Electron Transfer Followed by Proton Transfer
ArOHFlavonoid
BDEBond Dissociation Enthalpy
IPIonization Potential
PDEProton Dissociation Enthalpy
PAProton Affinity
ETEElectron Transfer Enthalpy
ΔrGGibbs Free Energy
σHammett sigma constants
Ffield/inductive effect
Rresonance effect
PPearson correlation coefficient

References

  1. Dizdaroglu, M.; Jaruga, P.; Birincioglu, M.; Rodriguez, H. Free radical-induced damage to DNA: Mechanisms and measurement. Free Radic. Biol. Med. 2002, 32, 1102–1115. [Google Scholar] [CrossRef]
  2. Fang, Y.Z.; Yang, S.; Wu, G. Free radicals, antioxidants, and nutrition. Nutrition 2002, 18, 872–879. [Google Scholar] [CrossRef]
  3. Maritim, A.C.; Sanders, R.A.; Watkins, J.B. Diabetes, oxidative stress, and antioxidants: A review. J. Biochem. Mol. Toxicol. 2003, 17, 24–37. [Google Scholar] [CrossRef] [PubMed]
  4. Pryor, W.A. Vitamin E and heart disease: Basic science to clinical intervention trials. Free Radic. Biol. Med. 2000, 28, 141–164. [Google Scholar] [CrossRef]
  5. Hahn, M.; Baierle, M.; Charão, M.F.; Bubols, G.B.; Gravina, F.S.; Zielinsky, P.; Arbo, M.D.; Cristina Garcia, S. Polyphenol-rich food general and on pregnancy effects: A review. Drug Chem. Toxicol. 2017, 40, 368–374. [Google Scholar] [CrossRef] [PubMed]
  6. Ross, J.A.; Kasum, C.M. Dietary flavonoids: Bioavailability, metabolic effects, and safety. Annu. Rev. Nutr. 2002, 22, 19–34. [Google Scholar] [CrossRef] [PubMed]
  7. Muhammad, D.; Hubert, J.; Lalun, N.; Renault, J.H.; Bobichon, H.; Nour, M.; Voutquenne-Nazabadioko, L. Isolation of flavonoids and triterpenoids from the fruits of Alphitonia neocaledonica and evaluation of their anti-oxidant, anti-tyrosinase and cytotoxic activities. Phytochem. Anal. 2015, 26, 137–144. [Google Scholar] [CrossRef] [PubMed]
  8. Trimech, I.; Weiss, E.K.; Chedea, V.S.; Marin, D.; Detsi, A.; Ioannou, E.; Roussis, V.; Kefalas, P. Evaluation of anti-oxidant and acetylcholinesterase activity and identification of polyphenolics of the invasive weed Dittrichia viscosa. Phytochem. Anal. 2015, 25, 421–428. [Google Scholar] [CrossRef] [PubMed]
  9. Tapas, A.R.; Sakarkar, D.M.; Kakde, R.B. Flavonoids as nutraceuticals: A review. Trop. J. Pharm. Res. 2008, 7, 1089–1099. [Google Scholar] [CrossRef]
  10. Gülçin, I. Antioxidant activity of food constituents: An overview. Arch. Toxicol. 2012, 86, 345–391. [Google Scholar] [CrossRef]
  11. Akhlaghi, M.; Bandy, B. Mechanisms of flavonoid protection against myocardial ischemia–reperfusion injury. J. Mol. Cell. Cardiol. 2009, 46, 309–317. [Google Scholar] [CrossRef] [PubMed]
  12. Olsen, H.T.; Stafford, G.I.; van Staden, J.; Christensen, S.B.; Jäger, A.K. Isolation of the MAOinhibitor naringenin from Mentha aquatica L. J. Ethnopharmacol. 2008, 117, 500–502. [Google Scholar] [CrossRef] [PubMed]
  13. Zbarsky, V.; Datla, K.P.; Parkar, S.; Rai, D.K.; Aruoma, O.I.; Dexter, D.T. Neuroprotective properties of the natural phenolic antioxidants curcumin and naringenin but not quercetin and fisetin in a 6-OHDA model of Parkinson’s disease. Free Radic. Res. 2005, 39, 1119–1125. [Google Scholar] [CrossRef] [PubMed]
  14. Youdim, M.B.; Bakhle, Y.S. Monoamine oxidase: Isoforms and inhibitors in Parkinson’s disease and depressive illness. Br. J. Pharmacol. 2006, 147, S287–296. [Google Scholar] [CrossRef] [PubMed]
  15. Chen, W.; Song, J.; Guo, P.; Cao, W.; Bian, J. Exploring a possible way to synthesize novel better antioxidants based on vitamin E: A DFT study. Bioorg. Med. Chem. Lett. 2006, 16, 5874–5877. [Google Scholar] [CrossRef] [PubMed]
  16. Klein, E.; Lukeš, V.; Ilčin, M. DFT/B3LYP study of tocopherols and chromans antioxidant action energetics. Chem. Phys. 2007, 336, 51–57. [Google Scholar] [CrossRef]
  17. Krygowski, T.M.; Stȩpień, B.T. Sigma-and pi-electron delocalization: Focus on substituent effects. Chem. Rev. 2005, 105, 3482–3512. [Google Scholar] [CrossRef]
  18. Hansch, C.; Leo, A.; Taft, R.W. A survey of Hammett substituent constants and resonance and field parameters. Chem. Rev. 1991, 91, 165–195. [Google Scholar] [CrossRef]
  19. Štefanič-Petek, A.; Krbavčič, A.; Šolmajer, T. QSAR of flavonoids: 4. Differential inhibition of aldose reductase and p56 lck protein tyrosine kinase. Croat. Chem. Acta 2002, 75, 517–529. [Google Scholar]
  20. Aparicio, S. A systematic computational study on flavonoids. Int. J. Mol. Sci. 2010, 11, 2017–2038. [Google Scholar] [CrossRef]
  21. Trouillas, P.; Fagnère, C.; Lazzaroni, R.; Calliste, C.; Marfak, A.; Duroux, J.L. A theoretical study of the conformational behavior and electronic structure of taxifolin correlated with the free radical-scavenging activity. Food Chem. 2004, 88, 571–582. [Google Scholar] [CrossRef]
  22. Leopoldini, M.; Russo, N.; Toscano, M. The molecular basis of working mechanisms of natural polyphenolic antioxidants. Food Chem. 2011, 125, 288–306. [Google Scholar] [CrossRef]
  23. Stepanić, V.; Trošelj, K.G.; Lučić, B.; Marković, Z.; Amić, D. Bond dissociation free energy as a general parameter for flavonoid radical scavenging activity. Food Chem. 2013, 141, 1562–1570. [Google Scholar] [CrossRef]
  24. Wright, J.S.; Johnson, E.R.; DiLabio, G.A. 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. [Google Scholar] [CrossRef] [PubMed]
  25. Nenadis, N.; Siskos, D. Radical scavenging activity characterization of synthetic isochroman-derivatives of hydroxytyrosol: A gas-phase DFT approach. Food Res. Int. 2015, 76, 506–510. [Google Scholar] [CrossRef] [PubMed]
  26. Vagánek, A.; Rimarčík, J.; Dropková, K.; Lengyel, J.; Klein, E. Reaction enthalpies of O–H bonds splitting-off in flavonoids: The role of non-polar and polar solvent. Comput. Theor. Chem. 2014, 1050, 31–38. [Google Scholar] [CrossRef]
  27. Vargas-Sánchez, R.D.; Mendoza-Wilson, A.M.; Balandrán-Quintana, R.R.; Torrescano-Urrutia, G.R.; Sánchez-Escalante, A. Study of the molecular structure and chemical reactivity of pinocembrin by DFT calculations. Comput. Theor. Chem. 2015, 1058, 21–27. [Google Scholar] [CrossRef]
  28. Wang, G.; Xue, Y.; An, L.; Zheng, Y.; Dou, Y.; Zhang, L.; Liu, Y. Theoretical study on the structural and antioxidant properties of some recently synthesised 2, 4, 5-trimethoxy chalcones. Food Chem. 2015, 171, 89–97. [Google Scholar] [CrossRef]
  29. 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. [Google Scholar] [CrossRef]
  30. Zheng, Y.Z.; Deng, G.; Chen, D.F.; Liang, Q.; Guo, R.; Fu, Z.M. Theoretical studies on the antioxidant activity of pinobanksin and its ester derivatives: Effects of the chain length and solvent. Food Chem. 2018, 240, 323–329. [Google Scholar] [CrossRef]
  31. Leopoldini, M.; Marino, T.; Russo, N.; Toscano, M. Antioxidant properties of phenolic compounds: H-atom versus electron transfer mechanism. J. Phys. Chem. A 2004, 108, 4916–4922. [Google Scholar] [CrossRef]
  32. Trouillas, P.; Marsal, P.; Siri, D.; Lazzaroni, R.; Duroux, J.L. A DFT study of the reactivity of OH groups in quercetin and taxifolin antioxidants: The specificity of the 3-OH site. Food Chem. 2006, 97, 679–688. [Google Scholar] [CrossRef]
  33. Leopoldini, M.; Rondinelli, F.; Russo, N.; Toscano, M. Pyranoanthocyanins: A theoretical investigation on their antioxidant activity. J. Agric. Food Chem. 2010, 58, 8862–8871. [Google Scholar] [CrossRef] [PubMed]
  34. Marković, Z.; Milenković, D.; Đorović, J.; Marković, J.M.D.; Stepanić, V.; Lučić, B.; Amić, D. PM6 and DFT study of free radical scavenging activity of morin. Food Chem. 2012, 134, 1754–1760. [Google Scholar] [CrossRef] [PubMed]
  35. Lengyel, J.; Rimarčík, J.; Vagánek, A.; Klein, E. On the radical scavenging activity of isoflavones: Thermodynamics of O–H bond cleavage. Phys. Chem. Chem. Phys. 2013, 15, 10895–10903. [Google Scholar] [CrossRef] [PubMed]
  36. Anouar, E.; Košinová, P.; Kozlowski, D.; Mokrini, R.; Duroux, J.L.; Trouillas, P. New aspects of the antioxidant properties of phenolic acids: A combined theoretical and experimental approach. Phys. Chem. Chem. Phys. 2009, 11, 7659–7668. [Google Scholar] [CrossRef]
  37. Silva, M.M.; Santos, M.R.; Caroço, G.; Rocha, R.; Justino, G.; Mira, L. Structure-antioxidant activity relationships of flavonoids: A re-examination. Free Radic. Res. 2002, 36, 1219–1227. [Google Scholar] [CrossRef]
  38. Zhang, H.Y.; Sun, Y.M.; Wang, X.L. Electronic effects on O–H proton dissociation energies of phenolic cation radicals: A DFT Study. J. Org. Chem. 2002, 67, 2709–2712. [Google Scholar] [CrossRef]
  39. Zhang, H.Y.; Sun, Y.M.; Chen, D.Z. O–H bond dissociation energies of phenolic compounds are determined by field/inductive effect or resonance effect? A DFT study and its implication. Mol. Inform. 2001, 20, 148–152. [Google Scholar] [CrossRef]
  40. Frisch, M.J.; Trucks, G.W.; Schlegel, H.B.; Scuseria, G.E.; Robb, M.A.; Cheeseman, J.R.; Scalmani, G.; Barone, V.; Mennucci, B.; Petersson, G.A.; Nakatsuji, H. Gaussian 09 (Revision B.01); Gaussian Inc.: Wallingford, CT, USA, 2013. [Google Scholar]
  41. Bartmess, J.E. Thermodynamics of the electron and the proton. J. Phys. Chem. 1994, 98, 6420–6424. [Google Scholar] [CrossRef]
  42. Rimarčík, 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 pphenylenediamine and tetracyano- p-phenylenediamine. J. Mol. Struct. Theochem. 2010, 925, 25–30. [Google Scholar] [CrossRef]
  43. Parker, V.D. 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. [Google Scholar] [CrossRef]
  44. Rodriguez-Lujan, I.; Huerta, R.; Elkan, C.; Cruz, C.S. Quadratic programming feature selection. J. Mach. Learn. Res. 2010, 11, 1491–1516. [Google Scholar]
Figure 1. The flavonoid basic structure and the chemical structure of naringenin.
Figure 1. The flavonoid basic structure and the chemical structure of naringenin.
Ijms 20 01450 g001
Figure 2. The Pearson correlation coefficients between the Hammett sigma constants (σm and σp)/F/R and the BDE/ ionization potential (IP)/ ionization potential (PA) in the gas (A), benzene (B) and water (C) phases.
Figure 2. The Pearson correlation coefficients between the Hammett sigma constants (σm and σp)/F/R and the BDE/ ionization potential (IP)/ ionization potential (PA) in the gas (A), benzene (B) and water (C) phases.
Ijms 20 01450 g002
Figure 3. Dependence of BDE (A), IP (B) and PA (C) on Hammett sigma constants for naringenin in gas, benzene and water phases.
Figure 3. Dependence of BDE (A), IP (B) and PA (C) on Hammett sigma constants for naringenin in gas, benzene and water phases.
Ijms 20 01450 g003
Table 1. O–H bond dissociation enthalpy (BDE) in kJ/mol obtained by the M062X/6-311+G ** method. The data in the form of the bold and underline represent the lowest BDEs in the molecule.
Table 1. O–H bond dissociation enthalpy (BDE) in kJ/mol obtained by the M062X/6-311+G ** method. The data in the form of the bold and underline represent the lowest BDEs in the molecule.
SubstituentsGasBenzeneWater
4′−OH5−OH7−OH4′−OH5−OH7−OH4′−OH5−OH7−OH
H366.5428391.1368.6423.2395.2362.2391.2393.4
NH2364.2425.3389.7366.7421.0393.6360.7389.9391.0
OH365.4426.1390.1367.9421.5393.8361.5390.1391.6
OMe365.2425.8390.4367.7421.3394.2362.4390.2392.1
Me367.0427.2390.9369.2422.7394.5362.1390.6392.5
F368.6429.1391.5370.8423.8396.1364.6391.6394.8
Cl369.3430.1393.9372.0423.9396.6366.7392.1395.5
CHO369.9430.9394.3372.3424.1398.4367.0392.7396.0
CF3371.0431.2394.8372.6424.7399.3367.1393.6396.9
CN372.1431.9396.2374.3425.5400.9367.3393.9398.5
NO2373.7432.4397.3375.5426.1402.3368.4394.3400.5
The bold and underline represent the lowest BDEs in the molecule.
Table 2. Ionization potential (IP) in kJ/mol obtained by the M062X/6-311+G ** method.
Table 2. Ionization potential (IP) in kJ/mol obtained by the M062X/6-311+G ** method.
SubstituentsGasBenzeneWater
H791.0678.2586.6
NH2786.9674.6580.6
OH797.8679.6587.6
OMe795.7680.8584.7
Me787.7676.8585.6
F813.7697.1595.5
Cl809.3695.6596.4
CHO810.4696.2599.2
CF3813.8698.0597.7
CN826.5707.1599.4
NO2832.4713.6605.4
Table 3. Proton affinity (PA) in kJ/mol obtained by the M062X/6-311+G** method. The data in the form of the bold and underline represent the lowest PAs in the molecule.
Table 3. Proton affinity (PA) in kJ/mol obtained by the M062X/6-311+G** method. The data in the form of the bold and underline represent the lowest PAs in the molecule.
SubstituentsGasBenzeneWater
4′−OH5−OH7−OH4′−OH5−OH7−OH4′−OH5−OH7−OH
H1406.31445.31373.3439.2468.8407.4162.3160.2138.4
NH21409.51447.71378.8442.2472.0409.8164.2161.9141.0
OH1405.71442.61363.1437.9464.1404.0162.0153.8137.9
OMe1405.51435.41365.8437.8462.3403.0161.9157.5137.6
Me1406.91447.31376.2440.4471.2408.2162.6161.6138.8
F1404.01423.61355.0437.8450.1393.1161.5150.0132.7
Cl1404.61424.01354.9437.4451.2393.2161.2152.1132.8
CHO1402.91418.91351.0436.4450.0392.8161.5152.7134.3
CF31401.11419.01352.3436.2448.5392.3160.2151.8133.4
CN1399.61407.01340.6433.6439.4384.3159.3148.9131.3
NO21397.81399.61334.9431.8433.0379.6158.6142.2127.7
The bold and underline represent the lowest PAs in the molecule.

Share and Cite

MDPI and ACS Style

Zheng, Y.-Z.; Deng, G.; Guo, R.; Chen, D.-F.; Fu, Z.-M. DFT Studies on the Antioxidant Activity of Naringenin and Its Derivatives: Effects of the Substituents at C3. Int. J. Mol. Sci. 2019, 20, 1450. https://doi.org/10.3390/ijms20061450

AMA Style

Zheng Y-Z, Deng G, Guo R, Chen D-F, Fu Z-M. DFT Studies on the Antioxidant Activity of Naringenin and Its Derivatives: Effects of the Substituents at C3. International Journal of Molecular Sciences. 2019; 20(6):1450. https://doi.org/10.3390/ijms20061450

Chicago/Turabian Style

Zheng, Yan-Zhen, Geng Deng, Rui Guo, Da-Fu Chen, and Zhong-Min Fu. 2019. "DFT Studies on the Antioxidant Activity of Naringenin and Its Derivatives: Effects of the Substituents at C3" International Journal of Molecular Sciences 20, no. 6: 1450. https://doi.org/10.3390/ijms20061450

APA Style

Zheng, Y. -Z., Deng, G., Guo, R., Chen, D. -F., & Fu, Z. -M. (2019). DFT Studies on the Antioxidant Activity of Naringenin and Its Derivatives: Effects of the Substituents at C3. International Journal of Molecular Sciences, 20(6), 1450. https://doi.org/10.3390/ijms20061450

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