Next Article in Journal
Insights into the Lysine Acetylome of the Haloarchaeon Haloferax volcanii during Oxidative Stress by Quantitative SILAC-Based Proteomics
Next Article in Special Issue
Full-Length Transcriptome Sequencing Reveals the Molecular Mechanism of Metasequoia glyptostroboides Seed Responding to Aging
Previous Article in Journal
Improved Photocatalytic and Antioxidant Activity of Olive Fruit Extract-Mediated ZnO Nanoparticles
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Unveiling the Intimate Mechanism of the Crocin Antioxidant Properties by Radiolytic Analysis and Molecular Simulations

Institut de Chimie Physique, Université Paris-Saclay, UMR8000 CNRS, Rue Michel Magat, F-91405 Orsay, France
*
Authors to whom correspondence should be addressed.
Antioxidants 2023, 12(6), 1202; https://doi.org/10.3390/antiox12061202
Submission received: 26 April 2023 / Revised: 25 May 2023 / Accepted: 30 May 2023 / Published: 1 June 2023
(This article belongs to the Special Issue Antioxidant Mechanisms and Redox Signalling in Seeds)

Abstract

:
The successive steps of the oxidation mechanism of crocin, a major compound of saffron, by the free OH radical are investigated by pulse radiolysis, steady-state (gamma) radiolysis methods, and molecular simulations. The optical absorption properties of the transient species and their reaction rate constants are determined. The absorption spectrum of the oxidized radical of crocin resulting from the H-abstraction presents a maximum of 678 nm and a band of 441 nm, almost as intense as that of crocin. The spectrum of the covalent dimer of this radical contains an intense band at 441 nm and a weaker band at 330 nm. The final oxidized crocin, issued from radical disproportionation, absorbs weaker with a maximum of 330 nm. The molecular simulation results suggest that the OH radical is electrostatically attracted by the terminal sugar and is scavenged predominantly by the neighbor methyl site of the polyene chain as in a sugar-driven mechanism. Based on detailed experimental and theoretical investigations, the antioxidant properties of crocin are highlighted.

1. Introduction

Saffron is a spice derived from the flower of the plant Crocus sativus, commonly known as the “saffron crocus”, which can be classified as a potent antioxidant plant [1]. It is widely cultivated in countries around the world such as Iran, Greece, Italy, Spain, India, China, and Japan [2]. After collection and dehydration, it is used mainly for medicinal purposes [3,4], for food coloring, and as a flavoring agent.
Saffron contains an impressive variety of plant compounds. The three main components of saffron are the crocin, the major one (30%), which accounts for the yellow pigmentation from the stigmas [5], the picrocrocin (5–15%), which gives the rusty, bittersweet flavor, and the safranal (up to 2.5%), the volatile oil, which lends the earthy fragrance to the spice [2]. The percentages of these compounds determine the quality and commercial grade of the saffron [6]. Moreover, saffron contains small amounts of vitamins, gums, proteins, amino acids, sugars, mineral matter, flavonoids, and other chemical compounds [7,8,9,10,11]. The chemical structure of crocin, that is, the diester formed from two saccharides gentiobioses and the dicarboxylic acid crocetin, is constituted of a conjugated chain of a polyene and of two terminal sugars, and it is presented in Figure 1.
The compounds of saffron are known to have strong antioxidant and radical scavenger properties against a variety of radical oxygen species (ROS) and pro-inflammatory cytokines. They give it medicinal properties and are used for a long time in traditional medicine for the treatment of different types of diseases [1].
Many studies mentioned that the health-promoting properties of saffron are attributed primarily to crocin, a unique carotenoid with a powerful antioxidant capacity [2]. Crocin scavenges free radicals, mainly OH radicals and superoxide anions, and so may defend cells against oxidative stress [1]. The crocin bleaching assay was also designed according to this important property of crocin as a basic element for the antioxidant activity of saffron [12,13]. The high molar absorptivity of crocin [14,15] allowed us to study crocin concentrations as low as 10 μM.
The radiolytic systems ensure a selective and specific generation of the oxygen radicals OH (and/or O2−• under aerated conditions), as well as of hydrated electrons eaq. The crocin bleaching efficiency can, therefore, be correlated with the yield of these radicals. The pulse-radiolytic data demonstrated a very rapid diffusion-controlled attack by both hydroxyl radicals (OH) or hydrated electrons (eaq), whereas superoxide anions (O2−•) did not react at all [16]. Rate constants of OH and eaq with crocin were obtained using either a competition approach [17] (k(OH+croc) = 1.7 × 1010 M−1 s−1, k(e−aq+croc) = 13 × 1010 M−1 s−1) [14] or a kinetic evaluation of the rapid bleaching by pulse radiolysis (k(OH+croc) = 3.33 × 1010 M−1 s−1 and k(e−aq+croc) = 10.5 × 1010 M−1 s−1, respectively) [14].
One quantum chemical calculation only was published on the isolated crocin molecule [18], and several classical force field simulations about interactions of crocin with biological systems were performed [19,20]. No computational investigations of the absorption spectra of crocin, or of its radicals, have been published. However, several computational studies have been devoted to the structure and antioxidant properties of various other carotenoids [21,22,23,24,25,26], and to the conformation studies of fatty acids radicals [27].
Understanding the chemistry of mechanisms, advantages, and limitations of the methods is critical for the valid assessment of antioxidant activity in specific samples or conditions [28].
The aim of the present research is to study the reactions arising from the oxidation by OH radicals in the water of the crocin as the major compound of saffron by giving the detailed mechanism of OH attack and the nature of the products. The reaction rate constants of formation and decay of the radicals and their full absorption spectrum are measured by using the fast pulse radiolysis method and by following the transients. Based on the knowledge of the radiolytic OH yield in water, the spectra may be calibrated in molar absorption coefficients. In addition, as the crocin contains several sites for OH radical attack, simulations of the absorption bands of the radicals corresponding to these sites are performed in order to assign the structure of the radical formed from the comparison between the simulated and experimental absorption bands. The important role of the environment of the sites for the H-abstraction by OH radical is demonstrated. Finally, further products of crocin oxidation are also identified by giving the detailed mechanism of their formation.

2. Materials and Methods

2.1. Solutions Preparation

Pure crocin was obtained from Sigma Aldrich Chemistry (Saint Quentin Fallavier, France). The nitrous oxide (N2O) was obtained from Air Liquide Industrial Gases Company (Limay, France). N2O is added to scavenge hydrated electrons and H radicals, enhancing the production of oxidizing OH radicals. Because of the well-known photosensitivity of the crocin, the solutions were freshly prepared in the dark before use, and the flasks were wrapped with aluminum films. The optical absorbance of samples was measured using a spectrophotometer Hewlett-Packard (Palo Alto, CA, USA). Given the difficulty of accessing very pure crocin, we checked the optical absorption properties of our samples.

2.2. Radiolysis Experiments

Radiolysis (in the steady-state or the pulse regime) is a method delivering homogeneously known amounts of radicals. The oxidation study was performed using the panoramic 60Co source (2 kGy h−1) of ICP for steady-state irradiation, and the picosecond laser-induced electron accelerator ELYSE (Institut de Chimie Physique, Orsay, France) for pulse radiolysis [29,30]. The electron energy was 8 MeV, and the pulse length was 7 ps. The optical path of the cuvette was 0.5 cm, and the room temperature was 22 °C. The time-resolved optical absorption spectra were detected by using a Hamamatsu streak camera (Shizuoka, Japan). The UV-visible absorbances at 270–500 nm and the visible absorbances at 450–740 nm were detected by two different measurements. The uncertainty on the absorbance values was ±5%. The solution was refreshed between the pulses by using a circulating pump. The dose per pulse was derived from the absorbance of the hydrated electron in pure water with G(eaq) = 3.3 × 10−7 mole J−1 at 5 ns and ε = 18,000 M−1 cm−1 at 660 nm [31,32].
The production of the oxidizing radicals OH arises from the radiolysis of the most abundant water molecules:
H2O vvvv > eaq (2.8), H+ (2.8), H (0.62), H2 (0.47), OH (2.8), H2O2 (0.73)
(in brackets are the yields in 10−7 mole J−1 units) [15].
These yields correspond to a scavenging factor of 10−7 s−1 (or 100 ns after the radical production and their mutual reactions in spurs). In N2O-saturated solutions, the crocin scavenges by Reaction (4) the OH radicals arising either primarily from the radiolysis of water (Reaction (1)) or from the scavenging by N2O of primary hydrated electrons and H radicals (Reactions (2) and (3)), respectively):
N2O + eaq + H+  →  N2 + OH
N2O + H  →   N2 + OH
croc + OH  →  croc (–H) + H2O
Reaction (3) plays a minor role because it is very slow and the yield of radical H is small.

2.3. Molecular Simulations

UV-visible absorption spectra are often calculated through geometry optimization and harmonic frequency calculations of both the ground and excited states. This is the method developed by Brémond et al. who were able to address in this way large organic molecules, with an explicit introduction of vibronic couplings [33,34].
Over the last few years, we have developed a different method based on molecular simulations of the species of interest: (i) the potential surface of the solute is investigated with the help of the classical MC (Monte Carlo) method [35]; (ii) at each MC step the electronic structure of the solute is calculated with the DFT (density functional theory) method and the lc-ωPBE functional [36]; (iii) the solvent is modeled with the help of the PCM (polarized continuum model) method, in its SMD (Solvent Model using Density) variant [37], (iv) the absorption spectrum is given by the convolution of a list of TDDFT (time dependent DFT) lines calculated on a sublist of MC configurations using the formula from the literature [38]. This sublist is built with uniform sampling, selecting every 100th step of the MC simulation. This method enabled us to obtain useful insights into the mechanisms of radiolytic processes in solution [39,40,41].
In the present simulations, we use two MC moves: (i) multi-stretch moves, operating collective displacements of all the atoms [32], with probability of 0.8, and (ii) torsional moves about all the chemical bonds with probability of 0.2. For critical species, we have performed simulations with 60,000 MC steps and discarded the first 20,000 steps for the calculation of the spectrum. For other species, we simulated with 20,000 MC steps only. The convolution of TDDFT lines uses a Gaussian function with a given fwhm parameter (full width at half maximum), to be discussed.
We also considered spin contamination in DFT results [36]. This phenomenon is due to the one-determinant DFT wave function [42]. In the present work, we generalize the spin screening procedure we already used for moderating this drawback of DFT [36]: in the list of TDDFT lines contributing to the absorption spectrum, all the lines with a <S2> value larger than the value in the ground state by some amount are discarded. This amount, the so-called spin screening parameter, will be discussed.

