Shamima
Begum
a,
Michele
Cianci
b,
Bo
Durbeej
c,
Olle
Falklöf
c,
Alfons
Hädener‡
a,
John R.
Helliwell
*a,
Madeleine
Helliwell
a,
Andrew C.
Regan
a and
C.
Ian F. Watt
a
aSchool of Chemistry, University of Manchester, Manchester, M13 9PL, UK. E-mail: john.helliwell@manchester.ac.uk
bEMBL, PETRA III, DESY, Notkestrasse 85, 22603 Hamburg, Germany
cDivision of Computational Physics, IFM, Linköping University, SE-58183 Linköping, Sweden
First published on 13th March 2015
The chemical basis of the blue-black to pink-orange color change on cooking of lobster, due to thermal denaturation of an astaxanthin–protein complex, α-crustacyanin, in the lobster carapace, has so far been elusive. Here, we investigate the relaxation of the astaxanthin pigment from its bound enolate form to its neutral hydroxyketone form, as origin of the spectral shift, by analyzing the response of UV-vis spectra of a water-soluble 3-hydroxy-4-oxo-β-ionone model of astaxanthin to increases in pH, and by performing extensive quantum chemical calculations over a wide range of chemical conditions. The enolization of astaxanthin is consistent with the X-ray diffraction data of β-crustacyanin (PDB code: 1GKA) whose crystals possess the distinct blue color. We find that enolate formation is possible within the protein environment and associated with a large bathochromic shift, thus offering a cogent explanation for the blue-black color and the response to thermal denaturation and revealing the chemistry of astaxanthin upon complex formation.
The crystal structure of β-crustacyanin reveals two protein sub-units linked by a pair of astaxanthin molecules,6 modeled therein as α-hydroxyketone tautomers 1. Terminal oxygenated rings are in close association, at one end of each astaxanthin, with the nitrogen of a histidine of one protein and, at the other end, with a water molecule and the acidic phenolic hydroxyl of a tyrosine residue of the second protein. These rings are almost coplanar with the polyene chains, which are slightly bowed with a closest approach of about 7 Å.
Two distinct themes have developed in explanation for the coloration. The first, and most thoroughly investigated, focuses on physical arrangements of the astaxanthins in the protein complexes. The angles of twist about the C6–C7 (and C6′–C7′) bonds of the astaxanthins are much reduced from those found in the 6-s-cis conformation of uncomplexed astaxanthin (ca. 50°), to a 6-s-trans conformation with end rings coplanar with the polyene chains in the protein bound form, permitting enhanced overlap along the conjugated array in the complexed forms. Exciton coupling between the paired astaxanthins in β-crustacyanin has also been considered. Quantum chemical calculations and experiment, however, suggest that neither conformational change nor exciton coupling can account for more than 30% of the observed shift, with the larger contribution arising from co-planarization of the end rings.8–10 Kuhn and Sörensen suggested in 1938 a second theme,11 which remained largely neglected, involving a reversible ionization of the astaxanthins upon complexation. The consequences of deprotonation of an α-hydroxycyclohexenone residue, and subsequent proton transfers, are shown in Scheme 1. However, Buchwald and Jencks12 stated that the proteins involved did not contain a site sufficiently basic to retain a proton.
The viability of enolate formation, however, also depends on the acidity of astaxanthin and although this remains unquantified, we note that much simpler compounds, also capable of forming an α-hydroxylated enolate in which the negative charge is distributed between two oxygen atoms, including 2,3-dihydroxy-2-propenal (triose reductone) and ascorbic acid,13,14 are relatively strong organic acids (pKa(1) = 5.2 and 4.1, respectively) with second ionizations (pKa(2) ≈ 13) also accessible in water. A priori rejection of the possibility of enolate formation within the protein thus seems premature. While solutions of astaxanthin in a strongly basic organic medium are known to develop deep blue colors,15 quantitative studies in aqueous media require disentangling the effects of aggregation16,17 from those of ionization. We have therefore examined the behavior of 3-hydroxy-4-oxo-β-ionone, 2, C13H18O3, as a water-soluble model that retains the trimethylated α-hydroxycyclohexenone of astaxanthin but replaces the rest of the molecule with a simpler α-β-unsaturated ketone at the 6-position of the ring. This permits formal conjugation of the ring carbonyl (at C4) with that on the side chain at C9 through the C6–C7 bond, which occupies the same steric environment as the corresponding bond in astaxanthin. This compound, readily synthesized from β-ionone,18,19 is shown by NMR spectroscopy to exist as its ketonic tautomeric form in both methanol and chloroform solution. Its enolization and enolate chemistry, expected to mimic that of astaxanthin (Scheme 1), is shown in Scheme 2. Our UV-vis experiments on an astaxanthin model compound and quantum chemical calculations support a mechanism by which the colors of crustacyanin protein complexes arise because astaxanthin is bound in its 1− form.
The DFT and MP2 calculations were performed with the GAUSSIAN 09 program.32 The CC2 and CASPT2 calculations, in turn, were performed with the TURBOMOLE 6.333,34 and MOLCAS 7.635 programs, respectively. Further details on the quantum chemical calculations are given in the ESI.†
Spectra of the aqueous solutions were also unchanged over long periods after acidification with dilute hydrochloric acid. Addition of aqueous NaOH solution, with cells maintaining a nitrogen atmosphere over the solution, enabled preparation and monitoring of basic solutions of 2 in water containing 0.7% (vol/vol) methanol. A typical sequence of spectra, following addition of 2 to a 0.2 M aqueous NaOH solution, is shown in Fig. 1a.
In Fig. 1a, trace A is the initial spectrum of 2 with λmax at 275 nm. Over about 45 minutes, this spectrum evolves to trace B, by decay of the peak at 275 nm and development of new maxima at 278, 356 and 454 nm (the longest-wavelength new maximum thus presents a bathochromic shift of no less than 180 nm relative to the parent material). The sequence displays a tight isosbestic point at 308 nm, and sections at any of the new maxima show that the absorbance changes are first order, with a rate constant kobs = 9.7 × 10−4 s−1. Spectra at this stage were stable for at least one hour, with slow decays associated with leakage of atmospheric oxygen into the cell (see further below). At this stage also, acidification of the solution to pH 4 by addition of the calculated amount of acetic acid quickly regenerated, with no detectable intermediate stages, the original spectrum (trace A). Furthermore, chromatography of an extract confirmed the recovery of 2 and showed only minor traces of other materials.
The spectroscopic changes, A → B (Fig. 1a), are consistent with the formation of either of the anions 2− and 2-dianion, or indeed of a rapidly equilibrating mixture thereof. Properly distinguishing between these possibilities led to examining spectra after equilibration in buffered basic solutions over as wide a pH range as possible. To avoid prolonged equilibration times and reproducibility issues due to slow infusion of atmospheric oxygen into the media, the solutions were equilibrated in 0.2 M NaOH, (spectrum B in Fig. 1a), and then the spectra were examined after adjustment of pH in both acidic and basic directions.
The results of increasing the basicity of the medium are shown in Fig. 1b. While the original trace A is included for comparison, the starting point for the chemistry here is trace B, corresponding to equilibrated 2− and 2-dianion species in a 0.2 M NaOH solution at pH 13.2. Traces C and D arose from sequential additions of highly concentrated NaOH, increasing the basicity firstly to pH 13.3 and then to pH 13.7. In both cases, spectra were corrected for small volume changes associated with the additions. In these spectra, the peaks at 278 and 454 nm increase with increasing basicity, whereas the peak at 356 nm decreases, with a new isosbestic point appearing at 395 nm. These changes are consistent with a situation in which, in trace B, the peak at 356 nm is associated with 2− and the peaks at 278 and 454 nm with the 2-dianion species. In a further experiment, 2 was equilibrated directly with 2 M NaOH (pH 14.5), which produced a spectrum that after scaling to match the 395 nm isosbestic point is shown as trace E in Fig. 1b. Despite the large variations in ionic strength, trace E falls into the pattern shown by traces B, C and D. A near superimposable trace was again produced when 1.5 M NaOH was used, indicating that at pH 14.5, 2 is almost fully converted to its 2-dianion form.
The results of decreasing the basicity of the medium are shown in Fig. 1c. Again, the starting point is the equilibrated state of 2− and 2-dianion species (trace B). Addition of measured amounts of neat trifluoroethanol, pKa = 12.4,37 then provided solutions with some buffering at pH > 12. In trace F (pH = 12.6), the original peak at 454 nm associated with 2-dianion decreased substantially, while the original peak at 356 nm associated with 2− increased slightly. Trace G (pH = 12.1) shows some small depletion of 2− and almost complete loss of the 2-dianion. For a more acidic medium yet, a portion of solid sodium bicarbonate, calculated to produce a solution buffered at pH 10.4, was added to a freshly equilibrated solution of 2 in 0.2 M aqueous NaOH. After small scaling of the spectrum to match the isosbestic point, this yielded trace H, in which it is clear that the 356 nm peak has dropped to about half its original value and there has been substantial growth of the 272 nm peak associated with 2 itself. Further acidifications continue the trend of reversion to the spectrum of 2.
The spectra and responses to pH changes presented in Fig. 1a–c are fully consistent with the scheme of deprotonations and proton transfers shown above in Scheme 2, where: (i) deprotonation of 2 at C3 followed by proton transfer from O3 to O4 yield 2−, with a sharing of the negative charge by conjugation between O3 and O9; (ii) reprotonation at O3 yields the ene-diol tautomer 2ed; and (iii) deprotonation of the enolic hydroxyl group of 2− produces 2-dianion. The spectra also permit initial estimates of the dissociation constants Ka(1) and Ka(2) defined in Scheme 2. From Fig. 1c, a comparison of traces H and F suggests that at pH 10.4, there is a ca. 50% conversion of 2 to 2−, with little further conversion to the 2-dianion. In that case, application of the Henderson–Hasselbach equation to the equal concentrations of 2 and 2− gives pKa(1) ≈ 10.4. Similar inspection of Fig. 1b, comparing traces B and E, suggests that at pH 13.2, the equilibrated solution contains near equal amounts of 2− and 2-dianion, with little residual un-ionized 2, so that pKa(2) ≈ 13.2. These values of the dissociation constants might be compared with the corresponding values of 4.12 (pKa(1)) and 11.51 (pKa(2)) for ascorbic acid.38
The preceding analysis providing estimates of Ka(1) and Ka(2) does not take into account the possibility that 2 can also be ionized to form its alkoxide, by deprotonation of the secondary alcohol (Scheme 3a), a process expected to have little effect on the UV-vis spectrum of 2. However, since this process involves equilibration of an oxygen acid, the associated relaxation in either acid or base in aqueous medium should be fast39 compared to that for the deprotonation leading to 2−, irrespective of thermodynamic acidities. The value of 10.4 for pKa(1) extracted above is not then simply that for dissociation of 2 as a carbon acid. Rather, any appreciable alkoxide formation would call for a correction according to the equation included in Scheme 3a, whereby the apparent ionization constant will be an underestimate of the true value, i.e., pKa(C3–H) < pKa(1).
For simple secondary aliphatic alcohols (16 < pKa(O3–H) < 16.5), extents of alkoxide formation would be small, even in strongly basic aqueous solutions, and the discrepancy between apparent and actual carbon acid dissociation constants would be undetectably small. The secondary alcohol in 2, however, is in a polar environment, and empirical additivity relationships (based on a very sparse data set) have yielded a relationship40 between acidity and Taft's polar substituent parameters (σ*) which, using σ* = 2.0 for –COCH = CH2,41 gives pKa(O3–H) ≈ 13.1 for 2. Even with this value, the difference between apparent and actual carbon acid dissociation constants would be small.
However, the uncertainty in the estimate of pKa(O3–H) is large, but while we have no direct experimental measurement, the kinetics of formation of 2− indicates that the value may be much higher than 13.1. This is also found by the quantum chemical calculations further described below. Earlier, we noted that the formation of 2− and 2-dianion mixture in 0.2 M NaOH show first-order behavior with kobs = 9.7 × 10−4 s−1. Assuming that the reaction is also first order in hydroxide, the second-order rate constant for the reaction is then 4.9 × 10−3 M−1 s−1. Comparison with the second-order rate constants for deprotonation (by hydroxide) at 25 °C of more familiar carbon acids such as nitromethane (pKa = 10.2) and acetylacetone (pKa = 8.9) of 2.76 × 10 M−1 s−1 (ref. 42) and 3.98 × 104 M−1 s−1,43 respectively, makes it very clear that carbon deprotonation of 2 is >103 slower than for a carbon acid of similar acidity. Nothing obvious in the structure of 2 accounts for the low reactivity, but if, under the conditions, 2 existed largely as its alkoxide and did not rearrange or react with hydroxide to yield 2−, the rate constant found here would need to be scaled upwards as shown in Scheme 3b to yield a true value, kC3–H, for the second-order rate constant for carbon deprotonation of 2 by hydroxide. If the correction factor kC3–H/kobs is indeed ca. 103, it is possible to work out an estimate of 10.0 for the true pKa(C3–H) value of 2 (cf. with the above estimate of 10.4 for the pKa(1) value of 2).
To identify the oxidation products, a large-scale reaction was carried out. Chromatography showed formation of a single new product. After acidification of the solution, a crystalline product was isolated, molecular formula C13H16O3, thus a simple oxidation product of 2. Revealingly, proton NMR spectroscopy confirmed loss of the ABX signal, associated with the C2 and C3 hydrogens of 2 and its replacement by a single one-proton signal at δ 6.04. Geminal methyl groups now show a single signal indicating equivalent environments. The spectrum also shows a broad 1H signal at δ 6.28, assigned to a hydrogen bonded OH, whose presence was confirmed by IR spectroscopy. The analytical data are consistent with this oxidation product being the enol 2,3-dehydro-3-hydroxy-4-oxo-β-ionone (3 in Scheme 3c), bearing the same relationship to 2 as astacene bears to astaxanthin (with which it often co-occurs1). The UV-vis spectra of 3 in neutral and basic aqueous solution were superimposable on traces J and I in Fig. 1d. Responses of spectra to changes in pH were immediate, and the behavior in borate buffers gave pKa = 9.9. The chemistry is summarized in Scheme 3c in which we suggest that atmospheric oxygen reacts with 2-dianion, yielding the enolate of 3 and probably hydrogen peroxide, as its anion.
In light of this discrepancy, we have quite exhaustively investigated whether factors that may have a negative impact on the accuracy of TD-DFT excitation energies are at play. This investigation, which is fully described in Sections 4.1–4.5 of the ESI† and presented in Tables S1 and S3–S9, includes a comparison of TD-DFT results with CC2 and CASPT2 ab initio results, an assessment of how molecular geometries and choice of basis set influence the calculated λmax, and additional benchmark calculations on the oxidation product of 2 (2,3-dehydro-3-hydroxy-4-oxo-β-ionone) and on 2,3-dihydroxy-2-propenal, which contains the minimum carbon framework in which conjugation between a carbonyl group and an ene-diol may occur. Overall, this investigation shows that TD-DFT calculations constitute a reliable tool for obtaining a qualitatively accurate assessment of how acid–base chemistry influences the UV-vis spectroscopic behavior of astaxanthin and models thereof, and strongly indicates that more elaborate ab initio calculations would support the same conclusions.
Describing the solvent using a hybrid cluster-continuum approach,45 we also calculated the pKa dissociation constants of 2 (Table S2, ESI†). Since most of these pKa values are available from experiment, the calculations serve both as a test of the experimental predictions and a validation of the computational accuracy of DFT for the present systems, with implications also for the excited-state calculations discussed below.
Firstly, it can be seen that the experimentally measured pKa of 10.4 for deprotonation of the hydroxyketone at C3 (pKa(C3–H) in Table S2 (ESI†) and Scheme 3; pKa(1) Scheme 2) to form the mono-anion is reproduced very well, and that the calculated estimates (ranging between 7.7 and 10.6) are not particularly sensitive to the number of explicit water molecules included in the modeling. Given that we propose a mechanism for the bathochromic shift in crustacyanin based on formation of the mono-anion of astaxanthin, it is also notable that this pKa is fairly small. From calculations on the ketone form lacking the C3 hydroxyl group (4-oxo-β-ionone), it is predicted that the smallness of pKa(C3–H) for the hydroxyketone is due precisely to the C3 hydroxyl group, because without this group, the calculated pKa(C3–H) values amount to ∼20 or even higher.
The calculations further predict that pKa(C3–H) for the hydroxyketone is considerably smaller than the pKa for deprotonation of the C3 hydroxyl group to form the alkoxide of the hydroxyketone (pKa(O3–H) in Table S2, ESI†). Thus, alkoxide formation should not impact on mono-anion formation. Although the estimates of pKa(O3–H) vary somewhat with the number of water molecules, the difference between pKa(O3–H) and pKa(C3–H) appears to converge to a value of around 10 in the calculations.
In addition to pKa values of the hydroxyketone, we have also calculated the pKa for deprotonation of the C4 hydroxyl group of the mono-anion (pKa(O4–H) in Table S2 (ESI†); pKa(2) in Scheme 2) to form the di-anion. This pKa has also been determined experimentally, and pleasingly, the calculated values approach (23.6 → 17.9 → 14.5) the experimental estimate (13.2) to within close to 1 pKa unit as the number of water molecules is increased.
Starting with the keto form of astaxanthin, it can be seen from Table 1 that all three density functionals (M06-HF, LC-ωPBE and LC-BLYP) give, in the absence of effects linked to aggregation phenomena, very similar λmax for the water and methanol solvents, and that the experimental value (476 nm) pertaining to a water–methanol mixture is reproduced with ∼20 nm (∼0.1 eV) accuracy. It can also be seen that the relatively small 6-31G(d,p) basis set gives similar λmax as the larger aug-cc-pVDZ basis set, which reflects the pronounced valence character of the strongly absorbing “Bu-like” excited state. Indeed, the diffuse functions included in the aug-cc-pVDZ basis set red shifts the λmax by no more than 13 nm (0.07 eV).
Method/basis setb | Solvent | Keto | Mono-enol | Di-enol | Mono-anion | Di-anion |
---|---|---|---|---|---|---|
a All calculations based on geometries optimized using B3LYP/6-31G(d,p) in combination with a CPCM description of the water or methanol solvent. b Basis set I = 6-31G(d,p); basis set II = aug-cc-pVDZ. c Solvent effects obtained using the cc-pVDZ basis set and added to gas-phase data obtained with the target aug-cc-pVDZ basis set. d Experimental value in a water–methanol mixture. e ε max given in units of dm3 mol−1 cm−1. | ||||||
M06-HF/I, λmax (ΔE, f) | H2O | 457 (2.71, 4.39) | 463 (2.68, 4.48) | 459 (2.70, 4.56) | 660 (1.88, 4.04) | 790 (1.57, 4.84) |
LC-ωPBE/I, λmax (ΔE, f) | H2O | 473 (2.62, 4.53) | 482 (2.57, 4.64) | 477 (2.60, 4.73) | 658 (1.88, 4.43) | 779 (1.59, 5.16) |
LC-BLYP/I, λmax (ΔE, f) | H2O | 463 (2.68, 4.59) | 470 (2.64, 4.70) | 465 (2.66, 4.79) | 643 (1.93, 4.45) | 754 (1.64, 5.19) |
M06-HF/II,cλmax (ΔE, f) | H2O | 470 (2.64, 4.31) | 476 (2.60, 4.37) | 471 (2.63, 4.44) | 666 (1.86, 4.17) | 816 (1.52, 4.67) |
LC-ωPBE/II,cλmax (ΔE, f) | H2O | 482 (2.57, 4.47) | 490 (2.53, 4.55) | 484 (2.56, 4.64) | 665 (1.86, 4.45) | 796 (1.56, 4.99) |
LC-BLYP/II,cλmax (ΔE, f) | H2O | 474 (2.62, 4.52) | 481 (2.58, 4.61) | 475 (2.61, 4.69) | 651 (1.90, 4.50) | 774 (1.60, 5.02) |
M06-HF/I, λmax (ΔE, f) | MeOH | 456 (2.72, 4.39) | 463 (2.68, 4.48) | 459 (2.70, 4.56) | 668 (1.86, 4.05) | 791 (1.57, 4.84) |
LC-ωPBE/I, λmax (ΔE, f) | MeOH | 473 (2.62, 4.52) | 481 (2.58, 4.64) | 477 (2.60, 4.73) | 666 (1.86, 4.45) | 779 (1.59, 5.16) |
LC-BLYP/I, λmax (ΔE, f) | MeOH | 463 (2.68, 4.59) | 470 (2.64, 4.70) | 465 (2.67, 4.79) | 651 (1.90, 4.46) | 754 (1.64, 5.19) |
Experimentdλmax (ΔE, εmaxe) | Mixture | 476 (2.61, 10000) | — | — | — | — |
Turning to the calculations on the other forms of astaxanthin also summarized in Table 1, but for which experimental data are not available, we first observe that neither single nor double enolization is predicted to significantly shift the λmax of the parent keto form, excluding tautomerism alone in the protein-bound astaxanthins as the origin of the much larger observed bathochromic shift. Notably, however, the calculations show that deprotonation of C3 followed by proton transfer between the oxygen atoms attached to C3 and C4 (or, equivalently, single enolization followed by deprotonation of the C3 hydroxyl group) to form the mono-anion, 1−, is associated with a red shift, lying uniformly between 177 and 212 nm (0.71–0.86 eV) for all methods, thus larger than the bathochromic shift in β-crustacyanin (λmax at 570 nm) and consistent with that in α-crustacyanin (λmax at 630 nm). For the dianion (1-dianion), the calculations yield a bathochromic shift of the order of 300 nm, bringing the absorption into the near-infrared (λmax ≈ 750–800 nm).
For this explanation to be valid chromophore–protein interactions should favor the enolate form over the α-hydroxyketone form. Thus, it would be of interest to investigate how specific short-range chromophore–protein interactions and long-range bulk electrostatic effects stabilize the various forms of astaxanthin inside the binding pocket of β-crustacyanin. While a detailed analysis of such interactions is beyond the scope of this work, it has previously been suggested, based on quantum chemical calculations,48 that the most proximal protein residues do not interact strongly with the neutral α-hydroxyketone form, which at least does not go against the idea that the protein shifts the equilibrium between enolate and α-hydroxyketone toward the former. Furthermore, this study also showed that the absorption of the α-hydroxyketone form is rather insensitive to short-range chromophore–protein interactions, which suggests that a more severe perturbation of the chromophore is needed to explain the bathochromic shift in the protein. Finally, preliminary B3LYP calculations that we have performed indicate that the histidine shifts the equilibrium between enolate and α-hydroxyketone by at least 8 kJ mol−1 toward the former species.
Astaxanthin is not unique amongst carotenoids in its ability to form colored complexes with the crustacyanin protein, but the overwhelming majority of color-forming carotenoids (including actinioerythrol, violerythrin, iso-norastaxanthin, astacene, 19,19′-dinorastacene, 4-hydroxy-4′-oxo-β-β-carotene, adonirubin, 7,8-didehydroastaxanthin and 15,15′-didehydroastaxanthin) possess structures with a fully conjugated polyene chain,2,49 linking oxygenated terminal rings with relatively acidic enolic or enolizable sites. Binding of carotenoids as an extended enolate within the protein may thus be a mechanism of some generality for formation of these spectacularly colored complexes in nature.
Footnotes |
† Electronic supplementary information (ESI) available: Complementary material for compound preparation and characterization and details and results of the quantum chemical calculations. See DOI: 10.1039/c4cp06124a |
‡ Pfaffenlohweg 29, CH-4125 Riehen, Switzerland. |
This journal is © the Owner Societies 2015 |