# Effect of Constituent Units, Type of Interflavan Bond, and Conformation on the Antioxidant Properties of Procyanidin Dimers: A Computational Outlook

Authors: Ana María Mendoza-Wilson; René Renato Balandrán-Quintana
Journal: Journal of Chemistry
 (2017)

---

## Abstract

Procyanidin (PC) dimers are powerful antioxidants, abundant in plant tissues, and also bioavailable. However, the role of the molecular structure of PCs on their antioxidant properties is still a controversial and not fully understood issue that needs to be addressed in a more specific way. The objective of this study was to analyze the effect of the constituent units, type of interflavan bond, and conformation on the antioxidant properties of PC dimers including PB3, PB4, PB5, PB6, PB7, and PB8, using the density functional theory (DFT) computational method. The analysis was performed in function of parameters that allow determining the ability of the molecules to transfer or to capture electrons, among which the chemical potential, bond dissociation enthalpy (BDE), gap energy, Fukui indices, and charge distribution of HOMO-LUMO orbitals. The factors that showed the most notable effects on the antioxidant properties of the PC dimers were the type of interflavan bond and the conformation. The antioxidant ability of the dimers PB3 and PB4 containing the interflavan bond C4–C8, in theirCompact conformation, was very similar to each other but greater than those of dimers PB5, PB6, PB7, and PB8 containing the C4–C6 interflavan bond. PB8 showed the lowest antioxidant ability.


---

## Body

## 1. Introduction

Procyanidins (PCs) are recognized as potent antioxidants due to their ability to inhibit the detrimental effects of oxidative enzymes, metals, and free radicals. These compounds are considered bioresources with possible applications in food, medicine, and pharmacology [1]. However, the relationship between the molecular structure of PCs and their antioxidant properties is still a controversial and not fully understood issue that needs to be addressed in a more specific way. This is in part due to the great structural diversity of the PCs regarding the degree of polymerization, type of constitutive units, type of interflavan bond, and conformations.The PCs are oligomeric or polymeric derivatives of flavan-3-ols, which are constituted by (+)-catechin and/or (−)-epicatechin units, mainly linked through C4–C8 or C4–C6 interflavan bonds. When the interflavan bond is positioned ahead of ring C, it is calledβ; but when this is oriented backward of ring C, it is appointed as α. The upper units of the chain are denominated extension units (EU) and the end unit is known as terminal unit (TU) [2]. Since there are several chiral centers (C2, C3, and C4) in these molecules, they can exist as a variety of stereoisomers. Due to the flexibility of ring C, the compact (Com) and extended (Ext) conformations are formed (Figure 1) [3]. In the conformationCom, the two rings B are opposite to each other, whereas they are facing each other in the conformationExt [4].Figure 1
General structure of procyanidin dimers with interflavan bond types C4–C8 and C4–C6. (a) PB3 (C4α–C8, catechin); PB4 (C4α–C8, catechin-epicatechin). (b) PB5 (C4β–C6, epicatechin); PB6 (C4α–C6, catechin); PB7 (C4β–C6, epicatechin-catechin); PB8 (C4α–C6, catechin-epicatechin).
(a) (b)The PCs are abundant in fruits, cereals, legumes, and their derived processed foods. Approximately 40% of the PCs ingested in the diet correspond to monomers, dimers, trimers, and tetramers [5, 6]. The dimers are not only abundant but also bioavailable, so that some common dimers in plan tissues as PB1 [(−)-epicatechin-(C4β→C8)-(+)-catechin] and PB2 [(−)-epicatechin-(C4β→C8)-(−)-epicatechin] have been preferentially studied. Other PC dimers, including PB3 [(+)-catechin-(C4α→C8)-(+)-catechin], PB4 [(+)-catechin-(C4α→C8)-(−)-epicatechin], PB5 [(−)-epicatechin-(C4β→C6)-(−)-epicatechin], PB6 [(+)-catechin-(C4α→C6)-(+)-catechin] PB7 [(−)-epicatechin-(C4β→C6)-(+)-catechin], and PB8 [(+)-catechin-(C4α→C6)-(−)-epicatechin], are less frequent in foods but are important components of fruits and cereals of high commercial value, among these grape, apple, and barley, as well as derived beverages, such as juice, red wine, and beer [7–9]. Together, these dimers afford contrasting structural characteristics related to the type of constituent units, type of interflavan bond, and conformation, whose study could provide insights concerning the role of molecular structure of PCs on their antioxidant properties. Previously, De Freitas et al. [7] performed an experimental-computational study that addressed the relationship between molecular structure and susceptibility to oxidation of those dimeric PCs. However, data relating the electron transfer and electron density changes were not taken into account.The antioxidant reactions operate through electron transfer from a donor molecule (PCs) to an acceptor one (enzyme, metal, and free radical). The electron transfer process can occur by different ways: either directly or in the form of hydrogen atoms; by combination with molecular oxygen or, even, by using electron transporters (NAD+, NADP+, FAD, FMN) [10]. The antioxidant reactions also can occur through nucleophilic or electrophilic attack.Since molecules reactivity depends on its structure, the first step to explore the antioxidant properties of compounds of interest is to know the most susceptible molecular sites for donating and accepting electrons, as well as the relative stabilities of the molecules.Computational chemistry is a useful tool for studying the reactivity properties of compounds, especially those with a great structural diversity like the PCs. Through a computational chemistry approach, diverse analysis can be executed which offer a better sight of the effect of structure on the antioxidant properties. In relation to the type of constituent units, type of interflavan bond, and conformation, parameters such as the dihedral angles of the interflavan bond, intramolecular hydrogen bonds between monomeric units, molecular surface area, and stability of the PCs may be of great relevance. Regarding the capacity for transferring electrons, the chemical potential is essential, as well as thermochemical parameters such as bond dissociation enthalpy (BDE), gap energy, Fukui indices, and charge distribution of HOMO-LUMO orbitals to determine sites of nucleophilic and electrophilic attack.The density functional theory (DFT) is among the computational methods which offer a high accuracy for simulating distinct conformations and exploring properties based on electron transfer, besides changes in electron density in complex molecular systems.The objective of this study was to analyze the effect of the constituent units, type of interflavan bond, and conformation on the antioxidant properties of PCs dimers including PB3, PB4, PB5, PB6, PB7, and PB8, employing the DFT method.


## 2. Computational Methods

