Drug-induced Cholestasis Research Articles

Development of a mechanistic biokinetic model for hepatic bile acid handling to predict possible cholestatic effects of drugs

Drug-induced liver injury (DILI) is a common reason for drug withdrawal from the market. An important cause of DILI is drug-induced cholestasis. One of the major players involved in drug-induced cholestasis is the bile salt efflux pump (BSEP; ABCB11). Inhibition of BSEP by drugs potentially leads to cholestasis due to increased (toxic) intrahepatic concentrations of bile acids with subsequent cell injury. In order to investigate the possibilities for in silico prediction of cholestatic effects of drugs, we developed a mechanistic biokinetic model for human liver bile acid handling populated with human in vitro data. For this purpose we considered nine groups of bile acids in the human bile acid pool, i.e. chenodeoxycholic acid, deoxycholic acid, the remaining unconjugated bile acids and the glycine and taurine conjugates of each of the three groups. Michaelis-Menten kinetics of the human uptake transporter Na+-taurocholate cotransporting polypeptide (NTCP; SLC10A1) and BSEP were measured using NTCP-transduced HEK293 cells and membrane vesicles from BSEP-overexpressing HEK293 cells. For in vitro-in vivo scaling, transporter abundance was determined by LC-MS/MS in these HEK293 cells and vesicles as well as in human liver tissue. Other relevant human kinetic parameters were collected from literature, such as portal bile acid levels and composition, bile acid synthesis and amidation rate. Additional empirical scaling was applied by increasing the excretion rate with a factor 2.4 to reach near physiological steady-state intracellular bile acid concentrations (80μM) after exposure to portal vein bile acid levels. Simulations showed that intracellular bile acid concentrations increase 1.7 fold in the presence of the BSEP inhibitors and cholestatic drugs cyclosporin A or glibenclamide, at intrahepatic concentrations of 6.6 and 20μM, respectively. This simplified model provides a tool for a first indication whether drugs at therapeutic concentrations might cause cholestasis by inhibiting BSEP.

Endoplasmic reticulum stress precedes oxidative stress in antibiotic-induced cholestasis and cytotoxicity in human hepatocytes.

Endoplasmic reticulum (ER) stress has been associated with various drug-induced liver lesions but its participation in drug-induced cholestasis remains unclear. We first aimed at analyzing liver damage caused by various hepatotoxic antibiotics, including three penicillinase-resistant antibiotics (PRAs), i.e. flucloxacillin, cloxacillin and nafcillin, as well as trovafloxacin, levofloxacin and erythromycin, using human differentiated HepaRG cells and primary hepatocytes. All these antibiotics caused early cholestatic effects typified by bile canaliculi dilatation and reduced bile acid efflux within 2h and dose-dependent enhanced caspase-3 activity within 24h. PRAs induced the highest cholestatic effects at non cytotoxic concentrations. Then, molecular events involved in these lesions were analyzed. Early accumulation of misfolded proteins revealed by thioflavin-T fluorescence and associated with phosphorylation of the unfolded protein response sensors, eIF2α and/or IRE1α, was evidenced with all tested hepatotoxic antibiotics. Inhibition of ER stress markedly restored bile acid efflux and prevented bile canaliculi dilatation. Downstream of ER stress, ROS were also generated with high antibiotic concentrations. The protective HSP27-PI3K-AKT signaling pathway was activated only in PRA-treated cells and its inhibition increased ROS production and aggravated caspase-3 activity. Overall, our results demonstrate that (i) various antibiotics reported to cause cholestasis and hepatocellular injury in the clinic can also induce such effects in in vitro human hepatocytes; (ii) PRAs cause the strongest cholestatic effects in the absence of cytotoxicity; (iii) cholestatic features occur early through ER stress; (iv) cytotoxic lesions are observed later through ER stress-mediated ROS generation; and (v) activation of the HSP27-PI3K-AKT pathway protects from cytotoxic damage induced by PRAs only.

Intrahepatic cholestasis of pregnancy or azithromycin-induced intrahepatic cholestasis

