Catechin and Epicatechin. What’s the More Reactive? ()
1. Introduction
Catechin and epicatechin are two isomeric flavanols. They differ from the R and S configuration linked to the two asymmetric carbons of 2-phenyl-3-chromanol (Figure 1) [1] [2] [3] [4]. In Figure 1, each ball represents an atom. Red is oxygen. A white ball illustrates the carbon. A small one describes hydrogen. Catechin and epicatechin are widespread in the plant kingdom through fruits such as grapes, apples and too many manufactured products such as wine, chocolate, and tea [2] [5]. With the consumption of these foods, humans benefit from antioxidant properties [6] [7] [8] [9] [10]. These delay cell damage and combat certain chronic diseases including cardiovascular and neurodegenerative diseases ones [6] [7] [11] [12] [13] [14]. These molecules also possess anti-cancer, anti-inflammatory and anti-viral properties [15] [16] [17] [18]. The wide spectrum of polyphenols’ biological activities gives them a great importance in therapeutic processes. Besides, only 2% of these alcohols consumed reach the plasma. This underperformance thus represents a problem of bioavailability [19]. Increasing their proportion in this organism would make them much more efficacious. For better control and improvement of polyphenols’ therapeutic effects, it’s advisable to discover ways to grow their concentration in plasma. The development of medicines incorporating polyphenols would be an asset in providing effective remedies for these illnesses that ruin human populations. Tests to manufacture catechin-specific anti-inflammatory drugs weren’t as conclusive as expected. The reasons are related to the failures of the very limited clinical tests because of their stability, their short half-life in the plasma and their low bioavailability [20] [21].
The properties exhibited by catechin and epicatechin emanate from their many hydroxyl groups [22]. The aromatic rings and their OH associated makes it possible to build several types of interactions. These permit polyphenols to scavenge hydroxyl radicals, inhibit reactions such as lipid peroxidation and
Figure 1. Structure of Catechin and Epicatechin.
prevent oxidation reactions. This process happens in a state of transition [22]; the latter promotes the occupation of their specific molecular orbitals [23]. These foster some CHO bond at three centres. However, despite the large number of studies on the biological mechanisms explaining the beneficial effects of polyphenols on human health, few quantitative data are available on the polyphenol’s reactivity. This research aims to fill this gap through the following question:
Which isomer is the most reactive between catechin and epicatechin?
Its response fulfills the gap in numerical data relating to the reactivity of the two compounds. But it’s difficult to identify all the factors involved in the instability of the two isomers. To correct this shortcoming, the work wants to address another question:
What are the most receptive sites of isomers?
The research plans to suggest avenues that will contribute to a better integration of molecular interactions and biological activities. It assimilates catechin and epicatechin to a polyphenol model. It evaluates their reactivity using physical and chemical quantities. These include the energy, infrared (IR) and ultraviolet (UV) spectroscopic parameters. They incorporate electronic properties such as the HOMO (High Occupied Molecular Orbital)-LUMO (Lowest Unoccupied Molecular Orbital) energy gap, chemical hardness and chemical potential. They comprise molecular electrostatic potential (MEP) and Fukui indices. These analyze the local and global reactivity. Furthermore, the acidity potential (pKa) specifies that of each hydroxyl group. Natural Bond Orbital (NBO) analysis allows following charge transfer and intramolecular interactions in catechin and epicatechin. This article discusses the calculations relating to these parameters. Previously, it presented the method and the material used to obtain them.
2. Materials and Methods
This part includes the descriptors of global and local reactivity. It explains those of acidity. It describes how to assess load transfers using NBO analysis. For the moment, it specifies the method of calculations. The bond lengths, the bond angles and the standard dihedral angle constitute the initial geometrical parameters. The geometry optimizations and the frequency calculations were carried out with the Density Functional Theory (DFT) method at the B3LYP/6-31 + G(d, p) level of theory [24] in gas and aqueous phases.
2.1. Calculation Method
The CPCM (Conductor-like Polarizable Continuum Model) calculation simulates the physiological environment. All calculations were performed with the GAUSSIAN 09 software [25]. Frequency calculations verify local and global minima. TD-DFT method permits to analyze the electron transitions in the ultraviolet (UV) range. The research determines the energies of the frontier molecular orbitals including the highest occupied one (HOMO) and the lowest unoccupied one (LUMO). It uses NBO calculation [26] examining the NBO. Gaussian’s calculations provide access to the frontier orbitals of a molecule. Their HOMO and LUMO contribute interpreting the reactivity of a molecule [27] [28].
2.2. Global Reactivity Descriptors
The energies of these two frontier orbitals lead to the determination of several reactivity parameters. According to Koopmans theorem [29], the ionization energy I and the electronic affinity A are directly related to the energy of the HOMO and the LUMO respectively.
(1)
(2)
The chemical electron potential μ and chemical hardness η are defined in terms of ionization energy and electron affinity [30]:
(3)
(4)
Also, global index of electrophilicity ω, introduced by Parr [31], is defined by the following formula:
(5)
Local reactivity index is also accessible by calculation. While the global reactivity parameters evaluate that of a molecule, Fukui [32] introduces indices (Fukui parameters) to describe the reactivity of each atom.
2.3. Local Reactivity Descriptors
These parameters denote either a nucleophilic attack (
), an electrophilic attack (
), or a radical attack (
). The following equations [33] help to determine them:
(6)
(7)
(8)
With
: the electron population of atom k in the neutral molecule.
: the electron population of atom k in the cationic molecule.
: the electron population of atom k in the anionic molecule.
The acidity represents an important parameter to integrate the intermolecular interactions. The calculations give access to the pKa.
2.4. Acidity
Catechin and epicatechin represent polyphenols with both five hydroxyl groups. They’re therefore likely to exchange one or more protons with other molecules, either during hydrogen bonding interactions or “antiradical” mechanisms [34]. The pKa of each hydroxyl group evaluates the proton transfer reactions between catechin and epicatechin and additional molecules. It helps classifying their acidity. The thermodynamic cycle below allows calculating it by considering the effect of solvation due to the presence of H2O [35].
The free energy of solvation
becomes:
(9)
The value of
equals −110.2 kcal/mol. Pliego and Riveros [36] determined it experimentally.
(10)
where [H2O] = 55.5 M. However, since the free energy of the hydronium ion constitutes a source of error. According to Pliego [35], the corrected formula for pKa is written as follows. On the other hand, this calculation of Pka leads to the NBO analysis.
(11)
2.5. NBO Analysis
NBO analysis [26] evaluates the charge transfer properties in molecules. It provides a deep understanding of intramolecular and intermolecular orbital interactions between occupied NBO donors and empty NBO acceptors [37] [38]. The charge transfer involves delocalization of electrons from an electron donor site Lewis-type orbital: lone pairs (n), natural bond (σ and π) to an electron acceptor one (anti-Lewis orbital: σ* and π*). This electron delocalization is with a decrease in the donor electronic density and an increase in the acceptor one. For each donor NBO (i) and acceptor NBO (j), the electron delocalization i → j is evaluated by the stabilization energy E(2) using second-order perturbation theory [39]. This latter is expressed as:
(12)
where
denotes the electronic density in the donor orbital,
designs a non-diagonal element of the Fock matrix,
and
correspond to the energies of the occupied i and empty j orbitals respectively. The higher
, the more redistribution of electrons between donor and acceptor confers stability. This affirmation introduces the presentation of the results and their analysis.
3. Results and Discussion
This section discusses those related to indicators of global and local responsiveness. Its analysis concerns acidity and NBO comprise results of calculations associated with spectral quantities. Catechin or epicatechin is composed of 35 atoms. Vibration mode analysis suggests 99 normal vibrations for both molecules.
3.1. Spectroscopic Parameters
These vibrations are divided into 34 stretching, 33 deformations and 32 torsions. The 34 stretching is repartitioned into five O–H, nine C–H, 13 C–C and seven C–O. This agrees with the structure of the two molecules (Figure 1). The analysis of the potential energy distribution (PED) by VEDA can efficiently assign the calculated normal modes of vibration. Table 1 presents the O–H elongation vibration frequencies of catechin and epicatechin. It indicates that the catechin and epicatechin O–H elongation frequencies vary from 3788 to 3845 cm−1. These agree with the experimental values of the O–H bonds of aromatic rings [40].
The first three O–H elongation frequencies stay at 3845 cm−1, 3833 cm−1, and 3830 cm−1 correspond to
,
, and
, respectively, of the two molecules. Their intensities are approximately the same. The ranking of the –H bonds by their elongation frequencies is:
for epicatechin and
for catechin.
changes from 3803 cm−1 in catechin to 3788 cm−1 in epicatechin. Their intensities of this vibration also change. They increase from 28.44 km/mole in the catechin to 61.35 km/mole. O30–H31 is affected by the asymmetry of carbon C12.
3.2. UV-Visible Spectral Analysis
The TD-DFT method at the level of the theory B3LYP/6-31 + G (d, p), in the gas phase evaluates the electron transitions of catechin or epicatechin in the ultraviolet and visible regions. It gives molecule excitation energy. This latter quantity permits comparing the reactivity of two molecules. Its higher value corresponds to the more stable compound. Table 2 shows the results of the TD-DFT calculations. It displays the first three excited states of both molecules as presented in Figure 2.
Table 1. Vibrational frequencies (cm−1), IR intensity (km/mol) and vibrational assignation (PED of Catechin and Epicatechin.
Figure 2. UV-visible spectra of catechin and epicatechin.
Table 2. Main transitions calculated for Catechin and Epicatechin.
For epicatechin, the first excited state wavelength is 262 nm. This is the most intense excitation. Its main component refers to the HOMO → LUMO transition with a weight of 61%. Its excitation energy is 4.7355 eV. The second one relates to the band with a wavelength equals 260.20 nm. It associates with the HOMO − 2 → LUMO + 1 transition. Its excitation energy is 4.7453 eV. The latter is greater than that of the first excited state. This result means that epicatechin is more reactive than catechin.
3.3. Global Reactivity Descriptors
To access the chemical reactivity of catechin or epicatechin, it’s necessary to determine their molecular border orbital HOMO and LUMO [27]. More, the first orbitals can explain the physical and chemical properties of molecules [23]. Here, they help to provide insight into intramolecular charge transfers; an electron-rich HOMO orbital acts as a donor; an electron-poor LUMO functions as an acceptor. Figure 3 and Figure 4 illustrate the HOMO − 3 to LUMO + 3 boundary orbitals of catechin and epicatechin. These molecular orbitals have a π character.
For HOMO, the π character is distributed over all the catechin and epicatechin. This means that the delocalization of electrons occurs through the molecules of catechin and epicatechin. Unlike LUMO, the π character of these molecules is concentrated on the B-cycle with a residual contribution of its oxygen atoms. The global reactivity indices are presented in Table 3. The smaller the energy gap (∆E) between HOMO and LUMO is, the more reactive is the molecule [27]. The energy gap of epicatechin (5.347 eV) is the lowest. It indicates that this molecule is most receptive. The global electrophilicity indices ω confirms this observation. Its epicatechin’s value is slightly higher than those of catechin. The epicatechin is more able to acquire electronic charges. The chemical hardness of epicatechin η (2.674 eV) is lower than that of catechin (2.720 eV). The epicatechin is softer than catechin. The Fukui indices
,
, and
represents the local reactivity descriptors.
Figure 3. Frontier orbitals of the catechin’s main calculated transitions.
Figure 4. Frontier orbitals of the epicatechin’s main calculated transitions.
Table 3. Global reactivity parameters of catechin and epicatechin.
3.4. Local Reactivity Descriptors
Fukui indices are used to assess the reactivity of atoms or functional groups (nucleophilic or electrophilic attack site) of a molecule [33]. This work is based on a natural population analysis (NPA) in gas and aqueous phase to estimate the electron population [41]. This leads to the Fukui indices, reported in Tables 4-6. The Molecular Electrostatic Potential (MEP) map displays the electron density of a molecule as a function of colour. It constitutes a descriptor of the local reactivity for a molecule. Colours indicate the electron density of areas of the molecule [42]. A red region illustrates an electron-rich site. A blue zone describes an electron-poor site. A green area is a neutral site. Figure 5 gives the plots of the catechin and epicatechin electrostatic potential.
These tables data show that the highest
values in both the aqueous and gas phases link to the hydrogen atoms H27, H29, H31, H33, and H35 of catechin and epicatechin. These constitute the preferential sites of nucleophilic attack. On the other hand, some are favoured on the oxygen atoms and the carbon atoms of the aromatic rings. The properties are preserved when the molecule passing from the gas phase of the aqueous phase.
Table 4. Global reactivity parameters of catechin and epicatechin (atoms 1 to 12).
Table 5. Global reactivity parameters of catechin and epicatechin (atoms 13 to 24).
catechinepicatechin
Figure 5. Molecular electrostatic potential of Catechin and Epicatechin.
Table 6. Global reactivity parameters of catechin and epicatechin (atoms 25 to 35).
This map indicates that the hydroxyl groups O28–H29, O32–H33, and O34–H35, are more depleted in electrons than those of O30–H31, O26–H27 and of the different Ci–Hj bonds. These hydrogen atoms of the Oi–Hj hydroxyl groups and Ci–Hj bonds of catechin and epicatechin are sites conducive to nucleophilic attack on the following decreasing order of reactivity : O28–H29 > O32–H33 > O34–H35 > O26–H27 > O30–H31 > Ci–Hj. The regions of the oxygen atoms and the aromatic rings are, on the other hand, favourable sites for electrophilic attacks.
3.5. Acidity
Table 7 shows the theoretical values of pKa for each hydroxyl group of the catechin and its isomer. This statistic varies from 6.77 to 21.77 for catechin and from 6.77 to 21.70 for epicatechin. The pka of the hydroxyl group’s O32–H33, O34–H35 and O26–H27 are respectively 9.92, 10.47, 11.47 for catechin. For epicatechin, they’re 10.56, 11.21, 1171. These values are close to the pka experimental value of the phenol’s hydroxyl; the latter is equal to 9.95 [43]. The pKa values of the hydroxyl group O28–H29 are 7.53 in catechin and 6.77 in epicatechin. These values remain inferior to 9.85. Therefore, the hydroxyl group O28–H29 of catechin or epicatechin is more acidic than the phenol’s hydroxyl group and those of catechin or epicatechin. The difference between the pKa of the hydroxyl group O30–H31 and that of the phenol is great. It shows that the hydroxyl group isn’t on an aromatic ring. Its acidic character is very weak compared to the other hydroxyl groups of catechin and epicatechin. The hydroxyl group acidity becomes: O28–H29 > O32–H33 > O34–H35 > O26–H27 > O30–H31. The hydroxyl groups O28–H29, O32–H33, O34–H35 are much more acidic. They’re more favourable to the anti-radical processes than the hydroxyl groups O26–H27 or O30–H31.
3.6. NBO Analysis
Table 8 and Table 9 collate the second-order perturbation energies E(2), the electron density (ED), the energy difference E(j) – E(i) of the donor NBO (i) and of the acceptor NBO (j) related to catechin and epicatechin. They show their elements of the Fock matrix F(i, j). These results indicate that for catechin and epicatechin, the main interactions are of two types. They’re intramolecular interactions of the types
and
where
,
,
Table 7. pKa values of catechin and epicatechin.
Table 8. The catechin’s NBO parameters.
Table 9. The epicatechin’s NBO parameters.
are: Lewis π orbital of the
double bond, the lone pair 2 of the Oε atoms, anti-Lewis π orbital of the
double bond. α, β, γ, δ, ε are atomic numbers. In the aromatic ring A, the stabilization energies of the
,
and
interactions are 25, 81, 25.41, 24.64 kcal/mol respectively for catechin and 25.74, 25.95, 24.71 kcal/mol for epicatechin.
In aromatic ring B, for the same intramolecular interaction
, the maximum stabilization energy is equal to 20.33 kcal/mol for catechin and 20.25 kcal/mol for epicatechin. Thus, the aromatic ring B is less stable than of the A one. Regarding the
interactions, the lowest stabilization energy is obtained with
. For catechin, it’s worthy 24.07 kcal/mol and 24.13 kcal/mol for epicatechin. The electronic density (ED) of the Lewis orbital
and
undergo a decreasing (ED < 2e) while those of the anti-Lewis orbital
increase (ED > 0).
These results show that intramolecular interactions are associated with the delocalization of π-elections from Lewis’s orbitals to anti Lewis’s orbitals. They also indicate that ring A is more stable than ring B in both molecules. The ring B is therefore more reactive than the ring A. Its reactivity influences that of the hydroxyl group O28–H29.
transition is the most active of the
transitions. Its delocalization energy (24.07 kcal/mol) is the lowest of the four catechin.
This finding is too valid for epicatechin. More, it contrasts the reactivity priority order of the hydroxyls. Its strong reactivity also enhances that of the hydroxyl O28–H29. Under these conditions, any electrophilic or nucleophilic attack takes place primarily and respectively on O28 and O29. In the event of unavailability, hydroxyl O26–H27 becomes priority; it attaches to ring B. The other two hydroxyls of ring A (O32–H33 = O34–H35) follow in order of reactivity. Hydroxyl O30–H31 loops the latter from the ring C. This aspect of the research leads to this article’s conclusion.
4. Conclusions
The research plans to compare the reactivity of two isomers, catechin and epicatechin, using the resources of theoretical chemistry. It’s carried out at the TD-DFT/B3LYP/6-31 + G(d, p) level. It’s interesting in the parameters of global and local responsiveness. The MEP and the Fukui indices make it possible to specify the latter. It also harnesses the distribution of potential energy through VEDA. It contrasts the acidity. It ends with an NBO analysis. These statistics probe that epicatechin remains the more reactive compound of the two. More, the Fukui indices and the analysis of vibration energies establish that hydroxyl group reactivity of two isomers varies in the following order: O28–H29 > O32–H33 = O34–H35 > O30–H31 > O26–H27. In other words, the O28–H29 site remains the least stable of all these OH. The molecular electrostatic potential confirms these results.
The acidity indicated that the hydroxyls O28–H29, O32–H33, O34–H35 remain the most acid in the sense of Brönsted. In other words, the hydroxyl O28–H29 the most favourable site to participate remains in anti-free radical processes. NBO analysis detected that the stabilization energies of the two isomers equal 25 81, 25.41, 24.64 kcal/mol. Besides, it presents the transition
,
, and
. Their energies remain equal to 25.74, 25.95, 24.71 kcal/mol for epicatechin. It appears that ring B is more reactive than A. Its reactivity reinforces that of O28–H29. This hydroxyl becomes the most receptive site of the catechin or epicatechin. As a result, the hydroxyls of ring A (O32–H33 = O34–H35) and of ring B (O26–H27) follow in decreasing order to their reactivity. This research furnishes the team’s next agenda.