
==== Front
Gut Microbes
Gut Microbes
Gut Microbes
1949-0976
1949-0984
Taylor & Francis

39264803
10.1080/19490976.2024.2387402
2387402
Version of Record
Research Article
Research Paper
Gut microbiota-mediated C-sulfonate metabolism impairs the bioavailability and anti-cholestatic efficacy of andrographolide
D. TANG ET AL.
GUT MICROBES
Tang Dafu a *
Hu Wanyu a *
Fu Bingxuan a
Zhao Xiaojie a
You Guoquan a
Xie Cong b
Wang Hong Yu a
Guo Xueni a
Zhang Qianbing c
Liu Zhongqiu a d
https://orcid.org/0000-0002-2012-7266
Ye Ling a
a NMPA Key Laboratory for Research and Evaluation of Drug Metabolism & Guangdong Provincial Key Laboratory of New Drug Screening, School of Pharmaceutical Sciences, Southern Medical University , Guangzhou, China
b Clinical Pharmacy Center, Nanfang Hospital, Southern Medical University , Guangzhou, China
c Cancer Research Institute, School of Basic Medical Sciences, Southern Medical University , Guangzhou, China
d International Institute for Translational Chinese Medicine, Guangzhou University of Chinese Medicine , Guangzhou, Guangdong, China
CONTACT Zhongqiu Liu liuzq@gzucm.edu.cn International Institute for Translational Chinese Medicine, Guangzhou University of Chinese Medicine, No. 232, Outer Ring East Road, Guangzhou, Guangdong 510006, China
Ling Ye yeling623@hotmail.com School of Pharmaceutical Sciences, Southern Medical University, No.1023, South Shatai Road, Baiyun District, Guangzhou, Guangdong 510515 China
* These authors contribute equally to this work.

12 9 2024
2024
12 9 2024
16 1 2387402Integra17 8 2024
Integra17 8 2024
15 10 2023
14 5 2024
29 7 2024
© 2024 The Author(s). Published with license by Taylor & Francis Group, LLC.
2024
The Author(s)
https://creativecommons.org/licenses/by/4.0/ This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited. The terms on which this article has been published allow the posting of the Accepted Manuscript in a repository by the author(s) or with their consent.

ABSTRACT

Cholestatic liver injury results from the accumulation of toxic bile acids in the liver, presenting a therapeutic challenge with no effective treatment available to date. Andrographolide (AP) has exhibited potential as a treatment for cholestatic liver disease. However, its limited oral bioavailability poses a significant obstacle to harnessing its potent therapeutic properties and restricts its clinical utility. This limitation is potentially attributed to the involvement of gut microbiota in AP metabolism. In our study, employing pseudo-germ-free, germ-free and strain colonization animal models, along with 16S rRNA and shotgun metagenomic sequencing analysis, we elucidate the pivotal role played by gut microbiota in the C-sulfonate metabolism of AP, a process profoundly affecting its bioavailability and anti-cholestatic efficacy. Subsequent investigations pinpoint a specific enzyme, adenosine-5’-phosphosulfate (APS) reductase, predominantly produced by Desulfovibrio piger, which catalyzes the reduction of SO42- to HSO3−. HSO3− subsequently interacts with AP, targeting its C=C unsaturated double bond, resulting in the formation of the C-sulfonate metabolite, 14-deoxy-12(R)-sulfo andrographolide (APM). Inhibition of APS reductase leads to a notable enhancement in AP bioavailability and anti-cholestatic efficacy. Furthermore, employing RNA sequencing analysis and farnesoid X receptor (FXR) knockout mice, our findings suggest that AP may exert its anti-cholestatic effects by activating the FXR pathway to promote bile acid efflux. In summary, our study unveils the significant involvement of gut microbiota in the C-sulfonate metabolism of AP and highlights the potential benefits of inhibiting APS reductase to enhance its therapeutic effects. These discoveries provide valuable insights into enhancing the clinical applicability of AP as a promising treatment for cholestatic liver injury.

GRAPHICAL ABSTRACT

KEYWORDS

Andrographolide
C-sulfonate metabolism
gut microbiota
oral bioavailability
anti-cholestatic effect
APS reductase
Desulfovibrio piger
FXR
National Natural Science Foundation of China 10.13039/501100001809 82274193 and 82074110 Guangdong Basic and Applied Basic Research Foundation 2024B1515020093 Science and Technology Innovation Project of Guangdong Medical Products Administration S2021ZDZ042 This work was supported by the National Natural Science Foundation of China [82274193 and 82074110], Guangdong Basic and Applied Basic Research Foundation [2024B1515020093] and the Science and Technology Innovation Project of Guangdong Medical Products Administration [S2021ZDZ042].
==== Body
pmc1. Introduction

Cholestasis represents a pivotal characteristic of cholestatic liver disease, and its untreated progression can lead to severe liver damage and end-stage liver diseases, thus constituting a significant public health concern.1,2 The report indicates that cholestasis occurs in approximately 35% of patients with chronic liver disease, with a higher prevalence in individuals diagnosed with primary biliary cholangitis and primary sclerosing cholangitis.3 Currently, therapeutic options for cholestatic liver disease are limited, with only ursodeoxycholic acid (UDCA) and obeticholic acid (OCA) approved by the Food and Drug Administration (FDA).4 Nevertheless, not all patients respond favorably to these treatments. UDCA is the first-line therapy in primary biliary cholangitis, the side-effects of UDCA treatment in most patients include bloating, weight gain, and sometimes self-reported thinning of hair.5,6 While OCA is currently approved for second-line therapy in patients with primary biliary cholangitis who have incomplete responses to UDCA or who are intolerant of UDCA, treatment with OCA is associated with a dose-dependent increase in pruritus, reported in roughly 40% of treated patients, which is the leading cause for drug discontinuation.7,8 S-adenosyl methionine (S-AMe) serves as a second-line treatment for cholestatic disease, primarily effective in the earlier stages of the disease.9 Hence, the development of novel therapeutic agents for cholestatic liver disease is of paramount importance. Natural products have demonstrated their potential as valuable sources of anti-cholestasis agents, offering a potential shortcut in the quest for effective treatments.10

Andrographolide (AP), the primary active compound found in Andrographis paniculata (Burm. f.) Nees, has been documented to exhibit protective and therapeutic effects against liver diseases,11 attributed to its well-known anti-inflammatory,12 anti-oxidative,13 and anti-fibrosis properties.14 Recent studies have revealed significant protective effects of AP against α-naphthylisothiocyanate (ANIT)-induced cholestatic liver injury in rats.15,16 Additionally, AP serves as a quality control component for Xiaoyan Lidan Formula,17 a classic Traditional Chinese Medicine formula known for its efficacy in alleviating various hepatobiliary disorders, including cholestasis, jaundice, and hepatitis in clinical practice.18 Except for Andrographis paniculata (Burm. f.) Nees (known as Chuanxinlian), this formula comprises two other Chinese herbs: Herba Rabdosia (the dried herb of Rabdosia serra (Maxim.) Hara, known as Xihuangcao), and Ramulus et Folium Picrasmae (the dried twigs and leaves of Picrasma quassioides (D. Don) Benn, known as Kumu).19 However, the clinical application of AP is severely hampered by its extremely low bioavailability, with rat absolute oral bioavailability ranging from 1.19% to 2.67%.20–22 Our previous study has illuminated that AP undergoes C-sulfonate metabolism in the rat intestine, resulting in the formation of a significant quantity of C-sulfonate metabolite 14-deoxy-12(R)-sulfo andrographolide (APM), which predominantly accounts for its poor oral bioavailability.22 Currently, the formulation of AP into solid lipid nanoparticles, a dosage form modification, offers a partial improvement in bioavailability, but AP still undergoes substantial metabolism into APM within the intestinal tract.23 Consequently, addressing the inhibition of metabolism and enhancement of AP’s bioavailability has emerged as an urgent challenge.

C-sulfonate metabolism denotes a process wherein sulfur atoms with lone pairs attack the β-carbon atoms of α, β-unsaturated carbonyl compounds to yield sulfonate compounds, although its precise mechanism remains elusive.22,24 This type of metabolism extends beyond AP and has been observed in various other compounds, including physalin A,25 nitecapone,26 2-pyridinethiol-1-oxide,27 and flumipropyn.28 Sulfonate metabolism can be mediated by host-derived sulfotransferases (SULTs) for O- and N-sulfonates,29 but likely not for C-sulfonates.30 It is highly probable that the gut microbiota contributes sulfite or bisulfite to mediate C-sulfonate metabolism, resulting in the generation of C-sulfonate metabolites.31

The gut microbiota serves as a pivotal “metabolic organ”, playing a crucial role in drug metabolism by harboring a diverse range of enzymes and metabolic pathways capable of both activating and inactivating drugs, which has significant implications for drug efficacy and safety.32,33 For instance, certain drugs undergo conversion into their active forms through gut microbiota-mediated metabolism, enhancing their pharmacological effects. For example, the gut microbiota facilitates the conversion of the prodrug sulfasalazine into its active anti-inflammatory component, 5-aminosalicylate,34 and sulindac into its active form, sulindac sulfide.35 Conversely, the gut microbiota can also metabolize drugs into inactive forms, thereby reducing their pharmacological activity. Specific gut bacteria, such as Lysinibacillus sphaericus, possess the ability to metabolize aspirin, decreasing its bioavailability.36 Additionally, microbiome-derived enzymes, such as gut commensal acetyltransferases, inactivate 5-aminosalicylic acid, consequently reducing its therapeutic efficacy.37 Moreover, the gut microbiota can diminish the effectiveness of the antidiabetic drug acarbose through microbial enzymes.38 This intricate interplay underscores the importance of considering microbiota-mediated metabolism in understanding drug responses and optimizing therapeutic outcomes. Consequently, further investigation is imperative to determine whether the gut microbiota is involved in the C-sulfonate metabolism of AP and its potential impact on AP’s bioavailability and efficacy.

The primary objective of our study was to examine the impact of gut microbiota-mediated C-sulfonate metabolism of AP on both its bioavailability and anti-cholestatic efficacy. Employing pseudo-germ-free, germ-free and strain colonization animal models, along with 16S rRNA and shotgun metagenomic sequencing analysis, we uncovered a significant revelation that the gut microbiota plays a pivotal role in the C-sulfonate metabolism of AP, subsequently influencing its bioavailability and anti-cholestatic effect. Furthermore, we identified adenosine-5’-phosphosulfate (APS) reductase as an enzyme primarily secreted by Desulfovibrio spp., particularly Desulfovibrio piger (D. piger), providing HSO3−, actively participating in the C-sulfonate metabolism of AP. Inhibiting APS reductase resulted in enhanced bioavailability and anti-cholestatic effect of AP. Moreover, utilizing RNA sequencing analysis and farnesoid X receptor (FXR) knockout mice, we delved into elucidating the potential mechanisms underlying AP’s anti-cholestatic effect through the activation of FXR pathway to promote bile acid efflux. In summary, our study elucidates the pivotal role of gut microbiota in the C-sulfonate metabolism of AP, emphasizing the potential advantages of APS reductase inhibitor to enhance its therapeutic effectiveness.

