# p-Hydroxyphenyl-pyranoanthocyanins: An Experimental and Theoretical Investigation of Their Acid-Base Properties and Molecular Interactions

Anna Vallverdu Queralt, Michal Biler, Emmanuelle Meudec, Christine Le Guerneve, Aude Vernhet, Jean Paul Mazauric, Jean-Luc Legras, Michele Loonis, Patrick Trouillas, Veronique Cheynier, et al.

► To cite this version:

Anna Vallverdu Queralt, Michal Biler, Emmanuelle Meudec, Christine Le Guerneve, Aude Vernhet, et al.: p-Hydroxyphenyl-pyranoanthocyanins: An Experimental and Theoretical Investigation of Their Acid-Base Properties and Molecular Interactions. International Journal of Molecular Sciences, 2016, 17 (11), 10.3390/ijms17111842. [https://doi.org/10.3390/ijms17111842](https://doi.org/10.3390/ijms17111842)

HAL Id: hal-01607617

https://hal.science/hal-01607617v1

Submitted on 27 May 2020

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons CC BY 4.0 - Attribution - International License
<figure><img src="image_2.png" /></figure>

Article

*p*-Hydroxyphenyl-pyranoanthocyanins:
An Experimental and Theoretical Investigation of
Their Acid—Base Properties and
Molecular Interactions

Anna Vallverdú-Queralt<sup>1,2</sup>, Michal Biler<sup>3,4</sup>, Emmanuelle Meudec<sup>1</sup>, Christine Le Guernevé<sup>1</sup>, Aude Vernhet<sup>1</sup>, Jean-Paul Mazauric<sup>1</sup>, Jean-Luc Legras<sup>1</sup>, Michèle Loonis<sup>5</sup>, Patrick Trouillas<sup>3,6</sup>, Véronique Cheynier<sup>1,*</sup> and Olivier Dangles<sup>5</sup>

Academic Editor: Nuno Mateus

Received: 26 September 2016; Accepted: 31 October 2016; Published: 5 November 2016

**Abstract:** The physicochemical properties of the wine pigments catechyl-pyranomalvidin-3-O-glucoside (PA1) and guaiacyl-pyranomalvidin-3-O-glucoside (PA2) are extensively revisited using ultraviolet (UV)-visible spectroscopy, dynamic light scattering (DLS) and quantum chemistry density functional theory (DFT) calculations. In mildly acidic aqueous solution, each cationic pigment undergoes regioselective deprotonation to form a single neutral quinonoid base and water addition appears negligible. Above pH = 4, both PA1 and PA2 become prone to aggregation, which is manifested by the slow build-up of broad absorption bands at longer wavelengths (λ ≥ 600 nm), followed in the case of PA2 by precipitation. Some phenolic copigments are able to inhibit aggregation of pyranoanthocyanins (PAs), although at large copigment/PA molar ratios. Thus, chlorogenic acid can dissociate PA1 aggregates while catechin is inactive. With PA2, both chlorogenic acid and catechin are able to prevent precipitation but not self-association. Calculations confirmed that the noncovalent dimerization of PAs is stronger with the neutral base than with the cation and also stronger than π-π stacking of PAs to chlorogenic acid (copigmentation). For each type of complex, the most stable conformation could be obtained. Finally, PA1 can also bind hard metal ions such as Al<fnref n="3" />⁺ and Fe<fnref n="3" />⁺ and the corresponding chelates are less prone to self-association.

**Keywords:** pyranoanthocyanins; acid-base; self-association; copigmentation; metal complexation; UV-visible absorption

<footnote n="1">¹ Institut National de la Recherche Agronomique (INRA), UMR1083 Sciences pour l’œnologie, 2 place Pierre Viala, 34060 Montpellier CEDEX, France; avallverdu@ub.edu (A.V.-Q.); meudec@supagro.inra.fr (E.M.); leguerneve@supagro.inra.fr (C.L.G.); vernhet@supagro.inra.fr (A.V.); mazauric@supagro.inra.fr (J.-P.M.); jean-luc.legras@inra.fr (J.-L.L.)</footnote>

<footnote n="2">² CIBER Fisiopatología de la Obesidad y la Nutrición (CIBEROBN), Instituto de Salud Carlos III, 28903 Madrid, Spain</footnote>

<footnote n="3">³ INSERM UMR 850, University of Limoges, School of Pharmacy, 2 rue du Dr. Marcland, F-87025 Limoges, France; michal.biler@seznam.cz (M.B.); patrick.trouillas@unilim.fr (P.T.)</footnote>

<footnote n="4">⁴ Department of Biophysics, Centre of the Region Haná for Biotechnological and Agricultural Research, Palacký University, Šlechtitelů 27, 78371 Olomouc, Czech Republic</footnote>

<footnote n="5">5 INRA, UMR408 SQPOV, University of Avignon, 84000 Avignon, France; michele.loonis@inra.fr (M.L.); olivier.dangles@univ-avignon.fr (O.D.)</footnote>

<footnote n="6">⁶ Regional Centre of Advanced Technologies and Materials, Department of Physical Chemistry, Faculty of Science, Palacký University Olomouc, tř. 17. listopadu 12, 77146 Olomouc, Czech Republic</footnote>

<footnote>* Correspondence: cheynier@supagro.inra.fr; Tel.: +33-04-9961-2298</footnote>
# 1. Introduction

Among flavonoids, anthocyanins are responsible for plant pigmentation and express a variety of colors from red to blue.

In highly acidic aqueous solution ($pH < 2$), the red flavylium ion $AH^+$ is the only anthocyanin form. As $pH$ increases, water addition at C2 (hydration) becomes thermodynamically favorable, leading to the colorless hemiketal (in equilibrium with low concentration of pale-yellow chalcones). Around $pH 4$, fast deprotonation of $AH^+$ competes with water addition, giving rise to the blue-purple quinonoid base A. Nevertheless, at equilibrium, A is in much lower concentration than the hemiketal, so that the solution is almost colorless. Around $pH 4$, deprotonation of $AH^+$ also takes place, giving rise to blue-purple quinonoid bases (A, kinetic products), which are in low concentration at equilibrium. However, in their natural medium, the anthocyanin colored forms are prone to form noncovalent assemblies, mainly by $\pi-\pi$ stacking interactions, either: (i) with themselves (self-association); or (ii) with the pool of colorless polyphenols especially, flavones, flavonols, and hydroxycinnamic acids (copigmentation) [1–3]. Copigmentation prevents or limits water addition, thereby maintaining higher concentrations of colored forms under mildly acidic to neutral conditions [4,5]. It is typically evidenced by a hyperchromic shift in the visible absorption band.

During wine aging, anthocyanins are progressively converted to specific pigments, in particular pyranomalcyanins with particular chromatic and physicochemical properties. Among them, *p*-hydroxyphenyl-pyranoanthocyanins result from reaction of anthocyanins with either vinylphenols [6] or hydroxycinnamic acids [7], the latter being the precursor of the former by decarboxylation. For instance, catechol and guaiacol are formed by decarboxylation of caffeic and ferulic acids, respectively, and react with malvidin-3-*O*-glucoside to form the catechyl-(PA1) and guaiacyl-(PA2) pyranomalvidin-3-*O*-glucoside derivatives (Figure 1). A detailed description of pyranoanthocyanin formation from malvidin 3-*O*-glucoside and vinylguaiacol has recently been provided, including the full characterization of the reaction medium by a combination of Liquid chromatography coupled to diode array detection and mass spectrometry (LC-DAD-MS) and accurate quantification of the major compounds by quantitative nuclear magnetic ressonance (qNMR) [8].

<figure><img src="image_3.png" /><figcaption>**Figure 1.** Structures of pyranomalvidins pigments catechyl-pyranomalvidin-3-*O*-glucoside (PA1) and guaiacyl-pyranomalvidin-3-*O*-glucoside (PA2).</figcaption></figure>

Controversial results have been reported on the physicochemical properties of pyranomalcyanins. Quijada-Morín et al. [9] investigated the physicochemical and chromatic features
of PA1 and PA2 in the pH range 2.0–4.5. A correct curve-fitting of the visible absorbance vs. pH curves was obtained assuming noncovalent dimerization of PA1 and two successive hydration reaction steps with no proton loss. In contrast, PA2 appeared virtually insensitive to water addition. Consistently, Oliveira et al. [10] proposed that homologues of vitisin A (R = CO₂H, Figure 1) do not undergo hydration in the pH range 0–8 but only deprotonation reactions. Cruz et al. [11] studied four flavanyl-pyranoanthocyanins and two catechyl-pyranoanthocyanins, including PA1, in hydroalcoholic solutions. With these compounds again, only deprotonation, not hydration, was observed. However, a slow and irreversible reaction of vitisin A with water in oxidative conditions leading to neutral pyranone-anthocyanins (oxovitisins) was demonstrated [12]. Remarkably, a complete characterization by NMR of epimeric water adducts of 10-acetyl-pyranoanthocyanins (R = COCH, Figure 1) was reported by Gomez-Alonso et al. [13]. Bakker and Timberlake [14] suggested that vitisins A and B (R = H, Figure 1) undergo water addition. Asenstorfer et al. [15] confirmed hydration of vitisin A but proposed that the competing deprotonation leading to the neutral quinonoid base prevails in mildly acidic solution.

There is also discrepancy about the values reported for the molar absorption coefficients of PAs. Namely, the ε<sub>max</sub> values of PA1 and PA2 at pH 0.8 were reported at 4500 and 3060 M<sup>-1</sup>·cm<sup>-1</sup>, respectively [9], whereas an earlier estimation for PA1 at pH 1.5 gave 12,200 M<sup>-1</sup>·cm<sup>-1</sup> [16].

The aim of this work is to revisit the chromatic and physicochemical properties of *p*-hydroxyphenyl-pyranoanthocyanins (PA1 and PA2) with a specific emphasis on: (a) proton transfer vs. hydration reactions; and (b) self-association vs. copigmentation. The experimental data by UV-visible spectroscopy and dynamic light scattering (DLS) are extensively analyzed and discussed with help of quantum chemical density functional theory (DFT) calculations.

## 2. Results and Discussion

### 2.1. Characterization of PA1 and PA2

After the purification steps, PA1 and PA2 were individually analyzed by LC-DAD/ESI-MS in positive-ion mode. PA1 (λ<sub>max</sub> = 506 nm) and PA2 (λ<sub>max</sub> = 506 nm) were detected at *m/z* 625 and *m/z* 639, respectively. Full ¹H and ¹³C NMR characterization confirmed the structures of both pigments (Table S1 and Table S2).

### 2.2. Structural Transformations

#### 2.2.1. UV-Visible Spectroscopy

The ε<sub>max</sub> values in acidified MeOH (0.1 HCl) were measured at 18 × 10³ and 42.5 × 10³ M<sup>-1</sup>·cm<sup>-1</sup> (λ<sub>max</sub> = 520 nm) for PA1 and PA2, respectively. Those ε<sub>max</sub> values are much higher than those previously reported in the literature, thereby attesting the high purity of the PA samples used in this work.

Small volumes of concentrated PA solutions in acidified MeOH (100% flavylium) were added to buffered aqueous solutions in the pH range 1–6 and the spectral properties were monitored over 15–30 min. Under such conditions, common anthocyanins undergo a typical color loss, reflecting water addition onto AH⁺ [17,18] during the period of time required to reach pseudo-equilibrium (without taking into account the slow *cis*-to-*trans* chalcone isomerization), i.e., from a few seconds (pH 2–3) to several minutes (pH 5–6). The apparent kinetics of hydration is largely governed by the fraction of AH⁺, the sole colored form to undergo water addition in solution. Unlike common anthocyanins, PA1 and PA2 do not clearly undergo such relatively fast color loss in agreement with previous investigations on various PAs [11,19]. Thus far, the sole PAs whose cationic form has been clearly shown to be electrophilic enough to add water (at C10) are the acetylpyranoanthocyanins (stemming from the addition of diacetyl-2,3-butanedione to anthocyanins) owing to the electron-withdrawing effect of the acetyl group [20]. By contrast, the visible spectra of PA1 and PA2 in the pH range 2.0–3.5 are essentially stable over time and identical to the spectra under highly acidic conditions (pH = 1), in which AH⁺
is the sole species. Above pH = 3.5, an instantaneous broadening and weakening of the visible band becomes observable with both PAs, which reflects proton loss. Therefore, a few seconds after mixing, the pH dependence of the visible absorbance at 505 nm (maximal absorption wavelength of AH⁺) can be safely analyzed with the assumption of a single proton loss. In an acetate buffer containing 12% EtOH, the first pKₐ of PA1 (concentration: 50 µM) and PA2 (ca. 25 µM) was estimated at 4.30 (±0.09) and 4.35 (±0.11), respectively. In the citrate-phosphate buffer (no EtOH), the first pKₐ of PA1 was estimated at 4.69 (±0.06) (Figure 2A). These values are consistent with those from the literature on PAs [11,18] and common anthocyanins [18].

<figure><img src="image_4.png" /><figcaption>**Figure 2.** (A) Plot of the visible absorbance of PA1 (48.4 µM, citrate-phosphate buffer) at 505 nm (collected immediately after addition of PA1 to buffer) as a function of pH. The solid line is the result of the curve-fitting assuming a single proton transfer: pKₐ = 4.69 (±0.06), rA = εA/εAH⁺ = 0.47 (±0.02), r = 0.993; (B) normalized spectra of pure cation (pH 1.88) and pure base (spectrum at pH 4.64 after subtraction of ca. 50% of the residual cation).</figcaption></figure>

The spectrum of the pure neutral base A (Figure 2B) was deduced from the spectrum obtained at pH ≈ pKₐ (50% flavylium + 50% base) and the spectrum of the pure AH⁺ form (pH ≈ 2). This deduced spectrum was used as the experimental control for calculations aimed at determining the distribution of the different tautomers (deprotonation of the acidic OH groups at C4B-OH, C4E-OH or C7A-OH) and their individual visible spectra (see below). Surprisingly, the experimental spectrum of A displays a main absorption band at λmax ≈ 480 nm, which means a hypsochromic shift of more than 20 nm with respect to AH⁺, and a shoulder in the range 530–540 nm. Its maximum molar absorption coefficient is circa half that of AH⁺. This uncommon hypsochromic shift is consistent with previous observations [11,19]. This could explain why hydration (formation of the colorless forms) and proton loss are not easily distinguished with PAs in absence of a careful kinetic investigation.

### 2.2.2. DFT-Based Rationalization

The AH⁺ form of PA1 exhibits a large band at 506 nm (Figure 2B), which is predicted by the TD-DFT calculation as corresponding to the first excited state S0 → S1 (at 496.3 nm, see Figure 3A).

This excited state is mainly described by the HOMO → LUMO electronic transition, both molecular orbitals (MO) being delocalized along an extended conjugated path (Figure 4A).
<figure><img src="image_5.png" /></figure>

<figure><img src="image_6.png" /></figure>

<figure><img src="image_7.png" /></figure>

<figure><img src="image_8.png" /><figcaption>**Figure 3.** Experimental and theoretical (TD-DFT) UV-visible absorption spectra of PA1 for: (A) AH⁺; and the three A forms resulting from deprotonation of: C4E-OH (B, base 1); C4B-OH (C, base 2); and C7A-OH (D, base 3).</figcaption></figure>

<figure><img src="image_9.png" /><figcaption>**Figure 4.** Molecular orbitals (MO) correlation diagram of PA1 for: (A) AH⁺; and (B) the A-form resulting from deprotonation of C7A-OH (base 3).</figcaption></figure>
The shoulder experimentally observed at 430 nm is ascribed to the second (S0 → S2) excited state theoretically lying at 446.4 nm (Figure 3A), which is mainly described by the HOMO-1 → LUMO electronic transition and mainly involving the E-ring, i.e., catechol moiety in the case of PA1 (Figure 4A). When increasing pH, the UV-visible absorption properties are modified, highlighting a clear hypsochromic shift (Figure 2B). The excited states of the deprotonated forms were evaluated for the three possible neutral quinonoid bases A. Concerning the first two deprotonated forms, i.e., proton loss from C4E-OH (base 1) or C4B-OH (base 2), a clear bathochromic shift is predicted for the first excited state (Figure 3B,C). In both cases, this excited state is still described by the HOMO → LUMO electronic transition. Whereas the LUMO is almost identical for AH⁺ and bases 1 and 2, deprotonation of C4E-OH or C4B-OH induces a destabilization of the HOMO (Figure S1A,B), therefore decreasing the energy gap and resulting in a bathochromic shift. By contrast, upon deprotonation of C7A-OH (base 3), the HOMO is dramatically displaced toward the A-ring, therefore strongly decreasing the overlap with the LUMO. As a consequence, the first excited state (S0 → S1) at 517.9 nm has relatively low oscillator strength compared to AH⁺ (0.2 and 0.9, respectively) and is slightly bathochromically shifted (as for bases 1 and 2) (Figure 3D). However, HOMO-1 is moved toward the C- and D-rings (Figure 4B), thus increasing the overlap with the LUMO and the oscillator strength of the S0 → S2 excited state at 454.1 nm. In conclusion, the observed hypsochromic shift is not an actual shift but an inversion in intensity between the first and second excited states. This is confirmed by the presence of a shoulder experimentally observed at 540 nm (Figure 3D). The agreement between the theoretical and experimental spectra supports the presence in solution of only one neutral quinonoid base corresponding to proton loss from C7A-OH. This is confirmed by the relative free enthalpies of the three quinonoid bases, showing that base 3 is energetically preferred by ca. 3 kcal·mol⁻¹ (Table 1).

<table><caption>**Table 1.** Relative free energies (kcal·mol⁻¹) of the quinonoid bases for pigments catechyl-pyranomalvidin-3-O-glucoside (PA1) and guaiacyl-pyranomalvidin-3-O-glucoside (PA2).</caption><thead><tr><th>Pyranoanthocyanin Form</th><th>PA1</th><th>PA2</th></tr></thead><tbody><tr><td>Base 1</td><td>4.28</td><td>3.60</td></tr><tr><td>Base 2</td><td>4.31</td><td>3.14</td></tr><tr><td>Base 3</td><td>0</td><td>0</td></tr></tbody></table>

According to the Boltzmann law, each minor quinonoid base (bases 1 and 2) represents only 1% of the total; their presence most probably contributes only slightly to the shoulder observed at 540 nm. The same rationalization also applies to PA2 (see Table S3 and Table S4).

## 2.3. Noncovalent Assemblies of Pyranoanthocyanins

### 2.3.1. Self-Association

Self-association of anthocyanins has been investigated by UV-visible spectroscopy, circular dichroism and NMR [21–23]. At high concentration (>0.1 mM), anthocyanins typically form left-handed aggregates, which results in relatively modest modifications in the visible spectra (band broadening with a slight decrease of absorption at $\lambda_{max}$). Interestingly, investigations by ¹H-NMR in acidic to mildly alkaline solutions [21,22] gave the following values for the dimerization constants of the colored forms of malvidin 3,5-diglucoside: AH⁺: 600, A: 900, A⁻: 200 M⁻¹. This means that, as expected, dimerization is stronger for the neutral form (pH 4–7), i.e., when electrostatic repulsion is minimized between the two chromophores.

By contrast, little data is available about the self-association of pyranoanthocyanins. One study was reported with the pyranoanthocyanin formed by addition of 8-vinylcatechin (itself a product of the aromatic electrophilic substitution of acetaldehyde with catechin) and malvidin 3-O-β-D-(6-O-p-coumaroyl)glucoside. At pH = 1.3, the AH⁺ form of this PA was shown to form a noncovalent dimer with a binding constant Kd of $3 \times 10^3$ M⁻¹ [11], which is only circa twice higher
than for the parent anthocyanin (estimated in the presence of 30% DMSO and thus probably lower than its value in aqueous solution) [23]. Although the $AH^+$ forms of PAs do not seem much more prone to self-association than the parent anthocyanins, under mildly acidic conditions, $AH^+...A$ or $A...A$ self-association is likely to be stronger due to the decrease in electrostatic repulsion.

In this work, a surprising observation is the slow decay of the main visible band taking place in 50 µM solutions of PA1 at pH > 3.5 (Figures 5 and 6).

<figure><img src="image_10.png" /><figcaption>**Figure 5.** (A) Changes in visible absorbance (Abs) following the addition of PA1 (48.4 µM) to a pH 4.39 acetate buffer + 12% EtOH; and (B) UV-visible spectra at different time points.</figcaption></figure>

<figure><img src="image_11.png" /><figcaption>**Figure 6.** (A) Changes in visible absorbance following the addition of PA1 (48.4 µM) to a pH 5.20 acetate buffer + 12% EtOH; and (B) UV-visible spectra at different time points.</figcaption></figure>

These spectral changes are not consistent with hydration, being too slow, non-first-order, faster as the pH increases, and accompanied by the build-up of broad absorption bands at longer wavelengths ($\lambda \ge 600$ nm). To the best of our knowledge, such observations have never been reported before. A reasonable explanation for the slow changes observed in the visible spectra of PA1 at pH > 3.5 is that self-association involving the neutral quinonoid base occurs. Mixed aggregates $AH^+...A$ may form first, followed by aggregates of pure base A at higher pH, which is consistent with the spectral differences between pH 4.4 (Figure 5) and pH 5.2 (Figure 6). The additional methyl group should make PA2 (25 µM) more prone to self-association and precipitation than PA1 (50 µM). As a direct consequence of more favored noncovalent interactions, precipitation rapidly occurred in PA2 solutions at pH > 4 (even in presence of 12% EtOH), despite relatively low concentration. By contrast, PA1 solutions remained free of precipitate over hours in the presence of 12% EtOH. Similar observations were made in dilute (10 µM) PA2 solutions. Unlike PA1, aggregation in dilute PA2 solution, which occurred relatively rapidly at pH > 4.5, was mainly manifested by a drop of visible absorbance and low-energy absorption bands were much weaker than with PA1 (Figure 7).
<figure><img src="image_12.png" /><figcaption>**Figure 7.** Visible spectra of PA2 (10 µM) in a pH 5.0 acetate buffer + 12% EtOH at different time (t) points.</figcaption></figure>

The higher propensity of A for self-association is corroborated by the concentration dependence of the first pK<sub>a</sub> value. Indeed, dilution by a factor of 2.5 shifted pK<sub>a</sub> value from 4.35 (±0.11) to 4.79 (±0.12) for PA2. This means that at high concentration, proton loss occurs at lower pH (acid-base equilibrium shifted toward the base) in agreement with A being more prone to self-association than AH<sup>+</sup>.

DFT calculations including dispersion correction confirmed the capacity of PA1 and PA2 to self-associate, both exhibiting negative energies of association (Table 2 and Table S3 for PA1 and PA2, respectively).

<table><caption>**Table 2.** Self-association energy (kcal·mol<sup>-1</sup>) of PA1, according to the charge state (A···A, AH<sup>+</sup>···A and AH<sup>+</sup>···AH<sup>+</sup>) of the two partners and according to the geometries of the different association complexes gathered in terms of structural similarities as three different groups (A, B and C). The association energies of the most stable conformers (for the three charged forms) are highlighted in bold.</caption><thead><tr><th>Group</th><th>Geometry</th><th>A···A</th><th>AH<sup>+</sup>···A</th><th>AH<sup>+</sup>···AH<sup>+</sup></th></tr></thead><tbody><tr><td rowspan="7">A</td><td>1</td><td>-28.0</td><td>-22.3</td><td>-10.1</td></tr><tr><td>2</td><td>-23.2</td><td>-25.3</td><td>-25.4</td></tr><tr><td>3</td><td>-22.8</td><td>-22.2</td><td>-11.4</td></tr><tr><td>4</td><td>-19.2</td><td>-19.9</td><td>-9.9</td></tr><tr><td>5</td><td>-18.2</td><td>-22.4</td><td>-14.9</td></tr><tr><td>6</td><td>-14.9</td><td>-15.9</td><td>-9.0</td></tr><tr><td>7</td><td>-12.7</td><td>-13.4</td><td>-7.9</td></tr><tr><td rowspan="2">B</td><td>8</td><td>-15.9</td><td>-21.6</td><td>-8.6</td></tr><tr><td>9</td><td>-15.6</td><td>-15.3</td><td>-9.5</td></tr><tr><td rowspan="2">C</td><td>10</td><td>-17.9</td><td>-8.5</td><td>-6.0</td></tr><tr><td>11</td><td>-15.2</td><td>-14.0</td><td>-10.0</td></tr></tbody></table>

The three possible self-associations (AH<sup>+</sup>···AH<sup>+</sup>, AH<sup>+</sup>···A, and A···A) appeared likely to occur. For each, several conformers were obtained (see Figure 8 and Table 2 for geometrical description and stabilizing energies, respectively), in which π-π stacking is the driving force of association.
<figure><img src="image_13.png" /><figcaption>**Figure 8.** Three representative conformers of PA1 self-association complexes of the three representative groups, exhibiting geometrical similarities: Group (A) interactions between A-, C-, and D-rings of both PA1 fragments; Group (B) interactions between A-, C-, and D-rings of one PA1 and the B-ring of the other PA1; and Group (C) interactions between A-, C-, and D-rings of one PA1 with the E-ring of the other PA1. For clarity purpose, one PA1 partner is in blue, whereas the other is both in red (pyranoanthocyanidin moiety) and in green (sugar moiety).</figcaption></figure>

Three types of conformers could be distinguished: Group A, exhibiting interactions between A-, C-, and D-rings of both PA1 units; Group B, exhibiting interactions between A-, C-, and D-rings of one PA1 and the B-ring of the other PA1; Group C, exhibiting interactions between A-, C-, and D-rings of one PA1 with the E-ring of the other PA1. In each group, electrostatic and H-bonding interactions influence orientation of both units with respect to each other. In particular, the strong electrostatic repulsion in $AH^+...AH^+$ complexes strongly distorted the geometries (compared to $AH^+...A$ and $A...A$), however without preventing the stabilizing $\pi-\pi$ stacking forces.

Self-association is here studied in terms of electronic energy, not Gibbs energy, mainly because entropy related to self-association cannot be accurately calculated [3]. Hence, if the most stable conformer can be regarded as the major conformer present in solution, all other conformers must also be considered, although with relative weight obtained from their relative stabilizing energy. Having this in mind, Table 2 shows that $A...A$ association is favored as the most stable $A...A$ conformer is more stable by ca. $3 \text{ kcal} \cdot \text{mol}^{-1}$ than the most stable $AH^+...A$ and $AH^+...AH^+$ conformers. Moreover, stable $A...A$ conformers appear more numerous, in particular when comparing to $AH^+...AH^+$ dimers. In this latter case indeed, only one $AH^+...AH^+$ dimer exhibits a stabilizing energy close to those of the $AH^+...A$ dimers. Hence, DFT calculations including dispersion forces confirm that $AH^+...AH^+$ ($\text{pH} < pK_a$) forms to a lesser extent than $AH^+...A$ ($\text{pH} \approx pK_a$) and $A...A$ ($\text{pH} > pK_a$), in agreement with the experimental observations.

The observation that PA2 is more prone to self-association and precipitation than PA1 was also confirmed by the dispersion corrected (DFT-D) calculations, as the most stable $A...A$ conformers exhibited association energies of -32.3 and -28.0 kcal·mol⁻¹, for PA2 and PA1, respectively (Table 2 and Table S3). As expected, replacing an OH group by an OMe group on an aromatic ring increases $\pi-\pi$ stacking interactions by slightly extending the $\pi$-conjugated system [3].

Noncovalent dimerization must be instantaneous (diffusion-controlled), which is in contrast to the slow spectral changes observed with PAs. It is thus very likely that PA self-association is not limited to the formation of noncovalent dimers or small-sized oligomers but leads to large aggregates. This was confirmed by DLS experiments. For PA1 (50 µM), an immediate increase in the relative scattering intensity ($I/I_0$) was observed at pH 4.4 (Figure 9A) and 5.2, which was related to the immediate formation of polydispersed aggregates (PI = 0.4) with an average hydrodynamic diameter (Dh) of ca. 150 nm. With such polydisperse dispersions, it must be kept in mind that the Dh value is only indicative and in favor of the largest particles. The average Dh and polydispersity index (PI) kept
increasing during the experiment, up to the formation of micron-sized aggregates and to a PI of 1, whereas I/I₀ reached a plateau-value after about 3–4 h. This indicated that no new aggregates formed at this stage. After 5 h, the decrease in I/I₀ was related to precipitation of micron-sized aggregates (Figure 9A). At pH = 3.5, there was no change in scattering intensity, meaning that no aggregation occurred. For PA2 (25 µM), precipitation rapidly occurred after mixing while interaction kinetic was significantly slowed down at lower concentration (10 µM).

<figure><img src="image_14.png" /></figure>

<figure><img src="image_15.png" /></figure>

**Figure 9.** Aggregation kinetics followed by dynamic light scattering (DLS) (observation angle = 90°). A rise in the normalized scattered intensity I/I₀ above 1 indicates the onset of aggregation. I₀, intensity scattered by the solvent; I, intensity scattered after PA1 addition. Dh: average hydrodynamic diameter of the aggregates, weighted according to the intensity. Y axis indicates I/I₀ or Dh value. (A) PA1 (50 µM) at pH 4.4 (+12% EtOH); and (B) PA1 (50 µM) with chlorogenic acid (100 equiv.) at pH 4.4 (+12% EtOH).

### 2.3.2. Copigmentation

Adding a large excess of chlorogenic acid to PA1 aggregates (50 µM PA1 solution ca. 30 min after preparation) led to conversion of the visible spectrum back to that of a fresh solution (i.e., before the onset of aggregation) (Figure 10).

<figure><img src="image_16.png" /></figure>

**Figure 10.** Visible spectra of PA1 (48.4 µM) in a pH 4.62 acetate buffer + 12% EtOH: immediately after addition (t = 0), after the development of self-association (t = 27.5), and after the subsequent addition of the copigment chlorogenic acid (Copigment, 100 equiv).
This suggests the dissociation of the PA1 aggregates due to the competitive formation of noncovalent PA1···chlorogenic acid complexes. A bathochromic shift of the $λ_{max}$ value by ca. 10 nm was observed, which is typical of copigmentation. Copigmentation is actually known to compete with anthocyanin self-association. Moreover, the time-dependent spectral changes that were observed due to PA1 self-association at pH value of 4.6 were markedly slowed down when the chlorogenic acid concentration was increased and totally abolished with 100 equiv. copigment (data not shown). DLS experiments confirmed that addition of chlorogenic acid (100 equiv.) inhibited self-association of PA1 and prevented the formation of large aggregates prone to precipitation (in acetate buffer containing 12% ethanol at pH = 4.4 and 5.2): the scattering intensity remained unchanged and only a slight Dh increase from 100 to 200 nm was observed (Figure 9B). Unlike chlorogenic acid, catechin was unable to reverse PA1 self-association in agreement with catechin having less affinity for anthocyanins than chlorogenic acid [3].

With PA2 (25 µM), both chlorogenic acid and catechin (20–100 equiv.) were able to prevent precipitation at pH > 4 but not self-association. As said above, under more dilute conditions (10 µM), precipitation is not observed, even in the absence of copigment, and self-association is slow at pH = 4.4. However, at pH ≥ 5, self-association is much faster (Figure 7) and its kinetics is not significantly influenced by any copigment. The DLS experiments showed that addition of chlorogenic acid to PA2 (pH = 4.4 and pH ≥ 5) did not induce any change in the scattering intensity, but the aggregate size kept increasing, meaning that chlorogenic acid can indeed inhibit precipitation but not self-association (data not shown).

Theoretical calculations confirmed the formation of copigmentation complexes between chlorogenic acid and both PA1 and PA2. For PA1, complexation with $AH^+$ appeared less favorable than with A (deprotonation of C7A-OH). For instance, the stabilizing energies of the most stable PA1···chlorogenic acid conformers were −18.0 and −24.5 kcal·mol⁻¹ for $AH^+$ and A, respectively (Table 3). This difference is weaker (1.3 kcal·mol⁻¹) with PA2 (Table S3).

<table><caption>**Table 3.** Association energies (kcal·mol⁻¹) of PA1 (in both $AH^+$ and A forms) with chlorogenic acid. Groups A, B, C and D are defined according to geometry similarity (see text and Figure 11). The association energies of the most stable conformers are highlighted in bold.</caption><thead><tr><th colspan="3">AH<sup>+</sup>···Chlorogenic Acid</th><th colspan="3">A···Chlorogenic Acid</th></tr><tr><th>Group</th><th>Geometry</th><th>Binding Energy</th><th>Group</th><th>Geometry</th><th>Binding Energy</th></tr></thead><tbody><tr><td rowspan="3">A</td><td>1</td><td>−13.2</td><td rowspan="3">A</td><td>1</td><td>−22.0</td></tr><tr><td>2</td><td>−13.2</td><td>2</td><td>−17.7</td></tr><tr><td>3</td><td>−12.5</td><td>3</td><td>−14.5</td></tr><tr><td rowspan="2">B</td><td>4</td><td><strong>−18.0</strong></td><td>B</td><td>4</td><td><strong>−24.5</strong></td></tr><tr><td>5</td><td>−12.9</td><td>C</td><td>5</td><td>−18.4</td></tr><tr><td>C</td><td>6</td><td>−12.9</td><td></td><td>6</td><td>−14.6</td></tr><tr><td>D</td><td>7</td><td>−12.9</td><td>D</td><td>7</td><td>−13.1</td></tr></tbody></table>

Copigmentation association is mainly driven by $π-π$ stacking with many possible orientations of chlorogenic acid with respect to pigment. Four groups can be distinguished in which the catechol moiety of chlorogenic acid binds the A-, C- and D-rings (Group A), the B-ring (Group B), the E-ring (Group C) or the sugar moiety (Group D) of PA (Figure 11, Figure S2 and Figure S3).

In Group D, H-bonding would mainly drive copigmentation; as expected, the corresponding conformers are less stable than those of Groups A, B, and C. In Groups A, B, and C, $π-π$ stacking is the driving force, intermolecular electrostatic interactions and H-bonding driving only the respective orientation of pigment and copigment. Copigmentation of PAs with chlorogenic acid appears less favorable than self-association. Indeed, when comparing geometries of noncovalent dimers and copigmentation complexes, $π-π$ stacking in noncovalent dimers is much more favored due to more extended $π$-conjugated systems facing each other. This is consistent with the observation that
chlorogenic acid has to be used in large excess (vs. pyranoanthocyanins) to prevent aggregation. Aggregation involves the formation of dimers and higher oligomers of high conformational variability. At high concentration, copigments can intercalate in small aggregates, thus slowing down or even inhibiting the formation of larger aggregates. Copigments can also intercalate in large aggregates and promote their dissociation.

<figure><img src="image_17.png" /><figcaption>**Figure 11.** Four representative conformers of copigmentation complexes between PA1 (blue) and chlorogenic acid (red) of the four representative groups, exhibiting geometrical similarities. The catechol moiety of chlorogenic acid is in interaction with the following moieties of PA1: Group (A) A-, C- and D-rings; Group (B) B-ring; Group (C) E-ring; and Group (D) glucose.</figcaption></figure>

### 2.3.3. Metal Complexation

Owing to its catechol group, PA1 is able to bind hard metal ions such as Al³⁺ and Fe³⁺. Under mildly acidic conditions, Al³⁺ binding and concomitant loss of the two catechol protons result in a strong bathochromic shift and a decrease in visible absorbance at λmax (Figure 12A). With Fe³⁺, a strong band broadening is also observed, which features additional ligand-to-metal charge transfer (Figure 12B).

<figure><img src="image_18.png" /><figcaption>**Figure 12.** Visible spectra of PA1 (48.4 µM) in a pH 4.62 acetate buffer + 12% EtOH: immediately before metal addition and after metal addition (10 equiv.) at two time points: (A) Al³⁺; and (B) Fe³⁺.</figcaption></figure>
Both bindings are fast, being completed within 10–20 s. It is noteworthy that the metal chelates thus formed seem less prone to self-association (a possible consequence of charge repulsion) as the spectral changes observed over 20–30 min are less intense than with metal-free PA1.

## 2.4. Spectral Shifts Related to Noncovalent Association

The spectral shifts characteristic of self-association and copigmentation were studied by TD-DFT calculations and analyzed according to the $\pi \rightarrow \pi^*$ electronic transitions constituting the first excited states (corresponding to $\lambda_{max}$). Only the spectral characteristics of the $AH^+$···chlorogenic acid and $A···A$ complexes will be detailed. In the former case, $\lambda_{max}$ is bathochromically shifted in agreement with the experiment. The low-energy band is assigned to the first excited state ($S0 \rightarrow S1$), which is mainly described by HOMO → LUMO electronic transition. As already observed [3], the spectral behavior of these noncovalent complexes is highly dependent on the conformer. In some cases, the charge transfer (CT) character is strong; in some others, it is very weak (Table 4 and Figure 13).

<table><caption>**Table 4.** Bathochromic shifts $\Delta\lambda$ in nm, molecular orbital (MO) descriptions (%) and charge transfer (CT) character of the PA1 ($AH^+$)···chlorogenic acid complex (same geometry numbering as in Table 3A).</caption><thead><tr><th>Geometry</th><th>Δλ</th><th>MO</th><th>CT</th></tr></thead><tbody><tr><td>1</td><td>8.03</td><td>HOMO → LUMO (0.84)</td><td>No</td></tr><tr><td>2</td><td>8.60</td><td>HOMO → LUMO (0.83)</td><td>Weak</td></tr><tr><td>3</td><td>12.03</td><td>HOMO → LUMO (0.83)</td><td>Yes</td></tr><tr><td>4</td><td>4.45</td><td>HOMO → LUMO (0.44)</td><td>Yes</td></tr><tr><td>5</td><td>0.04</td><td>HOMO → LUMO (0.63)</td><td>Yes</td></tr><tr><td>6</td><td>16.79</td><td>HOMO → LUMO (0.87)</td><td>No</td></tr><tr><td>7</td><td>7.86</td><td>HOMO → LUMO (0.85)</td><td>Yes</td></tr></tbody></table>

<figure><img src="image_19.png" /><figcaption>**Figure 13.** Molecular orbital description of the copigmentation complex of PA1 (flavylium ion) with chlorogenic acid: without (a); or with (b) CT character in the excited state.</figcaption></figure>

When CT is significant, the HOMO is on the A-, B-, C-, and E-rings of the pigment, whereas the LUMO is delocalized on the copigment (Figure 13B). The average value of the spectral shift is ca. 10 nm (Table 4). Copigmentation of the quinonoid base also induces significant bathochromic shifts, which can be even higher than for the cation (Table S5).

Self-association also leads to global bathochromic shifts. However in this case, the MO analysis is more complex. Indeed, with the $A···A$ dimer, several excited states must be considered, most of them
leading to bathochromic shifts with respect to the most intense band of the monomer. For instance, for the most stable conformer, the first four excited states correspond to bathochromic shifts of 11.5, 26.4, 41.9, and 57.9 nm, for S0 → S4 to S0 → S1, respectively (Table 5).

<table><caption>**Table 5.** Absorption wavelength λ in nm (oscillator strengths in brackets) of the first four excited states (S₀ → S₁) of the A···A PA1 dimer (same geometry numbering as in Table 2).</caption><thead><tr><th>States</th><th>1 *</th><th>2</th><th>3</th><th>4</th><th>5</th></tr></thead><tbody><tr><td>S₀ → S₁</td><td>451.2 (0.155)</td><td>462.1 (0.028)</td><td>452.9 (0.072)</td><td>454.1 (0.128)</td><td>468.3 (0.019)</td></tr><tr><td>S₀ → S₂</td><td>435.2 (0.093)</td><td>433.6 (0.469)</td><td>449.1 (0.146)</td><td>450.9 (0.047)</td><td>449.3 (0.301)</td></tr><tr><td>S₀ → S₃</td><td>419.7 (0.495)</td><td>410.7 (0.019)</td><td>411.3 (0.358)</td><td>420.5 (0.649)</td><td>410.1 (0.045)</td></tr><tr><td>S₀ → S₄</td><td>405.1 (1.052)</td><td>396.3 (0.213)</td><td>399.1 (0.636)</td><td>399.3 (0.747)</td><td>393.5 (0.861)</td></tr></tbody></table>

This leads to four visible absorption bands of different intensity and CT character (Figure 14).

<figure><img src="image_20.png" /><figcaption>**Figure 14.** Molecular orbital description of the A···A PA1 dimer: (A) Geometry 1 exhibiting slight CT in the first excited state (S0 → S1); and (B) Geometry 4 exhibiting strong CT in the fourth excited state (S0 → S4).</figcaption></figure>

Considering the multitude of conformers, the noncovalent dimer is optically characterized by a multitude of bathochromically shifted absorption bands. In the situation of small n-mers or large aggregates, this effect is expected to increase, i.e., multiplication of conformers and of bathochromically shifted excited states. This is consistent with the very broad low-energy visible absorption bands gradually appearing in the experimental PA spectra recorded in mildly acidic conditions (Figures 5 and 6).

### 3. Materials and Methods

#### 3.1. *Chemicals*

All solvents were of HPLC quality and all chemicals of p.a. grade. Vinylguaiacol and glycerol were purchased from Sigma (Saint-Louis, MO, USA). Malvidin-3-O-glucoside chloride, *p*-coumaric, chlorogenic and caffeic acids were purchased from Extrasynthèse (Genay, France). Purified deionized water (MilliQ purification system, Millipore, Molsheim, France) was used for the preparation of all solutions.

#### 3.2. *Production of Vinylcatechol*

Vinylcatechol was obtained from caffeic acid by bioconversion using the yeast strain *Saccharomyces cerevisiae* MTF3686 (internal collection). A two liters cell culture was prepared as described earlier [24]

<footnote>* The spectral shifts stated in the main text are considered with respect to the absorption wavelength of A, which is 393.3 nm (oscillator strength of 0.856).</footnote>
with some modifications: glycerol 10 g/L was used as a carbon source for the YNB media, and 0.5 mM of *p*-coumaric acid was added to induce the decarboxylase activity. Cells were recovered in exponential growth phase, washed twice with 50 mM citrate-phosphate buffer (pH = 3.5). Highly dense suspensions of resting cells were prepared in 1 L of citrate-phosphate buffer (pH = 3.5) containing 2% glucose and supplemented with 0.5 mM of caffeic acid as a bioconversion substrate. The medium was agitated at 27 °C during 4 days, then centrifuged (1000 × g) for 10 min. The supernatant was passed through a Sep-Pak® tC18 cartridge (10 g reversed-phase sorbent, Waters, Milford, MA, USA), eluted successively with H₂O and MeOH. The MeOH fraction containing vinylcatechol was mixed with water and concentrated under vacuum to remove MeOH, analysed by HPLC-DAD as described earlier [8] to determine the concentration as equivalent vinylguaiacol, and used for the synthesis immediately.

### 3.3. *Pyranoanthocyanin Synthesis*

Malvidin-3-O-glucoside (2 mM) and vinylguaiacol (2 mM) or vinylcatechol (2 mM equiv. vinylguaiacol) were prepared in MeOH/H₂O mixture (70/30). The pH of the mixture was adjusted with concentrated HCl to pH 3. The reactions were carried out during 5 days at 35 °C [8]. The reaction media were concentrated and re-dissolved in 1% HCOOH in H₂O.

### 3.4. *Purification of PA1 and PA2*

Semi-preparative purification of PA1 and PA2 was performed using a Varian LC 940 system (Agilent Technologies, Santa Clara, CA, USA) with a UV detector and fraction collector. After filtration on a 0.45 µm HVLP Millipore filter, the solution was injected onto a Microsorb 100—3 µm, C18 column (100 × 21.4 mm, Dynamax, Santa Clara, CA, USA). Elution was achieved at a flow rate of 10 mL/min with solvents A (0.1% aqueous HCl) and B (0.1% HCl in MeCN-H₂O (4:1)) with the following gradient: 0–6 min, linear 39%–41% B, 6–8 min, linear 41%–80% B, 8–9 min, isocratic 80% B, 9–10 min, linear 80%–39% B. The collected fractions were freeze-dried and stored at –80 °C until analysis.

### 3.5. *NMR Analysis*

All experiments were performed on an Agilent 500 MHz DD2 NMR spectrometer (Agilent Technologies, Santa Clara, CA, USA) operating at 500.05 MHz and 125.75 MHz for ¹H and ¹³C, respectively, and equipped with a 5 mm indirect detection Z-gradient probe. The chemical shifts were expressed with MeOH as the internal reference (δ = 3.31 and 49.1 ppm for ¹H and ¹³C, respectively). All data were collected with active sample temperature regulation. Assignments of proton and carbon signals were deduced from homonuclear ¹H correlation spectroscopy, rotating-frame overhauser spectroscopy, diffusion-ordered spectroscopy and heteronuclear ¹H/¹³C heteronuclear single quantum coherence and heteronuclear multiple bond correlation experiments. Data were processed and analyzed using the VNMRJ software (Agilent, Santa Clara, CA, USA).

### 3.6. *UV-Visible Spectroscopy*

Spectra were recorded with an Agilent 8453 diode array spectrophotometer fitted with a quartz cell (optical path length = 1 cm) equipped with a stirring magnet. A constant temperature of 25 °C in the cell was obtained by use of a water-thermostated bath.

Small volumes of concentrated PA solutions in acidified MeOH (0.1 M HCl, 100% flavylium) were added to a 0.2 M acetate buffer containing 12% EtOH or to a citrate-phosphate buffer (pH range 1–6). In the cell, the pigment concentration was ca. 50 µM for PA1 and 25 or 10 µM for PA2. The spectral monitoring was carried out over 15–30 min. Then, the samples were collected and kept overnight at room temperature to repeat the spectroscopic analysis. The experiments were also repeated with buffers containing a copigment (CP: chlorogenic acid or catechin) in large excess (20–100 equiv.) or a metal ion (PA1 + 10 equiv. AlCl₃ or Fe(NO₃)₃) in acetate buffer + 12% EtOH.
Alternatively, some samples without CP were combined at the end of the kinetic monitoring and the pH was adjusted to 4.5 or 5.2. To the resulting solutions, the copigment was added (CP-to-PA molar ratios = 20–100) and the UV-visible spectra were recorded.

### 3.7. Dynamic Light Scattering (DLS)

DLS experiments were conducted with a Malvern Autosizer 4800 HS (Malvern Instruments, Malvern, UK) equipped with a 30 mW He-Ne laser (λ = 633 nm) and APD detection. Measurements were carried out at an angle of 90° from the incident beam. Prior to experiments, all solvents were filtered on 0.2 µm microfiltration units (Millipore, Millex-GV, Molsheim France). Both PA1 and PA2 were prepared in 0.1 M HCl in MeOH. PA solutions (10 µL) were diluted in acetate buffer (pH = 3.5, 4.4 and 5.2, 2 mL) and introduced in the measurement cell at 25 °C. Each experiment was carried out in duplicate. Control measurements were performed on solvents. Addition of chlorogenic acid (100 equiv.) to PA1 and PA2 was performed in acetate buffer containing 12% ethanol at pH = 3.5, 4.4 and 5.2.

Sample evolution was followed by measuring the scattered intensity and the autocorrelation function G(t) of the scattered light. Each measurement was the average of 10 subruns. A rise in scattered intensity above control values indicated the onset of aggregation.

The autocorrelation function of the scattered light was first analyzed using the cumulant method, which gives an average value of the aggregate size in the dispersion. The logarithm of the normalized correlation function was fitted to a polynomial:

$$ \log (G(t)/B-1) = a + bt + ct^2 + \dots $$

where G(t) are the experimental correlation points, B is the baseline, and a, b, and c are the coefficients of the cumulant fit determined by a simple linear least-squares fitting procedure (a: intercept, b: slope giving access to the relaxation time of the signal). The diffusion coefficient D of the colloids was calculated with the following formula:

$$ 1/b = 2DK2 $$

where K = (4πn₀/λ)·sin(θ/2) is the scattering vector, n₀ the refractive index of the solvent, λ the laser wavelength and θ the scattering angle. The average hydrodynamic diameter of the particles Dh was then derived from D using the Stokes-Einstein equation and assuming spherical shapes:

$$ 1/b = kT/(3\pi\cdot\eta\cdot Dh) $$

where k is the Boltzmann constant, T the temperature and η the solvent viscosity.

Finally, the polydispersity index PI (0 < PI < 1), defined as c/2b, was estimated as a measure of the variance (squared standard deviation) of the particle size distribution.

### 3.8. Methods of Calculation for UV-Visible Properties of Single Compounds

Geometry optimization were achieved using density functional theory (DFT) calculations with the hybrid functional B3P86, which has appeared particularly adapted to assess both the thermodynamic and UV-visible absorption properties of polyphenols [20,25]. The Pople-type triple-ζ basis set, namely 6-31+G(d,p), allowed to establish reliable structure-property (UV-visible absorption) relationships for large series of polyphenols at a reasonable computational time [26]. The frequency analysis, performed at the same level of theory, confirmed the absence of imaginary frequencies in the ground state (GS). Excited singlet state (ES) energies were evaluated from the GS geometries, using the time dependent (TD) DFT formalism with the same functional and basis set. The allowed vertical π → π* electronic excitation energies were thus obtained, which consequently gave the absorption energies in the UV-visible range with the contribution of all one-electron transitions as well as their oscillator strength.
The TD-DFT calculations were performed both in the gas phase and in solvent, using the Polarizable Continuum Model (PCM) scheme. The integral equation formalism PCM (IEFPCM) method coupled to UA0 radii was used. In the PCM formalism, the solute is embedded in a shape-adapted cavity surrounded by the solvent implicitly described as a dielectric continuum characterized by its dielectric constant $\epsilon$ ($\epsilon = 32.6$ for MeOH). All calculations were performed with the Gaussian software package [27].

### 3.9. Methods of Calculation for Self-Association and Copigmentation

Self-association and copigmentation complexes were obtained from pure quantum chemistry calculations, using dispersion corrected DFT (DFT-D, D for "dispersion corrected"), which is a successful approach to evaluate noncovalent (mainly $\pi-\pi$ stacking) interactions in copigmentation complexes. However, pure quantum chemistry methods do not allow the explicit description of solvent, especially water, for such quite large molecular systems. As a consequence, the contribution of hydrogen bonding to stabilizing energies is overestimated. Moreover, only the stabilizing electronic energies (enthalpy changes), not the free enthalpy changes, were evaluated as calculation of entropic contributions in self-association and copigmentation is still an unsolved theoretical challenge [3]. In other words, this theoretical evaluation only provides trends and not absolute $\Delta G$ of association.

All geometries of the dimerization and copigmentation assemblies were obtained with the COSMO B3P86-D2/def2-SVP formalism in the ORCA package [3]. Electronic energies were calculated with the same functional using the larger basis set def2-TZVPP. For all noncovalent complexes, binding energies were obtained as the difference in energy between the complexes and the separated partners. The geometries were used to calculate theoretical UV-visible absorption properties (spectral shifts and electronic transitions) using time-dependent DFT with IEFPCM solvent model and $\omega B97XD/6-31+G(d,p)$ formalism (Gaussian software package). The $\omega B97XD$ functional appears crucial to correctly account for possible charge transfer in the excited state.

## 4. Conclusions

In conclusion, catechyl-pyranomalvidin-3-O-glucoside (PA1) and guaiacyl-pyranomalvidin-3-O-glucoside (PA2) are stable wine pigments that do not undergo water addition and thus simply behave as monoacids in mildly acidic solution. However, in the pH range 3.5-6 that is most important for application as food colorants, their intense color tends to evolve due to self-association of the neutral quinonoid base and eventual precipitation. If controlled by appropriate matrices, possibly including some copigments to avoid precipitation, this phenomenon, reported in this work for the first time, is very interesting as natural colors from red (monomer) to blue (aggregates) can thus be produced.

**Supplementary Materials:** Supplementary materials can be found at www.mdpi.com/1422-0067/17/11/1842/s1.

**Acknowledgments:** Anna Vallverdú-Queralt is grateful to the Alfonso Martín Escudero Foundation for the postdoctoral fellowship for carrying out research abroad. Funding for this work was provided by the Instituto de Salud Carlos III, ISCIII (CIBEROBN), the Quality Group from Generalitat de Catalunya (GC) 2014 SGR 773 and by INRA Departement CEPIA. Financial support from GIS IBiSA (Infrastructures en Biologie Santé et Agronomie), Région Languedoc Roussillon, and INRA CNOC for funding of the UPLC-MS and NMR equipment is also acknowledged.

**Author Contributions:** Patrick Trouillas, Véronique Cheynier and Olivier Dangles conceived and designed the experiments; Anna Vallverdú-Queralt and Michal Biler performed the experiments; Anna Vallverdú-Queralt and Michal Biler analyzed the data; Emmanuelle Meudec, Aude Vernhet, Christine Le Guernevé, Jean-Luc Legras, Jean-Paul Mazauric and Michèle Loonis contributed reagents/materials/analysis tools; and Anna Vallverdú-Queralt, Michal Biler, Patrick Trouillas, Véronique Cheynier and Olivier Dangles wrote the paper.

**Conflicts of Interest:** The authors declare no conflict of interest.

## References

1. Bayer, E.; Egeter, H.; Fink, A.; Nether, K.; Wegmann, K. Complex Formation and Flower Colors. *Angew. Chem. Int. Ed.* **1966**, *5*, 791–798. [CrossRef]
2. Brouillard, R. The in vivo expression of anthocyanin colour in plants. *Phytochemistry* **1983**, *22*, 1311–1323.
[CrossRef]

3. Trouillas, P.; Sancho-García, J.C.; de Freitas, V.; Gierschner, J.; Otyepka, M.; Dangles, O. Stabilizing and Modulating Color by Copigmentation: Insights from Theory and Experiment. *Chem. Rev.* **2016**, *116*, 4937–4982. [CrossRef] [PubMed]

4. Goto, T.; Kondo, T. Structure and Molecular Stacking of Anthocyanins. Flower Color Variation. *Angew. Chem. Int. Ed.* **1991**, *30*, 17–33. [CrossRef]

5. Yoshida, K. Metal Complex Pigment Involved in the Blue Sepal Color Development of Hydrangea. *J. Agric. Food Chem.* **2015**, *63*, 7630–7635.

6. Fulcrand, H.; Cameira dos Santos, P.J.; Sarni-Manchado, P.; Cheynier, V.; Favre-Bonvin, J.J. Structure of new anthocyanin-derived wine pigments. *Chem. Soc. Perkin Trans.* **1996**, *1*, 735–739. [CrossRef]

7. Schwarz, M.; Wabnitz, T.C.; Winterhalter, P. Pathway Leading to the Formation of Anthocyanin—Vinylphenol Adducts and Related Pigments in Red Wines. *J. Agric. Food Chem.* **2003**, *51*, 3682–3687. [CrossRef] [PubMed]

8. Vallverdú-Queralt, A.; Meudec, E.; Ferreira-Lima, N.; Sommerer, N.; Dangles, O.; Cheynier, V.; Guernevé, C.L. A comprehensive investigation of guaiacyl-pyranoanthocyanin synthesis by one-/two-dimensional NMR and UPLC-DAD-ESI-MSn. *Food Chem.* **2016**, *199*, 902–910. [CrossRef] [PubMed]

9. Quijada-Morín, N.; Dangles, O.; Rivas-Gonzalo, J.C.; Escribano-Bailón, M.T. Physico-Chemical and Chromatic Characterization of Malvidin 3-Glucoside-vinylcatechol and Malvidin 3-Glucoside-vinylguaiacol Wine Pigments. *J. Agric. Food Chem.* **2010**, *58*, 9744–9752. [CrossRef] [PubMed]

10. Oliveira, J.; Fernandes, V.; Miranda, C.; Santos-Buelga, C.; Silva, A.; Freitas, V.; Mateus, N. Color properties of four cyanidin-pyruvic acid adducts. *J. Agric. Food Chem.* **2006**, *54*, 6894–6903. [CrossRef] [PubMed]

11. Cruz, L.; Petrov, V.; Teixeira, N.; Mateus, N.; Pina, F.; de Freitas, V. Establishment of the Chemical Equilibria of Different Types of Pyranoanthocyanins in Aqueous Solutions: Evidence for the Formation of Aggregation in Pyranomalvidin-3-O-coumaroylglucoside-(+)-catechin. *J. Phys. Chem. B* **2010**, *114*, 13232–13240. [CrossRef] [PubMed]

12. He, J.; Oliveira, J.; Silva, A.M.S.; Mateus, N.; de Freitas, V. Oxovitisins: A new class of neutral pyranone-anthocyanin derivatives in red wines. *J. Agric. Food Chem.* **2010**, *58*, 8814–8819. [CrossRef] [PubMed]

13. Gómez-Alonso, S.; Blanco-Vega, D.; Gómez, V.; Hermosín-Gutiérrez, I. Synthesis, Isolation, Structure Elucidation, and Color Properties of 10-Acetyl-pyranoanthocyanins. *J. Agric. Food Chem.* **2012**, *60*, 12210–12223. [CrossRef] [PubMed]

14. Bakker, J.; Timberlake, C.F. Isolation, Identification, and Characterization of New Color-Stable Anthocyanins Occurring in Some Red Wines. *J. Agric. Food Chem.* **1997**, *45*, 35–43. [CrossRef]

15. Asenstorfer, R.E.; Jones, G.P. Charge equilibria and pK values of 5-carboxypyranomalvidin-3-glucoside (vitisin A) by electrophoresis and absorption spectroscopy. *Tetrahedron* **2007**, *63*, 4788–4792. [CrossRef]

16. Hakansson, A.E.; Pardon, K.; Hayasaka, Y.; de Sa, M.; Herderich, M. Structures and colour properties of new red wine pigments. *Tetrahedron Lett.* **2003**, *44*, 4887–4891. [CrossRef]

17. Pina, F. Chemical Applications of Anthocyanins and Related Compounds. A Source of Bioinspiration. *J. Agric. Food Chem.* **2014**, *62*, 6885–6897. [CrossRef] [PubMed]

18. Pina, F.; Melo, M.J.; Laia, C.A.T.; Parola, A.J.; Lima, J.C. Chemistry and Applications of Flavylium Compounds: A Handful of Colours. *Chem. Soc. Rev.* **2012**, *41*, 869–908. [CrossRef] [PubMed]

19. Oliveira, J.; Mateus, N.; de Freitas, V. Previous and recent advances in pyranoanthocyanins equilibria in aqueous solution. *Dyes Pigm.* **2014**, *100*, 190–200. [CrossRef]

20. Trouillas, P.; di Meo, F.; Gierschner, J.; Linares, M.; Sancho-García, J.C.; Otyepka, M. Optical properties of wine pigments: Theoretical guidelines with new methodological perspectives. *Tetrahedron* **2015**, *71*, 3079–3088. [CrossRef]

21. Hoshino, T. Self-Association of Flavylium Cations of Anthocyanidin 3,5-Diglucosides Studied by Circular Dichroism and ¹H NMR. *Phytochemistry* **1992**, *31*, 647–653. [CrossRef]

22. Hoshino, T. An Approximate Estimate of Self-Association Constants and the Self-Stacking Conformation of Malvin Quinonoidal Bases Studied by ¹H NMR. *Phytochemistry* **1991**, *30*, 2049–2055. [CrossRef]

23. Leydet, Y.; Gavara, R.; Petrov, V.; Diniz, A.M.; Parola, A.J.; Lima, J.C.; Pina, F. The Effect of Self-Aggregation on the Determination of the Kinetic and Thermodynamic Constants of the Network of Chemical Reactions in 3-Glucoside Anthocyanins. *Phytochemistry* **2012**, *83*, 125–135. [CrossRef] [PubMed]
24. Steyer, D.; Erny, C.; Claudel, P.; Riveill, G.; Karst, F.; Legras, J.-L. Genetic analysis of geraniol metabolism during fermentation. *Food. Microbiol.* **2013**, *33*, 228–234. [CrossRef] [PubMed]

25. Di Meo, F.; Sancho-García, J.C.; Dangles, O.; Trouillas, P. Highlights on anthocyanin pigmentation and copigmentation: A matter of flavonoid π-stacking complexation to be described by DFT-D. *J. Chem. Theory Comput.* **2012**, *8*, 2034–2043. [CrossRef] [PubMed]

26. Anouar, E.H.; Gierschner, J.; Duroux, J.L.; Trouillas, P. UV/Visible spectra of natural polyphenols: A time-dependent density functional theory study. *Food Chem.* **2012**, *131*, 79–89. [CrossRef]

27. 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.; et al. *Gaussian 09 (Revision A.02)*; Gaussian, Inc.: Wallingford, CT, USA, 2009.

<figure><img src="image_21.png" /></figure>

© 2016 by the authors; licensee MDPI, Basel, Switzerland. This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution (CC-BY) license (http://creativecommons.org/licenses/by/4.0/).