The ground-state molecular geometries, harmonic vibrational frequencies, and reactivity calculations were performed under the simulation of an aqueous medium, within the framework of DFT method implemented in the GAUSSIAN 09W computational package [11].The M05-2X/6-31G∗∗ model chemistry in combination with the SMD solvation method was selected for this study because they are considered appropriate for the type of variables to be analyzed. The M05-2X is a highly parameterized hybrid metaexchange-correlation functional, which yields satisfactory results for thermochemistry parameters and noncovalent interactions (hydrogen bonding) [12, 13]. On the other hand, the SMD model is applicable to any charged or uncharged solute molecule interacting with a continuum solvent. The elements present in the solutes included H, C, and O, which are part of the PCs. The SMD model employs a single set of parameters, among them the intrinsic atomic Coulomb radii and atomic surface tension coefficients, optimized over six electronic structure methods, within which are included M05-2X/6-31G∗ and M05-2X/6-31+G∗∗. With the use of 6-31G∗ basis set, the SMD model achieves mean errors of 0.6–1.0 kcal/mol in the solvation free energies of neutral molecules and mean errors around of 4 kcal/mol for ions [14].For a better understanding of the effect of the type of constituent units, type of interflavan bond, and conformation on the antioxidant properties of PC dimers, parameters such as the dihedral angles of interflavan bond, intramolecular hydrogen bonds between monomeric units, the molecular surface area, and stability of the PCs were determined.The dihedral angles of PC dimers were measured among C3–C4–C8–C9 and C3–C4–C6–C7 carbon atoms for C4–C8 and C4–C6 interflavan bonds, respectively.To know the molecular surface area of the PC dimers, a density cube file was generated from the optimized geometries in Gaussian 09; then the Chimera software 1.11rc [15] was used to read the cube data file and the surface area of PC dimers was determined at ±0.001 a.u. isodensity value.To estimate the stability of the PC dimers, the Gibbs free energy including the thermal free energy corrections (E0+Gcorr) was obtained as one of the thermochemical values, computed in the Gaussian 09 program from a vibrational frequency calculation in aqueous medium, using the rigid-rotor harmonic-oscillator approximation (RRHO) at 298 K [16]. The stability of the PC dimers was established as a relative value of free energy (ΔE), that is, as the difference in free energy between the lowest-energy dimer and the remaining dimers in their distinct conformations.The calculated reactivity properties included those parameters that define the chemical potential, such as ionization potential, electron affinity, hardness, electronegativity, and electrophilicity; also, the HOMO-LUMO gap energies and thermochemical parameters such as bond dissociation enthalpy (BDE), Fukui indices, and charge distribution of HOMO-LUMO orbitals were included.The parameters of the chemical potential were obtained from the differences of total electronic energies when adding or removing an electron, in relation to the neutral molecule, using a vertical approach that is without reoptimization of the molecular geometries [17]. The equations used to calculate these parameters are as follows:Ionization potential is(1)I=cationenergy-neutralmoleculeenergy.Electron affinity is(2)A=neutralmolecularenergy-anionenergy.Hardness is(3)η=12I-A.Electronegativity is(4)χ=12I+A.Electrophilicity is(5)ω=μ22η,where μ=-χ, the negative of electronegativity (χ) is equivalent to chemical potential (μ).The HOMO-LUMO gap energies were calculated as the difference in energy between the frontier orbitals HOMO and LUMO, as indicated in the following equation [18]:(6)gap=ELUMO-EHOMO.In order to know the ability of the PC dimers to transfer electrons in form of hydrogen atoms, the bond dissociation enthalpy (BDE) was calculated. This was determined as the difference of enthalpy between parent molecule (ArOH) and its phenoxyl radical (ArO•) formed after the hydrogen atom (H•) transfer from each phenolic or hydroxyl group (OH). The reaction (7) and (8) were applied [19]:(7)ArOH⟶ArO•+H•(8)ΔrH298K=ΣE0+Hcorrproducts-ΣE0+Hcorrreactants.The computed enthalpy of hydrogen atom (H•) in aqueous phase was −0.493267 Hartrees.The Fukui function is a local property derived from the electronic density and is defined by the following relation:(9)fr=∂ρr∂Nvr,whereρ represents the electronic density of the system under consideration [20].The condensed Fukui functions are found by taking the finite difference approximations from charge population analysis of atoms in molecules, depending on the direction of the electron transfer; for nucleophilic and electrophilic attack, the functions are described as follows [21]:(10)fk+=qkN+1-qkNnucleophilicattackfk-=qkN-qkN-1electrophilicattack,where qk is the gross charge of atomk in the molecule.The charge distributions of HOMO and LUMO orbitals were visualized from the optimized geometries of PC dimers in aqueous medium by GaussView 5.0.9 program.


## 3. Results and Discussion

### 3.1. Molecular Properties and Stability of the PC Dimers

