Skip to main content

Toxicity and toxicokinetics of the ethanol extract of Zuojin formula



Zuojin formula, a traditional Chinese medicine, comprises Coptis chinensis and Evodia rutaecarpa. In our previous study, the total alkaloid extract from Zuojin formula (TAZF) showed potent and improved efficacy. However, its safety and toxicokinetics remain unknown. The objective of this study was to evaluate the safety of repeated administrations of TAZF and investigate the internal exposure of the main components and its relationship with toxic symptoms.


Sprague–Dawley rats were orally administered TAZF at 0.4, 1.2 and 3.7 g/kg for 28 days, which was followed by a 14-day recovery period. The toxic effects were evaluated weekly by assessing body weight changes, food intake, blood biochemistry and haematological indices, organ weights and histological changes. A total of eight components were detected, including berberine, coptisine, epiberberine, palmatine, jatrorrhizine, columbamine, evodiamine, and rutaecarpine. The toxicokinetic profiles of the eight components were investigated after single and repeated administrations. Linear mixed effect models were applied to analyse the associations between internal exposure and toxic symptoms. Network pharmacology analysis was applied to explore the potential toxic mechanisms.


Compared with the vehicle group, the rats in the low- and medium-dose groups did not show noticeable abnormal changes, while rats in the high-dose group exhibited inhibition of weight gain, a slight reduction in food consumption, abdominal bloating and atrophy of the splenic white pulp during drug administration. The concentration of berberine in plasma was the highest among all compounds. Epiberberine was found to be associated with the inhibition of weight gain. Network pharmacology analysis suggested that the alkaloids might cause abdominal bloating by affecting the proliferation of smooth muscle cells. The benchmark dose lower confidence limits (based on body weight inhibition) of TAZF were 1.27 g/kg (male) and 1.91 g/kg (female).


TAZF has no notable liver or kidney toxicity but carries risks of gastrointestinal and immune toxicity at high doses. Alkaloids from Coptis chinensis are the main plasma components related to the toxic effects of TAZF.

Peer Review reports


Zuojin formula is a traditional Chinese medicine (TCM) comprised of Coptis chinensis (the dry roots of Coptis chinensis Franch.) and Evodia rutaecarpa (the dry mature fruits of Euodia rutaecarpa (Juss.) Benth.) at a ratio of 6:1. Zuojin formula is famous for its therapeutic effect against some digestive tract diseases, such as gastric acid reflux, stomachache, diarrhoea, and gastritis [1,2,3,4]. Zuojin formula can also exhibit antitumour and antidepressive effects [1, 5,6,7,8,9]. As a traditional formula with a long history, Zuojin Pill has been regarded as a safe and nontoxic drug, and there are few reports of adverse reactions.

However, Evodia rutaecarpa, one of the ingredients in Zuojin formula, has been found to have noticeable hepatotoxicity. Both water and alcohol extracts of Evodia rutaecarpa can lead to liver injury in a time- and dose-dependent manner [10,11,12,13]. Alkaloids, such as evodiamine and rutaecarpine, and limonoids in Evodia rutaecarpa are the main active components [14]. Since Zuojin Pill consists of these potential hepatotoxic components, it has a potential risk of toxicity.

In our previous study, the alkaloids in Coptis chinensis and Evodia rutaecarpa were enriched after being extracted and eluted with different concentrations of ethanol. The total alkaloid extract of Zuojin formula (TAZF) had better efficacy than Zuojin Pill, as TAZF administered at a dosage of 1.37 g/day (humans) or 0.12 g/kg (rats) exerted the same or better treatment effects against gastritis as 6 ~ 12 g/day Zuojin Pill (unpublished data). Our preliminary study found that the LD50 of TAZF in rats was over 5 g/kg, rating this formula as nontoxic (unpublished data). However, the toxicity of repeated TAZF treatment remains unknown. In the present study, we evaluated the safety and toxicokinetic profiles of TAZF in a 28-day oral toxicity study in rats. Furthermore, we intended to find the links between exposure and toxic symptoms using statistical modelling.


Materials, reagents and animals

Carbamazepine (CAS No.: 298–46-4), berberine (CAS No.: 2086-83-1), coptisine (CAS No.: 6020-18-4), epiberberine (CAS No.: 6873-09-2), palmatine (CAS No.: 3486-67-7), jatrorrhizine (CAS No.: 3621-38-3), columbamine (CAS No.: 483–34-1), evodiamine (CAS No.: 518–17-2) and rutaecarpine (CAS No.: 20575–76-2), all > 98% purity, were purchased from Shanghai Yuanye Bio-Technology Company. Acetonitrile and methanol (LC grade) were purchased from Fisher Chemical, USA. Formaldehyde (AR) was purchased from Beijing Tong Guang Fine Chemicals Company. Isoflurane was provided by HARVEYBIO Company.

TAZF was provided by the State Key Laboratory of Natural and Biomimetic Drugs, School of Pharmaceutical Sciences, Peking University. The concentrations of the compounds in TAZF were measured using high-performance liquid chromatography (HPLC), with berberine accounting for the majority of TAZF at a ratio of 14.73% (g/g), followed by palmatine (3.14%), coptisine (2.52%), evodiamine (1.41%), epiberberine (1.30%), rutaecarpine (1.27%), columbamine (0.90%) and jatrorrhizine (0.47%).

Repeated administration toxicity test

Animal administration and treatments

