Pharmacodynamics of Five Anthraquinones (Aloe-emodin, Emodin, Rhein, Chysophanol, and Physcion) and Reciprocal Pharmacokinetic Interaction in Rats with Cerebral Ischemia

(1) Background: Rhubarb anthraquinones—a class of components with neuroprotective function—can be used to alleviate cerebral ischemia reperfusion injury. (2) Methods: The three pharmacodynamic indicators are neurological function score, brain water content, and cerebral infarction area; UPLC-MS/MS was used in pharmacokinetic studies to detect plasma concentrations at different time points, and DAS software was used to calculate pharmacokinetic parameters in a noncompartmental model. (3) Results: The results showed that the pharmacodynamics and pharmacokinetics of one of the five anthraquinone aglycones could be modified by the other four anthraquinones, and the degree of interaction between different anthraquinones was different. The chrysophanol group showed the greatest reduction in pharmacodynamic indicators comparing with other four groups where the rats were administered one of the five anthraquinones, and there was no significant difference between the nimodipine group. While the Aloe-emodin + Physcion group showed the most obvious anti-ischemic effect among the groups where the subjects were administered two of the five anthraquinones simultaneously. Emodin, rhein, chrysophanol, and physcion all increase plasma exposure levels of aloe-emodin, while aloe-emodin lower their plasma exposure levels. (4) Conclusions: This experiment provides a certain preclinical basis for the study of anthraquinone aglycones against cerebral ischemia and a theoretical basis for the study of the mechanism of interaction between anthraquinones.


Introduction
Cerebral ischemic stroke is a major cause of severe acquired disability worldwide and the morbidity approximately account for 80% of global strokes [1,2]. However, the use of a thrombolysis agent is the only FDA-approved method for the treatment of acute ischemic stroke so far and there is a limited therapeutic window within 3 h, after which there is potential risk of hemorrhagic transformation [3][4][5][6].

Pharmacodynamic Results
The neurological function scores, brain water content, and infarct size of the sham operation group, the model group, the positive group, and the 16 experimental groups are shown in Table 1 and Figure 1. Comparing the model group with the sham operation group, the neurological function score and brain tissue water content were significantly increased and the ischemic brain tissue staining showed a large area of white cerebral infarction. The results indicated that the cerebral ischemia model was established successfully and the pharmacodynamic characteristics could be reflected with the three selected pharmacodynamic indicators. Compared with the model groups, the three pharmacodynamic indicators were significantly decreased in the positive group and the Aloe-emodin, Rhein, Chrysophanol, Aloe-emodin + Rhein, and Aloe-emodin + Physcion groups. Besides, the Aloe-emodin + Physcion group, which showed no significant pharmacodynamic differences with the positive group, possessed the strongest efficacy of the 15 experimental groups. Table 1. Pharmacodynamic index: the neurological function scores, the water content, and area of cerebral infarction. The results were expressed as means ± standard deviation. One-way ANOVA and bilateral inspection were used to analyze the data differences between the groups by some standards that ∆ , * p < 0.05; ∆∆ , ** p < 0.01. ∆ showed the comparison between the sham operation group, * showed the comparison between the model group.

Group
Neurological Function Score The Water Content (%)

Area of Cerebral Infarction (%)
Sham operation group -78. 40  and Figure 1. Comparing the model group with the sham operation group, the neurological function score and brain tissue water content were significantly increased and the ischemic brain tissue staining showed a large area of white cerebral infarction. The results indicated that the cerebral ischemia model was established successfully and the pharmacodynamic characteristics could be reflected with the three selected pharmacodynamic indicators. Compared with the model groups, the three pharmacodynamic indicators were significantly decreased in the positive group and the Aloeemodin, Rhein, Chrysophanol, Aloe-emodin + Rhein, and Aloe-emodin + Physcion groups. Besides, the Aloe-emodin + Physcion group, which showed no significant pharmacodynamic differences with the positive group, possessed the strongest efficacy of the 15 experimental groups. Table 1. Pharmacodynamic index: the neurological function scores, the water content, and area of cerebral infarction. The results were expressed as means ± standard deviation. One-way ANOVA and bilateral inspection were used to analyze the data differences between the groups by some standards that Δ , * p < 0.05; ΔΔ , ** p < 0.01. Δ showed the comparison between the sham operation group, * showed the comparison between the model group.

Pharmacokinetics Interaction Studies
According to the guidelines for the bioanalytical validation method issued by the FDA (Center for Drug Evaluation and Research [25]), we conducted methodology validation of selectivity, linearity, LOQ (limit of quantitation), accuracy, precision, extraction recovery, matrix effect, and stability. The methodological verification results are shown in Tables S1-S4 and Figure S2. The chromatographic method has proven to be linear (r > 0.997), precise (relative standard deviation RSD < 8%), and accurate (relative error RE ± 6%) in the concentration range of the experiment. The minimum detection limits for aloe-emodin, emodin, rhein, chrysophanol, and physcion were 0.026, 0.098, 0.098, 0.538, and 0.317 ng, respectively. The experimental conditions we selected can detect the blood concentration of five rhubarb aglycones rapidly, stably, exclusively, and high sensitively.
For each animal group the calculated mean plasma concentrations of the anthraquinone aglycones, plotted as a function of time, are shown in Figure 2. At all sampling times, the plasma level of aloe-emodin was increased by the co-administration of emodin, rhein, chrysophanol, or physcion; while plasma levels of emodin, rhein, and chrysophanol decreased when administered with aloe-emodin, emodin, rhein, chrysophanol, and physcion. Moreover, different pharmacokinetic parameters of the five anthraquinone aglycones were calculated using noncompartmental analysis, as presented in Table 2. Significant differences in AUC 0-t , AUC 0-∞ , AUMC 0-∞ , and CL/F values were observed when comparing the groups in which subjects were administered one of the anthraquinone aglycones (aloe-emodin, emodin, rhein, chrysophanol, or physcion) with that groups where subjects were administered two of the anthraquinone aglycones. According to the guidelines for the bioanalytical validation method issued by the FDA (Center for Drug Evaluation and Research [25]), we conducted methodology validation of selectivity, linearity, LOQ (limit of quantitation), accuracy, precision, extraction recovery, matrix effect, and stability. The methodological verification results are shown in Tables S1-S4 and Figure S2. The chromatographic method has proven to be linear (r > 0.997), precise (relative standard deviation RSD < 8%), and accurate (relative error RE ± 6%) in the concentration range of the experiment. The minimum detection limits for aloe-emodin, emodin, rhein, chrysophanol, and physcion were 0.026, 0.098, 0.098, 0.538, and 0.317 ng, respectively. The experimental conditions we selected can detect the blood concentration of five rhubarb aglycones rapidly, stably, exclusively, and high sensitively.
For each animal group the calculated mean plasma concentrations of the anthraquinone aglycones, plotted as a function of time, are shown in Figure 2. At all sampling times, the plasma level of aloe-emodin was increased by the co-administration of emodin, rhein, chrysophanol, or physcion; while plasma levels of emodin, rhein, and chrysophanol decreased when administered with aloe-emodin, emodin, rhein, chrysophanol, and physcion. Moreover, different pharmacokinetic parameters of the five anthraquinone aglycones were calculated using noncompartmental analysis, as presented in Table 2. Significant differences in AUC0-t, AUC0-∞, AUMC0-∞, and CL/F values were observed when comparing the groups in which subjects were administered one of the anthraquinone aglycones (aloe-emodin, emodin, rhein, chrysophanol, or physcion) with that groups where subjects were administered two of the anthraquinone aglycones.   Table 2. Pharmacokinetic parameters of 5 substances in different groups (mean ± SD; n = 10). The pharmacokinetic parameters were calculated according to the blood concentration time of each subject, and then the mean values and SD values of 10 individual pharmacokinetic parameters in each group were calculated. One-way ANOVA and bilateral test were used to compare the pharmacokinetic differences between the compatibility groups and their corresponding groups which the subjects were administered one of the aglycones (* p < 0.05; ** p < 0.01).

Group
Parameters

Discussion
In this experiment, the matrix effect and extraction recovery rate of five rhubarb aglycones and an internal standard were investigated under the conditions that methanol, acetonitrile, and methanol: acetonitrile (1:1) were used as precipitation solvent, respectively. The results showed that when pure methanol was used as a protein removal agent the impurity interference was minimal and the matrix effect were in the range of 75 to 115%; recovery rates were higher than 90%. Next, we compared four mobile phase systems-methanol-water (A), acetonitrile-water (B), methanol-0.1% formic acid water (C), and methanol-3 mmol/mL ammonium acetate (D). It was found that the organic phase acetonitrile did not possess peak shape following to inaccurate integral value; meanwhile methanol could avoid this phenomenon. Formic acid water can improve the detection of the quasi-ion peak [M + H] + of rhubarb aglycone. However, according to the fact that the five aglycones contain a plurality of phenolic hydroxyl groups [17,26], the five anthraquinones are more likely to exhibit a correspondingly strong signal quasi-ion peak [M − H] − in an alkaline environment. Therefore, methanol-3 mmol/mL ammonium acetate (D) was finally selected as the mobile phase system. Referring to published literature [27], the C max , t 1/2 , and AUC 0-t in thrombotic focal cerebral ischemia (TFCI)-induced rats were almost twice that in normal rats, and the CL value was significantly lower than in normal rats (p < 0.05), showing that the pharmacokinetic characteristics of the five anthraquinones were modified under pathological conditions. Thus, in order to match the treatment background, cerebral ischemia reperfusion model rats were selected as the research objects in this experiment. In order to provide some theoretical basis for the drug administration methods of the five drugs (single or combined), our research studied the pharmacokinetic characteristics and pharmacodynamics of the five anthraquinone glycosides in the cerebral ischemia reperfusion model rats and the pharmacokinetic interactions and pharmacodynamics of the five anthraquinone glycosides in the pairwise administration.
Considering the drug concentration-time curve, aloe-emodin exhibits different degrees of absorption double-peaks. This may be because, when administered orally, aloe-emodin was easily combined with glucose in the form of glycoside [16], and when the concentration of aloe-emodin in the blood reduced, a small amount of glycoside was hydrolyzed, so that aloe-emodin has a smaller peak. After all, aloe-emodin still exists in the form of a metabolite [28]. This phenomenon was consistent with previous studies; aloe-emodin met the two-compartment model [16]. The four aglycones (emodin, rhein, chrysophanol, and physcion) all prolonged aloe-emodin peak times, T max ; increased peak concentrations, C max ; reduced elimination rates; and extended t 1/2 . The increasing AUC 0-t and AUC 0-∞ showed increasing aloe-emodin plasma exposure level. Therefore, the absorption was increased, the elimination was slowed, and action time was prolonged. In addition, the decrease in apparent volume of distribution may cause the drug to concentrate on the pathological site, reducing the side effects of aloe-emodin. Comparing the effects of the four anthraquinones on the pharmacokinetics of aloe-emodin, physcion possessed the greatest influence on it, while rhein showed the least effect. Aloe-emodin, rhein, chrysophanol, and physcion all reduced the AUC 0-t value of Emodin, and the AUC 0-t value decreased to 67.1%, 72.1%, 59.7%, and 35.9%, respectively. However, they have different ways of changing the pharmacokinetics of emodin. Aloe-emodin increases the exposure level of emodin by shortening and increasing the apparent volume of distribution, V, while physcion shortens the t 1/2 and improves the clearance rate. The significantly shortened T max indicated that the absorption of emodin was faster and the effect was rapider with the co-administered with aloe-emodin, rhein, chrysophanol, or physcion. Aloe-emodin, emodin, chrysophanol, and physcion all significantly reduced the AUC 0-t , AUC 0-∞ , and C max of Rhein and significantly increased T max , V, and CL, which suggested that the drug-time curves become gentle and the absorption is slowed down. This result may create the reductive possibility of drug accumulation in a certain part to cause the side effects. After all, the drug s therapeutic effect or toxic side effects often depend on whether the drug dose or concentration exceeds a certain value. In general, aloe-emodin can decrease the AUC of the other four anthraquinones and increase the clearance rate and apparent distribution volume, indicating that aloe-emodin has an inhibitory effect on their utilization; meanwhile, emodin, rhein, and chrysophanol increase and reduce the absorption of aloe-emodin. A large number of experiments [29][30][31][32] have shown that both emodin and chrysophanol can inhibit the expression of TNF-α, IL-1β, and vascular adhesion molecule (VCAM-1), enhance the activity of antioxidant enzymes, inhibit the excessive production of NO, scavenge free radicals, and inhibit apoptosis after cerebral ischemic injury. Thus, many identical neuroprotective functions have demonstrated that these two components may have the same targets in the treatment of cerebral ischemia and exhibit competitive inhibition, which may also be the cause of chrysophanol can cause emodin performed reduced absorption, accelerated elimination, and extensive distribution. The effects of emodin, rhein, and chrysophanol on their mutual pharmacokinetic characteristics are similar.
Three pharmacodynamic indicators were considered as a whole. The anticerebral ischemia effect of chrysophanol was the most obvious when subjects were administered one of the five anthraquinones, while the effect of physcion was not obvious. According to the literature [33,34], chrysophanol can attenuate nitrosative/oxidative stress injury and inhibit the endoplasmic reticulum stress response induced by ischemia/reperfusion damage in the brain. In our study, chrysophanol significantly reduced brain water content and cerebral infarct size with no significant difference from positive pharmacodynamics. There was no significant difference in the neurological function scores between the physcion and the model group. However, the pharmacodynamic index, water content, and cerebral infarct size values were higher than other anthraquinone groups. The three pharmacodynamic indicators of aloe-emodin + rhein and aloe-emodin + physcion were significantly different from the model group and were statistically significant. Comprehensive evaluation of the values of the three indicators, the neurological function score, brain water content and cerebral infarct size values in the aloe-emodin + physcion group were smaller than other compatibility groups, indicating that the group showed the best effect among the compatibility groups. The three indicators of the aloe-emodin + physcion group were lower than those in the aloe-emodin group and physcion group, but, comparing with the model group, the deviation was not equal to the sum of reduction values in the corresponding groups in which the subjects were administered only one of the two anthraquinones. Aloe-emodin and physcion increased the reciprocal anti-ischemic effect, and their interaction was not a simple additive effect but a synergistic effect. Comparing the group with emodin and the compatibility groups with the combination of emodin and the other four anthraquinones, physcion reduced the two pharmacodynamic indicators values of cerebral water content and cerebral infarction area significantly. Therefore, when emodin was combined with other rhubarb aglycones to treat cerebral ischemic injury, physcion may be the best choose. The three indexes in the chrysophanol group were lower than the chrysophanol compatibility groups, and the cerebral infarction area was significantly different from the chrysophanol compatibility groups, indicating that the anticerebral ischemia effect of chrysophanol was better than the chrysophanol compatibility group.
Combined with pharmacodynamic results and pharmacokinetic results, it was found that the pharmacodynamic intensity of the drugs often related to the plasma concentration of the drug in plasma. In the pharmacokinetic study: the AUC 0-t of chrysophanol was significantly decreased under the influence of the other four anthraquinone aglycones. Aloe-emodin, emodin, rhein, and physcion reduced the AUC 0-t of chrysophanol by 15.5%, 48.4%, 4.1%, and 33.3%, respectively, and had similar trends with brain water content in pharmacodynamic studies. Brain water content: Emodin + Chrysophanol > Physcion + Chrysophanol > Aloe-emodin + Chrysophanol > Rhein + Chrysophanol > Chrysophanol.

Chemicals and Materials
The reference standards of aloe-emodin, emodin, rhein, chrysophanol, physcion, and 1,8-dihydroxyanthraquinone (internal standard (IS)) were purchased from the National Institutes for Food and Drug Control (Beijing, China). The chemical purity of these reference substances was higher than 98% as determined by HPLC analysis. The raw materials had a purity of ≥ 90%; aloe-emodin, emodin, rhein, chrysophanol, and physcion were obtained from Xi an Xin Rui Bio Technology Co., Ltd.  (Fairfield, OH, USA), and ammonium acetate, purchased from Fisher Company (Pittsburgh, PA, USA), were of mass spectrum grade. The ultrapure water was produced by the Millipore Milli-Q system (Bedford, MA, USA). When aloe-emodin, emodin, rhein, chrysophanol, and physcion were used, the corresponding concentration suspension was prepared with 0.5% sodium carboxymethyl cellulose.

Animals and Cerebral Ischemia Reperfusion Model in Rats
Male Sprague-Dawley (SD) Rats weighing 260-300g were purchased from the Hebei Province Experimental Animal Center (Hebei, China, SCXK2012-1-003). The study protocol is as described by the Experimental Animal Care and Ethics Committee of the First Affiliated Hospital, Henan University of Chinese Medicine (Henan, China, SYXK (Yu) 2015-0005). Before operation, all rats were exposed to 12 h light/12 h dark, temperature 22 ± 2 • C, relative humidity 50 ± 2% conditions for 7 days to adapt to the environment. All breeding and research on experimental animals strictly abide by the regulations on the administration of experimental animals in Henan province. In this study, we used a method previously known as middle cerebral artery occlusion [35][36][37][38] (MCAO) to model cerebral ischemia reperfusion in rats.
Briefly, after anesthesia using 10% chloral hydrate, we separated the left common carotid artery (CCA), external carotid artery (ECA) and internal carotid artery (ICA). A six-0 nylon suture monofilament with a rounded tip (Beijing, China, Beijing salon Biotechnology Co., Ltd.) was inserted through a small incision on the left common carotid artery and moved forward into the internal carotid artery until small resistance was felt, which signifies blocking the entrance of middle cerebral artery. Then ligate the proximal end of the internal carotid artery, suture the wound, and place the rats in a rat cage with a clean pad. At this point, the brain deficit model has been established. After 2 h, the filament was slowly withdrawn until there is resistance to allow reperfusion. The sham group exposed only the left side of the vessel without the insertion.

Grouping and Administration
After 7 days of acclimatization, 180 rats were randomly assigned to one of the following 18 groups (n = 10 each); Sham operation group, MCAO model group, nimodipine positive control group, and 15 experimental groups MCAO + (Aloe-emodin; Emodin; Rhein; Chrysophanol; Physcion; Aloe-emodin + Emodin; Aloe-emodin+ Rhein; Aloe-emodin+ Chrysophanol; Aloe-emodin+ Physcion; Emodin + Rhein; Emodin + Chrysophanol; Emodin + Physcion; Rhein + Chrysophanol; Rhein + Physcion; Chrysophanol + Physcion). Drug concentrations were as follows; Aloe-emodin 13.05 mg/kg, Emodin 34.65 mg/kg, Rhein 17.25 mg/kg, Chrysophanol 39.45 mg/kg, Physcion 26.4 mg/kg, and the positive drug nimodipine 10.8 mg/kg. All of the above drugs were given orally with the corresponding concentration of suspension prepared with 0.5%CMC-Na before use. The dosage was calculated at one percent of the rat s body weight. The sham operation group and the model group were given 0.5%CMC-Na of the corresponding volume. Rats in the experimental groups and the positive group were intragastrically administered drugs once a day for seven consecutive days. The rats were fasted but water-free at 12 h before the last dose, and then a cerebral apoplexy reperfusion model was established 2 h after the last administration.

Measurement and Evaluation
(1) The neurological function of rats was evaluated by Longa method [39][40][41]. Scoring criteria: 0 points: no symptoms of neurological deficit, normal activity; 1 point: cannot extend the opposite front claw completely; 2 points: turn to hemiplegia side when crawling; 3 points: when walking, the body falls to the side of hemiplegia; 4 points: unable to walk spontaneously, loss of consciousness; 5 points: death. One-to-three points were selected as successful samples of cerebral ischemia. (2) The scope and area of infarction after cerebral ischemia reperfusion were determined by TTC staining, and the infarction area of each group was calculated by Image analysis software. After staining, the nonischemic area was rose-red and the infarction area was white. (3) The dry and wet method was used to the measurement of brain water content on cerebral ischemia side. Cerebral ischemia reperfusion can induce significant increase of brain water content and lead to cerebral edema in rats. Moisture content of brain tissue = (wet brain weight − dry brain weight)/wet brain weight × 100%.

Apparatus and Operation Conditions
For the LC separation, an Ultimate3000 HPLC (Thermo Scientific, San Jose, CA, USA) equipped with a quaternary pump, a cooling autosampler, a column temperature compartment, and UV detector were utilized. An XBridgeTMC18 (2.1 mm × 150 mm, 5 µm) was employed at 35 • C. The mobile phase consisted of methanol (A) and 3 mmol/L ammonium acetate (B), which was delivered at a flow rate of 0.3 mL/min. The applied gradient elution is as follows; 30-60% A from 0 to 2 min, 60-95% A from 2 to 10 min, 95-60% A from 10 to 11 min, 60-30% A from 11 to 12 min, and 30% A from 12 to 13 min. The injection volume was 5 mL.
For MS, a Q-Orbitrap mass spectrometer (Thermo Scientific, San Jose, CA, USA) equipped with a heat electrospray ionization (HESI) operated in the full-scan mode. The sheath gas, auxiliary gas, and sweep gas pressure was 35 bar, 10 bar, and 10 bar, respectively. Positive and negative ions were simultaneous scanned and the spray voltage was set at +3.5 kV or −2.

Plasma Sample Preparation
Rat plasma samples were directly subjected to protein precipitation [42,43] with methanol. Plasma samples (100 mL) were spiked with 100 mL of IS (0.916 mg/mL) and 200 mL of methanol, vortexed for 15 min, and centrifuged at 15,000 r/min for 20 min, and then the supernatant was transferred to a clean test tube, evaporated to dryness at 40 • C under a flow of nitrogen gas, the residue was reconstituted in 100 mL of methanol-water (30:70, v/v) and subsequently, a 5 mL aliquot was injected for HPLC-MS/MS measurement. Because the concentration of rhein is too high, it is not within the range of our standard curve linear range. In the actual operation of measuring the blood concentration of rhein, we diluted the plasma samples after the above treatment with methanol at ratio 1:50 (v:v).

