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

Getting in the zone: Metabolite transport across liver zones

The liver has many functions including the regulation of nutrient and metabolite levels in the systemic circulation through efficient transport into and out of hepatocytes. To sustain these functions, hepatocytes display large functional heterogeneity. This heterogeneity is reflected by zonation of metabolic processes that take place in different zones of the liver lobule, where nutrient-rich blood enters the liver in the periportal zone and flows through the mid-zone prior to drainage by a central vein in the pericentral zone. Metabolite transport plays a pivotal role in the division of labor across liver zones, being either transport into the hepatocyte or transport between hepatocytes through the blood. Signaling pathways that regulate zonation, such as Wnt/β-catenin, have been shown to play a causal role in the development of metabolic dysfunction-associated steatohepatitis (MASH) progression, but the (patho)physiological regulation of metabolite transport remains enigmatic. Despite the practical challenges to separately study individual liver zones, technological advancements in the recent years have greatly improved insight in spatially divided metabolite transport. This review summarizes the theories behind the regulation of zonation, diurnal rhythms and their effect on metabolic zonation, contemporary techniques used to study zonation and current technological challenges, and discusses the current view on spatial and temporal metabolite transport.

Modeling sex-specific whole-body metabolic responses to feeding and fasting

Men generally favor carbohydrate metabolism, while women lean towards lipid metabolism, resulting in significant sex-based differences in energy oxidation across various metabolic states such as fasting and feeding. These differences are influenced by body composition and inherent metabolic fluxes, including increased lipolysis rates in women. However, understanding how sex influences organ-specific metabolism and systemic manifestations remains incomplete. To address these gaps, we developed a sex-specific, whole-body metabolic model for feeding and fasting scenarios in healthy young adults. Our model integrates organ metabolism with whole-body responses to mixed meals, particularly high-carbohydrate and high-fat meals. Our predictions suggest that differences in liver and adipose tissue nutrient storage and oxidation patterns drive systemic metabolic disparities. We propose that sex differences in fasting hepatic glucose output may result from the different handling of free fatty acids, glycerol, and glycogen. We identified a metabolic pathway, possibly more prevalent in female livers, redirecting lipids towards carbohydrate metabolism to support hepatic glucose production. This mechanism is facilitated by the TG-FFA cycle between adipose tissue and the liver. Incorporating sex-specific data into multi-scale frameworks offers insights into how sex modulates human metabolism.

A unified framework for prediction of liver steatosis dynamics in response to different diet and drug interventions

BACKGROUND & AIMS

Non-alcoholic fatty liver disease (NAFLD) is a common metabolic disorder, characterized by the accumulation of excess fat in the liver, and is a driving factor for various severe liver diseases. These multi-factorial and multi-timescale changes are observed in different clinical studies, but these studies have not been integrated into a unified framework. In this study, we aim to present such a unified framework in the form of a dynamic mathematical model.

METHODS

For model training and validation, we collected data for dietary or drug-induced interventions aimed at reducing or increasing liver fat. The model was formulated using ordinary differential equations (ODEs) and the mathematical analysis, model simulation, model formulation and the model parameter estimation were all performed in MATLAB.

RESULTS

Our mathematical model describes accumulation of fat in the liver and predicts changes in lipid fluxes induced by both dietary and drug interventions. The model is validated using data from a wide range of drug and dietary intervention studies and can predict both short-term (days) and long-term (weeks) changes in liver fat. Importantly, the model computes the contribution of each individual lipid flux to the total liver fat dynamics. Furthermore, the model can be combined with an established bodyweight model, to simulate even longer scenarios (years), also including the effects of insulin resistance and body weight. To help prepare for corresponding eHealth applications, we also present a way to visualize the simulated changes, using dynamically changing lipid droplets, seen in images of liver biopsies.

CONCLUSION

In conclusion, we believe that the minimal model presented herein might be a useful tool for future applications, and to further integrate and understand data regarding changes in dietary and drug induced changes in ectopic TAG in the liver. With further development and validation, the minimal model could be used as a disease progression model for steatosis.

Quantifying fat zonation in liver lobules: an integrated multiscale in silico model combining disturbed microperfusion and fat metabolism via a continuum biomechanical bi-scale, tri-phasic approach