2. Materials and methods

2.1. Chemicals and reagents

Andrographolide (AP, purity > 98%) and genistein (purity >98%, as an internal standard (IS)) were obtained from Chengdu Pufeide Co., Ltd (Chengdu, China). 14-deoxy-12(R)-sulfo andrographolide (APM) was synthesized and identified by our laboratory according to a previously described procedure24,39 and its purity was above 98%. APS reductase inhibitor bromocriptine mesylate (Inhibitor) was obtained from Shanghai Yuanye Biotechnology Co., Ltd (Shanghai, China). S-adenosyl methionine (S-AMe, as a positive drug), α-naphthylisothiocyanate (ANIT), sodium sulfate (Na2SO4), and sodium bisulfite (NaHSO3) were from Sigma-Aldrich Company (St. Louis, MO, USA). Ampicillin, vancomycin, neomycin sulfate, metronidazole, ethylenediaminetetraacetic acid disodium salt (EDTA-Na2), basic fuchsin were from Dalian Meilun Biotechnology Co., Ltd (Dalian, China). Pooled rat liver S9 fraction (RLS9) and pooled mouse liver S9 fraction (MLS9) were prepared based on the method used in previous study in our lab.40 Pooled human liver S9 fraction (HLS9) was from BD Biosciences (Woburn, MA). Other reagents were of analytical grade or better and were commercially available.

2.2. Animals

C57BL/6J male mice (18–22 g, 6–8 weeks old) and Sprague-Dawley (SD) male rats (180–220 g, 6–8 weeks old) were obtained from the Experimental Animal Center of Southern Medical University (Guangzhou, China) (License number: SCXK(Yue)-2016–0041). Germ-free C57BL/6J male mice (18–22 g, 6–8 weeks old) were obtained from the Experimental Animal Center of the First Affiliated Hospital of Sun Yat-sen University (Guangzhou, China) (License number: SYXK(Yue)-2020–0233). The wild-type (WT) mice and the FXR knockout (Fxr-/-) mice (18–22 g, 6–8 weeks old) were both on a C57BL/6J genetic background.41 All animals were housed in a specific pathogen-free (SPF) room with a temperature of 25–26°C and a 12-h light/dark cycle. Germ-free mice were raised in special sterile isolators. All animals underwent a one-week acclimation period to adapt to these conditions before commencing the experiments. During the experimental period, all animals were provided unrestricted access to both food and autoclaved water. This study was conducted in strict adherence to the “Guide for the Care and Use of Laboratory Animals” and received approval from the Institutional Animal Ethics Committee of Southern Medical University (Guangzhou, China).

2.3. Animals treatments

2.3.1. Antibiotics (ABX) treatment

C57BL/6J male mice and SD male rats were treated with antibiotics (ABX) to deplete gut microbiota. All animals were given autoclaved water supplemented with ampicillin (1 g/L) and received a daily oral gavage of an antibiotic mixture consisting of vancomycin (50 mg/kg), neomycin sulfate (100 mg/kg), and metronidazole (100 mg/kg) for a duration of 2 weeks. The establishment of the ABX (pseudo-germ-free) model has been successfully identified by the cultivation of viable bacteria (Supplementary Figure S1) and 16S rRNA sequencing.42

2.3.2. Desulfovibrio piger (D. piger) culture and treatment

D. piger (DSM 749) obtained from Ningbo MingZhou Bio Co.,Ltd. (Zhejiang, China) was cultured at 37°C in DSMZ Medium 63 (ShanDong Tuopu Bio-engineering Co., Ltd.) following the supplier’s instructions. The concentration of D. piger was measured at OD600 using a Tecan Infinite M1000 PRO Multifunctional Enzyme Labeler (Tecan, Switzerland). For D. piger species colonization experiment, bacterial cells were harvested by centrifugation and resuspended in PBS. Mice were orally administered 0.2 mL of the D. piger suspension (1 × 109 CFU/mL) daily for 2 weeks. The colonization of D. piger in mice was evaluated using real-time quantitative PCR.43

2.3.3. Pharmacodynamic experiments

To investigate the impact of gut microbiota on the anti-cholestatic effect of AP, Control and ABX mice were divided into five groups, respectively: Vehicle, AP (200 mg/kg), ANIT (50 mg/kg), ANIT+AP (200 mg/kg), and ANIT+S-AMe (100 mg/kg).

Similarly, to investigate the impact of C-sulfonate metabolism on the anti-cholestatic effect of AP, mice were randomly divided into two groups: p.o. group (oral administration) and i.p. group (intraperitoneal injection). The p.o. group was further divided into six subgroups: Vehicle, AP (400 mg/kg), ANIT (50 mg/kg), ANIT+AP (100 mg/kg), ANIT+AP (200 mg/kg), ANIT+AP (400 mg/kg). While the i.p. group was divided into five subgroups: Vehicle, AP (200 mg/kg), ANIT (50 mg/kg), ANIT+AP (100 mg/kg), ANIT+AP (200 mg/kg).

Additionally, to investigate the effect of APS reductase inhibitor on the anti-cholestatic effect of AP, mice were randomly divided into Control and Inhibitor groups, each further divided into five subgroups: Vehicle, AP (200 mg/kg), ANIT (50 mg/kg), ANIT+AP (200 mg/kg), and ANIT+S-AMe (100 mg/kg). The mice received APS reductase inhibitor bromocriptine mesylate by oral gavage (20 mg/kg, p.o.) 2 weeks before AP treatment and continued until the end of the experiment. Furthermore, to investigate the important role of FXR in the anti-cholestatic effect of AP, WT and Fxr-/- mice were, respectively, divided into the following groups: Vehicle, ANIT (50 mg/kg), and ANIT+AP (200 mg/kg).

In pharmacodynamic experiments, mice were administered AP at various dosages through different routes daily for a duration of 2 weeks. Additionally, S-AMe was administered at a dosage of 100 mg/kg via intraperitoneal injection (i.p.) daily for 2 weeks. A single dose of ANIT (50 mg/kg, p.o.) was administered 48 h before euthanizing the mice. At the conclusion of each experiment, the mice were euthanized, and blood samples were collected for biochemical analysis and tissue samples for histological examination and UPLC-MS/MS analysis.

2.3.4. Pharmacokinetics experiments

Control and ABX rats were orally administered a single dose of AP at 120 mg/kg. Blood samples (500 μL) were collected at predetermined time intervals (0, 5, 10, 20, 30, 40, 60, 90, 120, 180, 240, 360, and 480 min). In some cases, APS reductase inhibitor bromocriptine mesylate (20 mg/kg, p.o.) or Na2SO4 (200 mg/kg, p.o.) was given to rats 1 week before AP treatment.

Mice were orally administered APS reductase inhibitor bromocriptine mesylate (20 mg/kg, p.o.) or Na2SO4 (200 mg/kg, p.o.) or saline solution 1 week prior to AP treatment (a single dose of AP at 100 mg/kg, p.o.). Blood samples (50 μL) were collected at predetermined time points (0, 10, 20, 40, 60, 90, 120, 240, 360 min). In specific cases, mice were orally administrated with D. piger suspension (0.2 mL, 1 × 109 CFU/mL) or saline solution by oral gavage for 2 weeks before AP treatment (a single dose of AP at 100 mg/kg, p.o.).

The pharmacokinetic parameters of AP including peak time (Tmax), maximum concentration (Cmax), area under the curve from zero time to the last time point measured (AUC0-t), area under the curve from zero time to infinity (AUC0-∞), and elimination half-life (t1/2) were estimated using DAS 2.0 software (Shanghai, China). The relative oral bioavailability (Fr) of AP was calculated by the following equation:Fr(%)=AUC0−∞of ABX/Inhibitor/SO42−/D.pigerAUC0−∞of Control×100%

2.3.5. Tissue distribution experiments

Control and ABX rats were orally administered a single dose of AP at 120 mg/kg. Control, ABX and germ-free C57BL/6J mice were orally administered a single dose of AP at 100 mg/kg. After 90 min of AP administration, all animals were euthanized, and blood and tissue samples were collected for subsequent analysis.

Mice were orally administered APS reductase inhibitor bromocriptine mesylate (20 mg/kg, p.o.) or Na2SO4 (200 mg/kg, p.o.) or saline solution 1 week before AP treatment (a single dose of AP at 100 mg/kg, p.o.). After 360 min of AP treatment, all mice were euthanized, and tissue samples were collected for subsequent analysis. In specific cases, mice were orally administered D. piger suspension (0.2 mL, 1 × 109 CFU/mL) or saline solution 2 weeks before AP treatment (a single dose of AP at 100 mg/kg, p.o.).

2.4. Co-culturing AP with cecum samples or D. piger

The cecum contents from mice were collected and mixed with sterilized anaerobic medium in a ratio of 1.0 g to 10 mL, with gentle agitation. After filtered, the mixed medium containing the intestinal flora was placed in an anaerobic environment and preincubated at 37°C for 30 min before use. Mixed medium with or without APS reductase inhibitor bromocriptine mesylate (20 μg/mL) was incubated with AP (10 μg/mL), while autoclaved mixed medium without APS reductase inhibitor served as the control. Incubation was performed in a shaking incubator at 37°C for 6 h, with strict maintenance of anaerobic condition throughout the process.44 Similarly, D. piger (0.2 mL, 1 × 108 CFU/mL) with or without inhibitor, was incubated with AP for 0, 3, 6, and 12 h, with autoclaved D. piger without inhibitor serving as the control. The reaction was terminated by adding cold methanol contained IS. And then samples were analyzed by UPLC-MS/MS to detect the concentrations of AP and APM. Concentrations of APS reductase and HSO3− were also determined at the end of incubation.

2.5. Biochemical analysis and histological examination

Biochemical indicators of cholestatic liver injury, including alanine aminotransferase (ALT), aspartate aminotransferase (AST), total bile acid (TBA), and total bilirubin (TBIL), were assessed using assay kits obtained from Nanjing JianCheng Bioengineering Institute (Nanjing, China). The experimental procedures were conducted in accordance with the protocols provided by the manufacturer. Fresh liver tissues were fixed in 10% formalin, and prepared as 4-μm paraffin sections. These sections were stained with hematoxylin and eosin for examination of liver histopathology. The stained tissue sections were then observed and captured under an optical microscope for further analysis. Liver histological changes were quantified using a slightly modified scoring system published previously.45 The scoring system comprised four histological features that were evaluated semi-quantitatively: necrosis and fibrosis (0–4), steatosis (0–3), lobular inflammation (0–2), and hepatocellular ballooning (0–2). The liver necrosis and fibrosis areas were quantified by area fraction analysis (Image J software, V1.8.0, National Institutes of Health, USA).46

2.6. 16S rRNA sequencing