2.4. Thermochemistry

The thermochemical quantities, such as association free energies, barrier heights, and oxidation processes, have been performed with the Gaussian package [34], with the lc-ωPBE functional [33], and also ordinary DFT functionals such as B3LYP, and cam-B3LYP, and with the SMD method [34]. For the simulations, we used the small SDD (Stuttgart Dresden) Gaussian basis set. This small SDD basis has been provided with d polarization functions on C atoms (exponent 0.83) and O atoms (exp. 1.), and the resulting basis set will be called pSDD (polarized SDD) hereafter. For a few calculations, we also provided the pSDD basis set with diffuse functions on C atoms: s (0.07), p (0.08) and d (0.55, 0.15), and on H atoms: s (0.3, 0.1), and the resulting basis will be called pSDD+ hereafter. Thermochemical calculations were done with the large aug-cc-pvdz basis set [34]. This large basis will be called apvdz hereafter.
Association free energies are difficult quantities, and we propose approximate values with the help of the Wertz formula [43]:
ΔrGpcmcorr = ΔrGpcm − 0.5 (ΔrGgp − ΔrHgp) + RT
where ΔrGpcm is the crude PCM value of the association free energy and ΔrGgp and ΔrHgp are the gas phase values. The +RT term cancels the PV contribution to enthalpy in the Gaussian code. This method is efficient for covalent dimers [44].
For van der Waals dimers, we also included in the ΔrGpcm term of Formula (5) dispersion interactions, and correction of the basis set superposition error (BSSE) [34]. Structures were visualized and processed with the GaussView code [45].
Crocin displays two bulky sugar substituents (Figure 1), making DFT simulations unworkable. Alternative solutions could be the use of QMMM or DFTB methods [34], which would accelerate the quantum calculations, but make the simulations much too long. We, therefore, modeled crocin by the di-methanol ester of crocetin (Me2–crocetin), by replacing the two gentiobiose groups of crocin with two methyl groups (Figure 1). The specific behavior of these crocin sugar groups was considered separately. This approach will be firmly justified later.
We finally optimized the ω parameter of the lc-ωPBE functional, so as to reproduce for Me2–crocetin at 0 K the measured value of λmax = 441 nm for crocin. Optimization yields the value ω = 0.244 and λmax = 441 nm, f = 3.22.

3. Radiolytic Oxidation Results and Discussion

3.1. Solubility of the Crocin in Water

We observed the optical absorption spectra of the samples at various concentrations to be compared with the literature data. The absorptivity of crocin is very intense, and the optical path of the cell was 1 or 2 mm. The crocin spectrum is constituted of an intense band with a maximum of 441 nm, a shoulder at 460 nm, and two weak UV bands at 250 and 333 nm (Figure S1) [46]. No absorbance is detected beyond 550 nm. Two studies evaluated the molar absorption coefficient of crocin at the maximum of 441 nm: ε441(croc) = 1.369 × 105 M−1 cm−1 [14], and ε441(croc) = 1.335 × 105 M−1 cm−1 [13]. The average value of ε441(croc) = 1.35 × 105 M−1 cm−1 will be used in the present work.
The absorbance increases proportionally to the concentration up to [croc] = 0.08 mM. But, noteworthy, beyond 0.13 mM, the absorbance is almost constant (Figure S1, inset). The crocin solubility in water is, therefore, limited to this value.

3.2. Reaction of Crocin with OH Radical within 1 μs

Kinetics of the reaction. Figure 2a presents the pulse radiolysis results within 1 μs of the difference spectrum arising from the OH scavenging by the crocin at the concentration of 0.0385 mM. At 1 μs, all OH radicals formed in Reactions (1) and (2) have reacted (Reaction (3) is too slow to contribute). A bleaching band is increasing with time at 441 nm at the same wavelength as that of the crocin and with the same width, as already observed [16]. In addition, a new absorption band is observed with a maximum of 678 nm, which increases with time up to 1 μs. Noteworthy, two isosbestic points are observed at 375 and 510 nm. They indicate first that the reaction product arises directly and stoichiometrically from the crocin oxidation by the OH radical and that it must be assigned to the oxidized radical R = (croc (–H)). Secondly, this radical absorbs at these isosbestic points with the same absorptivity as the crocin, but less in the interval 375–510 nm where bleaching is observed (Figure 2). Around 333 nm, where the crocin spectrum contains a weak band, the absorbance is slightly enhanced in the differential radical spectrum. These optical features of the radical at 678, 441, and 333 nm were not yet mentioned in the literature. That at 678 nm will be crucial to conclude from molecular simulations (see Section 4.1) on the site of the H-abstraction by OH radicals.
Figure 2b presents the kinetics in the range of 1 μs of the negative differential absorbance (bleaching) at 441 nm where it increases, of the positive absorbances at 678 nm and 333 nm where it increases, and at the isosbestic wavelengths 375 and 510 nm where it remains zero. In addition, a bump of a rapid increase soon after the pulse, followed by a decrease at 50 ns, is observed on each curve. It is assigned to the formation, within the pulse, of the hydrated electrons which then react with N2O and yield OH radicals in <10 μs.
The transient differential spectra within 500–750 nm at 1 μs and the kinetics at 678 nm are compared in Figure S2 for various crocin concentrations between 0.023 and 0.13 mM. Similar features are found for these concentrations.
The pseudo-first-order test applied to the absorbance A678 increase at 0.13 mM is a straight line (Figure S3). The second order rate constant of the Reaction (4) of OH with crocin is thus kcroc+OH• = 3.4 × 1010 M−1 s−1. This value is close to the literature one (kcroc+OH• = 3.33 × 1010 M−1 s−1) [16] and confirms that the crocin is an efficient antioxidant to scavenge the radicals OH. In fact, the particularly high kcroc+OH• value between neutral species favors crocin in the competition with other molecules, even more concentrated, if their scavenging rate constant is much lower than kcroc+OH•. The crocin efficiency as an antioxidant is still more significant if it would prevent the deleterious effect of induction by OH of oxidation chain reactions of neighbor molecules with oxygen.

3.3. Calibration in Molar Absorption Coefficient of the Radical Optical Spectrum

At the dose of 97 Gy per pulse with a yield of G(OH) = Geaq− + GOH = 5.6 × 10−7 mol J−1, the concentration of OH radicals is [OH] = 5.4 × 10−5 M per pulse. We conclude that after reaction with OH radicals at 0.06 mM (Reaction (4)), the crocin at the initial concentration of 0.0385 mM is totally oxidized into the radical (croc (H)). The differential spectrum at 1 μs in Figure 2 is thus the difference between the spectrum of the radical and that of the initial crocin, whose initial absorbance for 1 mm optical path was A441nm(croc) = 0.52. From the stoichiometry indicated by the isosbestic points and the ratio of the respective absorbances of crocin at 441 nm and radical at 678 nm with the equal optical path, we derive the molar absorption coefficient at 678 nm of the radical: ε678nm(croc (–H)) = ε441nm(croc) × A678nm(croc (–H))/A441nm(croc) = 6230 M−1 cm−1. (Table 1 and Figure 3b).
At 441 nm, the transient differential absorption spectra of the radical allows bleaching relative to the initial intense absorption band of the crocin (Figure 2a). The negative absorbance of the bleaching at 441 nm and 1 μs is almost three times higher than the positive absorbance at 678 nm (Figure 2). However, it is clearly much smaller than the expected bleaching arising from the consumed crocin corresponding to the OH scavenging. Therefore, the radical also absorbs at around 441 nm, but less than crocin. It appears that, in addition to the supplementary band around 678 nm, the radical spectrum is very similar to that of the crocin, though less intense in the interval between the isosbestic wavelengths 375 and 510 nm (Figure 3a). A previous study had also reported the increase in bleaching at 440 nm, much less intense than expected from the OH scavenging, and had concluded just a partial OH scavenging, of a few percent only [15]. The simultaneous increasing absorbance arising from the radical, also at the same wavelength range as the crocin, had not been considered, and the conclusion deduced about the OH scavenging efficiency was wrong. At 441 nm, the differential molar coefficient Δε441nm = ε441nm (croc) − ε441nm(croc (–H)) = 18,200 M−1 cm−1, that is, ε441nm(croc (–H)) = 1.17 × 105 M−1 cm−1. At 333 nm, Δε333nm = ε333nm(croc (–H)) − ε333nm (croc) = 3000 M−1 cm−1, that is, ε333nm(croc (–H)) = 0.10 × 105 M−1 cm−1 (Table 1). The complete radical spectrum is obtained from the algebraic sum of the crocin spectrum and of the differential spectrum of the radical, both calibrated in molar absorption coefficients (Figure 3).

3.4. Reaction of the Crocin Oxidized Radical (R) within 1 to 500 μs

Reaction kinetics. At ≥1 μs, all the OH radicals have been consumed (by scavenging or dimerization), and the formation of the radical is over. Then, beyond 1 μs, it is observed that the absorbance at 678 nm of the crocin oxidized radical decreases with time. Figure 4 presents the spectra evolution beyond 1 μs for [croc] = 0.10 mM and the kinetics of this decay over 84 μs. For better accuracy, the curve summarizes the absorbance data at 0.8 μs (Figure 4b), and of the decay kinetics over 10 and 84 μs (insets) (without the initial rapid increase in the radical formation).
The absorbance at 678 nm completely vanishes at 500 μs (Figure S4), meaning that the product arising from the decay reaction of two radicals does not absorb at this wavelength and that the absorbance must be assigned exclusively to the (croc (–H)) radical.
The second-order test of the radical decay at 678 nm is given in Figure 4b, inset. The test in 1/At is a straight line, indicating that two radicals react with a rate constant kobs = 0.6 × 106 s−1. According to the molar absorption coefficient derived above ε678 (croc (–H))) = 6.23 × 103 M−1 cm−1, the second order rate constant of the radical decay is 2kR+R = = 3.7 × 109 M−1 s−1. This reaction could be assigned to a dimerization or to a disproportionation of the radical. However, this diffusion-controlled process during 80 μs corresponds better to the radical dimerization (see the Simulation section):
2(croc (–H))  →  (croc (–H))2

