Monitoring and Evaluating the Quality Consistency of Compound Bismuth Aluminate Tablets by a Simple Quantified Ratio Fingerprint Method Combined with Simultaneous Determination of Five Compounds and Correlated with Antioxidant Activities

A combination method of multi-wavelength fingerprinting and multi-component quantification by high performance liquid chromatography (HPLC) coupled with diode array detector (DAD) was developed and validated to monitor and evaluate the quality consistency of herbal medicines (HM) in the classical preparation Compound Bismuth Aluminate tablets (CBAT). The validation results demonstrated that our method met the requirements of fingerprint analysis and quantification analysis with suitable linearity, precision, accuracy, limits of detection (LOD) and limits of quantification (LOQ). In the fingerprint assessments, rather than using conventional qualitative “Similarity” as a criterion, the simple quantified ratio fingerprint method (SQRFM) was recommended, which has an important quantified fingerprint advantage over the “Similarity” approach. SQRFM qualitatively and quantitatively offers the scientific criteria for traditional Chinese medicines (TCM)/HM quality pyramid and warning gate in terms of three parameters. In order to combine the comprehensive characterization of multi-wavelength fingerprints, an integrated fingerprint assessment strategy based on information entropy was set up involving a super-information characteristic digitized parameter of fingerprints, which reveals the total entropy value and absolute information amount about the fingerprints and, thus, offers an excellent method for fingerprint integration. The correlation results between quantified fingerprints and quantitative determination of 5 marker compounds, including glycyrrhizic acid (GLY), liquiritin (LQ), isoliquiritigenin (ILG), isoliquiritin (ILQ) and isoliquiritin apioside (ILA), indicated that multi-component quantification could be replaced by quantified fingerprints. The Fenton reaction was employed to determine the antioxidant activities of CBAT samples in vitro, and they were correlated with HPLC fingerprint components using the partial least squares regression (PLSR) method. In summary, the method of multi-wavelength fingerprints combined with antioxidant activities has been proved to be a feasible and scientific procedure for monitoring and evaluating the quality consistency of CBAT.

A combination method of multi-wavelength fingerprinting and multi-component quantification by high performance liquid chromatography (HPLC) coupled with diode array detector (DAD) was developed and validated to monitor and evaluate the quality consistency of herbal medicines (HM) in the classical preparation Compound Bismuth Aluminate tablets (CBAT). The validation results demonstrated that our method met the requirements of fingerprint analysis and quantification analysis with suitable linearity, precision, accuracy, limits of detection (LOD) and limits of quantification (LOQ). In the fingerprint assessments, rather than using conventional qualitative "Similarity" as a criterion, the simple quantified ratio fingerprint method (SQRFM) was recommended, which has an important quantified fingerprint advantage over the "Similarity" approach. SQRFM qualitatively and quantitatively offers the scientific criteria for traditional Chinese medicines (TCM)/HM quality pyramid and warning gate in terms of three parameters. In order to combine the comprehensive characterization of multi-wavelength fingerprints, an integrated fingerprint assessment strategy based on information entropy was set up involving a super-information characteristic digitized parameter of fingerprints, which reveals the total entropy value and absolute information amount about the fingerprints and, thus, offers an excellent method for fingerprint integration. The correlation results between quantified fingerprints and quantitative determination of 5 marker compounds, including glycyrrhizic acid (GLY), liquiritin (LQ), isoliquiritigenin (ILG), isoliquiritin (ILQ) and isoliquiritin apioside (ILA), indicated that multi-component quantification could be replaced by quantified fingerprints. The Fenton reaction was Introduction According to statistics, approximately 10% of the world population suffers from gastric ulcers [1], which are caused or aggravated by etiological factors including infection of Helicobacter pylori, stress, frequent use of nonsteroidal anti-inflammatory drugs, smoking and alcohol consumption [2]. CBAT, a popular medicine used around the world, has been used for over 20 years in the treatment of such ulcers. CBAT includes not only three synthetic drugs, bismuth aluminate, heavy magnesium carbonate and sodium hydrogen carbonate, but also three medicinal herbs, Extract Licorice, Cortex Frangulae and Fructus Foeniculi (300, 25 and 10 mg in one tablet, respectively) [3]. Licorice, as a principal herb in CBAT, mainly contains triterpene saponins and flavonoids, which have been suggested to be main bioactive constituents [4,5]. Such major bioactive components as GLY, LQ, ILQ, ILG, ILA, etc. were often selected as marker compounds for the purpose of Licorice quality control [6][7][8].
Due to the mystery of chemical composition of TCM/HM that varies depending on a wide range of factors, such as the botanical species, the anatomical part of the plant used, the amount of sun the plant has been exposed to, the type of ground in which it is grown, harvesting time, cultivating region, and storage conditions [9], it is very difficult to maintain the quality consistency of TCM/HM at both manufacturer and batch level. In addition, TCM/HM are composed of hundreds of chemically different constituents, and multiple compounds often work synergistically in delivering therapeutic effects [10]. Consequently, a comprehensive and effective quality control method capable of guaranteeing their safety and efficacy in clinical applications is absolutely imperative. Now CBAT has been documented in the Chinese Pharmacopoeia (2010) and a quantitative standard with respect to its HM only involves GLY. Other reports concerning CBAT have mainly focused on analytical methods for bismuth [11][12][13][14][15] and have been confined to quantification of a single bioactive component [16,17]. The fingerprinting technique with an overall view of a complicated system has been proved to be of value in determining the identity, authenticity and batch-to-batch consistency of TCM/HM. Hence, such analytical profiles have been accepted and adopted as a strategy for TCM/HM quality assessment by the Food and Drug Administration (FDA) of the USA [18], European Medicines Agency (EMA) [19] and the State Food and Drug Administration (SFDA) of China [20,21]. Fingerprinting can be carried out using both chromatographic and spectroscopic techniques [22][23][24][25][26][27], and the chromatographic approach, especially HPLC, is regarded as one of the most potential and reliable means for TCM/HM quality control.
In the present study, we strategically established multi-wavelength chromatographic fingerprints and simultaneously determined 5 marker compounds to assess the HM quality in CBAT by HPLC-DAD, and the molecular structures of the quantification involved are shown in Fig. 1. In the fingerprint assessments, SQRFM, which has an outstanding quantified fingerprint advantage and exhibits a clear TCM/HM quality pyramid and warning gate using the combination of S 0 F , M F and α, was employed to assess the HM quality grades from a both qualitative and quantitative perspective. In order to combine the comprehensive characterization of multiwavelength ratio fingerprint profiles, an integrated fingerprint assessment method based on information entropy was set up. Because CBAT contains many components with antioxidant activities, correlation analysis between chromatographic fingerprint components and antioxidant activities of CBAT samples were performed and we obtained a robust calibration model by PLSR, which provided important medicinal efficacy information in vitro for CBAT quality control. It has been demonstrated that the present study offered an efficient way to monitor and evaluate the quality consistency of CBAT.

Theory of SQRFM
In the sample fingerprint vector (SFPV) x * ¼ ðx 1 ; x 2 ; . . . ; x n Þ and the reference fingerprint vector (RFPV) y * ¼ ðy 1 ; y 2 ; . . . ; y n Þ, x i and y i represent the peak area of the ith constituent in the sample and reference fingerprints, respectively. The 'Similarity', S F , between SFPV and RFPV can be expressed by Eq (1), as recommended by SFDA [28]. Moreover, the ratio similarity, S 0 F , in the ratio fingerprint profile as defined in Eq (2) can effectively cancel the effect of larger peaks masking smaller peaks [29]. A ratio fingerprint profile can be obtained by plotting x i /y i vs. the retention time, meanwhile, SFPV and RFPV are converted into X r * ¼ ðx 1 =y 1 ; x 2 =y 2 ; Á Á Á ; x n =y n Þ and Y r * ¼ ð1; 1; Á Á Á ; 1Þ, respectively. S 0 F not only clearly demonstrates the resemblance in the number and distribution of fingerprint peaks between SFPV and RFPV, but also takes into consideration the major and minor peaks in equal weights. The mean total content (M) of all the components in the ratio fingerprint profile can be calculated by Eq (3) and its corrected model, which can eliminate the error from cross-compensation among x i /y i values, is shown in Eq (4), where M F represents the corrected quantified ratio similarity. It should be noted that M F is also corrected for sample weight using a mass coefficient (f i ), which is defined as the ratio of the apparent weight of RFP (m R ) and the weight of the ith sample (m i ). The coefficient of variation (α) is defined in Eq (5) to reflect the dissimilarity between them. The approach to assess TCM/HM quality using the combination of S 0 F , M F and α in ratio fingerprint profiles is called SQRFM, by which TCM/HM are divided into 8 grades (Table 1). In the evaluation system, the TCM/HM quality pyramid and warning gate are clearly presented in Fig. 2. At the top of the pyramid is grade 1 which is the highest quality, so the proper criteria falling into grade 1-5 are certainly considered as qualified. The whole red region is termed the quality warning gate which reveals the characteristics of adulterants with S 0 F < 0.5, α > 0.50 and any M F .

Instruments
Chromatographic analysis was carried out on an Agilent 1100 HPLC series (Agilent Technologies, USA) equipped with a G1315B UV-vis DAD, a G1311A low pressure quaternary pump, a G1379A on-line degasser and a G1313A auto-sampler. System control and data analysis were carried out using an Agilent ChemStation workstation (Agilent, USA). Two columns, a Century SIL C 18 BDS (250 × 4.6 mm; 5.0 μm) (Century, China) and an Agilent poroshell 120SB C 18 (150 × 4.6 mm; 2.7 μm) (Agilent, USA), were used to optimize the separation conditions. The instrument for antioxidant activity assay was a 752-UV spectrophotometer (Gaomi Caihong, China).

Sample and Standard Solution Pretreatment
The tested samples were ground into powder with a mortar. Each powder sample (1.2 g) was accurately weighed and extracted three times with 25 ml of methanol for half an hour by refluxing. The combined extracts were concentrated in a rotary evaporator under vacuum, and then diluted to the mark with methanol in a 25 ml volumetric flask. The solutions were kept at 4˚C. Standard stock solutions of the 5 standards, ILQ, ILG, LQ, ILA and GHIA, were accurately weighed and then dissolved in methanol. A mixed standard solution was prepared by mixing the 5 individual standard stock solutions, and then six concentration levels of mixed standard solutions for the calibration curves were obtained by diluting the above mixed standard solution with methanol. These standard solutions were kept at 4˚C.
All the sample and standard solutions were filtered using a 0.22 μm nylon membrane prior to HPLC analysis.
All the sample fingerprint profiles were analyzed and evaluated by independently developed software 'Digitized Evaluation System for Super-Information Characteristics of TCM/HM Chromatographic Fingerprints 4.0' (Software certificate No. 0407573, China). In addition, SPSS 19.0, MATLAB 2009a, ORIGIN 8.5 and SIMCA 13.0 were also used for data processing and drawing.

Determination of Antioxidant Activity
We used the Fenton reaction [30] to determine the antioxidant activity of CBAT. Ferrous irons react with H 2 O 2 in an acidic medium to generate hydroxyl radicals (OH•) as shown in Eq (6); they possess outstanding oxidant ability and induce an advanced oxidation reaction, causing fading of the crystal violet. After adding the antioxidant to eliminate OH•, the degree of crystal violet fading was reduced. The scavenging ratio of OH• was indirectly deduced from the absorbance values of crystal violet.
The blank sample solution was prepared by adding successively 6.0 ml of 1.0 mmol/L FeSO 4 Á7H 2 O solution, 1.0 ml of 0.4 mmol/L crystal violet solution, and 0.09 ml of potassium hydrogen phthalate buffer solution (pH = 4.0) to a 25 ml volumetric flask and diluting to the mark with water, and then the absorbance (A 0 ) of the blank sample solution was measured at 588 nm. The reaction solution was prepared by adding 1.1 ml of 1% H 2 O 2 solution to the blank sample solution, then the absorbance (A b ) was measured.
The negative control samples were prepared by adding different volumes of methanol (0.10, 0.20, 0.25, 0.50 and 0.75 ml) to the blank sample solution, and the absorbance values (A B ) were measured. Finally the positive control samples were prepared by replacing methanol with corresponding volumes of sample solution as described in the section of 'Sample and Standard Solution Pretreatment'. These assay solutions were shaken and stored at 4˚C for 30 min, and then the absorbance values (A S ) were measured. The scavenging ratio (SR) of OH• was calculated using Eq (7):

Optimization of the Extraction and Chromatography Conditions
In order to obtain as much chemical information as possible and obtain as efficient a chromatographic separation as possible within a short analysis time, the extraction solvents, extraction times, analytical wavelengths, column types and mobile phases for CBAT samples were investigated. The index of the fingerprint information amount (I) [31] was applied as an objective function to optimize the extraction and chromatography conditions. The index I reflects the total signal strength, the signal homogenization and the information amount in a chromatogram, so the experimental conditions can be assessed by the I values of chromatograms and the higher the I value the better the condition. From S1 Fig., it was found that the I values of the methanol, 30 min, 250 nm, column 1 and mobile phases 2 were relatively higher than the other conditions, and thus a set of optimized experimental parameters was established as follows: methanol and 30 min were elected as the extraction solvent and extraction time, respectively, and the procedure 'Sample and Standard Solution Pretreatment' was adopted; 250 nm, column 1, i.e. the Century SIL C 18 BDS column (250 × 4.6 mm; 5 μm) as well as mobile phases 2, i.e. water-glacial acetic acid (A; 100:1, v/v) and acetonitrile-glacial acetic acid (B; 100:1, v/v) were used for chromatographic analysis. Because the maximum UV absorbance values of the marker compounds were observed at 250, 276, 330, 360 and 375 nm using a DAD detector, the above 5 wavelengths were selected for UV detection. Under the optimized extraction and chromatography conditions, ideal chromatograms were obtained which possessed not only abundant fingerprint information but also efficient separation characteristics including high signal responses, little interference and good peak shapes within 90 min (Fig. 3). Fig. 3 shows the typical chromatograms and ratio chromatograms of 27 batches of CBAT samples, as well as the chromatograms of the 5 marker standards. Chromatographic Fingerprint Analysis Methodology Validation of Fingerprint Analysis. All the sample solutions as described in the section of 'Sample and Standard Solution Pretreatment' were used to perform the following experiments. The repeatability of the method was assessed by analyzing six independently prepared samples (S1) using the same experimental procedure. The precision of the instrument was determined by loading the same sample solution (S1) six times consecutively. The stability of the sample was examined by analyzing the same sample solution (S1) after 0, 3, 6, 12 and 24 hours.
The relative retention time (RRT) and relative peak area (RPA) of each co-possessing fingerprint peak were calculated to estimate the repeatability, precision and stability, and the results showed that, for repeatability, the relative standard deviations (RSD) of RRT and RPA were less than 0.5 and 4.8%, respectively; for precision, the obtained values did not exceed 0.3 and 2.7%, respectively; for stability, the obtained values were less than 0.3 and 2.6%, respectively. Thus, the validation results demonstrated that the method satisfied the fingerprint analysis criteria for CBAT samples.
Evaluation of Multi-wavelength Fingerprints. Reference fingerprint profiles (RFP) in both chromatograms and ratio chromatograms (Fig. 3) were generated from the mean of 27 sample chromatograms at corresponding wavelengths. The fingerprint and ratio fingerprint profiles of CBAT samples at 5 wavelengths (250, 276, 330, 360 and 375 nm) contained 39, 31, 15, 15 and 9 co-possessing peaks, respectively. In order to carry out HM quality control, the chromatograms with comprehensive chemical information characterization were required. However, the obtained fingerprint and ratio fingerprint profiles clearly showed that it was impossible to achieve this at a single analytical wavelength because of the complexity of the HM composition. Although the most abundant fingerprint information was acquired and the largest responses of saponins were displayed at 250 nm, weaker responses of the remaining constituents were also observed. Therefore, an appropriately integrated fingerprint assessment method capable of synthesizing both abundant fingerprint information and intense signal responses at 5 wavelengths was absolutely essential.
In this study, a novel integrated fingerprint assessment method based on information entropy was set up. Information entropy is a super-information characteristic digitized parameter of fingerprints, which reveals the total entropy value and absolute information amount of a fingerprint system [32]. Integrated assessment procedures were carried out as follows. Firstly, sample fingerprint signals at 5 wavelengths were imported into the TCM/HM fingerprint software so as to acquire five sets of evaluation parameters (S 0 F , M F and α values, their RSD values were all below 3%, n = 3) of ratio fingerprints at 5 wavelengths (Table 2) and information entropy values (S j ) of ratio fingerprints at the 5 wavelength. Secondly, three assessment parameters were combined according to Eq (9)- (11), where e j represents the integrated weight at the jth wavelength and can be calculated based on Eq (8). S 0 Fj , M Fj and α j , as well as S 0 F , M 0 F and α 0 ( Table 2) represent three assessment parameters of the ratio fingerprints before and after integration, respectively. Finally, the integrated HM quality grades (Table 2) based on S 0 F , M 0 F and α 0 could be obtained according to TCM/HM quality criteria (Table 1) classified by SQRFM that could be used to comprehensively assess sample quality.
From Table 2, by observing the quality grades at 5 wavelengths and integrated wavelength, we found that: S1 and S13 had a constant quality grade at all wavelengths (grade 3 and 1, respectively); S17 and S25 had unqualified integrated grades (grade 7 and 8, respectively), while the remaining 25 samples had qualified integrated ones (grade 1-5); The integrated results exhibited some fluctuations and even greater differences compared with the quality grades at single wavelengths. For example, the integrated results of S3 and S22 were at medium level (grade 2 and 3, respectively) of grades (1, 2, 4 and 2, 3, 4, respectively) at single wavelengths. Interestingly, S17 and S25 exhibited unqualified integrated quality (grade 7 and 8, respectively) despite their having qualified grades (3, 5 and 4, respectively) at certain single wavelengths, conversely, S20 and S23 had qualified integrated quality (grade 3 and 5, respectively) in spite of their being unqualified (grade 6 and 6, respectively) at certain wavelengths. These grade changes might be attributed to differences in the fingerprint numbers and response strengths at different wavelengths, which illustrated that our evaluation strategy of combining multi-wavelength fingerprints was very comprehensive and essential to avoid any bias at a single wavelength.
In addition, the integrated S 0 F and α 0 values for 27 samples were all not below 0.89 and above 0.11, respectively, demonstrating that all samples had a similar chemical composition. Although they should be qualified and be in the range of grade 1-5 based on integrated S 0 F and α 0 from a qualitative perspective, in fact, S17 and S25 were judged as outliers (grade 7 and 8, respectively) in combination with integrated M 0 F (147.0 and 151.2%, respectively) from a quantitative perspective, indicating that qualitative evaluation was first performed (an acceptable S 0 F and α 0 should not be below 0.7 and above 0.3, respectively) and then further quantitative assessment (an acceptable M 0 F should be in the range of 70.0-130.0%) should not be ignored. Because the M 0 F can reflect the overall ingredient content in a sample, it has the potential to be associated with medicinal efficacy in clinics. In our opinion, SQRFM which includes both qualitative and quantitative assessments should be a better method for TCM/HM quality evaluation, and it should provide a vital measurement strategy for ensuring TCM/HM safety and efficacy.

Simultaneous Quantitative Analysis of the Five Marker Compounds
Methodology Validation of Quantitative Analysis. The linearity of the quantitative analysis method was assessed using a series of mixed standard solutions as described in the section of 'Sample and Standard Solution Pretreatment' under the optimized experimental conditions. The calibration curves were constructed by plotting the peak areas (y) vs. the injection masses (x, μg) of the 5 marker compounds. Their LOD and LOQ were determined by appropriately diluting the mixed standard solutions. The calibration curves, including the regression equations, correlation coefficients, linear ranges as well as LOD and LOQ, are displayed in Table 3. The correlation coefficients were all above 0.9996, indicating that the linear correlations of the 5 analytes were excellent between the peak area and injection mass over the tested ranges. According to the methodology validation experiments for fingerprint analysis, the peak areas of the 5 maker compounds were calculated for estimations of quantitative analysis repeatability, precision and stability, and the corresponding RSD values were found not to exceed 2.7, 1.2 and 1.6%, respectively. Recovery tests were performed using the standard addition method to evaluate the accuracy of the method and the mean recoveries of the 5 marker compounds were found to be between 96.0 and 104.0%. The above results demonstrated that the analysis method met quantitative requirements and was suitable for the simultaneous determination of ILQ, ILG, LQ, ILA and GLY in CBAT. Sample Analysis. The HM quality control in CBAT is carried out by determining the content of GLY in the Chinese Pharmacopoeia (2010). In fact, there are many other bioactive components besides GLY and, thus, the quantitative analysis of multiple components combined with chromatographic fingerprints was considered for quality control. The qualitative analysis of the 5 marker compounds in chromatograms were carried out by comparing the retention times and on-line UV spectra with those of marker standards to assign GLY, LQ, ILQ, ILG and ILA (Fig. 3). The developed HPLC method was applied to the simultaneous determination of the 5 marker compounds and the contents (mg/g) of the 5 marker compounds in 27 samples are summarized in S1 Table. From the obtained results, it was found that the contents of the 5 analytes varied markedly among the 27 samples, and this might be attributed to the raw herbs variability associated with a wide range of factors or variability in manufacturing processes.
Principal component analysis (PCA) is an important chemometrics method [33] capable of reducing dimensions for high dimensional variables and, thus, it is helpful for simplifying and analyzing complicated problems. In order to investigate the differentiating ability of marker contents, PCA was conducted using SIMCA software (version 13.0) and the contents of the 5 marker compounds were used as input data to construct two-dimensional matrices (27 × 5), with 27 and 5 representing the sample number and marker type, respectively. A two-component PCA model was obtained which cumulatively accounted for 94.7% of the variation. The total variance explained for the first principal component was 61.4% and that for the second principal component was 33.3%, which was visible at the left bottom of the PCA score plot (Fig. 4A). In the loading plot (Fig. 4B), the coordinate position of each marker reflected marker Table 3. Calibration plots, LOD and LOQ for the 5 compounds.