The 16S rRNA sequencing of fecal samples obtained from Control and ABX-treated mice was conducted with slight modifications following our previously described protocol.43 Total microbial genomic DNA was extracted from fecal samples using the PF Mag-Bind Stool DNA Kit (Omega Bio-tek, Georgia, USA). DNA quality and concentration were determined by 1.0% agarose gel electrophoresis using a NanoDrop® ND-2000 spectrophotometer (Thermo Scientific Inc., USA). The hypervariable region V3-V4 of the bacterial 16S rRNA gene was amplified using primer pairs 338F (5’-ACTCCTACGGGAGGCAGCAG-3’) and 806 R (5’-GGACTACHVGGGTWTCTAAT-3’) by the ABI GeneAmp® 9700 PCR thermocycler (ABI, CA, USA). Purified amplicons were pooled in equimolar amounts and paired-end sequenced on an Illumina PE300/PE250 platform (Illumina, San Diego, USA) according to the standard protocols by Majorbio Bio-Pharm Technology Co. Ltd. (Shanghai, China). Raw fastq files were de-multiplexed and quality-filtered using QIIME, and the sequences were matched against a high-quality 16S rRNA sequence from the Green Genes database after trimming the primer, barcode, and chimeras. Then, the optimized sequences were clustered into operational taxonomic units (OTUs) using UPARSE 11 with 97% sequence similarity level. Finally, we performed the data analysis.

2.7. Shotgun metagenomic sequencing

Shotgun metagenomic sequencing was performed using cecum samples from mice after AP treatment (400 mg/kg, oral, 14 d). DNA was extracted by a PF Mag-Bind DNA Kit (Omega Bio-tek, USA) and qualified by 1% Agarose gel electrophoresis. Extracted DNA was fragmented by Covaris M220 (Gene Company Limited, China) and screened at approximately 400 bp fragments for paired-end library construction. NEXTFLEX® Rapid DNA-Seq (Bioo Scientific, USA) was used to build paired-end library construction. The Illumina NovaSeq 6000 (Illumina, USA) sequencing platform was used for metagenomic sequencing. Fastp was used to perform quality cropping on the adapter sequences of the 3’ and 5’ ends of reads, preserving high-quality paired-end reads and single-end reads. BWA (version 0.7.17) was used to construct the reads with the host DNA sequence and remove contaminated reads with high alignment similarity. MEGAHIT (version 1.1.2) was used to splice and assemble the optimized sequence. SOAP aligner (version soap2.21release) was used to compare the high-quality reads of each sample with the non-redundant gene set (95% identity) and calculate the abundance information of genes in the corresponding samples. DIAMOND (version 2.0.11) was used to annotate species taxonomy, and KEGG function. Data was deposited into the NCBI Sequence Read Archive (SRA) database.

2.8. RNA sequencing analysis

RNA in liver tissue was isolated from mice treated with ANIT and ANIT+AP (200 mg/kg, i.p.) and purified by TRIzol (Invitrogen, CA). The isolated RNA was quantified using Qubit 2.0 (Thermo Fisher Scientific, USA) following the manufacturer’s protocol. Libraries were prepared using the NEBNext® UltraTM RNA Library Prep Kit for Illumina® (NEB, USA) and sequenced on a Novaseq 6000 (Illumina, USA) at Novogene’s deep sequencing facility. Mapped reads from each sample were assembled using StringTie (v1.3.3b) in a reference-based approach, and featureCounts was used to quantify mapped reads. Differential expression analysis of two conditions/groups was performed using the DESeq2 R package (1.20.0).

2.9. UPLC-MS/MS analysis

An ACQUITYTM UPLC-MS/MS system equipped with an ACQUITYTM UPLC BEH C18 Column (50 × 2.1 mm I.D., 1.7 μm, Waters, Milford, MA, USA) was utilized for the quantification of AP and APM as previously described.47 Briefly, the detection samples (50 μL) were mixed with 50 μL of IS (50 ng/mL). AP and APM were extracted by 400 μL cold methanol. Supernatants were dried under nitrogen stream. Then, the residue was dissolved with 50% methanol in water, centrifuged for 30 min at 13,000 rpm, and 10 µL of the supernatant was subjected to UPLC-MS/MS analysis. The mobile phases consisted of eluent A (water) and eluent B (methanol) with a gradient elution as follows: 0–0.8 min, 10% B; 0.8–1.2 min, 10–90% B; 1.2–5.0 min, 90% B; 5.0–5.5 min, 90–10% B; 5.5–7.0 min, 10% B. The flow rate was 0.3 mL/min and the injection volume of all samples was 10 μL. The quantification of AP, APM and IS was performed by UPLC-MS/MS using an MRM method with the transitions of m/z 349.1→m/z 287.2 for AP, m/z 413.2→m/z 287.2 for APM, m/z 269→m/z 133 for IS.

2.10. APS reductase and HSO3− concentration determination

For the determination of APS reductase concentration, the supernatant from mouse cecum contents and the D. piger strain were assessed using microorganism APS reductase elisa assay kit (Shanghai Senxingyan Biotechnology Co., Ltd., China), following the instructions provided with the kit. For the determination of HSO3− concentration, cecum contents from mice or rats were weighed and homogenized in five volumes of normal saline (based on tissue weight). The homogenate was then subjected to centrifugation at 13,000 rpm and 4°C for 15 min. The supernatant (or perfusate) was collected and mixed with 0.2% formaldehyde, 0.05% basic fuchsin, and 0.01% EDTA-Na2 in a ratio of 1:1:1:9 under acidic conditions. To establish a standard curve, different concentrations of NaHSO3 were used. Subsequently, the absorbance of the resulting solutions was measured at a wavelength of 560 nm after 40 min.

2.11. Real-time quantitative PCR (qPCR) analysis

Total RNA from frozen livers of Vehicle, AP (200 mg/kg, i.p.), ANIT (50 mg/kg), and ANIT+AP (200 mg/kg, i.p.) groups was extracted using the Animal Total RNA Isolation Kit (FOREGENE, China). cDNA was synthesized using a Color Reverse Transcription Kit (EZBioscience, Roseville, CA, USA). qPCR was performed using SYBR green/Rox qPCR master mix (EZBioscience, Roseville, CA, USA), according to the manufacturer’s instruction. The expression of target genes was normalized to β-actin, and fold changes were calculated using 2–ΔΔCT method. For bacterial identification, genomic DNA of fecal flora was isolated using the TIANamp Stool DNA Kit (TIANGEN, Beijing, China) and qPCR was performed. The relative abundance and expression of Desulfovibrio spp., D. piger and APS reductase were normalized to the amounts of total bacteria (16S). All primers were synthesized by Tsingke Biotechnology Co., Ltd. (Beijing, China) and were listed in Supplementary Table S1.

2.12. Western blot analysis

Liver tissue from Vehicle, AP (200 mg/kg, i.p.), ANIT (50 mg/kg), and ANIT+AP (200 mg/kg, i.p.) groups was homogenized in RIPA buffer (Thermo Scientific). Homogenized liver tissue was then incubated on ice for 1 h and then centrifuged at −4°C for 20 min at 13,000 rpm. The collected supernatant was quantified using BCA protein assay kit (Epizyme). Primary antibodies including anti-FXR (anti-NR1H4, bs-12,867 R, bioss), anti-MRP2 (anti-ABCC2, bs-1092 R, bioss), anti-MRP3 (anti-ABCC3, 39909S, CST), anti-MRP4 (anti-ABCC4, DF6921, Affinity), anti-OATP1 (anti-SLCO1C1, bs -11,436 R, bioss), NTCP Polyclonal Antibody (ABP53103, Abbkine) and GADPH (60004–1-Ig, proteintech) were incubated overnight at 4°C. The sections were incubated for 1 h at room temperature with a secondary antibody diluted 1:5000. The membrane was then washed three times with TBST for 5 min. To visualize the protein bands, the membrane was reacted with Immobilon Forte Western HRP (Millipore) substrate solution for 2 min and analyzed using Tanon5200 imager. Data analysis was performed using ImageJ software (V1.8.0, National Institutes of Health, USA).

2.13. Statistical analysis

Data analysis and graphing were conducted by GraphPad Prism v.8 (GraphPad, USA). Two-tailed unpaired Student’s t-test was used to evaluate the significant difference between two groups, while one-way ANOVA followed by Tukey’s multiple comparisons test was applied for three or more groups. Results were presented as mean ± standard error of the mean (SEM) or mean ± standard deviation (SD) and p ＜ 0.05 was considered statistically significant.

3. Results

3.1. Gut microbiota depletion enhances the anti-cholestatic effect of AP

To investigate the impact of gut microbiota on the anti-cholestatic effect of AP, the anti-cholestatic effect of AP (200 mg/kg, p.o.) was compared between the Control and ABX mice (Figure 1(a)). A significant reduction in both the quantity and diversity of gut microbiota was noted following ABX treatment (Figure 1(b-e) and Supplementary Figure S1), indicating a substantial depletion of the gut microbiome under the influence of ABX. An intrahepatic cholestasis model was induced by ANIT, as demonstrated by the increase in key biochemical indicators (ALT, AST, TBA and TBIL) in plasma, along with the presence of hepatic necrosis and inflammatory cell infiltration (Figure 1(f-j)). In Control mice with intact intestinal flora, oral administration of AP at a dosage of 200 mg/kg showed a weak protective effect against ANIT-induced liver injury (Figure 1(f-j)). Conversely, in ABX mice where intestinal flora were depleted, AP at the same dosage markedly ameliorated ANIT-induced elevation of plasma biochemical indicators, liver cell necrosis, and hepatic inflammatory cell infiltration (Figure 1(f-j)). These findings suggest that gut microbiota impairs the anti-cholestatic effect of AP. Figure 1. Gut microbiota depletion enhances the anti-cholestatic effect of AP. (a) Schematic representation of the procedure for administering antibiotics (ABX) and AP treatment in mice. (b) α-diversity analysis using the Shannon index at the operational taxonomic unit (OTU) level following 16S rRNA sequencing on fecal samples obtained from Control and ABX-treated mice (n = 6 per group). (c) β-diversity analysis using principal coordinates analysis (PCoA) on OTU level (n = 6 per group). (d) Number of shared OTUs between Control and ABX groups (n = 6 per group). (e) Percentage of community on the phylum level between Control and ABX groups (n = 6 per group). (f–i) Plasma levels of ALT, AST, TBA, and TBIL (n = 8 per group). (j) Representative H&E staining of the liver (magnification: 200×) and quantitative analysis of H&E staining score and necrosis area of the liver (n = 5 per group). Data are presented as mean ± SEM. *p < .05, **p < .01, ***p < .001, and ****p < 0.0001.

3.2. C-sulfonate metabolism mediated by gut microbiota impairs the anti-cholestatic effect of AP