Metabolic zonation refers to the spatial separation of metabolic functions along the sinusoidal axes of the liver. This phenomenon forms the foundation for adjusting hepatic metabolism to physiological requirements in health and disease (e.g., metabolic dysfunction-associated steatotic liver disease/MASLD). Zonated metabolic functions are influenced by zonal morphological abnormalities in the liver, such as periportal fibrosis and pericentral steatosis. We aim to analyze the interplay between microperfusion, oxygen gradient, fat metabolism and resulting zonated fat accumulation in a liver lobule. Therefore we developed a continuum biomechanical, tri-phasic, bi-scale, and multicomponent in silico model, which allows to numerically simulate coupled perfusion-function-growth interactions two-dimensionally in liver lobules. The developed homogenized model has the following specifications: (i) thermodynamically consistent, (ii) tri-phase model (tissue, fat, blood), (iii) penta-substances (glycogen, glucose, lactate, FFA, and oxygen), and (iv) bi-scale approach (lobule, cell). Our presented in silico model accounts for the mutual coupling between spatial and time-dependent liver perfusion, metabolic pathways and fat accumulation. The model thus allows the prediction of fat development in the liver lobule, depending on perfusion, oxygen and plasma concentration of free fatty acids (FFA), oxidative processes, the synthesis and the secretion of triglycerides (TGs). The use of a bi-scale approach allows in addition to focus on scale bridging processes. Thus, we will investigate how changes at the cellular scale affect perfusion at the lobular scale and vice versa. This allows to predict the zonation of fat distribution (periportal or pericentral) depending on initial conditions, as well as external and internal boundary value conditions.

Virtual clinical trials of anti-PD-1 and anti-CTLA-4 immunotherapy in advanced hepatocellular carcinoma using a quantitative systems pharmacology model

BACKGROUND

Hepatocellular carcinoma (HCC) is the most common form of primary liver cancer and is the third-leading cause of cancer-related death worldwide. Most patients with HCC are diagnosed at an advanced stage, and the median survival for patients with advanced HCC treated with modern systemic therapy is less than 2 years. This leaves the advanced stage patients with limited treatment options. Immune checkpoint inhibitors (ICIs) targeting programmed cell death protein 1 (PD-1) or its ligand, are widely used in the treatment of HCC and are associated with durable responses in a subset of patients. ICIs targeting cytotoxic T-lymphocyte-associated protein 4 (CTLA-4) also have clinical activity in HCC. Combination therapy of nivolumab (anti-PD-1) and ipilimumab (anti-CTLA-4) is the first treatment option for HCC to be approved by Food and Drug Administration that targets more than one immune checkpoints.

METHODS

In this study, we used the framework of quantitative systems pharmacology (QSP) to perform a virtual clinical trial for nivolumab and ipilimumab in HCC patients. Our model incorporates detailed biological mechanisms of interactions of immune cells and cancer cells leading to antitumor response. To conduct virtual clinical trial, we generate virtual patient from a cohort of 5,000 proposed patients by extending recent algorithms from literature. The model was calibrated using the data of the clinical trial CheckMate 040 (ClinicalTrials.gov number, NCT01658878).

RESULTS

Retrospective analyses were performed for different immune checkpoint therapies as performed in CheckMate 040. Using machine learning approach, we predict the importance of potential biomarkers for immune blockade therapies.

CONCLUSIONS

This is the first QSP model for HCC with ICIs and the predictions are consistent with clinically observed outcomes. This study demonstrates that using a mechanistic understanding of the underlying pathophysiology, QSP models can facilitate patient selection and design clinical trials with improved success.

Integrated experimental-computational analysis of a HepaRG liver-islet microphysiological system for human-centric diabetes research

Microphysiological systems (MPS) are powerful tools for emulating human physiology and replicating disease progression in vitro. MPS could be better predictors of human outcome than current animal models, but mechanistic interpretation and in vivo extrapolation of the experimental results remain significant challenges. Here, we address these challenges using an integrated experimental-computational approach. This approach allows for in silico representation and predictions of glucose metabolism in a previously reported MPS with two organ compartments (liver and pancreas) connected in a closed loop with circulating medium. We developed a computational model describing glucose metabolism over 15 days of culture in the MPS. The model was calibrated on an experiment-specific basis using data from seven experiments, where HepaRG single-liver or liver-islet cultures were exposed to both normal and hyperglycemic conditions resembling high blood glucose levels in diabetes. The calibrated models reproduced the fast (i.e. hourly) variations in glucose and insulin observed in the MPS experiments, as well as the long-term (i.e. over weeks) decline in both glucose tolerance and insulin secretion. We also investigated the behaviour of the system under hypoglycemia by simulating this condition in silico, and the model could correctly predict the glucose and insulin responses measured in new MPS experiments. Last, we used the computational model to translate the experimental results to humans, showing good agreement with published data of the glucose response to a meal in healthy subjects. The integrated experimental-computational framework opens new avenues for future investigations toward disease mechanisms and the development of new therapies for metabolic disorders.