Rationale:Azithromycin-induced liver injury has been rarely reported in adult individuals, let alone in a pregnant woman. Here, we describe the clinical features and outcomes of azithromycin-induced liver injury in a pregnant woman.Patient concerns:A 30-year-old pregnant woman presented with generalized pruritus and elevated serum bile acid level (123.6 μmol/L) on day 4 of azithromycin administration. A diagnosis of intrahepatic cholestasis of pregnancy was made, and cesarean section was performed immediately. Interestingly, the alanine aminotransferase level (ALT) reached 211.2 U/L on day 9 after azithromycin administration.Diagnosis:Therefore, drug-induced intrahepatic cholestasis was considered.Interventions:(1) Azithromycin withdrawal after the patient hospitalized. (2) Termination of pregnancy by cesarean section was performed inmmediately to protect the fetus. (3) Silymarin capsules and bifendate are used to protect the liver after liver enzymes elevation was discovered.Outcomes:The liver enzymes recovered within 4 weeks without any symptoms after treatment with silymarin capsules and bifendate, which helps reduce ALT level and protects the liver from further injury.Lessions:A pregnant woman developed azithromycin-induced intrahepatic cholestasis. Physicians should be aware of this side effect of azithromycin, which is widely prescribed.

Progressive and Preferential Cellular Accumulation of Hydrophobic Bile Acids Induced by Cholestatic Drugs Is Associated with Inhibition of Their Amidation and Sulfation.

Drug-induced intrahepatic cholestasis is characterized by cellular accumulation of bile acids (BAs), whose mechanisms remain poorly understood. The present study aimed to analyze early and progressive alterations of BA profiles induced by cyclosporine A, chlorpromazine, troglitazone, tolcapone, trovafloxacin, and tacrolimus after 4-hour, 24-hour, and 6-day treatments of differentiated HepaRG cells. In BA-free medium, the potent cholestatic drugs cyclosporine A, chlorpromazine, and troglitazone reduced endogenous BA synthesis after 24 hours, whereas the rarely cholestatic drugs tolcapone, trovafloxacin, and tacrolimus reduced BA synthesis only after 6 days. In the presence of physiologic serum BA concentrations, cyclosporine A, chlorpromazine, and troglitazone induced early and preferential cellular accumulation of unconjugated lithocholic, deoxycholic, and chenodeoxycholic acids that increased 8- to 12-fold and 47- to 50-fold after 24 hours and 6 days, respectively. Accumulation of these hydrophobic BAs resulted from strong inhibition of amidation, and in addition, for lithocholic acid reduction of its sulfoconjugation, and was associated with variable alterations of uptake and efflux transporters. Trovafloxacin also caused BA accumulation, especially after 6 days, whereas tolcapone and tacrolimus were still without effect. However, when exogenous BAs were added to the medium at cholestatic serum concentrations, a 6-day treatment with all drugs resulted in cellular BA accumulation with higher folds of chenodeoxycholic and lithocholic acids. At the tested concentration, tolcapone had the lowest effect. These results bring the first demonstration that major cholestatic drugs can cause preferential and progressive in vitro cellular accumulation of unconjugated toxic hydrophobic BAs and bring new insights into mechanisms involved in drug-induced cellular accumulation of toxic BAs.

A Dynamic Mathematical Model of Bile Acid Clearance in HepaRG Cells.

A dynamic model based on ordinary differential equations that describes uptake, basolateral and canalicular export of taurocholic acid (TCA) in human HepaRG cells is presented. The highly reproducible inter-assay experimental data were used to reliably estimate model parameters. Primary human hepatocytes were similarly evaluated to establish a mathematical model, but with notably higher inter-assay differences in TCA clearance and bile canaliculi dynamics. By use of the HepaRG cell line, the simultaneous TCA clearance associated to basolateral uptake, canalicular and sinusoidal efflux, was predicted. The mathematical model accurately reproduced the dose-dependent inhibition of TCA clearance in the presence and absence of the prototypical cholestatic drugs cyclosporine A (CsA) and chlorpromazine. Rapid inhibition of TCA clearance and recovery were found to be major characteristics of CsA. Conversely, the action of chlorpromazine was described by slow onset of inhibition relative to inhibition of TCA clearance by CsA. The established mathematical model, validated by the use of these 2 prototypical cholestatic drugs and the integration of bile canalicular dynamics, provides an important development for the further study of human hepatobiliary function, through simultaneous temporal and vectorial membrane transport of bile acids in drug-induced cholestasis.

Drug-induced cholestasis.

Cholestatic drug‐induced liver injury (DILI) can be a diagnostic challenge due to a large differential diagnosis, variability in clinical presentation, and lack of serologic biomarkers associated with this condition. The clinical presentation of drug‐induced cholestasis includes bland cholestasis, cholestatic hepatitis, secondary sclerosing cholangitis, and vanishing bile duct syndrome. The associate mortality of cholestatic DILI can be as high as 10%, and thus prompt recognition and removal of the offending agent is of critical importance. Several risk factors have been identified for drug‐induced cholestasis, including older age, genetic determinants, and properties of certain medications. Antibiotics, particularly amoxicillin/clavulanate, remain the predominant cause of cholestatic DILI, although a variety of other medications associated with this condition have been identified. In this review, we summarize the presentation, clinical approach, risk factors, implicated medications, and management of drug‐induced cholestatic liver injury. (Hepatology Communications 2017;1:726–735)

Validation of gene expression profiles from cholestatic hepatotoxicants in vitro against human in vivo cholestasis

Drug-induced liver injury remains the most common cause of acute liver failure and a frequently indicated reason for withdrawal of drugs. For the purpose of evaluating the relevance of liver cell models for assessing hepatotoxic risks in intact humans, we here aimed to benchmark 'omics-derived mechanistic data from three in vitro models for parenchymal liver function, intended for the investigation of drug-induced cholestasis, against 'omics data from cholestatic patients. Transcriptomic changes in HepG2 cells, primary mouse hepatocytes and primary human hepatocytes exposed to known cholestatic compounds were analyzed using microarrays. Some of the differentially expressed genes in HepG2 cells were also differentially expressed into the same direction in human cholestasis. The overlap between drug-induced transcriptomic responses in primary mouse hepatocytes and primary human hepatocytes appeared limited and no genes overlapping with in vivo cholestasis were found. Thereupon, a pathway for drug-induced cholestasis was used to map the drug-induced transcriptomic modifications involved in bile salt homeostasis. Indications of an adaptive response to prevent and reduce intracellular bile salt accumulation were observed in vivo as well as in the in vitro liver models. Furthermore, drug-specific changes were found, which may be indicative for their cholestatic properties. Furthermore, connectivity mapping was applied in order to investigate the predictive value of the in vitro models for in vivo cholestasis. This analysis resulted in a positive connection score for most compounds, which may indicate that for identifying cholestatic compounds the focus should be on gene expression signatures rather than on differentially expressed genes.

8-Methoxypsoralen disrupts MDR3-mediated phospholipids efflux and bile acid homeostasis and its relevance to hepatotoxicity.

Since its discovery in 1987, multidrug resistance 3 P-glycoprotein (MDR3) had recognized to play a crucial role in the translocation of phospholipids from the inner to outer leaflets of bile canalicular membranes. An increasing number of reports suggest that drug-mediated functional disruption of MDR3 is responsible for drug-induced cholestasis. 8-Methoxypsoralen (8-MOP) is used clinically to treat psoriasis, vitiligo and other skin disorders. However, psoralens safety for long-term use is a concern. In the current study, we evaluate 8-MOP's potential hepatotoxicity and effects on bile formation. Sprague Dawley (SD) rats were treated daily 200mg/kg or 400mg/kg of 8-MOP orally for 28 days. The result showed a prominent decrease in biliary phospholipids output, which associated with the down-regulation of MDR3. Elevated bile acid serum level and increased biliary bile acid outputs were observed in 8-MOP-treated groups. The disturbance of bile acid homeostasis was associated with changes in enzymes and proteins involved in bile acid synthesis, regulation and transport. Human liver cell line L02 was used to determine on the mRNA and protein levels of MDR3. Cells treated with 8-MOP reveled a decrease in fluorescent PC (phosphatidylcholine) secretion into the pseudocanaliculi (formed between adjacent cells) compared with untreated cells. Our investigation represent the first evidence that 8-MOP can induce cholestatic liver injury by disturbing MDR3-mediated phospholipids efflux and bile acid homeostasis.

Structure based classification for bile salt export pump (BSEP) inhibitors using comparative structural modeling of human BSEP