To investigate the influence of C-sulfonate metabolism on the anti-cholestatic effect of AP, we administered AP to mice via oral (p.o.) and intraperitoneal (i.p.) routes (Figure 2(a)). Compared to the ANIT group, oral treatment with AP at low dosage of 100 or 200 mg/kg showed week hepatoprotective effect, while at a high dosage of 400 mg/kg showed significant efficacy (Figure 2(b-f)). In contrast, intraperitoneal treatment with AP at low dosage of 100 mg/kg or 200 mg/kg both exhibited a superior anti-cholestatic effect, significantly reducing the severity of hepatic damage (Figure 2(b-f)). These results strongly indicate that intraperitoneal administration of AP has a more potent anti-cholestatic effect compared to oral administration. This suggests that intestinal C-sulfonate metabolism significantly weakens the anti-cholestatic efficacy of AP. Figure 2. Intestinal C-sulfonate metabolism impairs the anti-cholestatic effect of AP. (a) Schematic diagram illustrating oral (p.o.) and intraperitoneal (i.p.) administration of AP at various dosages in mice. (b–e) Plasma levels of ALT, AST, TBA, and TBIL (n = 6 per group). (f) Representative H&E staining of the liver (magnification: 200×) and quantitative analysis of H&E staining score and necrosis area of the liver (n = 5 per group). Data are presented as mean ± SEM. *p < .05, **p < .01, and ****p < .0001.

To explore the impact of gut microbiota on the C-sulfonate metabolism of AP, we employed ABX-treated rat, ABX-treated mouse and germ-free mouse models for comprehensive pharmacokinetic and tissue distribution analyses of AP. In ABX-treated rats, elevated levels of AP and reduced levels of APM were observed in plasma and tissues (Figure 3(a-b)), leading to a substantial 70.11% enhancement in the relative oral bioavailability of AP (Table 1). Similarly, in ABX-treated mice, the concentrations of AP and APM were notably increased and decreased, respectively, suggesting a reduction in the C-sulfonate metabolism of AP in mice with depleted gut microbiota (Figure 3(c)). Notably, in germ-free mice, the C-sulfonate metabolism of AP was also suppressed, consistent with the findings from ABX-treated rats and mice (Figure 3(c)). We further determined the concentrations of AP and APM in the small intestine and liver of Control and ABX-treated mice from previous pharmacological experiments and found that following ABX treatment, the concentration of AP significantly increased in the intestine and liver, while the level of its metabolite APM decreased significantly (Figure 3(d)). Additionally, in comparative pharmacological experiments of AP between oral administration and intraperitoneal injection, when AP exhibits better efficacy, its concentration in the small intestine and liver was correspondingly higher, while the concentration of its metabolite APM was lower (Figure 3(e)). Furthermore, our results indicated that the generation of the C-sulfonate metabolite was unaffected by the involvement of sulfotransferases (SULTs), as evidenced by the absence of detectable production of the metabolite APM in the classical incubation system of SULT-mediated metabolism (Suplementary Figure S2). Taken together, these findings convincingly establish the substantial role of gut microbiota in the C-sulfonate metabolism of AP, consequently diminishing the bioavailability and anti-cholestatic effect of AP. Figure 3. Gut microbiota mediates the C-sulfonate metabolism of AP. (a) Plasma pharmacokinetics of AP and APM following oral administration of AP (120 mg/kg) in rats (n = 4 per group). (b) Distribution of AP and APM in the duodenum (d), jejunum (j), ileum (i), colon (c), liver, and kidneys following oral administration of AP (120 mg/kg) in rats (n = 4 per group). (c) Distribution of AP and APM in the small intestine, colon, liver, and kidneys following oral administration of AP (100 mg/kg) in mice (n = 5 per group). (d) Concentrations of AP and APM in the small intestine and liver following AP treatment (200 mg/kg) in Control and ABX-treated mice (n = 8 per group). (e) Concentrations of AP and APM in the small intestine and liver following AP treatment at different dosages via oral (p.o.) or intraperitoneal (i.p.) administration routes in mice (n = 6 per group). Data are presented as mean ± SEM. *p < .05, **p < .01, ***p < .001, and ****p < .0001.

Table 1. Plasma pharmacokinetics of AP following oral gavage administration at a dose of 120 mg/kg in rats from the Control and ABX groups. The presented values are the mean ± SD (n = 4 per group). *p < 0.05, **p < 0.01, compared to Control group.

Parameters	Units	Control	ABX	
Tmax	min	50.00 ± 11.55	42.50 ± 12.58	
Cmax	μg/mL	0.69 ± 0.19	1.33 ± 0.25**	
t1/2	min	183.94 ± 55.70	138.82 ± 44.04	
AUC0-t	mg/L*min	141.73 ± 28.80	265.51 ± 63.63*	
AUC0-∞	mg/L*min	180.09 ± 42.46	306.35 ± 95.27	
Fr (%)	 	 	170.11	

3.3. D. piger-derived APS reductase mediates the C-sulfonate metabolism of AP

Our previous study demonstrated that large amounts of C-sulfonate metabolite APM were generated when AP was incubated with rat duodenal perfusate.22 Furthermore, we observed continued APM generation even after subjecting the duodenal perfusate to filtration and ultrafiltration (Suplementary Figure S3(a-b)). However, APM was not detectable when the perfusate underwent autoclaving treatment (Supplementary Figure S3(c)), coinciding with a significant reduction in the HSO3− content of the perfusate (Supplementary Figure S3(d)). Remarkably, the addition of HSO3− to the autoclaved perfusate resulted in the production of APM (Supplementary Figure S3(e)). These findings indicate that HSO3− is an essential ligand for the C-sulfonate metabolism of AP, independent of host-derived enzymes participation.

To investigate whether HSO3− originated from the intestine, we administered AP to mice and rats via the tail vein and orally, respectively. APM was undetectable in the intestine, liver, and kidneys after tail vein administration of AP, but a significant amount of APM was detected in all tissues after oral administration of AP (Supplementary Figure S4). This observation suggested that HSO3− primarily originates from the intestine. HSO3− is biosynthesized in vivo via two principal pathways: the host and the gut microbiota pathway. The host pathway involves the cysteine catabolic process, wherein cysteine is enzymatically converted into cysteine sulfinic acid by cysteine dioxygenase, subsequently undergoing decomposition to yield HSO3−.48 The gut microbiota pathway is the dissimilatory sulfate reduction route, wherein sulfate (SO42-) undergoes reduction to form adenosine 5’-phosphosulfate (APS), and APS, in turn, serves as a precursor for the production of HSO3− through the action of APS reductase, an enzyme produced by the gut microbiota (primarily by Desulfovibrio spp.).49 Our research findings demonstrated a significant decrease in HSO3− content in the cecum of ABX-treated rats, ABX-treated mice, and germ-free mice (Figure 4(a-b)), accompanied by a notable reduction in the abundance of Desulfovibrio spp. and APS reductase in ABX-treated mice (Figure 4(d)). Furthermore, administration of the APS reductase inhibitor bromocriptine mesylate, which did not induce damage to organism in our study (Supplementary Figure S5), resulted in a decrease in HSO3− content (Figure 4(c)) and APS reductase level (Figure 4(e)), while the administration of Na2SO4 led to a significant increase in HSO3− content (Figure 4(c)). Inhibition of APS reductase activity by bromocriptine mesylate significantly suppressed the C-sulfonate metabolism of AP, resulting in a remarkable 117.24% increase in AP bioavailability in rats and a 47.48% increase in mice (Figure 4(f-g)and Tables 2,3), accompanied by elevated AP levels in mouse tissues (Figure 4(h)). Conversely, treatment of Na2SO4 enhanced the C-sulfonate metabolism of AP, leading to a reduction in its bioavailability by 49.30% in rats and 70.58% in mice (Figure 4(f-g) and Tables 2,3), coupled with decreased AP levels in mouse tissues (Figure 4(h)). Similarly, a substantial increase in AP C-sulfonate metabolism was observed following the addition of Na2SO4 in an in vivo intestinal perfusion assay conducted in rats (Supplementary Figure S6). Figure 4. Gut microbiota-derived APS reductase mediates the C-sulfonate metabolism of AP. (a-b) Concentration of HSO3− in the cecum of ABX-treated rats, ABX-treated mice and germ-free mice (n = 5 per group). (c) Concentration of HSO3− in the cecum after oral administration of APS reductase inhibitor and Na2SO4 in mice (n = 4 per group). (d) Relative expression of Desulfovibrio spp. and APS reductase in the feces of Control and ABX-treated mice by qPCR (n = 8 per group). (e) Relative expression of Desulfovibrio spp. and APS reductase in the feces of mice after APS reductase inhibitor treatment by qPCR (n = 8 per group). (f) The effect of APS reductase inhibitor and Na2SO4 treatment on plasma pharmacokinetics of AP and APM after oral administration of AP (120 mg/kg) in rats (n = 4 for Control group and Inhibitor group, n = 5 for SO42- group). (g) The effect of APS reductase inhibitor and Na2SO4 treatment on plasma pharmacokinetics of AP and APM after oral administration of AP (100 mg/kg) in mice (n = 10 per group). (h) The effect of APS reductase inhibitor and Na2SO4 treatment on the distribution of AP and APM after oral administration of AP (100 mg/kg) in mice (n = 10 per group). Inhibitor refers to APS reductase inhibitor bromocriptine mesylate treatment and SO42- denotes Na2SO4 treatment. Data are presented as mean ± SEM. *p < .05, **p < .01, ***p < .001, and ****p < .0001.

Table 2. Plasma pharmacokinetics of AP following oral gavage administration at a dose of 120 mg/kg in rats from the Control, Inhibitor and SO42- groups.

Parameters	Units	Control	Inhibitor	SO42-	
Tmax	min	30.00 ± 24.50	35.00 ± 10.00	32.00 ± 4.47	
Cmax	μg/mL	1.42 ± 0.71	2.75 ± 0.71*	0.49 ± 0.19*	
t1/2	min	122.32 ± 23.81	109.32 ± 35.35	262.83 ± 194.05	
AUC0-t	mg/L*min	206.23 ± 59.18	451.39 ± 124.08*	81.66 ± 20.02**	
AUC0-∞	mg/L*min	221.37 ± 58.26	480.90 ± 120.03**	112.23 ± 35.42*	
Fr (%)	 	 	217.24	50.70	
The presented values are the mean ± SD (n = 4 for Control group and Inhibitor group, n = 5 for SO42- group). *p < .05, **p < .01, compared to control group.

Table 3. Plasma pharmacokinetics of AP following oral gavage administration at a dose of 100 mg/kg in mice from the Control, Inhibitor and SO42- groups.

Parameters	Units	Control	Inhibitor	SO42-	
Tmax	min	32.00 ± 32.59	27.00 ± 14.94	17.00 ± 4.83	
Cmax	μg/mL	0.50 ± 0.33	0.97 ± 0.37**	0.142 ± 0.03**	
t1/2	min	145.97 ± 67.38	145.86 ± 86.38	133.19 ± 104.25	
AUC0-t	mg/L*min	38.78 ± 9.13	54.34 ± 19.38*	11.93 ± 3.71****	
AUC0-∞	mg/L*min	46.29 ± 14.77	68.27 ± 30.94	13.62 ± 4.03****	
Fr (%)	 	 	147.48	29.42	
The presented values are the mean ± SD (n = 10 per group). *p < .05, **p < .01, ****p < .0001, compared to control group.

