Virtual Experiments Enable Exploring and Challenging Explanatory Mechanisms of Immune-Mediated P450 Down-Regulation

Utilizing virtual experiments to increase understanding of discrepancies involving in vitro-to-in vivo predictions of hepatic clearance

Predictions of xenobiotic hepatic clearance in humans using in vitro-to-in vivo extrapolation methods are frequently inaccurate and problematic. Multiple strategies are being pursued to disentangle responsible mechanisms. The objective of this work is to evaluate the feasibility of using insights gained from independent virtual experiments on two model systems to begin unraveling responsible mechanisms. The virtual culture is a software analog of hepatocytes in vitro, and the virtual human maps to hepatocytes within a liver within an idealized model human. Mobile objects (virtual compounds) map to amounts of xenobiotics. Earlier versions of the two systems achieved quantitative validation targets for intrinsic clearance (virtual culture) and hepatic clearance (virtual human). The major difference between the two systems is the spatial organization of the virtual hepatocytes. For each pair of experiments (virtual culture, virtual human), hepatocytes are configured the same. Probabilistic rules govern virtual compound movements and interactions with other objects. We focus on highly permeable virtual compounds and fix their extracellular unbound fraction at one of seven values (0.05–1.0). Hepatocytes contain objects that can bind and remove compounds, analogous to metabolism. We require that, for a subset of compound properties, per-hepatocyte compound exposure and removal rates during culture experiments directly predict corresponding measures made during virtual human experiments. That requirement serves as a cross-system validation target; we identify compound properties that enable achieving it. We then change compound properties, ceteris paribus, and provide model mechanism-based explanations for when and why measures made during culture experiments under- (or over-) predict corresponding measures made during virtual human experiments. The results show that, from the perspective of compound removal, the organization of hepatocytes within virtual livers is more efficient than within cultures, and the greater the efficiency difference, the larger the underprediction. That relationship is noteworthy because most in vitro-to-in vivo extrapolation methods abstract away the structural organization of hepatocytes within a liver. More work is needed on multiple fronts, including the study of an expanded variety of virtual compound properties. Nevertheless, the results support the feasibility of the approach and plan.

Involvement of NF-κB in the reversal of CYP3A down-regulation induced by sea buckthorn in BCG-induced rats

Previous studies reported that sea buckthorn (Hippophae rhamnoides L., Elaeagnaceae, HRP) exhibits hepatoprotective effects via its anti-inflammatory and antioxidant properties as well as its inhibitory effects on collagen synthesis. However, it is unclear whether this hepatoprotective effect is also achieved by regulating liver drug metabolism enzyme pathways. Herein, we examined the regulatory effect of HRP on cytochrome P450 3A (CYP3A) in rats with immune liver injury, and explored the molecular mechanism of its hepatoprotective effect. Rat models of immunological liver injury were induced by intravenous injections of Bacillus Calmette-Guerin (BCG; 125 mg kg-1; 2 wks). Specific protein levels were detected by ELISA or western blot, and CYP3A mRNA expression was detected by RT-PCR. High-performance liquid chromatography (HPLC) detected relative changes in CYP3A metabolic activity based on the rates of 1-hydroxylation of the probe drug midazolam (MDZ). BCG pretreatment (125 mg kg-1) significantly down-regulated liver CYP3A protein expression compared with the control, metabolic activity, and transcription levels while up-regulating liver NF-κB, IL-1β, TNF-α and iNOS. HRP intervention (ED50: 78 mg kg-1) moderately reversed NF-κB, inflammatory cytokines, and iNOS activation in a dose-dependent manner (P < 0.05), and suppressed CYP3A down-regulation (P < 0.05); thereby partially alleviating liver injury. During immune liver injury, HRP may reverse CYP3A down-regulation by inhibiting NF-κB signal transduction, and protect liver function, which involves regulation of enzymes transcriptionally, translationally and post-translationally. The discovery that NF-κB is a molecular target of HRP may initiate the development and optimization of a clinical therapeutic approach to mitigate hepatitis B and other immunity-related liver diseases.

Contrasting model mechanisms of alanine aminotransferase (ALT) release from damaged and necrotic hepatocytes as an example of general biomarker mechanisms