Previous studies have shown that PC dimers can acquire two energy-minimized structures, which are identified based on the dihedral angle. One of the optimized structures is characterized by adopting a conformationCom, while the others adopt a conformationExt [4, 22].The dihedral angles established among the C3–C4–C8–C9 or C3–C4–C6–C7 carbon atoms of the PC dimers optimized in aqueous medium are displayed in Table1. PB3 and PB4 have a C4α–C8 interflavan bond, so they can acquireCom andExt conformations; therefore the dihedral angles formed in both conformations are shown. Nonetheless, in PB5, PB6, PB7, and PB8, which possess C4α–C6 or C4β–C6 interflavan bonds, only one conformation can be adopted in its ground-state. In the case of PB5 and PB7, which have C4β–C6 interflavan bond (ahead of ring C), and (−)-epicatechin as extension unit, a conformationCom was acquired, while, in PB6 and PB8 with C4α–C6 interflavan bond (backward of ring C) and (+)-catechin as extension unit, a conformation Ext was adopted.Table 1
Molecular properties and stability (ΔE) of the procyanidin dimers optimized in aqueous medium (SMD model) using the M05-2X/6-31G∗∗ DFT method.
Procyanidin dimer Dihedral angle (degrees) Molecular surface area(Å2) Relative energy (ΔE kcal/mol) PB3Com 114.20 (C3–C4–C8–C9) 745.9 3.49 PB3Ext −35.39 (C3–C4–C8–C9) 678.5 6.49 PB4Com 102.08 (C3–C4–C8–C9) 744.4 0 PB4Ext −65.12 (C3–C4–C8–C9) 677.3 0.23 PB5 101.72 (C3–C4–C6–C7) 748.9 5.05 PB6 −59.32 (C3–C4–C6–C7) 751.4 3.12 PB7 101.37 (C3–C4–C6–C7) 745.3 4.63 PB8 −59.49 (C3–C4–C6–C7) 749.8 4.55 PB3 (C4α–C8, catechin); PB4 (C4α–C8, catechin-epicatechin); PB5 (C4β–C6, epicatechin); PB6 (C4α–C6, catechin); PB7 (C4β–C6, epicatechin-catechin); PB8 (C4α–C6, catechin-epicatechin). Δ E was calculated as the difference in free energy between the lowest-energy dimer and the remaining dimers.It is important to clarify that the dihedral angles obtained from the energy-minimized structures of the PC dimers do not come from calculations of torsional profiles from 0° to 360°. The way we check that the structures correspond to absolute minima or true minima was calculating the force constants and vibrational frequencies on the stationary points obtained after the optimization of PC dimers. With such data, the presence of imaginary frequencies (negative values) was discarded and it was verified that the forces representing the minima remained constant in all directions of the molecular plane. For the last point, four criteria that confirm convergence were taken into account [11] and these are as follows:(1)
The maximum component of the force must be below the cutoff value of 0.00045 in all directions of the molecular plane.(2)
The root-mean-square (RMS) of the forces must be below the defined tolerance of 0.0003.(3)
The calculated displacement for the next step must be smaller than the defined cutoff value of 0.0018.(4)
The RMS of the displacement for the next step must be below its cutoff value of 0.0012.The three-dimensional orientation of PC dimers induced by the adopted conformations (Com orExt) can lead to the formation of intramolecular hydrogen bonds (IHBs) between their constituent units. These bonds exert important effects on the interactions that the conformers may establish with other molecules such as oxidants or solvents. The analysis of interatomic bond distances of the PC dimers optimized in aqueous medium suggest that, in PB3Com, PB4Com, PB5, PB6, PB7, and PB8, IHBs were not formed. Instead, the conformationExt of PB3 and PB4 can favor the formation of IHBs between the catechol moieties (B rings) of the extension and terminal units, with values within a range from 2.0 to 2.53 Å (Figure 2).Figure 2
Optimized geometry of the procyanidin dimers in aqueous medium (SMD model) using the M05-2X/6-31G∗∗ DFT method.Com: compact conformation;Ext: extended conformation.With the purpose of corroborating the establishment of these IHBs, the O–H stretching frequencies of the catechol moieties were analyzed, and these appear around 3700 cm−1 in PB3Ext and PB4Ext. The literature indicates that, in orthosubstituted phenols, which can form IHBs, it is possible to locate the O–H stretch band in the region of 3670–3607 cm−1 in vapor phase and this band could be shifted to higher frequencies in dilute solutions. The abnormally high O–H stretching frequencies are attributed to the repulsion between O–H groups and protons of other groups [23]. In the case of PB3Com, PB4Com, PB5, PB6, PB7, and PB8, the O–H stretching frequencies including the catechol moieties were detected in the range of 3900–3800 cm−1 region, where the free O–H stretching in phenols can be observed.The above results are different to data reported in other studies, in which IHBs were favored in conformationCom but not in the conformationExt of distinct PC oligomers linked by C4β–C8 interflavan bonds, among which are the PB1 and PB2 dimers, PC1 trimer, and a B-type tetramer [24, 25]. These results highlight the importance of the extension unit monomer and the atoms involved in the interflavan bond (C4–C8 or C4–C6), as well as the bond orientation (α or β). Those that can favor the formation of IHBs are the joins C4–C8 but not the C4–C6. When the extension unit is (−)-epicatechin and the interflavan bond acquires the β orientation (C4β–C8), IHBs in conformationCom can be promoted, but if the extension unit is (+)-catechin and the interflavan bond acquires the α orientation (C4α–C8), IHBs in conformationExt could be favored.The conformationExt is characterized by a rearrangement where the two rings B are facing each other. Depending on the proximity of these catechol rings, an intramolecular π-π stacking arrangement can be adopted leading to the formation of IHBs, which reduces the total surface of the molecule [4]. Considering the repercussion that this could have on the accessibility of the solvent or oxidizing agents in the reactive sites of the molecules, the molecular surface area of PC dimers was measured.In Table1, it can be seen that PB5, PB6, PB7, and PB8 with C4–C6 interflavan bond showed the greatest molecular surface areas (745 to 751 Å2), whereas PB3 and PB4 with C4–C8 interflavan bond in their conformationExt showed the smallest molecular surface areas (677 to 678 Å2). Figure 3 depicts the total electronic density of the molecules of PB6, PB3Ext, and PB4Ext, in order to visualize the difference between the molecule with greater surface area (PB6) and the molecules with smaller surface area (PB3Ext and PB4Ext).Figure 3
Molecular surface area of the PB3Ext, PB4Ext, and PB6 procyanidin dimers, determined under simulation of an aqueous medium (SMD model) using the M05-2X/6-31G∗∗ DFT method.Ext: extended conformation.Regarding the stability of PC dimers, PB4 was the one that exhibited the most negative values of free energy in aqueous medium, mainly in the conformationCom followed by the conformationExt, indicating that it was the most stable dimer. The relative free energy (ΔE) between conformationsCom andExt of PB4 was 0.23 kcal/mol, suggesting that they could be exchanged easily (Table 1). However, in a study involving nuclear magnetic resonance and molecular modeling, it was reported that, in PB4, the conformationCom predominates in aqueous ethanol solution [7]. In another similar study, it was found that, in PB2 dimer, which has (−)-epicatechin as terminal unit such as PB4, even proportions ofCom andExt conformations were detected when it was dissolved in water and hydroalcoholic media [26].The relative stability of the dimers based on the lowest values ofΔE and taking PB4 as reference because of being the most stable was as follows (Table 1).PB4Com (0 kcal/mol) > PB4Ext (0.23 kcal/mol) > PB6 (3.12 kcal/mol) > PB3Com (3.49 kcal/mol) > PB8 (4.55 kcal/mol) > PB7 (4.63 kcal/mol) > PB5 (5.05) > PB3Ext (6.49 kcal/mol).As can be appreciated in Table1, PB3 in its conformationExt is the less stable dimer even establishing two IHBs between the B rings of both extension and terminal units. On the one hand, this can be attributed to the fact that the dihedral angle of the interflavan bond is small (−35.39) and generates steric hindrance among the remaining O–H groups by the proximity of the aromatic rings. On the other hand, it is important to take into account the fact that all molecules have an effective surface area, which is determined by the conformation, and this is more important for the interactions with solvent than the total surface area. In addition to the formation of IHBs, the establishment of intermolecular hydrogen bonds with the solvent can contribute to the stability of molecules in solution and marking the differences of solvation free energy among PC dimers. To this respect, it is important to mention that, in the SMD solvation model, the interactions of the solute and the surrounding solvent do not take into account the molecular structure of the solvent and that the dielectric response of the medium is uniform and linear at all positions outside the space that defines the solute. For this reason, the SMD model does not provide a completely valid description neither of the solvent structure nor of the solute-solvent interactions (hydrogen bonding, π-π stacking interactions). However, the SMD model allows estimating the geometry-dependent proportionality constants called atomic surface tensions that describe the solvent-accessible surface areas of the individual atoms of the solute [14]. On these bases, it is considered that the SMD solvation model could give an approximate and acceptable description of the relative stabilities of the PC dimers in their distinct conformations.