In order to enhance the understanding of the cause–effect relationship between Desulfovibrio spp., AP’s C-sulfonate metabolism and bioavailability, we carried out a series of experiments. For experiments involving the co-culturing of AP with cecum samples from mice, the results unveiled a notable conversion of AP into APM by the intestinal flora of mice (Figure 5(a)). This conversion was effectively impeded upon the addition of an APS reductase inhibitor (Figure 5(a)). Moreover, concomitant with this inhibition, there was a significant decrease in the concentrations of APS reductase and HSO3− (Figure 5(b)), suggesting that APS reductase produced by gut microbiota facilitates the C-sulfonate metabolism of AP. To further clarify the specific Desulfovibrio species responsible for APS reductase production and involvement in AP’s C-sulfonate metabolism, we conducted shotgun metagenomics analysis. We discovered four predominant Desulfovibrio species (D. piger, D. fairfieldensis, D. desulfuricans and D. vulgaris) in both Control and AP-treated mice. Subsequent analysis revealed that D. piger exhibited the highest abundance in both groups (Figure 5(c)). Notably, D. piger is one of the most prevalent and abundant Desulfovibrio species in the human body.50–52 Consequently, we proceeded to co-culture AP with D. piger and observed its time-dependent ability to metabolize AP into APM in vitro through the production of APS reductase and HSO3−. This activity was significantly inhibited by an APS reductase inhibitor (Figure 5(d-e)). To further validate the role of D. piger, we successfully transplanted it into mice. This resulted in a significant increase in the relative abundance of APS reductase and HSO3− concentration (Figure 5(f-g)), leading to a marked enhancement in the C-sulfonate metabolism of AP and subsequently reducing its bioavailability by 60.11% (Figure 5(h-i), Table 4). Furthermore, following colonization of D. piger in mice, there was a notable reduction in AP concentration in tissues such as the liver, along with an increase in the concentration of APM (Figure 5(i)). Figure 5. D. piger-derived APS reductase mediates the C-sulfonate metabolism of AP. (a) Concentrations of AP and APM following co-cultured of AP with cecum samples from mice in the absence or presence of APS reductase inhibitor (n = 5 per group). (b) Concentrations of APS reductase and HSO3− following co-cultured of AP with cecum samples from mice in the absence or presence of APS reductase inhibitor (n = 5 per group). (c) Relative abundance of four Desulfovibrio species obtained from shotgun metagenomic sequencing analysis of cecum samples in Control and AP-treated mice (n = 8 per group). (d) Concentrations of AP and APM following co-cultured of AP with D. piger in the absence or presence of APS reductase inhibitor (n = 5 per group). (e) Concentrations of APS reductase and HSO3− following co-cultured of AP with D. piger in the absence or presence of APS reductase inhibitor (n = 5 per group). (f) Relative abundance of Desulfovibrio spp. and D. piger in mouse fecal samples post-colonization of D. piger using qPCR (n = 12 per group). (g) Relative mRNA level of APS reductase and concentration of HSO3− in mouse fecal samples post-colonization of D. piger (n = 12 per group). (h) The effect of D. piger treatment on plasma pharmacokinetics of AP and APM after oral administration of AP (100 mg/kg) in mice (n = 12 per group). (i) The effect of D. piger treatment on the distribution of AP and APM after oral administration of AP (100 mg/kg) in mice (n = 12 per group). Data are presented as mean ± SEM. *p < .05, **p < .01, ***p < .001, and ****p < .0001.

Table 4. Plasma pharmacokinetics of AP following oral gavage administration at a dose of 100 mg/kg in mice from the Control and D. piger groups.

Parameters	Units	Control	D. piger	
Tmax	min	41.67 ± 39.27	35.00 ± 31.77	
Cmax	μg/mL	0.76 ± 0.44	0.36 ± 0.28*	
t1/2	min	135.90 ± 65.22	77.18 ± 38.61*	
AUC0-t	mg/L*min	54.56 ± 13.99	22.45 ± 7.65****	
AUC0-∞	mg/L*min	61.80 ± 15.00	24.65 ± 7.04****	
Fr (%)	 	 	39.89	
The presented values are the mean ± SD (n = 12 per group). *p < .05, ****p < .0001, compared to Control group.

These collective findings underscore the crucial role of APS reductase, produced by D. piger, in mediating the C-sulfonate metabolism of AP and modulating its bioavailability.

3.4. Inhibition of APS reductase enhances the anti-cholestatic effect of AP

To ascertain the impact of APS reductase-mediated C-sulfonate metabolism on the anti-cholestatic efficacy of AP, we conducted a treatment study in mice wherein AP was administered in combination with APS reductase inhibitor bromocriptine mesylate (Figure 6(a)). Our research findings revealed that the administration of AP (200 mg/kg, p.o.) alone did not exhibit significant alleviation of ANIT-induced cholestatic liver injury (Figure 6(b-f)). However, when AP was administered in combination with APS reductase inhibitor, a notable amelioration of liver injury as observed (Figure 6(b-f)). Furthermore, the APS reductase inhibitor effectively suppressed the C-sulfonate metabolism of AP, resulting in an increase in AP exposure and a corresponding decrease in APM within the small intestine and liver (Figure 6(g)). Based on these results, we conclude that the anti-cholestatic effect of AP can be enhanced through inhibition of APS reductase, which mediates the C-sulfonate metabolism of AP. Consequently, the co-administration of AP with APS reductase inhibitors holds promise as a robust and efficacious therapeutic strategy for cholestasis. Figure 6. APS reductase inhibitor enhances the anti-cholestasis effect of AP. (a) Schematic representation of the treatment procedure involving APS reductase inhibitor (bromocriptine mesylate, 20 mg/kg, p.o.) and AP (200 mg/kg, p.o.) in mice. (b–e) Plasma levels of ALT, AST, TBA, and TBIL (n = 8 per group). (f) Representative H&E staining of the liver (magnification: 200×) and quantitative analysis of H&E staining score and necrosis area of the liver (n = 5 per group). (g) Concentrations of AP and APM in the small intestine and liver following AP treatment in the control and APS reductase inhibitor-treated mice (n = 8 per group). Data are presented as mean ± SEM. **p < .01, ***p < .001, and ****p < .0001.

3.5. AP exerts anti-cholestatic effects via activating FXR-related pathway

To further investigate the mechanism by which AP alleviates cholestasis, we performed RNA-sequencing analysis to identify the potential targets and pathways of AP for the treatment of cholestasis. Previous studies have shown that AP exerts its anticholestatic effects mainly through the NF-κB and JNK/MAPK pathways.15,53,54 However, our analysis revealed no significant differences in the clustering of both NF-κB pathway genes (Vcam-1, Il-4, Il-6, Tnf, and Ptgs2) and JNK/MAPK pathway genes (Tnf, Il-1β, Timp-1, Mmp-2, Bax, and Bcl-2) (Supplementary Figure S7). Nevertheless, our findings showed a significant upregulation in the expression of Fxr (Nr1h4) (Figure 7(a-b)). Consistently, both qPCR and Western blot analyses confirmed the upregulation of FXR expression after AP treatment (Figure 7(c-d)). FXR plays a crucial role in regulating bile acid transporters, thus exerting an important function in anti-cholestatic processes.41,55,56 Subsequently, we employed both WT and Fxr-/- mice in our study to test whether AP exerts anti-cholestatic effects via activating FXR-related pathway (Figure 7(e)). In WT mice, AP exhibited a superior anti-cholestatic effect (Figure 7(f-j)). However, in Fxr-/- mice, the beneficial effect of AP was completely abolished at the same dosage (), highlighting the dependence of AP’s anti-cholestatic effect on FXR (Figure 7(f-j)). Furthermore, we observed that AP heightened the expression of FXR downstream bile acid transporters (Mrp2, Mrp3, Mrp4, Oatp1, and Ntcp, except for Bsep, Ost-α and Ost-β) (Figure 8(a-h)). Western blot analysis showed that the expression of MRP2, MRP4, and OATP1 was consistent with qPCR analysis (Figure 8(i-l)). In conclusion, these findings suggest that AP exerts its anti-cholestatic effect by activating FXR and thereby upregulating MRP2 and MRP4.Figure 7. AP exerts anti-cholestatic effects by agonizing FXR. (a) Volcano plot depicting the differentially expressed genes between the ANIT group and the ANIT+AP group, based on RNA sequencing analysis conducted on liver samples from mice treated with ANIT and ANIT+AP (n = 3 per group). (b) The read counts of Fxr obtained from RNA sequencing analysis on liver samples collected from mice treated with ANIT and ANIT+AP (n = 3 per group). (c) Relative mRNA level of Fxr in the liver of mice as determined by qPCR (n = 5 per group). (d) Relative protein level of FXR in the liver of mice as assessed by western blot analysis (n = 3 per group). (e) Schematic diagram illustrating administration of AP in WT and Fxr−/− ice. (f-i) Plasma levels of ALT, AST, TBA, and TBIL (n = 5 per group). (j) Representative H&E staining of the liver (magnification: 200×) and quantitative analysis of H&E staining score and necrosis area of the liver (n = 3 per group). Data are presented as mean ± SEM. *p < .05, **p < .01, ***p < .001, and ****p < .0001.

Figure 8. AP upregulates the relative mRNA and protein levels of FXR downstream bile acid transporters. (a–h) Relative mRNA levels of Mrp2, Mrp3, Mrp4, Oatp1, Ntcp, Bsep, Ost-α and Ost-β in the liver of mice by qPCR (n = 5 per group). (i–l) Relative protein levels of MRP2, MRP3, MRP4, OATP1 and NTCP in the liver of mice by western blot analysis (n = 3 per group). Data are presented as mean ± SEM. *p < .05, **p < .01, ***p < .001, and ****p < .0001.

4. Discussion

Cholestatic liver injury, characterized by the accumulation of toxic bile acids in the liver, presents a significant challenge in clinical management due to the lack of effective therapeutic options.57,58 Despite previous studies demonstrating the protective potential of AP against cholestatic liver injury,15,16 its efficacy may be hindered by low bioavailability.22 Our study focused on exploring the underlying mechanisms of AP’s low bioavailability, shedding light on the critical role played by gut microbiota in modulating its pharmacokinetics and anti-cholestatic effect.

One of the primary findings of our study, utilizing pseudo- and germ-free models, was the discovery that gut microbiota significantly influences the metabolism of AP, resulting in the formation of a C-sulfonate metabolite, APM. This transformative metabolic pathway is mediated by a specific gut microbial enzyme, APS reductase, which catalyzes the reduction of SO42- to HSO3−. HSO3− subsequently interacts with AP, targeting its C=C unsaturated double bond, thereby leading to the formation of the C-sulfonate metabolite. It is important to note that C-sulfonate metabolism constitutes a distinctive pathway, separate from the O- and N-sulfonate metabolism pathway governed by SULTs.29,30 Significantly, this study represents the first documentation of gut microbiota, especially D. piger-derived APS reductase’s pivotal involvement in directing the C-sulfonate metabolism of AP, thereby introducing a novel dimension to the already intricate network of drug metabolic pathways. In our study, we found that clearing gut microbiota or inhibiting APS reductase can significantly improve the bioavailability of AP. By depleting gut microbiota using ABX, the relative bioavailability of AP increased by 70%, while inhibiting APS reductase activity led to a relative bioavailability increase of 117% in rats. It is noteworthy that a multitude of pharmaceutical agents undergo metabolism not only by the host but also by the gut microbiota. This complexity is exemplified by various drugs, such as acarbose,38 aspirin,36 and levodopa,59 indicating the intricate interplay of drug metabolism within the gut microbiota. Notably, HSO3− is generated in vivo by both the host and the intestinal flora.48,60 While our research elucidated the vital role of intestinal flora-produced HSO3− in AP metabolism, it raises the intriguing question of whether host-derived HSO3− plays a role in this process, warranting further exploration and investigation.

