Simultaneous quantification of bioactive components in Chinese herbal spirits by ultra-high performance liquid chromatography coupled to triple-quadrupole mass spectrometry (UHPLC–QQQ–MS/MS)

Abstract

Background

The Chinese medicinal wine made from herbal medicines became prevalent among Chinese people. The Chinese herbal spirit is composed of several herbal extracts, and has the certain health functions, such as anti-fatigue and immune regulation. The quality evaluation of Chinese herbal spirit is greatly challenged by the enormous and complex components with great structural diversity and wide range of concentration distribution.

Methods

An ultra-high performance liquid chromatography coupled to triple quadrupole mass spectrometry (UHPLC-QQQ-MS/MS) with multiple reaction monitoring (MRM) method was developed to simultaneously determine forty-three bioactive components in the Chinese herbal spirits produced by year 2014 and 2018.

Results

Quantitative results showed that 11 components, i.e.., puerarin (5), purpureaside C (7), daidzin (8), echinacoside (9), acteoside (15), epimedin B (22), epimedin C (23), icariin (24), eugenol (27), chikusetsusaponin iva (30) and Z-ligustilide (40), significantly decreased along with the increasing years of storage, while 5 compounds, i.e.., geniposidic acid (1), protocatechuic acid (2), crustecdysone (14), daidzein (18) and icariside I (35), were basically stable in all samples across the years.

Concusion

The established method allowing to simultaneously determined 43 components with wide structural diversity and trace amounts will facilitate the quality control research of Chinese herbal spirits.

Introduction

Chinese medicinal wine, an alcoholic beverage, is commonly produced by soaking precious Chinese medicinal materials. It became prevalent among Chinese people due to its nourishing and tonic functions. A Chinese herbal spirit from Jing Brand Co., Ltd is one of the most popular medicinal wine in China, which is brewed with the faint-scented Xiaoqu liqueur and several authorized herbal extracts under the guidance of traditional Chinese medicine theory, including Dioscoreae rhizoma, Curculiginis rhizoma, Angelicae sinensis radix, Cistanches herba, Lycii fructus, Astragali radix, Epimedii folium, Cinnamomi cortex, and Caryophylli flos, etc. [1]. This Chinese herbal spirit has health-care functions of anti-fatigue and immune enhancement [1]. It was shown to strengthen the immune system in the Shen-yang deficient rats, which was associated with the activation of hypothalamic-pituitary-adrenal axis [2]. Furthermore, a clinic study indicated that the Chinese herbal spirit could relieve both physical fatigue and mental fatigue of the patients with the fatigued sub-health status [3].

According to enterprise criterion, total saponins, total flavonoids and icariin are current chemical markers for quality control of this herbal spirit. However, its quality evaluation is greatly challenged by the enormous and complex components with great structural diversity and wide range of concentration distribution. More than 150 ingredients were identified from the herbal spirit, including flavonoids, saponins, alkaloids, phenylethanoids, coumarins, anthraquinones and volatile oil, etc. To date, gas chromatography-mass spectrometry (GC-MS) [1, 4, 5] and liquid chromatography-mass spectrometry (LC-MS) [6], have been established to quantify the volatile flavoring substances and active ingredients in the spirit. However, the very limited quantity of ingredients were quantified. Ultra-high performance liquid chromatography coupled to triple quadrupole mass spectrometry (UHPLC-QQQ-MS) operated in multiple reaction monitoring (MRM) mode is an effective quantification method owing to its well-known high sensitivity and specificity, which could avoid the interference from the background matrix [7,8,9]. It has been successfully utilized to quantify bioactive components in the complex systems [10,11,12,13,14,15]. In the current study, therefore, an UHPLC-QQQ-MS method was developed to simultaneously quantify 43 bioactive components in the Chinese herbal spirit samples produced by year 2014 and 2018, and their concentrations in different peoduction years were also compared.

Materials and methods

Materials and chemicals