Sprague–Dawley (SD) rats (5 ~ 6 weeks old, male: 146.5 ± 7.9 g, female: 141.1 ± 9.1 g) were obtained from the Department of Laboratory Animal Science, Peking University Health Science Center. Animals were acclimated in environmentally controlled rooms (temperature of 20 ~ 25 °C; relative humidity of 50 ~ 55%; 12-h light/dark cycle) with food and water available ad libitum for approximately 1 week prior to the initiation of dosing. All experimental procedures were carried out following the National Institutes of Health Guide for the Care and Use of Laboratory Animals. Animal experimental protocols were approved by the Peking University Medical Ethical Review Committee (No. LA2018313).

Rats were randomly allocated to control and treated groups (16/group/sex), and administered the appropriate material via gavage for 28 days, which could support a 2-week clinical therapy regimen according to ICH guidelines [15,16,17]. Rats in the vehicle control group were treated with deionized water, while those in the treated groups were administered TAZF at dosages of 3.7 g/kg (high dose), 1.2 g/kg (medium dose), and 0.4 g/kg (low dose), which were close to 30-, 10-, and 3-fold the therapeutic dose in rats. The possible dosages of the main components are listed in Table 1.

Table 1 The dosage of main components in TAZF administered by gavage in rats

As shown in Fig. 1, blood was collected on the 7th, 14th, and 21st days for biochemical analysis. On the 29th day, twenty rats from each group were sacrificed under isoflurane anaesthesia. Blood was collected, and the main organs were weighed and fixed in formalin. Organ coefficients (OCs) (\(\mathrm{OCs}=\frac{\mathrm{organ}\ \mathrm{weight}}{\mathrm{body}\ \mathrm{weight}}\ast 100\%\)) were also calculated. The remaining 12 rats were observed for another 14 days (the recovery period).

Fig. 1
figure 1

Animal treatment schedule

Blood biochemistry measurements

Serum samples, separated from the blood collected on the 7th, 14th, 21st, and 28th days, were analysed with an AU2700 automatic biochemical analyser for the following biochemical indicators: liver function indices (levels of aspartate aminotransferase (AST), alanine aminotransferase (ALT), and alkaline phosphatase (ALP)), kidney function indices (urea nitrogen (BUN), creatinine (CRE)), blood lipid indices (total cholesterol (TC), triglycerides (TGs), high-density lipoprotein (HDL), low-density lipoprotein (LDL)), total protein (TP), albumin (Alb), blood glucose (Glu), total bilirubin (TBIL), creatine phosphokinase (CK)), and blood ion levels (sodium ions (Na+), potassium ions (K+) and chloride ions (Cl)).

Haematological examination

Blood samples collected on the last day of dosing and in the recovery period were analysed with a MEK-6318 K Automatic blood cell analyser for the following haematological indices: red blood cell count (RBC), haemoglobin (HGB), haematocrit (HCT), mean red blood cell volume (MCV), mean red blood cell haemoglobin (MCH), mean red blood cell haemoglobin concentration (MCHC), reticulocyte count (RET), white blood cell count (WBC) and classification (MONO%, LYMPH%, BASO%, EO%, NEUT%), platelet count (PLT) and prothrombin time (PT).

Histological observations

After sacrifice, the following main organs were fixed in formalin and made into pathological sections: brain, pituitary, thyroid, gastrointestinal tract, mesenteric lymph nodes, heart, liver, kidney, spleen, lung, testis, epididymis, ovary, and uterus. All sections were stained with haematoxylin and eosin (H&E) for microscopic examination.

Safety evaluation

