Zhen-jie Liu‡
abc,
Zhi-long Shi‡ab,
Can Tuab,
Hai-zhu Zhangab,
Dan Gaob,
Chun-yu Lib,
Qin Heb,
Rui-sheng Lib,
Yu-ming Guob,
Ming Niub,
Cong-en Zhangb,
Yong-shen Reng,
Han-shen Zhen*af,
Jia-bo Wang*be and
Xiao-he Xiao*dg
aCollege of Pharmacy, Chengdu University of Traditional Chinese Medicine, Chengdu 610075, China. E-mail: 8zhen@163.com
bChina Military Institute of Chinese Medicine, 302 Military Hospital, Beijing 100039, China. E-mail: pharm_sci@126.com
cGuangxi Botanical Garden of Medicinal Plants, Nanning 530023, China
dIntegrative Medical Center, 302 Military Hospital, Beijing 100039, China. E-mail: pharmacy302xxh@126.com; 8zhen@163.com
eState Key Laboratory Breeding Base of Systematic Research, Development and Utilization of Chinese Medicine Resources, Chengdu University of Traditional Chinese Medicine, Chengdu 610075, China
fCollege of Pharmacy, Guangxi University of Traditional Chinese Medicine, Nanning 530001, China
gSchool of Pharmaceutical Sciences, South-Central University for Nationalities, Wuhan, Hubei 430074, China
First published on 17th January 2017
Due to the inherent nature of the heterogeneity and complexity in the chemical composition of medicinal plants, there are significant challenges in the characterization and quantification of all active chemical constituents of botanicals, such as those used in traditional Chinese medicines (TCMs) and botanical drugs (BDs). The selection of one or even multiple marker chemical compounds for quality control may often be inadequate, as those marker compounds may not be associated with or fully represent the clinical effects of the botanical products. In this study, we proposed a novel activity-associated chemical analysis approach, named the effect-constituents index (ECI), to evaluate the quality of TCMs, using Salvia miltiorrhiza Bge. (SM), a commonly used TCM herb, as an example. The ECI is not only a quantitative analysis of multiple chemical constituents, but also includes an effect-based calibration of the relative bioactivity coefficient of each chemical constituent to express “the whole biological/pharmacological effect” of a TCM. Our chemical assay results showed significant variations in the contents (Ci) of nine main constituents, including cryptotanshinone and salvianolic acid B, in the SM samples. In addition, neither the concentrations of the individual constituents nor the sum of those constituents in the tested SM samples had good associations (r < 0.81) with the biopotency of SM itself, based on the activity data from an antiplatelet aggregation test. To calculate the ECI, cryptotanshinone and the other five constituents of high biopotency (above 660 U mmol−1) were selected, and their biopotency values were set as the calibration weight of each constituent; then, the sum of the products of content (Ci) and biopotency weight (Fi) were calculated as the ECI. By regression analysis, the ECI had either the highest correlation coefficient (r = 0.92, p < 0.001) with the biopotency of SM or the lowest residual compared to the chemical contents. Moreover, the ECI showed a good ability to distinguish and predict the biopotency-based quality grade, while the chemical markers alone did not. In conclusion, our study indicates that this effect-calibrated quantitative approach is a useful tool to associate the “quality” with the potential clinical effects of TCMs and BDs consisting of multiple active constituents or known chemical marker compounds.
Our group had established several biological assay methods for the quality control of TCMs, such as purgative activity assays, cardiac activity assays, antiviral activity assays, and antiplatelet aggregation assays.8–18 Despite the potential advantage of being able to show close association with clinical efficacy, biological assays still have some limitations, such as limited operability, high cost, low precision, and animal ethics issues, when live animals are used in some instances.2,19–22 However, those limitations are usually not limitations for chemical analyses. Hence, it might be attractive if we can combine the two types of approaches, chemical determination and biological assays, to form a novel methodology with added ability to control the quality of TCMs and BDs and demonstrate consistency.
As a result, we propose a novel efficacy-associated approach, named the effect-constituents index (ECI), to evaluate the quality of TCMs.2 The ECI is not only a determination of multiple chemical constituents but also an effect-based calibration using the relative bioactivity coefficient of each chemical constituent to express the whole effect of a TCM. To form the ECI, the chemical constituents of high biopotency are selected, and the biopotency values are set as the calibration weight of each constituent; then, the sum of the products of constituent content (Ci) and the biopotency weight (Fi) is calculated as the ECI, which equates to the “total bioactivity” of those constituents in a TCM or BD. Thus, the ECI would have a good association with clinical efficacy similar to a biological assay. Moreover, when the ECI has been constructed for a specific TCM or BD, the only instrumental analysis needed is the chemical determination of given constituents.
ECI was only clarified in theory, but it was not validated in practice. To test the proposed ECI approach, we chose the commonly used TCM Salvia miltiorrhiza Bge. (SM or Danshen in Chinese) as an example. SM is widely used in Asia, the United States and European countries for the prevention and treatment of cardiovascular and cerebrovascular diseases.23,24 The beneficial actions are attributable to anti-coagulant, anti-thrombotic, and especially anti-platelet aggregation effects.25,26 Platelet aggregation is one of the essential contributors to the onset and development of cardiovascular diseases. Some constituents isolated from SM, such as salvianolic acid B and cryptotanshinone, were proven to have antiplatelet aggregation activity.27–30 At present, quality control for SM is performed routinely with multi-constituent determination by high performance liquid chromatography, as stipulated in the China Pharmacopeia (ChP), United States Pharmacopeia (USP), and European Pharmacopeia (EP). For example, four chemicals, salvianolic acid B, tanshinone, tanshinone IIA, and cryptotanshinone, were determined as quality markers of SM in the ChP. However, we found that the sum of those four constituents could not distinguish the bioactivity-based quality grade of SM, indicating its limited ability to control the real quality associated with clinical efficacy. Therefore, there is an obvious need to establish a new quality control method for SM. The experimental flow chart of this study is shown in Fig. 1.
Clustering analysis was a statistic method used for determining the natural structure of object in multidimensional profiles. The clustering analysis can divide objects or individuals into specific clusters in a way that maximizes homogeneity within each cluster while also maximizing heterogeneity between clusters.31,32 We introduce the concept of cluster analysis to distinguish the difference of SM from various sources. The inherent characteristics of SM from various source could not be directly reflected by the data of contents. Hence, we introduced the concept of cluster analysis to dig the difference of chemical profile or to find out the rule of variation. Furthermore, it was suitable for us to compare chemical profile and biopotency of SM. Results of the cluster analysis showed that the SM separated into five clusters, as shown in Fig. 3. Cluster 1 included S04 and S05; cluster 2 included S01, S02, and S03; cluster 3 included S07 and S10; cluster 4 contained only S06; and cluster 5 contained S08, S09. The samples in cluster 1 were batches of SM with high liposoluble tanshinone and SalA contents. The samples in cluster 2 were batches of SM with low contents of all. The samples in cluster 3 had high contents of SalB. The S06 sample in cluster 4 was different from the others; it contained the highest content of ProA. The samples in cluster 5 had higher content of all the components except DSS. Although the SM were divided into 5 groups using chemical contents (Fig. 3), the groups were still different with biopotency grades of SM (Fig. 4D). Therefore, chemical profile could not accurately reflect bioactivity of SM.
The investigation results on the influence of DMSO in different concentrations are shown in Fig. 4B. The 0.5% DMSO solution appears to have no influence on platelet aggregation (Fig. 4B). To have a high maximum platelet aggregation rate of vehicle-treated platelet-rich plasma (PRP), the PRP should be adjusted to 4 × 108 platelets per mL (Fig. 4C). The anti-platelet aggregation effects of SM were examined using rat PRP, and the inhibition was calculated into biopotency. As shown in Fig. 4D, the measured biopotency was divided into 3 groups. The low group (S1, S2, S3) had a biopotency of 0.1–0.2 U; the middle group (S5, S6, S7, S10) had a biopotency of 0.2–0.3 U; and the high group (S4, S8, S9) had a biopotency of 0.3–0.4 U. The feasible limit of this bioassay was below 58.31%.
The correlation analysis results of the contents and biopotency are shown in Fig. 5. According to the results, all the contents and sum of contents in SM had different correlations (R = 0.2398–0.8046) with the biopotency of SM.
(1) |
Fig. 5 Correlation analysis results of the content of each constituent and biopotency (A). Correlation analysis results of the sum of contents and biopotency (B). |
Although the biopotency of the nine major constituents were measured and provided in Fig. 6, biopotency of the three constituents (DT1, T1, T2) were zero (Fig. 6B). Thus, we selected biopotency of 6 constituents (DSS, ProA, RosA, SalB, SalA, CT) to calculate ECI (Table 1). The ECI of SM was calculated by the formula in eqn (1). As shown in Table 1, the ECIs of SM were divided into 3 groups: the low group (S1, S2, S3) had ECI values of 1.0–1.6; the middle group (S5, S6, S7, S10) had ECI values of 1.7–2.7; and the high group (S4, S8, S9) had ECI values of 0.3–0.4. The ECI-based grades of SM were almost the same as the bio-based grades.
Samples | Fi ×Ci | ECI | |||||
---|---|---|---|---|---|---|---|
DSS | ProA | RosA | SalB | SalA | CT | ||
S1 | 0.0123 | 0.0003 | 0.0234 | 0.5739 | 0.0065 | 0.9147 | 1.5310 |
S2 | 0.0129 | 0.0003 | 0.0320 | 0.5257 | 0.0056 | 0.5206 | 1.0971 |
S3 | 0.0080 | 0.0004 | 0.0276 | 0.6316 | 0.0036 | 0.3367 | 1.0079 |
S4 | 0.0126 | 0.0004 | 0.0511 | 0.6350 | 0.0369 | 2.7722 | 3.5083 |
S5 | 0.0130 | 0.0003 | 0.0397 | 0.7791 | 0.0259 | 1.8137 | 2.6718 |
S6 | 0.0108 | 0.0006 | 0.0486 | 1.1550 | 0.0036 | 0.6969 | 1.9156 |
S7 | 0.0076 | 0.0003 | 0.0435 | 1.1107 | 0.0156 | 0.7689 | 1.9465 |
S8 | 0.0064 | 0.0004 | 0.0405 | 0.9456 | 0.0181 | 1.9023 | 2.9133 |
S9 | 0.0072 | 0.0004 | 0.0452 | 1.1656 | 0.0174 | 1.8595 | 3.0953 |
S10 | 0.0062 | 0.0003 | 0.0273 | 1.1350 | 0.0132 | 0.5556 | 1.7375 |
Fig. 7 Results of the correlation analysis (A); correlation analysis results of ECI and biopotency (B); residual of contents and residual of the effect-constituents index (ECI) (C). |
Results of the correlation analysis between the ECI and biopotency are shown in Fig. 7B. The ECI of SM has a high correlation (R = 0.9187) with biopotency.
The results of the residual analysis for all the contents and the ECI of SM are shown in Fig. 7C. The residual of the ECI was less than the others. The highest residual value was more than 0.06, and CT had a residual value of 0.028. However, the ECI only had a 0.013 residual value.
Fig. 8 Partial Least Squares Discriminant Analysis (PLS-DA) results of Salvia miltiorrhiza Bge. from different sources (A). Principal Component Analysis (PCA) results of contents (B). |
In addition, the ECI could overcome the limitations of the biological assay by providing good operability, high precision and low cost. The ECI, which is closely related to bioactivity, which in turn may indicate correlation with clinical efficacy, can be implemented by chemical assays after initial bioassay data were collected. The ECI was formed with the sum of products of the biopotency weight (Fi) of the six constituents and their content (Ci). In the present study, two types of constituents, water-soluble and liposoluble chemicals, were evaluated for anti-platelet aggregation effects, and the results were calculated into biopotency (Fi). Notably, CT, a liposoluble constituent, inhibited platelet aggregation significantly and had more than 50–200 folds higher biopotency than the water-soluble constituents (DSS, ProA, RosA, SalB and SalA) (Fig. 6B). Despite the Ci value of CT was approximately 30–80 folds lower than that of the water-soluble constituent SalB (Table S1†), the biopotency (Fi, U mg−1) of CT was over 50 times higher than that of SalB (Fig. 6B). Thus, the active constituents, e.g. CT, with lower concentrations but higher activity levels would contribute considerably to the total bioactivity of SM as the other active constituents with high concentrations but lower activity, e.g. SalB, did. This result indicated the necessity of activity-based calibration for the low concentration but highly active constituents in quality evaluation of SM and the other TCMs or BDs.
The value of Fi is reflected by the differences of the structure–activity relationships of the constituents. There was an obvious relationship between the molecular structure and the biological activity of these constituents. For water-soluble constituents, the ortho-dihydroxyl phenyl moiety (indicated in Fig. 6A) seems important to the anti-platelet aggregation activity. ProA or DSS has one ortho-dihydroxyl phenyl moiety, so biopotency of them were rough equal. Moreover, RosA, SalA and SalB have two, three and four of such moieties, respectively; and the biopotency of these three compounds rises consequently. However, the biopotency of these compounds rise exponentially with the numbers of ortho-dihydroxyl phenyl moiety (Fig. 6D). The results indicated that the higher degree of polymerization has, the stronger biopotency exhibits. For liposoluble constituents, there might be two structural domains, ring A and ring D (Fig. 6A), of close relationship with activity. When ring A is an unsaturated ring, these compounds have no activity, such as DT1 and T1. Meanwhile, as ring D is an unsaturated ring, such as T1 and T2, these components have no activity else. In contrast, CT has both unsaturated ring A and unsaturated ring D, and CT has strong activity. In addition, three-dimensional structure of these compounds was influenced by the factor that whether ring A and ring D were unsaturated or not (Fig. 6E). Therefore, these three-dimensional structures may affect the efficacy of compound by altering combining ability of the compound with target protein. The activity of tanshinone compounds were influenced by three-dimensional structure, and those examples were also found in evaluation of angiogenesis activity. For instance, CT could inhibit angiogenesis in vitro but T2 had no effect.29 Therefore, the results of the structure–activity relationship provides insights into the development of antiplatelet agents with improved effect to prevent platelet hyperactivity-related cardiovascular diseases.
The ECI, establishment based on data from bioassays and chemical assays, was an effect-calibrated quantitative approach and was essentially equivalent to bioassays in predicting clinically relevant activities, and yet using chemical assays only as the routine tests. Thus, this method will have better applicability and can be easily adopted by other research groups, industry, and regulatory agencies for evaluation of herbal products. The ECI is suitable for TCMs or BDs with multiple known active constituents that are also researched for their bioactivities. In this paper, nine constituents were selected, because they were either active constituents on anti-platelet aggregation, or quality markers of SM. After the bioassay, six chemicals, DDS, ProA, RosA, SalB, SalA and CT, were proven to have antiplatelet aggregation activity. Thus, we chose biopotency of the six constituents as calibration factors for calculating ECI. But once new active constituents of anti-platelet aggregation are found in SM, biopotency of them would be used as new calibration factors. Although the ECI is suitable for evaluating TCMs or BDs with multiple active constituents, it also has limitations regarding constituents with multiple activities, which require several bioassays to predict their overall activities and interactions of the whole products. Therefore, it is necessary to characterize the interactive antagonistic or synergistic effects of constituents, as well as applying multiple bioassays to study the whole extracts of the herbs and herbal combinations, to decipher their clinical effects. Based on the activities of the multiple active constituents, and their interactive effects from biological and pharmacological assays, the TCM herbs, and herbal preparations should be further evaluated with updated or new calibration factors for improved quality control of those herbal products.
Each standard was accurately weighed and dissolved in 70% methanol to produce a 0.2 mg mL−1 standard solution of RosA, SalB and a 0.05 mg mL−1 standard solution of DSS, ProA, SalA, DT1, CT, T1, and T2. Each solution was filtered through a 0.22 μm microporous membrane.
Dry powder (3 g) of SM was extracted with ultrasonication in 0.5 liters of 70% methanol. After 1 hour of extraction, the solution was filtered through filter paper (Whatmann filter 1), and the filtrate was filtered through a 0.22 μm microporous membrane before injection.
Clustered heatmaps were another very popular visualization tool in analyst. MetaboAnalyst 3.0 statistics software (MetaboAnalyst 3.0, Xia Lab@McGill University, Canada) was used to conduct the cluster heatmaps. The contents of the constituents were determined by HPLC using the chromatographic conditions. Based on the results, the contents in SM from different sources were selected for cluster analysis by the software MetaboAnalyst 3.0.
The SM test sample was extracted with 70% methanol. The extraction solution was concentrated to one-hundredth of its original volume in a rotary evaporator at 50 °C as stock solution, and this extract was diluted to the appropriate concentration with physiological saline as the test solution. Before the assay, each solution was filtered through a 0.22 μm microporous membrane.
To improve the precision and reliability of the bioassay, the effects of the experimental conditions relating to the experimental inducing agent, standard, platelet concentration, and DMSO concentration were investigated. The standardized bioassay was established based on the optimized experimental conditions.
After a preliminary screen and after considering the maximal platelet aggregation ratio and the response potency of the constituents by induction of adenosine diphosphate, arachidonic acid, collagen, and thrombin, we chose arachidonic acid as the test inducing agent of platelet aggregation for this bioassay. In addition, standard, platelet concentration and the influence of DMSO concentration were also investigated.
A quadra-channel aggregometer (AggRAM, Corporation, Helena, USA) was used to measure platelet aggregation in vitro according to the turbidimetric method. Blood, anticoagulated with 3.8% sodium citrate (1:9 citrate/blood, v/v), was withdrawn from male SD rats (anesthetized by sodium pentobarbital) from the abdominal aorta. Platelet-rich plasma (PRP) was obtained by centrifugation at 800 rpm for 15 min at 25 °C. Platelet-poor plasma (PPP) was prepared from the precipitated fraction of PRP by centrifugation at 3000 rpm for 10 min at 25 °C. The PRP was adjusted to 4 × 108 platelets per mL. Next, 175 μL of PRP was incubated at 37 °C for 1 min in a cuvette with 50 μL of test solution at a final concentration. 0.5% DMSO in distilled water was used as control. After incubation, platelet aggregation was induced by the addition of 25 μL of AA (250 μM). The maximum platelet aggregation rate was recorded within 5 min with continuous stirring at 37 °C. The light transmittance was calibrated with PPP. The percentage (%) of inhibition of platelet aggregation was calculated using eqn (2).
I% = [(X − Y)/X] × 100% | (2) |
The nine standard constituents were accurately weighed and respectively dissolved in 0.5% DMSO diluting water to produce a maximal concentration of solution. Then this solution was diluted to the appropriate concentration with physiological saline as the test solution. Then, the inhibition rate of nine constituents was measured on anti-platelet aggregation, and the method was the same as test of SM which was previously mentioned. Finally, the biopotencies of the nine standard constituents were calculated by the inhibition rate of anti-platelet aggregation.
The bioassay statistical method was established according to the inhibition rate of test samples and the reference. Constituents and SM were test samples, and CT was defined as reference. The biopotency was calculated according to the qualitative response parallel method. In the drug biological inspection, the dose was converted to logarithm, and the inhibition rate was converted into a probability unit, which made the logarithm dose and the reaction function appear in parallel. According to this principle, in the experiment, the dose and the inhibition rate of the reference group and the test group were transformed, as shown in Fig. S1.† The biopotency of CT was defined as 1000 U mg−1 (296350 U mmol−1). The biopotency was calculated according to the qualitative response parallel method.33
We estimated the effect calibration coefficient for each ingredient using the biopotency. The ECI computing formula is listed in eqn (2).
Footnotes |
† Electronic supplementary information (ESI) available: Tables S1–S3. See DOI: 10.1039/c6ra26281c |
‡ These authors contributed equally to this paper. |
This journal is © The Royal Society of Chemistry 2017 |