The bile salt export pump (BSEP) actively transports conjugated monovalent bile acids from the hepatocytes into the bile. This facilitates the formation of micelles and promotes digestion and absorption of dietary fat. Inhibition of BSEP leads to decreased bile flow and accumulation of cytotoxic bile salts in the liver. A number of compounds have been identified to interact with BSEP, which results in drug-induced cholestasis or liver injury. Therefore, in silico approaches for flagging compounds as potential BSEP inhibitors would be of high value in the early stage of the drug discovery pipeline. Up to now, due to the lack of a high-resolution X-ray structure of BSEP, in silico based identification of BSEP inhibitors focused on ligand-based approaches. In this study, we provide a homology model for BSEP, developed using the corrected mouse P-glycoprotein structure (PDB ID: 4M1M). Subsequently, the model was used for docking-based classification of a set of 1212 compounds (405 BSEP inhibitors, 807 non-inhibitors). Using the scoring function ChemScore, a prediction accuracy of 81% on the training set and 73% on two external test sets could be obtained. In addition, the applicability domain of the models was assessed based on Euclidean distance. Further, analysis of the protein–ligand interaction fingerprints revealed certain functional group-amino acid residue interactions that could play a key role for ligand binding. Though ligand-based models, due to their high speed and accuracy, remain the method of choice for classification of BSEP inhibitors, structure-assisted docking models demonstrate reasonably good prediction accuracies while additionally providing information about putative protein–ligand interactions.

The Causes and Clinical Spectrum of Drug-Induced Cholestatic Liver Injury

A variety of drugs, most notably anabolic steroids and some herbal and dietary supplements (HDS), are associated with drug-induced liver injury (DILI) and cholestasis. In the current review, we highlight the recent data on cholestatic DILI and discuss the potential for novel agents as therapeutic options. The clinical phenotype of cholestatic DILI is currently defined based on the “R” ratio calculated from the liver biochemistries at the time of initial presentation. Based on histology, cholestatic DILI may be categorized as acute cholestasis, chronic cholestasis, and cholestatic hepatitis. Histologic injuries include cholangiodestructive features and bile duct loss with the later presentation associated with prolonged recovery period and a poorer prognosis. Current management approaches include early discontinuation of the offending agent, close monitoring, symptom management, selective use of liver biopsy for prognostication, and use of hepatoprotective agents such as ursodeoxycholic acid. Novel treatment agents that are currently under evaluation for the treatment of cholestatic diseases may have a role in the treatment of cholestatic DILI. Cholestatic DILI is increasingly being appreciated, especially in the context of HDS with anabolic steroids. Recent studies report prolonged recovery period and poor prognosis when associated with bile duct loss.

Rat precision-cut liver slices predict drug-induced cholestatic injury

Drug-induced cholestasis (DIC) is one of the leading manifestations of drug-induced liver injury (DILI). As the underlying mechanisms for DIC are not fully known and specific and predictive biomarkers and pre-clinical models are lacking, the occurrence of DIC is often only reported when the drug has been approved for registration. Therefore, appropriate models that predict the cholestatic potential of drug candidates and/or provide insight into the mechanism of DIC are highly needed. We investigated the application of rat precision-cut liver slices (PCLS) to predict DIC, using several biomarkers of cholestasis: hepatocyte viability, intracellular accumulation of total as well as individual bile acids and changes in the expression of genes known to play a role in cholestasis. Rat PCLS exposed to the cholestatic drugs chlorpromazine, cyclosporine A and glibenclamide for 48 h in the presence of a 60 μM physiological bile acid (BA) mix reflected various changes associated with cholestasis, such as decrease in hepatocyte viability, accumulation and changes in the composition of BA and changes in the gene expression of Fxr, Bsep and Ntcp. The toxicity of the drugs was correlated with the accumulation of BA, and especially DCA and CDCA and their conjugates, but to a different extent for different drugs, indicating that BA toxicity is not the only cause for the toxicity of cholestatic drugs. Moreover, our study supports the use of several biomarkers to test drugs for DIC. In conclusion, our results indicate that PCLS may represent a physiological and valuable model to identify cholestatic drugs and provide insight into the mechanisms underlying DIC.

Synthesis, in vitro and in vivo evaluation of 3β-[18F]fluorocholic acid for the detection of drug-induced cholestasis in mice.