### 3.2. Chemical Potential, Gap Energy, and Stability of Anions and Cations

The chemical potential is defined by a set of parameters which measure the tendency to transfer or capture electrons and therefore describe the antioxidant properties of the molecules. The parameters that determine the ability to transfer electrons are the ionization potential (I) and hardness (η). I represents the amount of energy required to remove an electron of a molecule and η is a measure of the resistance to charge transference [18]. On the other hand, the ability to capture electrons is determined by the electron affinity (A), electronegativity (χ), and electrophilicity (ω). A is defined as the energy release when an electron is added to a neutral molecule, while χ indicates the tendency to attract electrons and ω establishes the decay of binding energy due to a maximum electron flow between a donor and a receptor [18].PCs, like all flavonoids, have atomic regions that allow them to easily transfer electrons but also have other atomic regions through which they react accepting electrons from other compounds. Molecules with improved antioxidant properties are those that readily transfer the electrons and have a low capacity to capture them, that is, having low values ofI, A, η, χ, and ω. It is known that a large value of A and ω may lead to prooxidant effects and a decrease in antioxidant ability [27, 28].As displayed in Table2, the chemical potential parameters fluctuated within a comparable range of values in all PC dimers: I (5.82–6.34 eV), A (0.28–0.40 eV), η (2.72–3.01 eV), χ (3.06–3.33 eV), ω (1.68–1.84 eV). These values are lower than those reported for quercetin [29], I (7.22 eV), A (0.76 eV), η (3.23 eV), χ (3.99 eV), ω (2.47 eV), which could mean that the PC dimers possess high antioxidant ability, since quercetin is considered one of the flavonoids with greater antioxidant activity.Table 2
Parameters of chemical potential and gap energy of the procyanidin dimers calculated by the M05-2X/6-31G∗∗ DFT method in aqueous medium (SMD model).
Procyanidin dimer I (eV) A (eV) η (eV) χ (eV) ω (eV) Gap (eV) PB3Com 5.84 0.28 2.78 3.06 1.68 6.10 PB3Ext 6.01 0.39 2.81 3.20 1.83 6.31 PB4Com 5.82 0.31 2.76 3.06 1.70 6.11 PB4Ext 5.84 0.40 2.72 3.12 1.79 6.27 PB5 5.95 0.34 2.81 3.15 1.77 6.37 PB6 5.99 0.29 2.85 3.15 1.74 6.27 PB7 5.92 0.28 2.82 3.09 1.70 6.34 PB8 6.34 0.32 3.01 3.33 1.84 6.32 PB3 (C4α–C8, catechin); PB4 (C4α–C8, catechin-epicatechin); PB5 (C4β–C6, epicatechin); PB6 (C4α–C6, catechin); PB7 (C4β–C6, epicatechin-catechin); PB8 (C4α–C6, catechin-epicatechin). I: ionization potential; A: electron affinity; η: hardness; χ: electronegativity; ω: electrophilicity.Notwithstanding the previously discussed, differences were observed among the analyzed PC dimers. In general terms, the dimers PB3 and PB4 in their conformationCom were the ones that showed the lowest values for all the variables associated with the chemical potential, while the dimer PB8 presented the highest values (Table 2).The band gap represents the conductive layer of electrons in a molecule and is determined by the energy difference between the valence orbitals HOMO and LUMO. Reactive molecular systems that easily transfer their electrons are associated with low gap values, whereas, in the opposite case, the less reactive molecules have higher gap values [17, 18]. The gap values obtained in our study, shown in Table 2, suggest that PB3 and PB4 in their conformationCom could be the dimers that more easily transfer their electrons, since the energy difference between HOMO and LUMO was lower (6.10–6.11 eV) in comparison to the PB5, PB6, PB7, and PB8 dimers (6.27–6.37 eV). These results corroborate those found for the chemical potential.


### 3.3. Bond Dissociation Enthalpy (BDE)

One of the most common ways of studying the ability of flavonoids to inhibit free radicals is through thermochemical parameters, among these the BDE, because the electron transfer in the form of hydrogen atoms is one of the preferential mechanisms followed by flavonoids [24, 30]. In Table 3, it can be seen that all the analyzed PC dimers showed the same behavior regarding the easiness of each phenolic group for donating its hydrogen atom; the lower the BDE value, the higher the easiness for transferring the hydrogen atom. On these bases, the following general sequence for the distinct OH groups of the dimers is proposed, in both the extension and terminal units.Table 3
Bond dissociation enthalpy (BDE) of the procyanidin dimers calculated in aqueous medium (SMD model) using the M05-2X/6-31G∗∗ DFT method.
Procyanidin dimer BDE (kcal/mol) Extension unit (EU) Terminal unit (TU) 3′-OH 4′-OH 5-OH 7-OH 3′-OH 4′-OH 5-OH 7-OH PB3Com 85.51 82.36 88.57 88.77 84.04 83.03 87.27 87.21 PB3Ext 84.60 81.87 88.53 89.92 86.11 82.42 86.51 85.78 PB4Com 86.27 83.03 88.28 88.79 85.70 82.42 86.96 87.10 PB4Ext 83.85 81.91 88.27 89.11 85.32 82.28 86.89 86.98 PB5 86.45 83.57 90.64 92.79 86.53 83.30 84.98 89.07 PB6 86.16 82.87 88.46 90.86 85.83 82.50 85.53 87.42 PB7 86.83 82.96 90.60 92.79 85.82 82.50 84.72 88.70 PB8 85.78 82.50 88.71 90.98 86.37 83.33 85.26 86.84 PB3 (C4α–C8, catechin); PB4 (C4α–C8, catechin-epicatechin); PB5 (C4β–C6, epicatechin); PB6 (C4α–C6, catechin); PB7 (C4β–C6, epicatechin-catechin); PB8 (C4α–C6, catechin-epicatechin).Easiness to transfer hydrogen atoms (BDE values) from the OH groups in the extension unit of PC dimers: 4′-OH (81.87–83.03 kcal/mol) > 3′-OH (83.85–86.83 kcal/mol) > 5-OH (88.27–90.64 kcal/mol) > 7-OH (88.77–92.79 kcal/mol).Easiness to transfer hydrogen atoms (BDE values) from the OH groups in the terminal unit of PC dimers: 4′-OH (82.28–83.33 kcal/mol) > 3′-OH (84.04–86.53 kcal/mol) > 5-OH (84.72–87.27 kcal/mol) > 7-OH (85.78–89.07 kcal/mol).The BDE values displayed above suggest that the 4′-OH and 3′-OH groups, which constitute the catechol moiety of ring B, show more facility to transfer their hydrogen atoms in both the extension and terminal units.The phenolic groups that show the greatest facility for transferring their hydrogen atoms (less BDE) are those which confer better antioxidant properties to the molecules. On these bases, it is suggested that the 4′-OH group is the most antioxidant site of the PC dimers.The BDE values of the phenoxyl radicals derived from the 4′-OH group of the different dimers were very close to each other; however, these tend to be slightly lower for PB3 and PB4, especially in the conformationExt, indicating greater antioxidant ability for such dimers compared to PB5, PB6, PB7, and PB8. These results are consistent with those reported by [7], who performed a kinetic study of oxidative susceptibility in the same dimers. The compounds with higher susceptibility to oxidation are those that may present better antioxidant properties by protecting other molecules from oxidation.