3.5. Calibration of the Optical Absorption Spectrum of the Radical Dimer

During the first μs, the absorbance at 678 nm and the bleaching at 441 nm both increase with a constant ratio because they belong to the same radical spectrum (Figure 2b). In contrast, within 1 and 420 μs, the absorbance at 678 nm decays to zero, whereas the 441 nm bleaching still increases, but more slowly, up to a plateau where ΔA441(420 μs) = −0.74 (Figure 5b). The kinetics at 441 nm is a convolution of decay of the radical absorbance, similar to that at 678 nm, and of an increase in a more intense radical dimer bleaching, arising by Reaction (5). Another small increase is observed in the positive 330 nm band. Note that, during the supplementary increase, the shape of the bleaching band changes with a maximum wavelength shifting from 441 nm to 460 nm. No isosbestic point is observed.
The concentration of (croc (–H))2 is half of that of the radical precursor, and the molar absorption coefficient of the dimer is thus double that of one monomer. Therefore, the molar absorption coefficient at 441 nm of the radical dimer ε441(croc (–H))2 is deduced from the ratio between the differential absorbances at 1 and 420 μs by the following equation:
441(croc) − ½ ε441(croc (–H))2]/[ε441(croc) − ε441(croc (–H))] = ΔA441(420 μs)/ΔA441(1 μs)
that is, ε441(croc (–H))2 = 1.84 × 105 M−1 cm−1 (Table 1). The complete dimer spectrum is obtained from the algebraic sum of the crocin spectrum and of the differential spectrum due to the dimer with the bleaching band at 420 μs, both calibrated in molar absorption coefficients (Figure 5b). The spectrum of the radical dimer is more intense (though half of it is less intense) than those of the crocin, or of the radical, and the shoulder at 460 nm is less noticeable. At 333 nm, ε333(croc (–H))2 = 0.31 × 105 M−1 cm−1 (Table 1).

3.6. Gamma Irradiation of Crocin

