Scott Q. Siler

In vitro to in vivo extrapolation and species response comparisons for drug-induced liver injury (DILI) using DILIsym™: a mechanistic, mathematical model of DILI

Drug-induced liver injury (DILI) is not only a major concern for all patients requiring drug therapy, but also for the pharmaceutical industry. Many new in vitro assays and pre-clinical animal models are being developed to help screen compounds for the potential to cause DILI. This study demonstrates that mechanistic, mathematical modeling offers a method for interpreting and extrapolating results. The DILIsym™ model (version 1A), a mathematical representation of DILI, was combined with in vitro data for the model hepatotoxicant methapyrilene (MP) to carry out an in vitro to in vivo extrapolation. In addition, simulations comparing DILI responses across species illustrated how modeling can aid in selecting the most appropriate pre-clinical species for safety testing results relevant to humans. The parameter inputs used to predict DILI for MP were restricted to in vitro inputs solely related to ADME (absorption, distribution, metabolism, elimination) processes. MP toxicity was correctly predicted to occur in rats, but was not apparent in the simulations for humans and mice (consistent with literature). When the hepatotoxicity of MP and acetaminophen (APAP) was compared across rats, mice, and humans at an equivalent dose, the species most susceptible to APAP was not susceptible to MP, and vice versa. Furthermore, consideration of variability in simulated population samples (SimPops™) provided confidence in the predictions and allowed examination of the biological parameters most predictive of outcome. Differences in model sensitivity to the parameters were related to species differences, but the severity of DILI for each drug/species combination was also an important factor.

Modeling Drug- and Chemical-Induced Hepatotoxicity with Systems Biology Approaches

We provide an overview of computational systems biology approaches as applied to the study of chemical- and drug-induced toxicity. The concept of “toxicity pathways” is described in the context of the 2007 US National Academies of Science report, “Toxicity testing in the 21st Century: A Vision and A Strategy.” Pathway mapping and modeling based on network biology concepts are a key component of the vision laid out in this report for a more biologically based analysis of dose-response behavior and the safety of chemicals and drugs. We focus on toxicity of the liver (hepatotoxicity) – a complex phenotypic response with contributions from a number of different cell types and biological processes. We describe three case studies of complementary multi-scale computational modeling approaches to understand perturbation of toxicity pathways in the human liver as a result of exposure to environmental contaminants and specific drugs. One approach involves development of a spatial, multicellular “virtual tissue” model of the liver lobule that combines molecular circuits in individual hepatocytes with cell–cell interactions and blood-mediated transport of toxicants through hepatic sinusoids, to enable quantitative, mechanistic prediction of hepatic dose-response for activation of the aryl hydrocarbon receptor toxicity pathway. Simultaneously, methods are being developing to extract quantitative maps of intracellular signaling and transcriptional regulatory networks perturbed by environmental contaminants, using a combination of gene expression and genome-wide protein-DNA interaction data. A predictive physiological model (DILIsym™) to understand drug-induced liver injury (DILI), the most common adverse event leading to termination of clinical development programs and regulatory actions on drugs, is also described. The model initially focuses on reactive metabolite-induced DILI in response to administration of acetaminophen, and spans multiple biological scales.

Exploring BSEP inhibition-mediated toxicity with a mechanistic model of drug-induced liver injury

Inhibition of the bile salt export pump (BSEP) has been linked to incidence of drug-induced liver injury (DILI), presumably by the accumulation of toxic bile acids in the liver. We have previously constructed and validated a model of bile acid disposition within DILIsym®, a mechanistic model of DILI. In this paper, we use DILIsym® to simulate the DILI response of the hepatotoxic BSEP inhibitors bosentan and CP-724,714 and the non-hepatotoxic BSEP inhibitor telmisartan in humans in order to explore whether we can predict that hepatotoxic BSEP inhibitors can cause bile acid accumulation to reach toxic levels. We also simulate bosentan in rats in order to illuminate potential reasons behind the lack of toxicity in rats compared to the toxicity observed in humans. DILIsym® predicts that bosentan, but not telmisartan, will cause mild hepatocellular ATP decline and serum ALT elevation in a simulated population of humans. The difference in hepatotoxic potential between bosentan and telmisartan is consistent with clinical observations. However, DILIsym® underpredicts the incidence of bosentan toxicity. DILIsym® also predicts that bosentan will not cause toxicity in a simulated population of rats, and that the difference between the response to bosentan in rats and in humans is primarily due to the less toxic bile acid pool in rats. Our simulations also suggest a potential synergistic role for bile acid accumulation and mitochondrial electron transport chain (ETC) inhibition in producing the observed toxicity in CP-724,714, and suggest that CP-724,714 metabolites may also play a role in the observed toxicity. Our work also compares the impact of competitive and noncompetitive BSEP inhibition for CP-724,714 and demonstrates that noncompetitive inhibition leads to much greater bile acid accumulation and potential toxicity. Our research demonstrates the potential for mechanistic modeling to contribute to the understanding of how bile acid transport inhibitors cause DILI.

Linking physiology to toxicity using DILIsym®, a mechanistic mathematical model of drug‐induced liver injury

The drug development industry faces multiple challenges in the realization of safe effective drugs. Computational modeling approaches can be used to support these efforts. One approach, mechanistic modeling, is new to the realm of drug safety. It holds the promise of not only predicting toxicity for novel compounds, but also illuminating the mechanistic underpinnings of toxicity. To increase the scientific communitys familiarity with mechanistic modeling in drug safety, this article seeks to provide perspective on the type of data used, how they are used and where they are lacking. Examples are derived from the development of DILIsym(®) software, a mechanistic model of drug-induced liver injury (DILI). DILIsym(®) simulates the mechanistic interactions and events from compound administration through the progression of liver injury and regeneration. Modeling mitochondrial toxicity illustrates the type and use of in vitro data to represent biological interactions, as well as insights on key differences between in vitro and in vivo conditions. Modeling bile acid toxicity illustrates a case in which the over-arching mechanism is well accepted, but many mechanistic details are lacking. Modeling was used to identify measurements predicted to strongly impact toxicity. Finally, modeling innate immune responses illustrates the importance of time-series data, particularly in the presence of positive and negative feedback loops, as well as the need for data from different animal species for better translation. These concepts are germane to most mechanistic models, although the details will vary. The use of mechanistic models is expected to improve the rational design of new drugs.

Mechanistic Modeling Reveals the Critical Knowledge Gaps in Bile Acid–Mediated DILI

Bile salt export pump (BSEP) inhibition has been proposed to be an important mechanism for drug‐induced liver injury (DILI). Modeling can prioritize knowledge gaps concerning bile acid (BA) homeostasis and thus help guide experimentation. A submodel of BA homeostasis in rats and humans was constructed within DILIsym, a mechanistic model of DILI. In vivo experiments in rats with glibenclamide were conducted, and data from these experiments were used to validate the model. The behavior of DILIsym was analyzed in the presence of a simulated theoretical BSEP inhibitor. BSEP inhibition in humans is predicted to increase liver concentrations of conjugated chenodeoxycholic acid (CDCA) and sulfate‐conjugated lithocholic acid (LCA) while the concentration of other liver BAs remains constant or decreases. On the basis of a sensitivity analysis, the most important unknowns are the level of BSEP expression, the amount of intestinal synthesis of LCA, and the magnitude of farnesoid‐X nuclear receptor (FXR)‐mediated regulation.

A Mechanistic Model of Drug‐Induced Liver Injury Aids the Interpretation of Elevated Liver Transaminase Levels in a Phase I Clinical Trial

Entolimod (CBLB502) is a Toll‐like receptor 5 agonist in development as a single‐dose countermeasure against total body irradiation. Efficacy can be assessed from animal studies, but the “Animal Rule” does not apply to safety assessment. Marked elevations of serum aminotransferases (exceeding 1,000 IU/l) were observed in some human subjects receiving Entolimod in a safety study, threatening its continued development. The percentage of total hepatocytes undergoing necrosis in these subjects was estimated using a mechanistic, multiscale, mathematical model (DILIsym). The simulations suggested that no subject in the safety study experienced more than a modest loss of hepatocytes (<5%), which was comparable to estimates from a study of healthy volunteers receiving treatment with heparins. The predicted hepatocyte loss with Entolimod was lower than that required to cause liver dysfunction or that is routinely excised from volunteers donating for autologous liver transplantation and did not likely represent a serious health risk.

Application of a Mechanistic Model to Evaluate Putative Mechanisms of Tolvaptan Drug-Induced Liver Injury and Identify Patient Susceptibility Factors

Tolvaptan is a selective vasopressin V2 receptor antagonist, approved in several countries for the treatment of hyponatremia and autosomal dominant polycystic kidney disease (ADPKD). No liver injury has been observed with tolvaptan treatment in healthy subjects and in non-ADPKD indications, but ADPKD clinical trials showed evidence of drug-induced liver injury (DILI). Although all DILI events resolved, additional monitoring in tolvaptan-treated ADPKD patients is required. In vitro assays identified alterations in bile acid disposition and inhibition of mitochondrial respiration as potential mechanisms underlying tolvaptan hepatotoxicity. This report details the application of DILIsym software to determine whether these mechanisms could account for the liver safety profile of tolvaptan observed in ADPKD clinical trials. DILIsym simulations included physiologically based pharmacokinetic estimates of hepatic exposure for tolvaptan and2 metabolites, and their effects on hepatocyte bile acid transporters and mitochondrial respiration. The frequency of predicted alanine aminotransferase (ALT) elevations, following simulated 90/30  mg split daily dosing, was 7.9% compared with clinical observations of 4.4% in ADPKD trials. Toxicity was multifactorial as inhibition of bile acid transporters and mitochondrial respiration contributed to the simulated DILI. Furthermore, simulation analysis identified both pre-treatment risk factors and on-treatment biomarkers predictive of simulated DILI. The simulations demonstrated that in vivo hepatic exposure to tolvaptan and the DM-4103 metabolite, combined with these 2 mechanisms of toxicity, were sufficient to account for the initiation of tolvaptan-mediated DILI. Identification of putative risk-factors and potential novel biomarkers provided insight for the development of mechanism-based tolvaptan risk-mitigation strategies.

Sandwich-Cultured Hepatocytes as a Tool to Study Drug Disposition and Drug-Induced Liver Injury.

Sandwich-cultured hepatocytes (SCH) are metabolically competent and have proper localization of basolateral and canalicular transporters with functional bile networks. Therefore, this cellular model is a unique tool that can be used to estimate biliary excretion of compounds. SCH have been used widely to assess hepatobiliary disposition of endogenous and exogenous compounds and metabolites. Mechanistic modeling based on SCH data enables estimation of metabolic and transporter-mediated clearances, which can be used to construct physiologically based pharmacokinetic models for prediction of drug disposition and drug-drug interactions in humans. In addition to pharmacokinetic studies, SCH also have been used to study cytotoxicity and perturbation of biological processes by drugs and hepatically generated metabolites. Human SCH can provide mechanistic insights underlying clinical drug-induced liver injury (DILI). In addition, data generated in SCH can be integrated into systems pharmacology models to predict potential DILI in humans. In this review, applications of SCH in studying hepatobiliary drug disposition and bile acid-mediated DILI are discussed. An example is presented to show how data generated in the SCH model were used to establish a quantitative relationship between intracellular bile acids and cytotoxicity, and how this information was incorporated into a systems pharmacology model for DILI prediction.

Systems pharmacology modeling of drug-induced hyperbilirubinemia: Differentiating hepatotoxicity and inhibition of enzymes/transporters

Elevations in serum bilirubin during drug treatment may indicate global liver dysfunction and a high risk of liver failure. However, drugs also can increase serum bilirubin in the absence of hepatic injury by inhibiting specific enzymes/transporters. We constructed a mechanistic model of bilirubin disposition based on known functional polymorphisms in bilirubin metabolism/transport. Using physiologically based pharmacokinetic (PBPK) model‐predicted drug exposure and enzyme/transporter inhibition constants determined in vitro, our model correctly predicted indinavir‐mediated hyperbilirubinemia in humans and rats. Nelfinavir was predicted not to cause hyperbilirubinemia, consistent with clinical observations. We next examined a new drug candidate that caused both elevations in serum bilirubin and biochemical evidence of liver injury in rats. Simulations suggest that bilirubin elevation primarily resulted from inhibition of transporters rather than global liver dysfunction. We conclude that mechanistic modeling of bilirubin can help elucidate underlying mechanisms of drug‐induced hyperbilirubinemia, and thereby distinguish benign from clinically important elevations in serum bilirubin.

Refining Liver Safety Risk Assessment: Application of Mechanistic Modeling and Serum Biomarkers to Cimaglermin Alfa (GGF2) Clinical Trials

Cimaglermin alfa (GGF2) is a recombinant human protein growth factor in development for heart failure. Phase I trials were suspended when two cimaglermin alfa‐treated subjects experienced concomitant elevations in serum aminotransferases and total bilirubin, meeting current US Food and Drug Administration criteria for a serious liver safety signal (i.e., “Hys Law”). We assayed mechanistic biomarkers in archived clinical trial serum samples which confirmed the hepatic origin of the aminotransferase elevations in these two subjects and identified apoptosis as the major mode of hepatocyte death. Using a mathematical model of drug‐induced liver injury (DILIsym) and a simulated population, we estimated that the maximum hepatocyte loss in these two subjects was <13%, which would not result in liver dysfunction sufficient to significantly increase serum bilirubin levels. We conclude that the two subjects should not be considered Hys Law cases and that mechanistic biomarkers and modeling can aid in refining liver safety risk assessment in clinical trials.