### 3.4. Fukui Indices and Charge Distribution in HOMO and LUMO Orbitals

The Fukui indices calculation was carried out using the Hirshfeld charge population analysis. Through this scheme it is possible to study well-defined atomic fragments, which when they integrated, they allow establishing the contribution of each fragment in the changes of electronic density, in the deformation of the bonds, electrostatic potential, and the transferability of net atomic charge at molecular level [31].During the reaction of the PCs against different radical species, oxidative enzymes, and metals, these can develop effective charges in the transition state, establishing sites susceptible to nucleophilic and electrophilic attack. According to the Fukui indices displayed in Table4, the preferential site for nucleophilic attack in all PC dimers was the C4′ atom (catechol group, ring B). For PB3Com, PB4Com, PB5, PB6, PB7, and PB8, this site was located in the extension unit (EU), while, for PB3Ext and PB4Ext, the site was located in the terminal unit (TU). These results are supported by the LUMO distribution (isovalue 0.02), which represents the molecular site for the positive charge density (Figure 4). The LUMO distribution also points out to other atoms that are potential sites for nucleophilic attack, but the charge density is mainly concentrated on ring B, either of the extension unit or of the terminal unit, which is where the C4′ atom is located.Table 4
Fukui indices of the procyanidin dimers calculated by the M05-2X/6-31G∗∗ DFT method using the Hirshfeld charge population analysis in an aqueous medium (SMD model).
Procyanidin dimer Atom/unit Nucleophilic attack Atom/unit Electrophilic attack PB3Com C4′ EU 0.095 O5 TU 0.085 PB3Ext C4′ TU 0.097 O4′ EU 0.104 PB4Com C4′ EU 0.095 O5 TU 0.088 PB4Ext C4′ TU 0.098 O5 TU 0.083 PB5 C4′ EU 0.105 O3′ EU 0.116 PB6 C4′ EU 0.096 O7 TU 0.056 PB7 C4′ EU 0.095 O7 TU 0.099 PB8 C4′ EU 0.097 O4′ EU 0.049 PB3 (C4α–C8, catechin); PB4 (C4α–C8, catechin-epicatechin); PB5 (C4β–C6, epicatechin); PB6 (C4α–C6, catechin); PB7 (C4β–C6, epicatechin-catechin); PB8 (C4α–C6, catechin-epicatechin). EU: extension unit; TU: terminal unit.Figure 4
Charge distribution in the LUMO orbital of the procyanidin dimers optimized in aqueous medium (SMD model) using the M05-2X/6-31G∗∗ DFT method.Com: compact conformation;Ext: extended conformation.The site for electrophilic attack was the most variable for the PC dimers. Table4 shows that in PB3Com, PB4Com, PB4Ext, PB6, and PB7 dimers, the electrophilic attack occurs preferentially either in the O5 or in O7, atoms that are part of the resorcinol moiety of ring A. The HOMO distribution in such molecules corroborates the electron density that is mainly concentered on ring A, especially in the terminal unit (Figure 5). In contrast, PB3Ext, PB5, and PB8 showed as a preferential site for electrophilic attack, either O3′ or O4′, which are part of the catechol moiety of ring B in the extension unit.Figure 5
Charge distribution in the HOMO orbital of the procyanidin dimers optimized in aqueous medium (SMD model) using the M05-2X/6-31G∗∗ DFT method.Com: compact conformation;Ext: extended conformation.The PB5 dimer proved to be the most susceptible for both nucleophilic and electrophilic attacks, since it presented the highest Fukui indices (0.105 and 0.116 eV, resp.). The PB3, PB4, PB6, PB7, and PB8 dimers showed susceptibilities for nucleophilic attack, similar to each other (Fukui indices 0.095–0.098 eV). However, the susceptibilities for electrophilic attack showed notable differences among these dimers (Fukui indices 0.049–0.116 eV), corresponding to the lowest value to the PB8 dimer.


## 3.1. Molecular Properties and Stability of the PC Dimers

Previous studies have shown that PC dimers can acquire two energy-minimized structures, which are identified based on the dihedral angle. One of the optimized structures is characterized by adopting a conformationCom, while the others adopt a conformationExt [4, 22].The dihedral angles established among the C3–C4–C8–C9 or C3–C4–C6–C7 carbon atoms of the PC dimers optimized in aqueous medium are displayed in Table1. PB3 and PB4 have a C4α–C8 interflavan bond, so they can acquireCom andExt conformations; therefore the dihedral angles formed in both conformations are shown. Nonetheless, in PB5, PB6, PB7, and PB8, which possess C4α–C6 or C4β–C6 interflavan bonds, only one conformation can be adopted in its ground-state. In the case of PB5 and PB7, which have C4β–C6 interflavan bond (ahead of ring C), and (−)-epicatechin as extension unit, a conformationCom was acquired, while, in PB6 and PB8 with C4α–C6 interflavan bond (backward of ring C) and (+)-catechin as extension unit, a conformation Ext was adopted.Table 1
Molecular properties and stability (ΔE) of the procyanidin dimers optimized in aqueous medium (SMD model) using the M05-2X/6-31G∗∗ DFT method.
Procyanidin dimer Dihedral angle (degrees) Molecular surface area(Å2) Relative energy (ΔE kcal/mol) PB3Com 114.20 (C3–C4–C8–C9) 745.9 3.49 PB3Ext −35.39 (C3–C4–C8–C9) 678.5 6.49 PB4Com 102.08 (C3–C4–C8–C9) 744.4 0 PB4Ext −65.12 (C3–C4–C8–C9) 677.3 0.23 PB5 101.72 (C3–C4–C6–C7) 748.9 5.05 PB6 −59.32 (C3–C4–C6–C7) 751.4 3.12 PB7 101.37 (C3–C4–C6–C7) 745.3 4.63 PB8 −59.49 (C3–C4–C6–C7) 749.8 4.55 PB3 (C4α–C8, catechin); PB4 (C4α–C8, catechin-epicatechin); PB5 (C4β–C6, epicatechin); PB6 (C4α–C6, catechin); PB7 (C4β–C6, epicatechin-catechin); PB8 (C4α–C6, catechin-epicatechin). Δ E was calculated as the difference in free energy between the lowest-energy dimer and the remaining dimers.It is important to clarify that the dihedral angles obtained from the energy-minimized structures of the PC dimers do not come from calculations of torsional profiles from 0° to 360°. The way we check that the structures correspond to absolute minima or true minima was calculating the force constants and vibrational frequencies on the stationary points obtained after the optimization of PC dimers. With such data, the presence of imaginary frequencies (negative values) was discarded and it was verified that the forces representing the minima remained constant in all directions of the molecular plane. For the last point, four criteria that confirm convergence were taken into account [11] and these are as follows:(1)
The maximum component of the force must be below the cutoff value of 0.00045 in all directions of the molecular plane.(2)
The root-mean-square (RMS) of the forces must be below the defined tolerance of 0.0003.(3)
The calculated displacement for the next step must be smaller than the defined cutoff value of 0.0018.(4)
The RMS of the displacement for the next step must be below its cutoff value of 0.0012.The three-dimensional orientation of PC dimers induced by the adopted conformations (Com orExt) can lead to the formation of intramolecular hydrogen bonds (IHBs) between their constituent units. These bonds exert important effects on the interactions that the conformers may establish with other molecules such as oxidants or solvents. The analysis of interatomic bond distances of the PC dimers optimized in aqueous medium suggest that, in PB3Com, PB4Com, PB5, PB6, PB7, and PB8, IHBs were not formed. Instead, the conformationExt of PB3 and PB4 can favor the formation of IHBs between the catechol moieties (B rings) of the extension and terminal units, with values within a range from 2.0 to 2.53 Å (Figure 2).Figure 2
Optimized geometry of the procyanidin dimers in aqueous medium (SMD model) using the M05-2X/6-31G∗∗ DFT method.Com: compact conformation;Ext: extended conformation.With the purpose of corroborating the establishment of these IHBs, the O–H stretching frequencies of the catechol moieties were analyzed, and these appear around 3700 cm−1 in PB3Ext and PB4Ext. The literature indicates that, in orthosubstituted phenols, which can form IHBs, it is possible to locate the O–H stretch band in the region of 3670–3607 cm−1 in vapor phase and this band could be shifted to higher frequencies in dilute solutions. The abnormally high O–H stretching frequencies are attributed to the repulsion between O–H groups and protons of other groups [23]. In the case of PB3Com, PB4Com, PB5, PB6, PB7, and PB8, the O–H stretching frequencies including the catechol moieties were detected in the range of 3900–3800 cm−1 region, where the free O–H stretching in phenols can be observed.The above results are different to data reported in other studies, in which IHBs were favored in conformationCom but not in the conformationExt of distinct PC oligomers linked by C4β–C8 interflavan bonds, among which are the PB1 and PB2 dimers, PC1 trimer, and a B-type tetramer [24, 25]. These results highlight the importance of the extension unit monomer and the atoms involved in the interflavan bond (C4–C8 or C4–C6), as well as the bond orientation (α or β). Those that can favor the formation of IHBs are the joins C4–C8 but not the C4–C6. When the extension unit is (−)-epicatechin and the interflavan bond acquires the β orientation (C4β–C8), IHBs in conformationCom can be promoted, but if the extension unit is (+)-catechin and the interflavan bond acquires the α orientation (C4α–C8), IHBs in conformationExt could be favored.The conformationExt is characterized by a rearrangement where the two rings B are facing each other. Depending on the proximity of these catechol rings, an intramolecular π-π stacking arrangement can be adopted leading to the formation of IHBs, which reduces the total surface of the molecule [4]. Considering the repercussion that this could have on the accessibility of the solvent or oxidizing agents in the reactive sites of the molecules, the molecular surface area of PC dimers was measured.In Table1, it can be seen that PB5, PB6, PB7, and PB8 with C4–C6 interflavan bond showed the greatest molecular surface areas (745 to 751 Å2), whereas PB3 and PB4 with C4–C8 interflavan bond in their conformationExt showed the smallest molecular surface areas (677 to 678 Å2). Figure 3 depicts the total electronic density of the molecules of PB6, PB3Ext, and PB4Ext, in order to visualize the difference between the molecule with greater surface area (PB6) and the molecules with smaller surface area (PB3Ext and PB4Ext).Figure 3
Molecular surface area of the PB3Ext, PB4Ext, and PB6 procyanidin dimers, determined under simulation of an aqueous medium (SMD model) using the M05-2X/6-31G∗∗ DFT method.Ext: extended conformation.Regarding the stability of PC dimers, PB4 was the one that exhibited the most negative values of free energy in aqueous medium, mainly in the conformationCom followed by the conformationExt, indicating that it was the most stable dimer. The relative free energy (ΔE) between conformationsCom andExt of PB4 was 0.23 kcal/mol, suggesting that they could be exchanged easily (Table 1). However, in a study involving nuclear magnetic resonance and molecular modeling, it was reported that, in PB4, the conformationCom predominates in aqueous ethanol solution [7]. In another similar study, it was found that, in PB2 dimer, which has (−)-epicatechin as terminal unit such as PB4, even proportions ofCom andExt conformations were detected when it was dissolved in water and hydroalcoholic media [26].The relative stability of the dimers based on the lowest values ofΔE and taking PB4 as reference because of being the most stable was as follows (Table 1).PB4Com (0 kcal/mol) > PB4Ext (0.23 kcal/mol) > PB6 (3.12 kcal/mol) > PB3Com (3.49 kcal/mol) > PB8 (4.55 kcal/mol) > PB7 (4.63 kcal/mol) > PB5 (5.05) > PB3Ext (6.49 kcal/mol).As can be appreciated in Table1, PB3 in its conformationExt is the less stable dimer even establishing two IHBs between the B rings of both extension and terminal units. On the one hand, this can be attributed to the fact that the dihedral angle of the interflavan bond is small (−35.39) and generates steric hindrance among the remaining O–H groups by the proximity of the aromatic rings. On the other hand, it is important to take into account the fact that all molecules have an effective surface area, which is determined by the conformation, and this is more important for the interactions with solvent than the total surface area. In addition to the formation of IHBs, the establishment of intermolecular hydrogen bonds with the solvent can contribute to the stability of molecules in solution and marking the differences of solvation free energy among PC dimers. To this respect, it is important to mention that, in the SMD solvation model, the interactions of the solute and the surrounding solvent do not take into account the molecular structure of the solvent and that the dielectric response of the medium is uniform and linear at all positions outside the space that defines the solute. For this reason, the SMD model does not provide a completely valid description neither of the solvent structure nor of the solute-solvent interactions (hydrogen bonding, π-π stacking interactions). However, the SMD model allows estimating the geometry-dependent proportionality constants called atomic surface tensions that describe the solvent-accessible surface areas of the individual atoms of the solute [14]. On these bases, it is considered that the SMD solvation model could give an approximate and acceptable description of the relative stabilities of the PC dimers in their distinct conformations.