A crocin solution at 0.057 mM was gamma-irradiated at an increasing dose of up to 211 Gy (Figure 6). The main optical absorption band at 441 nm and the small band at 230 nm of the crocin decrease in intensity while keeping the same shape. Simultaneously, a new band with a maximum of 330 nm increases. An isosbestic point is seen at 375 nm up to the dose of 97 Gy. Beyond 100 Gy, the isosbestic point disappears. The dose-dependence of the absorbance at these three wavelengths is presented in Figure 6a, inset.
Note that, under the conditions of the γ-irradiation, the end spectrum is observed in the range of several minutes, though the dose is of the same magnitude as in one electron pulse. Another main difference with pulse radiolysis is that the OH radicals are produced and scavenged steadily and that their concentration is much lower than in pulse radiolysis so that the OH scavenging is not in competition with their dimerization. Therefore, OH radicals are completely scavenged by crocin with the oxidation of crocin into the radical (croc (–H)), dimerizing into (croc (–H))2 (Reaction (6)) as a transient, then into the molecule (croc (–H2)) with 2 H-atoms abstracted, by a disproportionation reaction:
(croc (–H))2  →  (croc (–H2)) + croc
The radical (croc (–H)) resulting from the OH scavenging is also in low concentration. As well, it might first associate with the crocin into the complex (croc, croc (–H)) (Reaction (8)), before yielding the same oxidized molecule as above by direct disproportionation of this complex (Reaction (9)):
(croc (–H)) + croc  →  (croc, croc (–H))
2(croc, croc (–H))  →  3 croc + (croc (–2H))
The dose-dependent spectra are constituted of two components, the decreasing crocin (with the molar absorption coefficient ε441(croc) = 1.35 × 105 M−1 cm−1) and, after the disproportionation (9), the increasing reaction product assigned to the molecule (croc (–2H)). From the deconvolution of the dose-dependent spectra in γ-radiolysis up to 97 Gy (Figure 6b), the maximum of the (croc (–2H)) spectrum is at 330 nm. The molar absorption coefficients are ε330(croc (–2H)) = 0.25 × 105 M−1 cm−1 and ε441(croc (–2H)) = 0.11 × 105 M−1 cm−1. The isosbestic wavelength with the crocin is ε375(croc (–2H)) = 0.18 × 105 M−1 cm−1 (Table 1). The intensity of this spectrum is weak compared to those of crocin or of the radical dimer. This point will be discussed in Section 5 below. Beyond 100 Gy, crocin is partly exhausted, and the accumulated product (croc (–2H), also containing methyl sites, is now attacked by OH, therefore it plays the same role as an antioxidant. By this oxidation, a secondary product is formed with molar absorption coefficients at 330 and 375 nm lower than for croc (–2H) (Figure 6a, inset).
From the absorbance decrease at 97 Gy (ΔA441 = 0.49) and the differential molar absorption coefficient (Δε441 = ε441 (croc) − ε441 (croc (–2H)) = 1.24 × 105 M−1 cm−1), the experimental consumption yield of crocin in N2O-saturated solutions irradiated by γ-rays (Figure 6, inset) is Gexp(- croc)(γ) = 4.0 × 10−7 mole J−1 (Table 1). The hydrogen peroxide produced by the water radiolysis (Reaction (2)) is also able to oxidize slowly crocin during the γ-irradiation [14]. The calculated overall oxidation yield after the disproportionation would be equal to Gcalc(- croc)(γ) = Gcalc(croc (–2H))(γ) = ½ (Geaq + GOH + GH + 2GH2O2) = 3.84 × 10−7 mole J−1 (Scheme of Figure 7b and Table 1). This value is slightly lower than the experimental scheme above, possibly because of some simultaneous photodegradation of crocin. The agreement between Gcalc(- croc)(γ) and Gexp(- croc)(γ) and the existence of isosbestic points during the oxidation support the occurrence of the disproportionation Reaction (8). Instead, in pulse radiolysis, the dimerization is predominant and the oxidation by H2O2 is too slow and does not yet occur: Gcalc(croc(–2H))(pulse) = 1/2G(- croc) = 3.1 × 10−7 mole J−1 (Scheme of Figure 7 and Table 1).
The silver ion reduction is an adequate way to check the possible antioxidant properties of the radical (crocin (–H)). For that purpose, a mixed solution of 0.1 mM of crocin and 0.1 mM of silver perchlorate was irradiated by γ-rays. The silver cations scavenge the H radicals faster than N2O and are reduced into silver atoms (Reaction (9)).
Ag+ + H →  Ag0 + H+
It is observed (Figure 8a) that the absorption spectrum decreases at increasing doses with an isosbestic wavelength at 375 nm without Ag+. However, the bleaching yield of the crocin at 441 nm is much less in the mixed silver and crocin solution than in the solution without silver (Figure 6). Therefore, the radical (crocin (–H)) does not reduce Ag+ into a silver atom. The reduction potential of this radical E0(croc (–H2))/(croc (–H)) is thus much less negative than that of the radical of formate CO2−• [2], or of citrate (cit –(H)) [2], which did reduce the silver ions (E0(Ag2+/2Ag+) = −1.2 VNHE) and enhanced the antioxidant properties of their precursors.
Moreover, to explain that the observed yield with Ag+ Gexp(- croc)Ag+ = 2.1 × 10−7 mole J−1 is even less than the value without Ag+ (Gexp(- croc) = 4.0 × 10−7 mole J−1) (Table 1), we suggest that the silver atoms arising from the H scavenging do reduce the crocin into the radical (croc−•), which itself recombines with the oxidized radical (croc –(H)) into two crocin molecules (Scheme of Figure 8b):
Ag0 + croc    →    Ag+ + croc−•
croc−• + croc (–H) + H2O    →  2 croc + OH
The silver ions also catalyze the hydrogen peroxide decomposition. After complexation and disproportionation of the radicals (Reactions (8) and (9)), the calculated value Gcalc(- croc)Ag+ = Geaq + GOH• − GH• = 2.5 × 10−7 mole J−1 is in agreement with the experimental value (Table 1).

4. Thermochemistry of Crocin and Me2–Crocetin

4.1. Oxidation Free Energies

We considered all the possible radicals yielded by H abstraction from the model Me2–crocetin according to:
Me2–crocetin + OH  →   (Me2–crocetin (–H)) + H2O
The reaction free energies have been calculated at the lc-ωPBE (0.244)/apvdz/SMD level. Results for the H-abstraction from Ca and Cb carbons are given in Table 2 (the results for the symmetrical Ca′ and Cb′ carbons are obviously the same as for Ca and Cb). They have been obtained with the most stable conformer of each species. Results for abstraction from other species: –CH groups of the polyene chain and for the cation formation are given in Figure S5. It can be seen in Table 2 that oxidation is very efficient for the Ca and Cb carbon atoms, with reaction free energies close to −2 eV, and Ca is slightly more oxidable than Cb. Other carbon atoms, belonging to the polyene chain, are clearly less oxidable, with reaction free energies ranging between −1.6 eV (C3) and −1.2 eV (C4) (Figure S5).
In the following we note RCa and RCb the (Me2–crocetin (–H)) radicals after H-abstraction from the –CaH3 and –CbH3 groups.
Note that each crocin molecule contains four methyl groups and that their high affinity for OH radicals could explain the high value of the reaction rate constant per molecule measured above by pulse radiolysis. Moreover, according to the free energy (−2 eV) of the Me2–crocetin oxidation reaction to RCa or RCb (Table 2), and to the redox potential of OH radicals (E°(OH/H2O) ≈ 2.0 VNHE at pH 7), the redox potential of Me2–crocetin is E°(Me2–crocetin (–H)/Me2–crocetin) ≈ 0 VNHE. This means that as far as crocin may behave as Me2–crocetin, it is able also to play an indirect role of antioxidants in restoring, by H-atom transfer, oxidized neighbor molecules with positive redox potentials.

4.2. Transition States

Transition states (TS) have been calculated for the oxidation of β-carotene by the NO2 radical [20]. β-carotene has the same polyene chain as crocin, with four methyl groups, but very different end-of-chain groups. The authors found the following barriers to β-carotene oxidation: +0.45 eV for –CaH3 and +0.50 eV for –CbH3 (using the numbering of Figure 1). They also found lower barriers for the oxidation of the end groups. Calculations on a series of carotenoids show that oxidation barriers are always smaller on the end groups [22].
We found the calculation of some TS of Me2–-crocetin very delicate due to very small vibration frequencies. Nevertheless, we could obtain reliable values through the comparison of Me2–crocetin with two smaller systems: model 1, which has a much shorter polyene chain (Figure 9), and model 2, which has no COO groups (Figure S6). We also compared the small pSDD and large apvdz bases. The discussion is given in Section S4.2, results are shown in Tables S1 and S2 and Figure S7, and summarized in Figure 9.
TS free energies for the –CH groups are easily calculated. Figure 9 shows that the lower barriers for C3, C5, and C7 amount to ≈0.1 eV in model 1. In Table S2 we show the basis effect in model 1 and deduce that in Me2–crocetin these barriers are a little bit smaller, but remain positive. TS free energies for the –CH3 groups proved very delicate. In model 1, we had to force harmonic analysis with residual nuclear gradients ≈ 10−4 a.u., larger than the usual threshold and obtained negative values ≈ −0.4 eV. We could confirm this negative sign with the help of model 2 where both –COOCH3 groups have been replaced by two –CH3 groups (Figure S6 and Table S2). In model 2, the negative value −0.6 eV for –CaH3 was easily and securely obtained. Negative TS energies are not surprising: they mean that the free energy curve is attractive at longer distances and that the TS lies under the asymptote. Note that this negative value is due to the OH radical: Table S2 shows that the NO2 radical gives the positive free energy: +0.3 eV to this same TS in model 2, in agreement with the value of Ref. [20].
Structures of a few TS are shown in Figure S8. They explain the issues encountered by the –CH3 groups. The ”easy” TS of –CH groups is of a stretching type, with large imaginary frequencies, whereas the TS of the –CH3 groups are of a bending type, with OH pointing toward one O atom, provoking very small imaginary frequencies.
We conclude that the H abstraction from the methyl groups –CaH3 and –CbH3 is easily done, without barrier, and that this abstraction from the –CH groups of the polyene chain is scarcely possible at 22.5 °C. This is consistent with the reaction free energies of Table 2 and Figure S5. Besides, the H abstraction by the OH radical from a glucose molecule with the same method was also investigated. All these calculations yielded high potential barriers: +0.29 eV and +0.33 eV for the carbon and oxygen atoms of the –CH2OH group, and + 0.36 eV for –OH groups of the cycle. These numbers confirm that the sugar groups of crocin are much less oxidable than the methyl groups of the chain.

4.3. Conformation Analysis

It is well known from the literature that cis–trans isomerization can be critical for molecules with a polyene chain, such as carotenoids [24] or fatty acids [25]. In the following, we call all-trans the structure of Figure 1, with no cis–trans isomerization and (x,y) the cis–trans isomerization about the CxCy bond. We first investigated the cis–trans issue for the crocin and Me2–crocetin molecules, optimizing the seven structures corresponding to one cis–trans isomerization about one bond. In Table S3, we give the free energies and low energy optical absorption data for the all-trans and the first three cis–trans conformations lying under ≈0.1 eV. It can be seen that crocin and Me2–crocetin behave similarly as long as conformation and optical absorption are concerned.
Both molecules are mainly all-trans, with the lowest (6,7) cis–trans isomer 0.07 eV higher. Both other (2,3) and (8,8′) conformations also interfere. In Me2–crocetin a very constant 5–6 nm redshift, with respect to crocin, is observed in every conformation. This shift is very small. In both molecules, the optical data hardly depend on the conformation, with again a 5–6 nm shift between conformations.
For both radicals RCa and RCb, we used a more demanding method: (i) optimization of the 13 conformers involving one cis–trans isomerization about the 13 chemical bonds of each radical, (ii) 13 MC simulations at 1000 K, initialized with each of these optimized structures, and (iii) geometry optimizations and harmonic analysis from a sub-list of 130 MC conformations of these simulations. This procedure enables the test of multiple cis–trans isomerizations. The free energies and first excited states are given in Table S4.
The RCa radical also has a dominant all-trans conformation, but the bent conformers (5,6), (2,3), and (2′,3′) lie only 0.03, 0.5, and 0.08 eV higher. One intense transition at 460–470 nm can be seen, and two weak, lower energy transitions interfere. The RCb radical displays two quasi-degenerate conformers, the all-trans and the (6,7) bent ones. Two other bent conformers, (5,6) and (5,6) (6,7), lie only 0.025 and 0.05 eV higher. Therefore, a double cis–trans isomerization interferes. One intense transition at 410–420 nm can be seen, and three weak, lower energy transitions interfere.
The absorption behavior of the radicals is a critical issue. In Figure S9, the absorption lines of the optimized structures of Table S4 have been broadened in order to mimic absorption spectra. Figure S9 shows that cis–trans isomerization modifies the absorption spectra very slightly.

5. Molecular Simulations

The simulated optical absorption spectra may be directly compared with the observable time-resolved or final spectra of radiolysis. When fitting with the experimental ones, they give essential insights into the mechanisms of oxidation. It was mentioned in Section 2.3 that the extraction of absorption spectra from molecular simulations is conditioned by the fwhm parameter, that is, the broadness of the Gaussian function used for the convolution of TDDFT lines, and by the spin screening parameter, used for the empirical reduction of spin contamination in radical species. The choice of these parameters is discussed in Section S5 and Figure S10 in the case of the RCa radical. According to this discussion, three values of the fwhm parameter will be used: fwhm = 0.7 eV for the broadening of TDDFT lines of optimized structures, fwhm = 0.50 eV for the broadening of TDDFT lines of the intense bands, and for the weak, high-energy bands of the simulation, and fwhm = 0.20 eV for the weak, low energy bands of the simulation. The spin screening parameter will be 20%.

5.1. Absorption Spectrum of Me2–Crocetin

The absorption spectrum of crocin, obtained experimentally (Figure S1) is rather complex, with two weak, high-energy bands at 260 and 320 nm and an intense, dissymmetrical, and structured band at 441 nm. The simulated spectrum of Me2–crocetin, unconstrained and initialized with the all-trans conformer, is shown in Figure 10a, calculated with the pSDD basis, the fwhm value 0.50 eV, and 60,000 MC steps. This spectrum displays some of the features observed for crocin: the weak bands at 260 and 320 nm, and an intense band at 448 nm. However, this intense band is symmetrical and displays no structure. Figure 10b shows the same spectrum, calculated with the larger pSDD+ basis, but this only produces a slight redshift. This raises the question of whether the discrepancies between measurement and simulation are due to a defective simulation of Me2–crocetin, or to our modeling of crocin by Me2–crocetin.
The answer to this question is given by the conformation analysis: in fact, the isomerization free energies of Table S3: 0.07, 0.08, and 0.11 eV, are rather large, but the DFT calculations and harmonic analysis make them rather inaccurate. Actually, these free energies only suggest that the interference of bent conformers cannot be excluded. Therefore, we have performed constrained 20,000 steps MC simulations of all the bent conformers of Me2–crocetin, with only atom multi-stretch moves and no dihedral angle fluctuation. Typical structures of the most stable four conformers of Me2–crocetin according to Table S3, namely the all-trans conformer and the (2,3), (6,7), and (8,8′) bent conformers, are shown in Figure S11. Figure 11a presents their spectra. It can be seen that (i) the weak band at 250 nm can be attributed to the all-trans and (2,3) bent conformers, (ii) the weak band at 320 nm to the bent (6,7) and (8,8′) conformers, (iii) the all-trans, (2,3) bent and (6,7) bent conformers have very close maxima: 438, 444, and 440 nm, respectively, but the bent (8,8′) conformer at 451 nm undergoes a 13 nm red shift compared to the all-trans conformer. The conformers’ spectral features of Me2–crocetin (Figure 11a) provide a plausible interpretation of the shape of the experimental spectrum of crocin in Figure 11b. In particular, the striking shoulder on the red wing of the spectrum is assigned to the (8,8′) conformation.
Note that the lengthy simulated, unconstraint spectrum of Figure 10a, with λmax = 448 nm is close to the constraint, all-trans spectrum of Figure 11a with λmax = 438 nm. This unconstraint simulation has been initialized with the all-trans conformer and keeps it all the time. According to the large free energies of Table S3, this is quite normal. Therefore, the present unconstrained simulation, even if extended, cannot yield the structure of the observed spectrum. The right Monte Carlo treatment of these conformation issues would be the parallel tempering Monte Carlo method [32], with several parallel, coupled MC simulations at different temperatures. However, the analysis of Section 4.3 shows that parallel tempering is not necessary and that ordinary, unconstrained MC simulations of the main conformers provide a realistic investigation of the absorption spectra.
We conclude that Me2–crocetin is an adequate model of crocin. The present results for Me2–crocetin, including the interference of cis–trans conformers, emphasize the predominant influence of the polyene chain on the optical spectrum. Significantly, they may explain not only the spectrum of crocin but also the spectrum shape of β-carotene and of numerous carotenoids, mainly the same intense and dissymmetrical band [47].

5.2. Absorption Spectra of Me2–Crocetin Radicals

The spectra of all the possible oxidation products of Me2–crocetin were then simulated. Since the reaction barriers of Figure 9 show that the methyl groups are easily oxidable, the spectra of the RCa and RCb species are given in more detail. In Figure 10 the focus is made on the intense bands with fwhm = 0.50 eV and the small pSDD (a) and large pSDD+ (b) bases, and in Figure 12 we focus on the weak bands only, at higher (a) and lower (b) energies with the large basis only. We label “½ ½” the half-sum of the spectra of RCa and RCb species, assuming equal oxidation probabilities, and “5/6 1/6” another combination of the spectra assuming that RCa is predominant, these numbers will be explained. The spectra of the RC2, RC3, RC4, RC6, and RC7 species are given in Figure S12, together with that of the cation.
We now compare Figure 3a,b of the experimental spectrum and Figure 10 and Figure 12 of the simulated spectrum, respectively. Figure 3a shows that the oxidation product has an intense absorption band very similar to that of crocin, with a slight decrease in intensity. Figure 10a shows that the two possible products, RCa and RCb, have an intense band with very different maxima at 467 nm and 415 nm, respectively. From the large blue shift of 33 nm (0.20 eV) and the large absorption intensity decrease obtained for RCb alone, it appears that this radical is not predominant in contrast with RCa, yielding a slight decrease in intensity and a 20 nm redshift (0.10 eV). If RCa and RCb were equally formed (“½ ½” hypothesis), they would yield a larger decrease in intensity and a tiny 5 nm blue shift. Using the larger basis does not change this situation (Figure 10b).
Figure 12a (to be compared with Figure 3b) shows that the weak, high-energy band observed at 330 nm must be attributed to RCa only. The simulated band (in red) lies at 343 nm (error: 13 nm, 0.1 eV). Figure 12b shows that the observed weak, low-energy absorption band must be also attributed to RCa only. Note that the observed maximum lies at 678 nm and that the simulated values are 646 nm with the small basis (not shown, error: 0.1 eV) and 654 nm with the larger basis (error: 0.05 eV). Figure 12b also shows that the little bump at 620 nm, which is present in the measured spectrum, may be attributed to RCb alone. The simulated band lies at 592 nm (error: 0.10 eV). This bump appears to be the only direct interference of RCb in the spectrum. In brief, Figure 10 and Figure 12 suggest that RCa is predominant, and RCb observable but rare, in contradiction with thermochemistry which rather implies equal probabilities of the two radicals and the ½ ½ spectra of Figure 10 and Figure 12. The predominant formation of RCa now requests supplementary arguments.
We first note that the real crocin molecule offers to the OH radical a polyene chain and two bulky sugar substituents of comparable sizes. The probability that OH first meets the polyene chain amounts to ≈1/3 only. Therefore, and according to thermochemistry, the oxidation probabilities of RCa and RCb are the same: ½ × 1/3 = 1/6.
We then consider that, arriving on a sugar substituent, the OH radical might undergo an electrostatic attraction by the sugar. Actually, the partial charges of the oxygen atoms of the gentiobiose unit (between −0.9 and −1.1 a.u.) are higher than the partial charge of oxygen in the water molecules (−0.7 a.u.), and the charge of the H atom of the OH radical (+0.4 a.u.) is higher than that in the water molecule (+0.34) (Figure S13). However, note also that due to thermodynamical barriers and despite this electrostatic attraction, the OH radical cannot oxidize the sugar moieties and will rather migrate from site to site on their surface. This favors the encounter and oxidation of the neighbor –CaH3, rather than –CbH3 (Figure 13). The sugar-driven probabilities of RCa formation becomes: 2/3 + 1/6 = 5/6 and of RCb: 0 + 1/6 = 1/6. The corresponding spectra are labeled “5/6 1/6” in Figure 10 and Figure 12.
According to this discussion, the formation probabilities of the RCa and RCb products should lie between the extreme cases: (1/2 1/2) (no interference of the sugar substituents) and (5/6, 1/6) (sugar-driven oxidation of –CaH3). Obviously, this probability 5/6 has been roughly estimated and simply says that the RCa product is predominant. Figure 10 and Figure 12 suggest that this is the case. Note also that the cation, with its spectrum on Figure S12, is clearly not observed.

5.3. Absorption Spectra of Covalent Dimers

The RCa and RCb species may undergo covalent dimerization, through singlet pairing of their two unpaired electrons. In Figure 5b the dimer displays an intense band, with a slight redshift and an absorbance enhancement, with respect to crocin, and also a weaker band at 320 nm, with a ~50% absorbance enhancement with respect to the weak band of crocin at 250 nm. However, the corresponding dimers are too large for the present method of molecular simulation. We could perform simplified simulations, nevertheless, with dimers made of frozen monomers, taken in the optimized structures of the dimers, shown in Figure S14. This yields the spectra of Figure 14, where the spectrum of frozen Me2–crocetin, namely of its structure at 0 K, is added for comparison.
In Figure 14 we also extended to the dimers our two hypotheses about the RCa and RCb monomers: if the monomers have the probabilities (½ ½), then the dimers must have the probabilities (1/2)2 = 0.25 for CaCa, 2 × (1/2)2 = 0.5 for CaCb and 0.25 for CbCb; if the monomers have the probabilities (5/6 1/6), then the dimers must have the probabilities (5/6)2 = 0.69 for CaCa, 2 × (5/6) × (1/6) = 0.28 for CaCb and (1/6)2 = 0.03 for CbCb. Figure 14a shows that dimerization modifies the intense band of Me2–crocetin, with a little redshift and an intensity enhancement, in agreement with the measurement of Figure 5b. Figure 14b shows that the observed band at 330 nm is reproduced by the calculations. The ratio between the intensity of this band and that of the band at 250 nm of Me2–crocetin amounts to ≈2.6 for the (½ ½) combination of the monomers, to ≈1 for their (5/6 1/6) combination, to be compared to the 1.5 factor of Figure 5b. This suggests that the predominance of the RCa product, as concluded from the results at shorter times, is still observed in the spectra of their dimers, and that the 5/6 probability of this product is probably overestimated. Actually, the spectra of the dimers would deserve more rigorous simulations.

5.4. Absorption Spectra of Products of Radical’s Disproportionation

In γ-radiolysis (at a weak dose rate), the final products (Figure 6b) have a spectrum much different from that of radical dimers. One significant difference in γ radiolysis, compared to pulse radiolysis, is the low concentration of radicals and the formation of van der Waals complexes, (crocin, crocin (–H)) between radicals and non-oxidized crocin (Reactions (7) and (8)).
Investigation of such complexes is at the limit of our present methods because of the size of the systems involved. We first considered the complexes (Me2–crocetin, Me2–crocetin(–H)). The optimized structures of the complexes of the RCa and RCb radicals are shown in Figure S15. The complexes display a parallel arrangement of the polyene chains, thus achieving mutual solvation of these hydrophobic moieties. The formation free energies of these complexes are very negative (−0.5 eV) (Table 2). This means that if such a complex forms, then it does not separate. Unfortunately, these structures are disappointing because the encounter of two such complexes can provoke either disproportionation, if the H transfer occurs, or dimerization as well, if one radical center, CaH2 or CbH2, of one of the complexes, comes close to the radical center of the other complex. However, this possibility of dimerization should be reduced by bulky substituents. With this in mind, we introduced glucose moieties and optimized a few structures of the (Glu2-crocetin, Glu2-crocetin(–H)) complexes. At this stage, we do not claim to have a general view of such complexes. One typical structure involving the RCa radical is shown in Figure S15. Note that this structure is ruled by parallel polyene chains, like in the Me2–crocetin complex, and also by hydrogen bonds between neighbor glucose moieties.
If now two such complexes are encountered, it seems clear that a parallel arrangement is hindered by the substituents, and that crossed structures, with an X shape or a T shape, will be rather favored. Then, the dimerization vs. disproportionation alternative occurs if, in the resulting dimer of complexes, both radicals are in contact inside the structure. If these radicals are both of the predominant RCa type, then the crossed shape will forbid dimerization, and allow H atom transfer and disproportionation. If one of these radicals is of the rare RCb type, dimerization cannot be excluded but will be rare as well.
It can be seen in Table 2 that disproportionation involves two steps: first transfer of one H atom from one radical to the other one, yielding a primary, twice oxidized species, then cyclization of this species. In the following, we note Me2–crocetin(–2H) CxCy the product of disproportionation, where the Cx and Cy carbon atoms have been oxidized. It can be seen also that in most cases the primary species is a closed shell molecule with alternate single and double CC bonds between two =CH2 groups. This is the case for the Me2–crocetin(–2H) CaCa′, CaCb′, or CbCb′ species, see Figure S16. This is due to an odd number of bonds between the two oxidized carbon atoms (Figure 1). Table 2 shows that the formation free energies of these species amount to about −0.3 eV. In one case only, Me2–crocetin(–2H)CaCb, the primary species is not a closed shell molecule, because of an even number (6) of CC bonds. This species is a diradical, made of two –CH2 groups, for which the triplet state can be calculated by DFT, but the singlet state cannot. In Table 2 we consider that in this case, the triplet state is more stable than the singlet, thus suggesting that the formation free energy of this singlet species is positive, and the species inaccessible. In any case both –CH2 groups can come close to each other and fuse as a –CH2–CH2– bond, leaving a cyclic, closed shell molecule (Figure S16). Note that cyclization is also possible for the unexpected Me2–crocetin(–2H) CaCb system. It is clear that the primary twice oxidized radicals are transients and cannot be observed in γ radiolysis.
In Figure 15 the absorption spectra of the four possible cyclic products are shown: Me2–crocetin(–2H) CaCa′, CaCb, CaCb′, and CbCb′, and of three of their possible combinations. The spectra of the three accessible transients, primary products are also shown in Figure S17.
The mechanism of the crocin oxidation by γ-radiolysis may be derived from the best fitting between the experimental absorption spectrum of the product in Figure 6b and the simulated spectrum of Figure 15 and Figure S17. Figure 15a shows that, in contrast with Figure 6b, the cyclic CaCb (or the symmetric Ca′Cb′) molecule has an intense absorption spectrum similar to that of Me2–crocetin (with a 15 nm blue shift). The three other accessible cyclic products, CaCb′, CaCa′, and CbCb′, have rather close spectra very different from that of Me2–crocetin and with only weak bands in the 400–450 nm zone. If we consider the 1/3 1/3 1/6 1/6 combination, assuming that the RCa and RCb species have the same probability and that the CaCb cyclic is also present, this spectrum would have still an intense band at 425 nm, which is clearly not observed. The 0 1/2 1/4 1/4 combination makes the same assumption on the RCa and RCb species but discards the CaCb cycle. This spectrum, with a weak band at 321 nm and a still weaker band at 407 nm is in fair agreement with the observed spectrum. The 0 1/2 5/12 1/12 combination also discards the CaCb cycle and reflects the probabilities 5/6 and 1/6 for RCa and RCb. Actually, the corresponding spectrum, with a weak band at 323 nm and a still weaker band at 414 nm, is also very close to the observed spectrum because the spectra of the involved cyclic products are close to each other. Figure 6 shows that after γ radiolysis, the absorbance of the sample is drastically reduced, up to ≈1/9 of that of crocin at 441 nm. Figure 15b shows that the simulations well reproduce this fact, with the same factor of 1/9. Figure S17 presents the spectra of the transient primary species and their 0 1/2 1/4 1/4 combination. It can be seen that all the spectra are different from those of the final cyclic molecules and that their combined spectrum is clearly not observed.
In brief, this section provides the interpretation of the observed spectrum after γ radiolysis (weak dose rate). The Me2–crocetin(–2H) CaCb product is not observed. The three accessible cyclic products are observed through their superposition, with a weak band at ≈320 nm, to be compared with the observed one at ≈330 nm, and a weaker low energy band at ≈410 nm to be compared with the observed one at 430 nm. The drastic absorption coefficient collapse after γ-radiolysis is well reproduced, but the spectra cannot help discriminate the assumptions about the probabilities of the RCa and RCb species.
The critical point is that the RCa species is predominant and has a radical center close to a bulky substituent. This prevents the close approach of two reactive –CaH2 groups to each other, and thus from a covalent dimerization. Note that this is true only if (crocin, crocin (–H) radical) complexes predominantly exist, with high steric requirements. In pulse radiolysis, such complexes do not exist, the approach of the radicals is ruled by chance and the sugars may step aside and give way to dimerization.

6. Conclusions

To understand the detailed mechanism of the antioxidant properties of crocin, observations of transients to the end-products of the crocin oxidation were performed by the radiolysis method, in pulse or steady-state regime, and crossed with demanding molecular simulations.
The high value of the rate constant of the crocin oxidation by the OH radical confirms the efficient direct antioxidant properties of crocin, thus protecting the neighbor molecules from oxidation, and from the deleterious effects of possible induced chain reactions with oxygen. Moreover, the thermochemical calculation on the model Me2–crocetin indicates that by its redox potential, crocin may behave also as an indirect antioxidant in restoring, through an H-atom transfer, neighbor molecules from their oxidized radical. Note that the present radiolysis and simulation studies are restricted to the oxidative attack by one OH radical per molecule. Moreover, the stable products of the reaction still possess several similar oxidable sites that enhance the antioxidant power of this kind of molecule.
The simulated spectra, which reproduce fairly well the spectra obtained by radiolysis, gave essential insights into the molecular structures and the detailed mechanisms of oxidation. The observed optical absorption spectrum of crocin, mostly constituted of an intense asymmetrical band, is interpreted with the help of conformation analysis. This analysis emphasizes the predominant influence of the polyene chain conformers on the optical spectrum. Significantly, they may also explain the similar ones of numerous carotenoid spectra, mainly the same intense and asymmetrical band at around 440 nm.
The optical spectrum of the radical resulting from the crocin oxidation by OH radical is constituted of a specific band at 678 nm in addition to the main asymmetrical band, similar to crocin at 441 nm but slightly less intense. The thermochemical calculations and the molecular simulations of the optical spectra provide an assignment of the site of the OH radical attack on crocin. It is important to note that this OH attack, contrary to other carotenoids, occurs only on the polyene chain. Interestingly, the sugar sites of crocin do not scavenge but just attract electrostatically OH radicals, and thus play an important role in favoring the H-abstraction from methyl groups of the chain, predominantly from the close neighbor ones (sugar-driven mechanism).
After a pulse (of high dose rate), the radicals dimerize into a stable covalent dimer. During the gamma irradiation (of lower dose rate), the radicals first are complexed by excess crocin molecules before their disproportionation into a twice-oxidized crocin and a crocin molecule. From molecular simulations, the oxidized molecules may display various cyclic conformations.

Supplementary Materials

The following supporting information, containing 17 Figures and 4 Tables can be downloaded at: https://www.mdpi.com/article/10.3390/antiox12061202/s1.

Author Contributions

Conceptualization, M.M., P.A. and J.B.; methodology, S.A.G., D.A. and P.A.; software, P.A.; validation, M.M., J.B. and P.A.; formal analysis, S.A.G., M.M., J.B. and P.A.; investigation, S.A.G.; writing—original draft preparation, S.A.G., J.B., P.A. and M.M.; writing—review and editing, J.B., P.A. and M.M.; supervision, J.B., P.A. and M.M.; project administration, M.M.; funding acquisition, M.M. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

The data that support the findings of this study are available on request from the corresponding authors (M.M. and P.A.).

Acknowledgments

The authors thank Jean-Philippe Larbre and Pierre Jeunesse of the picosecond pulse radiolysis facility ELYSE, for their assistance in the time-resolved measurements and Dominik Domin for technical assistance in the computations.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Mashmoul, M.; Azlan, A.; Khaza’ai, H.; Yusof, B.N.M.; Noor, S.M. Saffron: A natural potent antioxidant as a promising anti-obesity drug. Antioxidants 2013, 2, 293–308. [Google Scholar] [CrossRef] [PubMed]
  2. Ramadan, A.; Soliman, G.; Mahmoud, S.S.; Nofal, S.M.; Abdel-Rahman, R.F. Evaluation of the safety and antioxidant activities of Crocus sativus and Propolis ethanolic extracts. J. Saudi Chem. Soc. 2012, 16, 13–21. [Google Scholar] [CrossRef]
  3. Kalalinia, F.; Ghasim, H.; Farzad, S.A.; Pishavar, E.; Ramezani, M.; Hashemi, M. Comparison of the effect of crocin and crocetin, two major compounds extracted from saffron, on osteogenic differentiation of mesenchymal stem cells. Life Sci. 2018, 208, 262–267. [Google Scholar] [CrossRef] [PubMed]
  4. Wang, M.Z.; Gao, J.; Chu, Y.; Niu, J.; Chen, M.; Shang, Q.; Peng, L.H.; Jiang, Z.H. Synthesis of crocetin derivatives and their potent inhibition in multiple tumor cells proliferation and inflammatory property of macrophage. BMC Compl. Med. Ther. 2020, 20, 29. [Google Scholar] [CrossRef]
  5. Kabiri, M.; Rezadoost, H.; Ghassempour, A. A comparative quality study of saffron constituents through HPLC and HPTLC methods followed by isolation of crocins and picrocrocin. LWT 2017, 84, 1–9. [Google Scholar] [CrossRef]
  6. Kyriakoudi, A.; Chrysanthou, A.; Mantzouridou, F.; Tsimidou, M.Z. Revisiting extraction of bioactive apocarotenoids from Crocus sativus L. dry stigmas (saffron). Anal. Chim. Acta. 2012, 755, 77–85. [Google Scholar] [CrossRef]
  7. Abdullaev, F.I. Biological effects of saffron. BioFactors 1993, 4, 83–86. [Google Scholar]
  8. Rios, J.L.; Recio, M.C.; Giner, R.M.; Manez, S. An update review of saffron and its active constituents. Phytother. Res. 1996, 10, 189–193. [Google Scholar] [CrossRef]
  9. Winterhalter, P.; Straubinger, M. Saffron—Renewed interest in an ancient spice. Food Rev. Int. 2000, 16, 39–59. [Google Scholar] [CrossRef]
  10. Cerdá-Bernad, D.; Valero-Cases, E.; Pastor, J.; Frutos, M. Saffron bioactives crocin, crocetin and safranal: Effect on oxidative stress and mechanisms of action. Crit. Rev.Food Sci. Nutr. 2020, 62, 3232–3249. [Google Scholar] [CrossRef]
  11. Chengtao, S.; Shivraj, H.N.; Yiting, Z.; Luping, Q.; El-Seedi, H.R.; Daglia, M.; Guoyin, K. Novel insight into utilization of flavonoid glycosides and biological properties of saffron (Crocus sativus L.) flower by-products. J. Agric. Food Chem. 2020, 68, 10685–10696. [Google Scholar]
  12. Schmidt, M.; Betti, G.; Hensel, A. Saffron in pytotherapy: Pharmacology and clinical uses. Wien. Med. Wochenschr. 2007, 157, 315–319. [Google Scholar] [CrossRef]
  13. Liu, N.; Yang, Y.; Mo, S.; Liao, J.; Jin, J. Calcium antagonistic effects of chinese crude drugs: Preliminary investigation and evaluation by 45ca. Appl. Radiat. Isot. 2005, 63, 151–155. [Google Scholar] [CrossRef]
  14. Friend, J.; Mayer, A.M. The enzymic destruction of carotenoids by isolated chloroplasts. Biochim. Biophys. Acta 1960, 4, 422–429. [Google Scholar] [CrossRef]
  15. Weber, F.; Laskawy, G.; Grosch, W. Co-oxidation of carotene and crocin by soyabean lipoxygenase isoenzymes. Z. Lebensmit.-Untersuch. Forsch. 1974, 155, 142–150. [Google Scholar] [CrossRef]
  16. Bors, W.; Saran, M.; Michel, C. Radical intermediates involved in the bleaching of the carotenoid crocin. Hydroxyl radicals, superoxide anions and hydrated electrons. Int. J. Radiat. Biol. Relat. Stud. Phys. Chem. Med. 1982, 41, 493–501. [Google Scholar] [CrossRef]
  17. Buxton, G.V. An overview of the radiation chemistry of liquids. In Radiation Chemistry from Basics to Applications in Material and Life Sciences; Spotheim- Maurizot, M., Mostafavi, M., Douki, T., Belloni, J., Eds.; EDP Sciences-L’Act. Chim.: Paris, France, 2008; pp. 3–16. [Google Scholar]
  18. Akhtari, K.; Hassanzadeh, K.; Fakhraei, B.; Fakhraei, N.; Hassanzadeh, H.; Zarei, S.A. A density functional study of the reactivity descriptors and antioxidant behavior of crocin. Comp. Theor. Chem. 2013, 1013, 123–129. [Google Scholar] [CrossRef]
  19. Kordzadeh, A.; Saadatabadi, A.R.; Adi, A. Investigation on penetration of saffron components through lipid bilayer bound to spike protein of SARS-CoV-2 using steered molecular dynamics simulation. Heliyon 2020, 6, e05681. [Google Scholar] [CrossRef]
  20. Saffari, B.; Amininasab, M. Crocin inhibits the fibrillation of human α synuclein and disassembles maturefFibrils: Experimental findings and mechanistic insights from molecular dynamics Simulation. ACS Chem. Neurosci. 2021, 12, 4037–4057. [Google Scholar] [CrossRef]
  21. Martinez, A.B. A. antiradical power of carotenoids and vstamin E: Testing the hydrogen atom transfer mechanism. J. Phys.Chem. B 2008, 112, 16945–16951. [Google Scholar] [CrossRef]
  22. Ceron-Carrasco, J.P.; Bastida, A.; Requena, A.; Zuniga, J.; Miguel, B. A theoretical study of the reaction of β-carotene with the nitrogen dioxide radical in solution. J. Phys. Chem. B 2010, 114, 4366–4372. [Google Scholar] [CrossRef] [PubMed]
  23. Cerezo, J.Z.; Bastida, A.; Requena, A.; Ceron-Carrasco, J.P.; Eriksson, L.A. Antioxidant properties of β-carotene isomers and their role in photosystems: Insights from ab initio simulation. J. Phys. Chem. A 2012, 116, 3498–3506. [Google Scholar] [CrossRef] [PubMed]
  24. Liao, F.X.; Hu, C.H. A density functional theory study for the role of end groups on the antioxidative potency of carotenoids. Theor. Chem. Acc. 2013, 132, 1357. [Google Scholar] [CrossRef]
  25. Monego, D.L.; Barcellos da Rosa, M.; do Nascimento, P.C. Applications of computational chemistry to the study of the antiradical activity of carotenoids: A review. Food Chem. 2017, 217, 37–44. [Google Scholar] [CrossRef] [PubMed]
  26. Sandhiya, L.; Zipse, H. Conformation-dependent antioxidant properties of β-carotene. Org. Biomol. Chem. 2022, 20, 152–162. [Google Scholar] [CrossRef]
  27. Tzeng, Y.Z.; Hu, C.H. Radical-induced cis-trans isomerization of fatty acids: A theoretical study. J. Phys. Chem. A 2014, 118, 4554–4564. [Google Scholar] [CrossRef]
  28. Christodoulou, M.C.; Orellana Palacios, J.C.; Hesami, G.; Jafarzadeh, S.; Lorenzo, J.M.; Domínguez, R.; Moreno, A.; Hadidi; M. Spectrophotometric methods for Measurement of antioxidant activity in food and pharmaceuticals. Antioxidants 2022, 11, 2213. [Google Scholar] [CrossRef]
  29. Belloni, J.M.; Monard, H.; Gobert, F.; Larbre, J.-P.; Demarque, A.; De Waele, V.; Lampre, I.; Marignier, J.-L.; Mostafavi, M.; Bourdon, J.-C.; et al. ELYSE. A picosecond electron accelerator for pulse radiolysis research. Nucl. Instrum. Meth. Phys. Res. Sect. A 2005, 539, 527–539. [Google Scholar] [CrossRef]
  30. Marignier, J.-L.; De Waele, V.; Monard, H.; Gobert, F.; Larbre, J.-P.; Demarque, A.; Mostafavi, M.; Belloni, J. Time-resolved spectroscopy at the picosecond laser-triggered electron accelerator ELYSE. Radiat. Phys. Chem. 2006, 75, 1024–1033. [Google Scholar] [CrossRef]
  31. Torche, F.; Marignier, J.-L. Direct evaluation of the molar absorption coefficient of hydrated electron by the isosbestic point method. J. Phys. Chem. B 2016, 120, 7201–7206. [Google Scholar] [CrossRef]
  32. Wang, F.; Schmidhammer, U.; Larbre, J.P.; Zong, Z.; Marignier, J.-L.; Mostafavi, M. Time-dependent yield of the hydrated electron and the hydroxyl radical in D2O: A picosecond pulse radiolysis study. Phys. Chem. Chem. Phys. 2018, 20, 15671–15679. [Google Scholar] [CrossRef]
  33. Jacquemin, D.; Brémond, E.; Planchat, A.; Ciofini, I.; Adamo, C. TD-DFT vibronic couplings in anthraquinones : From basis set and functional benchmarks to applications for industrial dyes. J. Chem. Theory Comput. 2011, 7, 1882–1892. [Google Scholar] [CrossRef]
  34. Brémond, E.; Alberto, M.E.; Russo, N.; Ciofini, I.; Adamo, C. Photophysical properties of NIR-emitting fluorescence probes: Insights from TD-DFT. Phys. Chem. Chem. Phys. 2013, 15, 10019. [Google Scholar] [CrossRef]
  35. Ungerer, P.; Tavitian, B.; Boutin, A. Applications of Molecular Simulations in the Oil and Gas Industry: Monte-Carlo Methods; Editions Technip: Paris, France, 2005. [Google Scholar]
  36. Vydrov, O.A.; Scuseria, G.E. Assessment of a long-range corrected functional. J. Chem. Phys. 2006, 125, 234109. [Google Scholar] [CrossRef]
  37. Frisch, M.J.; Trucks, G.W.; Schlegel, H.B.; Scuseria, G.E.; Robb, M.A.; Cheeseman, J.R.; Scalmani, G.; Barone, V.; Petersson, G.A.; Nakatsuji, H.; et al. Gaussian 16, Rev B.01; Gaussian, Inc.: Wallingford, CT, USA, 2016. [Google Scholar]
  38. Lever, A.B.; Gorelski, S.I. Ruthenium complexes of non-innocent Ligands: Aspects of charge transfer spectroscopy. Struct. Bond. 2004, 107, 77–114. [Google Scholar]
  39. Wang, F.; Pernot, P.; Marignier, J.-L.; Archirel, P.; Mostafavi, M. Mechanism of (SCN)2.- Formation and decay in neutral and basic KSCN solutions upon irradiation from a pico to microsecond range. J. Phys. Chem. B 2019, 123, 6599–6608. [Google Scholar] [CrossRef] [PubMed]
  40. Archirel, P.; Houée Levin, C.; Marignier, J.-L. Radiolytic oxidation of two inverse dipeptides, methionine-valine and valine-methionine: A Joint experimental and computational study. J. Phys. Chem. B 2019, 123, 9087–9097. [Google Scholar] [CrossRef] [PubMed]
  41. Léonard, C.; Le Quéré, F.; Adjei, D.; Denisov, S.A.; Mostafavi, M.; Archirel, P. Oxidation of silver cyanide Ag(CN)2- by the OH radical: From ab initio calculations to molecular simulations, and to experiment. J. Phys. Chem. A 2020, 124, 10787–10798. [Google Scholar] [CrossRef]
  42. Ipatov, A.; Cordova, F.; Doriol, F.J.; Casida, M.E. Excited states spin contamination in time-dependent density-functional theory for molecules with open-shell ground state. J. Mol. Struct. THEOCHEM 2009, 914, 60–73. [Google Scholar] [CrossRef]
  43. Wertz, D.H. Relationship between the gas phase entropies of molecules and their entropies of solvation in water. J. Am. Chem. Soc. 1980, 102, 5316–5322. [Google Scholar] [CrossRef]
  44. Lattach, Y.; Archirel, P.; Remita, S. Influence of chemical functionalities of a molecular imprinted polymers on its sensing properties: Electrochemical measurements and semiempirical DFT calculations. J. Phys. Chem. B. 2012, 116, 1467–1481. [Google Scholar] [CrossRef]
  45. Dennington, R.D.; Keith, T.A.; Millam, J.M. GaussView 6.1.1; Semichem Inc.: Shawnee, KS, USA, 2019. [Google Scholar]
  46. Kuhn, R.; Wang, Y. Synthese des Tetradecacetylcrocins und verwandter Verbindungen. Ber. Der Dtsch. Chem. Ges. 1939, 72, 871–878. [Google Scholar] [CrossRef]
  47. Cvetković, D.; Marković, D. Stability of carotenoids toward UV-irradiation in hexane solution. J. Serb. Chem. Soc. 2008, 73, 15–27. [Google Scholar] [CrossRef]
Figure 1. (Top) The all-trans crocin molecule. (Bottom) The all-trans Me2–crocetin molecule used as a model for the molecular simulations. The numbering of the carbon atoms is indicated on the polyene chain (between the dotted lines).
Figure 1. (Top) The all-trans crocin molecule. (Bottom) The all-trans Me2–crocetin molecule used as a model for the molecular simulations. The numbering of the carbon atoms is indicated on the polyene chain (between the dotted lines).
Antioxidants 12 01202 g001
Figure 2. (a) Transient difference spectra within 1 μs after the pulse for ([croc] = 0.0385 mM). (b) Kinetics at 333, 375, 441, 510 et 678 nm. Dose = 97 Gy. Optical path: 0.5 cm.
Figure 2. (a) Transient difference spectra within 1 μs after the pulse for ([croc] = 0.0385 mM). (b) Kinetics at 333, 375, 441, 510 et 678 nm. Dose = 97 Gy. Optical path: 0.5 cm.
Antioxidants 12 01202 g002
Figure 3. (a) Normalized optical absorption spectra in the molar absorption coefficient of the crocin (in black) and of the oxidized radical of crocin (in pink). (b) Zoom of both spectra between 500 and 750 nm, with an isosbestic point at 510 nm.
Figure 3. (a) Normalized optical absorption spectra in the molar absorption coefficient of the crocin (in black) and of the oxidized radical of crocin (in pink). (b) Zoom of both spectra between 500 and 750 nm, with an isosbestic point at 510 nm.
Antioxidants 12 01202 g003
Figure 4. (a) Kinetics at 678 nm within 84 μs for [croc] = 0.13 mM, reconstituted according to the absorbance at 1 μs (Figure S3) and to the signals over 10 and 100 μs (insets). (b) Transient spectra within 84 μs for ([croc] = 0.13 mM. Inset: Second-order test.
Figure 4. (a) Kinetics at 678 nm within 84 μs for [croc] = 0.13 mM, reconstituted according to the absorbance at 1 μs (Figure S3) and to the signals over 10 and 100 μs (insets). (b) Transient spectra within 84 μs for ([croc] = 0.13 mM. Inset: Second-order test.
Antioxidants 12 01202 g004
Figure 5. (a) Time-resolved spectra within 20 and 418 μs. The spectrum at 0.86 μs (in pink) of Figure 2a is added for comparison. Inset: Kinetics at 441 nm and 678 nm. Data before 0.86 μs are taken from Figure 2b. (b) Normalized optical absorption spectra in t molar absorption coefficient of the crocin (in black) and of the radical dimer (croc (–H))2 (in red).
Figure 5. (a) Time-resolved spectra within 20 and 418 μs. The spectrum at 0.86 μs (in pink) of Figure 2a is added for comparison. Inset: Kinetics at 441 nm and 678 nm. Data before 0.86 μs are taken from Figure 2b. (b) Normalized optical absorption spectra in t molar absorption coefficient of the crocin (in black) and of the radical dimer (croc (–H))2 (in red).
Antioxidants 12 01202 g005
Figure 6. (a) Dose-dependent optical absorption spectra of a γ−irradiated crocin solution. Inset: Variation of the absorbances at 330, 375, and 441 nm versus the dose. [croc] = 0.057 mM. Dose rates = 3.5, 6.0, 9.8, or 21.1 Gy/min. (b) Spectra of the crocin (in black) and of its final oxidation product (in green) arising by γ-irradiation and obtained from the deconvolution, both calibrated in molar absorption coefficients.
Figure 6. (a) Dose-dependent optical absorption spectra of a γ−irradiated crocin solution. Inset: Variation of the absorbances at 330, 375, and 441 nm versus the dose. [croc] = 0.057 mM. Dose rates = 3.5, 6.0, 9.8, or 21.1 Gy/min. (b) Spectra of the crocin (in black) and of its final oxidation product (in green) arising by γ-irradiation and obtained from the deconvolution, both calibrated in molar absorption coefficients.
Antioxidants 12 01202 g006
Figure 7. Schemes of the pulse and gamma radiolysis of aqueous solutions of crocin saturated by N2O. The yield values are in 10−7 mole J−1 units.
Figure 7. Schemes of the pulse and gamma radiolysis of aqueous solutions of crocin saturated by N2O. The yield values are in 10−7 mole J−1 units.
Antioxidants 12 01202 g007
Figure 8. (a) Absorption spectra evolution with the γ-dose of a mixed solution of 0.084 mM of crocin and 0.1 M of silver perchlorate saturated by N2O. (b) Scheme of the radiolysis of the solutions. The yield values are in 10−7 mole J−1 units.
Figure 8. (a) Absorption spectra evolution with the γ-dose of a mixed solution of 0.084 mM of crocin and 0.1 M of silver perchlorate saturated by N2O. (b) Scheme of the radiolysis of the solutions. The yield values are in 10−7 mole J−1 units.
Antioxidants 12 01202 g008
Figure 9. Model 1 for the evaluation of transition states for H atom abstraction by the OH radical at the B3LYP/apvdz/SMD level. Barrier heights (eV) and imaginary frequencies (cm−1), written with a negative sign. Bracketed numbers for Ca and Cb are approximate but the negative sign is ensured (see text).
Figure 9. Model 1 for the evaluation of transition states for H atom abstraction by the OH radical at the B3LYP/apvdz/SMD level. Barrier heights (eV) and imaginary frequencies (cm−1), written with a negative sign. Bracketed numbers for Ca and Cb are approximate but the negative sign is ensured (see text).
Antioxidants 12 01202 g009
Figure 10. Absorption spectra of Me2 crocetin, of its RCa and RCb radicals and their sum with two coefficients, with fwhm = 0.50 eV, spin screening 20% and (a) the small pSDD and (b) the larger pSDD+ gaussian bases.
Figure 10. Absorption spectra of Me2 crocetin, of its RCa and RCb radicals and their sum with two coefficients, with fwhm = 0.50 eV, spin screening 20% and (a) the small pSDD and (b) the larger pSDD+ gaussian bases.
Antioxidants 12 01202 g010
Figure 11. (a) Absorption spectrum of a few Me2–crocetin conformers, yielded by MC simulations with all-trans and bent constraints; with fwhm = 0.50 eV. (b) Measured absorption spectrum of crocin (Figure S1) with structure propositions according to the Me2–crocetin results of Figure 11a.
Figure 11. (a) Absorption spectrum of a few Me2–crocetin conformers, yielded by MC simulations with all-trans and bent constraints; with fwhm = 0.50 eV. (b) Measured absorption spectrum of crocin (Figure S1) with structure propositions according to the Me2–crocetin results of Figure 11a.
Antioxidants 12 01202 g011
Figure 12. Simulated spectra of the RCa and RCb species with the large pSDD+ basis. The focus is made on weak bands at (a) higher energies with fwhm = 0.50 eV, and (b) at lower energies with fwhm = 0.20 eV and spin screening 20%.
Figure 12. Simulated spectra of the RCa and RCb species with the large pSDD+ basis. The focus is made on weak bands at (a) higher energies with fwhm = 0.50 eV, and (b) at lower energies with fwhm = 0.20 eV and spin screening 20%.
Antioxidants 12 01202 g012
Figure 13. Despite the electrostatic attraction, the OH radical cannot oxidize the sugar moieties of crocin. This sugar-driven mechanism favors the oxidation of the neighbor –CaH3, rather than –CbH3.
Figure 13. Despite the electrostatic attraction, the OH radical cannot oxidize the sugar moieties of crocin. This sugar-driven mechanism favors the oxidation of the neighbor –CaH3, rather than –CbH3.
Antioxidants 12 01202 g013
Figure 14. Absorption spectra of the covalent dimers through simulation of frozen monomers, with fwhm = 0.70 eV. The focus is made on (a) the intense bands, and (b) the weak, high-energy bands. The labels (½ ½) and (5/6 1/6) refer to the monomer probabilities.
Figure 14. Absorption spectra of the covalent dimers through simulation of frozen monomers, with fwhm = 0.70 eV. The focus is made on (a) the intense bands, and (b) the weak, high-energy bands. The labels (½ ½) and (5/6 1/6) refer to the monomer probabilities.
Antioxidants 12 01202 g014
Figure 15. Absorption spectra of (a) the four possible cyclic products of the radical’s disproportionation, and (b) three of their combinations, with fwhm = 0.50 eV.
Figure 15. Absorption spectra of (a) the four possible cyclic products of the radical’s disproportionation, and (b) three of their combinations, with fwhm = 0.50 eV.
Antioxidants 12 01202 g015
Table 1. Radiolytic yields (in 10−7 mole J−1) and molar absorption coefficients (in 105 M−1 cm−1) of the crocin oxidation products ([croc] = 0.13 mM).
Table 1. Radiolytic yields (in 10−7 mole J−1) and molar absorption coefficients (in 105 M−1 cm−1) of the crocin oxidation products ([croc] = 0.13 mM).
SolutionCrocin
(Pulse Radiolysis)
Crocin
(γ-Rays)
Crocin + Ag+
(γ-Rays)
Gexp(- croc)-4.02.1
Gcalc(- croc)6.2 3.84 2.5
Gcalc(croc –(H)6.2 --
ε678nm(croc –(H))
ε441nm(croc –(H))
ε333nm(croc –(H))
0.0623
1.17
0.10
--
Gcalc(croc –(H))23.1 --
ε333nm(croc –(H))2
ε441nm(croc –(H))2
ε678nm(croc –(H))2
0.31
1.84
0

-

-
Gcalc(croc (–2H))-3.842.5
ε333nm(croc (–2H))
ε441nm(croc (–2H))
-0. 26
0.13

-
Table 2. Molar reaction free energies (in eV) for processes occurring in the oxidation of Me2–crocetin by the OH radical. R = (Me2–crocetin (–H)) radical, with one oxidized carbon atom. The carbon atoms Ca and Cb are shown in Figure 1.
Table 2. Molar reaction free energies (in eV) for processes occurring in the oxidation of Me2–crocetin by the OH radical. R = (Me2–crocetin (–H)) radical, with one oxidized carbon atom. The carbon atoms Ca and Cb are shown in Figure 1.
ReactionrG (eV)
Me2–crocetin + OH→ RCa+ H2O−1.98
→ RCb−1.94
RCa + RCa→ dimer CaCa-−1.09
RCa + RCb→ dimer CaCb-−1.05
RCb + RCb→ dimer CbCb-−1.12
RCa + RCa→ Me2–crocetin (–2H) CaCa′ 1A+ Me2–crocetin−0.26
→ Me2–crocetin (–2H) CaCa′ cyclic−0.48
→ Me2–crocetin (–2H) CaCb• 3a+0.30
→ Me2–crocetin (–2H) CaCb• 1A≥+0.30
→ Me2–crocetin (–2H) CaCb cyclic−1.28
→ Me2–crocetin (–2H) CaCb′ 1A−0.26
→ Me2–crocetin (–2H) CaCb′ cyclic−0.39
RCa + RCb→ Me2–crocetin (–2H) CaCa′ 1A+ Me2–crocetin−0.30
→ Me2–crocetin (–2H) CaCa′ cyclic−0.52
→ Me2–crocetin (–2H) CaCb• 3a+0.26
→ Me2–crocetin (–2H) CaCb• 1A≥+0.26
→ Me2–crocetin (–2H) CaCb cyclic−1.32
→ Me2–crocetin (–2H) CaCb′ 1A−0.30
→ Me2–crocetin (–2H) CaCb′ cyclic−0.43
→ Me2–crocetin (–2H) CbCb′ 1A−0.29
→ Me2–crocetin (–2H) CbCb′ cyclic−1.04
RCb + RCb→ Me2–crocetin (–2H) CaCb• 3a+ Me2–crocetin+0.22
→ Me2–crocetin (–2H) CaCb• 1A≥+0.22
→ Me2–crocetin (–2H) CaCb cyclic−1.36
→ Me2–crocetin (–2H) CaCb′ 1A−0.33
→ Me2–crocetin (–2H) CaCb′ cyclic−0.47
→ Me2–crocetin (–2H) CbCb′ 1A−0.33
→ Me2–crocetin (–2H) CbCb′ cyclic−1.08
RCa + Me2–crocetin(RCa, Me2–crocetin)-−0.52
RCb + Me2–crocetin(RCb, Me2–crocetin)-−0.52
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Gharib, S.A.; Archirel, P.; Adjei, D.; Belloni, J.; Mostafavi, M. Unveiling the Intimate Mechanism of the Crocin Antioxidant Properties by Radiolytic Analysis and Molecular Simulations. Antioxidants 2023, 12, 1202. https://doi.org/10.3390/antiox12061202

AMA Style

Gharib SA, Archirel P, Adjei D, Belloni J, Mostafavi M. Unveiling the Intimate Mechanism of the Crocin Antioxidant Properties by Radiolytic Analysis and Molecular Simulations. Antioxidants. 2023; 12(6):1202. https://doi.org/10.3390/antiox12061202

Chicago/Turabian Style

Gharib, Sarah Al, Pierre Archirel, Daniel Adjei, Jacqueline Belloni, and Mehran Mostafavi. 2023. "Unveiling the Intimate Mechanism of the Crocin Antioxidant Properties by Radiolytic Analysis and Molecular Simulations" Antioxidants 12, no. 6: 1202. https://doi.org/10.3390/antiox12061202

APA Style

Gharib, S. A., Archirel, P., Adjei, D., Belloni, J., & Mostafavi, M. (2023). Unveiling the Intimate Mechanism of the Crocin Antioxidant Properties by Radiolytic Analysis and Molecular Simulations. Antioxidants, 12(6), 1202. https://doi.org/10.3390/antiox12061202

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop