Population pharmacokinetic model to generate mechanistic insights in bile acid homeostasis and drug-induced cholestasis.

Bile acids (BA) fulfill a wide range of physiological functions, but are also involved in pathologies, such as cholestasis. Cholestasis is characterized by an intrahepatic accumulation of BAs and subsequent spillage to the systemic circulation. The aim of the present study was to develop physiologically based kinetic (PBK) models that would provide a tool to predict dose-dependent BA accumulation in humans upon treatment with a Bile Salt Export Pump (BSEP) inhibitor. We developed a PBK model describing the BA homeostasis using glycochenodeoxycholic acid as an exemplary BA. Population wide distributions of BSEP abundances were incorporated in the PBK model using Markov Chain Monte Carlo simulations, and alternatively the total amount of BAs was scaled empirically to describe interindividual differences in plasma BA levels. Next, the effects of the BSEP inhibitor bosentan on the BA levels were simulated. The PBK model developed adequately predicted the in vivo BA dynamics. Both the Markov Chain Monte Carlo simulations based on a distribution of BSEP abundances and empirical scaling of the total BA pool readily described the variations within and between data in human volunteers. Bosentan treatment disproportionally increased the maximum BA concentration in individuals with a large total BA pool or low BSEP abundance. Especially individuals having a large total BA pool size and a low BSEP abundance were predicted to be at risk for rapid saturation of BSEP and subsequent intrahepatic BA accumulation. This model provides a first estimate of personalized safe therapeutic external dose levels of compounds with BSEP-inhibitory properties.

Omics-based responses induced by bosentan in human hepatoma HepaRG cell cultures.

Bosentan is well known to induce cholestatic liver toxicity in humans. The present study was set up to characterize the hepatotoxic effects of this drug at the transcriptomic, proteomic, and metabolomic levels. For this purpose, human hepatoma-derived HepaRG cells were exposed to a number of concentrations of bosentan during different periods of time. Bosentan was found to functionally and transcriptionally suppress the bile salt export pump as well as to alter bile acid levels. Pathway analysis of both transcriptomics and proteomics data identified cholestasis as a major toxicological event. Transcriptomics results further showed several gene changes related to the activation of the nuclear farnesoid X receptor. Induction of oxidative stress and inflammation were also observed. Metabolomics analysis indicated changes in the abundance of specific endogenous metabolites related to mitochondrial impairment. The outcome of this study may assist in the further optimization of adverse outcome pathway constructs that mechanistically describe the processes involved in cholestatic liver injury.