Abstract

Yangyin Tongnao (YYTN) prescription is used as a traditional Chinese herbal formula, and it has antioxidant activity that mainly contributes in the treatment of cardiovascular and cerebrovascular diseases. However, the compounds related to its antioxidant activity are still unknown. In the present study, the fingerprints of YYTN extracts under different extraction conditions were obtained by high performance liquid chromatography (HPLC) to identify the common peaks to all the samples processed. A 2,2-diphenyl-1-picrylhydrazyl (DPPH) radical scavenging activity assay and ferric reducing antioxidant power (FRAP) assay were carried out to evaluate the antioxidant activity of the extracts. Spectrum-effect relationship between HPLC fingerprints and antioxidant activity of YYTN was assessed by Pearson product-moment correlation coefficient (PPMCC) and multiple linear regression analysis (MLRA). The results showed that peaks 5, 6, 13, 15, and 24 of the fingerprints were closely connected to antioxidant activity. Five peaks were identified: vanillic acid (P5), puerarin (P7), ferulic acid (P13), daidzein (P21), and formononetin (P23). Our study successfully established the spectrum-effect relationship between HPLC fingerprints and antioxidant activity of YYTN, which provided a general method for establishing quality standards with a combination of chromatography and antioxidant activity.

1. Introduction

Yangyin Tongnao prescription is a traditional Chinese herbal formula composed of six kinds of Chinese herbal medicines, Radix Rehmanniae, Radix Astragali, Puerariae Lobatae Radix, Ligustici Chuanxiong Rhizoma, Dendrobii Caulis, and Hirudo. However, four kinds of these Chinese herbal medicines except Dendrobii Caulis and Hirudo have been studied as a simplified formula in a lot of research [1, 2]. The proportions of herbal composition of Yangying Tongnao (YYTN) simplified prescription are shown in Table 1 [3]. The efficacy of YYTN has been evaluated by clinical trials and animal experiments. According to the Chinese medicine, it has the effects of nourishing yin and replenishing qi (deficiency of yin or qi is a kind of syndrome; deficiency of yin usually appears hot and bothered, and deficiency of qi has a feeling of weariness). In the Western medicine, it could promote blood circulation, especially in the treatment of cardiovascular and cerebrovascular diseases [4]. Previous researches have proved that YYTN had the functions of inhibiting neuronal apoptosis, improving haematological disorders [5], and is related with antioxidant and anti-inflammatory properties [6]. As a traditional Chinese herbal formula, YYTN can also improve the neurological deficits after focal cerebral ischemic infarction in rats [3]. At the same time, it was found that YYTN has protective effect on oxidative stress-induced cerebral ischemia injury in vivo through animal experiments [7]. However, there does not seem to be any report about the relationship between active ingredients and pharmaceutical efficacy of YYTN.

Antioxidants play a crucial role in preventing various serious diseases caused by oxidative stress damage. At present, several methods have been applied to evaluate the antioxidant activity, including assays for ferric reducing antioxidant power (FRAP) [810], 2, 2-diphenyl-1-picrylhydrazyl (DPPH) radical scavenging [1113], and hydroxyl radical scavenging [14, 15]. Therefore, a growing number of researchers are paying attention to the search for natural compounds from plants or vegetal drugs with antioxidant properties [1618]. However, the compounds of YYTN responsible for its antioxidant activity are unclear.

Nowadays, spectrum-effect is used to study the relationship between components and efficacy of traditional Chinese medicines. For example, the main vasorelaxant components of rosemary have been found by using the chemical spectrum-effect relationship [19]. HPLC fingerprint and bioactivity verification were used to investigate the potential correlation between the chemical components and pharmaceutical effects [20, 21]. It is a representative method to find the main active components in complex traditional Chinese herbal formulas based on the study of spectrum-activity relationship.

There are three major categories of data processing methods for spectrum-effect relationships, such as correlation analysis, cluster analysis, and regression analysis. Correlation analysis is a statistical method for studying the degree of closeness between variables [22], including grey relational analysis (GRA) and Pearson product-moment correlation coefficient (PPMCC). Cluster analysis is one of the important techniques in data mining and exploratory data analysis [23], including hierarchical clustering analysis (HCA) and non-hierarchical clustering analysis (N-HCA). Regression analysis refers to a statistical analysis method to determine the quantitative relationship between two or more variables [24], including multiple linear regression analysis (MLRA) and partial least squares (PLS). Therefore, GRA, PPMCC, PLS, and MLRA are applied to study the spectrum-effect relationship in the correlation analysis between the components and the efficacy of Chinese formulas [2529].

In our study, HPLC fingerprints were combined with the data of the 2,2-diphenyl-1-picrylhydrazyl (DPPH) radical scavenging and ferric reducing antioxidant power (FRAP) assays to study the spectrum-effect relationship of YYTN using PPMC and MLRA. It lays the foundation for the pharmacokinetics substantial basis of YYTN.

2. Materials and Methods

2.1. Materials and Reagents

Radix Rehmanniae, Radix Astragali, Puerariae Lobatae Radix, and Chuanxiong Rhizoma were purchased from Traditional Chinese Medicine Pieces Co., Ltd. of Zhejiang Chinese Medical University (Hangzhou, China) and identified by Professor Shengwu Huang, Zhejiang Chinese Medical University (Hangzhou, China). HPLC-grade methanol and acetonitrile were obtained from Tedia Company, Inc. (Fairfield, OH, USA). Water was purified with a Millipore purification system (Millipore Co., Ltd, Billerica, MA, USA). DPPH was purchased from Shanghai Biochemical Co., Ltd. (Shanghai, China). FRAP kit was purchased from Nanjing Institute of Bioengineering (Nanjing, China). The reference standards of vanillic acid (≥99.0%), puerarin (≥98.0%), ferulic acid (≥98.0%), daidzein (≥99.0%), and formononetin (≥98.0%) were purchased from Shanghai Yuanye Co., Ltd. (Shanghai, China).

2.2. HPLC Fingerprints
2.2.1. Chromatographic Conditions

HPLC analysis was performed using a Waters 2695 HPLC equipped with 2498 UV detector. The chromatographic separation was performed on a Hypersil BDS C18 column (4.6 mm × 300 mm, 5 μm) maintained at 25°C. The mobile phase consisted of acetonitrile (A) and 0.5% formic acid water (B) with the following optimized gradient elution: 0 min: 2% A; 10 min: 10% A; 29 min: 12% A; 35 min: 15% A; 50 min: 25% A; 65 min: 60% A; 70 min: 90% A. The flow rate was 1.0 mL/min with an injection volume of 10 μL, and the detection wavelength was set at 280 nm.

2.2.2. Preparation of Standard Solutions

The mixed reference standard solution was prepared by adding 3 mg of vanillic acid, 36 mg of puerarin, 2.4 mg of ferulic acid, 6 mg of daidzein, and 1.8 mg of formononetin in a 10 mL volumetric flask and dissolved in 50% methanol to obtain final concentrations of 0.3, 3.6, 0.24, 0.6, and 0.18 mg/mL, respectively. Then the solution was filtered through a 0.22 µm microporous membrane, and 50% methanol was added to final concentrations of 0.15, 0.075, 0.05, 0.0375, and 0.03 mg/mL of vanillic acid; 1.8, 0.9, 0.6, 0.45, and 0.36 mg/mL of puerarin; 0.12, 0.06, 0.04, 0.03, and 0.024 mg/mL of ferulic acid; 0.3, 0.15, 0.1, 0.075, and 0.06 mg/mL of daidzein; 0.09, 0.045, 0.03, 0.0225, and 0.018 mg/mL of formononetin for HPLC analysis.

2.2.3. Preparation of Sample Solutions

The four herbs were crushed into crude grains. Accurately weighed 1.00 g of Radix Rehmanniae, 1.00 g of Radix Astragali, 1.20 g of Puerariae Lobatae Radix, and 0.67 g of Chuanxiong Rhizoma crude grains were mixed, and aqueous extracts were prepared using ultrasound-assisted extraction according to Tables 2 and 3. All the solutions were filtered through a 0.22 µm microporous membrane and stored at 4°C until analysis.

2.2.4. Validation of Methodology

Method precision was evaluated by five replicate injections of one sample, while repeatability was evaluated by five samples prepared independently from the same initial conditions. For the stability study, sample solutions were stored at room temperature and analyzed at 0, 2, 4, 8, 12, and 24 h, respectively.

2.2.5. Peak Identification and Evaluation of Fingerprints

The reference standard solutions were injected into the HPLC system, and the chromatographic peaks in the sample solutions were identified according to the retention time of the reference substances.

Each sample solution and each standard solution was analyzed with the HPLC optimized condition to obtain the HPLC fingerprints. The fingerprints were matched automatically by the Similarity Evaluation System for Chromatographic Fingerprint of Traditional Chinese Medicine [30] (Version 2012A). Then the reference fingerprint was formed by the system using the median method from the comparison of nine samples, and the similarity values between the reference and nine samples were calculated using this software.

2.3. Antioxidant Activity Determination
2.3.1. DPPH Assay

DPPH radical scavenging activity was determined according to Zhou et al. [31] with slight modifications. Briefly, 100 µL of 0.4 mM DPPH anhydrous ethanol solution and 100 µL of each sample solution at various concentrations were added to the 96 microplate. After incubation for 30 min at room temperature in the dark, the absorbance of the mixture was measured at 517 nm (SpectraMax Plus384 Absorbance Microplate Reader, Molecular Devices, USA). DPPH radical scavenging activity was calculated according to the following equation:where Ai was the absorbance of sample + DPPH; Aj was the absorbance of sample + anhydrous ethanol; A0 was the absorbance of anhydrous ethanol + DPPH.

2.3.2. FRAP Assay

The total antioxidant activity was evaluated by FRAP assay [32]. The FRAP solution was prepared by mixing 1.50 mL of 300 mmol/L acetate buffer (pH = 3.6) with 0.15 mL of 10 mmol/L TPTZ solution in 40 mmol/L HCl and 0.15 mL of 20 mmol/L FeCl3 solution, incubated for 30 min at 37°C in the dark. A fresh 100 mM FeSO4 solution was prepared from FeSO4.7H2O standard material. Then 5 μL of the sample and 180 μL of FRAP solution were added to the 96 microplate and allowed to react in the dark for 6 min at room temperature, and the absorbance was measured at 593 nm (SpectraMax Plus384 Absorbance Microplate Reader, Molecular Devices, USA). Aqueous standard solutions of 0.15, 0.3, 0.6, 0.9, 1.2, and 1.5 mmol/L of FeSO4 were used to generate a calibration curve. The final results were expressed as the equivalent concentration of the FeSO4 standard solution.

2.4. Spectrum-Effect Relationship
2.4.1. Pearson Product-Moment Correlation Coefficient

The PPMCC is a statistical method to study the degree of closeness and direction of change between the two variables [33]. In this study, PPMCC was used to calculate the relationship between the peak area values of twenty-four common peaks to all the samples (P1 to P24) in HPLC fingerprints and the antioxidant activities.

2.4.2. Multiple Linear Regression Analysis

To eliminate the problem of multicollinearity that may exist among multiple independent variables, this study used the stepwise regression module in multiple linear regression in SPSS statistical software [34] to analyze the peaks of HPLC fingerprints, IC50 values of DPPH radical scavenging percentage, and the ferric reducing antioxidant power, respectively. Therefore, characteristic peaks with significant correlation among multiple independent variables were selected.

2.5. Statistical Analysis

The IC50 was calculated using the following URL: https://www.aatbio.com/tools/ic50-calculator. MLRA and PPMCC were performed by SPSS 25.0 software.

3. Results

3.1. Results of HPLC Fingerprints
3.1.1. Method Validation

Method validation for HPLC fingerprint results showed that precision of the same sample solution was within the range of 0.19–1.37% for relative retention time and 0.75–2.67% for the average peak area of the common peaks to all the samples. Repeatability of the experiment appeared within the range of 0.15–1.04% for the relative retention time and 0.57–2.48% for relative peak area of common peaks. Sample stability was 0.11–1.32% for relative retention time and 0.69–2.74% for relative peak area of the common peaks to all the samples. All results indicated that the method of HPLC analysis was valid and satisfactory.

3.1.2. HPLC Fingerprints

The HPLC fingerprints of the extracts of YYTN and their reference fingerprints are shown in Figure 1, and sample and mixed reference solution are shown in Figures 2 and 3, respectively. Twenty-four peaks with good separation and large peak areas were selected as common peaks to all the samples. The average retention time and peak areas of samples and the coefficient of variation (CV) of twenty-four peak areas are listed in Table 4. Values of CV were all more than 8.96%, which indicated that the content of each component of the samples under different extraction conditions varies greatly, especially the unknown compound represented by peak 2.

3.1.3. Quantitative Analysis of YYTN

After comparing with the reference standard, five of the twenty-four common peaks were identified, which were vanillic acid (P5), puerarin (P7), ferulic acid (P13), daidzein (P21), and formononetin (P23). As shown in Table 5, quality control was performed on the contents of five compounds in nine samples.

3.1.4. Similarity Analysis of Fingerprints

Each sample fingerprint of YYTN was analyzed by comparing with the reference fingerprint to evaluate the similarity among these samples. Results of similarity values between the reference and the nine samples are shown in Table 6.

3.2. Results of the Antioxidant Activity Test
3.2.1. DPPH Assay

The IC50 values of DPPH radical scavenging activity test are listed in Table 7. The results revealed that nine extracts of YYTN solution exhibited the IC50 values that appeared within the range of 3.023–4.396 mg/mL. As shown in Table 7, S5 exhibited the most DPPH radical scavenging activity.

3.2.2. FRAP Assay

According to the FRAP kit, the absorbance value of the FeSO4 standard was plotted on the abscissa, and the concentration corresponding to each absorbance value was plotted on the ordinate. The regression equation was as follows: Y = 3.4887X − 0.0291, R2 = 0.9994. The results revealed that the nine extracts of YYTN solution varied from 5.331 to 8.245 mmol/L equivalent to FeSO4. As shown in Table 7, S9 exhibited the highest ferric reducing antioxidant power.

3.3. Spectrum-Effect Relationship Results
3.3.1. Pearson Product-Moment Correlation Coefficient

PPMCC analysis is performed with peak areas of twenty-four peaks as independent variables and IC50 values of DPPH radical scavenging percentage and FRAP values as dependent variables. Results of PPMCC are shown in Table 8. The PPMCC represents the correlation between the peak areas and the IC50 values of radical scavenging percentage and the FRAP values. If the PPMCC value was positive, it indicates that there was a positive correlation between independent and dependent variables. Conversely, if the PPMCC value is negative, then there is a negative correlation. The larger the absolute value of PPMCC, the better the correlation.

3.3.2. Multiple Linear Regression Analysis

Using SPSS software, the peak areas of twenty-four common peaks were selected as independent variables and expressed as P1P24, and IC50 values of DPPH radical scavenging percentage and FRAP values were selected as dependent variables and expressed as Y1 and Y2, respectively. By stepwise regression analysis, the multiple regression equation of IC50 values of DPPH radical scavenging percentage variable was expressed as Y1 = 0.308 − 0.004 × P5 − 0.001 × P6 + 0.001 × P9 − 0.002 × P13 0.004 × P15 + 0.006 × P16 − 0.001 × P21 + 0.011 × P24, and the equation of FRAP was expressed as Y2= −1.742 + 0.011 × P5 + 0.008 × P6 × P9 + 0.001 × P13 − 0.026 × P15 + 0.024 × P16 − 0.001 × P21 + 0.044 × P24. These two multivariate linear regression equations represent peaks that were strongly correlated with IC50 values of DPPH radical scavenging percentage and FRAP values. In the regression process, some peaks with no obvious correlation and multicollinearity were eliminated. Finally, the correlation peak was chosen to get the regression equation.

4. Discussion

YYTN is composed of Radix Rehmanniae, Radix Astragali, Puerariae Lobatae Radix, and Chuanxiong Rhizoma. Radix Rehmanniae had been found to have anticancer effects through promoting the activation of natural killer cells and dendritic cells [35, 36]. Radix Astragali has antitumor activity in vitro and in vivo [37]. It also exerts pharmacological activities on the cardiovascular, immune, respiratory, and hepatic systems [38]. Puerariae Lobatae Radix has been shown to be effective in the treatment of hypertension, cerebral ischemia, and diabetes [3941]. Chuanxiong Rhizoma has a protective effect on diabetic nephropathy and has anticancer and antioxidant pharmacological activities [42, 43]. Therefore, YYTN, as the traditional Chinese herbal formula that combines these four Chinese medicines, has been extensively applied in the treatment of brain injury [3]. However, there are few reports on the HPLC fingerprint and pharmacological efficacy of YYTN.

From Figure 1 and Table 4, there were certain differences in the peak areas between the nine samples under different extraction conditions. However, as shown in Table 6, the similarity value of nine sample fingerprints and reference fingerprints ranged from 0.907 to 0.999, which indicated that the quality of the YYTN composed of Radix Rehmanniae, Radix Astragali, Puerariae Lobatae Radix, and Chuanxiong Rhizoma was relatively consistent and stable. We identified five compounds in YYTN by comparing Figures 2 and 3: P5 was vanillic acid, P7 was puerarin, P13 was ferulic acid, P21 was daidzein, and P23 was formononetin. Subsequently, we carried out the quantitative determination of these compounds’ content. The results are shown in Table 5: the contents of known components in YYTN were puerarin > formononetin > daidzein > vanillic acid > ferulic acid, indicating that puerarin was the most abundant in YYTN. The result of this experiment was consistent with that of previous research of our laboratory.

As shown in Table 7, DPPH radical scavenging percentage and FRAP values can reflect the antioxidant activity [44]. The results showed that YYTN had a good antioxidant activity. The literature showed that the antioxidant activity increased with the decrease of IC50 of DPPH radical scavenging percentage and the increase of FRAP values [45, 46]. According to the PPMCC shown in Table 8, P4, P6, P7, P8, P9, P10, P12, P16, P21, P23, and P24 are positively correlated with DPPH scavenging percentage and FRAP values. P6, P10, and P12 were negatively correlated with IC50 of DPPH radical scavenging percentage, while these peaks were positively correlated with FRAP. The order of correlation between peaks and IC50 of DPPH radical scavenging percentage followed P12 > P10 > P6, and FRAP values followed P6 > P12 > P10. MLRA could eliminate the problem of multicollinearity that may exist among multiple independent variables. Therefore, MLRA analysis showed fewer peaks related to antioxidant activity than PPMCC. By MLRA, P5, P6, P9, P13, P15, P16, P21, and P24 were related to IC50 of DPPH radical scavenging percentage and FRAP values. Among them, P5, P6, P13, P15, and P24 were negatively correlated with IC50 of DPPH radical scavenging percentage, and P5, P6, P13, and P24 were positively correlated with FRAP, which meant that P5, P6, P13, P15, and P24 may have significant antioxidant activity. P5 and P13 were identified by HPLC as vanillic acid and ferulic acid, but P6, P15, and P24 were not identified in this study. Therefore, we will further identify these three unknown components and prove whether these five peaks have a greater contribution to the antioxidant activity.

Because of the complexity of the components of traditional Chinese medicines, the chemical fingerprint cannot accurately evaluate the related information between the components and their efficacy [47, 48]. In this study, the fingerprints of YYTN were established, and no related reports were found in the literature. It was concluded that this study laid a foundation for the establishment of the fingerprint of YYTN. Subsequently, the antioxidant activity of YYTN was tested, which showed that YYTN had a certain extent of antioxidant activity and provided a certain therapeutic idea for the treatment of cardiovascular and cerebrovascular diseases because YYTN has been shown to be effective for cerebral ischemia and cardiac ischemia in animal experiments [49, 50]. Finally, the components of YYTN related to antioxidant activity were determined by correlation analysis, which may be significant to the research of the therapeutic material basis of YYTN.

Mathematical modeling is used in spectrum-effect relationship of traditional Chinese medicines [34]. For example, the two methods (PPMCC and MLRA) we used in this work indicated that both could establish the spectral correlation model to determine the chromatographic peaks which are closely related to the antioxidant activity. The two antioxidant methods mentioned in this work were also used in a lot of research [5153]. The results showed that the methods were reliable for the determination of antioxidant activity in vitro. Of course, we will verify this result further in both cell and animal experiments in vivo. Our study not only established the spectrum-effect relationship, but also quantitative analysis of the active components in YYTN compared to other traditional Chinese medicine spectrum-effect correlations [20, 34].

In general, the main advantage of this research is that the relationship between HPLC fingerprints and the antioxidant activity of YYTN was established. Main antioxidant components of YYTN are further determined and would contribute to the quality standards of YYTN and to the study of the mechanism of other antioxidant components.

5. Conclusion

In the present study, we established the spectrum-effect relationship between the HPLC fingerprints and scavenging activity for DPPH and the total antioxidant activity for FRAP of YYTN. The HPLC analysis was used to build fingerprints of YYTN which contained up to twenty-four common peaks, and the similarity values of these fingerprints were evaluated by similarity analysis. The results showed that the similarity values of nine samples were more than 0.907, and puerarin was the most abundant chemical constituent in YYTN. The results of the spectrum-effect relationship indicate that peaks 5, 6, 13, 15, and 24 may be the main components responsible for the antioxidant activity of YYTN. This study shows that the fingerprinting method has been validated; it could provide a reliable and practical method for the consistency of quality of traditional Chinese medicine and other herbal preparations.

Data Availability

The data used to support the findings of this study are available from the corresponding author upon request.

Conflicts of Interest

The authors declare no conflicts of interest.

Authors’ Contributions

Yangyang Zhang and Li Yu contributed equally to this work. Yangyang Zhang performed the experiments. Li Yu designed the experiments and drafted the manuscript. Xixi Zhao analyzed the data. Haofang Wan and Yu He helped to revise the manuscript. Haitong Wan and Jiehong Yang conceived the study and revised the manuscript.

Acknowledgments

This research was funded by the National Natural Science Foundation of China (Grant nos. 81904083, 81973560, and 81630105), key projects of the Natural Science Foundation of Zhejiang Province (Grant no. LZ18H270001), and the Foundation of Zhejiang Chinese Medical University (Grant no. 2018ZJ11).