IntroductionDrug-induced cholestasis is a liver disorder that might be caused by interference of drugs with the hepatobiliary bile acid transporters. It is important to identify this interference early on in drug development. In this work, Positron Emission Tomography (PET)-imaging with a 18F labeled bile acid analogue was introduced to detect disturbed hepatobiliary transport of bile acids.Methods3β-[18F]fluorocholic acid ([18F]FCA) was prepared by nucleophilic substitution of a mesylated precursor with [18F]fluoride, followed by deprotection with sodium hydroxide. Transport of [18F]FCA was assessed in vitro using CHO-NTCP, HEK-OATP1B1, HEK-OATP1B3 transfected cells and BSEP & MRP2 membrane vesicles. Investigation of [18F]FCA metabolites was performed with primary mouse hepatocytes. Hepatobiliary transport of [18F]FCA was evaluated in vivo in wild-type, rifampicin and bosentan pretreated FVB-mice by dynamic μPET scanning.ResultsRadiosynthesis of [18F]FCA was achieved in a moderate radiochemical yield (8.11 ± 1.94%; non-decay corrected; n = 10) and high radiochemical purity (>99%). FCA was transported by the basolateral bile acid uptake transporters NTCP, OATP1B1 and OATP1B3. For canalicular efflux, BSEP and MRP2 are the relevant bile acid transporters. [18F]FCA was found to be metabolically stable. In vivo, [18F]FCA showed fast hepatic uptake (4.5 ± 0.5 min to reach 71.8 ± 1.2% maximum % ID) and subsequent efflux to the gallbladder and intestines (93.3 ± 6.0% ID after 1 hour). Hepatobiliary transport of [18F]FCA was significantly inhibited by both rifampicin and bosentan.ConclusionA 18F labeled bile acid analogue, [18F]FCA, has been developed that shows transport by NTCP, OATP, MRP2 and BSEP. [18F]FCA can be used as a probe to monitor disturbed hepatobiliary transport in vivo and accumulation of bile acids in blood and liver during drug development.

Predicting Drug-Induced Cholestasis with the Help of Hepatic Transporters-An in Silico Modeling Approach.

Cholestasis represents one out of three types of drug induced liver injury (DILI), which comprises a major challenge in drug development. In this study we applied a two-class classification scheme based on k-nearest neighbors in order to predict cholestasis, using a set of 93 two-dimensional (2D) physicochemical descriptors and predictions of selected hepatic transporters’ inhibition (BSEP, BCRP, P-gp, OATP1B1, and OATP1B3). In order to assess the potential contribution of transporter inhibition, we compared whether the inclusion of the transporters’ inhibition predictions contributes to a significant increase in model performance in comparison to the plain use of the 93 2D physicochemical descriptors. Our findings were in agreement with literature findings, indicating a contribution not only from BSEP inhibition but a rather synergistic effect deriving from the whole set of transporters. The final optimal model was validated via both 10-fold cross validation and external validation. It performs quite satisfactorily resulting in 0.686 ± 0.013 for accuracy and 0.722 ± 0.014 for area under the receiver operating characteristic curve (AUC) for 10-fold cross-validation (mean ± standard deviation from 50 iterations).

Assessing the Risk of Drug-Induced Cholestasis Using Unbound Intrahepatic Concentrations.

Inhibition of the bile salt export pump (BSEP) has been recognized as a key factor in the development of drug-induced cholestasis (DIC). The risk of DIC in humans has been previously assessed using in vitro BSEP inhibition data (IC50) and unbound systemic drug exposure under assumption of the "free drug hypothesis." This concept, however, is unlikely valid, as unbound intrahepatic drug concentrations are affected by active transport and metabolism. To investigate this hypothesis, we experimentally determined the in vitro liver-to-blood partition coefficients (Kpuu) for 18 drug compounds using the hepatic extended clearance model (ECM). In vitro-in vivo translatability of Kpuu values was verified for a subset of compounds in rat. Consequently, unbound intrahepatic concentrations were calculated from clinical exposure (systemic and hepatic inlet) and measured Kpuu data. Using these values, corresponding safety margins against BSEP IC50 values were determined and compared with the clinical incidence of DIC. Depending on the ECM class of a drug, in vitro Kpuu values deviated up to 14-fold from unity, and unbound intrahepatic concentrations were affected accordingly. The use of in vitro Kpuu-based safety margins allowed separation of clinical cholestasis frequency into three classes (no cholestasis, cholestasis in ≤2%, and cholestasis in >2% of subjects) for 17 out of 18 compounds. This assessment was significantly superior compared with using unbound extracellular concentrations as a surrogate for intrahepatic concentrations. Furthermore, the assessment of Kpuu according to ECM provides useful guidance for the quantitative evaluation of genetic and physiologic risk factors for the development of cholestasis.

Nimesulide and 4'-Hydroxynimesulide as Bile Acid Transporters Inhibitors Are Contributory Factors for Drug-Induced Cholestasis.

