Metabolomic Analysis Reveals the Therapeutic Effects of MBT1805, a Novel Pan-Peroxisome Proliferator-Activated Receptor Agonist, on α-Naphthylisothiocyanate-Induced Cholestasis in Mice

Background and Aims: Therapeutic drugs that are used to treat cholestatic liver disease are limited; however, the results of clinical trials on primary biliary cholangitis treatment targeting peroxisome proliferator-activated receptors (PPARs) are encouraging. In this study, we aimed to identify the effects of MBT1805, a novel balanced PPARα/γ/δ agonist, on cholestasis induced by α-naphthylisothiocyanate (ANIT) and elucidate the underlying mechanisms through untargeted and bile acid-targeted metabolomic analysis.

Methods: Levels of serum biochemical indicators (transaminase, aspartate transaminase, alkaline phosphatase, and total bilirubin) and liver histopathology were analyzed to evaluate the therapeutic effects of MBT1805 on ANIT-induced cholestasis in C57BL/6 mice. Untargeted and bile acid-targeted metabolomic analysis of liver tissues was performed using ultrahigh-performance liquid chromatography-triple quadrupole mass spectrometry (UPLC-MC/MC). qRT-PCR and Western blot analysis were carried out to measure the expression of key enzymes and transporters regulating bile acid synthesis, biotransformation, and transport.

Results: MBT1805 significantly improved abnormal levels of liver biochemical indicators and gallbladder enlargement induced by ANIT. Histopathological analysis showed that MBT1805 effectively relieved ANIT-induced necrosis, vacuolation, and inflammatory infiltration. Untargeted metabolomic analysis identified 27 metabolites that were involved in the primary biliary acid biosynthesis pathway. In addition, bile acid-targeted metabolomics showed that MBT1805 could alleviate the abnormal bile acid content and composition induced by ANIT. Furthermore, qRT-PCR and Western blot results confirmed that MBT1805 could effectively regulate bile acid synthesis, biotransformation, and transport which helps relieve cholestasis.

Conclusions: MBT1805 is a potential candidate drug for cholestasis, with a balanced PPARα/γ/δ activation effect.

Introduction

Impaired bile flow and accumulation of toxic bile acids are the main causes of cholestasis injury, which can lead to liver tissue destruction, inflammation infiltration, and subsequent fibrosis and cirrhosis (Arrese and Trauner, 2003). Patients who are not properly treated may develop end-stage liver disease. Once the end-stage liver disease develops, liver transplantation is the only treatment option since effective drugs are still lacking.

Primary biliary cholangitis (PBC) and primary sclerosing cholangitis (PSC) are the main chronic cholestasis disease. Ursodeoxycholic acid (UDCA) was the only treatment for PBC until the approval of obeticholic acid (OCA), which is approved as a monotherapy or combined with UDCA (Trivedi et al., 2016). Approximately 40% of patients with PBC have an inadequate response to UDCA treatment (Corpechot et al., 2011), whereas OCA has some adverse effects such as pruritus, fatigue, and discomfort (Kowdley et al., 2018). Although dose-dependent reductions in alkaline phosphatase (ALP) levels were observed in trials of PSC treatment with UDCA, long-term effects, such as cirrhosis manifestation, acute decompensation of cirrhosis, liver transplantation, and death events, need further assessment (Olsson et al., 2005).

Nuclear receptors (NRs) effectively regulate gene expression related to the elimination of toxic biliary constituents accumulating in cholestasis and are potential targets of PBC treatment (Halilbasic et al., 2013). Peroxisome proliferator-activated receptors (PPARs) belong to the NR family and comprise three isoforms: PPARα, PPARβ/δ, and PPARγ. PPARα agonists can treat hyperlipidemia, whereas PPARγ agonists effectively manage type 2 diabetes (Hong et al., 2018). However, some adverse effects have restricted the use of PPAR agonists, such as weight gain and fluid retention. Dual or pan-PPAR agonists are expected to reduce the side effects of single receptor activation. For instance, elafibranor, a dual PPARα/δ agonist, has completed phase-2 clinical trials and has reached the primary and secondary endpoints. Fibrates, such as fenofibrate and bezafibrate, have also been used to treat PBC with encouraging outcomes (Honda et al., 2013; Corpechot et al., 2018; Ghonem et al., 2020).

MBT1805 is a novel drug candidate built on resveratrol scaffolds and it is the most balanced PPAR-α/γ/δ agonist with moderate activity (EC50: 8.46, 11.94, and 11.5 μM, α/γ/δ, respectively). The synthesis of MBT1805 has been elaborately described in an early study (Li et al., 2010). Previous studies showed that activation of PPARα alleviates cholestasis injury caused by α-naphthylisothiocyanate (ANIT) (Fang et al., 2017; Dai et al., 2017; Zhao et al., 2017; Dai et al., 2018). The animal model of cholestasis induced by ANIT used in this research was widely used in previous studies, especially in the exploration of novel treatment drugs for cholestasis (Pollheimer and Fickert, 2015; Fang et al., 2017; Wang H. et al., 2017; Wu et al., 2017).

Metabolomics focuses on the study of endogenous small molecules with relative molecular weights of <1,000 Da (Nicholson et al., 1999), mainly using nuclear magnetic resonance (NMR) and mass spectrometry. Untargeted metabolomics intends to comprehensively analyze all measurable metabolites that are widely used in mechanism analysis and biomarker discovery, whereas targeted metabolomics focuses on the analysis of specific metabolic components that are suitable for pharmacokinetic studies (Roberts et al., 2012; Gertsman and Barshop, 2018).