The Chinese herbal spirit samples produced by year 2014 (n = 20) and 2018 (n = 20) were kindly provided by Jing brand Co. Ltd (Hubei, China). The voucher specimens were deposited at room temperature and shielded from light at Institute of Chinese Medical Sciences, University of Macau, Macao. Forty-three reference compounds listed in Table 1 were purchased from Baoji herbest Bio-Tech Co. Ltd (Baoji, China). Their chemical structures were shown in Additional file 1: Figure S1, and the purities of reference standards were over 98 % as confirmed by HPLC-UV. HPLC-grade acetonitrile and methanol were purchased from Merck (Darmstadt, Germany). Formic acid was obtained from Aladdin Industrial Inc. (Shanghai, China). Deionized water was purified using a Millipore Milli-Q purification system (Bedford, MA, USA).

Table 1 Mass parameters, calibration curves, LOD, LOQ, precision and recovery of 43 investigated analytes by LC-MS/MS with MRM mode

Sample preparation and standard solution preparation

An aliquot of 1mL of the spirit was diluted with the equivalent volume of acetonitrile, and vortexed for 1 min. The mixture was centrifuged at 14,800 rpm for 20 min. After the centrifugation, the supernatant was filtered through a 0.22 μm filter (PVDF Millex-GV, 13 mm, Millipore) prior to the quantitative analysis. All of 43 reference standards were accurately weighed and dissolved in methanol to prepare individual stock solutions at the concentrations ranging from 0.58 to 2.29 mg mL^− 1. The mixed standard solution was prepared by mixing appropriate volumes of the individual stock solutions and further diluted to a series of proper concentrations with methanol.

LC-MS/MS analysis

Forty bathes of the spirit samples were analyzed by an ACQUITY™ UPLC system (Waters, Milford, MA, US) coupled with a Xevo TQD triple-quadrupole tandem mass spectrometry (QQQ-MS/MS, Waters Co., Manchester, UK). Chromatographic separation was implemented on an ACQUITY UPLC BEH C₁₈ column (150 × 2.1 mm, 1.7 μm). The mobile phases were consisted of 0.1 % (v/v) aqueous formic acid solution (phase A) and acetonitrile containing 0.1 % formic acid (phase B) at the flow-rate of 0.3 mL min⁻¹ using a gradient elution program as follows: 5–40 % B at 0–12 min, 40–100 % B at 12–16 min, isocratic 100 % B for 2 min, and the re-equilibrated by 5 % B for 3 min. The sample injection volume was set at 2 µL. The temperatures of column and injector were set at 35 ℃ and 8 ℃, respectively.

Data acquisition was performed by a Xevo TQD QQQ-MS equipped with an electrospray ionization (ESI) using MRM mode. MS was operated in either the positive (ESI+) or negative mode (ESI-) to obtain satisfactory MS response for 43 investigated compounds due to their different properties. The MS and MS/MS spectra of each compound were acquired using the mixed standard solution. The optimized MS parameters were as follows: capillary voltage, 3.5 kV (positive ion mode) and − 3.0 kV (negative ion mode); source temperature, 140 ℃; desolvation gas flow and temperature, 650 L h⁻¹ and 350 ℃; cone flow, 50 L h⁻¹. The ion transitions, cone voltage, and collision energy for each compound were optimized and shown in Table 1. All instrumentations were synchronized and controlled by Waters Masslynx software (version, 4.1).

Method validation

To evaluate sensitivity and precision of the established UHPLC-QQQ-MS/MS method, the linearity, limit of detection (LOD), limit of quantitation (LOQ), precision and recovery of 43 analytes were tested. The calibration curve of each compound was constructed by plotting the peak areas against the concentrations using the mixed standard solution at a series of concentrations. The precision was examined by calculating intra- and inter-day variations of each analyte using the mixed standards for five replicates within a day and three consecutive days. The LOD and LOQ for each analyte were estimated at signal-to-noise ratio (S/N) of about 3 and 10, respectively. The accuracy of the established method was assessed by spike recovery experiments. A known amount (equal to the content for each analyte in the sample) of the mixed standards was spiked into the random spirit sample (S-2018-08). The sample was prepared with six replicates and analyzed by the method mentioned above. The recovery (%) was calculated as the following equation:

$${\text{Recovery }}\left( \% \right)\; = \;{1}00\; \times \;\left( {{\text{detected amount}}\; - \;{\text{original amount}}} \right)/{\text{ spiked amount}}$$

Statistical analysis

The concentrations of 43 analytes in the spirit samples were presented as mean ± standard deviation (SD). The difference between groups was assessed by student t-test using a GraphPad Prism package (version 6.0, San Diego, CA, USA), and a p-value of less than 0.05 was considered statistically significant. Orthogonal partial least squares-discriminant analysis (OPLS-DA), a supervised multiple regression, was conducted to discriminate the spirit samples manufactured in different years according to the levels of investigated analytes by SIMCA-P software (version 14.1, Umetrics, Umeå, Sweden).

Results and discussion

LC-MS/MS method development

The quality control research related to this Chinese herbal spirit mainly focused on the determination of volatile components by GC-MS [1, 4, 5], and very limited number of nonvolatile bioactive components were quantified by LC-MS [6]. Due to consisting of multiple herbal extracts and the faint-scented Xiaoqu liqueur, the herbal spirit is a very complex matrix containing the numerous organic and inorganic compounds with wide range of concentrations. Thus, the simultaneous quantitation of forty-three compounds with various chemical types encounters the great challenge in short running time using UHPLC. MRM is a highly specific technique for quantifying the targeted analyte, regardless of baseline chromatographic separation. The targeted analyte in chromatographic co-elution could be accurately quantified if they have different MS or MS/MS characteristics. However, the co-eluted analytes may cause the potential mutual ionization suppression in ESI, leading to the low MS response. Therefore, it is also necessary to optimize the chromatographic conditions, including column and mobile phase, to achieve the high sensitivity and fast separation in LC-MS/MS analysis. Three UHPLC columns, such as BEH C₁₈ column, BEH HILIC column and HSS T3 C₁₈ column, were examined. As results, an ACQUITY BEH C₁₈ column was most suitable for the separation of the targeted compounds in the sample owing to the best resolution and the most peak capacity. Furthermore, several types of mobile phases, including methanol/water and acetonitrile/water system supplemented with various modifiers, were tested. The results shown that 0.1 % aqueous formic acid solution / acetonitrile with 0.1 % formic acid was the optimum mobile phases to obtain the chromatogram with the best resolution.

It is critical to design ion transition of each analyte, including precursor ions and their corresponding product ions, in MRM analysis. The full scan was used to select the precursor ions using their reference standards, and the dominated fragment ion in daughter scan was chosen as the corresponding product ion (Fig. 1). Taking sagittatoside A (32, MW = 676.24), a main component derived from Epimedii folium, as an example, the protonated ion [M + H]⁺ m/z 677.33 was presented with the highest abundance in full scan spectrum, the adduct ion [M + Na]⁺ (m/z 699.34), was also observed with less intensity. While, the product ion m/z 369.24 was dominated in the corresponding daughter scan of [M + H]⁺ (Fig. 1 a, b). The mass difference between parent and product ions was 308.09, which corresponded to the loss of one glucose and one rhamnose units [15]. Herein, the ion transition (m/z 677.33 → 369.24) was selected to quantify sagittatoside A in the liqueur by MRM. Likewise, the ion transition (m/z 417.36 → 255.22) was optimized to determine daidzin (8, Fig. 1 c, d). To achieve maximum signal, both positive and negative ion modes were tested and compared. The results indicated that all compounds, except Chikusetsusaponin iva (30), shown the higher sensitivities in the positive ion mode than in the negative mode. For example, MRM chromatograms of daidzein (18) and puerarin (5) in positive ion mode were remarkably higher than that in negative mode (Fig. 2a, b). However, the higher MS intensity of chikusetsusaponin iva was observed in negative mode (Fig. 2c). Furthermore, the transition of chikusetsusaponin iva in the negative mode were further optimized. We found that multiple ion monitoring (MIM), i.e. m/z 793.57 → 793.57, was more suitable for the detection of chikusetsusaponin iva, compared to other transitions, such as m/z 793.57 → 631.67 (Fig. 2d).