Microphysiological systems (MPS) are powerful tools to unravel biological knowledge underlying disease. MPS provide a physiologically relevant, human-based in vitro setting, which can potentially yield better translatability to humans than current animal models and traditional cell cultures. However, mechanistic interpretation and extrapolation of the experimental results to human outcome remain significant challenges. In this study, we confront these challenges using an integrated experimental-computational approach. We present a computational model describing glucose metabolism in a previously reported MPS integrating liver and pancreas. This MPS supports a homeostatic feedback loop between HepaRG/HHSteC spheroids and pancreatic islets, and allows for detailed investigations of mechanisms underlying type 2 diabetes in humans. We show that the computational model captures the complex dynamics of glucose-insulin regulation observed in the system, and can provide mechanistic insight into disease progression features, such as insulin resistance and β-cell dynamics. Furthermore, the computational model can explain key differences in temporal dynamics between MPS and human responses, and thus provides a tool for translating experimental insights into human outcome. The integrated experimental-computational framework opens new avenues for future investigations toward disease mechanisms and the development of new therapies for metabolic disorders.

Digital twin predicting diet response before and after long-term fasting

Today, there is great interest in diets proposing new combinations of macronutrient compositions and fasting schedules. Unfortunately, there is little consensus regarding the impact of these different diets, since available studies measure different sets of variables in different populations, thus only providing partial, non-connected insights. We lack an approach for integrating all such partial insights into a useful and interconnected big picture. Herein, we present such an integrating tool. The tool uses a novel mathematical model that describes mechanisms regulating diet response and fasting metabolic fluxes, both for organ-organ crosstalk, and inside the liver. The tool can mechanistically explain and integrate data from several clinical studies, and correctly predict new independent data, including data from a new study. Using this model, we can predict non-measured variables, e.g. hepatic glycogen and gluconeogenesis, in response to fasting and different diets. Furthermore, we exemplify how such metabolic responses can be successfully adapted to a specific individual’s sex, weight, height, as well as to the individual’s historical data on metabolite dynamics. This tool enables an offline digital twin technology.

Fasting and diet are central components of prevention against cardiovascular disease. Unfortunately, there is little consensus regarding which diet schemes are optimal. This is partially because different clinical studies contribute with different non-connected pieces of knowledge, which have not been fully integrated into a useful and interconnected big picture. In principle, mathematical models describing meal responses could be used for such an integration. However, today’s models still lack critical mechanisms, such as protein metabolism and a dynamic glycogen regulation. Herein, we present a) a new expanded model structure including these mechanisms; b) a set of parameters which can simultaneously describe a wide array of complementary estimation data, in both healthy and diabetic populations; c) a personalisation-script, which allows these generic parameters to be tuned to an individual/sub-population, using demographics (age, weight, height, diabetes status) and historic metabolic data. We exemplify how this personalisation can be used to predict new independent data, including a new clinical study, where a qualitatively new prediction is validated: that an oral protein tolerance test gives a clear response in plasma glucose, after, but not before, a 48h fasting period. Our combined model, parameters, and fitting script lay the foundation for an offline digital twin.

Applications of Quantitative Systems Pharmacology (QSP) in Drug Development for NAFLD and NASH and Its Regulatory Application

Nonalcoholic steatohepatitis (NASH) is a widely prevalent disease, but approved pharmaceutical treatments are not available. As such, there is great activity within the pharmaceutical industry to accelerate drug development in this area and improve the quality of life and reduce mortality for NASH patients. The use of quantitative systems pharmacology (QSP) can help make this overall process more efficient. This mechanism-based mathematical modeling approach describes both the pathophysiology of a disease and how pharmacological interventions can modify pathophysiologic mechanisms. Multiple capabilities are provided by QSP modeling, including the use of model predictions to optimize clinical studies. The use of this approach has grown over the last 20 years, motivating discussions between modelers and regulators to agree upon methodologic standards. These include model transparency, documentation, and inclusion of clinical pharmacodynamic biomarkers. Several QSP models have been developed that describe NASH pathophysiology to varying extents. One specific application of NAFLDsym, a QSP model of NASH, is described in this manuscript. Simulations were performed to help understand if patient behaviors could help explain the relatively high rate of fibrosis stage reductions in placebo cohorts. Simulated food intake and body weight fluctuated periodically over time. The relatively slow turnover of liver collagen allowed persistent reductions in predicted fibrosis stage despite return to baseline for liver fat, plasma ALT, and the NAFLD activity score. Mechanistic insights such as this that have been derived from QSP models can help expedite the development of safe and effective treatments for NASH patients.

A biologically based model to quantitatively assess the role of the nuclear receptors liver X (LXR), and pregnane X (PXR) on chemically induced hepatic steatosis

Hepatic steatosis is characterized by the intracellular increase of free fatty acids (FFAs) in the form of triglycerides in hepatocytes. This hepatic adverse outcome can be caused by many factors, including exposure to drugs or environmental toxicants. Mechanistically, accumulation of lipids in the liver can take place via several mechanisms such as de novo synthesis and/or uptake of FFAs from serum via high fat content diets. De novo synthesis of FFAs within the liver is mediated by the liver X receptor (LXR), and their uptake into the liver is mediated through the pregnane X receptor (PXR). We investigated the impact of chemical exposure on FFAs hepatic content via activation of LXR and PXR by integrating chemical-specific physiologically based pharmacokinetic (PBPK) models with a quantitative toxicology systems (QTS) model of hepatic lipid homeostasis. Three known agonists of LXR and/or PXR were modeled: T0901317 (antagonist for both receptors), GW3965 (LXR only), and Rifampicin (PXR only). Model predictions showed that T0901317 caused the most FFAs build-up in the liver, followed by Rifampicin and then GW3965. These modeling results highlight the importance of PXR activation for serum FFAs uptake into the liver while suggesting that increased hepatic FAAs de novo synthesis alone may not be enough to cause appreciable accumulation of lipids in the liver under normal environmental exposure levels. Moreover, the overall PBPK-hepatic lipids quantitative model can be used to screen chemicals for their potential to cause in vivo hepatic lipid content buildup in view of their in vitro potential to activate the nuclear receptors and their exposure levels.

A Spatial Model of Hepatic Calcium Signaling and Glucose Metabolism Under Autonomic Control Reveals Functional Consequences of Varying Liver Innervation Patterns Across Species

Rapid breakdown of hepatic glycogen stores into glucose plays an important role during intense physical exercise to maintain systemic euglycemia. Hepatic glycogenolysis is governed by several different liver-intrinsic and systemic factors such as hepatic zonation, circulating catecholamines, hepatocellular calcium signaling, hepatic neuroanatomy, and the central nervous system (CNS). Of the factors regulating hepatic glycogenolysis, the extent of lobular innervation varies significantly between humans and rodents. While rodents display very few autonomic nerve terminals in the liver, nearly every hepatic layer in the human liver receives neural input. In the present study, we developed a multi-scale, multi-organ model of hepatic metabolism incorporating liver zonation, lobular scale calcium signaling, hepatic innervation, and direct and peripheral organ-mediated communication between the liver and the CNS. We evaluated the effect of each of these governing factors on the total hepatic glucose output and zonal glycogenolytic patterns within liver lobules during simulated physical exercise. Our simulations revealed that direct neuronal stimulation of the liver and an increase in circulating catecholamines increases hepatic glucose output mediated by mobilization of intracellular calcium stores and lobular scale calcium waves. Comparing simulated glycogenolysis between human-like and rodent-like hepatic innervation patterns (extensive vs. minimal) suggested that propagation of calcium transients across liver lobules acts as a compensatory mechanism to improve hepatic glucose output in sparsely innervated livers. Interestingly, our simulations suggested that catecholamine-driven glycogenolysis is reduced under portal hypertension. However, increased innervation coupled with strong intercellular communication can improve the total hepatic glucose output under portal hypertension. In summary, our modeling and simulation study reveals a complex interplay of intercellular and multi-organ interactions that can lead to differing calcium dynamics and spatial distributions of glycogenolysis at the lobular scale in the liver.

Citrus Flavanone Narirutin, In Vitro and In Silico Mechanistic Antidiabetic Potential

Citrus fruits and juices have been studied extensively for their potential involvement in the prevention of various diseases. Flavanones, the characteristic polyphenols of citrus species, are the primarily compounds responsible for these studied health benefits. Using in silico and in vitro methods, we are exploring the possible antidiabetic action of narirutin, a flavanone family member. The goal of the in silico research was to anticipate how narirutin would interact with eight distinct receptors implicated in diabetes control and complications, namely, dipeptidyl-peptidase 4 (DPP4), protein tyrosine phosphatase 1B (PTP1B), free fatty acid receptor 1 (FFAR1), aldose reductase (AldR), glycogen phosphorylase (GP), alpha-amylase (AAM), peroxisome proliferator-activated receptor gamma (PPAR-γ), alpha-glucosidase (AGL), while the in vitro study looked into narirutin’s possible inhibitory impact on alpha-amylase and alpha-glucosidase. The results indicate that the studied citrus flavanone interacted remarkably with most of the receptors and had an excellent inhibitory activity during the in vitro tests suggesting its potent role among the different constituent of the citrus compounds in the management of diabetes and also its complications.

Application of Magnetic Resonance Imaging in Liver Biomechanics: A Systematic Review