Interpretations of elevated blood levels of alanine aminotransferase (ALT) for drug-induced liver injury often assume that the biomarker is released passively from dying cells. However, the mechanisms driving that release have not been explored experimentally. The usefulness of ALT and related biomarkers will improve by developing mechanism-based explanations of elevated levels that can be expanded and elaborated incrementally. We provide the means to challenge the ability of closely related model mechanisms to generate patterns of simulated hepatic injury and ALT release that scale (or not) to be quantitatively similar to the wet-lab validation targets, which are elevated plasma ALT values following acetaminophen (APAP) exposure in mice. We build on a published model mechanism that helps explain the generation of characteristic spatiotemporal features of APAP hepatotoxicity within hepatic lobules. Discrete event and agent-oriented software methods are most prominent. We instantiate and leverage a small constellation of concrete model mechanisms. Their details during execution help bring into focus ways in which particular sources of uncertainty become entangled with cause-effect details within and across several levels. We scale ALT amounts in virtual mice directly to target plasma ALT values in individual mice. A virtual experiment comprises a set of Monte Carlo simulations. We challenge the sufficiency of four potentially explanatory theories for ALT release. The first of the tested model theories failed to achieve the initial validation target, but each of the three others succeeded. Results for one of the three model mechanisms matched all target ALT values quantitatively. It explains how ALT externalization is the combined consequence of lobular-location-dependent drug-induced cellular damage and hepatocyte death. Falsification of one (or more) of the model mechanisms provides new knowledge and incrementally shrinks the constellation of model mechanisms. The modularity and biomimicry of our explanatory models enable seamless transition from mice to humans.

Propagation of Pericentral Necrosis During Acetaminophen-Induced Liver Injury: Evidence for Early Interhepatocyte Communication and Information Exchange

Acetaminophen (APAP)-induced liver injury is clinically significant, and APAP overdose in mice often serves as a model for drug-induced liver injury in humans. By specifying that APAP metabolism, reactive metabolite formation, glutathione depletion, and mitigation of mitochondrial damage within individual hepatocytes are functions of intralobular location, an earlier virtual model mechanism provided the first concrete multiattribute explanation for how and why early necrosis occurs close to the central vein (CV). However, two characteristic features could not be simulated consistently: necrosis occurring first adjacent to the CV, and subsequent necrosis occurring primarily adjacent to hepatocytes that have already initiated necrosis. We sought parsimonious model mechanism enhancements that would manage spatiotemporal heterogeneity sufficiently to enable meeting two new target attributes and conducted virtual experiments to explore different ideas for model mechanism improvement at intrahepatocyte and multihepatocyte levels. For the latter, evidence supports intercellular communication via exosomes, gap junctions, and connexin hemichannels playing essential roles in the toxic effects of chemicals, including facilitating or counteracting cell death processes. Logic requiring hepatocytes to obtain current information about whether downstream and lateral neighbors have triggered necrosis enabled virtual hepatocytes to achieve both new target attributes. A virtual hepatocyte that is glutathione-depleted uses that information to determine if it will initiate necrosis. When a less-stressed hepatocyte is flanked by at least two neighbors that have triggered necrosis, it too will initiate necrosis. We hypothesize that the resulting intercellular communication-enabled model mechanism is analogous to the actual explanation for APAP-induced hepatotoxicity at comparable levels of granularity.

Simulation enabled search for explanatory mechanisms of the fracture healing process

A significant portion of bone fractures fail to heal properly, increasing healthcare costs. Advances in fracture management have slowed because translation barriers have limited generation of mechanism-based explanations for the healing process. When uncertainties are numerous, analogical modeling can be an effective strategy for developing plausible explanations of complex phenomena. We demonstrate the feasibility of engineering analogical models in software to facilitate discovery of biomimetic explanations for how fracture healing may progress. Concrete analogical models—Callus Analogs—were created using the MASON simulation toolkit. We designated a Target Region initial state within a characteristic tissue section of mouse tibia fracture at day-7 and posited a corresponding day-10 Target Region final state. The goal was to discover a coarse-grain analog mechanism that would enable the discretized initial state to transform itself into the corresponding Target Region final state, thereby providing an alternative way to study the healing process. One of nine quasi-autonomous Tissue Unit types is assigned to each grid space, which maps to an 80×80 μm region of the tissue section. All Tissue Units have an opportunity each time step to act based on individualized logic, probabilities, and information about adjacent neighbors. Action causes transition from one Tissue Unit type to another, and simulation through several thousand time steps generates a coarse-grain analog—a theory—of the healing process. We prespecified a minimum measure of success: simulated and actual Target Region states achieve ≥ 70% Similarity. We used an iterative refinement protocol to explore many combinations of Tissue Unit logic and action constraints. Workflows progressed through four stages of analog mechanisms. Similarities of 73–90% were achieved for Mechanisms 2–4. The range of Upper-Level similarities increased to 83–94% when we allowed for uncertainty about two Tissue Unit designations. We have demonstrated how Callus Analog experiments provide domain experts with a fresh medium and tools for thinking about and understanding the fracture healing process.

Translation barriers have limited the generation of mechanism-based explanations of fracture healing processes. Those barriers help explain why, to date, biological therapeutics have had only a minor impact on fracture management. Alternative approaches are needed, and we present one that is intended to help develop incrementally better mechanism-based explanations of fracture healing phenomena. We created virtual Callus Analogs to simulate how the histologic appearance of a mouse fracture callus may transition from day-7 to day-10. Callus Analogs use software-based model mechanisms, and simulation experiments enable challenging and improving those model mechanisms. During execution, model mechanism operation provides a coarse-grain explanation (a theory) of a four-day portion of the healing process. Simulated day-10 callus histologic images achieved 73–94% Similarity to a corresponding day-10 fracture callus image, thus demonstrating feasibility. Simulated healing provides an alternative perspective on the actual healing process and an alternative way of thinking about plausible fracture healing mechanisms. Our working hypothesis is that the approach can be extended to cover more of the healing process while making features of simulated and actual fracture healing increasingly analogous. The methods presented are intended to be extensible to other research areas that use histologic analysis to investigate and explain tissue level phenomena.

The Virtual Anemia Trial: An Assessment of Model-Based In Silico Clinical Trials of Anemia Treatment Algorithms in Patients With Hemodialysis

In silico approaches have been proposed as a novel strategy to increase the repertoire of clinical trial designs. Realistic simulations of clinical trials can provide valuable information regarding safety and limitations of treatment protocols and have been shown to assist in the cost‐effective planning of clinical studies. In this report, we present a blueprint for the stepwise integration of internal, external, and ecological validity considerations in virtual clinical trials (VCTs). We exemplify this approach in the context of a model‐based in silico clinical trial aimed at anemia treatment in patients undergoing hemodialysis (HD). Hemoglobin levels and subsequent anemia treatment were simulated on a per patient level over the course of a year and compared to real‐life clinical data of 79,426 patients undergoing HD. The novel strategies presented here, aimed to improve external and ecological validity of a VCT, significantly increased the predictive power of the discussed in silico trial.

Competing Mechanistic Hypotheses of Acetaminophen-Induced Hepatotoxicity Challenged by Virtual Experiments

Acetaminophen-induced liver injury in mice is a model for drug-induced liver injury in humans. A precondition for improved strategies to disrupt and/or reverse the damage is a credible explanatory mechanism for how toxicity phenomena emerge and converge to cause hepatic necrosis. The Target Phenomenon in mice is that necrosis begins adjacent to the lobule's central vein (CV) and progresses outward. An explanatory mechanism remains elusive. Evidence supports that location dependent differences in NAPQI (the reactive metabolite) formation within hepatic lobules (NAPQI zonation) are necessary and sufficient prerequisites to account for that phenomenon. We call that the NZ-mechanism hypothesis. Challenging that hypothesis in mice is infeasible because 1) influential variables cannot be controlled, and 2) it would require sequential intracellular measurements at different lobular locations within the same mouse. Virtual hepatocytes use independently configured periportal-to-CV gradients to exhibit lobule-location dependent behaviors. Employing NZ-mechanism achieved quantitative validation targets for acetaminophen clearance and metabolism but failed to achieve the Target Phenomenon. We posited that, in order to do so, at least one additional feature must exhibit zonation by decreasing in the CV direction. We instantiated and explored two alternatives: 1) a glutathione depletion threshold diminishes in the CV direction; and 2) ability to repair mitochondrial damage diminishes in the CV direction. Inclusion of one or the other feature into NZ-mechanism failed to achieve the Target Phenomenon. However, inclusion of both features enabled successfully achieving the Target Phenomenon. The merged mechanism provides a multilevel, multiscale causal explanation of key temporal features of acetaminophen hepatotoxicity in mice. We discovered that variants of the merged mechanism provide plausible quantitative explanations for the considerable variation in 24-hour necrosis scores among 37 genetically diverse mouse strains following a single toxic acetaminophen dose.