Subscriber access provided by NEW MEXICO STATE UNIV
Article
Head to head comparison of uHPLC-DAD versus qNMR for the quantitative analysis of the silymarin complex in Silybum marianum fruit extracts Antigoni Cheilari, Sonja Sturm, Daniel Intelmann, Christoph Seger, and Hermann Stuppner J. Agric. Food Chem., Just Accepted Manuscript • DOI: 10.1021/acs.jafc.5b05494 • Publication Date (Web): 23 Jan 2016 Downloaded from http://pubs.acs.org on January 25, 2016
Just Accepted “Just Accepted” manuscripts have been peer-reviewed and accepted for publication. They are posted online prior to technical editing, formatting for publication and author proofing. The American Chemical Society provides “Just Accepted” as a free service to the research community to expedite the dissemination of scientific material as soon as possible after acceptance. “Just Accepted” manuscripts appear in full in PDF format accompanied by an HTML abstract. “Just Accepted” manuscripts have been fully peer reviewed, but should not be considered the official version of record. They are accessible to all readers and citable by the Digital Object Identifier (DOI®). “Just Accepted” is an optional service offered to authors. Therefore, the “Just Accepted” Web site may not include all articles that will be published in the journal. After a manuscript is technically edited and formatted, it will be removed from the “Just Accepted” Web site and published as an ASAP article. Note that technical editing may introduce minor changes to the manuscript text and/or graphics which could affect content, and all legal disclaimers and ethical guidelines that apply to the journal pertain. ACS cannot be held responsible for errors or consequences arising from the use of information contained in these “Just Accepted” manuscripts.
Journal of Agricultural and Food Chemistry is published by the American Chemical Society. 1155 Sixteenth Street N.W., Washington, DC 20036 Published by American Chemical Society. Copyright © American Chemical Society. However, no copyright claim is made to original U.S. Government works, or works produced by employees of any Commonwealth realm Crown government in the course of their duties.
Page 1 of 28
Journal of Agricultural and Food Chemistry
Head to Head Comparison of uHPLC-DAD versus qNMR for the Quantitative Analysis of the Silymarin Complex in Silybum marianum Fruit Extracts
Antigoni Cheilari1, Sonja Sturm1*, Daniel Intelmann2, Christoph Seger1, Hermann Stuppner1
Affiliation Institute of Pharmacy, Department of Pharmacognosy, CCB – Centrum of Chemistry and Biomedicine, University of Innsbruck, Innrain 80/82, 6020 Innsbruck, Austria Bionorica research GmbH, Mitterweg 24, 6020 Innsbruck, Austria
* Correspondence: Dr. Sonja Sturm, Institute of Pharmacy, Department of Pharmacognosy, CCB – Centrum of Chemistry and Biomedicine, University of Innsbruck, Innrain 80/82, 6020 Innsbruck, Austria Email:
[email protected], Tel: 0043 512 507 58408, Fax: 0043 512 507 58499
ACS Paragon Plus Environment
Journal of Agricultural and Food Chemistry
1
Abstract
2
Quantitative NMR spectroscopy (qNMR) is known as excellent alternative to chromatography based
3
mixture analysis. NMR spectroscopy is a non-destructive method, needs only limited sample
4
preparation and can be readily automated. A head-to-head comparison of qNMR to an uHPLC-DAD
5
based quantitative analysis of six flavonolignan congeners (silychristin, silydianin, silybin A, silybin B,
6
isosilybin A, isosilybin B) of the Silybum marianum silymarin-complex is presented. Both assays
7
showed similar performance characteristics (linear range, accuracy, precision, LOQs) with analysis
8
times below 30 min / sample. The assays were applied to industrial S. marianum extracts (AC
9
samples) and to extracts locally prepared from S. marianum fruits (PL samples). Assay comparison by
10
Bland-Altman plots (relative method bias AC samples = -0.1%, 2S range ± 5.1%, relative method bias
11
PL samples = -0.3%, 2S range ± 7.8%) and Passing-Bablok regression analysis (slope and intercept for
12
AC and PL samples not significantly different from 1.00 and 0.00, respectively; Spearman’s coefficient
13
of rank correlation >0.99) did show, that qNMR and uHPLC-DAD can be used interchangeably to
14
quantitate flavonolignans in the silymarin complex.
15 16 17 18 19 20
Key words
21
Silybum marianum, Asteraceae, silymarin, flavonolignans, uHPLC-DAD, qNMR, quantitation
22
ACS Paragon Plus Environment
Page 2 of 28
Page 3 of 28
23
Journal of Agricultural and Food Chemistry
Introduction
24 25
NMR spectroscopy based analysis of mixtures is an appealing alternative to chromatographic analysis
26
– both in the qualitative (metabolic profiling)1,2 and the quantitative (qNMR) context.3 The analytical
27
concept of qNMR, usually performed using 1H NMR resonances, can be considered an untargeted
28
approach (any analyte combination can be quantified if appropriate undisturbed signals can be
29
identified) and offers some advantages over chromatographic analyses.4 It is non-destructive, sample
30
preparation is of limited complexity and generally compatible to chromatography, can be automated,
31
and provides – without the need of chromatographic separation or additional analytical devices
32
(hyphenation) – both qualitative and quantitative analyte information. Aside quantitative
33
information derived from carefully selected single NMR signals the complete 1H NMR spectrum of
34
the targeted compound allows a deep insight into its chemical structure. A set of 1H NMR resonances
35
can be utilized to confirm the identity of a targeted entity in a mixture – e.g. a methine 1H-NMR
36
resonance is used for quantitation and the methyl group 1H-NMR signal pattern is used for
37
qualification. In addition, an unbiased and holistic view of the total sample composition (alike the
38
total ion current in LC-MS) is obtained and allows checking for sample adulterations, unusual signal
39
pattern, solvent residuals and many other features.4,5,6
40 41
From a practical point of view, quantitative and qualitative phytochemical analysis in research,
42
production, and quality control is always constrained by limited availability of resources including
43
sample material, access to instrumentation, and investigation time. Hence replacing technologies
44
with low sample throughput, high personnel demand, and extensive sample workup requirements by
45
alternative approaches is a central demand in such processes. Due to the advantages of NMR
46
spectroscopy over chromatographic assays as listed above qNMR can be understood as such an
47
alternative. However, although both qNMR and HPLC-UV are well established and frequently applied
48
technologies in this context, head-to-head comparisons of quantitative assays have hardly been
49
performed.
50 51
The group around Verpoorte7 did report on artemisinin quantitation in eight samples by qNMR,
52
HPLC-ELSD, HPLC-MS, and TLC regarding qNMR as gold standard. They reported qNMR easy to use
53
and found all methods quantitatively comparable. Pieri et al. did compare NMR derived quantitative
54
results to findings from chromatographic assays in comprehensive Thymus vulgaris metabolite
55
profiling,8 in a targeted analysis of cynaropicrin from artichoke leaf extracts,9 and in an investigation
56
devoted to quantitate steviol glycosides from Stevia rebaudiana.10 Pauli and co-workers did devote a
57
recent publication to the analysis of green tea catechins by qNMR and LC-MS/MS.11 These authors
ACS Paragon Plus Environment
Journal of Agricultural and Food Chemistry
58
used the qNMR/HiFSA (1H-NMR iterative full spin analysis) approach4 for NMR signal selection and
59
validated their approach against the NIST reference material SRM3255. Seven analytes were
60
quantified in two extracts and the SRM with LC-MS/MS and qNMR results showing high congruency.
61 62
Milk thistle (Silybum marianum (L) GAERTN.) preparations have been used for centuries to treat a
63
variety of illnesses, particularly liver maladies.12 The effective biological activities of the plant are
64
credited to silymarin, an isomeric mixture of at least six flavonolignans (Figure 1), namely silychristin
65
(1), silydianin (2), silybin A (3), silybin B (4), isosilybin A (5) and isosilybin B (6), and additionally
66
taxifolin, a flavonoid. The fruits contain up to 9% of those compounds, which mainly accumulate in
67
their external cover.13,14 The hepatoprotective activity of silymarin can be explained by the
68
antioxidant properties due to the phenolic nature of flavonolignans, stimulating liver cells
69
regeneration and cell membrane stabilization to prevent hepatotoxic agents from entering
70
hepatocytes.15 Numerous studies have demonstrated the efficacy of silymarin as a hepatoprotective,
71
but also as an anti-inflammatory, antioxidant, antifibrotic and cancer chemopreventive agent.15-18
72 73
In herbal supplement preparations the amount of silymarin in S. marianum extracts varies from 55-
74
85% with the exact composition still not completely unravelled.19 Up to now a number of HPLC-UV
75
methods have been published to support quality control measures; some however with either rather
76
long analysis times (30 min – 100 min)20-26 or limited validation data presentation.21,26 In addition LC-
77
MS/MS instrumentation was used to characterize S. marianum extract constituents27-29 or to follow
78
them in bio-analytical relevant matrices like plasma.30
79 80
With these reports at hands it was envisioned to perform a head-to-head comparison of an HPLC
81
platform with qNMR. Since there is still the need to prove the equivalence of qNMR to the
82
phytoanalytical “gold-standard” methods due to its inherent orthogonality to HPLC, the
83
establishment of a thoroughly validated HPLC method was considered as prerequisite in such
84
approach and not as an “add-on” to prove the validity on NMR derived quantitative data.
85
Consequently, it was decided to establish a state-of-the-art HPLC assay utilizing a modern stationary
86
phase. Due to the chosen research object (the S. marianum derived silymarin complex) MS/MS
87
detection was deemed unnecessary and UV/VIS detection was seen sufficient.
88 89
Materials and Methods
90
Chemicals and Reference Compounds
91
Methanol and formic acid were of HPLC grade and purchased from Merck (Darmstadt, Germany) and
92
Sigma-Aldrich (Vienna, Austria) respectively. Methanol-d4 (99.90% containing 0.03% TMS) was from
ACS Paragon Plus Environment
Page 4 of 28
Page 5 of 28
Journal of Agricultural and Food Chemistry
93
Euriso-Top (Saint-Aubin, France). Water for HPLC was produced on site using an Arium611UV
94
(Sartorius Stedim Biotech, Göttingen, Germany) water purification system. Reference compounds 1 -
95
6 were purchased from PhytoLab (Vestenbergsgreuth, Germany) and anthracene from Sigma Aldrich
96
(Vienna, Austria). Purity of all standard compounds was ≥98.5% (determined by HPLC and NMR).
97 98
Plant material and sample preparation
99
Bionorica Research GmbH (Innsbruck, Austria) provided dry S. marianum fruit extracts (extraction
100
solvent acetone, AC) and dry S. marianum fruits (PL). Voucher specimens of these materials were
101
deposited at the Institute of Pharmacy/Pharmacognosy. AC samples were weighted (10.0 ± 0.1 mg)
102
into 1.5 ml polyethylene microcentrifuge tubes (Eppendorf, Hamburg, Germany). 1 mL methanol-d4
103
containing anthracene (0.8 mg/mL) as internal standard (IS) was added and after vigorous mixing (1
104
min) on a Vortex mixer (VWR, Vienna, Austria) samples were incubated at 24°C under mild agitation
105
overnight to complete solvatization. After centrifugation (10 min at 24 oC and 20817 g) 650 μL of the
106
supernatants was transferred to NMR tubes, 300 µl to HPLC vials. Each sample was prepared in
107
triplicate. All samples were stored at 4°C until analysis. PL samples were ground and 100.0 ± 0.1 mg
108
were extracted five times with 1 mL methanol in the Ultra Sonic bath for 15 minutes. After each
109
extraction cycle, samples were centrifuged as given for the AC samples. Centrifugation supernatants
110
were combined, evaporated, and redissolved in 2 mL of methanol-d4 fortified with the IS (0.8
111
mg/mL). After centrifugation (5 min at 24 oC and 20817 g) 650 μL of the supernatants was transferred
112
to NMR tubes, 300 µl to HPLC vials. Each sample was prepared in triplicate. All samples were stored
113
at 4°C until analysis.
114 115
NMR experimental parameters
116
NMR spectra were acquired at 300 K with a Bruker Avance II 600 spectrometer (Bruker Biospin,
117
Rheinstetten, Germany) equipped with a 5 mm Prodigy Cryo-probehead with Z-gradient.
118
Experiments were performed in automation mode, using a Bruker BACS-60 sample changer operated
119
by Bruker IconNMR. Data acquisition and processing were done with Bruker Topspin 3.1. 1H-NMR
120
spectra were acquired using the Bruker zg0pr pulse program with the following settings: relaxation
121
delay (d1) = 15 s; flip angle = 45o; acquisition time = 2.99 s; FID (free induction decay) data points = 64
122
K; spectral width = 20 ppm; number of scans = 96. For experiments using presaturation, the
123
transmitter offset was manually set in order to achieve optimal suppression of the residual water
124
signal. In all cases the acquired FIDs were Fourier transformed to yield spectra with 128 K data points
125
(zero filling). Manual phase correction and automatic polynomial baseline correction were always
126
used. Chemical shift values were referenced to the calibration standard (TMS) signal. Inversion-
127
recovery experiments were performed using the Bruker t1ir pulse program, with standard acquisition
ACS Paragon Plus Environment
Journal of Agricultural and Food Chemistry
128
parameters, setting the τ (tau) parameter to 0.01, 0.05, 0.1, 0.25, 0.5, 1, 2, 4, 8, and 15 seconds. T1
129
(longitudinal relaxation) values were calculated using the T1 relaxation routine (Topspin 3.1).
130 131
HPLC-DAD and LC-MS conditions
132
HPLC analyses were performed with a HP1260 system equipped with binary pump, auto sampler,
133
column thermostat, and photodiode array detector (Agilent, Waldbronn, Germany). An Agilent
134
Poroshell 120 EC-C18 column (3.0 x 50 mm, 2.7 μm particle size) was used as stationary phase and
135
water (A) / methanol (B), both containing 0.1% formic acid, as mobile phase solvents. Separation of
136
the analytes (1 µl injection volume) was achieved by applying a solvent gradient of 0 min 80% A; 3
137
min 70% A; 6 min 60% A; 8 min 60% A; 12 min 55% A; 13 min 2% A, 17 min 2% A delivered at a flow
138
rate of 0.6 ml/min. The gradient program was followed by a re-equilibration period of 10 min at 80%
139
A. The column oven temperature was set to 50°C, the detection wavelength to 286 nm. HPLC-MS/MS
140
experiments for HPLC-DAD peak assignment confirmation were performed on an amaZon SL ion-trap
141
mass spectrometer (Bruker-Daltonics, Bremen, Germany) hyphenated to the Agilent HP1260 system
142
operated under the given chromatography conditions. MS parameters: ESI alternating mode; spray
143
voltage = 4.5 kV; dry gas = 8 L/min at 250°C; nebulizer = 30 psi; full scan mode = m/z 100-1500.
144
MS/MS experiments were performed in the automatic mode.
145 146
Calibration and Validation
147
Both uHPLC and NMR methods were validated according to the ICH guidelines ‘’Validation of
148
Analytical Procedures: Text and Methodology Q2 (R1)’’31 for linearity, limit of detection and
149
quantitation, peak purity, accuracy, precision and repeatability. Standard stock solutions were
150
prepared by dissolving compounds 1 – 6 in methanol-d4 containing the internal standard (IS) (0.8
151
mg/mL). From these stock solutions nine calibrator levels (2000 μg/mL, 1000 μg/mL, 700 μg/mL, 500
152
μg/mL, 310 μg/mL, 200 μg/mL, 45 μg/mL, 25 μg/mL, 10 μg/mL) were prepared. Calibration curves
153
were obtained by plotting the peak areas of analytes 1 – 6 versus the concentrations of 1 – 6 in
154
uHPLC and by plotting the respective ratios of areas and concentrations of 1 – 6 and IS in NMR. The
155
regression parameters intercept, slope and correlation coefficient (R2) were calculated by linear
156
regression analysis. The limit of detection (LOD) and the limit of quantitation (LOQ) for each analyte
157
were calculated from the regression models including only the lowest three dilution levels. From the
158
obtained regression equation data the LOD was calculated as three times the standard deviation (SD)
159
of the y-intercept divided by the slope, whereas the LOQ was calculated as ten times the SD of the y-
160
intercept divided by the slope. Accuracy was determined by recovery experiments spiking a dry
161
extract (AC 01) and a fruit batch (PL 01) with known amounts of selected analytes prior to sample
162
workup. All samples were prepared in triplicate. By analysing the variability (triplicate analysis) of all
ACS Paragon Plus Environment
Page 6 of 28
Page 7 of 28
Journal of Agricultural and Food Chemistry
163
recovery samples on three consecutive days intra-day and inter-day precision of the assays was
164
assessed. Quantitative results have not been corrected for residual traces of moisture or solvents.
165 166
Statistics
167
Data analysis for validation and calculation of analyte concentrations was done within Microsoft
168
Excel 2011 (Redmond, WA, USA) spread sheets. Statistical data analysis (comparative Bland-Altman
169
plots, Passing Bablok regression, rank correlation analysis) was done with MedCalc for Windows,
170
version 15.8 (MedCalc Software, Ostend, Belgium)
171 172
Results and Discussion
173
Sample preparation for NMR and HPLC
174
Two sample types were used in this investigation. Beside industry based plant extracts (extraction
175
solvent acetone, AC) typically used as formulation basis in commercial milk thistle medications (>50%
176
silymarin), laboratory prepared fruit extracts (PL) were additionally investigated to introduce broader
177
sample variability into the investigation. Due to the distinctively different matrix properties, tailored
178
sample preparation protocols for AC and PL had to be established. In a series of preliminary 1H-NMR
179
experiments several deuterium labeled solvents including methanol-d4, aceton-d6, DMSO-d6,
180
acetonitrile-d3, and deuterium oxide as well as mixtures thereof were tested as potential extraction
181
solvents for both AC and PL samples and as solvents for pure reference materials. The obtained
182
spectra showed that methanol-d4 was the solvent of choice enabling sufficient solubility of reference
183
compounds and extracts. A good separation of NMR key resonances and a clean baseline in the
184
relevant spectral regions was obtained. Anthracene was chosen as internal standard for qNMR due
185
to its solubility and stability in methanol, its low volatility and its aromatic nature resulting neither in
186
1
187
were directly re-dissolved in methanol-d4, the preparation of PL extracts was optimized to get a
188
maximum of secondary metabolites into solution. Different extraction solvents (ethanol, methanol,
189
chloroform, hexane and mixtures of them with water) and extraction procedures (sonication,
190
Eppendorf shaker) were evaluated. Sonication supported quantitative extraction of the ground plant
191
material with methanol as solvent was found to be most efficient.
H NMR signal overlap with signals of 1 - 6 nor in HPLC co-elution with 1 - 6. Whereas AC extracts
192 193
Validation of quantitative NMR analysis
194
The application of quantitative 1H NMR (qNMR) requires that at least one non-overlapping signal for
195
each molecule to be quantified is available for integration.5 The research efforts were concentrated
196
to the principal components of silymarin – silybin A (1), silybin B (2), isosilybin A (3), isosilybin B (4),
197
silychristin (5), silydianin (6) (Figure 1). All compounds are isobaric (MW 482) with compounds 3-6
ACS Paragon Plus Environment
Journal of Agricultural and Food Chemistry
Page 8 of 28
198
even representing isobaric pairs of diastereoisomers (3/4 and 5/6) with very similar 1H and 13C NMR
199
spectra. Silybins (3/4) and isosilybins (5/6) have the same trans-conformation in the C ring of the
200
flavonoid skeleton and opposite conformation and substitution pattern in the dioxan scaffold with
201
1
202
Pauli and co-workers have shown that applying quantum mechanical driven HiFSA even almost
203
overlapping proton signals can be used for analyte quantitation – but this approach was not
204
extended to S. marianum extracts.32 Having the complexity of an extract compared to a reference
205
compound standard mixture in mind and since the HiFSA approach depends on software only
206
commercially available; it was decided not to take up this approach but to quantitate 3,4 and 5,6 as
207
compound pairs.
H-NMR shift differences (Δδ) not exceeding 0.02 ppm at 600 MHz 1H-NMR resonance frequency.
208 209
Careful analysis of the 1H-NMR spectra of the single compounds and of mixtures in the extraction
210
solvent led to a series of methine resonances suitable for analyte quantitation in mixtures (Table 1).
211
For analyte pairs 3,4 and 5,6 aromatic doublets (JHH ≈ 2.1 Hz) between 7.08 ppm and 7.15 ppm have
212
been chosen, for 1 an OCHCH doublet at 5.55 ppm (JHH ≈ 6.3 Hz) and for 2 a broadened methine
213
singlet at 2.89 ppm (Figure 2). The peak purity of these signals was assured by the aid of homo- and
214
heteronuclear 2D NMR experiments (COSY, HSQC, HMBC), which enabled unambiguous assignment
215
of these resonances in the proton scaffold of compounds 1 – 6. The patterns of these key signals
216
were identical in pure compounds and matrix containing sample materials, an influence of
217
satellite signals in the signal pairs 3,4 and 5,6 can be excluded, since the distance between these
218
resonances does not match the anticipated 1J(13C,1H) ≈ 120 – 150 Hz.
13
C
219 220
A crucial parameter in qNMR is the relaxation delay, which enables equilibrium of the magnetization
221
between pulses when accumulating co-added FIDs. It depends on the relaxation properties of the
222
various nuclei in the mixture and is determined by inversion-recovery T1 experiments for all
223
resonances of interest. The highest T1-value was calculated to be 4.7 s for the selected signal of
224
anthracene; all Silybum constituents showed shorter T1-values between 1.54 s and 2.22 s (Table 1).
225
In order to obtain a full relaxation of all nuclei of interest a relaxation delay of 5*T1 = 23.5 s was
226
calculated for a flip angle of 90o. After optimization experiments a flip angle of 45o and a relaxation
227
delay of 15 s (acquisition time 2.99 s) were chosen and resulted in a 99.9% recovery of the internal
228
standard signal with the longest T1. The total analysis time including sample transfer, sample
229
temperature equilibration, and sample shimming was about 35 min with a net acquisition time of
230
about 29 min.
231
ACS Paragon Plus Environment
Page 9 of 28
Journal of Agricultural and Food Chemistry
232
NMR method validation covered assay linearity, limit of detection and quantitation, accuracy,
233
precision, and repeatability. Based on the calibration models, linearity of the method was confirmed
234
for the concentration range of 10 µg/mL to 2000 μg/mL for compound 2, 3, and 6 and 25 µg/mL to
235
2000 μg/mL for 1, 4, and 5 (Table 2). In these ranges the correlation coefficient of the calibration
236
function was better than 0.9986 for all analytes. Limits of detection (LOD) and limits of quantitation
237
(LOQ) determined by an independent approach were found to range from 0.5 µg/mL to 3.4 µg/mL
238
and 1.6 µg/mL to 10.4 μg/mL, respectively. The discrepancy between found LOQs and the lower limit
239
of the linear range does prove, that for analyte concentration ranges between these figures of merit
240
a separate calibration model must be implemented. Consequently in this investigation the LLOQs
241
(lower limits of quantitation) were set to the lowest calibrator concentration in the linear range. The
242
signal to noise ratios for these concentrations did exceed 20 units.
243 244
Accuracy was determined in recovery experiments with three different concentrations of standard
245
compounds spiked to AC samples or added to the plant material (PL) prior to extraction. All results
246
were within the range of 90.7% – 106.4% of the amount fortified with no observable bias between
247
the methods of the materials (Table 4). Hence for both assays the analytical analyte recovery is
248
quantitative and the LOQs established for the solvent-based calibrators can be used for the AC and
249
PL samples. Precision of the assay was determined by preparation and analysis of fortification
250
samples on three days. Intra-day precision was better than 7.7% over all analytes and materials;
251
inter-day precision was found to be better than 8.6% (Table 5). There was no precision difference
252
between the materials AC and PL, result scatter for analyte 2, present in rather low concentrations
253
(Table 6 and Table 7), was slightly increased compared to the other analytes.
254 255
Validation of quantitative uHPLC analysis
256
To allow head-to-head comparison of the silymarin complex qNMR assay with a more conventional
257
analytical technique an uHPLC-DAD assay was established. As to be expected, the quite similar
258
physicochemical properties of the diastereomeric analyte pairs 3/4 and 5/6 were the major hurdle in
259
the development of a sufficient selective method. In a first round of experiments several silica based
260
reverse phase materials (RP-12, RP-18, phenyl-, phenyl-hexyl) as stationary phases and different
261
mobile phase conditions were evaluated. The use of acetonitrile in the mobile phase was clearly not
262
appropriate, but methanol alone was also not giving satisfactory separation. Acidic mobile phase
263
additives were improving peak shape; optimal results were achieved with water to methanol
264
gradient elution using formic acid as modifier. Conventional stationary phase materials as the
265
endcapped RP-18 column (Zorbax Eclipse) were found insufficient. Even under optimized and lengthy
266
chromatographic conditions signal overlap at silydianin (2) led to a remarkable method bias with 2-3-
ACS Paragon Plus Environment
Journal of Agricultural and Food Chemistry
267
fold concentrations in the chromatographic assay compared to qNMR (data not shown). Taking into
268
account the selectivity of the chosen wavelength (286 nm) and the possibility of differences in
269
UV/VIS detection response factors (molar extinction coefficients) it is likely that an unaddressed
270
silydianin isomer or silydianin degradation product33 was co-eluting with silydianin (2) in this method.
271
Only replacing the fully porous Zorbax column by a superficially porous particle stationary phase
272
(Poroshell C18) led to baseline separation of all analytes including silydianin (2). To prove satisfactory
273
analyte separation in chromatography in addition to spiking experiments and comparison of
274
retention times and UV-spectra to reference materials, HPLC-MS/MS analysis was utilized to confirm
275
the peak assignment of 2 (tR = 6.40 min) as silydianin (Figure 3) since this analyte yields a distinctive
276
fragmentation pattern in MS/MS analysis (m/z = 151 and 169). The peak eluting shortly after 2 (tR =
277
6.60 min) showed an identical mass spectrum (negative ESI base peak with m/z = 481) and a similar
278
fragmentation pattern to compound 1 with m/z = 355 as most abundant fragment ion. Thus, this
279
peak could be tentatively assigned to silychristin B, a diastereomer of silychristin (1).29 The analyte
280
peak eluting between 1 and 2 (tR = 6.20 min) showed a significant different behaviour to compound 1
281
in mass spectrometry with a negative ESI base peak at m/z = 383 and – in contrast to the other
282
analytes – a favoured ionisation in the positive ESI modus (m/z = 385). Hence this analyte is most
283
likely not a flavonolignan derivative. The HPLC-DAD assay was validated according to the ICH-
284
guidelines.31 Linearity was confirmed from 10 -1000 µg/mL for 1 and 4 and from 25-2000 µg/mL for
285
2, 3, 5 and 6, correlation coefficients were higher than 0.999 in all cases (Table 3). LODs ranged from
286
1.0 µg/mL to 6.0 µg/mL and LOQs from 2.0 µg/mL to 19.0 µg/mL respectively (S/N > 20). The results
287
of accuracy experiments were comparable to qNMR data with recovery rates from 91.8 – 110.0 %
288
(Table 4). Intra-day variations less than 7.4% and within days variability did not exceed 5.9% (Table
289
5).
290 291
Sample analysis
292
The validated uHPLC-DAD and qNMR methods were subsequently applied for the quantitation of
293
nine batches of AC extracts and ten batches of PL extracts (Tables 6 and 7). AC 01 – 09 showed
294
relatively small variations in the contents of compounds 1 – 6, with AC 02 showing a slightly higher
295
amount of silybins (3,4) and AC 07 showing a higher content in isosilybin derivatives (5,6) and
296
silychristin (1). The ratio of 3 to 4 was relatively constant in AC 01 – 09 and was determined as 1 : 1.7
297
(± 3.3%) by HPLC-DAD, the ratio of 6 to 5 was 1 : 2.8 (± 0.5%). In all cases the total content of
298
flavonolignans was 59.8 ± 2.4% of the total extract. The content of all other flavonolignans, except 5
299
and 6 was significantly decreased. PL extracts derived from plant raw material showed distinctive
300
higher concentration variations and – as to be expected – overall lower relative amounts of 1 – 6
301
(silymarin complex content 9.4 ± 2.4%). Whereas the ratio of 3 to 4 in PL 01 – 10 was with 1 : 1.6 (±
ACS Paragon Plus Environment
Page 10 of 28
Page 11 of 28
Journal of Agricultural and Food Chemistry
302
12.1%) comparable to the AC materials, the ratio of 6 to 5 was found to be significantly different with
303
1 : 1.4 (± 4.3%).
304 305
The comparison of uHPLC and qNMR data by statistical means did prove that both techniques are
306
equivalent regarding the quantitative assessment of S. marianum flavonolignan congeners present in
307
the “silymarin” compound complex. Passing-Bablok correlation analysis of qNMR (Table 6, analytes 1,
308
2, 3/4, 5/6) and uHPLC (Table 7, analytes 1, 2 and sums of 3+4 and 5+6) derived analyte
309
concentrations in AC and PL samples (separate analysis due to different concentration units) did
310
result in ideal correlation equations (Figure 4) with intercepts statistically not different from zero (AC:
311
intercept = 2.91, 95% CI = -1.02 to 18.66; PL: intercept = -4.61, 95% CI = -12.75 to 1.35) and slopes
312
statistically not different from one (AC: slope = 0.99, 95% CI = 0.98 to 1.01; PL: slope = 1.02, 95% CI =
313
1.00 to 1.07). Calculating the rank correlation coefficient (Speaman’s rho) confirmed the visually
314
impressive correlation with rho > 0.99 for both analyses (Figure 4). Relative Bland-Altman plot based
315
method comparison analysis (Figure 4) did confirm the results of the Passing-Bablok correlation
316
analysis and further unveiled that neither AC nor PL samples analysis did show any significant
317
concentration dependent bias. This does prove, that minor matrix constituents neighbouring the
318
NMR signals of compounds 1 and 5+6 (Figure 2) do not significantly perturbate the qNMR based
319
analyte quantification. The 2S confidence intervals of the AC and PL Bland-Altman plots (2S = 5.1%
320
and 7.8% respectively) were in good agreement with the inter-day RSD data gathered in the
321
validation process of the assays (Table 5) and the RSD data obtained from repeated analysis of the
322
samples (Table 6 and Table 7). No signs of additional sources of imprecision (e.g. bias related effects
323
from matrix interferences) were detected. If however, a compound-by-compound quantification is
324
desired, the presented uHPLC assay must be seen superior to the qNMR method due to partial signal
325
overlap in the latter one.
326 327
This presentation aims to prove, that qNMR and uHPLC analysis, although in principal orthogonal
328
instrumental analysis technologies, can be interchangeably used to quantitate secondary metabolites
329
from raw material (fruit extracts = PL) and enriched industry based extracts (acetone dry extracts =
330
AC). In this particular case study sample preparation protocols were identical for qNMR and uHPLC,
331
extracts solved in deuterated NMR solvents can be easily applied to chromatographic assays. Key
332
figures of merit of the assay validation – linear range, LOQ / LOD, repeatability were found to be
333
more the less identical for both assays – this is remarkably different from a recent report, were LC-
334
MS/MS analyses showed a distinctively higher RSD than qNMR.11 Both assays were applied to
335
nineteen sample batches – nine AC samples and ten PL samples. Identical quantitative results were
336
obtained, no statistical differences were found between the methods. From an economic point of
ACS Paragon Plus Environment
Journal of Agricultural and Food Chemistry
337
view, both approaches are comparable, too when used in routine analysis. Whilst NMR equipment is
338
usually more expensive than chromatography instrumentation, costs per sample are lower in NMR
339
due to reduced solvent use and less need for consumables. Also measurement times were
340
comparable, although not been fully optimized both techniques are capable of running several dozen
341
samples per day. With all the care needed judging such a limited investigation the obtained results
342
provide confidence that – proper assay design and assay validation provided – the concept to use
343
(automatable) qNMR and uHPLC interchangeably for quantitative analysis of primary and secondary
344
metabolites from rather crude mixtures, e.g. from fermentation processes, from extraction plants, or
345
within raw material entry control will inspire further investigations in this topic.
346 347
Acknowledgements
348
The authors acknowledge the substantial support by Bionorica SE, Neumarkt, Germany and Bionorica
349
Research GmbH, Innsbruck, Austria.
350
ACS Paragon Plus Environment
Page 12 of 28
Page 13 of 28
351
Journal of Agricultural and Food Chemistry
References
352 353 354 355 356 357 358 359 360 361 362 363
(1) Seger, C.; Sturm S. Analytical aspects of plant metabolite profiling platforms: current standings and future aims. J. Proteome Res. 2007, 6, 480-497. (2) Holmes, E.; Tang, H.; Wang, Y.; Seger C. The assessment of plant metabolite profiles by NMRbased methodologies. Planta Med. 2006, 72, 771-785. (3) Pauli, G. F.; Gödecke, T.; Jaki, B.U.; Lankin, D. C. Quantitative 1H NMR. Development and potential of an analytical method: an update. J. Nat. Prod. 2012, 75, 834-851. (4) Simmler, C.; Napolitano, J. G.; McAlpine, J. B.; Chen, S. N.; Pauli, G. F. Universal quantitative NMR analysis of complex natural samples. Curr. Opin. Biotechnol. 2014, 25 ,51-59. (5) Pauli, G. F.; Jaki, B. U.; Lankin, D. C. Quantitative 1H NMR: development and potential of a method for natural products analysis. J. Nat. Prod. 2005, 68, 133-149. (6) Maggio, R. M.; Calvo, N. L.; Vignaduzzo, S. E.; Kaufman, T. S. Pharmaceutical impurities and
364
degradation products: uses and applications of NMR techniques. J. Pharm. Biomed. Anal. 2014,
365
101, 102-122.
366 367
(7) Liu, N. Q.; Choi, Y. H; Verpoorte, R.; van der Kooy, F. Comparative quantitative analysis of artemisinin by chromatography and qNMR. Phytochem. Anal. 2010, 21, 451-456.
368
(8) Pieri, V.; Sturm, S.; Seger, C; Franz, C.; Stuppner, H. 1H NMR-based metabolic profiling and target
369
analysis: a combined approach for the quality control of Thymus vulgaris. Metabolomics 2012, 8,
370
335-346.
371 372 373
(9) Pieri, V.; Stuppner, H. Quantification of cynaropicrin in artichoke leaf extracts by ¹H NMR spectroscopy. Planta Med, 2011, 77, 1756-1758. (10) Pieri, V.; Belancic, A.; Morales, S.; Stuppner H. Identification and quantification of major steviol
374
glycosides in Stevia rebaudiana purified extracts by 1H NMR spectroscopy. J. Agric. Food Chem.
375
2011, 59, 4378-4384
376
(11) Napolitano, J. G.; Gödecke, T.; Lankin, D. C.; Jaki, B. U.; McAlpine, J. B.; Chen, S. N.; Pauli, G. F.
377
Orthogonal analytical methods for botanical standardization: determination of green tea
378
catechins by qNMR and LC-MS/MS. J. Pharm. Biomed. Anal. 2014, 93, 59-67.
379 380 381 382 383
(12) Flora, K.; Hahn, M.; Rosen, H.; Benner, K. Milk thistle (Silybum marianum) for the therapy of liver disease. Am. J. Gastroenterol. 1998, 93, 139-143. (13) Morazzoni, P.; Bombardelli, E. Silybum marianum (Carduus marianus). Fitoterapia 1995, 66, 342. (14) Madrid, E.; Corchete, P. Silymarin secretion and its elicitation by methyl jasmonate in cell
384
cultures of Silybum marianum is mediated by phospholipase D-phosphatidic acid. J. Exp. Bot.
385
2010, 61, 747-754.
ACS Paragon Plus Environment
Journal of Agricultural and Food Chemistry
386 387 388 389
(15) Fraschini, F.; Demartini, G.; Esposti, D. Pharmacology of Silymarin. Clin. Drug. Invest. 2002, 22, 51-65. (16) Abenavoli, L.; Capasso, R.; Milic, N.; Capasso F. Milk thistle in liver diseases: past, present, future. Phytother. Res. 2010, 24, 1423–1432.
390
(17) AbouZid, S. Silymarin, natural flavonolignans from milk thistle. In: Phytochemicals - A Global
391
Perspective of Their Role in Nutrition and Health; Rao, V. Ed.; InTech: Rijeka, Croatia, 2012; pp
392
255-272.
393
(18) Burczynski, F. J.; Yan, J.; Gong, Y.; Nguyen, D.; Wang, G.; Burczynski, S. D.; Smith, H. J.; Gong, Y.
394
The hepatoprotective effect of diltiazem and silymarin. Nat. Prod. Chem. Res. 2013, 1, 111.
395
(19) Kroll, D. J.; Shaw, H. S.; Oberlies, N. H. Milk thistle nomenclature: Why it matters in cancer
396 397 398 399
research and pharmacokinetic studies. Integr. Cancer Ther. 2007, 6, 110–119. (20) Kuki, A.; Nagy, L.; Deak, G.; Nagy, M., Zsuga, M.; Keki, S. Identification of Silymarin Constituents: An Improved HPLC―MS Method. Chromatographia 2012, 75, 175-180. (21) Saleh, I. A.; Vinatoru, M.; Mason, T. J.; Abdel-Azim, N. S.; Aboutabl, E. A.; Hammouda, F. M.
400
Ultrasonic-assisted extraction and conventional extraction of Silymarin from Silybum marianum
401
seeds; a comparison. Res. J. Pharm. Biol. Chem. Sci. 2015, 6, 709-717.
402 403 404
(22) Eklund, L.; Simon, J. P.; Ballenger, J. High Performance Liquid Chromatography of flavonolignans in commercial milk thistle supplements. Bios 2009, 80, 164-169. (23) Lee, J. I.; Narayan, M.; Barrett, J. S. Analysis and comparison of active constituents in commercial
405
standardized silymarin extracts by liquid chromatography-electrospray ionization mass
406
spectrometry. J. Chromatogr. B 2007, 845, 95-103.
407 408
(24) Kvasnicka, F.; Bíba, B.; Sevcík, R.; Voldrich, M.; Krátká, J. Analysis of the active components of silymarin. J. Chromatogr A. 2003, 990, 239-245.
409
(25) Omar, A.; Hadad, G.; Badr J. First detailed quantification of silymarin components in the leaves of
410
Silybum marianum cultivated in Egypt during different growth stages. Acta Chromtogr. 2012, 24,
411
463-474.
412 413 414
(26) Wianowska, D.; Wiśniewski, M. Simplified procedure of silymarin extraction from Silybum marianum L. Gaertner. J. Chromatogr. Sci. 2015, 53, 366-372. (27) Wang, K.; Zhang, H.; Shen, L.; Du, Q.; Li, J. Rapid separation and characterization of active
415
flavonolignans of Silybum marianum by ultra-performance liquid chromatography coupled with
416
electrospray tandem mass spectrometry. J. Pharm. Biomed. Anal. 2010, 53, 1053-1057.
417
(28) Lee, J. I.; Hsu, B. H.; Wu, D.; Barrett, J. S. Separation and characterization of silybin, isosilybin,
418
silydianin and silychristin in milk thistle extract by liquid chromatography-electrospray tandem
419
mass spectrometry. J. Chromatogr. A 2006, 1116, 57-68.
ACS Paragon Plus Environment
Page 14 of 28
Page 15 of 28
420
Journal of Agricultural and Food Chemistry
(29) Shibano, M.; Lin, A. S.; Itokawa, H.; Lee, K. H. Separation and characterization of active
421
flavonolignans of Silybum marianum by liquid chromatography connected with hybrid ion-trap
422
and time-of-flight mass spectrometry (LC-MS/IT-TOF). J. Nat. Prod. 2007, 70, 1424-1428.
423
(30) Brinda, B. J.; Zhu, H. J.; Markowitz, J. S. A sensitive LC-MS/MS assay for the simultaneous analysis
424
of the major active components of silymarin in human plasma. J. Chromatogr. B 2012, 902, 1-9.
425
(31) N.N. (2006) ICH Harmonized Tripartite Guideline Validation of Analytical Procedures: Text and
426
Methodology Q2(R1); http://www.ich.org/products/guidelines/quality/article/quality-
427
guidelines.html (accessed August 23, 2015)
428
(32) Napolitano, J. G.; Lankin, D. C.; Graf, T. N.; Friesen, J. B.; Chen, S. N.; McAlpine J. B.; Oberlies, N.
429
H; Pauli, G.F. HiFSA Fingerprinting applied to isomers with near-identical NMR spectra: The
430
Silybin/Isosilybin case. J. Org. Chem. 2013, 78, 2827-2839.
431 432 433
(33) Bilia, A.R.; Bergonzi, M.C.; Gallori, S.; Mazzi, G.; Vincieri, F. F. J. Pharm. Biomed. Anal. 2002, 30, 613-624. (34) Passing, H.; Bablok, W. A new biometrical procedure for testing the equality of measurements
434
from two different analytical methods. Application of linear regression procedures for method
435
comparison studies in Clinical Chemistry, Part I. J. Clin. Chem. Clin. Biochem. 1983, 21, 709-720.
436
(35) Dewitte, K.; Fierens, C.; Stöckl, D.; Thienpont, L. M. Application of the Bland-Altman plot for
437
interpretation of method-comparison studies: a critical investigation of its practice. Clin. Chem.
438
2002, 48, 799-801.
439
ACS Paragon Plus Environment
Journal of Agricultural and Food Chemistry
440
Figures Legends
441 442
Figure 1: Structures of major S. marianum flavonolignans (1 - 6). 1H NMR resonances of highlighted
443
protons were used for qNMR based analyte quantitation.
444 445
Figure 2: Representative 1H-NMR spectra of selected S. marianum AC and PL extracts recorded under
446
qNMR conditions. Highlighted inserts show NMR peaks (top in sample AC03, bottom in PL01) used
447
for analyte quantitation, numbers correspond to analytes 1 – 6.
448 449
Figure 3: Representative S. marianum uHPLC-DAD chromatograms of AC and PL extracts recorded at
450
286 nm. Numbers correspond to analytes 1 – 6.
451 452
Figure 4: Passing-Bablok correlation (left-hand column) and Bland-Altman plot (right-hand column)
453
analysis of AC (top row) and PL (bottom row) samples.
454
ACS Paragon Plus Environment
Page 16 of 28
Page 17 of 28
Journal of Agricultural and Food Chemistry
Tables
Table 1: 1H-NMR signal information (methanol-d4, 600 MHz) for NMR resonances used within the qNMR experiments for silymarin flavonolignans quantitation.
Compound
δH (ppm)
multiplicity, JHH (Hz)
T1 (sec)
1
5.55
d, 6.28
1.5
2
2.89
br s, -
2.0
3
7.10
d, 2.02
2.1
4
7.08
d, 2.02
2.1
5
7.15
d, 2.08
2.2
6
7.13
d, 2.08
2.2
ACS Paragon Plus Environment
Journal of Agricultural and Food Chemistry
Page 18 of 28
Table 2: qNMR based calibration function parameters for 1-6, including regression equations, correlation coefficients (R2), linearity range, LOD and LOQ.
Linearity range
LOD
(μg/mL)
(μg/mL)
0.9986
25-2000
2.9
8.8
y = 0.1476x - 0.0016
0.9998
10-2000
3.4
10.4
3
y = 0.1882x - 0.0002
0.9999
10-2000
0.6
1.8
4
y = 0.1878x - 0.0005
0.9994
25-2000
0.5
1.6
5
y = 0.1848x - 0.0004
1.0000
25-2000
2.3
7.0
6
y = 0.1626x + 0.0046
0.9997
10-2000
1.4
4.3
2
Y (peak area) = kcA+d
R
1
y = 0.1570x + 0.0189
2
Compound
LOQ (μg/mL)
ACS Paragon Plus Environment
Page 19 of 28
Journal of Agricultural and Food Chemistry
Table 3: uHPLC-DAD based calibration function parameters for 1-6, including regression equations, correlation coefficients (R2), linearity range, LOD and LOQ.
Ret. Time Compound (min)
Linearity
2
Y (peak area) = kcA+ d
R
LOD
LOQ
range (μg/mL) (μg/mL) (μg/mL)
1
5.86
y= 2374.1x – 8.4163
0.9993
10-2000
1.0
3.0
2
6.60
y= 1886.7x – 22.4990
0.9997
25-2000
4.0
12.0
3
8.38
y= 2667.6x – 26.5870
0.9994
25-2000
6.0
19.0
4
8.74
y= 2528.0x – 5.2708
0.9992
10-2000
1.0
3.0
5
10.01
y= 2710.3x – 60.6340
0.9994
25-2000
1.0
2.0
6
10.38
y= 2598.4x – 63.9100
0.9993
25-2000
1.0
3.0
ACS Paragon Plus Environment
Journal of Agricultural and Food Chemistry
Page 20 of 28
Table 4: qNMR and uHPLC-DAD assay accuracy for compounds 2, 3, and 4. Inter-day recovery values (n = 9) expressed in percentage amount added (mean ± RSD).
Compound
NMR Recovery
Compound for
uHPLC Recovery
for NMR
(% ± RSD)
uHPLC
(% ± RSD)
Amount added (µg) AC 500.0 250.0
AC 104.9 ± 2.7
2
91.2 ± 6.7
98.9 ± 3.1 2
100.7 ± 6.8
100.0
90.7 ± 5.1
103.1 ± 6.1
500.0
98.4 ± 0.7
97.2 ± 2.1
250.0
3
100.0
106.4 ± 2.2 101.8 ± 0.6
PL 500.0 250.0
3
103.6 ± 1.5 PL
95.3 ± 2.2 2
98.4 ± 2.0
91.2 ± 5.4
101.3 ± 3.6 2
94.7 ± 4.5
100.0
103.6 ± 3.3
102.1 ± 4.2
500.0
103.1 ± 0.9
106.5 ± 3.0
250.0 100.0
4
105.5 ± 1.0
4
100.8 ± 1.9
ACS Paragon Plus Environment
91.8 ± 2.4 110.0 ± 0.8
Page 21 of 28
Journal of Agricultural and Food Chemistry
Table 5: qNMR and uHPLC-DAD assay precision for compounds 2, 3, and 4. Intra-day (n = 3 on each day) and inter-day (n = 9) precision for both types of extracts.
RSD (%) Compound
Day 1 Day 2 Day 3
Day 1-3
NMR 2 (AC)
7.7
3.4
6.6
8.6
3 (AC)
1.4
0.9
2.0
1.4
2 (PL)
3.5
3.5
5.2
4.0
4 (PL)
2.3
1.4
2.9
2.3
uHPLC 2 (AC)
3.3
6.0
7.4
5.3
3 (AC)
3.8
1.4
1.8
2.9
2 (PL)
3.3
5.2
2.2
5.9
4 (PL)
2.5
2.9
2.8
3.0
ACS Paragon Plus Environment
Journal of Agricultural and Food Chemistry
Page 22 of 28
Table 6: Quantitative qNMR results for compounds 1 – 6 in S. marianum samples AC 01-09 and PL 0110 (RSD in parenthesis, n=3).
Samples
1
2
3+4
5+6
AC (µg / 10 mg AC material) AC 01
1388.5 (2.1)
43.1 (9.6)
3718.5 (1.5)
676.7 (1.0)
AC 02
1458.7 (0.67) 41.2 (6.8)
3854.7 (0.5)
659.3 (0.7)
AC 03
1593.8 (0.9)
46.0 (5.8)
3311.2 (0.6)
735.2 (1.7)
AC 04
1577.5 (0.9)
40.2 (6.0)
3410.4 (0.4)
729.0 (1.5)
AC 05
1666.4 (1.1)
41.0 (5.8)
3583.9 (0.8)
749.8 (3.3)
AC 06
1683.2 (1.0)
37.2 (1.7)
3619.8 (0.5)
756.5 (1.6)
AC 07
1869.9 (0.9)
39.6 (6.3)
3721.6 (0.9)
843.0 (2.5)
AC 08
1703.5 (0.9)
42.6 (5.5)
3411.3 (1.2)
778.1 (4.6)
AC 09
1752.4 (0.2)
47.5 (6.1)
3284.8 (1.8)
788.5 (0.6)
PL (µg / 100 mg PL) PL 01
372.3 (4.4)
34.0 (4.3)
860.7 (2.8)
221.5 (7.1)
PL 02
106.0 (6.5)
224.3 (6.3)
221.6 (8.2)
169.5 (5.3)
PL 03
114.3 (8.9)
169.0 (9.4)
241.1 (5.9)
181.4 (1.8)
PL 04
158.6 (9.0)
249.9 (4.2)
306.9 (1.0)
207.5 (3.2)
PL 05
144.0 (8.6)
260.0 (4.9)
311.0 (2.7)
201.9 (8.7)
PL 06
128.1 (7.1)
287.0 (6.6)
211.3 (4.6)
214.0 (3.6)
PL 07
151.0 (8.9)
301.2 (5.5)
301.4 (5.2)
201.2 (5.6)
PL 08
104.0 (2.3)
217.0 (3.7)
247.8 (6.7)
182.8 (7.2)
PL 09
148.0 (6.9)
273.2 (7.5)
281.8 (3.5)
196.4 (6.7)
PL 10
137.0 (6.5)
321.7 (6.1)
280.9 (5.5)
200.3 (7.8)
ACS Paragon Plus Environment
Page 23 of 28
Journal of Agricultural and Food Chemistry
Table 7: Quantitative uHPLC-DAD results for compounds 1-6 in S. marianum samples AC 01-09 and PL 01-10 (RSD in parenthesis, n=3) Sample
1
2
AC 01
1387.2 (1.7)
43.5 (6.9)
AC 02
1421.4 (4.2)
AC 03
3
4
3+4
5
6
5+6
AC (µg / 10 mg AC material) 1463.9 (1.9) 2293.1 (1.8) 3757.0
466.2 (1.5)
172.7 (1.8)
638.9
41.7 (7.5)
1527.1 (3.8)
2379.7 (3.9)
3906.8
477.7 (4.2)
183.8 (2.8)
661.5
1576.8 (0.8)
48.3 (3.5)
1275.8 (1.4)
2120.6 (1.2)
3396.4
528.0 (0.6)
176.7 (1.8)
704.7
AC 04
1566.5 (2.8)
41.6 (4.2)
1298.0 (2.9)
2140.1(3.0)
3438.1
527.2 (2.9)
181.8 (3.6)
709.0
AC 05
1644.0 (2.5)
40.6 (6.6)
1332.0 (2.9)
2207.0 (3.0)
3539.0
550.2 (2.7)
189.1 (3.7)
739.3
AC 06
1670.1 (5.1)
38.6 (6.2)
1389.6 (5.0)
2290.8 (5.1)
3680.4
561.8 (4.8)
193.2 (4.9)
754.8
AC 07
1868.7 (7.7)
41.6 (6.7)
1400.2 (7.4)
2378.6 (7.5)
3778.8
607.5 (7.1)
214.4 (5.4)
821.9
AC 08
1719.5 (2.1)
41.3 (3.6)
1279.1 (1.9)
2174.4 (2.0)
3453.5
559.9 (1.5)
196.4 (1.8)
756.4
AC 09
1710.1 (2.2)
48.6 (3.2)
1251.9 (2.2)
2141.4 (2.2)
3393.3
557.6 (2.1)
194.7 (2.2)
752.3
PL 01
382.1 (1.44)
30.5 (9.7)
PL (µg / 100 mg plant material) 325.8 (1.3) 536.9 (1.9) 862.8
153.3 (0.8)
70.4 (0.8)
223.7
PL 02
113.5 (7.89)
228.1 (5.4)
82.2 (5.6)
135.1 (6.2)
217.3
96.6 (4.2)
77.6 (2.8)
174.1
PL 03
118.8 (6.18)
162.1 (7.9)
98.0 (7.0)
139.0 (8.3)
237.0
105.0 (3.5)
84.0 (4.5)
189.0
PL 04
158.3 (1.91)
253.3 (5.9)
108.2 (8.6)
191.2 (1.5)
299.4
116.4 (0.6)
89.2 (0.5)
205.6
PL 05
146.1 (5.75)
243.5 (7.5)
103.3 (5.0)
190.8 (4.0)
294.1
123.5 (3.0)
89.9 (2.5)
213.4
PL 06
120.9 (8.12)
290.3 (5.4)
72.7 (8.4)
129.4 (8.2)
202.1
125.1 (4.2)
96.2 (2.9)
221.3
PL 07
157.5 (9.09)
300.7 (3.7)
107.2 (7.8)
184.6 (8.0)
291.8
120.7 (7.2)
87.5 (6.5)
208.2
PL 08
102.2 (2.69)
211.0 (8.7)
116.2 (6.3)
136.7 (3.9)
252.9
101.0 (6.2)
77.3 (6.6)
178.3
PL 09
141.0 (8.44)
264.2 (5.2)
106.0 (7.4)
170.0 (8.0)
276.0
109.0 (7.6)
86.0 (8.4)
195.0
PL 10
142.8 (8.51)
316.1 (9.7)
103.2 (9.0)
169.7 (8.6)
272.9
117.1 (4.3)
84.4 (4.8)
201.5
ACS Paragon Plus Environment
Journal of Agricultural and Food Chemistry
Page 24 of 28
Figures
OH O HO
H
O
HO
OH OH O silydianin (2)
silychristin (1) H H
H O
O
CH2OH
H OCH3
O OH
H
HO
O
O
CH2OH OCH3
O
OH
OH
OH O
H
OH
OH O
silybin A (3)
silybin B (4) OH
H H HO
OCH3
H
O
OH OH O
HO
OH
HO
OCH3 CH2OH
O
O
O
O
O OH
H
OCH3 CH2OH
HO
O
O
H
OH O
OH
H
O OH
H
OH O isosilybin A (5)
isosilybin B (6)
Figure 1
ACS Paragon Plus Environment
OCH3 CH2OH
Page 25 of 28
Journal of Agricultural and Food Chemistry
3,4 2
1 5,6
7.14
7.12
7.10
7.08
5.58
5.56
5.54
5.52
2.92
2.90
2.88
2.86
AC03 IS
PL01
8
6
4
2
0
[ppm] 3,4
1
2
5,6
7.14
7.12
7.10
7.08
5.58
5.56
5.54
5.52
Figure 2
ACS Paragon Plus Environment
2.92
2.90
2.88
2.86
Journal of Agricultural and Food Chemistry
Figure 3
ACS Paragon Plus Environment
Page 26 of 28
Page 27 of 28
Journal of Agricultural and Food Chemistry
4000
(uHPLC-qNMR)/Mean of methods (%)
6
AC – samples
Analyte concentration qNMR (µg / 10 mg extract)
y = 0.99 x + 2.91; rho = 0.996
3000
2000
1000
AC - samples
2
0
Mean = - 0.1 % -2
-4
-2SD = -5.2% -6
0 0
1000
2000
3000
4000
0
1000
Analyte concentration uHPLC (µg / 10 mg extract)
2000
3000
4000
Mean of methods (µg / 10 mg extract)
(uHPLC-qNMR)/Mean of methods (%)
10
1000
PL – samples y = 1.02 x – 4.61; rho = 0.994
Analyte concentration qNMR (µg / 100 mg plant)
+2SD = 5.0%
4
800 600 400 200
PL - samples +2SD = 7.5%
5
0
Mean = - 0.3 %
-5 -2SD = -8.1% -10
0 0
200
400
600
800
1000
0
200
Analyte concentration uHPLC (µg / 100 mg plant)
Figure 4
ACS Paragon Plus Environment
400
600
Mean of methods (µg / 100 mg plant material)
800
1000
Journal of Agricultural and Food Chemistry
”TOC” graphic
ACS Paragon Plus Environment
Page 28 of 28