Our study also demonstrated a remarkable enhancement in AP’s anti-cholestatic effect upon inhibition of APS reductase produced by Desulfovibrio spp. (particularly D. piger). D. piger is one of the most common and abundant Desulfovibrio species in the intestinal tract of both mice and humans.50 Our investigation revealed that D. piger metabolizes AP through the production of APS reductase in vitro and in vivo, thereby impeding the bioavailability of AP. Therefore, it is highly probable that D. piger affects AP’s C-sulfonate metabolism in the human body. This finding suggests a potential strategy for optimizing the therapeutic application of AP in patients suffering from cholestatic liver injury. APS reductase, predominantly produced by sulfate-reducing bacteria including Desulfovibrio spp.,49 serves as a pivotal enzyme within the sulfate reduction metabolic pathway, playing a vital role for the growth and development of both plants and animals.61 It is notable that APS reductase has garnered attention as a validated target for the development of novel antitubercular agents.62 Several compounds with inhibitory effects on APS reductase have been investigated biologically, including bromocriptine, flupirtine, and mexiletine.63,64 These compounds have already found clinical utility in the treatment of conditions such as Parkinson’s disease, acute pain, and migraines, respectively65–67 In our study, the co-administration of AP with the APS reductase inhibitor bromocriptine significantly enhanced the anti-cholestatic effect of AP. Thus, combining these inhibitors with AP holds substantial promise for enhancing the clinical efficacy of AP in the treatment of cholestatic liver injury. Additionally, strategies aimed at modifying gut microbiota composition and subsequently reducing APS reductase activity will also enhance the clinical efficacy of AP in the treatment of cholestatic liver injury. These strategies involve utilizing probiotics like Lactobacillus and Bifidobacterium, consuming prebiotic-rich foods to promote beneficial bacteria growth, limiting sulfur intake to hinder sulfate-reducing bacteria, and in some cases, employing targeted antibiotic therapy.68–71 Our shotgun metagenomic analysis revealed that AP treatment did not induce significant alterations in the expression of APS reductase or inhibit the growth of APS reductase-producing gut microbiota, despite its ability to influence the abundance and composition of gut microbes (Supplementary Figure S8(a-f)). The Xiaoyanlidan formula, with AP as its quality control component, exhibits promising clinical efficacy in addressing cholestatic liver injury.17,18 Further investigation is warranted to ascertain whether specific components within it have the potential to inhibit APS reductase directly or impede the proliferation of APS reductase-producing bacterial strains.

Furthermore, our investigation provided intriguing insights into the mechanisms underlying the anti-cholestatic effects of AP. While previous studies suggested that AP exerted its anti-cholestatic effects primarily through the NF-κB and JNK/MAPK pathways,15,53,54 our transcriptomic results revealed no significant effect of AP on these pathways. Instead, we observed upregulation of FXR following AP treatment, confirmed by qPCR and WB analysis. Further, we compared the anti-cholestatic effects of AP in wild-type mice and FXR knockout mice, observing significant effects in wild-type mice but no protective effect in knockout mice. FXR, a nuclear receptor, acts as a key regulator of bile acid synthesis, transport, and metabolism.55 Specifically, FXR activation regulates the expression of bile acid transporters, such as BSEP, MRP2/3/4.72 BSEP and MRP2 facilitate the secretion of bile acids into bile, while MRP3 and MRP4 responsible for transporting bile acids out of hepatocytes into the bloodstream.73 We found that AP upregulated the gene and protein levels of MRP2 and MRP4 in our pharmacological model, consistent with a previous study.15 Our results also indicated that AP upregulated OATP1 expression, while the expression of OATP1 regulated by FXR may vary in different contexts.74,75 Traditional Chinese medicines like Yin Chen Hao Decoction and Gardenia jasminoides Ellis have demonstrated their ability to ameliorate cholestatic liver injury through FXR-mediated mechanisms.76,77 Currently, commercially available dosage forms of AP include tablets, capsules, and dripping pills. The inherent qualities of AP, such as its well-tolerated nature, high safety profile, absence of side effects or adverse reactions, and cost-effectiveness,78 position it as a potent therapeutic agent for cholestasis. Previous studies have underscored a close association between cholestatic diseases and the gut microbiota, with specific focus on Enterococcus, Lactobacillus, and Escherichia bacteria, which are deemed pivotal in their onset and progression.79,80 Nonetheless, our shotgun metagenomic analysis unveiled that AP did not suppress the expression of these bacteria (Supplementary Figure S8(g-i)), despite eliciting changes in the composition of the gut microbiota. However, whether AP exerts its anti-cholestatic effects by intervening with other bacteria requires further investigation.

In conclusion, our study has unraveled the intricate relationship between gut microbiota especially Desulfovibrio piger and the C-sulfonate metabolism of AP, shedding light on a novel pathway for optimizing its therapeutic effects in cholestatic liver injury (Figure 9). The identification of APS reductase as a key enzyme and the potential for its inhibition to enhance AP bioavailability hold promise for future drug development efforts. Furthermore, the engagement of FXR-mediated MRP2/4 in AP’s mechanism of action enhances our comprehension of its therapeutic potential. These findings offer valuable insights and open new avenues for improving the clinical application of AP as a promising treatment for cholestatic liver injury.Figure 9. The APS reductase secreted by gut microbiota, particularly Desulfovibrio piger, mediates the C-sulfonate metabolism of AP, thereby impairing its bioavailability and anti-cholestatic efficacy.

5. Conclusion

The C-sulfonate metabolism of AP is facilitated by APS reductase, which is produced by gut microbiota, particularly Desulfovibrio piger. Targeting APS reductase leads to in an increase of the bioavailability and anti-cholestatic efficacy of AP.

Supplementary Material

Supplemental Material

Acknowledgments

The Fxr knockout mice utilized in this study were generously provided by Professor Haiping Hao from China Pharmaceutical University.

Disclosure statement

No potential conflict of interest was reported by the author(s).

Authors contributions

D.T. and W.H. conducted the experiments with assistance from B.F., X.Z., G.Y., C.X., H.W., X.G., Q.Z. and L.Y.; D.T. and G.Y. analyzed data; L.Y. and Z.L. supervised the project; L.Y., Z.L., D.T. and W.H. wrote and revised the manuscript.

Data availability statement

The 16S rRNA and shotgun metagenomic sequencing datasets are available in the NCBI SRA database under the accession numbers PRJNA1073304 and PRJNA1073453, respectively. Additionally, the RNA sequencing data have been deposited in the NCBI SRA database under the accession number PRJNA1073703. The remaining data are available within the article, Supplementary Materials, or Source data file.

Supplementary material

Supplemental data for this article can be accessed online at https://doi.org/10.1080/19490976.2024.2387402.
==== Refs
References