In this study, we aimed to evaluate the role of MBT1805 for treating ANIT-induced cholestatic liver injury and revealed the possible therapeutic mechanisms through untargeted and bile acid-targeted metabolomics.

Materials and Methods

Animals and Treatments

Six-week-old male C57BL/6 mice were purchased from Charles River Laboratories (Beijing, China). Mice were housed in a breeding room at the Experimental Animal Center of the Translational Medical Research Institute, First Hospital of Jilin University, Changchun, China. They were maintained at 25°C ± 1 C, 60%–70% humidity, and a light/dark cycle of 12 h with libitum access to feed and water. After adaptive feeding of normal growth chow from Xietong (Jiangsu, China) for 1 wk, mice were divided into groups: control, no treatment (n = 10); model, treatment only with ANIT (n = 10); UDCA-60, treatment with 60 mg kg⁻¹ UDCA and ANIT (n = 10). The dose of UDCA was referred to the research of Wang et al. (Wang L. et al., 2017): MBT1805-10, treatment with 10 mg kg⁻¹ MBT1805 and ANIT (n = 10); MBT1805-20, treatment with 20 mg kg⁻¹ and ANIT (n = 10); and MBT1805-30, treatment with 30 mg kg⁻¹ MBT1805 and ANIT (n = 10). UDCA and MBT1805 were administrated by gavage for 7 days. All mice were euthanized 24 h after the administration of 75 mg kg⁻¹ ANIT. ANIT and MBT1805 were supplied by Beijing, JK HuaYuan Med Tech Company LTD (Beijing, China).

Biochemical Analysis

Blood samples were collected from the heart, placed into 4 ml centrifuge tubes, and left at 25°C for 30 min. Serum samples were centrifuged at 4,000 rpm for 10 min at 4°C. Serum alanine aminotransferase (ALT), aspartate aminotransferase (AST), ALP, and total bilirubin (Tbil) levels were measured using assay kits (Nanjing Jiancheng Bioengineering Institute, Nanjing, China), according to the manufacturer’s instructions. Triglyceride and total protein were detected using assay kits (Nanjing Jiancheng Bioengineering Institute, Nanjing, China).

Histopathological Evaluation

The liver tissue was cut into several parts, one of which was fixed with 4% paraformaldehyde solution. Fixed liver tissues of the above experiments were dehydrated in a serial concentration of alcohol and xylene followed by paraffin embedding. 4 μm of each sample was cut and stained with hematoxylin and eosin. Other liver tissues were stored in liquid nitrogen for subsequent RNA, protein extraction, and metabolomic analysis.

Metabolomics Study

Sample Preparation and Operation Conditions

Liver tissue samples stored in liquid nitrogen were thawed on ice. For untargeted metabolomics, 50 ± 2 mg of tissue was mixed with cold steel balls and homogenized at 30 Hz for 3 min. Therefore, 1 ml of 70% methanol with standard internal extract (Metware Biotechnology, Wuhan, China) was added to the homogenized tissue and mixed by stirring for 5 min. The supernatant was collected after centrifugation at 12,000 rpm for 10 min at 4°C and stored at −20°C overnight. After centrifugation at 12,000 rpm for 3 min at 4°C, 200 μl of the supernatant was collected for ultrahigh performance liquid chromatography-triple quadrupole mass spectrometry (UPLC-MC/MC).

For bile acid-targeted metabolomics, 20 mg of tissue was ground with a ball mill, mixed with 200 μl methanol, and then stored at −20 C to precipitate protein. Extracts were dried by evaporation and reconstituted in 100 μl of 50% methanol (v/v) for UPLC-MC/MC. The instrumentation and operation conditions are described in the appendix, section1 (Supplementary Appendix S1). To identify metabolites, we used the self-built widely targeted metabolome database MWDB (Metware Biotechnology, Wuhan, China).

Metabolomic Data Analysis

Principal component analysis (PCA) and partial least squares discrimination analysis (PLS-DA) of metabolomic data were performed using SIMCA-P 14.1 (Umetrics, Umea, Sweden). All variables were mean-centered and scaled to a Pareto variance before PCA. Variable influence on the projection (VIP) was used to identify variables responsible for group separation, and only those with VIP ≥2.0 were selected for further analysis. Compounds with significant changes among groups (p < 0.05) were considered biomarkers and used for enrichment and pathway analysis with MetaboAnalyst 5.0.

Quantitative Analysis of Hepatic Bile Acids

UPLC-MS/MS was used to detect the levels of various bile acids in mouse livers. The instrumentation and operation conditions are described in the appendix, section 2 (Supplementary Appendix S2).

qRT-PCR

All RNA from liver tissues was isolated using the Promega SV Total Isolation System (Promega, Madison, WI, United States) and quantified using a microplate reader. cDNA was synthesized using the PrimerScript first strand cDNA Synthesis Kit (Takara, Shiga, Japan). The mRNA expressions of key enzymes (CYP27A1, CYP7A1, and CYP8b1), transport (NTCP, bile salt excretion pump (BSEP), MRP2, MRP3, and MRP4), and metabolic enzymes (CYP) were detected using the Agilent Mx3005 P Real-Time PCR System (Agilent, Santa Clara, CA, United States) with the SYBR Green PCR kit. The thermal conditions were as follows: 95°C for 60 s, followed by 45 cycles at 95°C for 15 s and 60°C for 60 s. β-Actin was used as an internal control. All primers used for qRT-PCR were synthesized by Sangon Biotech (Shanghai, China) and are presented in appendix, section 3 (Supplementary Appendix S3).

Western Blot Analysis

Protein extracted from rat liver tissue was prepared using radio-immunoprecipitation assay lysis buffer (Beyotime, Shanghai, China). Protein concentrations were detected using a bicinchoninic acid protein assay (BCA) kit. Protein extract was loaded and electrophoresed through sodium dodecyl sulfate-polyacrylamide gel electrophoresis (SDS-PAGE) and transferred onto a polyacrylamide difluoride membrane (Immobilon®-P^SQ). After blocking with 5% skim milk in Tris-buffered saline containing 0.1% Tween-20 (TBST) for 1 h at room temperature, membranes were incubated overnight at 4 C with primary antibodies and washed with 0.1% TBST Tris-buffered saline. The membranes were subsequently incubated with HRP-conjugated secondary antibodies for 30 min and then washed with 0.1% TBST Tris-buffered saline. Then, enhanced chemiluminescent kit (ECL) (New Cell & Molecular Biotech Co., Ltd) was used to detect the protein-antibody complexes using ChemiDoc XRS+ system (Bio-Rad, United States). All primary and secondary antibodies used in this experiment were supplied in appendix, section 4 (Supplementary Appendix S4).

Statistical Analysis

Data were expressed as means ± standard error. Statistical differences among groups were determined by Welch and Brown-Forsythe one-way analysis of variance and the post hoc test of Dunnett T3. Statistical analysis was performed by SPSS 23.0. The statistical significance was set at p < 0.05.

Results

Biochemical Analysis and Pathological Assessment

Compared with the control, all serum indicators were markedly increased in the model group (Figures 1A–D). UDCA and MBT1805 were effective in reducing ALT, AST, and Tbil, and MBT1805 was better. In addition, ALT and AST levels significantly reduced in the MBT1805-treated group in a dose-dependent manner (Figures 1A,B). Reduction of ALP was only observed in high dose of MBT1805 compared with the model group. ALP level was lower in the UDCA-treated group than that of the MBT1805-treated group (Figure 1C). Compared with the control, the gallbladder of the model was abnormally enlarged and filled with bile acids, whereas that of the MBT1805-treated groups appeared normal without any liver tissue congestion (Figure 2A).

FIGURE 1

FIGURE 1. Serum parameters. (A) Alanine transaminase (ALT); (B) Aspartate aminotransferase (AST); (C) Alkaline phosphatase (ALP); (D) Total bilirubin (Tbil). #p < 0.05 and ##p < 0.01 compared with the control group; *p < 0.05 and **p < 0.01 compared with the ANIT group.

FIGURE 2

FIGURE 2. Appearance of gallbladders and pathology of liver tissues. (A) Appearance of gallbladders in different groups. (B) Pathology of liver tissues with magnification of 100×. (C) Pathology of liver tissues with magnification of 200×.

Histopathological evaluation revealed significant hepatic tissue changes, such as necrosis, inflammatory infiltration, and vacuolization, in the model. In contrast, the UDCA- and MBT1805-treated groups exhibited less inflammatory infiltration and no necrotic tissues (Figures 2B,C).

Multivariate Statistical Analysis

In total, 793 metabolites were recognized after the combination of ESI⁺/ESI⁻ detection signals. PCA identified differences among the control, model, and MBT1805-30 groups (Supplementary Appendix S5). We found that the sample named treat6 in the MBT1805-treated group was an outlier (confidence interval, >95%) in PCA and was, thus, removed from the following analysis. The clustering of the model was significantly different from that of the control or MBT1805-30; however, no differences were identified between the latter two groups (Figure 3A). Therefore, further analysis is needed. PLS-DA separated the control, model, and MBT1805-treated groups. R²Y and Q² were 0.906 and 0.832, respectively (Figure 3B), indicating the high accuracy of predictions. A permutation test was performed with 100 iterations. The validation plots had an intercept of Q² < 0.5 (Figure 3C), and VIP ≥2.0 was used as a threshold to identify potential metabolites (Figure 3D).

FIGURE 3

FIGURE 3. PCA and PLS-DA. (A) PCA score of control (n = 6), ANIT (n = 6), and MBT1805 30 mg/kg treatment (treat, n = 5). (B) Score plots of PLS-DA. (C) PLS-DA model maturation test with 100 (Q2 intercept =0.048). (D) VIP values with threshold, VIP ≥ 2.

Identification of Potential Metabolites

Among the 793 detected compounds, 52 were selected from the control, model, and MBT1805-30 as candidates for Welch and Brown-Forsythe one-way ANOVA. Of these, variables that significantly differed among the groups (p < 0.05) were considered candidate biomarkers. Table 1 summarizes the 27 compounds identified by their molecular weight, formula. p values of Welch and Brown-Forsythe one-way analysis of variance and the post hoc test of Dunnett T3 are also presented in Table 1.

TABLE 1

TABLE 1. Metabolites with significant differences among groups.

Pathway Analysis

The heat map revealed changes in various types of substances among the control, model, and MBT1805-treated groups (Figure 4A). Upregulation is in red and downregulation is in green. Control, ANIT, and MBT1805-treated groups were significantly separated due to expression difference. Bile acids in the ANIT group were upregulated while reduced in the MBT1805-treated group. Figure 4B shows the heatmap of bile acids in bile acid-targeted metabolomics. MBT1805 significantly reversed ANIT-induced abnormal contents of bile acids. As shown in Figure 4C, the top four pathways were pyrimidine metabolism, glycine, serine and threonine metabolism, glutathione metabolism, and primary bile acid biosynthesis. Untargeted and bile acid-targeted metabolomics indicated that MBT1805 was effective in improving cholestasis and probably achieved by regulating bile acid metabolism.

FIGURE 4

FIGURE 4. Heat map and pathways analysis among the control, ANIT, and MBT1805 treated groups. (A) Heat map of potential metabolites. Red in gradient presented the increases. Green in gradient presented the decreases. Samples with green were the control group (n = 6), samples with red were the ANIT group (n = 6), and samples with purple were the MBT1805-treated group. (B) Heat map of bile acids of bile acid-targeted metabolomic. (C) Bubble map of control, ANIT, and MBT1805 30 mg/kg groups in pathway analysis.

Quantitative Analysis of Bile Acids

Metabolomic analysis revealed that the bile acid biosynthesis pathway was an essential pathway altered by MBT1805. In agreement with the results of untargeted metabolomics, the levels of conjugated bile acids, taurine-conjugated bile acids (TCA, TCDCA, TUDCA, and Tβ-MCA), and glycine-conjugated bile acids (GCA and GUDCA) were significantly elevated in the model compared with the control but significantly reduced in MBT1805-treated groups (Figure 5). In summary, bile acid metabolism was disrupted by ANIT, whereas MBT1805 treatment significantly reversed the abnormality of bile acid homeostasis.

FIGURE 5

FIGURE 5. Bile acid levels of different groups. Analysis of 21 common bile acids (A). Analysis of several hydrophobic bile acids, Deoxycholic acid, DCA (B); Glycolithocholic acid-3-sulfate, GLCA-3S (C); Chenodeoxycholic acid-3-β-D-glucuronide, CDCA-3Gln (D). Control, control group; ANIT, model group; MBT1805, MBT1805 30 mg/kg treated group; #p < 0.05, ##p < 0.01, vs. control; *p < 0.05, vs. ANIT.

Expression Levels of Key Enzymes and Transporters That Regulate Bile Acid Synthesis, Detoxification, and Transport

Primary bile acid, including CA and CDCA, is synthesis by classical and alternative pathways, which are mediated by CYP7A1 and CYP27A1, respectively (Wahlström et al., 2016). In ANIT-induced cholestasis, the synthesis of bile acid increased, especially through an alternative pathway. However, mRNA and protein expressions were inconsistent with CYP 27A1 (Figures 6A,B). We supposed that a decrease in the mRNA level of CYP27A1 was an adaptive reaction when cholestasis had formed. MBT1805 significantly reduced the expression of CYP27A1; thus, the synthesis of bile acid was interrupted. In the cholestasis model, transport of bile acid from hepatocytes to bile duct and to systemic circulation was increased, because the expressions of BSEP and multidrug resistance proteins (MRP2, MRP3, and MRP4) were evaluated. BESP and MRP2 mediate transport of bile acid from hepatocytes to the bile duct, while MRP3 and MRP4 mediate transport to the systemic circulation (Jia et al., 2018). MBT1805 decreased bile acid transport to circulation; therefore, serum indicators significantly improved. Bile acid reabsorption mediated by sodium taurocholate cotransporting polypeptide (NTCP) was inhibited in the cholestasis model and was restored in the MBT1805-treated group. Hepatic biotransformation, especially phase Ⅱ metabolic reaction, is crucial in reducing the toxicity of hydrophobic bile acid. MBT1805 could also regulate the expressions of UDP-glucuronosyltransferases, which contributed to detoxication.

FIGURE 6

FIGURE 6. Influence of MBT1805 on the expression of bile acid synthesis enzymes and transports. Data are mean ± SD, n = 5. #p < 0.05 and ##p < 0.01 compared with control group; *p < 0.05 and **p < 0.01 compared with ANIT group. (A,C) Western blot; (B,D) qRT-PCR.

Discussion

Accumulation of toxic bile acids due to destruction or obstruction of bile ducts is the main cause of cholestasis. UDCA and OCA are currently recommended treatment options for cholestatic liver injury (Gao et al., 2020). However, up to 40% of patients with PBC have an inadequate response to UDCA treatment, whereas OCA can aggravate pruritus, which is a common symptom of cholestasis (Roberts et al., 2020; Gomez et al., 2021). PPAR agonists, used as drugs for hyperlipidemia and type2 diabetes, show some encouraging results in PBC treatment (Jones et al., 2017; Corpechot et al., 2018) and, thus, are considered potential therapeutic targets.

MBT1805 was built on resveratrol scaffolds and has a balanced PPARα/γ/δ activation effect. In this study, we aimed to assess the therapeutic effects of MBT1805 and explore the underlying mechanisms of action against ANIT-induced cholestasis. ANIT causes pathological necrosis of hepatocytes and intrahepatic bile ducts, periportal inflammation, and fibrosis in mice (Mariotti et al., 2018). We found that serum levels of biochemical indicators of liver injury (ALT and AST) significantly increased in the model compared with the control and that MBT1805 could significantly reduce the levels in a dose-dependent manner (Figures 1A,B). A similar pattern was observed for Tbil (Figure 1C). ALT, AST, and Tbil levels in the MBT1805-treated group were lower than those in the UDCA-treated group. But ALP level (Figure 1C) was lower in the UDCA-treated group. Additionally, reduced ALP was observed only in the high dose of the MBT1805-treated group. We speculated that the difference of PPARs distribution and affinity to agonist between hepatocytes and bile duct cells might provide an explanation. PPARα (Yang et al., 2019) and PPARγ (Xiang et al., 2021) have been reported to have protective effects in cholestasis by regulating bile acid metabolism in hepatocytes. In cholangiocytes, studies mainly focused on PPARγ, which played an important role in maintaining immune tolerance (Marra et al., 2005; Zhang et al., 2018). Therefore, hepatocytes and cholangiocytes could respond inconsistently to pan-PPAR agonist. Histopathological observations revealed that bile acid was highly accumulated in the gallbladder of the model, whereas the levels in the MBT1805-treated groups were similar to those in the control (Figure 2A). Besides, necrosis of the liver tissues, vacuolation, and inflammatory infiltration were observed in the model; however, no necrosis was found in the MBT1805-treated groups, whereas inflammatory infiltration was markedly reduced. Overall, the serological and histopathological results suggested that MBT1805 was an effective anticholestasis compound.

To investigate the related mechanisms of MBT1805 against cholestasis, we employed metabolomic analysis using UPLC-MS/MS, which mainly focuses on substances with molecular weights <1,000 Da. We detected 793 compounds after combining ESI⁺/ESI⁻ detection patterns and deleting any duplicates. PCA distinguished the model from the control and MBT1805-treated groups and revealed an overlap of the latter two groups. Metabolites with VIP ≥2 were selected after PLS-DA of the control, model, and MBT1805-treated groups. In total, 27 compounds were identified as significant markers (Table 1). In the model, the bile acid metabolism increased, and the phospholipid metabolism decreased. Fang et al. reported that elevated levels of lysophosphatidylcholine 18:0 activate the NF-κB/IL-6/STAT3 signaling pathway in ANIT-induced hepatotoxicity and aggravate inflammatory damage. In ANIT-induced injury, the ANIT-glutathione complex is secreted into the bile mediated by MRP2, exhausting glutathione that plays an essential role in eliminating free radicals and protecting against oxidative stress injury (Christoph G. Dietrich, 2001). We found that MBT1805 significantly reduced the levels of lysophosphatidylcholine 18:0 and oxidative glutathione and increased the level of glutathione, which could contribute to decreased hepatic inflammatory infiltration and reactive oxygen species production.

Pathway analysis identified that the top four pathways were pyrimidine metabolism, glycine, serine and threonine metabolism, glutathione metabolism, and primary bile acid biosynthesis, indicating that MBT1805 probably influences bile acid metabolism and plays an important role in cholestasis. To further explore the effects of MBT1805 on bile acid metabolism in cholestasis, we carried out bile acid target metabolomics and identified more than 20 bile acids. Cholestasis is characterized by the intrahepatic accumulation of excessive bile acids, which probably induces hepatic parenchymal cell death, bile duct proliferation, liver inflammation, and fibrosis (Cai and Boyer, 2021). Hydrophobic bile acid retention promotes oxidative stress that causes mitochondrial dysfunction, leading to liver damage (Arduini et al., 2012). In the present study, metabolomics revealed that levels of bile acids, including CA, GCA, β/ωMCA, and Tβ-MCA, were modulated by ANIT and reversed by MBT1805. In particular, some hydrophobic and toxic acid derivatives such as DCA, GLCA-3S, and CDCA-3Gln were reduced to undetectable levels in the MBT1805-treated group, as shown in Figures 5B–D. Therefore, MBT1805 effectively reduced the accumulation of toxic bile acids accumulation and their damage to hepatocytes.

Hepatic bile acid metabolism is a complex biochemical process. Hepatic cholestasis occurs when bile acid metabolism processes are abnormal, including bile acid synthesis, biotransformation, and transport. qRCR and Western blot were performed to verify the expression of key enzymes and transporters regulating bile acid synthesis, biotransformation, and transport. Herein, we found an increase of bile acid synthesis in ANIT-induced cholestasis, particularly alternative pathway of bile acid synthesis. CYP27A1 is the rate-limiting enzyme of alternative pathway (Wahlström et al., 2016), whose expression significantly increases in the model group (Figure 6A). However, results of mRNA and protein expression levels were inconsistent (Figures 6A,B), and we speculated that a decrease of mRNA expression could be an adaptive response when cholestasis had been already formed. BESP and MRP2 mediate bile acid transport from hepatocytes to bile duct, while MRP3 and MRP4 mediate bile acid transport from hepatocytes to the systemic circulation (Jia et al., 2018). In the ANIT-induced cholestasis model, transport from hepatocytes to bile duct and to systemic circulation was enhanced, and reabsorption of bile acid mediated by NTCP was inhibited, while MBT1805 significantly reversed the abnormal expression of bile acid transport. Glucuronosyltransferase enzymes play an important role in hepatic phase II metabolic reactions, which mediate the conversion of hydrophobic bile acids to hydrophilic bile acid, thus reducing bile acid toxicity (Burchell, 2003). MBT1805 could enhance bile acid biotransformation (Figures 6C,D), which contributed to relieve the toxicity of bile acid.

Hepatic steatosis accompanies cholestasis usually. We have evaluated hepatic triacylglycerol to describe hepatic cytoarchitecture. As shown in appendix (Supplementary Figure S6A), triacylglycerol level in the ANIT-treated group elevated, without significant difference. In addition, the decrease of triacylglycerol in MBT1805-treated also had no statistical difference. We speculated that ANIT-induced cholestasis was a model of acute liver injury. Steatosis could be present if damage persisted. In addition, mRNA expressions of PPARα, PPARβ/δ, and PPARγ were detected. As we see in Supplementary Figure S6B, mRNA expressions of PPARα, PPARβ/δ, and PPARγ were significantly elevated.

In summary, MBT1805 significantly improved cholestasis induced by ANIT. Hepatic nontargeted and bile acid-targeted metabolomic analysis suggested that MBT1805 could alleviate cholestasis through regulating bile acid metabolic pathway. Results of qRT-PCR and Western blot indicated that MBT1805 significantly inhibited synthesis of bile acid and enhanced phase II metabolic reactions.

Conclusion

Our study focused on the anticholestatic effects of the novel PPARα/γ/δ agonist MBT1805. Serum biochemistry and liver histopathology indicated that MBT1805 effectively protected against cholestasis. Untargeted metabolomics combined with bile acid-targeted metabolomics revealed that the abnormal primary bile acid biosynthesis caused by ANIT was recovered by MBT1805. We also verified the key enzymes and transporters involved in the bile acid metabolism pathway, showing that MBT1805 significantly reversed the abnormal synthesis, metabolism, and transport of bile acids. In summary, MBT1805 is a potential candidate drug for cholestasis.

Data Availability Statement

The original contributions presented in the study are included in the article/Supplementary Material, further inquiries can be directed to the corresponding authors.

Ethics Statement

The animal study was reviewed and approved by the Animal Ethics Committee of the First Hospital of Jilin University.

Author Contributions

CW mainly completed the experiment and wrote this article with the help of YS and XW. MBT1805 used in this research was supplied by BZ, and FP completed the experiment about proper doses in the early stage. XJ and JN supervised the overall design of the experiment and provided funding and equipment.

Funding

This study was supported by the National Natural Science Foundation of China (grants No. 81970519), Program for JLU Science and Technology, Innovative Research Team (2017TD-08), and the Fundamental Research Funds for the Central Universities.

Conflict of Interest

Author BZ was employed by Beijing JK HuaYuan Med Tech Company LTD.

The remaining authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Publisher’s Note

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations or those of the publisher, the editors, and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.

Acknowledgments

The authors wish to express their gratitude to Metware Biotechnology and Editage for their contributions to metabolites detection and manuscript polish. Hongqin Xu from the First Hospital of Jilin University was especially appreciated for her guidance in statistical analysis.

Supplementary Material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fphar.2021.732478/full#supplementary-material

References

Arduini, A., Serviddio, G., Tormos, A. M., Monsalve, M., and Sastre, J. (2012). Mitochondrial Dysfunction in Cholestatic Liver Diseases. Front. Biosci. (Elite Ed.) 4, 2233–2252. doi:10.2741/539

PubMed Abstract | CrossRef Full Text | Google Scholar

Arrese, M., and Trauner, M. (2003). Molecular Aspects of Bile Formation and Cholestasis. Trends Mol. Med. 9 (12), 558–564. doi:10.1016/j.molmed.2003.10.002

PubMed Abstract | CrossRef Full Text | Google Scholar

Burchell, B. (2003). Genetic Variation of Human UDP-Glucuronosyltransferase: Implications in Disease and Drug Glucuronidation. Am. J. Pharm. 3 (1), 37–52. doi:10.2165/00129785-200303010-00006

CrossRef Full Text | Google Scholar

Cai, S. Y., and Boyer, J. L. (2021). The Role of Bile Acids in Cholestatic Liver Injury. Ann. Transl Med. 9 (8), 737. doi:10.21037/atm-20-5110

PubMed Abstract | CrossRef Full Text | Google Scholar

Corpechot, C., Chazouillères, O., and Poupon, R. (2011). Early Primary Biliary Cirrhosis: Biochemical Response to Treatment and Prediction of Long-Term Outcome. J. Hepatol. 55 (6), 1361–1367. doi:10.1016/j.jhep.2011.02.031

CrossRef Full Text | Google Scholar

Corpechot, C., Chazouillères, O., Rousseau, A., Le Gruyer, A., Habersetzer, F., Mathurin, P., et al. (2018). A Placebo-Controlled Trial of Bezafibrate in Primary Biliary Cholangitis. N. Engl. J. Med. 378 (23), 2171–2181. doi:10.1056/NEJMoa1714519

PubMed Abstract | CrossRef Full Text | Google Scholar

Dai, M., Hua, H., Lin, H., Xu, G., Hu, X., Li, F., et al. (2018). Targeted Metabolomics Reveals a Protective Role for Basal PPARα in Cholestasis Induced by α-Naphthylisothiocyanate. J. Proteome Res. 17 (4), 1500–1508. doi:10.1021/acs.jproteome.7b00838

PubMed Abstract | CrossRef Full Text | Google Scholar

Dai, M., Yang, J., Xie, M., Lin, J., Luo, M., Hua, H., et al. (2017). Inhibition of JNK Signalling Mediates PPARα-dependent Protection Against Intrahepatic Cholestasis by Fenofibrate. Br. J. Pharmacol. 174 (18), 3000–3017. doi:10.1111/bph.13928

CrossRef Full Text | Google Scholar

Dietrich, C. G., Ottenhoff, R., de Waart, D. R., and Oude Elferink, R. P. (2001). Role of MRP2 and GSH in Intrahepatic Cycling of Toxins. Toxicology 167, 73–81. doi:10.1016/s0300-483x(01)00459-0

PubMed Abstract | CrossRef Full Text | Google Scholar

Fang, Z. Z., Tanaka, N., Lu, D., Jiang, C. T., Zhang, W. H., Zhang, C., et al. (2017). Role of the Lipid-Regulated NF-κB/IL-6/STAT3 Axis in Alpha-Naphthyl Isothiocyanate-Induced Liver Injury. Arch. Toxicol. 91 (5), 2235–2244. doi:10.1007/s00204-016-1877-6

PubMed Abstract | CrossRef Full Text | Google Scholar

Gao, L., Wang, L., Woo, E., He, X., Yang, G., Bowlus, C., et al. (2020). Clinical Management of Primary Biliary Cholangitis-Strategies and Evolving Trends. Clin. Rev. Allergy Immunol. 59 (2), 175–194. doi:10.1007/s12016-019-08772-7

PubMed Abstract | CrossRef Full Text | Google Scholar

Gertsman, I., and Barshop, B. A. (2018). Promises and Pitfalls of Untargeted Metabolomics. J. Inherit. Metab. Dis. 41 (3), 355–366. doi:10.1007/s10545-017-0130-7

CrossRef Full Text | Google Scholar

Ghonem, N. S., Auclair, A. M., Hemme, C. L., Gallucci, G. M., de la Rosa Rodriguez, R., Boyer, J. L., et al. (2020). Fenofibrate Improves Liver Function and Reduces the Toxicity of the Bile Acid Pool in Patients with Primary Biliary Cholangitis and Primary Sclerosing Cholangitis Who Are Partial Responders to Ursodiol. Clin. Pharmacol. Ther. 108 (6), 1213–1223. doi:10.1002/cpt.1930

PubMed Abstract | CrossRef Full Text | Google Scholar

Gomez, E., Garcia Buey, L., Molina, E., Casado, M., Conde, I., Berenguer, M., et al. (2021). Effectiveness and Safety of Obeticholic Acid in a Southern European Multicentre Cohort of Patients with Primary Biliary Cholangitis and Suboptimal Response to Ursodeoxycholic Acid. Aliment. Pharmacol. Ther. 53 (4), 519–530. doi:10.1111/apt.16181

PubMed Abstract | CrossRef Full Text | Google Scholar

Halilbasic, E., Baghdasaryan, A., and Trauner, M. (2013). Nuclear Receptors as Drug Targets in Cholestatic Liver Diseases. Clin. Liver Dis. 17 (2), 161–189. doi:10.1016/j.cld.2012.12.001

PubMed Abstract | CrossRef Full Text | Google Scholar

Honda, A., Ikegami, T., Nakamuta, M., Miyazaki, T., Iwamoto, J., Hirayama, T., et al. (2013). Anticholestatic Effects of Bezafibrate in Patients with Primary Biliary Cirrhosis Treated with Ursodeoxycholic Acid. Hepatology 57 (5), 1931–1941. doi:10.1002/hep.26018

PubMed Abstract | CrossRef Full Text | Google Scholar

Hong, F., Xu, P., and Zhai, Y. (2018). The Opportunities and Challenges of Peroxisome Proliferator-Activated Receptors Ligands in Clinical Drug Discovery and Development. Int. J. Mol. Sci. 19 (8), 2189. doi:10.3390/ijms19082189

PubMed Abstract | CrossRef Full Text | Google Scholar

Jia, W., Xie, G., and Jia, W. (2018). Bile Acid-Microbiota Crosstalk in Gastrointestinal Inflammation and Carcinogenesis. Nat. Rev. Gastroenterol. Hepatol. 15 (2), 111–128. doi:10.1038/nrgastro.2017.119

PubMed Abstract | CrossRef Full Text | Google Scholar

Jones, D., Boudes, P. F., Swain, M. G., Bowlus, C. L., Galambos, M. R., Bacon, B. R., et al. (2017). Seladelpar (MBX-8025), A Selective PPAR-δ Agonist, in Patients with Primary Biliary Cholangitis with an Inadequate Response to Ursodeoxycholic Acid: a Double-Blind, Randomised, Placebo-Controlled, Phase 2, Proof-Of-Concept Study. Lancet Gastroenterol. Hepatol. 2 (10), 716–726. doi:10.1016/s2468-1253(17)30246-7

PubMed Abstract | CrossRef Full Text | Google Scholar

Kowdley, K. V., Luketic, V., Chapman, R., Hirschfield, G. M., Poupon, R., Schramm, C., et al. (2018). A Randomized Trial of Obeticholic Acid Monotherapy in Patients with Primary Biliary Cholangitis. Hepatology 67 (5), 1890–1902. doi:10.1002/hep.29569

PubMed Abstract | CrossRef Full Text | Google Scholar

Li, W., He, X., Shi, W., Jia, H., and Zhong, B. (2010). Pan-PPAR Agonists Based on the Resveratrol Scaffold: Biological Evaluation and Docking Studies. ChemMedChem 5 (12), 1977–1982. doi:10.1002/cmdc.201000360

PubMed Abstract | CrossRef Full Text | Google Scholar

Mariotti, V., Strazzabosco, M., Fabris, L., and Calvisi, D. F. (2018). Animal Models of Biliary Injury and Altered Bile Acid Metabolism. Biochim. Biophys. Acta Mol. Basis Dis. 1864 (4 Pt B), 1254–1261. doi:10.1016/j.bbadis.2017.06.027

PubMed Abstract | CrossRef Full Text | Google Scholar

Marra, F., DeFranco, R., Robino, G., Novo, E., Efsen, E., Pastacaldi, S., et al. (2005). Thiazolidinedione Treatment Inhibits Bile Duct Proliferation and Fibrosis in a Rat Model of Chronic Cholestasis. World J. Gastroenterol. 11 (32), 4931–4938. doi:10.3748/wjg.v11.i32.4931

CrossRef Full Text | Google Scholar

Nicholson, J. K., Lindon, J. C., and Holmes, E. (1999). 'Metabonomics': Understanding the Metabolic Responses of Living Systems to Pathophysiological Stimuli via Multivariate Statistical Analysis of Biological NMR Spectroscopic Data. Xenobiotica 29 (11), 1181–1189. doi:10.1080/004982599238047

PubMed Abstract | CrossRef Full Text | Google Scholar

Olsson, R., Boberg, K. M., de Muckadell, O. S., Lindgren, S., Hultcrantz, R., Folvik, G., et al. (2005). High-dose Ursodeoxycholic Acid in Primary Sclerosing Cholangitis: A 5-year Multicenter, Randomized, Controlled Study. Gastroenterology 129 (5), 1464–1472. doi:10.1053/j.gastro.2005.08.017

PubMed Abstract | CrossRef Full Text | Google Scholar

Pollheimer, M. J., and Fickert, P. (2015). Animal Models in Primary Biliary Cirrhosis and Primary Sclerosing Cholangitis. Clin. Rev. Allergy Immunol. 48 (2-3), 207–217. doi:10.1007/s12016-014-8442-y

PubMed Abstract | CrossRef Full Text | Google Scholar

Roberts, L. D., Souza, A. L., Gerszten, R. E., and Clish, C. B. (2012). Targeted Metabolomics. Curr. Protoc. Mol. Biol. 98 (1), 1–34. doi:10.1002/0471142727.mb3002s98

CrossRef Full Text | Google Scholar

Roberts, S. B., Ismail, M., Kanagalingam, G., Mason, A. L., Swain, M. G., Vincent, C., et al. (2020). Real-World Effectiveness of Obeticholic Acid in Patients with Primary Biliary Cholangitis. Hepatol. Commun. 4 (9), 1332–1345. doi:10.1002/hep4.1518

PubMed Abstract | CrossRef Full Text | Google Scholar

Trivedi, P. J., Hirschfield, G. M., and Gershwin, M. E. (2016). Obeticholic Acid for the Treatment of Primary Biliary Cirrhosis. Expert Rev. Clin. Pharmacol. 9 (1), 13–26. doi:10.1586/17512433.2015.1092381

PubMed Abstract | CrossRef Full Text | Google Scholar

Wahlström, A., Sayin, S. I., Marschall, H. U., and Bäckhed, F. (2016). Intestinal Crosstalk Between Bile Acids and Microbiota and its Impact on Host Metabolism. Cell Metab 24 (1), 41–50. doi:10.1016/j.cmet.2016.05.005

PubMed Abstract | CrossRef Full Text | Google Scholar

Wang, H., Fang, Z. Z., Meng, R., Cao, Y. F., Tanaka, N., Krausz, K. W., et al. (2017a). Glycyrrhizin and Glycyrrhetinic Acid Inhibits Alpha-Naphthyl Isothiocyanate-Induced Liver Injury and Bile Acid Cycle Disruption. Toxicology 386, 133–142. doi:10.1016/j.tox.2017.05.012

PubMed Abstract | CrossRef Full Text | Google Scholar

Wang, L., Wu, G., Wu, F., Jiang, N., and Lin, Y. (2017b). Geniposide Attenuates ANIT-Induced Cholestasis through Regulation of Transporters and Enzymes Involved in Bile Acids Homeostasis in Rats. J. Ethnopharmacol 196, 178–185. doi:10.1016/j.jep.2016.12.022

CrossRef Full Text | Google Scholar

Wu, J. S., Li, Y. F., Li, Y. Y., Dai, Y., Li, W. K., Zheng, M., et al. (2017). Huangqi Decoction Alleviates Alpha-Naphthylisothiocyanate Induced Intrahepatic Cholestasis by Reversing Disordered Bile Acid and Glutathione Homeostasis in Mice. Front. Pharmacol. 8, 938. doi:10.3389/fphar.2017.00938

PubMed Abstract | CrossRef Full Text | Google Scholar

Xiang, J., Yang, G., Ma, C., Wei, L., Wu, H., Zhang, W., et al. (2021). Tectorigenin Alleviates Intrahepatic Cholestasis by Inhibiting Hepatic Inflammation and Bile Accumulation via Activation of PPARγ. Br. J. Pharmacol. 178 (12), 2443–2460. doi:10.1111/bph.15429

CrossRef Full Text | Google Scholar

Yang, S., Wei, L., Xia, R., Liu, L., Chen, Y., Zhang, W., et al. (2019). Formononetin Ameliorates Cholestasis by Regulating Hepatic SIRT1 and PPARα. Biochem. Biophys. Res. Commun. 512 (4), 770–778. doi:10.1016/j.bbrc.2019.03.131

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhang, H., Leung, P. S. C., Gershwin, M. E., and Ma, X. (2018). How the Biliary Tree Maintains Immune Tolerance. Biochim. Biophys. Acta Mol. Basis Dis. 1864 (4 Pt B), 1367–1373. doi:10.1016/j.bbadis.2017.08.019

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhao, Q., Yang, R., Wang, J., Hu, D. D., and Li, F. (2017). PPARα Activation Protects Against Cholestatic Liver Injury. Sci. Rep. 7 (1), 9967. doi:10.1038/s41598-017-10524-6

PubMed Abstract | CrossRef Full Text | Google Scholar