MRI-based biomechanical studies can provide a deep understanding of the mechanisms governing liver function, its mechanical performance but also liver diseases. In addition, comprehensive modeling of the liver can help improve liver disease treatment. Furthermore, such studies demonstrate the beginning of an engineering-level approach to how the liver disease affects material properties and liver function. Aimed at researchers in the field of MRI-based liver simulation, research articles pertinent to MRI-based liver modeling were identified, reviewed, and summarized systematically. Various MRI applications for liver biomechanics are highlighted, and the limitations of different viscoelastic models used in magnetic resonance elastography are addressed. The clinical application of the simulations and the diseases studied are also discussed. Based on the developed questionnaire, the papers' quality was assessed, and of the 46 reviewed papers, 32 papers were determined to be of high-quality. Due to the lack of the suitable material models for different liver diseases studied by magnetic resonance elastography, researchers may consider the effect of liver diseases on constitutive models. In the future, research groups may incorporate various aspects of machine learning (ML) into constitutive models and MRI data extraction to further refine the study methodology. Moreover, researchers should strive for further reproducibility and rigorous model validation and verification.

Computational Hypothesis: How Intra-Hepatic Functional Heterogeneity May Influence the Cascading Progression of Free Fatty Acid-Induced Non-Alcoholic Fatty Liver Disease (NAFLD)

Non-Alcoholic Fatty Liver Disease (NAFLD) is the most common type of chronic liver disease in developed nations, affecting around 25% of the population. Elucidating the factors causing NAFLD in individual patients to progress in different rates and to different degrees of severity, is a matter of active medical research. Here, we aim to provide evidence that the intra-hepatic heterogeneity of rheological, metabolic and tissue-regenerating capacities plays a central role in disease progression. We developed a generic mathematical model that constitutes the liver as ensemble of small liver units differing in their capacities to metabolize potentially cytotoxic free fatty acids (FFAs) and to repair FFA-induced cell damage. Transition from simple steatosis to more severe forms of NAFLD is described as self-amplifying process of cascading liver failure, which, to stop, depends essentially on the distribution of functional capacities across the liver. Model simulations provided the following insights: (1) A persistently high plasma level of FFAs is sufficient to drive the liver through different stages of NAFLD; (2) Presence of NAFLD amplifies the deleterious impact of additional tissue-damaging hits; and (3) Coexistence of non-steatotic and highly steatotic regions is indicative for the later occurrence of severe NAFLD stages.

Virtual metabolic human dynamic model for pathological analysis and therapy design for diabetes

A virtual metabolic human model is a valuable complement to experimental biology and clinical studies, because in vivo research involves serious ethical and technical problems. I have proposed a multi-organ and multi-scale kinetic model that formulates the reactions of enzymes and transporters with the regulation of hormonal actions at postprandial and postabsorptive states. The computational model consists of 202 ordinary differential equations for metabolites with 217 reaction rates and 1,140 kinetic parameter constants. It is the most comprehensive, largest, and highly predictive model of the whole-body metabolism. Use of the model revealed the mechanisms by which individual disorders, such as steatosis, β cell dysfunction, and insulin resistance, were combined to cause diabetes. The model predicted a glycerol kinase inhibitor to be an effective medicine for type 2 diabetes, which not only decreased hepatic triglyceride but also reduced plasma glucose. The model also enabled us to rationally design combination therapy.

•

A standard of virtual metabolic human dynamic models is proposed

•

It integrates the three scales of molecules, organs, and whole body

•

It gets insight into pathological mechanisms of type 1 and type 2 diabetes

•

It enables the computer-aided design of medication treatment for diabetes

TCTP participates in hepatic metabolism by regulating gene expression involved in insulin resistance

Translationally controlled tumor protein (TCTP) has various cellular functions and molecular interactions, many related to its growth-promoting and antiapoptotic properties. Recently, TCTP expression was reported to increases in insulin-resistant mice fed with high-fat diet. TCTP is a multifunctional protein, but its role in liver metabolism is unclear. Here, we investigated the function and mechanism of TCTP in HepG2 cells. Knock-down of TCTP led to 287 differentially expressed genes (DEGs) that were highly associated with cellular apoptosis and signal response, TNF and NF-κB signaling pathways, glycolysis/gluconeogenesis, insulin resistance, FoxO and insulin signaling pathways, adipocytokine and AMPK signaling pathways. shTCTP downregulated the expression of the key gluconeogenesis enzyme phosphoenolpyruvate carboxykinase (PCK1). Furthermore, TCTP regulated the alternative splicing of genes enriched in the phospholipid biosynthetic process and glycerophospholipid metabolism. We further showed that shTCTP down-regulated the intracellular levels of triglyceride and total cholesterol. Our results showed that TCTP regulates the liver cell transcriptome at both the transcriptional and alternative splicing levels. The TCTP regulatory network predicts the biological functions of TCTP in glucose and lipid metabolism, and also insulin resistance, which may be associated with liver metabolism and diseases such as nonalcoholic fatty liver disease.

Computational Modeling of Fructose Metabolism and Development in NAFLD

Non-alcohol fatty liver disease (NAFLD) is a common disorder that has increased in prevalence 20-fold over the last three decades. It covers a spectrum of conditions resulting from excess lipid accumulation in the liver without alcohol abuse. Among all the risk factors, over-consumption of fructose has been repeatedly reported in both clinical and experimental studies to be highly associated with the development of NAFLD. However, studying in vivo systems is complicated, time consuming and expensive. A detailed kinetic model of fructose metabolism was constructed to investigate the metabolic mechanisms whereby fructose consumption can induce dyslipidaemia associated with NAFLD and to explore whether the pathological conditions can be reversed during the early stages of disease. The model contains biochemical components and reactions identified from the literature, including ~120 parameters, 25 variables, and 25 first order differential equations. Three scenarios were presented to demonstrate the behavior of the model. Scenario one predicts the acute effects of a change in carbohydrate input in lipid profiles. The results present progressive triglyceride accumulations in the liver and plasma for three diets. The rate of accumulation was greater in the fructose diet than that of the mixed or glucose only models. Scenario two explores the variability of metabolic reaction rate within the general population. Sensitivity analysis reveals that hepatic triglyceride concentration is most sensitive to the rate constant of pyruvate kinase and fructokinase. Scenario three tests the effect of one specific inhibitor that might be potentially administered. The simulations of fructokinase suppression provide a good model for potentially reversing simple steatosis induced by high fructose consumption, which can be corroborated by experimental studies. The predictions in these three scenarios suggest that the model is robust and it has sufficient detail to present the kinetic relationship between fructose and lipid in the liver.

A mechanobiological mathematical model of liver metabolism

The liver plays a complex role in metabolism and detoxification, and better tools are needed to understand its function and to develop liver-targeted therapies. In this study, we establish a mechanobiological model of liver transport and hepatocyte biology to elucidate the metabolism of urea and albumin, the production/detoxification of ammonia, and consumption of oxygen and nutrients. Since hepatocellular shear stress (SS) can influence the enzymatic activities of liver, the effect of SS on the urea and albumin synthesis are empirically modeled through the mechanotransduction mechanisms. The results demonstrate that the rheology and dynamics of the sinusoid flow can significantly affect liver metabolism. We show that perfusate rheology and blood hematocrit can affect urea and albumin production by changing hepatocyte mechanosensitive metabolism. The model can also simulate enzymatic diseases of the liver such as hyperammonemia I, hyperammonemia II, hyperarginemia, citrollinemia, and argininosuccinicaciduria, which disrupt the urea metabolism and ammonia detoxification. The model is also able to predict how aggregate cultures of hepatocytes differ from single cell cultures. We conclude that in vitro perfusable devices for the study of liver metabolism or personalized medicine should be designed with similar morphology and fluid dynamics as patient liver tissue. This robust model can be adapted to any type of hepatocyte culture to determine how hepatocyte viability, functionality, and metabolism are influenced by liver pathologies and environmental conditions.

Emergent Properties of the HNF4α-PPARγ Network May Drive Consequent Phenotypic Plasticity in NAFLD

Non-alcoholic fatty liver disease (NAFLD) is the most common form of chronic liver disease in adults and children. It is characterized by excessive accumulation of lipids in the hepatocytes of patients without any excess alcohol intake. With a global presence of 24% and limited therapeutic options, the disease burden of NAFLD is increasing. Thus, it becomes imperative to attempt to understand the dynamics of disease progression at a systems-level. Here, we decoded the emergent dynamics of underlying gene regulatory networks that were identified to drive the initiation and the progression of NAFLD. We developed a mathematical model to elucidate the dynamics of the HNF4α-PPARγ gene regulatory network. Our simulations reveal that this network can enable multiple co-existing phenotypes under certain biological conditions: an adipocyte, a hepatocyte, and a "hybrid" adipocyte-like state of the hepatocyte. These phenotypes may also switch among each other, thus enabling phenotypic plasticity and consequently leading to simultaneous deregulation of the levels of molecules that maintain a hepatic identity and/or facilitate a partial or complete acquisition of adipocytic traits. These predicted trends are supported by the analysis of clinical data, further substantiating the putative role of phenotypic plasticity in driving NAFLD. Our results unravel how the emergent dynamics of underlying regulatory networks can promote phenotypic plasticity, thereby propelling the clinically observed changes in gene expression often associated with NAFLD.

Liver Bioreactor Design Issues of Fluid Flow and Zonation, Fibrosis, and Mechanics: A Computational Perspective