1. Hirschfield GM, Dyson JK, Alexander GJM, Chapman MH, Collier J, Hübscher S, Patanwala I, Pereira SP, Thain C, Thorburn D, et al. The British society of gastroenterology/UK-PBC primary biliary cholangitis treatment and management guidelines. Gut. 2018;67 (9 ):1568–27. doi:10.1136/gutjnl-2017-315259.29593060
2. Carey EJ, Lindor KD. Current pharmacotherapy for cholestatic liver disease. Expert Opin Pharmacother. 2012;13 (17 ):2473–2484. doi:10.1517/14656566.2012.736491.23094715
3. Bortolini M, Almasio P, Bray G, Budillon G, Coltorti M, Frezza M, Okolicsanyi L, Salvagnini M, Williams R. Multicentre survey of the prevalence of intrahepatic cholestasis in 2520 consecutive patients with newly diagnosed chronic liver disease. Drug Invest. 1992;4 (4 ):83–89. doi:10.1007/BF03258368.
4. Floreani A, Mangini C. Primary biliary cholangitis: old and novel therapy. Eur J Intern Med. 2018;47 :1–5. doi:10.1016/j.ejim.2017.06.020.28669591
5. Lleo A, Wang GQ, Gershwin ME, Hirschfield GM. Primary biliary cholangitis. Lancet. 2020;396 (10266 ):1915–1926. doi:10.1016/s0140-6736(20)31607-x.33308474
6. Hirschfield GM, Beuers U, Corpechot C, Invernizzi P, Jones D, Marzioni M, Schramm C. EASL clinical practice guidelines: the diagnosis and management of patients with primary biliary cholangitis. J Hepatol. 2017;67 (1 ):145–172. doi:10.1016/j.jhep.2017.03.022.28427765
7. Levy C, Manns M, Hirschfield G. New treatment paradigms in primary biliary cholangitis. Clin Gastroenterol Hepatol. 2023;21 (8 ):2076–2087. doi:10.1016/j.cgh.2023.02.005.36809835
8. Hasegawa S, Yoneda M, Kurita Y, Nogami A, Honda Y, Hosono K, Nakajima A. Cholestatic liver disease: current treatment strategies and new therapeutic agents. Drugs. 2021;81 (10 ):1181–1192. doi:10.1007/s40265-021-01545-7.34142342
9. Wunsch E, Raszeja-Wyszomirska J, Barbier O, Milkiewicz M, Krawczyk M, Milkiewicz P. Effect of S-adenosyl-L-methionine on liver biochemistry and quality of life in patients with primary biliary cholangitis treated with ursodeoxycholic acid. A prospective, open label pilot study. J Gastrointestin Liver Dis. 2018;27 (3 ):273–279. doi:10.15403/jgld.2014.1121.273.icz.30240471
10. Huang D, Zhao YY, Wang RM, Li W, Yuan FY, Yan XL, Yang X, Tang G-H, Yin S, Bi H-C, et al. Natural product-based screening led to the discovery of a novel PXR agonist with anti-cholestasis activity. Acta Pharmacol Sin. 2022;43 (8 ):2139–2146. doi:10.1038/s41401-021-00793-3.34931017
11. Qin X, Wang X, Tian M, Dong Z, Wang J, Wang C, Huang Q. The role of andrographolide in the prevention and treatment of liver diseases. Phytomed: Int J Phytother Phytopharmacology. 2023;109 :154537. doi:10.1016/j.phymed.2022.154537.
12. Lee KC, Chang HH, Chung YH, Lee TY. Andrographolide acts as an anti-inflammatory agent in lps-stimulated RAW264.7 macrophages by inhibiting STAT3-mediated suppression of the nf-κB pathway. J Ethnopharmacol. 2011;135 (3 ):678–684. doi:10.1016/j.jep.2011.03.068.21497192
13. Huang Z, Wu Z, Zhang J, Wang K, Zhao Q, Chen M, Yan S, Guo Q, Ma Y, Ji L, et al. Andrographolide attenuated mct-induced HSOS via regulating NRF2-initiated mitochondrial biogenesis and antioxidant response. Cell Biol Toxicol. 2023;39 (6 ):3269–3285. doi:10.1007/s10565-023-09832-7.37816928
14. Lee TY, Chang HH, Wen CK, Huang TH, Chang YS. Modulation of thioacetamide-induced hepatic inflammations, angiogenesis and fibrosis by andrographolide in mice. J Ethnopharmacol. 2014;158 (Pt A ):423–430. doi:10.1016/j.jep.2014.10.056.25446592
15. Khamphaya T, Chansela P, Piyachaturawat P, Suksamrarn A, Nathanson M, Weerachayaphorn J. Effects of andrographolide on intrahepatic cholestasis induced by alpha-naphthylisothiocyanate in rats. Eur J Pharmacol. 2016;789 :254–264. doi:10.1016/j.ejphar.2016.07.032.27475677
16. Wang L, Cao F, Zhu L, Liu P, Shang Y, Liu W, Dong X, Bao H-D, Gong P, Wang Z-Y, et al. Andrographolide impairs alpha-naphthylisothiocyanate-induced cholestatic liver injury in vivo. J Nat Med. 2019;73 (2 ):388–396. doi:10.1007/s11418-018-01275-3.30617707
17. Zhang K, Wang M, Yao Y, Huang T, Liu F, Zhu C, Lin C. Pharmacokinetic study of seven bioactive components of Xiaoyan Lidan formula in cholestatic and control rats using UPLC-MS/MS. Biomed Pharmacother = Biomedecine & Pharmacotherapie. 2021;139 :111523. doi:10.1016/j.biopha.2021.111523.33831838
18. Shen B, Yin H, Zhu XP, Xiao WX, Zhou HJ, Xiao G-Y, Zhou H-J. Laparoscopic cholecystectomy combined with xiaoyanlidan tablets for treatment of acute cholecystitis patients: curative effect and impact on prognosis and serological indicators. World Chin J Digestology. 2018;26 (25 ):1523–1528. doi:10.11569/wcjd.v26.i25.1523.
19. Zhang R, Huang T, Zhang Q, Yao Y, Liu C, Lin C, Zhu C. Xiaoyan lidan formula ameliorates α-naphthylisothiocyanate-induced intrahepatic cholestatic liver injury in rats as revealed by non-targeted and targeted metabolomics. J Pharm Biomed Anal. 2020;179 :112966. doi:10.1016/j.jpba.2019.112966.31735473
20. Zeng B, Wei A, Zhou Q, Yuan M, Lei K, Liu Y, Song J, Guo L, Ye Q. Andrographolide: a review of its pharmacology, pharmacokinetics, toxicity and clinical trials and pharmaceutical researches. Phytotherapy Res: PTR. 2022;36 (1 ):336–364. doi:10.1002/ptr.7324.
21. Chen H, Huang C, Li C, Lin A, Huang Y, Wang T, Yao H-T, Lii C-K. Bioavailability of andrographolide and protection against carbon tetrachloride-induced oxidative damage in rats. Toxicol Appl Pharmacol. 2014;280 (1 ):1–9. doi:10.1016/j.taap.2014.07.024.25110055
22. Ye L, Wang T, Tang L, Liu W, Yang Z, Zhou J, Zheng Z, Cai Z, Hu M, Liu Z, et al. Poor oral bioavailability of a promising anticancer agent andrographolide is due to extensive metabolism and efflux by P-glycoprotein. J Pharm Sci. 2011;100 (11 ):5007–5017. doi:10.1002/jps.22693.21721007
23. Yang T, Sheng HH, Feng NP, Wei H, Wang ZT, Wang CH. Preparation of andrographolide-loaded solid lipid nanoparticles and their in vitro and in vivo evaluations: characteristics, release, absorption, transports, pharmacokinetics, and antihyperlipidemic activity. J Pharm Sci. 2013;102 (12 ):4414–4425. doi:10.1002/jps.23758.24166599
24. He X, Li J, Gao H, Qiu F, Hu K, Cui X, Yao X. Identification of a rare sulfonic acid metabolite of andrographolide in rats. Drug Metab Dispos. 2003;31 (8 ):983–985. doi:10.1124/dmd.31.8.983.12867485
25. Liu H, Wang K, Xia G, Wang K, Chai L, Donkor PO, Ding L, Qiu F. Two sulfonate metabolites of physalin a in rats. Xenobiotica. 2018;48 (1 ):11–17. doi:10.1080/00498254.2016.1271961.28001472
26. Wikberg T, Taskinen J. Identification of major metabolites of the catechol-O-methyltransferase inhibitor entacapone in rats and humans. Drug Metab Dispos. 1993;21 (2 ):81–92. doi:10.1007/bf03190186.8095232
27. Adams MD, Wedig JH, Jordan RL, Smith LW, Henderson R, Borzelleca JF. Urinary excretion and metabolism of salts of 2-pyridinethiol-1-oxide following intravenous administration to female Yorkshire pigs. Toxicol Appl Pharmacol. 1976;36 (3 ):523–531. doi:10.1016/0041-008x(76)90231-3.941152
28. Yoshino H, Matsunaga H, Kaneko H, Yoshitake A, Nakatsuka I, Yamada H. Metabolism of N-[4-chloro-2-fluoro-5-[(1-methyl-2-propynyl)oxy]phenyl]-3,4,5,6-tetrahydrophthalimide (S-23121) in the rat: I. Identification of a new, sulphonic acid type of conjugate. Xenobiotica. 2008;23 (6 ):609–619. doi:10.3109/00498259309059399.
29. Meerman MJHN. Sulfation and glucuronidation as competing pathways in the metabolism of hydroxamic acids: the role of N,O-sulfonation in chemical carcinogenesis of aromatic amines. Environ Health Perspect. 1983;49 :27–32. doi:10.1289/ehp.834927.6339226
30. Huo ZP, Feng XC, Wang Y, Tian YT, Qiu F. Sulfite as the substrate of C-sulfonate metabolism of α, β-unsaturated carbonyl containing andrographolide: analysis of sulfite in rats’ intestinal tract and the reaction kinetics of andrographolide with sulfite. Chin J Nat Med. 2021;19 (9 ):706–712. doi:10.1016/s1875-5364(21)60094-8.34561083
31. Carbonero F, Benefiel AC, Alizadeh-Ghamsari AH, Gaskins HR. Microbial pathways in colonic sulfur metabolism and links with health and disease. Front Physiol. 2012;3 :448. doi:10.3389/fphys.2012.00448.23226130
32. Parvez MM, Basit A, Jariwala PB, Gáborik Z, Kis E, Heyward S, Redinbo MR, Prasad B. Quantitative investigation of irinotecan metabolism, transport, and gut microbiome activation. Drug Metab Dispos. 2021;49 (8 ):683–693. doi:10.1124/dmd.121.000476.34074730
33. Haiser HJ, Seim KL, Balskus EP, Turnbaugh PJ. Mechanistic insight into digoxin inactivation by Eggerthella lenta augments our understanding of its pharmacokinetics. Gut Microbes. 2014;5 (2 ):233–238. doi:10.4161/gmic.27915.24637603
34. Peppercorn MA, Goldman P. The role of intestinal bacteria in the metabolism of salicylazosulfapyridine. J Pharmacol Exp Ther. 1972;181 (3 ):555–562.4402374
35. Strong HA, Renwick AG, George CF, Liu YF, Hill MJ. The reduction of sulphinpyrazone and sulindac by intestinal bacteria. Xenobiotica. 1987;17 (6 ):685–696. doi:10.3109/00498258709043976.3630204
36. Zhao R, Coker OO, Wu J, Zhou Y, Zhao L, Nakatsu G, Bian X, Wei H, Chan AWH, Sung JJY, et al. Aspirin reduces colorectal tumor development in mice and gut microbes reduce its bioavailability and chemopreventive effects. Gastroenterology. 2020;159 (3 ):969–983.e964. doi:10.1053/j.gastro.2020.05.004.32387495
37. Mehta RS, Mayers JR, Zhang Y, Bhosle A, Glasser NR, Nguyen LH, Ma W, Bae S, Branck T, Song K, et al. Gut microbial metabolism of 5-ASA diminishes its clinical efficacy in inflammatory bowel disease. Nat Med. 2023;29 (3 ):700–709. doi:10.1038/s41591-023-02217-7.36823301
38. Balaich J, Estrella M, Wu G, Jeffrey PD, Biswas A, Zhao L, Korennykh A, Donia MS. The human microbiome encodes resistance to the antidiabetic drug acarbose. Nature. 2021;600 (7887 ):110–115. doi:10.1038/s41586-021-04091-0.34819672
39. Hong L. Synthesis method of andrographolide sodium bisulfite. 2007. CN101362740.
40. Zhang Q, Zhu L, Gong X, Ruan Y, Yu J, Jiang H, Wang Y, Qi X, Lu L, Liu Z, et al. Sulfonation disposition of Acacetin: In Vitro and in vivo. J Agric Food Chem. 2017;65 (24 ):4921–4931. doi:10.1021/acs.jafc.7b00854.28540728
41. Sinal CJ, Tohkin M, Miyata M, Ward JM, Lambert G, Gonzalez FJ. Targeted disruption of the nuclear receptor FXR/BAR impairs bile acid and lipid homeostasis. Cell. 2000;102 (6 ):731–744. doi:10.1016/s0092-8674(00)00062-3.11030617
42. Zarrinpar A, Chaix A, Xu ZZ, Chang MW, Marotz CA, Saghatelian A, Knight R, Panda S. Antibiotic-induced microbiome depletion alters metabolic homeostasis by affecting gut signaling and colonic metabolism. Nat Commun. 2018;9 (1 ):2872. doi:10.1038/s41467-018-05336-9.30030441
43. Ye L, Hou Y, Hu W, Wang H, Yang R, Zhang Q, Feng Q, Zheng X, Yao G, Hao H, et al. Repressed Blautia-acetate immunological axis underlies breast cancer progression promoted by chronic stress. Nat Commun. 2023;14 (1 ):6160. doi:10.1038/s41467-023-41817-2.37789028
44. Pan L, Yu H, Fu J, Hu J, Xu H, Zhang Z, Bu M, Yang X, Zhang H, Lu J, et al. Berberine ameliorates chronic kidney disease through inhibiting the production of gut-derived uremic toxins in the gut microbiota. Acta Pharm Sin B. 2023;13 (4 ):1537–1553. doi:10.1016/j.apsb.2022.12.010.37139409
45. Kleiner DE, Brunt EM, Van Natta M, Behling C, Contos MJ, Cummings OW, Ferrell LD, Liu Y-C, Torbenson MS, Unalp-Arida A, et al. Design and validation of a histological scoring system for nonalcoholic fatty liver disease. Hepatology. 2005;41 (6 ):1313–1321. doi:10.1002/hep.20701.15915461
46. Li D, Chen Y, Wan M, Mei F, Wang F, Gu P, Zhang X, Wei R, Zeng Y, Zheng H, et al. Oral magnesium prevents acetaminophen-induced acute liver injury by modulating microbial metabolism. Cell Host Microbe. 2024;32 (1 ):48–62.e49. doi:10.1016/j.chom.2023.11.006.38056458
47. Yang R, Hu W, Xie C, Tang D, Zhao X, Fu B, Wu J, Ye L, Liao R. A sensitive liquid chromatography-mass spectrometry method for determination of 14-deoxy-12(R)-sulfo andrographolide concentration in rat plasma and its application to a pharmacokinetic study. Curr Drug Metab. 2022;23 (11 ):905–911. doi:10.2174/1386207326666221020110217.36278441
48. Kohl JB, Mellis AT, Schwarz G. Homeostatic impact of sulfite and hydrogen sulfide on cysteine catabolism. Br J Pharmacol. 2019;176 (4 ):554–570. doi:10.1111/bph.14464.30088670
49. Kushkevych I. Dissimilatory sulfate reduction in the intestinal sulfate-reducing bacteria. Studia Biologica. 2016;10 (N 1 ):197–228. doi:10.30970/sbi.1001.560.
50. Nakao K, Tanaka K, Ichiishi S, Mikamo H, Shibata T, Watanabe K. Susceptibilities of 23 desulfovibrio isolates from humans. Antimicrob Agents Chemother. 2009;53 (12 ):5308–5311. doi:10.1128/aac.00630-09.19786606
51. Chen YR, Jing QL, Chen FL, Zheng H, Chen LD, Yang ZC. Desulfovibrio is not always associated with adverse health effects in the Guangdong gut microbiome project. Peer J. 2021;9 :e12033. doi:10.7717/peerj.12033.34466295
52. Warren YA, Citron DM, Merriam CV, Goldstein EJ. Biochemical differentiation and comparison of desulfovibrio species and other phenotypically similar genera. J Clin Microbiol. 2005;43 (8 ):4041–4045. doi:10.1128/jcm.43.8.4041-4045.2005.16081948
53. Zhang Y, Pan Y, Lin C, Zheng Y, Sun H, Zhang H, Wang J, Yuan M, Duan T, Du Q, et al. Bile acids evoke placental inflammation by activating Gpbar1/NF-κB pathway in intrahepatic cholestasis of pregnancy. J Mol Cell Biol. 2016;8 (6 ):530–541. doi:10.1093/jmcb/mjw025.27402811
54. Lee TY, Lee KC, Chang HH. Modulation of the cannabinoid receptors by andrographolide attenuates hepatic apoptosis following bile duct ligation in rats with fibrosis. Apoptosis. 2010;15 (8 ):904–914. doi:10.1007/s10495-010-0502-z.20446039
55. Lee FY, Lee H, Hubbert ML, Edwards PA, Zhang Y. FXR, a multipurpose nuclear receptor. Trends Biochem Sci. 2006;31 (10 ):572–580. doi:10.1016/j.tibs.2006.08.002.16908160
56. Yuan ZQ, Li KW. Role of farnesoid X receptor in cholestasis. J Dig Dis. 2016;17 (8 ):501–509. doi:10.1111/1751-2980.12378.27383832
57. Alamoudi JA, Li W, Gautam N, Olivera M, Meza J, Mukherjee S, Alnouti Y. Bile acid indices as biomarkers for liver diseases II: the bile acid score survival prognostic model. World J Hepatol. 2021;13 (5 ):543–556. doi:10.4254/wjh.v13.i5.543.34131469
58. Gossard AA, Talwalkar JA. Cholestatic liver disease. Med Clin N Am. 2014;98 (1 ):73–85. doi:10.1016/j.mcna.2013.09.002.24266915
59. van Kessel SP, Frye AK, El-Gendy AO, Castejon M, Keshavarzian A, van Dijk G, van Dijk G, El Aidy S. Gut bacterial tyrosine decarboxylases restrict levels of levodopa in the treatment of Parkinson’s disease. Nat Commun. 2019;10 (1 ):310. doi:10.1038/s41467-019-08294-y.30659181
60. Kushkevych I, Cejnar J, Treml J, Dordević D, Kollar P, Vítězová M. Recent advances in metabolic pathways of sulfate reduction in intestinal bacteria. Cells. 2020;9 (3 ):698. doi:10.3390/cells9030698.32178484
61. Duarte AG, Santos AA, Pereira IA. Electron transfer between the QmoABC membrane complex and adenosine 5’-phosphosulfate reductase. Biochim Biophys Acta. 2016;1857 (4 ):380–386. doi:10.1016/j.bbabio.2016.01.001.26768116
62. Bhave DP, Hong JA, Lee M, Jiang W, Krebs C, Carroll KS. Spectroscopic studies on the [4Fe-4S] cluster in adenosine 5’-phosphosulfate reductase from mycobacterium tuberculosis. J Biol Chem. 2011;286 (2 ):1216–1226. doi:10.1074/jbc.M110.193722.21075841
63. Palde PB, Bhaskar A, Pedró Rosa LE, Madoux F, Chase P, Gupta V, Spicer T, Scampavia L, Singh A, Carroll KS, et al. First-in-class inhibitors of sulfur metabolism with bactericidal activity against non-replicating M. tuberculosis. ACS Chem Biol. 2016;11 (1 ):172–184. doi:10.1021/acschembio.5b00517.26524379
64. Cosconati S, Hong JA, Novellino E, Carroll KS, Goodsell DS, Olson AJ. Structure-based virtual screening and biological evaluation of mycobacterium tuberculosis adenosine 5’-phosphosulfate reductase inhibitors. J Med Chem. 2008;51 (21 ):6627–6630. doi:10.1021/jm800571m.18855373
65. Jansen EN, Meerwaldt JD. Long-term treatment with high-dosage bromocriptine in advanced Parkinson’s disease. Adv Neurol. 1987;45 :539–542.3103396
66. Gupta AK, Kohli M, Pandey PK, Dwivedi PD, Singh VP. Comparative study of analgesia of ketorolac, tramadol, and flupirtine in the management of third molar surgery. Natl J Maxillofac Surg. 2022;13 (2 ):262–268. doi:10.4103/njms.NJMS_261_20.36051803
67. Marmura MJ, Passero FC Jr., Young WB. Mexiletine for refractory chronic daily headache: a report of nine cases. Headache. 2008;48 (10 ):1506–1510. doi:10.1111/j.1526-4610.2008.01234.x.18793209
68. Nanda Kumar NS, Balamurugan R, Jayakanthan K, Pulimood A, Pugazhendhi S, Ramakrishna BS. Probiotic administration alters the gut flora and attenuates colitis in mice administered dextran sodium sulfate. J Gastroenterol Hepatol. 2008;23 (12 ):1834–1839. doi:10.1111/j.1440-1746.2008.05723.x.19120873
69. Dordević D, Jančíková S, Vítězová M, Kushkevych I. Hydrogen sulfide toxicity in the gut environment: meta-analysis of sulfate-reducing and lactic acid bacteria in inflammatory processes. J Adv Res. 2021;27 :55–69. doi:10.1016/j.jare.2020.03.003.33318866
70. Kushkevych I, Kotrsová V, Dordević D, Buňková L, Vítězová M, Amedei A. Hydrogen sulfide effects on the survival of lactobacilli with emphasis on the development of inflammatory bowel diseases. Biomolecules. 2019;9 (12 ):752. doi:10.3390/biom9120752.31756903
71. Barton LL, Ritz NL, Fauque GD, Lin HC. Sulfur cycling and the intestinal microbiome. Dig Dis Sci. 2017;62 (9 ):2241–2257. doi:10.1007/s10620-017-4689-5.28766244
72. Sánchez B. Bile acid-microbiota crosstalk in gastrointestinal inflammation and carcinogenesis: a role for bifidobacteria and lactobacilli? Nat Rev Gastroenterol Hepatol. 2018;15 (4 ):205. doi:10.1038/nrgastro.2018.23.
73. Keppler D. Cholestasis and the role of basolateral efflux pumps. Z Gastroenterol. 2011;49 (12 ):1553–1557. doi:10.1055/s-0031-1281755.22139880
74. Liu M, Zhu D, Wen J, Ding W, Huang S, Xia C, Zhang H, Xiong Y. Berberine promotes OATP1B1 expression and rosuvastatin uptake by inducing nuclear translocation of FXR and LXRα. Front Pharmacol. 2020;11 :375. doi:10.3389/fphar.2020.00375.32292349
75. Qin S, Tian J, Zhao Y, Wang L, Wang J, Liu S, Meng J, Wang F, Liu C, Han J, et al. Gardenia extract protects against intrahepatic cholestasis by regulating bile acid enterohepatic circulation. J Ethnopharmacol. 2024;319 (Pt 1 ):117083. doi:10.1016/j.jep.2023.117083.37634748
76. Wei C, Qiu J, Wu Y, Chen Z, Yu Z, Huang Z, Yang K, Hu H, Liu F. Promising traditional Chinese medicine for the treatment of cholestatic liver disease process (cholestasis, hepatitis, liver fibrosis, liver cirrhosis). J Ethnopharmacol. 2022;297 :115550. doi:10.1016/j.jep.2022.115550.35863612
77. Fang S, Wang T, Li Y, Xue H, Zou J, Cai J, Shi R, Wu J, Ma Y. Gardenia jasminoides Ellis polysaccharide ameliorates cholestatic liver injury by alleviating gut microbiota dysbiosis and inhibiting the TLR4/NF-κB signaling pathway. Int J Biol Macromol. 2022;205 :23–36. doi:10.1016/j.ijbiomac.2022.02.056.35176320
78. Songvut P, Rangkadilok N, Pholphana N, Suriyo T, Panomvana D, Puranajoti P, Akanimanee J, Satayavivad J. Comparative pharmacokinetics and safety evaluation of high dosage regimens of andrographis paniculata aqueous extract after single and multiple oral administration in healthy participants. Front Pharmacol. 2023;14 :1230401. doi:10.3389/fphar.2023.1230401.37663270
79. Isaacs-Ten A, Echeandia M, Moreno-Gonzalez M, Brion A, Goldson A, Philo M, Patterson AM, Parker A, Galduroz M, Baker D, et al. Intestinal microbiome-macrophage crosstalk contributes to cholestatic liver disease by promoting intestinal permeability in mice. Hepatology. 2020;72 (6 ):2090–2108. doi:10.1002/hep.31228.32168395
80. Leung H, Xiong L, Ni Y, Busch A, Bauer M, Press AT, Panagiotou G. Impaired flux of bile acids from the liver to the gut reveals microbiome-immune interactions associated with liver damage. NPJ Biofilms Microbiomes. 2023;9 (1 ):35. doi:10.1038/s41522-023-00398-0.37286586
