A multiscale modelling approach to assess the impact of metabolic zonation and microperfusion on the hepatic carbohydrate metabolism

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.

Pregnancy alters fatty acid metabolism, glucose regulation, and detoxification of the liver in synchrony with biomechanical property changes

Pregnancy places a metabolic burden on the body including the liver, which is responsible for ensuring adequate nutrition for the maternal and fetal systems. To gain a better understanding of liver adaptation, this study investigates metabolic shifts occurring in livers of pregnant rats. Metabolic capacities of the livers of pregnant and non-pregnant female Wistar rats were assessed using comprehensive metabolic models. Kinetic metabolic models were generated for each animal based on protein abundance data from proteomics analysis allowing for a subject-specific assessment of hepatic metabolic functions. Data are available via ProteomeXchange with identifier PXD050758. Additionally, tissue stiffness, viscosity, and water diffusion obtained from magnetic resonance imaging and elastography were correlated with metabolic capabilities to study the relationship between metabolic function and biophysical properties. Proteome profiling revealed differences in protein expression in the livers of pregnant and non-pregnant animals. Functional analysis showed significant variations in metabolic capacities. Livers of pregnant rats had reduced capacities in carbohydrate and fatty acid metabolism, along with altered urea synthesis. Additionally, there were associations between metabolic functions and biophysical properties highlighting potential links between changes in liver structure and metabolic capacities during pregnancy. In summary, our work reveals extensive hepatic metabolic changes in pregnant rats. The liver adapts its metabolic capacities to ensure whole-body metabolic homeostasis but may struggle to counteract nutritional challenges, such as hypoglycemia. The study, employing a personalized approach combining proteomics, kinetic modeling, and advanced imaging, sheds light on the intricate interplay between hepatic adaptations and medical imaging markers, providing a foundation for further investigations into the implications for maternal and fetal health.

Kiphynet: an online network simulation tool connecting cellular kinetics and physiological transport

INTRODUCTION

Human metabolism is sustained by functional networks that operate at diverse scales. Capturing local and global dynamics in the human body by hierarchically bridging multi-scale functional networks is a major challenge in physiological modeling.

OBJECTIVES

To develop an interactive, user-friendly web application that facilitates the simulation and visualization of advection-dispersion transport in three-dimensional (3D) microvascular networks, biochemical exchange, and metabolic reactions in the tissue layer surrounding the vasculature.

METHODS

To help modelers combine and simulate biochemical processes occurring at multiple scales, KiPhyNet deploys our discrete graph-based modeling framework that bridges functional networks existing at diverse scales. KiPhyNet is implemented in Python based on Apache web server using MATLAB as the simulator engine. KiPhyNet provides the functionality to assimilate multi-omics data from clinical and experimental studies as well as vascular data from imaging studies to investigate the role of structural changes in vascular topology on the functional response of the tissue.

RESULTS

With the network topology, its biophysical attributes, values of initial and boundary conditions, parameterized kinetic constants, biochemical species-specific transport properties such as diffusivity as inputs, a user can use our application to simulate and view the simulation results. The results of steady-state velocity and pressure fields and dynamic concentration fields can be interactively examined.

CONCLUSION

KiPhyNet provides barrier-free access to perform time-course simulation experiments by building multi-scale models of microvascular networks in physiology, using a discrete modeling framework. KiPhyNet is freely accessible at   http://pallab.cds.iisc.ac.in/kiphynet/ and the documentation is available at   https://deepamahm.github.io/kiphynet_docs/ .

Recreating metabolic interactions of the tumour microenvironment

Tumours are heterogeneous tissues containing diverse populations of cells and an abundant extracellular matrix (ECM). This tumour microenvironment prompts cancer cells to adapt their metabolism to survive and grow. Besides epigenetic factors, the metabolism of cancer cells is shaped by crosstalk with stromal cells and extracellular components. To date, most experimental models neglect the complexity of the tumour microenvironment and its relevance in regulating the dynamics of the metabolism in cancer. We discuss emerging strategies to model cellular and extracellular aspects of cancer metabolism. We highlight cancer models based on bioengineering, animal, and mathematical approaches to recreate cell-cell and cell-matrix interactions and patient-specific metabolism. Combining these approaches will improve our understanding of cancer metabolism and support the development of metabolism-targeting therapies.

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.

Strategic Replication of the Hepatic Zonation In Vitro Employing a Biomimetic Approach

The varied functions of the liver are dependent on the metabolic heterogeneity exhibited by the hepatocytes within the liver lobule spanning the porto-central axis. This complex phenomenon plays an important role in maintaining the physiological homeostasis of the liver. Standard in vitro culture models fail to mimic this spatial heterogeneity of hepatocytes, assuming a homogeneous population of cells, which leads to inaccurate translation of results. Here, we demonstrate the development of an in vitro model of hepatic zonation by mimicking the microarchitecture of the liver using a 3D printed mini bioreactor and decellularized liver matrix to provide the native microenvironmental cues. There was a differential expression of hypoxic and metabolic markers across the developed mini bioreactor, showing the establishment of gradients of oxygen, Wnt/β-catenin pathway, and other metabolic pathways. The model also showed the establishment of zone-dependent toxicity on treatment with acetaminophen. The developed model would thus be a promising avenue in the field of tissue engineering for understanding the liver physiology and pathophysiology and for drug screening to evaluate the potential of new pharmaceutical interventions.

Unveiling the power of microenvironment in liver regeneration: an in-depth overview

The liver serves as a vital regulatory hub for various physiological processes, including sugar, protein, and fat metabolism, coagulation regulation, immune system maintenance, hormone inactivation, urea metabolism, and water-electrolyte acid-base balance control. These functions rely on coordinated communication among different liver cell types, particularly within the liver's fundamental hepatic lobular structure. In the early stages of liver development, diverse liver cells differentiate from stem cells in a carefully orchestrated manner. Despite its susceptibility to damage, the liver possesses a remarkable regenerative capacity, with the hepatic lobule serving as a secure environment for cell division and proliferation during liver regeneration. This regenerative process depends on a complex microenvironment, involving liver resident cells, circulating cells, secreted cytokines, extracellular matrix, and biological forces. While hepatocytes proliferate under varying injury conditions, their sources may vary. It is well-established that hepatocytes with regenerative potential are distributed throughout the hepatic lobules. However, a comprehensive spatiotemporal model of liver regeneration remains elusive, despite recent advancements in genomics, lineage tracing, and microscopic imaging. This review summarizes the spatial distribution of cell gene expression within the regenerative microenvironment and its impact on liver regeneration patterns. It offers valuable insights into understanding the complex process of liver regeneration.

On the relationship between metabolic capacities and in vivo viscoelastic properties of the liver

The liver is the central metabolic organ. It constantly adapts its metabolic capacity to current physiological requirements. However, the relationship between tissue structure and hepatic function is incompletely understood; this results in a lack of diagnostic markers in medical imaging that can provide information about the liver's metabolic capacity. Therefore, using normal rabbit livers, we combined magnetic resonance elastography (MRE) with proteomics-based kinetic modeling of central liver metabolism to investigate the potential role of MRE for predicting the liver's metabolic function in vivo. Nineteen New Zealand white rabbits were investigated by multifrequency MRE and positron emission tomography (PET). This yielded maps of shear wave speed (SWS), penetration rate (PR) and standardized uptake value (SUV). Proteomic analysis was performed after the scans. Hepatic metabolic functions were assessed on the basis of the HEPATOKIN1 model in combination with a model of hepatic lipid-droplet metabolism using liquid chromatography-mass spectrometry. Our results showed marked differences between individual livers in both metabolic functions and stiffness properties, though not in SUV. When livers were divided into 'stiff' and 'soft' subgroups (cutoff SWS = 1.6 m/s), stiff livers showed a lower capacity for triacylglycerol storage, while at the same time showing an increased capacity for gluconeogenesis and cholesterol synthesis. Furthermore, SWS was correlated with gluconeogenesis and PR with urea production and glutamine exchange. In conclusion, our study indicates a close relationship between the viscoelastic properties of the liver and metabolic function. This could be used in future studies to predict non-invasively the functional reserve capacity of the liver in patients.

Transitional phase Gd-EOB-DTPA-enhanced MRI: Visual assessment of hepatic function and fibrosis based on uptake rate of Gd-EOB-DTPA

PURPOSE

The aim of the study is to investigate the feasibility for hepatic function and fibrosis visual assessment using transitional phase imaging based on the uptake process of Gd-EOB-DTPA.

MATERIALS AND METHODS

We retrospectively selected 105 consecutive patients who underwent Gd-EOB-DTPA enhanced MRI examination at 1.5 Tesla for intrahepatic lesion evaluation from June 2020 to June 2022. Data were classified into two groups defined by the signal intensity (SI) difference in the hepatic vein against that of the hepatic parenchyma at transitional phase as follows: High and Iso-SI group: hepatic vein SI equal to or greater than the hepatic parenchymal SI; and Low-SI group: hepatic vein SI lower than hepatic parenchymal SI. We evaluated whether significant differences in ALBI score, FIB-4, APRI and LSR of hepatobiliary phase between two groups. We measured cut-off values between two groups in all items according to receiver operating characteristic analysis. Furthermore, the inter-reader reproducibility of the visual assessment on transitional phase images between two readers was evaluated using the ICC.

RESULTS

The visual assessment results were as follows: High and Iso- and Low-SI groups included 48, 57, patients, respectively. Significant differences were observed in ALBI score, FIB-4, APRI and LSR between two groups. The cut-off values of ALBI score, FIB-4, APRI and LSR were -2.69, 2.28, 0.49 and 2.15. ICC of transitional phase image visual assessment between two readers was 0.86.

CONCLUSIONS

Hepatic function and fibrosis might be assessed by visual assessment of transitional phase images in Gd-EOB-DTPA enhanced MRI.

Simulation of the crosstalk between glucose and acetaminophen metabolism in a liver zonation model

The liver metabolizes a variety of substances that sometimes interact and regulate each other. The modeling of a single cell or a single metabolic pathway does not represent the complexity of the organ, including metabolic zonation (heterogeneity of functions) along with liver sinusoids. Here, we integrated multiple metabolic pathways into a single numerical liver zonation model, including drug and glucose metabolism. The model simulated the time-course of metabolite concentrations by the combination of dynamic simulation and metabolic flux analysis and successfully reproduced metabolic zonation and localized hepatotoxicity induced by acetaminophen (APAP). Drug metabolism was affected by nutritional status as the glucuronidation reaction rate changed. Moreover, sensitivity analysis suggested that the reported metabolic characteristics of obese adults and healthy infants in glucose metabolism could be associated with the metabolic features of those in drug metabolism. High activities of phosphoenolpyruvate carboxykinase (PEPCK) and glucose-6-phosphate phosphatase in obese adults led to increased APAP oxidation by cytochrome P450 2E1. In contrast, the high activity of glycogen synthase and low activities of PEPCK and glycogen phosphorylase in healthy infants led to low glucuronidation and high sulfation rates of APAP. In summary, this model showed the effects of glucose metabolism on drug metabolism by integrating multiple pathways into a single liver metabolic zonation model.

LKB1 acts as a critical brake for the glucagon-mediated fasting response

As important as the fasting response is for survival, an inability to shut it down once nutrients become available can lead to exacerbated disease and severe wasting. The liver is central to transitions between feeding and fasting states, with glucagon being a key initiator of the hepatic fasting response. However, the precise mechanisms controlling fasting are not well defined. One potential mediator of these transitions is liver kinase B1 (LKB1), given its role in nutrient sensing. Here, we show LKB1 knockout mice have a severe wasting and prolonged fasting phenotype despite increased food intake. By applying RNA sequencing and intravital microscopy, we show that loss of LKB1 leads to a dramatic reprogramming of the hepatic lobule through robust up-regulation of periportal genes and functions. This is likely mediated through the opposing effect that LKB1 has on glucagon pathways and gene expression. Conclusion: Our findings show that LKB1 acts as a brake to the glucagon-mediated fasting response, resulting in "periportalization" of the hepatic lobule and whole-body metabolic inefficiency. These findings reveal a mechanism by which hepatic metabolic compartmentalization is regulated by nutrient-sensing.

A Graph-Based Framework for Multiscale Modeling of Physiological Transport

Trillions of chemical reactions occur in the human body every second, where the generated products are not only consumed locally but also transported to various locations in a systematic manner to sustain homeostasis. Current solutions to model these biological phenomena are restricted in computability and scalability due to the use of continuum approaches in which it is practically impossible to encapsulate the complexity of the physiological processes occurring at diverse scales. Here, we present a discrete modeling framework defined on an interacting graph that offers the flexibility to model multiscale systems by translating the physical space into a metamodel. We discretize the graph-based metamodel into functional units composed of well-mixed volumes with vascular and cellular subdomains; the operators defined over these volumes define the transport dynamics. We predict glucose drift governed by advective-dispersive transport in the vascular subdomains of an islet vasculature and cross-validate the flow and concentration fields with finite-element-based COMSOL simulations. Vascular and cellular subdomains are coupled to model the nutrient exchange occurring in response to the gradient arising out of reaction and perfusion dynamics. The application of our framework for modeling biologically relevant test systems shows how our approach can assimilate both multi-omics data from in vitro-in vivo studies and vascular topology from imaging studies for examining the structure-function relationship of complex vasculatures. The framework can advance simulation of whole-body networks at user-defined levels and is expected to find major use in personalized medicine and drug discovery.

Role and Regulation of Wnt/β-Catenin in Hepatic Perivenous Zonation and Physiological Homeostasis

Metabolic heterogeneity or functional zonation is a key characteristic of the liver that allows different metabolic pathways to be spatially regulated within the hepatic system and together contribute to whole body homeostasis. These metabolic pathways are segregated along the portocentral axis of the liver lobule into three hepatic zones: periportal, intermediate or midzonal, and perivenous. The liver performs complementary or opposing metabolic functions within different hepatic zones while synergistic functions are regulated by overlapping zones, thereby maintaining the overall physiological stability. The Wnt/β-catenin signaling pathway is well known for its role in liver growth, development, and regeneration. In addition, the Wnt/β-catenin pathway plays a fundamental and dominant role in hepatic zonation and signals to orchestrate various functions of liver metabolism and pathophysiology. The β-catenin protein is the central player in the Wnt/β-catenin signaling cascade, and its activation is crucial for metabolic patterning of the liver. However, dysregulation of Wnt/β-catenin signaling is also implicated in different liver pathologies, including those associated with metabolic syndrome. β-Catenin is preferentially localized in the central region of the hepatic lobule surrounding the central vein and regulates multiple functions of this region. This review outlines the role of Wnt/β-catenin signaling pathway in controlling the different metabolic processes surrounding the central vein and its relation to liver homeostasis and dysfunction.

CARDIOKIN1: Computational Assessment of Myocardial Metabolic Capability in Healthy Controls and Patients With Valve Diseases

Supplemental Digital Content is available in the text.

Background:

Many heart diseases can result in reduced pumping capacity of the heart muscle. A mismatch between ATP demand and ATP production of cardiomyocytes is one of the possible causes. Assessment of the relation between myocardial ATP production (MV_ATP) and cardiac workload is important for better understanding disease development and choice of nutritional or pharmacologic treatment strategies. Because there is no method for measuring MV_ATP in vivo, the use of physiology-based metabolic models in conjunction with protein abundance data is an attractive approach.

METHOD:

We developed a comprehensive kinetic model of cardiac energy metabolism (CARDIOKIN1) that recapitulates numerous experimental findings on cardiac metabolism obtained with isolated cardiomyocytes, perfused animal hearts, and in vivo studies with humans. We used the model to assess the energy status of the left ventricle of healthy participants and patients with aortic stenosis and mitral valve insufficiency. Maximal enzyme activities were individually scaled by means of protein abundances in left ventricle tissue samples. The energy status of the left ventricle was quantified by the ATP consumption at rest (MV_ATP[rest]), at maximal workload (MV_ATP[max]), and by the myocardial ATP production reserve, representing the span between MV_ATP(rest) and MV_ATP(max).

Results:

Compared with controls, in both groups of patients, MV_ATP(rest) was increased and MV_ATP(max) was decreased, resulting in a decreased myocardial ATP production reserve, although all patients had preserved ejection fraction. The variance of the energetic status was high, ranging from decreased to normal values. In both patient groups, the energetic status was tightly associated with mechanic energy demand. A decrease of MV_ATP(max) was associated with a decrease of the cardiac output, indicating that cardiac functionality and energetic performance of the ventricle are closely coupled.

Conclusions:

Our analysis suggests that the ATP-producing capacity of the left ventricle of patients with valvular dysfunction is generally diminished and correlates positively with mechanical energy demand and cardiac output. However, large differences exist in the energetic state of the myocardium even in patients with similar clinical or image-based markers of hypertrophy and pump function.

Registration:

URL: https://www.clinicaltrials.gov; Unique identifiers: NCT03172338 and NCT04068740.

Diagnostic accuracy of the attenuation value in abdominal contrast enhanced dynamic multi-detector-row computed tomography for esophageal varices in patients with liver cirrhosis

Purpose

To investigate whether the attenuation value obtained by subtracting the CT value obtained from abdominal dynamic contrast enhanced (ADCE)-MDCT imaging of the equilibrium phase from the value obtained from that of the portal phase in hepatic parenchyma is useful in distinguishing normal liver from liver cirrhosis (LC) and to predict the development of esophageal varices (EVs) in patients with LC.

Materials and methods

We assigned 72 outpatients to group A (n = 45; normal liver) and group B (n = 27; LC), who underwent ADCE-MDCT. The attenuation value and CT value of the hepatic parenchymal portal and equilibrium phase were compared, and the correlation between attenuation value and biomarkers (ALB, T-bil, PLT, FIB-4, APRI, and AAR) was investigated. Furthermore, we investigated differences in the attenuation value, FIB-4, APRI, and AAR between the two subgroups of group B [without EVs (group a) and with EVs (group b)]. We performed receiver operating characteristic (ROC) analysis of the attenuation value, FIB-4, APRI, and, AAR for subgroup a vs b and evaluated the diagnostic accuracy.

Results

Significant differences were observed between groups A and B in all items. The attenuation value correlated with ALB, T-bil, PLT, FIB-4, and APRI. Only attenuation value showed a significant difference between groups a and b. The best cut-off attenuation value, FIB-4, APRI, and AAR for predicting EVs, according to ROC analysis was 13.4 HU, 6.8, 1.9, and 1.5.

Conclusions

Attenuation value can be useful to quantitatively classify normal liver and LC and to predict EVs in patients with LC.

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.

Functional Consequences of Metabolic Zonation in Murine Livers: Insights for an Old Story

BACKGROUND AND AIMS

Zone-dependent differences in expression of metabolic enzymes along the portocentral axis of the acinus are a long-known feature of liver metabolism. A prominent example is the preferential localization of the enzyme, glutamine synthetase, in pericentral hepatocytes, where it converts potentially toxic ammonia to the valuable amino acid, glutamine. However, with the exception of a few key regulatory enzymes, a comprehensive and quantitative assessment of zonal differences in the abundance of metabolic enzymes and, much more important, an estimation of the associated functional differences between portal and central hepatocytes is missing thus far.

APPROACH AND RESULTS

We addressed this problem by establishing a method for the separation of periportal and pericentral hepatocytes that yields sufficiently pure fractions of both cell populations. Quantitative shotgun proteomics identified hundreds of differentially expressed enzymes in the two cell populations. We used zone-specific proteomics data for scaling of the maximal activities to generate portal and central instantiations of a comprehensive kinetic model of central hepatic metabolism (Hepatokin1).

CONCLUSIONS

The model simulations revealed significant portal-to-central differences in almost all metabolic pathways involving carbohydrates, fatty acids, amino acids, and detoxification.

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

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.

Spatial sorting enables comprehensive characterization of liver zonation

The mammalian liver is composed of repeating hexagonal units termed lobules. Spatially resolved single-cell transcriptomics revealed that about half of hepatocyte genes are differentially expressed across the lobule, yet technical limitations impeded reconstructing similar global spatial maps of other hepatocyte features. Here, we show how zonated surface markers can be used to sort hepatocytes from defined lobule zones with high spatial resolution. We apply transcriptomics, miRNA array measurements and mass spectrometry proteomics to reconstruct spatial atlases of multiple zonated features. We demonstrate that protein zonation largely overlaps with mRNA zonation, with the periportal HNF4α as an exception. We identify zonation of miRNAs such as miR-122, and inverse zonation of miRNAs and their hepatocyte target genes, highlighting potential regulation of protein levels through zonated mRNA degradation. Among the targets we find the pericentral Wnt receptors Fzd7 and Fzd8 and the periportal Wnt inhibitors Tcf7l1 and Ctnnbip1. Our approach facilitates reconstructing spatial atlases of multiple cellular features in the liver and other structured tissues.

Spatial heterogeneity in the mammalian liver

Hepatocytes operate in highly structured repeating anatomical units termed liver lobules. Blood flow along the lobule radial axis creates gradients of oxygen, nutrients and hormones, which, together with morphogenetic fields, give rise to a highly variable microenvironment. In line with this spatial variability, key liver functions are expressed non-uniformly across the lobules, a phenomenon termed zonation. Technologies based on single-cell transcriptomics have constructed a global spatial map of hepatocyte gene expression in mice revealing that ~50% of hepatocyte genes are expressed in a zonated manner. This broad spatial heterogeneity suggests that hepatocytes in different lobule zones might have not only different gene expression profiles but also distinct epigenetic features, regenerative capacities, susceptibilities to damage and other functional aspects. Here, we present genomic approaches for studying liver zonation, describe the principles of liver zonation and discuss the intrinsic and extrinsic factors that dictate zonation patterns. We also explore the challenges and solutions for obtaining zonation maps of liver non-parenchymal cells. These approaches facilitate global characterization of liver function with high spatial resolution along physiological and pathological timescales.

Dynamic Metabolic Zonation of the Hepatic Glucose Metabolism Is Accomplished by Sinusoidal Plasma Gradients of Nutrients and Hormones

Being the central metabolic organ of vertebrates, the liver possesses the largest repertoire of metabolic enzymes among all tissues and organs. Almost all metabolic pathways are resident in the parenchymal cell, hepatocyte, but the pathway capacities may largely differ depending on the localization of hepatocytes within the liver acinus-a phenomenon that is commonly referred to as metabolic zonation. Metabolic zonation is rather dynamic since gene expression patterns of metabolic enzymes may change in response to nutrition, drugs, hormones and pathological states of the liver (e.g., fibrosis and inflammation). This fact has to be ultimately taken into account in mathematical models aiming at the prediction of metabolic liver functions in different physiological and pathological settings. Here we present a spatially resolved kinetic tissue model of hepatic glucose metabolism which includes zone-specific temporal changes of enzyme abundances which are driven by concentration gradients of nutrients, hormones and oxygen along the hepatic sinusoids. As key modulators of enzyme expression we included oxygen, glucose and the hormones insulin and glucagon which also control enzyme activities by cAMP-dependent reversible phosphorylation. Starting with an initially non-zonated model using plasma profiles under fed, fasted and diabetic conditions, zonal patterns of glycolytic and gluconeogenetic enzymes as well as glucose uptake and release rates are created as an emergent property. We show that mechanisms controlling the adaptation of enzyme abundances to varying external conditions necessarily lead to the zonation of hepatic carbohydrate metabolism. To the best of our knowledge, this is the first kinetic tissue model which takes into account in a semi-mechanistic way all relevant levels of enzyme regulation.