## 3.2. Chemical Potential, Gap Energy, and Stability of Anions and Cations

The chemical potential is defined by a set of parameters which measure the tendency to transfer or capture electrons and therefore describe the antioxidant properties of the molecules. The parameters that determine the ability to transfer electrons are the ionization potential (I) and hardness (η). I represents the amount of energy required to remove an electron of a molecule and η is a measure of the resistance to charge transference [18]. On the other hand, the ability to capture electrons is determined by the electron affinity (A), electronegativity (χ), and electrophilicity (ω). A is defined as the energy release when an electron is added to a neutral molecule, while χ indicates the tendency to attract electrons and ω establishes the decay of binding energy due to a maximum electron flow between a donor and a receptor [18].PCs, like all flavonoids, have atomic regions that allow them to easily transfer electrons but also have other atomic regions through which they react accepting electrons from other compounds. Molecules with improved antioxidant properties are those that readily transfer the electrons and have a low capacity to capture them, that is, having low values ofI, A, η, χ, and ω. It is known that a large value of A and ω may lead to prooxidant effects and a decrease in antioxidant ability [27, 28].As displayed in Table2, the chemical potential parameters fluctuated within a comparable range of values in all PC dimers: I (5.82–6.34 eV), A (0.28–0.40 eV), η (2.72–3.01 eV), χ (3.06–3.33 eV), ω (1.68–1.84 eV). These values are lower than those reported for quercetin [29], I (7.22 eV), A (0.76 eV), η (3.23 eV), χ (3.99 eV), ω (2.47 eV), which could mean that the PC dimers possess high antioxidant ability, since quercetin is considered one of the flavonoids with greater antioxidant activity.Table 2
Parameters of chemical potential and gap energy of the procyanidin dimers calculated by the M05-2X/6-31G∗∗ DFT method in aqueous medium (SMD model).
Procyanidin dimer I (eV) A (eV) η (eV) χ (eV) ω (eV) Gap (eV) PB3Com 5.84 0.28 2.78 3.06 1.68 6.10 PB3Ext 6.01 0.39 2.81 3.20 1.83 6.31 PB4Com 5.82 0.31 2.76 3.06 1.70 6.11 PB4Ext 5.84 0.40 2.72 3.12 1.79 6.27 PB5 5.95 0.34 2.81 3.15 1.77 6.37 PB6 5.99 0.29 2.85 3.15 1.74 6.27 PB7 5.92 0.28 2.82 3.09 1.70 6.34 PB8 6.34 0.32 3.01 3.33 1.84 6.32 PB3 (C4α–C8, catechin); PB4 (C4α–C8, catechin-epicatechin); PB5 (C4β–C6, epicatechin); PB6 (C4α–C6, catechin); PB7 (C4β–C6, epicatechin-catechin); PB8 (C4α–C6, catechin-epicatechin). I: ionization potential; A: electron affinity; η: hardness; χ: electronegativity; ω: electrophilicity.Notwithstanding the previously discussed, differences were observed among the analyzed PC dimers. In general terms, the dimers PB3 and PB4 in their conformationCom were the ones that showed the lowest values for all the variables associated with the chemical potential, while the dimer PB8 presented the highest values (Table 2).The band gap represents the conductive layer of electrons in a molecule and is determined by the energy difference between the valence orbitals HOMO and LUMO. Reactive molecular systems that easily transfer their electrons are associated with low gap values, whereas, in the opposite case, the less reactive molecules have higher gap values [17, 18]. The gap values obtained in our study, shown in Table 2, suggest that PB3 and PB4 in their conformationCom could be the dimers that more easily transfer their electrons, since the energy difference between HOMO and LUMO was lower (6.10–6.11 eV) in comparison to the PB5, PB6, PB7, and PB8 dimers (6.27–6.37 eV). These results corroborate those found for the chemical potential.


## 3.3. Bond Dissociation Enthalpy (BDE)

One of the most common ways of studying the ability of flavonoids to inhibit free radicals is through thermochemical parameters, among these the BDE, because the electron transfer in the form of hydrogen atoms is one of the preferential mechanisms followed by flavonoids [24, 30]. In Table 3, it can be seen that all the analyzed PC dimers showed the same behavior regarding the easiness of each phenolic group for donating its hydrogen atom; the lower the BDE value, the higher the easiness for transferring the hydrogen atom. On these bases, the following general sequence for the distinct OH groups of the dimers is proposed, in both the extension and terminal units.Table 3
Bond dissociation enthalpy (BDE) of the procyanidin dimers calculated in aqueous medium (SMD model) using the M05-2X/6-31G∗∗ DFT method.
Procyanidin dimer BDE (kcal/mol) Extension unit (EU) Terminal unit (TU) 3′-OH 4′-OH 5-OH 7-OH 3′-OH 4′-OH 5-OH 7-OH PB3Com 85.51 82.36 88.57 88.77 84.04 83.03 87.27 87.21 PB3Ext 84.60 81.87 88.53 89.92 86.11 82.42 86.51 85.78 PB4Com 86.27 83.03 88.28 88.79 85.70 82.42 86.96 87.10 PB4Ext 83.85 81.91 88.27 89.11 85.32 82.28 86.89 86.98 PB5 86.45 83.57 90.64 92.79 86.53 83.30 84.98 89.07 PB6 86.16 82.87 88.46 90.86 85.83 82.50 85.53 87.42 PB7 86.83 82.96 90.60 92.79 85.82 82.50 84.72 88.70 PB8 85.78 82.50 88.71 90.98 86.37 83.33 85.26 86.84 PB3 (C4α–C8, catechin); PB4 (C4α–C8, catechin-epicatechin); PB5 (C4β–C6, epicatechin); PB6 (C4α–C6, catechin); PB7 (C4β–C6, epicatechin-catechin); PB8 (C4α–C6, catechin-epicatechin).Easiness to transfer hydrogen atoms (BDE values) from the OH groups in the extension unit of PC dimers: 4′-OH (81.87–83.03 kcal/mol) > 3′-OH (83.85–86.83 kcal/mol) > 5-OH (88.27–90.64 kcal/mol) > 7-OH (88.77–92.79 kcal/mol).Easiness to transfer hydrogen atoms (BDE values) from the OH groups in the terminal unit of PC dimers: 4′-OH (82.28–83.33 kcal/mol) > 3′-OH (84.04–86.53 kcal/mol) > 5-OH (84.72–87.27 kcal/mol) > 7-OH (85.78–89.07 kcal/mol).The BDE values displayed above suggest that the 4′-OH and 3′-OH groups, which constitute the catechol moiety of ring B, show more facility to transfer their hydrogen atoms in both the extension and terminal units.The phenolic groups that show the greatest facility for transferring their hydrogen atoms (less BDE) are those which confer better antioxidant properties to the molecules. On these bases, it is suggested that the 4′-OH group is the most antioxidant site of the PC dimers.The BDE values of the phenoxyl radicals derived from the 4′-OH group of the different dimers were very close to each other; however, these tend to be slightly lower for PB3 and PB4, especially in the conformationExt, indicating greater antioxidant ability for such dimers compared to PB5, PB6, PB7, and PB8. These results are consistent with those reported by [7], who performed a kinetic study of oxidative susceptibility in the same dimers. The compounds with higher susceptibility to oxidation are those that may present better antioxidant properties by protecting other molecules from oxidation.


