D.
Clara‡
a,
C. K.
Pezzei‡
a,
S. A.
Schönbichler
a,
M.
Popp
b,
J.
Krolitzek
b,
G. K.
Bonn
ac and
C. W.
Huck
*a
aInstitute of Analytical Chemistry and Radiochemistry, CCB – Center for Chemistry and Biomedicine, Leopold-Franzens University, Innrain 80/82, 6020 Innsbruck, Austria. E-mail: Christian.W.Huck@uibk.ac.at; Tel: +43 512 507 57304
bBionorica SE, Kerschensteinerstrasse 12-15, 92318 Neumarkt/Oberpfalz, Germany
cADSI – Austrian Drug Screening Institute, Innrain 66a, 6020 Innsbruck, Austria
First published on 26th August 2015
Near-infrared diffuse reflectance (NIR) and attenuated-total-reflectance mid-infrared (ATR-IR) spectroscopy techniques in hyphenation with multivariate analysis were utilized to determine the antioxidant capacity of ground Sambuci flos samples. Folin–Ciocalteu (FC), ferric ion reducing antioxidant power (FRAP), cupric reducing antioxidant capacity (CUPRAC), 2,2-diphenyl-picrylhydrazyl (DPPH) and 2,2′-azino-bis-(3-ethylbenzothiazoline-6-sulfonic acid) diammonium salt (ABTS) were optimized and performed as reference methods. To remove systematic errors several spectral pretreatments like 1st and 2nd derivative Savitzky–Golay, standard normal variate (SNV) or multiplicative scatter correction (MSC) were applied. Cross-validations and test-set validations were performed for all assays. The quality parameters, standard error of prediction (SEP) and the ratio performance deviation (RPD), were calculated. An acceptable quality of the calibration can be confirmed for ATR-IR spectroscopy (e.g. for the CUPRAC assay: R2: 0.85, RPDcorr: 2.68, SECV: 0.13% GAE for cross-validation; R2: 0.81, RPDcorr: 2.20, SEP: 0.15% GAE for test-set validation). Surprisingly all models calculated for NIR spectroscopy were of poor quality and point to unpredictability of the antioxidative capacity. Further investigations of extracts by high performance liquid chromatography (HPLC) with a diode array detector (DAD) coupled to mass spectroscopy (MS) were performed to analyze the principal compounds. Thus, rutin and chlorogenic acid were confirmed to be the main components in the samples. This study demonstrates that ATR-IR spectroscopy is suitable to determine the antioxidative capacity in ground Sambuci flos samples and can be used for quality control.
Antioxidants can lead to inhibition of the aforementioned effects and can be useful as preventers regarding oxidative molecules. Fruits, vegetables and medicinal herbs are the most important sources of antioxidant compounds.5 Herbal medicine is still the most used medicine in comparison to synthetic drugs due to its better cultural acceptability and lower side effects.5 In this study, spectroscopic methods to determine the antioxidant power of ground Sambucus nigra flowers (Sambuci flos) were developed. Sambuci flos belongs to the lineage of Adoxaceae (prior Caprifoliaceae), to the subgenus European elderberry.6,7 The most important antioxidants in Sambuci flos are flavonoids (flavonol glycosides, up to 3%) and phenolic acids (ca. 3%).8
Sambuci flos has a diaphoretic effect for the treatment of cold, chills and fever.8 Of great importance are the antiviral activity9 and the antibacterial activity against hospital pathogens.10
The antioxidant capacity of plants, food and chemical compounds can be measured by various assays.11 The various chemical mechanisms of these assays can be divided into two basic groups, the HAT (hydrogen atom transfer) or the SET (single electron transfer) mechanism.11,12 The currently used assays Folin–Ciocalteu (FC), ferric ion reducing antioxidant power (FRAP), cupric reducing antioxidant capacity (CUPRAC), 2,2-diphenyl-picrylhydrazyl (DPPH) and 2,2′-azino-bis-(3-ethylbenzothiazoline-6-sulfonic acid) diammonium salt (ABTS) function via the SET mechanism and the last two assays (DPPH and ABTS) additionally exhibit the HAT mechanism.11 These assays are destructive, time and chemical consuming. Near-infrared (NIR) and attenuated-total-reflectance infrared (ATR-IR) spectroscopy techniques measure the interaction of light with the material in a nondestructive way.13 NIR and ATR-IR spectroscopies in hyphenation with partial least squares (PLS) regression have already proved in a case study their capability to replace the aforementioned assay methods.14
:
9.5 with ethanol/H2O (50/50, v/v). 300 μl of diluted sample extract, 1600 μl H2O and 1100 μl FRAP reagent were mixed together. After incubation at 60 °C for 45 minutes and cooling for five minutes the absorbance was measured at 593 nm.
:
1
:
1 (v/v/v). This mixture (CUPRAC reagent) was allowed to rest at least for one hour. 50 μl sample extract, 775 μl H2O and 2250 μl CUPRAC reagent were mixed and incubated at 60 °C for 20 minutes. The absorbance was detected at 450 nm after the samples were cooled down for five minutes.
All measurements were performed in triplicate against a blank which contained the same amount of H2O instead of the sample in FC, FRAP and CUPRAC and the same amount of ethanol/H2O (50/50, v/v) instead of the sample in DPPH and ABTS.
All calibrations were performed using gallic acid as a reference compound in the ranges of 35–350 mg l−1, 2.3–18.4 mg l−1, 30–300 mg l−1, 14–140 mg l−1 and 12–120 mg l−1 for FC, FRAP, CUPRAC, DPPH and ABTS, respectively. The antioxidative capacities were calculated as % GAE (gallic acid equivalent) to enable a comparison. The measurements of one assay lasted three weeks on average. Ten replicate measurements at high, middle and low concentration levels were performed and the relative standard deviation (RSD) was calculated for each level to provide the various precisions.
000 cm−1 to 4000 cm−1. The absolute wavelength accuracy was ±2 cm−1 and the relative reproducibility was 0.2 cm−1. Every ground sample was measured in triplicate with 64 scans and the resolution was set to 8 cm−1. All measurements were performed at 22 °C within three days.
For spectral pretreatments and multivariate data analysis like the principal component analysis (PCA) with NIPALS (nonlinear iterative partial least squares) algorithm, partial least square regression (PLS) and all statistical analyses, The Unscrambler 10.2 chemometric software (Camo, Oslo, Norway) was used.
The spectra were transformed log
1/R and the sample average of the sample measurements (three for NIR and five for ATR-IR) was determined. The ATR-IR spectra were first normalized between one and zero then averaged and subsequently treated with 1st or 2nd derivative Savitzky–Golay.24 The NIR spectra were pretreated with 1st or 2nd derivative Savitzky–Golay,24 standard normal variate (SNV)25 and multiplicative scatter correction (MSC).26 Wavenumber regions with low information content or higher noise were removed. For all derivative calculations the optimal polynomial order and window size for smoothing were evaluated. For each assay a PLS regression model was calculated. To get the significant wavelength regions for every model, the uncertainty test function of The Unscrambler 10.2 software was executed. For both spectroscopic methods test-set and cross-validation models were performed. For the test-set validation, two thirds of the samples representing the reference data (calibration set) and the other one third that generates the test set (validation set) were chosen.
To evaluate each PLS calibration method, the quality parameters, standard error of prediction (SEP), standard error of cross-validation (SECV), R2 and the ratio performance to deviation (RPD), were computed and compared. The SEP is the standard deviation of differences between reference data and the results in the validation set.27 The squared coefficient of determination (R2) describes the squared quotient of the covariance of the reference data and the predicted data, divided by the product of their standard variations.28 The RPD is defined as the ratio of the standard deviation of the reference data to the SEP or SECV.28 The best PLS and the optimum number of PCs were selected according to the SECV/SEC or SEP/SEC ratio near to one, the highest R2 and the highest RPD value. SEC is the standard error of calibration.
Solvent A was water with 0.5% formic acid and B acetonitrile. A gradient program was executed in the following steps (min/% B): 0/5, 15/16, 30/25, 40/35, 41/95 and 47/95. The flow rate was 1 ml min−1 and the injection volume was 10 μl. The detection wavelength for rutin was 260 nm and for chlorogenic acid 320 nm. The mass spectrometry measurements were executed using ESI in negative-ion mode. The scan range was 200–1000 m/z and the scan speed was 830 u s−1. The interface voltage was 0.95 kV. The temperature of the desolvation line (DL) was 250 °C and that of the heat block was 200 °C. Drying gas flow was set at 15 l min−1 and nebulizing gas flow was 1.5 l min−1. The unfrozen extracts of samples 15 and 29 were centrifuged once again at 3500 rpm for 30 minutes at 22 °C and 200 μl were filled in a HPLC vial inlet. Both samples were measured by HPLC-DAD-MS.
The results of FC, FRAP and CUPRAC assays showed high Pearson correlations among each which are displayed in Table 1.
| FC | FRAP | CUPRAC | DPPH | ABTS | |
|---|---|---|---|---|---|
| FC | 1.00 | 0.94 | 0.95 | 0.85 | 0.88 |
| FRAP | 1.00 | 0.94 | 0.86 | 0.84 | |
| CUPRAC | 1.00 | 0.86 | 0.85 | ||
| DPPH | 1.00 | 0.84 | |||
| ABTS | 1.00 |
O stretching, non-conjugated), 1606 cm−1 (C
C stretching and C
O stretching, conjugated) and 1027 cm−1 (C–O, C–C stretch of carbohydrates and flavonoids).
For the PLS methods (test-set and cross-validation) of all assays the 2nd derivative Savitzky–Golay with a polynomial order of three and a window size of 23 (see Fig. 5) showed the best results. To compare the PLS models the cross-validation results were used because all samples are included in the test set. In the test-set validation the validation samples vary which complicates a direct comparison. The cross-validation gave the best results for the CUPRAC assay and its regression coefficient plot of the model is displayed in Fig. 3. Most prominent regression coefficients appeared at 1187 cm−1, 1123 cm−1 and 986 cm−1 and are assigned to aromatic CH in plane bending vibrations. The band at 814 cm−1 was assigned to C–O–C symmetric stretching vibrations of saturated heterocycles like sugars and to aromatic CH out of plane bending. The band at 1726 cm−1 was assigned to carbonyl stretching. Band assignments were realized according to Pretsch, Bühlmann and Badertscher.30
![]() | ||
| Fig. 4 CUPRAC predicted vs. reference plot and explained variance plot calculated with ATR-IR spectra applying cross-validation. | ||
![]() | ||
| Fig. 5 CUPRAC predicted vs. reference plot and explained variance plot calculated with ATR-IR spectra applying test-set validation. | ||
The quality parameters of the PLS models are listed in Table 3.
Example of the ATR-results with CUPRAC as the reference method: both validations were performed applying 2nd derivative Savitzky–Golay, polynomial order of three and a window size of 23 (Fig. 4 and Fig. 5).
Looking at the different assays executed in this study, the Folin–Ciocalteu (FC) assay showed the best precision and accuracy among the used wet chemical methods as shown in Table 2. Another advantage of this assay was its simple preparation procedure.
| Inter-day precision of different concentrations | Intra-day precision of different concentrations | Accuracy | |||||
|---|---|---|---|---|---|---|---|
| Low | Middle | High | Low | Middle | High | ||
| FC | 0.72 | 0.73 | 0.60 | 0.51 | 0.41 | 0.38 | 99.27 |
| FRAP | 1.48 | 1.70 | 1.25 | 0.66 | 1.15 | 0.66 | 97.40 |
| CUPRAC | 1.37 | 1.12 | 0.96 | 0.59 | 0.65 | 0.38 | 98.68 |
| DPPH | 2.53 | 1.46 | 1.11 | 1.56 | 0.94 | 0.82 | 86.61 |
| ABTS | 3.01 | 2.72 | 2.48 | 2.06 | 1.74 | 1.40 | 95.90 |
However, the results of the ATR-IR PLS models did not show the best RPD and R2 values for the FC assay. The various models resulted in the following order with decreasing RPD and R2 values: CUPRAC > ABTS > FC > DPPH > FRAP. All RPD values of the test-set validation were a bit lower than the corresponding values in cross-validation, because the cross-validation includes all samples for validation, compared with the test-set validation the results of which depend on the selected validation samples (see Table 3).
| R 2 | SECV | SEC | SECV/SEC | LV | RPD | RPDcorr | RPDPrec | RBE | Range | PrecMean | |
|---|---|---|---|---|---|---|---|---|---|---|---|
| ATR-IR cross-validation | |||||||||||
| FC | 0.84 | 0.15 | 0.11 | 1.38 | 6 | 2.54 | 2.62 | 10.16 | 3 | 4.1–6.0 | 0.68 |
| FRAP | 0.82 | 0.09 | 0.07 | 1.23 | 5 | 2.29 | 2.47 | 6.11 | 7 | 1.6–2.6 | 1.48 |
| CUPRAC | 0.85 | 0.13 | 0.10 | 1.25 | 5 | 2.56 | 2.68 | 8.72 | 4 | 2.8–4.6 | 1.15 |
| DPPH | 0.83 | 0.07 | 0.05 | 1.27 | 5 | 2.43 | 2.56 | 7.74 | 5 | 0.9–1.8 | 1.70 |
| ABTS | 0.84 | 0.04 | 0.04 | 1.24 | 5 | 2.47 | 3.85 | 3.23 | 36 | 0.9–1.5 | 2.73 |
| R 2 | SEP | SEC | SEP/SEC | LV | RPD | RPDcorr | RPDPrec | RBE | Range | PrecMean | |
|---|---|---|---|---|---|---|---|---|---|---|---|
| ATR-IR test-set validation | |||||||||||
| FC | 0.84 | 0.17 | 0.14 | 1.25 | 6 | 2.27 | 2.32 | 10.25 | 2 | 4.1–6.1 | 0.68 |
| FRAP | 0.80 | 0.09 | 0.08 | 1.16 | 5 | 2.35 | 2.53 | 6.27 | 7 | 1.6–2.6 | 1.48 |
| CUPRAC | 0.81 | 0.15 | 0.13 | 1.14 | 5 | 2.12 | 2.20 | 8.29 | 3 | 2.8–4.6 | 1.15 |
| DPPH | 0.84 | 0.07 | 0.06 | 1.17 | 5 | 2.39 | 2.51 | 8.04 | 5 | 0.9–1.9 | 1.70 |
| ABTS | 0.80 | 0.05 | 0.04 | 1.04 | 5 | 2.24 | 3.42 | 2.97 | 34 | 0.9–1.5 | 2.73 |
The best results for the ATR-IR models were achieved with the 2nd derivative Savitzky–Golay. A polynomial order of three and a window size of 23 were used to perform the PLS models. Performing the FC and CUPRAC PLS models no sample was removed from the sample set. Compared to the other PLS models two, three and four samples were removed from the sample set of FRAP, DPPH and ABTS, respectively. Except for the PLS model of the FRAP assay, all developed modes showed a RPD value higher than 2.5 in cross-validation. In test-set validation only ABTS reached a good RDP value about 2.66. According to Williams et al. RPD values over 2.5 are acceptable.28 In both validation methods a R2 value higher than 0.8 was achieved. Five LVs were needed for almost all PLS models except for FC (LV = 6). As displayed in Fig. 3 wavelengths between 800 cm−1 and 1800 cm−1 were important to perform the models. The bands over 2800 cm−1 (O–H stretching of water, carbohydrates and flavonoids and C–H symmetric and asymmetric stretching) were not important to perform the models.
The NIR models performed really poor and the cross-validation results are depicted in Table 4. The number of LVs needed for all PLS models varied between five and seven. Only for the FC PLS model acceptable values for R2 and RPD were obtained, but the explained variance plot showed that the calculated model was of low quality because the first two LVs resulted in explained variances lower than 8%. A reason for the worse prediction power of NIR could be that the antioxidant power is carried mainly by two single antioxidative compounds (rutin and chlorogenic acid) in the plant material which belongs to two different structural classes. Former studies revealed a higher predictive power for the determination of single compounds with ATR-IR in plant matrices,33–36 which is advantageous in this application. In ATR-IR, fundamental vibrations are present exhibiting less overlap than the combination and overtone bands of NIR. This makes NIR more sensitive to matrix variations,35 which are usually severe in natural products. Also in earlier publications ATR-IR spectroscopy revealed better results than NIR spectroscopy for the quantification of single compounds in plant materials.33,34 Comparing the standard errors obtained from the studies of Navarro et al.,33 Schönbichler et al.,34 Krähmer et al.35 and Baranska et al.,36 which dealt with quantification of single compounds in plant materials, the NIR exhibited values with a factor 1.1–2.7 times higher than the ATR-IR in the cross-validation results. Even if there is evidence pointing to the advantages of ATR-IR in single compound prediction, it has not been reported until now that ATR-IR gives good results for an application and the NIR fails completely in the same application. The exact reason for this difference in performance of the two spectroscopic methods in this study could not be found. On top of that, our previous work comparing ATR-IR and NIR determining the antioxidative capacity of primulae flos powder revealed better results for NIR.14 The good results for NIR in the case of primulae flos can be explained by the contribution of mainly one class of substances to the antioxidative capacity. Here the broad bands of NIR seem to have an advantage for predicting a sum parameter of those substances. In the case of Sambuci flos mainly two substances contribute to the antioxidative capacity, and they belong to different structural classes. The authors assume that for NIR either the interference of chlorogenic acid with rutin or the interference of other matrix compounds causes problems in the PLS regression models. Probably ATR-IR can distinguish better between those two substances and matrix components.
| R 2 | SECV/SEC | LV | RPD | RPDcorr | RBE | Cumulative explained variance for the first LVs in % | |||
|---|---|---|---|---|---|---|---|---|---|
| 1st LV | 2nd LV | 3rd LV | |||||||
| FC | 0.88 | 1.30 | 7 | 2.93 | 3.06 | 4 | 2 | 5 | 60 |
| FRAP | 0.77 | 1.25 | 7 | 2.10 | 2.24 | 6 | 3 | 5 | 17 |
| CUPRAC | 0.76 | 1.33 | 7 | 2.09 | 2.15 | 3 | 2 | 7 | 30 |
| DPPH | 0.78 | 1.33 | 6 | 2.02 | 2.09 | 3 | 5 | 2 | 25 |
| ABTS | 0.79 | 1.20 | 5 | 1.92 | 2.34 | 18 | 5 | 15 | 40 |
Footnotes |
| † Electronic supplementary information (ESI) available. See DOI: 10.1039/c5ay01314c |
| ‡ These authors contributed equally. |
| This journal is © The Royal Society of Chemistry 2016 |