Nimesulide (NIM) is a classic nonsteroidal anti-inflammatory drug. However, some patients treated with NIM experienced cholestatic liver injury. For this reason, we investigated the potential mechanism underlying NIM-induced cholestasis by using in vivo and in vitro models. Oral administration of 100 mg/kg/day NIM to Wistar rats for 5 days increased the levels of plasma total bile acids, alkaline phosphatase, alanine aminotransferase, and aspartate aminotransferase by 1.49-, 1.31-, 1.60-, and 1.29-fold, respectively. In sandwich-cultured rat hepatocytes, NIM and 4'-hydroxynimesulide (M1) reduced the biliary excretion index of d8-taurocholic acid (d8-TCA) and 5 (and 6)-carboxy-2',7'-dichlorofluorescein in a concentration-dependent manner, indicating the inhibition of the efflux transporters bile salt export pump and multidrug resistance-associated protein 2, respectively. In suspended rat hepatocytes, NIM and M1 inhibited the uptake transporters of d8-TCA for Na+-taurocholate cotransporting polypeptide at IC50 values of 21.3 and 25.0 μM, respectively, and for organic anion-transporting proteins at IC50 values of 45.6 and 39.4 μM, respectively. By contrast, nitro-reduced NIM and the further acetylated metabolite did not inhibit or only marginally inhibited these transporters at the maximum soluble concentrations. Inhibitory effects of NIM and M1 on human bile acid transporters were also confirmed using sandwich-cultured human hepatocytes. These data suggest that the inhibition of bile acid transporters by NIM and M1 is one of the biologic mechanisms of NIM-induced cholestasis.

Evaluation of transcriptomic signature as a valuable tool to study drug-induced cholestasis in primary human hepatocytes.

Primary human hepatocyte (PHH) sandwich cultures from five different donors were daily exposed to cyclosporine A (CsA), ibuprofen (IBU), chlorpromazine (CPZ), amiodarone (AMI) and paracetamol (APAP) at their respective Cmax (total) for short-term (1-3 days) and long-term treatment (14days). Whole genome mRNA profiles (34,693 genes in total) were conducted using an Illumina microarray platform. The impact of compound treatments on gene signatures involved in liver differentiation, cholestasis and in bile acid homeostasis was evaluated. Notably, PHH from the five donors showed a highly comparable phenotype of terminally differentiated hepatocytes. As expected, PHH exposed to 100µM APAP showed no signs of hepatotoxicity both after short- and long-term treatment. CsA at 0.7µM, IBU at 100µM, AMI at 2.5µM and CPZ at 0.1-0.2µM presented, in line with their cholestatic syndromes reported at therapeutic doses, transcriptomic signatures of cholestasis in PHH cultures; deregulation of genes involved in bile acid homeostasis further confirmed this finding. The strength of the cholestasis signature obtained after treatment with CsA, IBU and AMI could be directly related to the basal expression of the respective drug metabolizing enzymes in the various PHH cultures from different individuals. Our data show that the PHH model system combined with transcriptomics carries the future promise to identify individual gene expression profiles predictive of increased cholestasis risk. As the present work suggests possible correlation between mRNA levels of ADME relevant genes and a transcriptomic signature of cholestasis, particular focus on this research question could be the emphasis of additional data collection.

Mechanism-based risk assessment strategy for drug-induced cholestasis using the transcriptional benchmark dose derived by toxicogenomics.

Cholestasis is one of the major causes of drug-induced liver injury (DILI), which can result in withdrawal of approved drugs from the market. Early identification of cholestatic drugs is difficult due to the complex mechanisms involved. In order to develop a strategy for mechanism-based risk assessment of cholestatic drugs, we analyzed gene expression data obtained from the livers of rats that had been orally administered with 12 known cholestatic compounds repeatedly for 28 days at three dose levels. Qualitative analyses were performed using two statistical approaches (hierarchical clustering and principle component analysis), in addition to pathway analysis. The transcriptional benchmark dose (tBMD) and tBMD 95% lower limit (tBMDL) were used for quantitative analyses, which revealed three compound sub-groups that produced different types of differential gene expression; these groups of genes were mainly involved in inflammation, cholesterol biosynthesis, and oxidative stress. Furthermore, the tBMDL values for each test compound were in good agreement with the relevant no observed adverse effect level. These results indicate that our novel strategy for drug safety evaluation using mechanism-based classification and tBMDL would facilitate the application of toxicogenomics for risk assessment of cholestatic DILI.

