Computational experiments reveal plausible mechanisms for changing patterns of hepatic zonation of xenobiotic clearance and hepatotoxicity

From Seeing to Simulating: A Survey of Imaging Techniques and Spatially-Resolved Data for Developing Multiscale Computational Models of Liver Regeneration

Liver regeneration, which leads to the re-establishment of organ mass, follows a specifically organized set of biological processes acting on various time and length scales. Computational models of liver regeneration largely focused on incorporating molecular and signaling detail have been developed by multiple research groups in the recent years. These modeling efforts have supported a synthesis of disparate experimental results at the molecular scale. Incorporation of tissue and organ scale data using noninvasive imaging methods can extend these computational models towards a comprehensive accounting of multiscale dynamics of liver regeneration. For instance, microscopy-based imaging methods provide detailed histological information at the tissue and cellular scales. Noninvasive imaging methods such as ultrasound, computed tomography and magnetic resonance imaging provide morphological and physiological features including volumetric measures over time. In this review, we discuss multiple imaging modalities capable of informing computational models of liver regeneration at the organ-, tissue- and cellular level. Additionally, we discuss available software and algorithms, which aid in the analysis and integration of imaging data into computational models. Such models can be generated or tuned for an individual patient with liver disease. Progress towards integrated multiscale models of liver regeneration can aid in prognostic tool development for treating liver disease.

Zonation of hepatic fat accumulation: insights from mathematical modelling of nutrient gradients and fatty acid uptake

Intrinsic of non-alcoholic fatty liver diseases is an aberrant accumulation of triglycerides (steatosis), which occurs inhomogeneously within lobules. To improve our understanding of the mechanisms involved in this zonation patterning, we developed a mathematical multicompartment model of hepatic fatty acid metabolism accompanied by blood flow simulations. A model analysis determines the influence of the uptake process of fatty acids, the porto-central gradient of plasma fatty acid concentration, and the oxygen supply via blood on the zonation of triglyceride accumulation. From this theoretical perspective, the plasma oxygen gradient, but not the fatty acid gradient, leads the way to a zonated triglyceride accumulation by its decisive role in oxidative processes. In addition, the uptake mechanism of fatty acids seems to be fundamental for a pericentral dominance of steatosis. However, the mechanism of cellular fatty acid uptake from the blood is still under debate. Our theoretical approach supports the transporter-mediated uptake mechanism and reveals that the maximal velocity of fatty acid uptake affects the switching between a periportal and a pericentral triglyceride accumulation. Further research on hepatic fatty acid uptake is needed to push forward our understanding of aberrant triglyceride accumulation in diet-induced steatosis.

A Computational Model of Hepatic Energy Metabolism: Understanding Zonated Damage and Steatosis in NAFLD