Fig. 1

Full scan spectra (left), the corresponding daughter scan spectra and the proposed fragmentation pattern (right) of sagittatoside A (32, a) and daidzin (18, b)

Fig. 2

The MRM chromatograms of daidzein (18, a), puerarin (5, b) and Chikusetsu saponin iva (30, c) in both positive and negative ion modes; The MRM chromatograms of Chikusetsu saponin iva (d) with different transitions (m/z 793.57→793.57 and m/z 793.57→631.67) in the negative ion mode 

Additionally, cone voltage (CV) and collision energy (CE), the important factors that affect the sensitivity of UHPLC-QQQ-MS/MS analysis, were also optimized for each analyte using reference standard. CV and CE were optimized from 10 V to 50 V with a step of 10 V and from 5 V to 50 V with a step of 5 V, respectively. Taking sagittatoside A as an example, the signal intensity of product ion m/z 369.24 increased along with CV from 10 V to 30 V or CE from 5 V to 15 V, then decreased with the increasing voltages (Fig. 3a). Therefore, 30 V of CV and 15 V of CE were chosen for the quantification of sagittatoside A. Likewise, 30 V of CV and 20 V of CE were used to determine the daidzin (Fig. 3b). In a similar manner, the mass spectrometry conditions of all analytes were optimized and listed in Table 1. Under the optimized LC-MS/MS conditions, 43 analytes were well separated and detected in 18 min. The representative MRM chromatograms of the mixed standards and the Chinese herbal spirit are shown in Fig. 4.

Fig. 3

The intensity distributions of the fragment ions derived from sagittatoside A (m/z 677.33→369.24, a) and daidzin (m/z 417.36→255.22, b) with the cone voltage ranging from 10 V to 50 V (left), and collision energy ranging from 5 to 50 V (right) in positive ion mode

Fig. 4

The MRM chromatograms of glycosides (left) and the remaining analytes (right) in the mixed standards (a) and the Chinese herbal spirit sample (S-2018-08, b)

Method validation

The established LC-MS/MS method was validated by a series of experiments, including linearity, sensitivity, precision, and accuracy. As shown in Table 1, the calibration curves of all analytes exhibited good linear regression (R² ≥ 0.9911) within the wide dynamic range. The LODs and LOQs of analytes were less than 8.54 and 27.7 ng mL^− 1, respectively. Other than that of geniposide was 23.1 and 46.2 ng mL^− 1, respectively, which was higher than other compounds. The overall intra-day and inter-day variations were lower than 4.89 % and 5.81 %, respectively. Additionally, the developed method had the acceptable accuracy with recoveries ranging from 90.3 to 111.7 %. Taken together, the proposed LC-MS/MS method is sensitive, precise and accurate for the simultaneous determination of these 43 compounds in the Chinese herbal spirit produced in different years.

Quantification of 43 compounds in the Chinese herbal spirit