Compound
Calibration equation a y and x were, respectively, the peak areas and masses (μg) of the analytes. b LOD was defined as the mass for which signal-to-noise ratio was 3 and the LOQ was defined as the mass for which the signal-to-noise ratio was 10. weights in principal components and, so, the farther away from the coordinates-origin the marker position, the greater the correlation between the principal component and the marker. It was clear that GLY had the greatest correlativity with PC1, the same as LQ with PC2, and LQ and ILA had a greater correlativity with PC1, the same as GLY and ILA with PC2, while ILG and ILQ had less effect on PC1 and PC2. In the score plot, 27 batches of samples from the same manufacturer were distinctly divided into three clusters marked by group 1, 2 and 3, respectively. Three samples judged as unqualified (S17 and S25) and close to unqualified (S21) products by the integrated fingerprint assessment method fell into the same cluster (group 1), which was obviously different from the other two groups. In addition, the samples in group 2 had a positive correction with GLY, and the GLY contents (in the range of 19.992-23.387 mg/g) in group 2 were all higher than those (in the range of 10.377-19.65 mg/g) in group 3. Also, the samples in group 1 exhibited a positive correlation with LQ and ILA, hence the contents of LQ and ILA (in the ranges of 10.354-12.877 mg/g LQ and 3.669-4.839 mg/g ILA, respectively) in group 1 were all higher than those (in the ranges of 6.225-9.022 mg/g LQ and 3.639-1.599 mg/g ILA, respectively) in other two groups of samples. Consequently, it was found that the three products with the worst grades were clustered in the same group (group 1), mainly because of their high contents of LQ and ILA; the difference between group 2 and group 3 was mainly attributed to the contents of GLY. Therefore, we were able to conclude that the obtained two-component PCA model displayed better discriminating ability among the 27 samples from the same manufacturer.
Correlation between Corrected Quantified Ratio Similarity and Quantitative Analysis. In order to explore the relationship between the corrected quantified ratio similarity M F values of ratio fingerprints and quantitative results of 5 marker compounds, data processing and analysis were performed as follows. Firstly, the ith marker compound contents (z i ) in 27 samples were converted into percentage values (P i ) of the marker contents out of its average content ( z i ) of 27 samples as defined in Eq (12), and it was noteworthy that each sample weight should be corrected using a mass coefficient (f i ) when calculating P i ; Secondly, the mean of P i values of the 5 marker compounds in each sample (i.e. P 5C , Table 2) was calculated according to Eq (13), which was termed the logarithm-exponent mean.
Finally, the P 5C values of 27 samples (y) were, respectively, plotted vs. the M F values (x) at 6 wavelengths including integrated wavelength, 250, 276, 330, 360 and 375 nm (Fig. 5). From Fig. 5, it was apparent that there were better linear correlations between P 5C and M F values at 6 wavelengths, because their correlation coefficients (R) were successively 0.9224, 0.8552, 0.8565, 0.8927, 0.7871 and 0.9227, which were all above 0.8500 except for 0.7871 at 360 nm and, especially at an integrated wavelength, reached a more satisfactory value of 0.9224. This indicated that the M F values at the 6 wavelengths were basically consistent with the logarithm-exponent mean of the quantitative results for the 5 marker compounds. We concluded that multi-component quantification can be substituted by quantified fingerprint analysis for the purpose of TCM/HM quality control, which can be potentially applied to TCM/HM practical production due to its simplicity and economic advantage.