In non-alcoholic fatty liver disease (NAFLD), lipid build-up and the resulting damage is known to occur more severely in pericentral cells. Due to the complexity of studying individual regions of the sinusoid, the causes of this zone specificity and its implications on treatment are largely ignored. In this study, a computational model of liver glucose and lipid metabolism is presented which treats the sinusoid as the repeating unit of the liver rather than the single hepatocyte. This allows for inclusion of zonated enzyme expression by splitting the sinusoid into periportal to pericentral compartments. By simulating insulin resistance (IR) and high intake diets leading to the development of steatosis in the model, we identify key differences between periportal and pericentral cells accounting for higher susceptibility to pericentral steatosis. Secondly, variation between individuals is seen in both susceptibility to steatosis and in its development across the sinusoid. Around 25% of obese individuals do not show excess liver fat, whilst 16% of lean individuals develop NAFLD. Furthermore, whilst pericentral cells tend to show higher lipid levels, variation is seen in the predominant location of steatosis from pericentral to pan-sinusoidal or azonal. Sensitivity analysis was used to identify the processes which have the largest effect on both total hepatic triglyceride levels and on the sinusoidal location of steatosis. As is seen in vivo, steatosis occurs when simulating IR in the model, predominantly due to increased uptake, along with an increase in de novo lipogenesis. Additionally, concentrations of glucose intermediates including glycerol-3-phosphate increased when simulating IR due to inhibited glycogen synthesis. Several differences between zones contributed to a higher susceptibility to steatosis in pericentral cells in the model simulations. Firstly, the periportal zonation of both glycogen synthase and the oxidative phosphorylation enzymes meant that the build-up of glucose intermediates was less severe in the periportal hepatocyte compartments. Secondly, the periportal zonation of the enzymes mediating β-oxidation and oxidative phosphorylation resulted in excess fats being metabolised more rapidly in the periportal hepatocyte compartments. Finally, the pericentral expression of de novo lipogenesis contributed to pericentral steatosis when additionally simulating the increase in sterol-regulatory element binding protein 1c (SREBP-1c) seen in NAFLD patients in vivo. The hepatic triglyceride concentration was predicted to be most sensitive to inter-individual variation in the activity of enzymes which, either directly or indirectly, determine the rate of free fatty acid (FFA) oxidation. The concentration was most strongly dependent on the rate constants for β-oxidation and oxidative phosphorylation. It also showed moderate sensitivity to the rate constants for processes which alter the allosteric inhibition of β-oxidation by acetyl-CoA. The predominant sinusoidal location of steatosis meanwhile was most sensitive variations in the zonation of proteins mediating FFA uptake or triglyceride release as very low density lipoproteins (VLDL). Neither the total hepatic concentration nor the location of steatosis showed strong sensitivity to variations in the lipogenic rate constants.

Proposing a scientific confidence framework to help support the application of adverse outcome pathways for regulatory purposes

An adverse outcome pathway (AOP) describes the causal linkage between initial molecular events and an adverse outcome at individual or population levels. Whilst there has been considerable momentum in AOP development, far less attention has been paid to how AOPs might be practically applied for different regulatory purposes. This paper proposes a scientific confidence framework (SCF) for evaluating and applying a given AOP for different regulatory purposes ranging from prioritizing chemicals for further evaluation, to hazard prediction, and ultimately, risk assessment. The framework is illustrated using three different AOPs for several typical regulatory applications. The AOPs chosen are ones that have been recently developed and/or published, namely those for estrogenic effects, skin sensitisation, and rodent liver tumor promotion. The examples confirm how critical the data-richness of an AOP is for driving its practical application. In terms of performing risk assessment, human dosimetry methods are necessary to inform meaningful comparisons with human exposures; dosimetry is applied to effect levels based on non-testing approaches and in vitro data. Such a comparison is presented in the form of an exposure:activity ratio (EAR) to interpret biological activity in the context of exposure and to provide a basis for product stewardship and regulatory decision making.

Zonation of hepatic fatty acid metabolism - The diversity of its regulation and the benefit of modeling

A pronounced heterogeneity between hepatocytes in subcellular structure and enzyme activities was discovered more than 50years ago and initiated the idea of metabolic zonation. In the last decades zonation patterns of liver metabolism were extensively investigated for carbohydrate, nitrogen and lipid metabolism. The present review focuses on zonation patterns of the latter. We review recent findings regarding the zonation of fatty acid uptake and oxidation, ketogenesis, triglyceride synthesis and secretion, de novo lipogenesis, as well as bile acid and cholesterol metabolism. In doing so, we expose knowledge gaps and discuss contradictory experimental results, for example on the zonation pattern of fatty acid oxidation and de novo lipogenesis. Thus, possible rewarding directions of further research are identified. Furthermore, recent findings about the regulation of metabolic zonation are summarized, especially regarding the role of hormones, nerve innervation, morphogens, gender differences and the influence of the circadian clock. In the last part of the review, a short collection of models considering hepatic lipid metabolism is provided. We conclude that modeling, despite its proven benefit for understanding of hepatic carbohydrate and ammonia metabolisms, has so far been largely disregarded in the study of lipid metabolism; therefore some possible fields of modeling interest are presented.

Agent-based modeling: a systematic assessment of use cases and requirements for enhancing pharmaceutical research and development productivity