Integrative "-Omics" Analysis in Primary Human Hepatocytes Unravels Persistent Mechanisms of Cyclosporine A-Induced Cholestasis.

Cyclosporine A (CsA) is an undecapeptide with strong immunosuppressant activities and is used a lot after organ transplantation. Furthermore, it may induce cholestasis in the liver. In general, the drug-induced cholestasis (DIC) pathway includes genes involved in the uptake, synthesis, conjugation, and secretion of bile acids. However, whether CsA-induced changes in the cholestasis pathway in vitro are persistent for repeated dose toxicity has not yet been investigated. To explore this, primary human hepatocytes (PHH) were exposed to a subcytotoxic dose of 30 μM CsA daily for 3 and 5 days. To investigate the persistence of induced changes upon terminating CsA exposure after 5 days, a subset of PHH was subjected to a washout period (WO-period) of 3 days. Multiple -omics analyses, comprising whole genome analysis of DNA methylation, gene expression, and microRNA expression, were performed. The CsA-treatment resulted after 3 and 5 days, respectively, in 476 and 20 differentially methylated genes (DMGs), 1353 and 1481 differentially expressed genes (DEGs), and in 22 and 29 differentially expressed microRNAs (DE-miRs). Cholestasis-related pathways appeared induced during CsA-treatment. Interestingly, 828 persistent DEGs and 6 persistent DE-miRs but no persistent DMGs were found after the WO-period. These persistent DEGs and DE-miRs showed concordance for 22 genes. Furthermore, 29 persistent DEGs changed into the same direction as observed in livers from cholestasis patients. None of those 29 DEGs which among others relate to oxidative stress and lipid metabolism are yet present in the DIC pathway or cholestasis adverse outcome pathway (AOP) thus presenting novel findings. In summary, we have demonstrated for the first time a persistent impact of repeated dose administration of CsA on genes and microRNAs related to DIC in the gold standard human liver in vitro model with PHH.

Lipoprotein-X in cholestatic patients causes xanthomas and promotes foam cell formation in human macrophages

Lipoprotein-X (Lp-X) is an abnormal phospholipid-rich lipoprotein found in patients with cholestatic liver disease. Some patients exhibit skin xanthomas and severe hyperlipidemia. We investigated whether Lp-X induces foam cell formation in human-derived macrophages. To compare the atherogenic properties of Lp-X and modified LDL, we isolated Lp-X from 2 patients who had drug-induced cholestasis and xanthoma striata in the interphalangeal folds. We prepared oxidized LDL and acetylated LDL from healthy volunteers for the positive control experiments. When human monocyte-derived macrophages were incubated with these lipoproteins, the isolated Lp-X induced more prominent lipid accumulation than oxidized LDL or acetylated LDL. One case underwent liver biopsy, with the bile ducts showing marked damage, fulfilling the criteria for vanishing bile duct syndrome. The other case was clinically diagnosed as drug-induced hypersensitivity syndrome. In both cases, Lp-X levels decreased markedly and the xanthomas disappeared completely after the improvement of cholestasis. This study indicates that Lp-X induces foam cell formation in human-derived macrophages. Our findings strongly suggest that persistently elevated Lp-X may cause xanthomas.

A dynamic mathematical model of bile acid clearance in HepaRG cells

A dynamic model based on ordinary differential equations that describes uptake, basolateral and canalicular export of taurocholic acid (TCA) in human HepaRG cells is presented. The highly reproducible inter-assay experimental data were used to reliably estimate model parameters. Primary human hepatocytes were similarly evaluated to establish a mathematical model, but with notably higher inter-assay differences in TCA clearance and bile canaliculi dynamics. By use of the HepaRG cell line, the simultaneous TCA clearance associated to basolateral uptake, canalicular and sinusoidal efflux, was predicted. The mathematical model accurately reproduced the dose-dependent inhibition of TCA clearance in the presence and absence of the prototypical cholestatic drugs cyclosporine A (CsA) and chlorpromazine. Rapid inhibition of TCA clearance and recovery were found to be major characteristics of CsA. Conversely, the action of chlorpromazine was described by slow onset of inhibition relative to inhibition of TCA clearance by CsA. The established mathematical model, validated by the use of these 2 prototypical cholestatic drugs and the integration of bile canalicular dynamics, provides an important development for the further study of human hepatobiliary function, through simultaneous temporal and vectorial membrane transport of bile acids in drug-induced cholestasis.