The validated method was successfully applied to quantify 43 compounds in the liqueur samples. As shown in Table 2 and 35 of 43 analytes were detected and quantified in the liqueur samples produced in both 2014 and 2018. Eight analytes, including scopolin (4), salvianolic acid A (19), quercetin (20), rhamnocitrin (29), geniposide (31), desmethyl icaritin (38), icaritin (41), tanshinone IIA (42), were not detected in all samples. Eight compounds, including eugenol (27), puerarin (5), icariin (24), epimedin C (23), epimedin B (22), chikusetsusaponin iva (30), daidzin (8) and crustecdysone (14) were identified as major components in the liqueur with contents more than 1 µg mL^− 1. Among them, eugenol (27) is a major components with the highest abundance in the liqueur samples. While, daidzein (18), formononetin (25), baohuoside II (26) and icariside I (35) were trace components with the content of less than 0.1 µg mL^− 1. Additionally, 26 analytes were observed to be significantly different between the Chinese herbal spirits produced in year 2014 and 2018 (p < 0.05).

Table 2 The contents of 43 investigated components in the Chinese herbal spirits (ND = not detected)

In order to further visualize the difference between samples, OPLS-DA, a supervised multivariate data analysis, was further constructed to characterize the differences between groups on the basis of the levels of 35 analytes determined. As illustrated in Fig. 5a, the spirit samples produced at 2014 and 2018 were unambiguously segregated into two tight clusters (R²X = 0.826, R²Y  = 0.96, Q² = 0.887), suggesting that the great alteration in the investigated compounds were presented. The cumulative values of R²X, R²Y, and Q² were close to optimal value of 1.0, indicating the established models with excellent predictive capability and fitness [16]. To identify the differentiated components that contribute most to the group separation, the differentiated compounds were selected by S-plot derived from the constructed OPLS-DA. Eleven compounds with variable importance in the projection (VIP) of more than 1 (VIP > 1) and p value of less than 0.05 were highlighted in the S-plot (Fig. 5b). The group separation of the samples with the different production years could also be clearly classified in OPLS-DA score plot according to the contents of these 11 highlighted components (Fig. 5c). Compared to the spirit produced in 2018, the contents of the highlighted compounds, puerarin (5), purpureaside C (7), daidzin (8), echinacoside (9), acteoside (15), epimedin B (22), epimedin C (23), icariin (24), eugenol (27), chikusetsusaponin iva (30) and Z-ligustilide (40), significantly decreased in the samples produced in year 2014 (Fig. 6a). Geniposidic acid (1), protocatechuic acid (2), crustecdysone (14), daidzein (18) and icariside I (35) around the origin of S-plot were also highlighted as the stable compounds, which shown the comparable contents with less variation between samples from two production years (Fig. 6b). Additionally, a heatmap according to the relative contents of 35 detedcted analytes was constructed to display the changes in the analyte levels between groups (Fig. 6c).

Fig. 5

OPLS-DA score plots (a) and the corresponding S-plot (b) based on the contents of 43 analytes investigated in the Chinese herbal spirit samples produced in year 2014 and 2018. OPLS-DA score plot (c) according to the contents of 11 highlighted differentiated ions across the samples. The differentiated ions (VIP > 1 and p < 0.05) were marked in red. The compound number in the S-plot represent in the same manner as in Table 1

Fig. 6

The content comparisons of the highlighted analytes (a) and the compounds without difference (b). Heatmap (c) visualizing the differences in the contents of the investigated analytes across the samples. The compound number in the S-plot represent in the same manner as in Table 1

Conclusions

In this study, a rapid and sensitive UHPLC-QQQ-MS with MRM mode was developed to simultaneously quantify 43 bioactive components in Chinese herbal spirits. Quantitative results showed that 11 components, i.e.., puerarin (5), purpureaside C (7), daidzin (8), echinacoside (9), acteoside (15), epimedin B (22), epimedin C (23), icariin (24), eugenol (27), chikusetsusaponin iva (30) and Z-ligustilide (40), significantly decreased along with the increasing years of storage, while 5 compounds, i.e.., geniposidic acid (1), protocatechuic acid (2), crustecdysone (14), daidzein (18) and icariside I (35), were basically stable in all samples across the years. The established method allowing to simultaneously determined 43 components with wide structural diversity and trace amounts will facilitate the quality control research of Chinese herbal spirits.