A crisis continues to brew within the pharmaceutical research and development (R&D) enterprise: productivity continues declining as costs rise, despite ongoing, often dramatic scientific and technical advances. To reverse this trend, we offer various suggestions for both the expansion and broader adoption of modeling and simulation (M&S) methods. We suggest strategies and scenarios intended to enable new M&S use cases that directly engage R&D knowledge generation and build actionable mechanistic insight, thereby opening the door to enhanced productivity. What M&S requirements must be satisfied to access and open the door, and begin reversing the productivity decline? Can current methods and tools fulfill the requirements, or are new methods necessary? We draw on the relevant, recent literature to provide and explore answers. In so doing, we identify essential, key roles for agent-based and other methods. We assemble a list of requirements necessary for M&S to meet the diverse needs distilled from a collection of research, review, and opinion articles. We argue that to realize its full potential, M&S should be actualized within a larger information technology framework—a dynamic knowledge repository—wherein models of various types execute, evolve, and increase in accuracy over time. We offer some details of the issues that must be addressed for such a repository to accrue the capabilities needed to reverse the productivity decline. © 2013 Wiley Periodicals, Inc.

Systems toxicology from genes to organs

This unique overview of systems toxicology methods and techniques begins with a brief account of systems thinking in biology over the last century. We discuss how systems biology and toxicology continue to leverage advances in computational modeling, informatics, large-scale computing, and biotechnology. Next, we chart the genesis of systems toxicology from previous work in physiologically based models, models of early development, and more recently, molecular systems biology. For readers interested in further details this background provides useful linkages to the relevant literature. It also lays the foundations for new ideas in systems toxicology that could translate laboratory measurements of molecular responses from xenobiotic perturbations to adverse organ level effects in humans. By providing innovative solutions across disciplinary boundaries and highlighting key scientific gaps, we believe this chapter provides useful information about the current state, and valuable insight about future directions in systems toxicity.

Individualized, discrete event, simulations provide insight into inter- and intra-subject variability of extended-release, drug products

OBJECTIVE

Develop and validate particular, concrete, and abstract yet plausible in silico mechanistic explanations for large intra- and interindividual variability observed for eleven bioequivalence study participants. Do so in the face of considerable uncertainty about mechanisms.

METHODS

We constructed an object-oriented, discrete event model called subject (we use small caps to distinguish computational objects from their biological counterparts). It maps abstractly to a dissolution test system and study subject to whom product was administered orally. A subject comprises four interconnected grid spaces and event mechanisms that map to different physiological features and processes. Drugs move within and between spaces. We followed an established, Iterative Refinement Protocol. Individualized mechanisms were made sufficiently complicated to achieve prespecified Similarity Criteria, but no more so. Within subjects, the dissolution space is linked to both a product-subject Interaction Space and the GI tract. The GI tract and Interaction Space connect to plasma, from which drug is eliminated.

RESULTS

We discovered parameterizations that enabled the eleven subject simulation results to achieve the most stringent Similarity Criteria. Simulated profiles closely resembled those with normal, odd, and double peaks. We observed important subject-by-formulation interactions within subjects.

CONCLUSION

We hypothesize that there were interactions within bioequivalence study participants corresponding to the subject-by-formulation interactions within subjects. Further progress requires methods to transition currently abstract subject mechanisms iteratively and parsimoniously to be more physiologically realistic. As that objective is achieved, the approach presented is expected to become beneficial to drug development (e.g., controlled release) and to a reduction in the number of subjects needed per study plus faster regulatory review.

Zonation related function and ubiquitination regulation in human hepatocellular carcinoma cells in dynamic vs. static culture conditions

Background

Understanding hepatic zonation is important both for liver physiology and pathology. There is currently no effective systemic chemotherapy for human hepatocellular carcinoma (HCC) and its pathogenesis is of special interest. Genomic and proteomic data of HCC cells in different culture models, coupled to pathway-based analysis, can help identify HCC-related gene and pathway dysfunctions.

Results

