Mechanistic Modelling of Drug-Induced Liver Injury: Investigating the Role of Innate Immune Responses

Modeling and Simulation of Acetaminophen Pharmacokinetics and Hepatic Biomarkers After Overdoses of Extended-Release and Immediate-Release Formulations in Healthy Adults Using the Quantitative Systems Toxicology Software Platform DILIsym

Acetaminophen (APAP) has been formulated as immediate-, modified-, and extended-release tablets (APAP-IR, -MR, and -ER, respectively). However, there was concern that APAP-MR previously available in Europe could form a bezoar after a large overdose, leading to delayed absorption and atypical pharmacokinetics (PK) compared to APAP-IR, and that current treatment guidelines developed for APAP overdose to prevent severe hepatotoxicity are inappropriate for APAP-MR. In contrast, APAP-ER caplets available in the United States are designed with an IR layer and an erodible ER layer. Using modeling and simulation, predicted PK and hepatotoxicity biomarkers following various acute overdose and repeated supratherapeutic ingestion (RSTI) scenarios with APAP-IR and APAP-ER were compared to investigate the differences between these two formulations. The existing APAP-IR representation within DILIsym v8A, a quantitative systems toxicology model of drug-induced liver injury, was updated, and an APAP-ER model was developed, using newly acquired in vitro (e.g., tiny-TIMsg) and clinical data. The model and simulated populations (SimPops) representing healthy adults were extensively validated, before simulating PK and three clinically useful hepatic biomarkers after various overdose scenarios. On average, APAP exposure after acute overdose and RSTI in healthy adults was predicted to be slightly lower for APAP-ER compared to APAP-IR, partially due to lower APAP absorption for APAP-ER, while not markedly impacting the expected time course of APAP plasma concentrations. Similar hepatic biomarker profiles were predicted for both APAP formulations. Based on these results, the APAP overdose consensus treatment guidelines updated in 2023 are not further impacted by this report.

Critical Roles of the Sphingolipid Metabolic Pathway in Liver Regeneration, Hepatocellular Carcinoma Progression and Therapy

The sphingolipid metabolic pathway, an important signaling pathway, plays a crucial role in various physiological processes including cell proliferation, survival, apoptosis, and immune regulation. The liver has the unique ability to regenerate using bioactive lipid mediators involving multiple sphingolipids, including ceramide and sphingosine 1-phosphate (S1P). Dysregulation of the balance between sphingomyelin, ceramide, and S1P has been implicated in the regulation of liver regeneration and diseases, including liver fibrosis and hepatocellular carcinoma (HCC). Understanding and modulating this balance may have therapeutic implications for tumor proliferation, progression, and metastasis in HCC. For cancer therapy, several inhibitors and activators of sphingolipid signaling, including ABC294640, SKI-II, and FTY720, have been discussed. Here, we elucidate the critical roles of the sphingolipid pathway in the regulation of liver regeneration, fibrosis, and HCC. Regulation of sphingolipids and their corresponding enzymes may considerably influence new insights into therapies for various liver disorders and diseases.

Investigating bile acid-mediated cholestatic drug-induced liver injury using a mechanistic model of multidrug resistance protein 3 (MDR3) inhibition

Inhibition of the canalicular phospholipid floppase multidrug resistance protein 3 (MDR3) has been implicated in cholestatic drug-induced liver injury (DILI), which is clinically characterized by disrupted bile flow and damage to the biliary epithelium. Reduction in phospholipid excretion, as a consequence of MDR3 inhibition, decreases the formation of mixed micelles consisting of bile acids and phospholipids in the bile duct, resulting in a surplus of free bile acids that can damage the bile duct epithelial cells, i.e., cholangiocytes. Cholangiocytes may compensate for biliary increases in bile acid monomers via the cholehepatic shunt pathway or bicarbonate secretion, thereby influencing viability or progression to toxicity. To address the unmet need to predict drug-induced bile duct injury in humans, DILIsym, a quantitative systems toxicology model of DILI, was extended by representing key features of the bile duct, cholangiocyte functionality, bile acid and phospholipid disposition, and cholestatic hepatotoxicity. A virtual, healthy representative subject and population (n = 285) were calibrated and validated utilizing a variety of clinical data. Sensitivity analyses were performed for 1) the cholehepatic shunt pathway, 2) biliary bicarbonate concentrations and 3) modes of MDR3 inhibition. Simulations showed that an increase in shunting may decrease the biliary bile acid burden, but raise the hepatocellular concentrations of bile acids. Elevating the biliary concentration of bicarbonate may decrease bile acid shunting, but increase bile flow rate. In contrast to competitive inhibition, simulations demonstrated that non-competitive and mixed inhibition of MDR3 had a profound impact on phospholipid efflux, elevations in the biliary bile acid-to-phospholipid ratio, cholangiocyte toxicity, and adaptation pathways. The model with its extended bile acid homeostasis representation was furthermore able to predict DILI liability for compounds with previously studied interactions with bile acid transport. The cholestatic liver injury submodel in DILIsym accounts for several processes pertinent to bile duct viability and toxicity and hence, is useful for predictions of MDR3 inhibition-mediated cholestatic DILI in humans.

Applications of In Silico Models to Predict Drug-Induced Liver Injury

Drug-induced liver injury (DILI) is a major cause of the withdrawal of pre-marketed drugs, typically attributed to oxidative stress, mitochondrial damage, disrupted bile acid homeostasis, and innate immune-related inflammation. DILI can be divided into intrinsic and idiosyncratic DILI with cholestatic liver injury as an important manifestation. The diagnosis of DILI remains a challenge today and relies on clinical judgment and knowledge of the insulting agent. Early prediction of hepatotoxicity is an important but still unfulfilled component of drug development. In response, in silico modeling has shown good potential to fill the missing puzzle. Computer algorithms, with machine learning and artificial intelligence as a representative, can be established to initiate a reaction on the given condition to predict DILI. DILIsym is a mechanistic approach that integrates physiologically based pharmacokinetic modeling with the mechanisms of hepatoxicity and has gained increasing popularity for DILI prediction. This article reviews existing in silico approaches utilized to predict DILI risks in clinical medication and provides an overview of the underlying principles and related practical applications.

The immunological mechanisms and therapeutic potential in drug-induced liver injury: lessons learned from acetaminophen hepatotoxicity

Acute liver failure caused by drug overdose is a significant clinical problem in developed countries. Acetaminophen (APAP), a widely used analgesic and antipyretic drug, but its overdose can cause acute liver failure. In addition to APAP-induced direct hepatotoxicity, the intracellular signaling mechanisms of APAP-induced liver injury (AILI) including metabolic activation, mitochondrial oxidant stress and proinflammatory response further affect progression and severity of AILI. Liver inflammation is a result of multiple interactions of cell death molecules, immune cell-derived cytokines and chemokines, as well as damaged cell-released signals which orchestrate hepatic immune cell infiltration. The immunoregulatory interplay of these inflammatory mediators and switching of immune responses during AILI lead to different fate of liver pathology. Thus, better understanding the complex interplay of immune cell subsets in experimental models and defining their functional involvement in disease progression are essential to identify novel therapeutic targets for the treatment of AILI. Here, this present review aims to systematically elaborate on the underlying immunological mechanisms of AILI, its relevance to immune cells and their effector molecules, and briefly discuss great therapeutic potential based on inflammatory mediators.

Quantitative Systems Toxicology and Drug Development: The DILIsym Experience

DILIsym^® is a Quantitative Systems Toxicology (QST) model that has been developed over the last decade by a public-private partnership to predict the liver safety liability in new drug candidates. DILIsym integrates the quantitative abilities of parent and relevant metabolites to cause oxidative stress, mitochondrial dysfunction, and alter bile acid homeostasis. Like the prediction of drug-drug interactions, the data entered into DILIsym are assessed in the laboratory in human experimental systems, and combined with estimates of liver exposure to predict the outcome. DILIsym is now frequently used in decision-making within the pharmaceutical industry and its modeling results are increasingly included in regulatory communications and NDA submissions. DILIsym can be used to identify dominant mechanisms underlying liver toxicity and this information is increasingly being used to identify patient-specific risk factors, including certain disease states. DILIsym is also increasingly used to optimize the interpretation of liver injury biomarkers. DILIsym provides an example of how QST modeling can help speed the delivery of safer new drugs to the patients who need them.

A computational model of liver tissue damage and repair

Drug induced liver injury (DILI) and cell death can result from oxidative stress in hepatocytes. An initial pattern of centrilobular damage in the APAP model of DILI is amplified by communication from stressed cells and immune system activation. While hepatocyte proliferation counters cell loss, high doses are still lethal to the tissue. To understand the progression of disease from the initial damage to tissue recovery or death, we computationally model the competing biological processes of hepatocyte proliferation, necrosis and injury propagation. We parametrize timescales of proliferation (α), conversion of healthy to stressed cells (β) and further sensitization of stressed cells towards necrotic pathways (γ) and model them on a Cellular Automaton (CA) based grid of lattice sites. 1D simulations show that a small α/β (fast proliferation), combined with a large γ/β (slow death) have the lowest probabilities of tissue survival. At large α/β, tissue fate can be described by a critical γ/β* ratio alone; this value is dependent on the initial amount of damage and proportional to the tissue size N. Additionally, the 1D model predicts a minimum healthy population size below which damage is irreversible. Finally, we compare 1D and 2D phase spaces and discuss outcomes of bistability where either survival or death is possible, and of coexistence where simulated tissue never completely recovers or dies but persists as a mixture of healthy, stressed and necrotic cells. In conclusion, our model sheds light on the evolution of tissue damage or recovery and predicts potential for divergent fates given different rates of proliferation, necrosis, and injury propagation.

Mechanistic Investigations Support Liver Safety of Ubrogepant

Small-molecule calcitonin gene-related peptide (CGRP) receptor antagonists have demonstrated therapeutic efficacy for the treatment of migraine. However, previously investigated CGRP receptor antagonists, telcagepant and MK-3207, were discontinued during clinical development because of concerns about drug-induced liver injury. A subsequent effort to identify novel CGRP receptor antagonists less likely to cause hepatotoxicity led to the development of ubrogepant. The selection of ubrogepant, following a series of mechanistic studies conducted with MK-3207 and telcagepant, was focused on key structural modifications suggesting that ubrogepant was less prone to forming reactive metabolites than previous compounds. The potential for each drug to cause liver toxicity was subsequently assessed using a quantitative systems toxicology approach (DILIsym) that incorporates quantitative assessments of mitochondrial dysfunction, disruption of bile acid homeostasis, and oxidative stress, along with estimates of dose-dependent drug exposure to and within liver cells. DILIsym successfully modeled liver toxicity for telcagepant and MK-3207 at the dosing regimens used in clinical trials. In contrast, DILIsym predicted no hepatotoxicity during treatment with ubrogepant, even at daily doses up to 1000 mg (10-fold higher than the approved clinical dose of 100 mg). These predictions are consistent with clinical trial experience showing that ubrogepant has lower potential to cause hepatotoxicity than has been observed with telcagepant and MK-3207.

The DILI-sim Initiative: Insights into Hepatotoxicity Mechanisms and Biomarker Interpretation

The drug‐induced liver injury (DILI)‐sim Initiative is a public‐private partnership involving scientists from industry, academia, and the US Food and Drug Administration (FDA). The Initiative uses quantitative systems toxicology (QST) to build and refine a model (DILIsym) capable of understanding and predicting liver safety liabilities in new drug candidates and to optimize interpretation of liver safety biomarkers used in clinical studies. Insights gained to date include the observation that most dose‐dependent hepatotoxicity can be accounted for by combinations of just three mechanisms (oxidative stress, interference with mitochondrial respiration, and alterations in bile acid homeostasis) and the importance of noncompetitive inhibition of bile acid transporters. The effort has also provided novel insight into species and interpatient differences in susceptibility, structure‐activity relationships, and the role of nonimmune mechanisms in delayed idiosyncratic hepatotoxicity. The model is increasingly used to evaluate new drug candidates and several clinical trials are underway that will test the model's ability to prospectively predict liver safety. With more refinement, in the future, it may be possible to use the DILIsym predictions to justify reduction in the size of some clinical trials. The mature model could also potentially assist physicians in managing the liver safety of their patients as well as aid in the diagnosis of DILI.

Assessing Effects of BHV-0223 40 mg Zydis Sublingual Formulation and Riluzole 50 mg Oral Tablet on Liver Function Test Parameters Utilizing DILIsym

For patients with amyotrophic lateral sclerosis who take oral riluzole tablets, approximately 50% experience alanine transaminase (ALT) levels above upper limit of normal (ULN), 8% above 3× ULN, and 2% above 5× ULN. BHV-0223 is a novel 40 mg rapidly sublingually disintegrating (Zydis) formulation of riluzole, bioequivalent to conventional riluzole 50 mg oral tablets, that averts the need for swallowing tablets and mitigates first-pass hepatic metabolism, thereby potentially reducing risk of liver toxicity. DILIsym is a validated multiscale computational model that supports evaluation of liver toxicity risks. DILIsym was used to compare the hepatotoxicity potential of oral riluzole tablets (50 mg BID) versus BHV-0223 (40 mg BID) by integrating clinical data and in vitro toxicity data. In a simulated population (SimPops), ALT levels > 3× ULN were predicted in 3.9% (11/285) versus 1.4% (4/285) of individuals with oral riluzole tablets and sublingual BHV-0223, respectively. This represents a relative risk reduction of 64% associated with BHV-0223 versus conventional riluzole tablets. Mechanistic investigations revealed that oxidative stress was responsible for the predicted ALT elevations. The validity of the DILIsym representation of riluzole and assumptions is supported by its ability to predict rates of ALT elevations for riluzole oral tablets comparable with that observed in clinical data. Combining a mechanistic, quantitative representation of hepatotoxicity with interindividual variability in both susceptibility and liver exposure suggests that sublingual BHV-0223 confers diminished rates of liver toxicity compared with oral tablets of riluzole, consistent with having a lower overall dose of riluzole and bypassing first-pass liver metabolism.

Understanding Idiosyncratic Toxicity: Lessons Learned from Drug-Induced Liver Injury

Idiosyncratic adverse drug reactions (IADRs) encompass a diverse group of toxicities that can vary by drug and patient. The complex and unpredictable nature of IADRs combined with the fact that they are rare makes them particularly difficult to predict, diagnose, and treat. Common clinical characteristics, the identification of human leukocyte antigen risk alleles, and drug-induced proliferation of lymphocytes isolated from patients support a role for the adaptive immune system in the pathogenesis of IADRs. Significant evidence also suggests a requirement for direct, drug-induced stress, neoantigen formation, and stimulation of an innate response, which can be influenced by properties intrinsic to both the drug and the patient. This Perspective will provide an overview of the clinical profile, mechanisms, and risk factors underlying IADRs as well as new approaches to study these reactions, focusing on idiosyncratic drug-induced liver injury.

Prediction of Safety Margin and Optimization of Dosing Protocol for a Novel Antibiotic using Quantitative Systems Pharmacology Modeling

Elevations of liver enzymes have been observed in clinical trials with BAL30072, a novel antibiotic. In vitro assays have identified potential mechanisms for the observed hepatotoxicity, including electron transport chain (ETC) inhibition and reactive oxygen species (ROS) generation. DILIsym, a quantitative systems pharmacology (QSP) model of drug-induced liver injury, has been used to predict the likelihood that each mechanism explains the observed toxicity. DILIsym was also used to predict the safety margin for a novel BAL30072 dosing scheme; it was predicted to be low. DILIsym was then used to recommend potential modifications to this dosing scheme; weight-adjusted dosing and a requirement to assay plasma alanine aminotransferase (ALT) daily and stop dosing as soon as ALT increases were observed improved the predicted safety margin of BAL30072 and decreased the predicted likelihood of severe injury. This research demonstrates a potential application for QSP modeling in improving the safety profile of candidate drugs.

Systems Toxicology Approach to Identifying Paracetamol Overdose

Paracetamol (acetaminophen (APAP)) is one of the most commonly used analgesics in the United Kingdom and the United States. However, exceeding the maximum recommended dose can cause serious liver injury and even death. Promising APAP toxicity biomarkers are thought to add value to those used currently and clarification of the functional relationships between these biomarkers and liver injury would aid clinical implementation of an improved APAP toxicity identification framework. The framework currently used to define an APAP overdose is highly dependent upon time since ingestion and initial dose; information that is often highly unpredictable. A pharmacokinetic/pharmacodynamic (PK/PD) APAP model has been built in order to understand the relationships between a panel of biomarkers and APAP dose. Visualization and statistical tools have been used to predict initial APAP dose and time since administration. Additionally, logistic regression analysis has been applied to histology data to provide a prediction of the probability of liver injury.

Advancing Predictive Hepatotoxicity at the Intersection of Experimental, in Silico, and Artificial Intelligence Technologies

Adverse drug reactions, particularly those that result in drug-induced liver injury (DILI), are a major cause of drug failure in clinical trials and drug withdrawals. Hepatotoxicity-mediated drug attrition occurs despite substantial investments of time and money in developing cellular assays, animal models, and computational models to predict its occurrence in humans. Underperformance in predicting hepatotoxicity associated with drugs and drug candidates has been attributed to existing gaps in our understanding of the mechanisms involved in driving hepatic injury after these compounds perfuse and are metabolized by the liver. Herein we assess in vitro, in vivo (animal), and in silico strategies used to develop predictive DILI models. We address the effectiveness of several two- and three-dimensional in vitro cellular methods that are frequently employed in hepatotoxicity screens and how they can be used to predict DILI in humans. We also explore how humanized animal models can recapitulate human drug metabolic profiles and associated liver injury. Finally, we highlight the maturation of computational methods for predicting hepatotoxicity, the untapped potential of artificial intelligence for improving in silico DILI screens, and how knowledge acquired from these predictions can shape the refinement of experimental methods.