Tissue engineering, with the goal of repairing or replacing damaged tissue and organs, has continued to make dramatic science-based advances since its origins in the late 1980’s and early 1990’s. Such advances are always multi-disciplinary in nature, from basic biology and chemistry through physics and mathematics to various engineering and computer fields. This review will focus its attention on two topics critical for tissue engineering liver development: (a) fluid flow, zonation, and drug screening, and (b) biomechanics, tissue stiffness, and fibrosis, all within the context of 3D structures. First, a general overview of various bioreactor designs developed to investigate fluid transport and tissue biomechanics is given. This includes a mention of computational fluid dynamic methods used to optimize and validate these designs. Thereafter, the perspective provided by computer simulations of flow, reactive transport, and biomechanics responses at the scale of the liver lobule and liver tissue is outlined, in addition to how bioreactor-measured properties can be utilized in these models. Here, the fundamental issues of tortuosity and upscaling are highlighted, as well as the role of disease and fibrosis in these issues. Some idealized simulations of the effects of fibrosis on lobule drug transport and mechanics responses are provided to further illustrate these concepts. This review concludes with an outline of some practical applications of tissue engineering advances and how efficient computational upscaling techniques, such as dual continuum modeling, might be used to quantify the transition of bioreactor results to the full liver scale.

Enzymatic and nonenzymatic protein acetylations control glycolysis process in liver diseases

Impaired glycolysis has pathologic effects on the occurrence and progression of liver diseases, and it appears that glycolysis is increased to different degrees in different liver diseases. As an important post-translational modification, reversible lysine acetylation regulates almost all cellular processes, including glycolysis. Lysine acetylation can occur enzymatically with acetyltransferases or nonenzymatically with acetyl-coenzyme A. Accompanied by the progression of liver diseases, there seems to be a temporal and spatial variation between enzymatic and nonenzymatic acetylations in the regulation of glycolysis. Here, we summarize the most recent findings on the functions and targets of acetylation in controlling glycolysis in the different stages of liver diseases. In addition, we discuss the differences and causes between enzymatic and nonenzymatic acetylations in regulating glycolysis throughout the progression of liver diseases. Then, we review these new discoveries to provide the potential implications of these findings for therapeutic interventions in liver diseases.-Li, J., Wang, T., Xia, J., Yao, W., Huang, F. Enzymatic and nonenzymatic protein acetylations control glycolysis process in liver diseases.

Harnessing Human Microphysiology Systems as Key Experimental Models for Quantitative Systems Pharmacology

Two technologies that have emerged in the last decade offer a new paradigm for modern pharmacology, as well as drug discovery and development. Quantitative systems pharmacology (QSP) is a complementary approach to traditional, target-centric pharmacology and drug discovery and is based on an iterative application of computational and systems biology methods with multiscale experimental methods, both of which include models of ADME-Tox and disease. QSP has emerged as a new approach due to the low efficiency of success in developing therapeutics based on the existing target-centric paradigm. Likewise, human microphysiology systems (MPS) are experimental models complementary to existing animal models and are based on the use of human primary cells, adult stem cells, and/or induced pluripotent stem cells (iPSCs) to mimic human tissues and organ functions/structures involved in disease and ADME-Tox. Human MPS experimental models have been developed to address the relatively low concordance of human disease and ADME-Tox with engineered, experimental animal models of disease. The integration of the QSP paradigm with the use of human MPS has the potential to enhance the process of drug discovery and development.

Stem cell biomanufacturing under uncertainty: A case study in optimizing red blood cell production

As breakthrough cellular therapy discoveries are translated into reliable, commercializable applications, effective stem cell biomanufacturing requires systematically developing and optimizing bioprocess design and operation. This article proposes a rigorous computational framework for stem cell biomanufacturing under uncertainty. Our mathematical tool kit incorporates: high‐fidelity modeling, single variate and multivariate sensitivity analysis, global topological superstructure optimization, and robust optimization. The advantages of the proposed bioprocess optimization framework using, as a case study, a dual hollow fiber bioreactor producing red blood cells from progenitor cells were quantitatively demonstrated. The optimization phase reduces the cost by a factor of 4, and the price of insuring process performance against uncertainty is approximately 15% over the nominal optimal solution. Mathematical modeling and optimization can guide decision making; the possible commercial impact of this cellular therapy using the disruptive technology paradigm was quantitatively evaluated. © 2017 American Institute of Chemical Engineers AIChE J, 64: 3011–3022, 2018

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 multiscale modelling approach to assess the impact of metabolic zonation and microperfusion on the hepatic carbohydrate metabolism

The capacity of the liver to convert the metabolic input received from the incoming portal and arterial blood into the metabolic output of the outgoing venous blood has three major determinants: The intra-hepatic blood flow, the transport of metabolites between blood vessels (sinusoids) and hepatocytes and the metabolic capacity of hepatocytes. These determinants are not constant across the organ: Even in the normal organ, but much more pronounced in the fibrotic and cirrhotic liver, regional variability of the capillary blood pressure, tissue architecture and the expression level of metabolic enzymes (zonation) have been reported. Understanding how this variability may affect the regional metabolic capacity of the liver is important for the interpretation of functional liver tests and planning of pharmacological and surgical interventions. Here we present a mathematical model of the sinusoidal tissue unit (STU) that is composed of a single sinusoid surrounded by the space of Disse and a monolayer of hepatocytes. The total metabolic output of the liver (arterio-venous glucose difference) is obtained by integration across the metabolic output of a representative number of STUs. Application of the model to the hepatic glucose metabolism provided the following insights: (i) At portal glucose concentrations between 6–8 mM, an intra-sinusoidal glucose cycle may occur which is constituted by glucose producing periportal hepatocytes and glucose consuming pericentral hepatocytes, (ii) Regional variability of hepatic blood flow is higher than the corresponding regional variability of the metabolic output, (iii) a spatially resolved metabolic functiogram of the liver is constructed. Variations of tissue parameters are equally important as variations of enzyme activities for the control of the arterio-venous glucose difference.

Glucose homeostasis is one of the central liver functions. The liver extracts glucose from the blood when plasma glucose levels are high and produces glucose when plasma glucose levels are low. To fulfill this function the liver is organized in smallest functional units, the sinusoidal tissue units (STUs). These STUs consist of a single sinusoid surrounded by linear arranged hepatocytes. Liver zonation describes the spatial separation of metabolic pathways along the STUs. As blood flows through the sinusoid the plasma nutrient and hormone composition changes and in conjunction with the heterogeneous endowment of metabolic enzymes this leads to big differences in the metabolic performance of hepatocytes depending on their position within the sinusoid. This makes liver zonation and blood flow two central determinants for the functional output of the liver. In this work we present a tissue model of hepatic carbohydrate metabolism that combines liver zonation and microperfusion within the STU. We show that structural properties, enzymatic properties and regional bloodflow are equally important for the understanding of liver functionality. With our work we provide a true multi-scale model bridging the scale from the cellular to the tissue level.

Computational Modeling in Liver Surgery

The need for extended liver resection is increasing due to the growing incidence of liver tumors in aging societies. Individualized surgical planning is the key for identifying the optimal resection strategy and to minimize the risk of postoperative liver failure and tumor recurrence. Current computational tools provide virtual planning of liver resection by taking into account the spatial relationship between the tumor and the hepatic vascular trees, as well as the size of the future liver remnant. However, size and function of the liver are not necessarily equivalent. Hence, determining the future liver volume might misestimate the future liver function, especially in cases of hepatic comorbidities such as hepatic steatosis. A systems medicine approach could be applied, including biological, medical, and surgical aspects, by integrating all available anatomical and functional information of the individual patient. Such an approach holds promise for better prediction of postoperative liver function and hence improved risk assessment. This review provides an overview of mathematical models related to the liver and its function and explores their potential relevance for computational liver surgery. We first summarize key facts of hepatic anatomy, physiology, and pathology relevant for hepatic surgery, followed by a description of the computational tools currently used in liver surgical planning. Then we present selected state-of-the-art computational liver models potentially useful to support liver surgery. Finally, we discuss the main challenges that will need to be addressed when developing advanced computational planning tools in the context of liver surgery.

Metabonomics screening of serum identifies pyroglutamate as a diagnostic biomarker for nonalcoholic steatohepatitis

OBJECTIVE

A key step in managing non-alcoholic fatty liver disease (NAFLD) is to differentiate nonalcoholic steatohepatitis (NASH) from simple steatosis (SS).

METHOD

Serum samples were collected from three groups: NASH patients (N=21), SS patients (N=38) and healthy controls (N=31). High performance liquid chromatography-mass spectrometry (HPLC-MS) was used to analyse the metabolic profile of the serum samples. The acquired data were processed by multivariate principal component analysis (PCA) and orthogonal partial least-squares-discriminant analysis (OPLS-DA) to identify novel metabolites. The potential biomarkers were quantitatively determined and their diagnostic power was further validated.

RESULTS

A total of 56 metabolites were capable of distinguishing NASH from SS samples based on the OPLS-DA model. Pyroglutamate was found to be the most promising factor in distinguishing the NASH from SS groups. With an optimal cut-off value of 4.82mmol/L, the sensitivity and specificity of the diagnosis of NASH were 72% and 85%, respectively. The area under the receiver operating characteristic (AUROC) of the pyroglutamate levels of NASH versus SS patients was more than those of tumor necrosis factor-α, adiponectin and interleukin-8.

CONCLUSION

These data suggest that pyroglutamate may be a new and useful biomarker for the diagnosis of NASH.