Standard Solution Preparation and Sample Quality Control
Five reference products (aloe-emodin, emodin, rhein, chrysophanol, and physcion) of rhubarb were accurately weighed, and methanol was used as solvent to prepare stock solutions with appropriate concentration. The five stock solutions were moderately mixed to afford a final mixed standard solution.
A series of working solutions of these five analytes were freshly prepared by diluting the mixed standard solution with methanol gradiently. 1,8-Dihydroxyanthraquinone, an IS, was prepared at a concentration of 0.916 µg/mL by methanol. Calibration samples were prepared by spiking 100 µL working solutions of the corresponding concentrations; 100 µL methanol and 100 µL of IS solution to 100 µL of blank rat plasma. Then, calibration samples were pretreated according to the "4.3.2 Plasma sample preparation". Low-, medium-, and high-quality control (QC) samples were also prepared in the same way; the medium concentrations were of aloe-emodin 93.2 ng/mL; emodin 77.6 ng/mL; rhein 82.8 ng/mL; chrysophanol 85.6 ng/mL; and physcion 86.4 ng/mL, and were used to measure the precision and accuracy of our methods. All stock solutions were refrigerated at 4 • C.

Administration and Sample Collection
After seven days of acclimatization, 150 rats were randomly divided into groups of 10 each. And those subjects were fasted but free to drink water for 12 h before building the cerebral ischemia model. The rats were administered when the cerebral ischemia model established for 12 h. The administered composition and dose of the 15-group pharmacokinetic study were the same as those of the 15 experimental groups in the pharmacodynamic study (4.2.1). Blood was collected from the tail vein before and after administration at 0.083, 0.25, 0.75, 1, 2, 4, 6, 8, 12, 24 h, respectively, and collected in the heparin centrifuge tube. The blood was then centrifuged in a refrigerated centrifuge for 20 min (6000 r/min) to obtain plasma. Plasma was stored in a −80 • C refrigerator for cryopreservation.

Statistical Methods
SPSS19.0 analysis software was adopted for data analysis. The results were expressed in terms of means ± standard deviation and one-way ANOVA was used to analyze the data differences between the groups by some standards; p < 0.05 indicates statistical significance, p < 0.01 indicates significant difference.

Conclusions
In each group rats were administered with one of the five anthraquinone aglycones; the anticerebral ischemia effect of chrysophanol was the most obvious. In the compatibility groups, the anticerebral ischemia effect of aloe-emodin + physcion group was the most obvious. After co-administration with Aloe-emodin, the anticerebral ischemia effect of emodin, rhein, and chrysophanol was weakened. Therefore, when treating cerebral ischemic injury, it was not recommended to administer aloe-emodin with emodin, rhein, and chrysophanol. The pharmacokinetic interaction between the five anthraquinones aglycones was intricate. Emodin, rhein, chrysophanol, and physcion increased the absorption and slowed the elimination of aloe-emodin, while aloe-emodin reduced the AUC 0-t values of the above four anthraquinones. Emodin, chrysophanol, and rhein have a similar pharmacokinetic effect on each other. It was speculated that the three anthraquinones had the same targets in anticerebral ischemia. In this study, the pharmacodynamics and pharmacokinetics of five rhubarb aglycones were systematically studied, and their pharmacodynamics and pharmacokinetic interaction were explored, which provided a theoretical basis for the study of the mechanism of action.
Supplementary Materials: The following are available online, Figure S1. Chemical formula for five anthraquinones (Rhein, Emodin, Chrysophanol, Aloe-emodin, Physcion) and IS. Figure S2. Six-channel chromatogram of five anthraquinones (Aloe-emodin, Aloe, Rhein, Chrysophanol, Pphyscion). Table S1. Calibration curves for the 5 anthraquinones in rat plasma. Table S2. Intra-day and inter-day precision and accuracy for determination of the 5 anthraquinones in rat plasma. Table S3. Matrix effects and extraction recovery of the 5 anthraquinones in rat plasma. Table S4. Stability of the 5 anthraquinones in rat plasma.  : 81274139; 30873265). Both of these laboratories contributed to the study design, research, and interpretation of data, and the writing, reviewing, and approving of the publication.