Correlation Analysis between Quantified Fingerprints and Antioxidant Activities
CBAT contains many components with antioxidant activities, and the sample antioxidant capacity should take account of all the antioxidant components. In the present paper, the EC 50 values of CBAT samples were used as the measurement parameter of antioxidant activity for 27 samples (Table 4). In order to explore the relationship between the quantified fingerprints and the antioxidant activities, attempts were made to correlate the EC 50 to HPLC fingerprints at 250 nm and, as shown in S2 Fig., there were 28 out of 39 peaks in the fingerprints exhibiting negative correlation and 11 peaks had positive correlation with the EC 50 . PLSR, a well-known multivariate regression method [34,35], was also applied to explore the spectrum-effect relationship by taking into consideration the EC 50 and the peak areas of 39 co-possessing peaks at 250 nm as two groups of variables. The regression characteristics were used to adjust a model to the measured data with the purpose of quantifying the relationship between two groups of variables. The appropriately adjusted model could be used to describe the relationship between the two variables and to predict new variables. The PLSR model was validated by means of full cross-validation for the purpose of rationality. From the relative errors (Table 4) of the predicted values from full cross-validation, two samples (S17 and S25) were identified as singular points, and thus they were eliminated when final mathematical model was constructed. Interestingly, two singular points, i.e. S17 and S25, were unqualified samples (grade 7 and 8, respectively) judged by the integrated fingerprint assessment method, indicating that the integrated quality grades on the basis of SQRFM were in agreement with the antioxidant activities of samples in vitro, while SQRFM as an TCM/HM fingerprint evaluation method and integrated multi-wavelength fingerprints as an assessment strategy were found to be both scientific and powerful. The obtained calibration model using the PLSR method is expressed by Eq (14) EC 50 ¼ 0:9621 À 0:0005x 1 þ 0:00112x 4 À 0:00115x 5 þ 0:00054x 6 À 0:00199x 7 À0:00016x 12 À 0:00018x 15 À 0:00013x 16 þ 0:00042x 17 À 0:00015x 18 þ 0:00013x 20 þ0:00016x 21 À 0:00019x 22 À 0:0009x 23 þ 0:00115x 24 À 0:00059x 27 À 0:00046x 28 À0:00021x 29 À 0:00188x 30 þ 0:00044x 34 À 0:00019x 35 þ 0:00039x 38 The above model was applied to predict the antioxidant activities of 27 batches of samples. The predicted EC 50 values (Table 4) could be calculated using the areas of 22 peaks in the fingerprints correlated with EC 50 , and the absolute values of the relative errors of the predicted results (Table 4) were all less than 5% except for two singular samples (S17 and S25), indicating that the calibration model was robust. The calibration curve was also constructed by plotting the predicted EC 50 (y) vs. the measured EC 50 (x) of 25 batches of CBAT samples (except for S17 and S25), as shown in Fig. 6. We found that the calibration curve had a slope of 0.93818 (standard error, S = 0.04974) and an intercept of 0.04148 (S = 0.03382) at the left top of the plot, further indicating that the predicted EC 50 values were coincident with the measured ones and the calibration model had an excellent predictive ability. Therefore, correlation analysis between the quantified fingerprints and the antioxidant activities can provide important medicinal efficacy information in vitro for HM quality control in CBAT.

Conclusions
In the present study, a combination method of integrating multi-wavelength fingerprints and quantifying multiple bioactive components using HPLC-DAD was successfully developed and applied to 27 CBAT samples for the purpose of HM quality control. SRQFM for TCM/HM quality assessment was recommended, because it contains both qualitative and quantitative similarity evaluations and, thus, can overcome the deficiency of a single qualitative criterion and reflect the genuine quality of TCM/HM. The quality of 27 CBAT samples from the same manufacturer were well differentiated based on the integrated fingerprint assessment method of SQRFM, and two samples (S17 and S25) were judged as unqualified products (grade 7 and 8, respectively) and the remaining 25 samples were in the range grade 1-5. Moreover the correlation analysis between fingerprints and antioxidant activities was carried out using PLSR and the established model had excellent predictive ability, providing important information for CBAT quality control. This study offers a scientific, sensitive and comprehensive analytical strategy for the quality control of CBAT practical production.
Supporting Information