We identified zonation-related expression profiles contributing to selective phenotypes of HCC, by integrating relevant experimental observations through gene set enrichment analysis (GSEA). Analysis was based on gene and protein expression data measured on a human HCC cell line (HepG2/C3A) in two culture conditions: dynamic microfluidic biochips and static Petri dishes. Metabolic activity (HCC-related cytochromes P450) and genetic information processing were dominant in the dynamic cultures, in contrast to kinase signaling and cancer-specific profiles in static cultures. That, together with analysis of the published literature, leads us to propose that biochips culture conditions induce a periportal-like hepatocyte phenotype while standard plates cultures are more representative of a perivenous-like phenotype. Both proteomic data and GSEA results further reveal distinct ubiquitin-mediated protein regulation in the two culture conditions.

Conclusions

Pathways analysis, using gene and protein expression data from two cell culture models, confirmed specific human HCC phenotypes with regard to CYPs and kinases, and revealed a zonation-related pattern of expression. Ubiquitin-mediated regulation mechanism gives plausible explanations of our findings. Altogether, our results suggest that strategies aimed at inhibiting activated kinases and signaling pathways may lead to enhanced metabolism-mediated drug resistance of treated tumors. If that were the case, mitigating inhibition or targeting inactive forms of kinases would be an alternative.

Relational grounding facilitates development of scientifically useful multiscale models

We review grounding issues that influence the scientific usefulness of any biomedical multiscale model (MSM). Groundings are the collection of units, dimensions, and/or objects to which a variable or model constituent refers. To date, models that primarily use continuous mathematics rely heavily on absolute grounding, whereas those that primarily use discrete software paradigms (e.g., object-oriented, agent-based, actor) typically employ relational grounding. We review grounding issues and identify strategies to address them. We maintain that grounding issues should be addressed at the start of any MSM project and should be reevaluated throughout the model development process. We make the following points. Grounding decisions influence model flexibility, adaptability, and thus reusability. Grounding choices should be influenced by measures, uncertainty, system information, and the nature of available validation data. Absolute grounding complicates the process of combining models to form larger models unless all are grounded absolutely. Relational grounding facilitates referent knowledge embodiment within computational mechanisms but requires separate model-to-referent mappings. Absolute grounding can simplify integration by forcing common units and, hence, a common integration target, but context change may require model reengineering. Relational grounding enables synthesis of large, composite (multi-module) models that can be robust to context changes. Because biological components have varying degrees of autonomy, corresponding components in MSMs need to do the same. Relational grounding facilitates achieving such autonomy. Biomimetic analogues designed to facilitate translational research and development must have long lifecycles. Exploring mechanisms of normal-to-disease transition requires model components that are grounded relationally. Multi-paradigm modeling requires both hyperspatial and relational grounding.

Enabling clearance predictions to emerge from in silico actions of quasi-autonomous hepatocyte components

We demonstrate the feasibility of using in silico hepatocyte cultures (ISHCs) to provide predictions of the intrinsic clearance (CL) of compounds in hepatocyte cultures. We compare results with predictions obtained using a multiple linear regression method. Our expectation is that the method can be extended to predict in vivo clearance of new compounds in humans. Within ISHCs, mobile "compounds" carry information describing referent compound properties. We used an iterative refinement protocol for ISHC refinement and development of parameterization methods. Quasi-autonomous "hepatocytes" and their components (including "transporters" and "enzymes") use a small, event-specific subset of compound properties to decide how to interact with encountered compounds each simulation cycle. The probability of occurrence for each event is specified by a rule that uses a subset of compound properties known to influence that event in vitro. ISHC experiments mimic in vitro counterparts. In silico clearance is measured the same as in vitro clearance and is used to predict a corresponding CL value. For 39 of 73 compounds having calculated CL S.D.s, 79% of ISHC predictions and 23% of regression predictions were within CL ± 2 S.D. For all 73 compounds, 38% of ISHC predictions and 32% of regression predictions were within a factor of 2 of the referent CL values. ISHC details during simulations stand as a mechanistic hypothesis of how clearance phenomena emerge during in vitro experiments.