The benchmark dose (BMD) and benchmark dose lower confidence limits (BMDLs) were calculated to assess the safety of TAZF. Toxic symptoms and indices or sensitive changes were regarded as the response. BMD analysis was performed with the Bayesian BMD platform ( [18] and the benchmark response (BMR) was set to 0.1.

Toxicokinetics test in the satellite groups

Animals in each satellite group

SD rats (5 ~ 6 weeks, 3/group/sex) were randomly allocated to three treatment groups and given the same dosage as those in the repeated toxicity experiment for 28 days. On the 1st and 28th days, blood samples were collected at 0.25, 0.5, 1, 2, 4, 6, 10 and 24 h after administration. Eight components in TAZF were detected: berberine, coptisine, epiberberine, palmatine, jatrorrhizine, columbamine, evodiamine and rutaecarpine. Carbamazepine was chosen as the internal standard (IS).

Plasma sample preparation

The blood samples were centrifuged to obtain plasma. Then, 300 μL of acetonitrile containing 40 ng of IS was added to 100 μL of plasma. After vortexing for 1 min, the mixture was centrifuged at 17,000 rpm for 15 min at 4 °C. The supernatant of each sample was injected into a UPLC–MS/MS system for determination. The injection volume was 1 μL.

UPLC–MS analysis

UPLC–MS analysis was performed on a Thermo Fisher Scientific Ultimate 3000 UPLC system (San Jose, USA) and an API 4000 QTRAP mass spectrometer (Foster City, USA) equipped with an ESI source. The mobile phase consisted of 2 solvents: (A) water containing 0.1% formic acid and (B) 100% acetonitrile. The elution gradient of the mobile phase is shown in Table 2. For MS, multiple reaction monitoring (MRM) in positive ion mode was used as the scanning mode for the main components and IS in the plasma samples. The dwell time (DT), collision energy (CE), quadrupole 1 pre-rod bias (Q1) and quadrupole 3 pre-rod bias (Q3) are listed in Table 3. The other parameters were optimized as follows: the curtain gas was 10 psi, ion source gas 1 and gas 2 were both 55 psi, the ion spray voltage was 5500 V, and the interface temperature was 600 °C.

Table 2 The gradient of mobile phase A and B in UPLC-MS
Table 3 Acquisition parameters in MRM mode

Analytical method validation

Standard curves for berberine were constructed with various concentrations: 0.5, 1, 2, 4, 10, 20, 50, 100, 200 and 400 ng/ml. The concentrations of the other components were 0.5, 1, 2, 4, 5, 7.5, 10, 25, 50 and 100 ng/ml due to their lower concentrations in plasma. The regression coefficient (r) of each equation was calculated, and the values should be over 0.99.

Recovery was assessed by analysing quality control (QC) samples at three concentrations: 1.5, 100 and 300 ng/ml for berberine and 1.5, 25 and 75 ng/ml for the other compounds. Recovery was determined by the following 2 equations, and RE% should be within 20%.

$$\mathrm{recovery}\ \mathrm{rate}=\frac{\mathrm{Measured}\ \mathrm{value}}{\mathrm{QC}\ \mathrm{concentration}}\times 100\%$$
$$\mathrm{RE}\%=\frac{\left|\mathrm{Measured}\ \mathrm{value}-\mathrm{QC}\ \mathrm{concentration}\right|}{\mathrm{QC}\ \mathrm{concentration}}\times 100\%$$

Toxicokinetic data analysis

The toxicokinetic parameters were calculated by PKSolver software [19], and the AUC0-t values after 28 days of repeated administration were used to evaluate the internal exposure.

Associations between internal exposure to the components in TAZF and changes in rat body weight

Linear mixed effect models were applied to explore the relationship between internal exposure to the components and changes in the body weights (BWs) of the rats. The models of internal exposure and changes in BW can be described as follows:

Basic model:

log(BW) = gender + weight0 + days + (1| days)


Total model:

log(BW) = gender + weight0 + days + component1 + (component1| days) + component2 + (component2| days) + …


Specific model:

log(BW) = gender + weight0 + days + component + (component| days)


where “ log(BW) ” represents the logarithm of the BWs of the rats in each satellite group; “ gender ”, “ weight0 ” and “ days ” represent the gender, basic BW and days of treatment, respectively, “ component ”, “ component1 ”, “ component2 ”, etc., represent the internal exposure to the components in plasma (berberine, epiberberine, etc.), and “ (1| days) ” and “ (component| days) ” represent the random effects in the models. The Akaike information criterion (AIC) values from the models were compared. When the AIC values of the component were notably lower than those of the basic model, this component was considered to influence BW changes.

Network pharmacology analysis

Network pharmacology analysis was performed to investigate the potential toxic mechanisms of TAZF. The components detected in plasma were searched in the Comparative Toxicogenomics Database (CTD, for related target genes. Studies have indicated that drug-induced paralytic ileus might cause abdominal bloating [20]. Therefore, “Ileus” and “intestinal obstruction” were chosen as disease/toxic effect keywords, and the related genes in the DisGeNET database ( and MalaCards database ( were gathered. Genes that interacted with both components and diseases were selected to perform protein–protein interaction (PPI) network analysis with STRING ( Gene Ontology enrichment analysis was performed using the clusterProfiler package in R 4.1.2.

Statistical analysis

The results are presented as the mean ± S.D.. R software (4.1.2) was applied to perform t tests and ANOVA. Pictures were constructed in GraphPad Prism 7. A value of P < 0.05 was considered statistically significant.


Repeated toxicity study

Clinical observations

Compared with the vehicle control group, no apparent changes were found in rats in the low- or medium-dose groups. In the high-dose group, asthma and lung rales were observed and lasted for 3 to 4 days after treatment with TAZF for approximately 1 week. Some deaths (male: 5/16; female: 7/16) occurred after several days of abdominal distention and an obvious inhibition of BW gain. Anatomical examinations showed that the intestines were full of air, with the accumulation of deep dark mixtures in the terminal caecum. Asthma and lung rales were also observed in rats receiving the high dose on the 7th day of the recovery period. Fortunately, no rats died during this period.

Figure 2 shows the changes in rat BW and food intake. Compared with the vehicle control group, male rats in the high-dose group (Fig. 2A) showed apparent BW gain and food intake reduction from the second week of dosing. The same decline was also observed in females (Fig. 2B). Food intake in the treated groups was also less than that in the vehicle control group (Fig. 2C and D).

Fig. 2
figure 2

Body weight changes in rats (A male, B female) and food intake (C male, D female) during and after oral administration of TAZF for 28 days. Compared with the vehicle control group, *: P < 0.05; **: P < 0.01; ***: P < 0.001

Organ coefficients

Compared with the vehicle control group, the OCs of the kidneys, lungs, spleens, hearts, thymus, ovary, and uterus did not show significant differences (Tables S1 and S2). As shown in Fig. 3A, the liver OCs of the rats in the high-dose group was significantly increased, and this difference disappeared after stopping drug administration. The brain and adrenal gland in the male rats in the high-dose group (Fig. 3B and C) had noticeable weight gain. Male reproductive organs, such as the testis (Fig. 3D) and epididymis (Fig. 3E), were also heavier in the rats in the high-dose group than in those in the control group.

Fig. 3
figure 3

Organ coefficients of the liver (A), brain (B), adrenal gland (C), testis (D) and epididymis (E) in rats during (n = 10) and after (n = 6) TAZF administration. Compared with the vehicle control group, *** P < 0.001; **: P < 0.01; *: P < 0.05

Blood biochemistry

Figure 4 presents the results of blood chemistry during drug administration. Compared with the vehicle control group, liver function indices (ALT, AST, ALP and TBIL levels) did not show significant toxic changes. The decreases in TC and TGs and the increase in HDL started as early as the first week, suggesting a blood lipid decreasing effect. Although BUN and CRE levels increased in several groups, the changes were in the normal reference range for SPF rats [21], which meant that the administration of TAZF was not harmful to the kidney. The results of the plasma protein indices and blood ion levels are presented in Figure S1.

Fig. 4
figure 4

Results of blood biochemistry during TAZF administration. Liver function indices: (A) ALT; (B) AST; (C) ALP; (D) TBIL. Blood lipid indices: (E) TC; (F) TGs; (G) HDL; (H) LDL. Kidney function indices: (I) BUN; (J) Cre. Compared with the vehicle control group, ***: P < 0.001; **: P < 0.01; *: P < 0.05

The serum biochemical levels during the recovery period are presented in Figure S2. Most changes, except those of TGs and LDL, which continued to decrease, disappeared after stopping drug administration.

Haematological examination

The results of the haematological examination are shown in Tables S3 and S4. Compared with the vehicle control group, the female rats had some changes in RBC and WBC after administration of TAZF. During treatment, the MCV and MCHC in the female rats in the high-dose group increased, while the MCH decreased. The classification of WBCs changed along with an increase in MONO% and a decrease in EO% during the recovery period, indicating that TAZF may affect the distribution of immune cells in peripheral blood.

Histological examination

All of the main organs of the rats in the control and high-dose groups were inspected by histological examination. As shown in Fig. 5, apparent atrophy of splenic white pulp was observed in rats receiving the high dose of TAZF, along with a decrease in lymphocyte density and proliferation of sinus tissue cells. Such changes occurred in 60% of the male (6/10) and 44% of the female rats (4/9) in the high-dose group. No other significant changes were found in the other organs, including the kidney, brain, liver, thymus, lung, male gonad, and digestive tract. These histological changes are presented in Figures S3 and S4.

Fig. 5
figure 5

Microphotographs of splenic white pulp atrophy in rats in the high-dose group. (A) control group-100×; (B) high-dose group-100×; (C) control group-200×; (D) high-dose group-200×. The border of the splenic white pulp in the control group was clear and the shape was complete, while the splenic white pulp from the rats in the treatment groups had unclear borders with an incomplete shape. The splenic white pulp area and lymphocyte density decreased, and sinus tissue cells proliferated

Safety evaluation

BW inhibition was chosen as a measurable toxic index to calculate the benchmark dose of TAZF. As shown in Table 4, the BMDLs of TAZF were 1.27 g/kg for males and 1.91 g/kg for females.

Table 4 BMD and BMDL values of body weights changes in rats

Toxicokinetics test

Analytical method validation

The linearity regression and linear range results are listed in Table S5. The calibration curves of the components presented satisfactory linearity with r values greater than 0.99. Table S6 presents the sensitivity and recovery results. All the RE% values of the components, except for 1.2 ng/ml evodiamine, showed an acceptable value under 20%, suggesting that these results were accurate.

Toxicokinetic profiles

The levels of evodiamine and rutaecarpine in plasma were under the detection limit (0.5 ng/ml). Figure S5 presents the toxicokinetic profiles of the other six components (berberine, epiberberine, coptisine, palmatine, jatrorrhizine and columbamine) after single and repeated administrations of TAZF. As shown in Fig. 6 and Table 5, after a 28 days of administration, the systemic exposure to berberine, coptisine, palmatine, jatrorrhizine, and columbamine was much lower than that after a single administration. The concentration of berberine in plasma was the highest, suggesting that berberine may play an essential role in the efficacy of TAZF. Most components had higher AUC0-t values in males than in females, and the former also showed significant dose-dependent relationships with the exception of berberine and columbamine on the 28th day. The 28-day AUCs of epiberberine in the high-dose group were 1 to 2 times greater than those of the first day (Fig. 6A and B), indicating that repeated administration might affect the elimination process. The Cmax of most components increased less than proportionally with the increase in dose. Most components had the highest plasma concentration in the medium-dose group on the first day (Fig. 6C and D). Other toxicokinetic parameters are presented in Table S7.

Fig. 6
figure 6

AUC0-t((A), (B)) and Cmax (ng/ml) ((C), (D)) of berberine, epiberberine, coptisine, palmatine, jatrorrhizine and columbamine on the first and 28th day of treatment in male ((A), (C)) and female ((B), (D)) rats

Table 5 Toxicokinetic parameters of 6 components in rat plasma on the first and 28th day treated with TAZF

Relationship between internal exposure to the components and changes in BW

Table 6 shows the AIC values of the linear mixed effect models. Except for the total model, only the model of epiberberine had a lower AIC value than the basic model, indicating that the internal exposure of epiberberine was responsible for the changes in BW in the treated groups. The coefficient of this model was < 0, suggesting that as epiberberine exposure increased, BW increased more slowly.

Table 6 AIC values of linear mixed effect models

Network pharmacology analysis

A total of 271 genes were included in the TAZF component target set, and 136 were included in the ileus-related gene set. Fourteen genes were found in both sets: ACHE, ALB, GPT, NGF, CCL2, HMOX1, IL10, IL1B, IL6, MMP9, PTGS2, TNF, CTNNB1, and MYC. The PPI network shown in Fig. 7 suggests that these genes are related. Figure 8 presents the results of GO analysis. These results showed that target genes might affect smooth muscle cell proliferation and that the cell membrane could be the action site.

Fig. 7
figure 7

Protein–protein interaction network of target genes from network pharmacology analysis

Fig. 8
figure 8

Results of molecular function (A), cellular component (B) and biological process (C) from GO analysis


In this study, we tried to assess the safety of TAZF via a repeated toxicity experiment. Although TAZF does not show acute toxicity in rats (unpublished data), it can cause a few toxic symptoms after repeated administration. These toxic symptoms appeared in rats in the high-dose (3.7 g/kg) group. The most severe symptom was abdominal bloating, which might lead to death. BW gain inhibition, food intake reduction, and transitory asthma were also observed in our study. Furthermore, atrophy of the splenic white pulp appeared in rats in the high-dose group.

Although Evodia rutaecarpa has been regarded as a hepatotoxic ingredient, TAZF treatment did not cause any hepatotoxic symptoms, except an increase in the liver organ coefficient in the high-dose group. The concentrations of evodiamine and rutaecarpine, the potential hepatotoxic components of Evodia rutaecarpa, were both low in the plasma. This may explain why TAZF is not hepatotoxic. However, the internal exposure of berberine in this study was much higher than that in other studies [22,23,24,25,26]. A previous study administered pure berberine at a dose of 90 mg/kg to male rats and found that the AUC0–60 was 305.8 ng/ml*h, while the Cmax was 28.2 ng/ml [23]. The AUC0–24 and Cmax in our study were 1222.27 ng/ml*h and 221.6 ng/ml, respectively, when rats were administered TAZF with a dose of berberine of 58.92 mg/kg. Female SD rats were orally administered pure berberine (90 mg/kg), and these AUC0–36 and Cmax values were 88.8 ng/ml*h and 29.2 ng/ml, respectively, while the AUC0–24 and Cmax in our study were 121.57 ng/ml*h and 148.77 ng/ml, respectively, when the female rats were administered 58.92 mg/kg berberine in TAZF. Berberine and other alkaloids, including jatrorrhizine and palmatine, can be metabolized by CYP3A4, CYP1A2 and CYP2D6 [27,28,29,30]. These enzymes can also be inhibited by epiberberine and berberine [31,32,33,34]. Thus, there might be drug–drug interactions (DDIs) between these alkaloids, which can alter the internal exposure of a single component. More studies are needed to clarify this process.

A sex difference was also observed for the plasma concentrations of the compounds. The gender difference effects were significant regarding the AUCs rather than the Cmax values [23, 25], indicating that metabolic factors might not be the primary reason for these differences. There were several peaks in the concentration-time profiles, and these multiple peaks played a major role in the sex differences in the pharmacokinetic profiles. Multiple blood concentration peaks from berberine have been reported in male rats [35]. A possible cause of multiple peaks is enterohepatic circulation, as enterohepatic circulation of berberine in rats has been reported [36]. However, no systematic research on enterohepatic circulation and sex differences among these compounds has been performed thus far. Much work needs to be done to investigate the multiple peaks and gender differences.

BW inhibition was chosen as the BMR toxic effect because the blood biochemistry and haematological examination did not exhibit noticeable toxic changes. However, BW can be easily affected by food consumption. Figure 2 reveals the connection between TAZF dose, body weight and food consumption, suggesting that BW could be influenced by a decrease in food consumption. Additionally, BMDLs were calculated based on the dose of TAZF, which contains multiple compounds, and applied to evaluate the safety overall. Although some compounds, such as berberine, exhibited a noticeable difference in exposure (AUC) based on sex, the influence of these compounds on the toxicity of TAZF could be complicated. These might be reasons why the AUC of a single compound showed a sex difference, while the BMDL values did not. Further studies on the relationship between multi-compound exposure and toxic symptoms are needed.

Smooth muscle cells are some of the main cells that make up the digestive tract and are related to peristalsis [20]. Network pharmacology analysis also suggested that these components could affect the proliferation of smooth muscle cells. The effect of berberine on smooth muscle cells has been confirmed by previous studies [37], and may contribute to abdominal bloating symptoms. Berberine is an effective antidiarrhoeal drug [2, 38, 39]. A previous study found that administration of berberine at dosages of 50 and 100 mg/kg for 14 days can improve diarrhoea-predominant irritable bowel syndrome in rats [40]. This effect was also observed in rats with thyroid hormone-induced diarrhoea after 7 days of 60 mg/kg berberine oral gavage [41]. The high-dose TAZF group in the present study received a similar dose of berberine, but this compound was administered for a longer period (28 days), resulting in higher internal exposure to alkaloids. This high level of berberine might contribute to the adverse effects on the intestine found in this study.

The atrophy of splenic white pulp observed in the high-dose TAZF group might also be related to the alkaloids in TAZF. Berberine has been found to have anti-inflammatory activity and immunosuppressive effects [42, 43]. Recent studies showed that mice administered berberine via intraperitoneal injection at a dose of 1 mg/kg had reduced levels of CD4+ Th and CD4+CXCR5+ Tfh cells and an increased level of Foxp3+ Tregs [44, 45]. Berberine administered to mice via intraperitoneal injection (10 mg/kg) had a strong immunosuppressive effect and reduced spleen weight in addition to changes in the classification of immune cells in the spleen [46]. These findings may partly explain the atrophy of splenic white pulp observed in our study. Further studies are still needed.


In the present study, we evaluated the safety of TAZF via 28 days of repeated administration. The BMDL of TAZF was found to be 1.27 g/kg (male) and 1.91 g/kg (female). TAZF has no noticeable liver or kidney toxicity but carries risks of gastrointestinal and immune toxicity at a high dose. Alkaloids from Coptis chinensis are the main plasma components related to TAZF’s toxic effects.

Availability of data and materials

All data generated or analysed during this study are included in this published article and its supplementary information files.


  1. Zhu C, Shi W, Chen Y, Zhu X. Clinical application and experimental research progress of Zuojin pill. J Pract Tradit Chin Med. 2019;35(5):626–9.

    Google Scholar 

  2. Guan M, Xu W, Li H. Advances in modern pharmacology of Zuojin pills. Acta Chin Med Pharmacol. 2020;48(5):78–81.

    Google Scholar 

  3. Wang Y, Ma C. Research Progress on treatment of gastric ulcer with Zuojin pill. J Pract Tradit Chin Int Med. 2021;35(5):103–6.

    Google Scholar 

  4. Wang J, Zhang T, Zhu L, Ma C, Wang S. Anti-ulcerogenic effect of Zuojin pill against ethanol-induced acute gastric lesion in animal models. J Ethnopharmacol. 2015;173:459–67.

    Article  CAS  PubMed  Google Scholar 

  5. Ye N, Liu X, Li Q. Experimental research progress of Zuojinwan and its active constituents in anti gastric cancer. China J Tradit Chin Med Pharm. 2014;29(4):1162–4.

    CAS  Google Scholar 

  6. Zhang J, Li Q. Progress of antitumor mechanism of Zuojin pills. J Changchun Univ Chin Med. 2018;34(3):609–12.

    Google Scholar 

  7. Wang QS, Ding SL, Mao HP, Cui YL, Qi XJ. Antidepressant-like effect of ethanol extract from Zuojin pill, containing two herbal drugs of Rhizoma Coptidis and Fructus Evodiae, is explained by modulating the monoaminergic neurotransmitter system in mice. J Ethnopharmacol. 2013;148(2):603–9.

    Article  CAS  PubMed  Google Scholar 

  8. Wei Y, Wang R, Ren S, Liu X, Jing M, Li R, et al. Zuojin pill ameliorates inflammation in indomethacin-induced gastric injury via inhibition of MAPK pathway. J Ethnopharmacol. 2021;275:114103.

    Article  CAS  PubMed  Google Scholar 

  9. Wang K, Miao X, Kong F, Huang S, Mo J, Jin C, et al. Integrating network pharmacology and experimental verification to explore the mechanism of effect of Zuojin pills in pancreatic Cancer treatment. Drug Design Dev Ther. 2021;15:3749–64.

    Article  Google Scholar 

  10. Zhang S, Liu S, Zhong Z, Chen B, Wei H. The relationship of“dose -time-toxicity”of different processed products of Evodia rutaecarpa ( Juss.) Benth on hepatotoxicity in mice. Lishizhen Med Mater Med Res. 2018;29(4):881–4.

    Google Scholar 

  11. Sun R, Huang W, Lv L. Study on the “dose-time-toxicity” relationship of Evodia rutaecarpa volatile oil single administration on liver toxicity in mice. Pharmacol Clin Chin Mater Med. 2012;38(3):55–8.

    Google Scholar 

  12. Huang W, Sun R, Lv L, Bao Z, Wang H. Study on the “dose-time-toxicity” relationship of single administration of Evodia rutaecarpa water extract on liver toxicity in mice. Pharmacol Clin Chin Mater Med. 2012;28(5):96–9.

    Google Scholar 

  13. Huang W, Li X, Sun R. Study on the “dose-time-toxicity” relationship of multiple administration of aqueous extract of Evodia rutaecarpa on liver toxicity in mice. China J Chin Mater Med. 2012;37(15):2223–7.

    Google Scholar 

  14. Liu L, Zhang X, Xu J, Zhang J, Zhang T, Chen C, et al. Chemical components and pharmacological action for Euodiae Fructus and predictive analysis on its Q-marker. Chin Tradit Herbal Drugs. 2020;51(10):2689–701.

    Google Scholar 

  15. ICH M3 (R2) Nonclinical Safety Studies for the Conduct of Human Clinical Trials and Marketing Authorization for Pharmaceuticals. 2010.

  16. CFDA [Z]GPT3-1 Long-term Toxicity Research of Traditional Chinese Medicine and Natural Drugs. 2005.

  17. ICH M3(R2) Guidance for Industry M3 Nonclinical Safety Studies for the Conduct of Human Clinical Trials for Pharmaceutical.  1997.

  18. Shao K, Shapiro AJ. A web-based system for Bayesian benchmark dose estimation. Environ Health Perspect. 2018;126(1):017002.

    Article  PubMed  PubMed Central  Google Scholar 

  19. Zhang Y, Zhou J, Huo M. A pharmacokinetic data analysis program based on Excel -- Writing and verification of pkslover 1.0. J Mathematical Med. 2007;20:58–61.

    CAS  Google Scholar 

  20. Qiang Q, Sheng L, Zhang B, Wang W. Classification of intestinal obstruction and its diagnostic methods. Western J Tradit Chin Med. 2015;28(10):140–2.

    Google Scholar 

  21. Fan J, Qin J, Liu C, Tang X, Tian H, Liu Y, et al. Establishment of normal reference range of blood hematological and biochemical indicators of SPF wistar rats. J Public Health Prev Med. 2010;21(6):50–2.

    Google Scholar 

  22. Buchanan B, Meng Q, Poulin MM, Zuccolo J, Azike CG, Gabriele J, et al. Comparative pharmacokinetics and safety assessment of transdermal berberine and dihydroberberine. PLoS One. 2018;13(3):e0194979.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  23. Li G, Yang F, Liu M, Su X, Zhao M, Zhao L. Development and application of a UPLC-MS/MS method for simultaneous determination of fenofibric acid and berberine in rat plasma: application to the drug-drug pharmacokinetic interaction study of fenofibrate combined with berberine after oral administration in rats. Biomed Chromatogr. 2016;30(7):1075–82.

    Article  CAS  PubMed  Google Scholar 

  24. Hu Y, Fan S, Liao X, Chen C, Su L, Li X. Pharmacokinetics, excretion of 8-cetylberberine and its main metabolites in rat urine. J Pharm Biomed Anal. 2017;132:195–206.

    Article  CAS  PubMed  Google Scholar 

  25. Liu M, Su X, Li G, Zhao G, Zhao L. Validated UPLC-MS/MS method for simultaneous determination of simvastatin, simvastatin hydroxy acid and berberine in rat plasma: application to the drug-drug pharmacokinetic interaction study of simvastatin combined with berberine after oral administration in rats. J Chromatogr B. 2015;1006:8–15.

    Article  CAS  Google Scholar 

  26. Lu T, Liang Y, Song J, Xie L, Wang GJ, Liu XD. Simultaneous determination of berberine and palmatine in rat plasma by HPLC-ESI-MS after oral administration of traditional Chinese medicinal preparation Huang-Lian-Jie-Du decoction and the pharmacokinetic application of the method. J Pharm Biomed Anal. 2006;40(5):1218–24.

    Article  CAS  PubMed  Google Scholar 

  27. Guo Y, Li F, Ma X, Cheng X, Zhou H, Klaassen CD. CYP2D plays a major role in berberine metabolism in liver of mice and humans. Xenobiotica. 2011;41(11):996–1005.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  28. Zhou H, Shi R, Ma B, Ma Y, Wang C, Wu D, et al. CYP450 1A2 and multiple UGT1A isoforms are responsible for jatrorrhizine metabolism in human liver microsomes. Biopharm Drug Dispos. 2013;34(3):176–85.

    Article  CAS  PubMed  Google Scholar 

  29. Xiao W, Deng Z, Lai C, Lu H, Huang M, Wen Y, et al. Inhibitory effect of ketoconazole, quinidine and 1-aminobenzotriazole on pharmacokinetics of l-tetrahydropalmatine and its metabolite in rats. Xenobiotica. 2021;51(4):447–54.

    Article  CAS  PubMed  Google Scholar 

  30. Vrba J, Papouskova B, Pyszkova M, Zatloukalova M, Lemr K, Ulrichova J, et al. Metabolism of palmatine by human hepatocytes and recombinant cytochromes P450. J Pharm Biomed Anal. 2015;102:193–8.

    Article  CAS  PubMed  Google Scholar 

  31. Zhao Y, Hellum BH, Liang A, Nilsen OG. The in vitro inhibition of human CYP1A2, CYP2D6 and CYP3A4 by tetrahydropalmatine, neferine and berberine. Phytother Res. 2012;26(2):277–83.

    Article  CAS  PubMed  Google Scholar 

  32. Lee SY, Jang H, Lee JY, Ma JY, Oh SJ, Kim SK. Inhibitory effects of Hwang-Ryun-Hae-Dok-Tang on cytochrome P450 in human liver microsomes. Xenobiotica. 2015;45(2):131–8.

    Article  CAS  PubMed  Google Scholar 

  33. Zhao X, Zhang JJ, Wang X, Bu XY, Lou YQ, Zhang GL. Effect of berberine on hepatocyte proliferation, inducible nitric oxide synthase expression, cytochrome P450 2E1 and 1A2 activities in diethylnitrosamine- and phenobarbital-treated rats. Biomed Pharmacother. 2008;62(9):567–72.

    Article  CAS  PubMed  Google Scholar 

  34. Chen N, Yang XY, Guo CE, Bi XN, Chen JH, Chen HY, et al. The oral bioavailability, excretion and cytochrome P450 inhibition properties of epiberberine: an in vivo and in vitro evaluation. Drug Design Dev Ther. 2018;12:57–65.

    Article  CAS  Google Scholar 

  35. Deng Y, Liao Q, Li S, Bi K, Pan B, Xie Z. Simultaneous determination of berberine, palmatine and jatrorrhizine by liquid chromatography-tandem mass spectrometry in rat plasma and its application in a pharmacokinetic study after oral administration of coptis-evodia herb couple. J Chromatogr B Analyt Technol Biomed Life Sci. 2008;863(2):195–205.

    Article  CAS  PubMed  Google Scholar 

  36. Qu Q, Zhang M. Pharmacokinetics of berberine. Chin Acad Med Magazine Organ. 2002;4:31–5.

    Google Scholar 

  37. Xu Z, Tang Z, Lin Y. Effect ofberberine on the function ofmyosin and on the contractility of intestinal and gastric smooth muscle. Prog Modern Biomed. 2009;9:3216–61.

    CAS  Google Scholar 

  38. Mirhadi E, Rezaee M, Malaekeh-Nikouei B. Nano strategies for berberine delivery, a natural alkaloid of Berberis. Biomed Pharmacother. 2018;104:465–73.

    Article  CAS  PubMed  Google Scholar 

  39. Hua X, Rui J, Ren B. Effects of traditional Chinese medicine Zuojin pill on gastrointestinal system. Drug Eval Res. 2017;40(2):190–5.

    Google Scholar 

  40. Lu Y, Huang J, Zhang Y, Huang Z, Yan W, Zhou T, et al. Therapeutic effects of Berberine hydrochloride on stress-induced diarrhea-predominant irritable bowel syndrome rats by inhibiting neurotransmission in colonic smooth muscle. Front Pharmacol. 2021;12:596686.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  41. Cheng ZF, Zhang YQ, Liu FC. Berberine against gastrointestinal peptides elevation and mucous secretion in hyperthyroid diarrheic rats. Regul Pept. 2009;155(1–3):145–9.

    Article  CAS  PubMed  Google Scholar 

  42. Ehteshamfar SM, Akhbari M, Afshari JT, Seyedi M, Nikfar B, Shapouri-Moghaddam A, et al. Anti-inflammatory and immune-modulatory impacts of berberine on activation of autoreactive T cells in autoimmune inflammation. J Cell Mol Med. 2020;24(23):13573–88.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  43. Zou K, Li Z, Zhang Y, Zhang HY, Li B, Zhu WL, et al. Advances in the study of berberine and its derivatives: a focus on anti-inflammatory and anti-tumor effects in the digestive system. Acta Pharmacol Sin. 2017;38(2):157–67.

    Article  CAS  PubMed  Google Scholar 

  44. Vita AA, Aljobaily H, Lyons DO, Pullen NA. Berberine delays onset of collagen-induced arthritis through T cell suppression. Int J Mol Sci. 2021;22(7):3522.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  45. Zhou BG, Liu FC, Zhao HM, Zhang XY, Wang HY, Liu DY. Regulatory effect of Zuojin pill on correlation with gut microbiota and Treg cells in DSS-induced colitis. J Ethnopharmacol. 2020;262:113211.

    Article  CAS  PubMed  Google Scholar 

  46. Mahmoudi M, Zamani Taghizadeh Rabe S, Balali-Mood M, Karimi G, Memar B, Rahnama M, et al. Immunotoxicity induced in mice by subacute exposure to berberine. J Immunotoxicol. 2016;13(2):255–62.

    Article  CAS  PubMed  Google Scholar 

Download references


The authors acknowledge the funding provided by the National Key Research and Development Program of China.


This work was supported by the National Key Research and Development Program of China (2018YFC1704506).

Author information

Authors and Affiliations



Shuo Wang: Investigation, Data curation, Formal analysis, Methodology, Visualization, Writing – original draft, Writing – review & editing; Tao Zhang: Methodology, Investigation, Writing – review & editing; Xiaoyan Liu: Investigation, Writing – review & editing; Zheng Yang: Investigation; Ludi Li: Investigation; Danping Shan: Investigation; Yadong Gao: Investigation; Yingzi Li: Investigation; Yanying Li: Investigation; Youbo Zhang: Conceptualization, Writing – review & editing; Qi Wang: Conceptualization, Supervision, Resources, Writing – review & editing. The author(s) read and approved the final manuscript.

Corresponding authors

Correspondence to Youbo Zhang or Qi Wang.

Ethics declarations

Ethics approval and consent to participate

Animal experimental procedures were approved by Peking University Medical Ethical Review Committee (No. LA2018313). Animals were cared for in accordance with the principles of the Guide to the Care and Use of Experimental Animals (Committee on Care and Use of Laboratory Animals, 1985) and were carried out in compliance with the ARRIVE guidelines.

Consent for publication

Not applicable.

Competing interests

The authors declare that there are no conflicts of interest.

Additional information

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Supplementary Information

Additional file 1: Figure S1.

Results of plasma protein indexes and blood ion levels during the administration. Compared with vehicle control group ***: P < 0.001; **: P < 0.01; *: P < 0.05. Figure S2. Results of blood chemistry after the oral administration of TAZF. Compared with vehicle control group, ***: P < 0.001; **: P < 0.01; *: P < 0.05. Figure S3. Microphotographs of liver, kidney and lung in control group and high-dose group. No obvious toxic change was found. Each figure is a representative photomicrograph from a rat in each group. Figure S4. Microphotographs of brain and digestive tract in control group and high-dose group. No obvious toxic change was found. Each figure is a representative photomicrograph from a rat in each group. Figure S5. Mean plasma concentration profiles of components in rat plasma after 1-day and 28-day oral administration of TAZF. Table S1. OCs of lung, spleen and heart in male rats during (n = 10) and after (n = 6) administration. Table S2. OCs of lung, spleen and heart in female rats during (n = 10) and after (n = 6) administration. Table S3. Results of hematological examination in male rats during (n = 10) and after (n = 6) administration. Table S4. Results of haematological examination in female rats during (n = 10) and after (n = 6) administration. Table S5. Linearity and linear range of the components. Table S6. Results of sensitivity and recovery. Table S7. Pharmacokinetic parameters of 6 components on the first and 28th day treated day in male rats. Table S8. Pharmacokinetic parameters of 6 components on the first and 28th day treated day in female rats.

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit The Creative Commons Public Domain Dedication waiver ( applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Wang, S., Zhang, T., Liu, X. et al. Toxicity and toxicokinetics of the ethanol extract of Zuojin formula. BMC Complement Med Ther 22, 220 (2022).

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: