Abstract

Zuogui Pill is a kidney-yin-tonifying formula in traditional Chinese medicine that is widely used to manage osteoporosis with kidney-yin-deficiency in China. Herein, an efficient and accurate high-performance liquid chromatography-tandem mass spectrometry method was developed to determine the concentrations of 5 bioactive compounds in rat plasma following oral administration of Zuogui Pill. Because drug absorption and distribution differ under physiological and pathological conditions, the established method was used to quantify blood components and dynamic change in osteoporotic rats with different syndrome types. Moreover, integrated pharmacokinetic study was conducted to describe the overall pharmacokinetic characteristics of traditional Chinese medicine. The results showed that the absorption, distribution, and metabolism of Zuogui Pill varied widely under different states. The bioavailability of most active components showed significant advantages in osteoporotic rats with kidney-yin-deficiency, which corresponds to the opinion that Zuogui Pill has the effect of nourishing kidney-yin. It is hoped that this finding could interpret the pharmacodynamic substances and mechanism of Zuogui Pill in the treatment of osteoporosis with kidney-yin-deficiency.

1. Introduction

Osteoporosis is a progressive bone remolding disease characterized by osteopenia and bone structure deterioration. The essence of osteoporosis is an imbalance between bone formation and bone resorption, which can cause skeletal fragility, increased fracture risk, and post-fracture mortality [1]. Recent evidence shows that the incidence of osteoporosis has been increasing year by year. Global osteoporosis-induced hip fracture patients are expected to increase by 550% by 2050 [2, 3]. Hormone replacement therapy (HRT) is currently proven to be effective in preventing bone loss and lowering the risk of fractures in osteoporotic patients [4, 5]. Nonetheless, long-term HRT has been linked to an increased risk of a variety of diseases [6, 7]. Concerns about the adverse effects of HRT have heightened interest in developing alternative strategies for osteoporosis treatment, and traditional Chinese medicine (TCM) has attracted considerable interest due to its superior curative effect and lower side effects. The TCM theory points out that the prevalence of osteoporosis is basically kidney deficiency; clinically, it mainly manifests as kidney-yang-deficiency and kidney-yin-deficiency syndrome [810]. As a result, kidney-tonifying is a common and traditional strategy for treating bone diseases such as osteoporosis with TCM herbal formulas that have kidney-tonifying activity.

Zuogui Pill is a classic kidney-yin-tonifying formula which was developed by a famous TCM doctor Jingyue Zhang during the Ming dynasty (1563–1640 AD). The formula consists of Rehmanniae Radix Praeparata, Dioscoreae Rhizoma, Lycii Fructus, Corni Fructus, Cyathulae Radix, Cuscutae Semen, Cervi Cornus Colla, and Testudinis Carapacis et Plastri Colla. The majority of TCM in this formula have similar and synergistic kidney-yin-tonifying functions and thus have the potential to benefit bone health [11, 12]. Many bioactive ingredients in Zuogui Pill, including aucubin, loganin, chlorogenic acid, and morroniside can promote osteoblast differentiation, reduce calcitonin levels, and inhibit osteoblast apoptosis and osteoclast proliferation [1317]. Meanwhile, protocatechuic aldehyde has been shown to be effective in the prevention and treatment of osteoporosis [18, 19]. Therefore, we selected the 5 bioactive components of Zuogui Pill (Figure 1) for further study.

Mounting evidence shows that differences in body conditions can influence drug absorption and distribution [20, 21]. From the perspectives of serum pharmacochemistry research ideas, only the components that can be absorbed are important for clinical effects [22]. Thus, studying pharmacokinetic properties in rats under pathological and normal conditions may reveal the mechanism of Zuogui Pill in the treatment of osteoporosis. Up to now, the methods for qualitative and quantitative detection of chlorogenic acid, loganin, and morroniside in Zuogui Pill have been established [23, 24]. There is no pharmacokinetic study on the active ingredients of Zuogui Pill in rat plasma. In this study, we developed an efficient and rapid HPLC-MS/MS method for detecting 5 components in plasma. Following that, the developed method was successfully used to compare the pharmacokinetic characteristics of Zuogui Pill after oral administration in normal rats, osteoporotic rats, and osteoporotic rats with kidney-yang-deficiency and kidney-yin-deficiency syndrome. Due to the complex composition of TCM compound, any single component poses limitations in describing the overall pharmacokinetic characteristics of TCM. Consequently, a study of integrated pharmacokinetics in each syndrome was conducted for systematic evaluation [25].

2. Materials and Methods

2.1. Chemicals and Reagents

Rehmanniae Radix Praeparata, Dioscoreae Rhizoma, Lycii Fructus, Corni Fructus, Cyathulae Radix, Cuscutae Semen, Cervi Cornus Colla, and Testudinis Carapacis were purchased from Huadong Pharmaceutical (Hangzhou, China) and confirmed by Professor Rusong Zhang (School of Pharmaceutical Science, Zhejiang Chinese Medical University).

Aucubin, loganin, protocatechuic aldehyde, chlorogenic acid, and tinidazole (internal standard, IS) (purity > 98%) were obtained from the National Institutes for Food and Drug Control (Beijing, China). Morroniside (purity > 98%) was procured from Beijing Yanxinlv Biological Technology (Beijing, China). Methanol, acetonitrile, ammonium formate, and formic acid (HPLC-grade) were obtained from Merck (Darmstadt, Germany). Hydrocortisone injections were provided by Jinyao Amino Acid Co. Ltd. (Tianjin, China). Thyroxine tablets were supplied by Great Wall Pharmaceutical Co. Ltd. (Shanghai, China). Ultrapure water was purchased from Wahaha Co. Ltd. (Hangzhou, China).

2.2. Preparation of Zuogui Pill Extract

The Zuogui Pill extract was prepared using a previously described method [26]. The raw materials, which included Rehmanniae Radix Praeparata (24 g), Dioscoreae Rhizoma (12 g), Lycii Fructus (12 g), Corni Fructus (12 g), Cyathulae Radix (9 g), Cuscutae Semen (12 g), Cervi Cornus Colla (12 g), and Testudinis Carapacis (12 g), were sliced and soaked in a 10-fold volume of 50% ethanol for 60 min. After that, a 30 min ultrasonic extraction was performed twice at 45°C. The extract solution was filtered and condensed to a 4 g/ml concentration.

2.3. Animals and Establishment of Osteoporosis Models

Female Sprague-Dawley rats (180–220 g) were purchased from Shanghai Laboratory Animal Centre (Shanghai, China). The rats were maintained in a standard environment (22 ± 1°C, 60 ± 10% relative humidity, and a 12 h light/dark cycle) for one week to acclimate and fed water and food ad libitum. Twenty-four rats were divided randomly into three groups and treated through ovary enucleation to establish the model of osteoporosis [27]. Ten weeks after operation, osteoporosis models of kidney-yang-deficiency and kidney-yin-deficiency were established by intramuscular injection of hydrocortisone (25 mg/kg) and intragastric administration of thyroxine suspension (160 mg/kg) for 10 days, respectively. The femur of two rats from each group was then stained with hematoxylin and eosin (H&E) to assess histopathological change. Also, the validation of osteoporosis kidney-yang-deficiency and kidney-yin-deficiency models was judged by detecting the level of cyclic adenosine monophosphate (cAMP) and cyclic guanosine monophosphate (cGMP) in rat plasma [28]. Laboratory Animal Management and Ethics Committee of Zhejiang Chinese Medical University approved all animal treatments (ZSLL-2014-48).

2.4. Preparation of Standards and Quality Controls

The standards of 5 active compounds and IS were weighed and solubilized in methanol to 1.0 mg/mL. To make working solutions, stock solutions were diluted with methanol, and the IS working solutions were diluted in a methanol : acetonitrile ratio of 1 : 1 at a final concentration of 50 ng/mL.

To obtain calibration standards, 10 μL working solutions were added in 90 μL blank plasma. The final concentrations were 1–1000 ng/mL for chlorogenic acid, 20–20000 ng/mL for loganin, 4–4000 ng/mL for morroniside, 0.5–500 ng/mL for aucubin, and 0.05–50 ng/mL for protocatechuic aldehyde.

Quality control (QC) samples were prepared at the following concentrations: 2, 40, and 800 ng/mL for chlorogenic acid, 40, 800, and 16000 ng/mL for loganin, 8, 160, and 3200 ng/mL for morroniside, 1, 20, and 400 ng/mL for aucubin, and 0.1, 2, and 40 ng/mL for protocatechuic aldehyde.

2.5. Sample Preparation

One hundred microliters of plasma was mixed with 100 μL of the IS solution, followed by 700 μL of methanol:acetonitrile (1 : 1). After 2 min of vortex mixing, the mixture was spun for 15 min at 8000 rpm and 4°C. The supernatants were collected and dried in a high-purity N2 environment. The residues were redissolved in 100 μL methanol and centrifuged for 15 min at 13000 rpm and 4°C. Finally, an aliquot of 10 μL supernatants was used to conduct the analysis.

2.6. Instrument and Analytical Conditions

An Agilent 1290 HPLC system (Agilent, California, USA) was used in conjunction with an AB Sciex Q-TRAP 5500 MS/MS system (AB Sciex, Foster City, USA) equipped with an ESI source to measure the 5 components.

The LC separation was performed on a Waters Atlantis C18 column (3.0 mm × 150 mm, 3 μm) with a mobile phase of 5 mmol/L ammonium formate in 0.05% formic acid aqueous solution (A) and acetonitrile (B) flowing at a rate of 0.3 mL/min. The gradient elution program was performed as follows: 0-1 min, 90% A; 1–5 min, 90–20% A; 5–10 min, 20% A; 10–10.1 min, 20–90% A; and 10.1–12 min, 90% A. The column temperature was 20°C.

The mass spectrometer was set to detect in a negative mode under the following conditions: ion spray voltage, 4500 V; curtain gas, 40 psi; source temperature, 650°C; and collision gas, medium. The analytes were quantified using multiple reaction monitoring, and the optimized parameters are listed in Table 1.

2.7. Method Validation

The methods were validated in accordance with US Food and Drug Administration guidelines for bioanalytical method validation [29].

2.7.1. Specificity

Chromatograms of blank biological samples from six rats, plasma samples containing standards, and plasma samples after administration of Zuogui Pill extract were compared to assess selectivity.

2.7.2. Linearity and Sensitivity

Calibration curves were constructed by drawing peak area ratios of 5 compounds to IS against nominal concentrations using least squares regression (weighted factor: 1/x2). The lower limit of quantification (LLOQ) was defined as the lowest level on calibration curves with a signal-to-noise ratio (S/N) = 10 used to demonstrate method sensitivity.

2.7.3. Accuracy and Precision

Six batches of QC samples at three levels and LLOQ were used to determine accuracy and precision. To determine intra-day accuracy and precision, samples were analysed on a single assay day, while inter-day accuracy and precision were determined over three consecutive days.

2.7.4. Extraction Recovery and Matrix Effect

To determine extraction recovery, peak areas of analytes in plasma samples were compared to post-extraction plasma containing pure analyte standards at corresponding concentrations. The matrix effect was calculated by comparing the peak areas of analytes added to blank plasma after extraction to those of standards liquefied in methanol at comparable concentrations.

2.7.5. Stability

QC samples at three levels using five replicates were used to evaluate the stability under the following conditions: 4 h at room temperature (25°C), 24 h in autosampler (4°C), and three freeze-thaw cycles (−20°C to 25°C).

2.8. Pharmacokinetic Study and Statistical Analysis

Four groups of rats were orally administered with Zuogui Pill extract at a dose of 1.5 mL/100 g. Then, 0.3 mL blood samples were collected from orbital veins at 0, 0.08, 0.17, 0.25, 0.5, 1, 2, 4, 6, 8, 12, and 24 h after administration, placed in an anticoagulant tube, and centrifuged for 15 min at 5000 rpm and 4°C. The plasma was prepared and stored at −20°C until analysis.

Phoenix WinNonlin 6.4 software was used to estimate pharmacokinetic parameters using noncompartmental analysis. Furthermore, SPSS 26 software was used to compute comparisons of the pharmacokinetic data among the four groups according to a one-way ANOVA. The results are presented as mean ± SD.

2.9. Integrated Pharmacokinetic Study

Integrated pharmacokinetic study was conducted based on the area under the concentration-time curve from zero to infinity (AUC0-∞) of each component [30, 31]. The ratio of AUC0-∞ of 5 components to total AUC0-∞ was defined as the weight coefficient () in the integrated concentration using equations (1) and (2). Formula (3) was used for obtaining the integrated concentration of the 5 components, where represent the concentration of the five components, and the integrated pharmacokinetic parameters were further calculated.

3. Results and Discussion

3.1. Evaluation of Osteoporosis Models

H&E staining revealed that the subchondral bone and trabecula were intact and arranged regularly in the normal rats. Osteocytes in trabecular bone were clearly visible, and the nuclei were centrally located and large. The bone marrow was rich in hematopoietic cells, with relatively few adipocytes in the medullary cavity and normal morphology, while the trabecular bone area and the number of hematopoietic cells were significantly reduced in the osteoporosis group. In addition, the volume of adipocytes increased and even fused into vesicles (Figure 2), both of which confirmed the successful establishment of the osteoporosis model.

The contents of cAMP and cGMP in rat plasma are covered in Table 2. After modelling, the plasma content of cAMP in the osteoporotic rats with kidney-yang-deficiency significantly decreased (), while that of cGMP increased () compared with the osteoporosis group. Also, the levels of plasma cAMP and cGMP showed the opposite trend in kidney-yin-deficiency rats. The increase or decrease of the cAMP/cGMP ratio indicated that the two types of osteoporosis models were successfully established.

3.2. Operational Condition Optimization

Considering the low drug concentration in plasma due to its first-pass effect, QTOF-MS was used for analysis in the pre-experiment. However, the sensitivity of QTOF-MS cannot meet the quantitative requirements; therefore, Q-TRAP MS was employed for detection in the end.

Different types of columns were tested, including Agilent ZORJBAX Extend C18 column (4.6 mm × 150 mm, 5 μm), Agilent ZORJBAX SB C18 column (4.6 mm × 150 mm, 5 μm) and Waters Atlantis C18 column (3.0 mm × 150 mm, 3 μm), and it turned out that Waters Atlantis C18 column (3.0 mm × 150 mm, 3 μm) could present a better peak shape and isolation. Acetonitrile, methanol, and several volatile acid-base additives in the aqueous phase such as formic acid, ammonia, ammonium formate, and ammonium acetate were evaluated to optimize the chromatography conditions. Finally, acetonitrile and 5 mmol/L ammonium formate in 0.05% formic acid aqueous solution were chosen as the mobile-phase system.

Because the response value of analytes to be measured was higher in the negative ion ionization state than in the positive mode, the former was used for the assay. Figure 3 depicts the ion mass spectrums.

3.3. Concentrations of the 5 Compounds in Zuogui Pill

The contents of 5 analytes were evaluated according to the method described in Section 2.6. The contents of chlorogenic acid, loganin, morroniside, aucubin, and protocatechuic aldehyde in Zuogui Pill were 33.2, 479, 411, 16.2, and 5.61 μg/mL, respectively.

3.4. Method Validation
3.4.1. Specificity

Figure 4 shows chromatograms of different plasma, including blank plasma, blank plasma with standards, and plasma sample from a rat 30 min after oral administration. The results demonstrated that endogenous plasma components did not influence the test substance or internal standard. Also, the IS response in the blank sample was no more than 5% of the average IS responses of the calibrators and QCs.

3.4.2. Linearity and Sensitivity

Analyte calibration curves were linear over their respective ranges, with correlation coefficients (R) greater than 0.9990 for all sequences. The data in Table 3 demonstrated that the linearity and LLOQ were suitable for assay with acceptable precision and accuracy.

3.4.3. Accuracy and Precision

Table 4 shows the results of the inter- and intra-day accuracy and precision of the 5 components. The precision was less than 15%, while the accuracy was 2.3% to 18.6% for intra-day assay and 2.5% to 18.0% for inter-day assay. All the data showed that the method was reliable and accurate.

3.4.4. Extraction Recovery and Matrix Effect

The extraction recovery of the analytes at different selected levels ranged from 65.0% to 91.9%, and the average matrix effect at the same concentrations ranged from 76.3% to 94.4% (Table 5). The absence of matrix components in plasma results in a striking variation in the mass spectrum response to analytes and IS, which is consistent with the analytical requirements of the biological sample.

3.4.5. Stability

The stability analysis of the 5 compounds revealed no obvious degradation of components detected in rat plasma under the above storage conditions with RSD below 15% (Table 6).

3.5. Pharmacokinetic Analysis

The validated method was used for pharmacokinetic assays of the 5 active compounds in plasma following oral administration of Zuogui Pill. Figure 5shows the mean plasma concentrations versus time profiles of 5 compounds and integration under various syndrome types.Table 7 displays the pharmacokinetic variables of the 5 compounds as well as the integrated pharmacokinetic variables.

Although the dosage of Zuogui Pill in each group was the same, the pharmacokinetic parameters of the four groups were dramatically different. The elimination half-life (t1/2) of chlorogenic acid, loganin, morroniside, and aucubin in the osteoporosis group was shorter than that in the normal group, indicating that they were eliminated faster in the osteoporosis state. Clearance (CL) for the 5 compounds was increased in osteoporotic rats. All the components were absorbed and reached their peak plasma concentration in 1 h. Aside from aucubin, the maximum plasma concentration (Cmax) of the remaining four components was higher in normal rats than in osteoporotic rats. Furthermore, in the osteoporosis group, the AUC0–t and AUC0–∞ of all ingredients were reduced, particularly in protocatechuic aldehyde (). By comparing the pharmacokinetic parameters in normal and osteoporosis groups, we found that the bioavailability of the 5 components in normal rats was higher than that in osteoporotic rats, while the removal rates of the 5 components in normal rats were lower.

The clearance of all analytes was statistically increased () in kidney-yang-deficiency rats, while the AUC0–t and AUC0–∞ of chlorogenic acid, loganin, aucubin, and protocatechuic aldehyde were significantly decreased (), suggesting that the above components were poorly absorbed and rapidly metabolized in this state. Except for protocatechuic aldehyde, Cmax of other compounds increased remarkably () in the kidney-yin-deficiency group when compared with the osteoporosis group. Increases in the value of mean residence time (MRT), AUC0–t, and AUC0–∞ of loganin, morroniside, and aucubin were observed () in the kidney-yin-deficiency group. The results demonstrated that the three components had better absorption and a long duration, which may lead to a significant accumulation of active components in vivo. Thus, loganin, morroniside, and aucubin are likely to be the medicinal compositions in the treatment of osteoporosis with kidney-yin-deficiency in Zuogui Pill extract as reported in the literature [32].

The above changes in pharmacokinetic behaviours could be due to the fact that different types of osteoporosis potentially influence the absorption, distribution, and metabolism of drugs. One of the factors influencing the rate and quantity of drug absorption in osteoporosis is the change in gut microbiota and metabolite [33, 34]. The phenomenon that all the components showed double or multiple peaks in the mean plasma concentration-time profiles could be attributed to Zuogui Pill reabsorption in plasma and enterohepatic circulation following oral administration [35, 36].

3.6. Integrated Pharmacokinetic Analysis

TCM compounds are formulated in accordance with the composition theory of a TCM prescription, which contains complex chemical components that potentially interact with one another to influence their pharmacokinetic behaviours [37]. The results of our assay revealed significant differences in the pharmacokinetic parameters of 5 compounds, demonstrating that any single ingredient poses limitations in describing pharmacokinetic behaviour in rats. The integrated pharmacokinetic parameters, on the other hand, were more practical and reasonable. According to Figure 5 and Table 6, it was obvious that the bioavailability and effective drug duration of integrated ingredient showed significant advantages in the state of kidney-yin-deficiency by comparison with other groups, which was consistent with the TCM theory that Zuogui Pill has the effect of nourishing kidney-yin [38].

4. Conclusions

The presently developed HPLC-MS/MS method is efficient and effective in detecting the 5 analytes in plasma. Therefore, it was used to compare the pharmacokinetics of 5 components in different syndrome types of osteoporosis. The evaluation of the pharmacokinetic behaviours of the research objects may aid in better understanding the mechanism of the pharmacological action and provide references on the clinical application of Zuogui Pill.

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 that there are no conflicts of interest regarding the publication of this paper.

Authors’ Contributions

Jiawei Qiu and Yaoyao Zhu contributed equally to this study.

Acknowledgments

This study was supported by the National Natural Science Foundation of China (81973868 and 81273772) and Natural Science Foundation of Zhejiang Province (LY17H270011).