## 3.4. Fukui Indices and Charge Distribution in HOMO and LUMO Orbitals

The Fukui indices calculation was carried out using the Hirshfeld charge population analysis. Through this scheme it is possible to study well-defined atomic fragments, which when they integrated, they allow establishing the contribution of each fragment in the changes of electronic density, in the deformation of the bonds, electrostatic potential, and the transferability of net atomic charge at molecular level [31].During the reaction of the PCs against different radical species, oxidative enzymes, and metals, these can develop effective charges in the transition state, establishing sites susceptible to nucleophilic and electrophilic attack. According to the Fukui indices displayed in Table4, the preferential site for nucleophilic attack in all PC dimers was the C4′ atom (catechol group, ring B). For PB3Com, PB4Com, PB5, PB6, PB7, and PB8, this site was located in the extension unit (EU), while, for PB3Ext and PB4Ext, the site was located in the terminal unit (TU). These results are supported by the LUMO distribution (isovalue 0.02), which represents the molecular site for the positive charge density (Figure 4). The LUMO distribution also points out to other atoms that are potential sites for nucleophilic attack, but the charge density is mainly concentrated on ring B, either of the extension unit or of the terminal unit, which is where the C4′ atom is located.Table 4
Fukui indices of the procyanidin dimers calculated by the M05-2X/6-31G∗∗ DFT method using the Hirshfeld charge population analysis in an aqueous medium (SMD model).
Procyanidin dimer Atom/unit Nucleophilic attack Atom/unit Electrophilic attack PB3Com C4′ EU 0.095 O5 TU 0.085 PB3Ext C4′ TU 0.097 O4′ EU 0.104 PB4Com C4′ EU 0.095 O5 TU 0.088 PB4Ext C4′ TU 0.098 O5 TU 0.083 PB5 C4′ EU 0.105 O3′ EU 0.116 PB6 C4′ EU 0.096 O7 TU 0.056 PB7 C4′ EU 0.095 O7 TU 0.099 PB8 C4′ EU 0.097 O4′ EU 0.049 PB3 (C4α–C8, catechin); PB4 (C4α–C8, catechin-epicatechin); PB5 (C4β–C6, epicatechin); PB6 (C4α–C6, catechin); PB7 (C4β–C6, epicatechin-catechin); PB8 (C4α–C6, catechin-epicatechin). EU: extension unit; TU: terminal unit.Figure 4
Charge distribution in the LUMO orbital of the procyanidin dimers optimized in aqueous medium (SMD model) using the M05-2X/6-31G∗∗ DFT method.Com: compact conformation;Ext: extended conformation.The site for electrophilic attack was the most variable for the PC dimers. Table4 shows that in PB3Com, PB4Com, PB4Ext, PB6, and PB7 dimers, the electrophilic attack occurs preferentially either in the O5 or in O7, atoms that are part of the resorcinol moiety of ring A. The HOMO distribution in such molecules corroborates the electron density that is mainly concentered on ring A, especially in the terminal unit (Figure 5). In contrast, PB3Ext, PB5, and PB8 showed as a preferential site for electrophilic attack, either O3′ or O4′, which are part of the catechol moiety of ring B in the extension unit.Figure 5
Charge distribution in the HOMO orbital of the procyanidin dimers optimized in aqueous medium (SMD model) using the M05-2X/6-31G∗∗ DFT method.Com: compact conformation;Ext: extended conformation.The PB5 dimer proved to be the most susceptible for both nucleophilic and electrophilic attacks, since it presented the highest Fukui indices (0.105 and 0.116 eV, resp.). The PB3, PB4, PB6, PB7, and PB8 dimers showed susceptibilities for nucleophilic attack, similar to each other (Fukui indices 0.095–0.098 eV). However, the susceptibilities for electrophilic attack showed notable differences among these dimers (Fukui indices 0.049–0.116 eV), corresponding to the lowest value to the PB8 dimer.


## 4. Conclusions

The type of constituent units did not show a well-defined trend in relation to the effect that is exerted on the ability of PC dimers to transfer or capture electrons. Based on the analysis of the chemical potential and BDE, it is suggested that the factors with the most notable effects on the antioxidant properties of the PC dimers were the type of interflavan bond (C4–C8 or C4–C6) and the conformation (Com orExt). Under this context, the antioxidant ability of the dimers PB3 and PB4 containing the interflavan bond C4–C8, in their conformationCom, was very similar to each other but greater than those of the dimers PB5, PB6, PB7, and PB8 containing the C4–C6 interflavan bond. It was also observed that PB8 showed the lowest antioxidant ability of all the analyzed PC dimers. Regarding the Fukui indices, the main differences between dimers were established on the base of their susceptibility to electrophilic attack, and the determinant factors were the conformation, together with the type and orientation of the interflavan bond. PB5, which contains the C4β–C6 interflavan bond and acquires a conformationCom in its ground-state, showed greater susceptibility to electrophilic attack, while PB8, which contains the C4α–C6 interflavan bond and acquires a conformationExt in its ground-state, showed the lowest susceptibility to electrophilic attack.



---

*Source: 3535148-2017-02-21.xml*