
Citation: Marcella Noorman, Richard Allen, Cynthia J. Musante, H. Thomas Banks. Analysis of compartments-in-series models of liver metabolism as partial differential equations: the effect of dispersion and number of compartments[J]. Mathematical Biosciences and Engineering, 2019, 16(3): 1082-1114. doi: 10.3934/mbe.2019052
[1] | Norbert Hungerbühler . Optimal control in pharmacokinetic drug administration. Mathematical Biosciences and Engineering, 2022, 19(5): 5312-5328. doi: 10.3934/mbe.2022249 |
[2] | H. Thomas Banks, W. Clayton Thompson, Cristina Peligero, Sandra Giest, Jordi Argilaguet, Andreas Meyerhans . A division-dependent compartmental model for computing cell numbers in CFSE-based lymphocyte proliferation assays. Mathematical Biosciences and Engineering, 2012, 9(4): 699-736. doi: 10.3934/mbe.2012.9.699 |
[3] | Christine K. Yang, Fred Brauer . Calculation of $R_0$ for age-of-infection models. Mathematical Biosciences and Engineering, 2008, 5(3): 585-599. doi: 10.3934/mbe.2008.5.585 |
[4] | Glenn Ledder . Incorporating mass vaccination into compartment models for infectious diseases. Mathematical Biosciences and Engineering, 2022, 19(9): 9457-9480. doi: 10.3934/mbe.2022440 |
[5] | Rebeccah E. Marsh, Jack A. Tuszyński, Michael Sawyer, Kenneth J. E. Vos . A model of competing saturable kinetic processes with application to the pharmacokinetics of the anticancer drug paclitaxel. Mathematical Biosciences and Engineering, 2011, 8(2): 325-354. doi: 10.3934/mbe.2011.8.325 |
[6] | Avner Friedman, Wenrui Hao . Mathematical modeling of liver fibrosis. Mathematical Biosciences and Engineering, 2017, 14(1): 143-164. doi: 10.3934/mbe.2017010 |
[7] | Sudhi Sharma, Victorita Dolean, Pierre Jolivet, Brandon Robinson, Jodi D. Edwards, Tetyana Kendzerska, Abhijit Sarkar . Scalable computational algorithms for geospatial COVID-19 spread using high performance computing. Mathematical Biosciences and Engineering, 2023, 20(8): 14634-14674. doi: 10.3934/mbe.2023655 |
[8] | Maryam Al-Yahyai, Fatma Al-Musalhi, Ibrahim Elmojtaba, Nasser Al-Salti . Mathematical analysis of a COVID-19 model with different types of quarantine and isolation. Mathematical Biosciences and Engineering, 2023, 20(1): 1344-1375. doi: 10.3934/mbe.2023061 |
[9] | Süleyman Cengizci, Aslıhan Dursun Cengizci, Ömür Uğur . A mathematical model for human-to-human transmission of COVID-19: a case study for Turkey's data. Mathematical Biosciences and Engineering, 2021, 18(6): 9787-9805. doi: 10.3934/mbe.2021480 |
[10] | Yi Jiang, Kristin M. Kurianski, Jane HyoJin Lee, Yanping Ma, Daniel Cicala, Glenn Ledder . Incorporating changeable attitudes toward vaccination into compartment models for infectious diseases. Mathematical Biosciences and Engineering, 2025, 22(2): 260-289. doi: 10.3934/mbe.2025011 |
Non-alcholic fatty liver disease (NAFLD) encompasses a spectrum of liver disease that begins as noninflammatory build-up of extra fat in liver cells (steatosis), can progress to non-alcoholic steatohepatitis (NASH) which is the occurrence of inflammation and liver damage in addition to the steatosis and eventually lead to cirrhosis [5,8,45]. Globally, an estimated 11–46% of people suffer from NAFLD [5,8,41]. Approximately 25–40% of patients with NAFLD will progress to NASH [8] with 40–50% of patients developing fibrosis and 20% developing cirrhosis [6]. The increasing prevalence of NAFLD has been linked to the rise of obesity [6] and is strongly associated with insulin resistance and type 2 diabetes mellitus [6,8,41,45]. Considered the hepatic manifestation of metabolic syndrome [8,17,45], NAFLD is currently the most common cause of chronic liver disease [5,8] and NASH is the second leading cause of cirrhosis in adults waiting for liver transplants in the United States [5,41,43], expected to become the first leading cause by 2030 [5,6]. It has also been suggested that a recent increase in the prevalence of hepatocellular carcinoma is caused by NAFLD [5,88].
As the prevalence of NAFLD continues to increase, the cause and pathogenesis of the disease remains an area of intensive study. Though our understanding of this disease has progressed substantially and continues to grow, there is much still to uncover. The pathogenesis of NAFLD is multifactorial and extremely complex, comprised of environmental and genetic factors, the specific contributions of which are not yet known [5]. Further, though the metabolic processes involved in steatosis have been investigated for decades, there is still much unknown about why lipids accumulate where they do in the liver.
In the continuing investigation of this disease, mathematical and computational models will serve an important role as they give us the ability to explore the complex relationships between various components in a low cost setting. In this paper, we will look at a model of the liver metabolism and examine certain aspects of it in order to inform future development of such models.
The liver is the main metabolic organ in the human body and, as such, plays a central role in the regulation of key metabolites including glucose and free fatty acid (FFA). Most of the metabolic processes that occur in the liver take place in the hepatocyte, the main cell type of the liver. These hepatocytes are tightly organized into tesselating columns embedded into liver lobules which make up each of the four lobes of the liver. The lobules have hexagonal cross sections, each corner of which harbors a portal triad which is comprised of a bile duct, hepatic artery (which supplies the hepatocytes with oxygenated blood), and a portal vein (which supplies the hepatocytes with nutrient filled blood). The blood exits the lobule through the central vein which is located at the center of the hexagonal structures. Numerous capillaries, called sinusoids, connect the periportal blood vessels of the portal triad with the central vein (see Figure 1 as adopted from [7]). These sinusoids are surrounded by the columns of hepatocytes and are the mechanism through which the hepatocytes have access to the blood provided by the portal triad.
Hepatocytes vary in enzyme expression and function depending upon where they lie along the sinusoid, a phenomena called zonation [40,48,49,50,53,54,55]. This heterogeneity among hepatocytes is, at least in part, the result of differing concentrations of metabolites and signaling molecules along the sinusoid. For example, as blood flows through the sinusoid, oxygen diffuses and is taken up into the hepatocytes resulting in higher levels of oxygen on the periportal end of the sinusoid (near the portal triad) and lower levels of oxygen on the pericentral end of the sinusoid (near the central vein). This oxygen gradient was discovered decades ago and has long since been considered a significant driving force of zonation [14,51,52,56,57,58,89]. Some studies have also indicated that the hormones insulin and glucagon play an important role in the regulation of zonated gene expression [50,53,72]. More recently, other signaling pathways have also been proposed and studied as potential mechanisms for zonation [25,46,87].
Metabolic zonation is not the only form of heterogeneity seen across the sinusoid. Most liver diseases show a varying amount of damage across the liver microstructure [50]. In NAFLD, studies have shown that, for adult patients, steatosis is typically most severe in pericentral cells [2,18,19,23,46,91]. However not all patients show such an aberrant accumulation of lipids on a particular part of the sinusoid [19,23]. Further, steatosis in pediatric patients with NAFLD tends to occur most severely in periportal cells [91]. The cause of this "zonated steatosis" in NAFLD is still not fully understood, however the heterogeneity of hepatocytes likely results in certain liver cells being more vulnerable to lipid build up than others [7]. Further, in patients who have progressed to NASH, inflammation and fibrosis tends to be most severe on the pericentral end of the sinusoid [1,18,23,46].
Although it is well known that lipids tend to accumulate heterogeneously across the sinusoid (both in NAFLD and in other liver diseases), the effect of metabolic zonation on such zonated steatosis has received limited attention experimentally. To a great extent, this is because investigating changes in individual regions of the sinusoid is extremely time consuming and complex [7]. Computational models of liver metabolism that are spatially distributed give us the ability to explore relationships between metabolic zonation and zonated steatosis without the high cost of experimentation. Since NAFLD is induced by a change in metabolism, that is, the build up of lipids in the liver, it is likely that metabolic zonation plays a major role in the pathology of this disease. Thus, the development of computational models of liver metabolism that include the spatial effects of zonation are of utmost importance in the continuing effort to understand and find treatment for NAFLD.
Computational models have been used previously to study the effect of zonation on detoxification processes such as the detoxification of ammonia [67], xenobiotics [39,82,84], and drugs [70,74], as well as more generally in pharmacokinetics [4,69]. However, few models currently exist that explicitly consider the heterogeneous expression of enzymes which defines metabolic zonation.
Existing mechanistic models of liver metabolism that include spatial distribution can be roughly divided into two categories depending on the types of equations used: those which use partial differential equations (PDEs) and those which use ordinary differential equations (ODEs). As will be shown later on, models from different categories may not always be as different from each other as may appear at first glance. The PDE models can further be divided into models derived using the theory of porous media and models that utilize advection or advection-dispersion equations.
The advantages of deriving a system using the theory of porous media is in the ability to more accurately describe the process of blood perfusion through the liver lobule. This is of particular interest since steatosis can affect hepatic hemodynamics [34,61,65,85]. In a series of publications [75,76,77], Ricken et al. develop a multi scale approach that does just this, coupling the blood perfusion through the liver lobule with a description of hepatic cell metabolism. While metabolic zonation is not explicitly included, they are still able to simulate some zonated steatosis.
Models using advection or advection-dispersion equations to represent the blood flow in the liver have been around for decades finding applications in many areas such as hepatic elimination, drug clearance, and capillary tissue exchange [11,12,29,69,78]. These models have generally been referred to as parallel-tube and distributed models (advective flow) [11,12,37] and dispersion models (advective-dispersive flow) [79]. A rather succinct but complete discussion pertaining to mass balance and mass transport involving compartmental modeling, local concentration gradients, advection and advection/dispersion in moving fluids can be found in Chapter 4 of [10]. In particular see Sections 4.4.1 and 4.4.2 of that reference along with a number of other excellent similar presentations and discussions in Chapter 5 of [80], Chapters 15–17 of [59] as well as the by now classical text of 68 [68].
Similar models to the dispersion model have been used in the past decade to study aspects of hepatic metabolism as well [21,24]. In 2006, 24 et al. [24] used advection-dispersion-reaction equations to represent the concentrations of metabolites in the sinusoid, coupled through metabolic transport to mass-balance equations for the metabolites in the tissue. Enzymatic zonation is included for central carbohydrate pathways, but not for lipid metabolism. In 2008, 21 et al. [21] used a model similar to the one given in [24] to illustrate how one may embed such a model of cellular liver metabolism into a Bayesian framework for parameter estimation purposes. However, zonation is excluded from the model and the system is only considered at steady state. Use of these types of PDEs to model the liver microstructure has appeared to lag in the past decade as compartmental models have started to show prevalence. An explanation for this may be due to the appearance of these PDEs as being more difficult to use in comparison to compartmental models which are composed of ODEs [4,42].
Compartmental models for the liver that include spatial distribution were originally inspired by "tanks-in-series" models from chemical engineering [42] and can be considered "compartments-in-series" models. These models are comprised of systems of ODEs derived by considering a series of compartments that are connected together, typically through blood flow. While early models of this type had each compartment represent both blood and tissue [42] so that there is only a single series of compartments, modern models tend to separate blood and tissue into different compartments with a transport between them [6,7,13,67,81]. These models have a series of hepatic blood compartments that are connected through blood flow along with a series of tissue compartments which are only connected to their corresponding hepatic blood compartment (i.e., there is no explicit connection between tissue compartments). These models often treat the series of compartments as being lined up along the portocentral axis of the liver sinusoid, and as such treat this as the repeating unit of the liver (rather than the hepatocyte). Recently, compartmental models that include spatial distribution have been applied more specifically to hepatic metabolism.
Schleicher et al. [81] used a three compartment model to simulate hepatic lipid metabolism. The hepatic blood flow is a basic advective flow that transports concentrations of metabolites in the hepatic blood from the periportal to the pericentral end of the sinusoid. This three compartment model is used to look at how the plasma oxygen gradient, plasma FA gradient, as well as the process for FA uptake affect the zonation of steatosis under a high-fat diet. Despite the fact that zonated enzyme expressions are not included in this model, model simulations do show a zonated steatosis.
A much more complex compartmental model was used by Berndt et al. [13] to study glucose metabolism and how it is affected by zonated enzyme expressions, metabolite and hormone gradients in the sinusoid, and blood perfusion. This model includes the blood in the sinusoid, accompanying space of Disse, as well as the adjacent layers of hepatocytes allowing the authors to take into account both morphological and systemic parameters. Each component (sinusoid, space of Disse, and hepatocytes) is partitioned into
Among the more advanced recent compartmental models for liver metabolism is the one presented by Ashworth et al. in [6,7]. This
As NAFLD becomes more and more prevalent, the ability to perform accurate in silico experiments is an integral part of uncovering important mechanisms and targets for these conditions. Compartments in series models are an attractive tool that allow one to investigate how metabolic zonation, hepatic blood gradients of metabolites and hormones, and even blood perfusion, affect liver metabolism at various levels of complexity. In particular, these models boast the ability to investigate various aspects of liver metabolism with sharp detail while keeping computational times reasonable. However, one may question if the simplification to a number of compartments rather than a continuum may result in the loss of information or misrepresentation of results. In what follows, we will examine a model similar to the one presented in [6,7] in order to begin attempting to answer this question. The two main goals of this paper are (ⅰ) to investigate the effect that the number of compartments has on model simulations and (ⅱ) to investigate the effect of dispersion on the system.
We present the models of the liver microstructure under consideration here. Our models closely resembles that of [6,7], with the main differences being how concentrations in the systemic blood are modeled as well as how we view the equations in the hepatic blood. We chose this model because it was similar to one used by one (MN) of us in an earlier computational effort during a summer project at Pfizer and because it seems sufficiently detailed to illustrate the ideas of our focus here: how the number of compartments and inclusion of dispersion in the flow affect simulation results. Even though it does not treat a number of important factors treated in other liver modeling efforts [39,75,76,77,84] such as hematocrit, transfer of oxygen between blood and hepatocyte, Hep to Hep transfer, varying sinusoidal velocities, etc., we believe the model is sufficient for our purposes here.
We do model the concentrations of a number of metabolites, hormones, and oxygen in the liver sinusoid. Equations are separated into those for the hepatocytes (referred to as hepatic equations), for hepatic blood (i.e., blood/plasma in the sinusoid), and systemic blood (i.e., blood/plasma in the rest of the body). While oxygen, hormones, and metabolites are all modeled in the hepatic and systemic blood, only metabolites are modeled in the hepatocyte. However, the concentration of hormones and oxygen in the hepatic blood do affect the metabolic processes that occur in the hepatocyte. A full list of all the concentrations modeled is included in Table 1 along with what kinds of equations they have and their abbreviated variable name that will be used to denote them throughout this paper. Figure 2 (see Figure 1 of [6]) gives a schematic of the overall structure of the model. We will use the following notation convention: variables in the hepatocyte will have a subscript
Eqns in H, HB, SB | Eqns in HB, SB | Eqns in H | |||
Glucose | | Insulin | | Glycogen | |
Triglycerides | | Glucagon | | Glucose-6-phosphate | |
Fatty acids | | Oxygen | | Glycerol-3-phosphate | |
Glycerol | | Acetyl-CoA | | ||
Lactate | | Inorganic phosphate | | ||
Guanosine Tri/Di-phosphate | | ||||
Uridine Tri/Di-phosphate | | ||||
Adenosine Tri/Di/Mono-phosphate | |
The equations for molecule concentrations in the hepatocytes are given in Table 3 below.
Molecule | Equation |
Glycogen | |
Glucose-6-phosphate | |
Glucose | |
Glycerol-3-phosphate | |
Lactate | |
Acetyl-CoA | |
Fatty acid | |
Triglyceride | |
Glycerol | |
Guanosine Tri-phosphate | |
Guanosine Di-phosphate | |
Uridine Tri-phosphate | |
Uridine Di-phosphate | |
Inorganic phosphate | |
Adenosine Tri-phosphate | |
Adenosine Di-phosphate | |
Adenosine Mono-phosphate |
In the model presented in [6,7], the rate of hepatic metabolic processes are modeled through the use of Hill functions. Rather than modeling every individual enzyme, sections of the metabolic pathways are often modeled by a single function, thus depending upon a number of intermediate enzymes rather than a single one. Each substrate and allosteric activator/inhibitor is represented by a Hill function and the rate of the process is calculated to be the product of these Hill functions. The hepatic metabolic processes included in the model are given in Tables 4–6 along with the function for the rate of that process.
Metabolic process | Abbrev. | Rate |
Glucokinase | | |
Glucose-6-phosphatase | | |
Glycogen Synthase | | |
Glycogen Phosphorylase | | |
Phosphoenolpyruvate Carboxykinase | | |
Fructose Biphosphatase | | |
Phosphofructokinase | | |
| ||
Pyruvate Kinase | |
Metabolic process | Abbrev. | Rate |
Pyruvate Oxidation | | |
| | |
Lipogenesis | | |
Lipolysis | | |
Triglyceride Synthesis | | |
Glycerol Kinase | |
Metabolic process | Abbrev. | Rate |
Adenosine Tri-phosphate Synthesis | | |
Uridine Diphosphate Kinase | | |
Guanosine Diphosphate Kinase | | |
Adenosine Kinase | | |
Additional Adenosine Tri-phosphate Use | | |
Control of Cellular Phosphate Levels | |
The transport of metabolites between plasma and hepatocyte are modeled similarly. For uni-directional transports, the rate is calculated using a Hill function where the plasma molecule is treated as the substrate and the hepatic molecule as the product. For bi-directional transports, the rate is calculated by a "Hill-like" function that depends on the difference in concentration of the molecule between the plasma and hepatocyte. Table 2 shows the rates of transport of metabolites between plasma and hepatocyte.
{Transport} | Rate |
Glycerol | |
Triglyceride | |
Fatty acid | |
Glucose (GLUT2) | |
Lactate | |
The regulation of fatty acid (FA) transport and uptake into the liver is a complex and often debated topic [62]. There is strong evidence for protein-mediated FA uptake with a number of membrane proteins being implicated [3,62]. In particular, in the liver such proteins include fatty acid translocase FAT/CD36 [20,62,63], liver fatty acid binding protein LFABP [66,73], and fatty acid transport proteins, specifically FATP2, FATP3, and FATP5 [3,15,16,32,33,35,62,64]. Passive diffusion also contributes to hepatic FA uptake, however this process is minimal in comparison to facilitated uptake [62].
In [6,7], the FFA transport is composed of a uni-directional, insulin dependent term and a bi-directional, insulin independent term. We omit the insulin dependent term here as, under physiological conditions, the hepatic FA transporters that appear to be the major players in hepatic FA uptake (LFABP, FATP2, FATP3, and FATP5) [3,16,66] are not known to be regulated by insulin. The greatest evidence for insulin regulated hepatic FA uptake is in FAT/CD36 mediated FA uptake.
FAT/CD36 is a fatty acid transporter that is regulated by muscle contraction and/or insulin and is an important protein implicated in FA uptake by skeletal muscle cells, cardiomyocytes, and adipocytes [20]. Under physiological conditions, FAT/CD36 is weakly expressed in hepatocytes suggesting FA uptake is largely FAT/CD36 independent [20]. However, there is evidence that hepatic mRNA and protein levels of FAT/CD36 at the plasma membrane of hepatocytes increases in patients with NAFLD [20,62,63].
Of the FATP family of transporters, the only one known to be regulated by insulin is FATP1 [3]. While FATP1 does have some expression in the liver, this transport is more significant for skeletal muscle, heart, and adipose tissue [3,16]. In fact, studies have shown that while mice lacking FATP1 show a decrease in insulin-stimulated long chain FA uptake in adipocytes and skeletal muscle, they actually have an increase in long chain FA uptake in the liver [90].
The zonation of the enzyme expressions are modeled through the metabolic processes (pathways) they are involved in by utilizing information about the enzymes involved in the process as well as the location of the hepatocyte along the sinusoid. Hepatocytes are assigned values
Each metabolic process
vMx=(1+zxkM)vMb | (2.1) |
where
It is in the representation and inclusion of metabolic processes in the systemic blood (referred to as the body compartment in [6,7]) that the model presented here truly differs from the model in [6,7] (excluding the addition of dispersion discussed in Section 2.3.2). These changes were made in order to improve biological detail as well as to simplify the equations here. Since the main focus of this model is on liver metabolism, we leave these equations as simple as possible while keeping the concentrations at realistic levels. The main issue we came across was that the meal inputs were too significant of a driver of the system and resulted in metabolite concentrations falling too quickly in the absence of these inputs. To address this, we adjust the equations for the concentrations of metabolites and hormones in the systemic blood. A summary of the most significant changes is given below. Plots showing selected metabolite and hormone levels in the systemic blood before and after these changes are given in Figure 3. Equations for the concentrations in the systemic blood are given in Table 7 (note,
Molecule} | Equation |
Glycerol | |
Triglyceride | |
Fatty acid | |
Glucose | |
Lactate | |
Insulin | |
Glucagon | |
Oxygen | |
In the original model, glucose uptake in the rest of the body is represented by a function
˜SgCup=v[gCSB]KM+[gCSB]. | (2.2) |
While there is insulin independent glucose uptake in the body such as in the brain, the uptake of glucose in other tissues (such as muscle and adipose) does depend upon insulin [22,83]. For this reason, we formulated a second term for the rate of insulin dependent glucose uptake in the body using Hill functions. Our resulting function for the rate of glucose uptake in the rest of the body is given by
SgCup=v1[gCSB]n1Kn1M1+[gCSB]n1+v2([InsSB]nIKnIMI+[InsSB]nI)([gCSB]n2Kn2M2+[gCSB]n2). | (2.3) |
The parameters for this function were estimated using a weighted least squares scheme with data from [28] and predictions were tested using data from [86] (see Figure 4). Finally, to get
We note that the data used to estimate the parameters for
With no adipose compartment in the model, there is also no triglyceride storage in adipose tissue. In [6], plasma triglycerides are treated as adipose triglycerides for lipolysis in adipose tissue. To account for the fact that the size of the plasma triglyceride pool is much smaller than those stored in adipose tissue, rate constants that had a slow dependence upon plasma glucagon and insulin concentration were used. In our simulations, however, this adjustment to the rate constants did not prevent the plasma triglyceride pool from emptying during fasting periods (i.e., periods without meal inputs). To alleviate this, we take the rate of lipolysis in adipose tissue to be driven by insulin action alone. Figure 5 gives schematic representations of how lipolysis in adipose tissue is represented in the model from [6] (Figure 5a) and in the model given here (Figure 5b).
Following the work of Periwal et al. [60,71], we model the rate of lipolysis in adipose tissue with the following system
d[FASB]dt=L(IA)−cf[FASB] | (2.4a) |
dIAdt=cI([InsSB]−IA−Ib) | (2.4b) |
L(IA)=l0+l21+(IA/IA2)α | (2.4c) |
where
Once the parameters for (2.4) were estimated, parameters for the transport as well as
In [6], the meal inputs are given by a high-powered sine function (i.e.,
The DE used to simulate meal ingestion is a two-compartment model where each compartment represents part of the gut and is given by
{dG1dt=−k1G1+αGinfs,G1(0)=αGindG2dt=k1G1−k2G2,G2(0)=0. | (2.5) |
The rate of change in plasma glucose from the ingestion of meals is then given by
In this paper, we will look at modeling the hepatic blood flow two different ways. One is modeling the blood flow as just a pure advective flow, which leads us to an advection-reaction (AR) equation for the concentrations of metabolites in the hepatic blood. The other is using Taylor's axial dispersion model which results in advection-dispersion-reaction (ADR) equations for the concentrations of metabolites in the hepatic blood. We perform simulations using both the AR and ADR equations and compare the results to determine whether or not the addition of dispersion in the blood flow, while more biologically accurate, has a significant effect upon the model solutions. The equations for both the AR system and the ADR system are given in Table 8. We note that we leave the equations for the hormones and oxygen the same for both models as they just give a basic gradient.
Mol. | AR Equation | ADR Equation |
| | |
| | |
| | |
| | |
| | |
| | |
| | |
| | |
In [6,7], the authors explain that they partition the hepatocytes around the sinusoid and blood flowing through the sinusoid into compartments depending upon where they lie on the sinusoid (see Figure 8). In this view, one might think of the model as multiple two compartment models that are connected by the blood flow through the hepatic blood compartments. This blood flow is modeled through a basic rate-in minus rate-out term. That is, supposing we've partitioned the hepatic blood into
d[CHB]idt=vbf([CHB]i−1−[CHB]i)−γHTCi | (2.6) |
for
vbf=QHn | (2.7) |
per second, where
d[CHB]idt=QH1/n([CHB]i−1−[CHB]i)−γHTCi | (2.8) |
and observe that
∂[CHB]i∂x≈QH1/n([CHB]i−1−[CHB]i) | (2.9) |
for
∂[CHB]∂t=QH∂[CHB]∂x−γHTC([CHB],[CH]). | (2.10) |
We note that
In [6,7], when viewing the model as a compartments-in-series model, the authors treat the blood flow in the systemic blood (referred to in [6,7] as the body compartment) similarly. If
d[CSB]dt=QHγB([CHB]n−[CSB])+S | (2.11) |
where
∂[CHB]∂t=QH∂[CHB]∂x−γHTC([CHB],[CH]),0<x<1,t>0, | (2.12) |
[CHB](0,t)=[CSB](t),t>0, | (2.13) |
[CHB](x,0)=fHB,0<x<1, | (2.14) |
where
d[CSB]dt=QBγB([CHB](1,t)−[CSB](t))+S, | (2.15) |
with initial condition
While an AR equation such as the one given in Section 2.3.1 is a good model for simple flows, the actual dynamics of blood flow through a capillary are much more complex than just a pure advective flow. Here, we include the effect of the dispersion that occurs due to the flow by modeling the blood flow using Taylor's axial dispersion model [36]. This model has been used previously for blood flow through the liver sinusoid in other works such as [21,24]. Here, we add a dispersion or diffusion term [10] to the AR model of Section 2.3.1 in order to get an advection-dispersion-reaction (ADR) equation in order to compare the two and thus determine whether or not a dispersion term makes a significant difference for the model presented here (one of the primary goals of this presentation).
The addition of a dispersion term, and thus the addition of a second spatial derivative, requires boundary conditions at both the periportal and pericentral ends of the sinusoid in order to remain well-posed (in contrast, the AR equation only requires boundary conditions at one end). To accommodate this, we assign homogeneous Neumann conditions at the pericentral end to indicate that the concentration of the metabolite leaving the sinusoid is the same as the concentration entering the systemic blood. The concentration of metabolite
∂[CHB]∂t=QH∂[CHB]∂x−DC∂2[CHB]∂x2−γHTC([CHB],[CH]),0<x<1,t>0, | (2.16) |
[CHB](0,t)=[CSB](t),t>0, | (2.17) |
∂[CHB]∂x(1,t)=0,t>0, | (2.18) |
[CHB](x,0)=fHB,0<x<1, | (2.19) |
where
Determining the appropriate values for
If
vbf=average rate of blood flow in the livervolume of blood compartment. | (2.20) |
In general, the rate of blood flow in the liver is given by the volume of blood in the liver divided by the time it takes the blood to go through the liver. Since the volume of a blood compartment should be the volume of blood in the liver divided by
vbf=ntime it takes blood to go through the liver. | (2.21) |
Thus, from (2.7), we have that, under this scaling, the time it takes blood to go through the liver is
In order to have an unscaled advection equation, we would need to know the velocity of the blood in the liver. Let
v=distance traveled through livertime it takes blood to go through the liver=L1/QH | (2.22) |
or
v=QHL. | (2.23) |
Now, while our model is of the liver microstructure (that is of a liver sinusoid and the surrounding cells), the concentrations of metabolites our model gives is actually representative of the accumulation from all sinusoids in the liver. On the other hand, in [24], the authors model the production from a single sinusoid using Taylor's axial dispersion and using physical measurements for the length of the sinusoid, velocity of blood, etc. Assuming then that our model gives us the concentration from a single sinusoid times the number of sinusoids in the liver, one can obtain that we should have
L=vQH=LsinQHτ=1.23mm. | (2.24) |
Now that we have our scaling factor
DiC=1.013×10−4×(MW)−0.46cm2s | (2.25) |
to compute the diffusion coefficient
DsinC=DiC+(vsinds2)248DiC | (2.26) |
to calculate
DC=DsCL2. | (2.27) |
The values used here are the same as those used in Chalhoub, et al. We do not believe simulating shorter sinusoids (e.g. 0.25 mm) with higher transit times (e.g., 1 second) will significantly affect the results. Doing such a simulation would cause the scaling factor L to be 1.67 rather than 1.23 and hence means that we would be dividing dispersion coefficients by 2.7889 rather than 1.5129. Since the unscaled dispersion coefficients themselves are already on the order of
Discretization in space is carried out using a first-order upwind difference for the advection term and a second-order central difference for the dispersion term. The resulting ODE's are then solved using ode15s from the Matlab ODE suite.
First, we look at how varying the mesh size
Largest | Largest | Largest | |||
Variable | change | Variable | change | Variable | change |
| |||||
| |||||
| |||||
| |||||
| |||||
| |||||
| |||||
| |||||
The differences seen due to the mesh size (or number of compartments) can be explained by the addition of numerical dispersion in the hepatic blood equations through the first-order upwind advection discretization. It can be shown that while the first-order upwind scheme is first-order accurate to the advection equation
ut+aux=0, | (3.1) |
it is actually second-order accurate to the modified advection-dispersion equation
˜ut+a˜ux=12ah˜uxx, | (3.2) |
where
Secondly, we want to investigate how the intentional inclusion of dispersion affects the system. That is, we wish to compare the AR simulation results with the ADR simulation results. In an attempt to limit the effect of numerical dispersion from the advection term while keeping computation time fairly short, we use
Largest | Largest | Largest | |||
Variable | change | Variable | change | Variable | change |
| |||||
| |||||
| |||||
| |||||
| |||||
| |||||
| |||||
| |||||
Figure 10 depicts the metabolite levels across the sinusoid at the time for which the solutions to the AR and ADR versions of the model have the largest difference for each particular metabolite (as measured by the discrete
Largest | Largest | Largest | |||
Variable | change | Variable | change | Variable | change |
Lastly, under this discretization, it appears that including the dispersion term actually increases computation time exponentially. Simulations were run using various mesh sizes (from
tADR=0.1720e0.3119tAR | (3.3) |
where
In this paper, we investigate the effect the number of compartments used in a model of the liver microstructure (similar to the one presented in [6,7]) has on the simulated concentrations as well as how including dispersion affects the system. The model of the liver microstructure from [6,7] was adjusted in order to keep metabolic levels realistic when the system is not being driven by meal inputs. This adjustment included, most influentially, the addition of insulin dependent glucose uptake in the body excluding the liver, rate of lipolysis in adipose tissue depending solely upon insulin action, and representation of meal inputs as the solution to a two-compartment model of ingestion in the gut. The model in [6,7] is presented as a compartments-in-series model, however we present it here as a system of advection-reaction equations that give the dynamics of concentrations in the hepatic blood, coupled to mass-balance equations that provide dynamics for the concentrations in the hepatic tissue. The in-flow boundary condition for the hepatic blood equations is given by the solution to an ODE that yields a basic representation of metabolism in the rest of the body which depends upon the out-flow concentrations, hence creating a feedback loop. The reason for this difference in representation is the fact that the equations for concentrations in the hepatic blood in the compartments-in-series model given in [6,7] are equivalent to advection-reaction equations that have been spatially discretized using a upwind forward difference and a mesh size of
Numerical simulations were carried out using the upwind forward difference discretization for the advection term in order to investigate how changing the number of compartments
We also address how intentionally included dispersion affects the system. This is of interest because it is biologically more accurate and has been included in other models of the liver microstructure. Most notably, 13 et al. [13] include a dispersive affect in the blood flow of their compartments-in-series model. When looked at closely, one sees that this dispersion is in fact modeled through a second-order central difference discretization for a standard dispersion term (that is, a discretization for
In the past, arguments have been made for compartments-in-series models over PDE models for hepatic elimination often citing the mathematical simplicity of ODEs compared to PDEs [4,42]. However, many of these compartments-in-series models are equivalent to a spatially discretized PDE, especially in how they've been used thus far to model liver metabolism. Thus, considering these models strictly as ODEs is potentially limiting the information and insights these systems may be able to give us. In particular, if our goal is to understand why and how certain spatially dependent processes occur, such as metabolic zonation and the phenomena of zonated steatosis in liver disease, then it makes sense to consider systems of PDEs which give us more options when it comes to discretizing spatially in numerical schemes. Further, PDEs lends themselves more appropriately in this problem to theoretical analyses that can unveil properties of the system such as the functional framework and stability. Compartments-in-series models are very useful computational tools and considering these systems from the perspective of a PDE does not diminish that. Rather, it provides more options and tools with which to analyze the system, both numerically and theoretically.
This research was supported in part by the Air Force Office of Scientific Research (HTB, MN) under grant number AFOSR FA9550-18-1-0457 and in part by a CRSC/Lord Corporation Fellowship (MN).
All authors declare no conflicts of interest in this paper.
[1] | M. F. Abdelmalek and A. M. Diehl, Nonalcoholic fatty liver disease as a complication of insulin resistance, Med. Clin. North Am., 91 (2007), 1125-1149. |
[2] | L. A. Adams, P. Angulo and K. D. Lindor, Nonalcoholic fatty liver disease, CMAJ, 172 (2005), 899-905. |
[3] | C. M. Anderson and A. Stahl, SLC27 fatty acid transport proteins, Mol. Aspects Med., 34 (2013), 516-528. |
[4] | Y. G. Anissimov and M. S. Roberts, A compartmental model of hepatic disposition kinetics: 1. Model development and application to linear kinetics, J. Pharmacokinet. Pharmacodyn., 29 (2002), 131-156. |
[5] | J. P. Arab, M. Arrese and M. Trauner, Recent insights into the pathogenesis of nonalcoholic fatty liver disease, Annu. Rev. Pathol. Mech. Dis., 13 (2018), 321-350. |
[6] | W. B. Ashworth, N. A. Davis and I. D. L. Bogle, A computational model of hepatic energy metabolism: understanding zonated damage and steatosis in NAFLD, PLoS Comput. Biol., 12 (2016), e1005105. |
[7] | W. B. Ashworth, C. Perez-Galvan, N. A. Davies and I. D. L. Bogle, Liver function as an engineering system, AIChE J., 62 (2016), 3285-3297. |
[8] | B. A. Banini and A. J. Sanyal, Nonalcoholic fatty liver disease: epidemiology, pathogenesis, natural history, diagnosis, and current treatment options, Clin. Med. Insights Ther., 8 (2016), 75-84. |
[9] | H. T. Banks, Modeling and Control in the Biomedical Sciences, Lecture Notes in Biomathematics, Vol 6, Springer-Verlag, New York, 1975. |
[10] | H.T. Banks and H.T. Tran, Mathematical and Experimental Modling of Physical and Biological Processes, Chapman & Hall/CRC Press, Taylor & Francis Group, Boca Raton, 2009. |
[11] | L. Bass, S. Keiding, K. Winkler and N. Tygstrup, Enzymatic elimination of substrates flowing through the intact liver, J. Theor. Biol., 61 (1976), 393-409. |
[12] | L. Bass, P. Robinson and A. J. Bracken, Hepatic elimination of flowing substrates: the distributed model, J. Theor. Biol., 72 (1978), 161-184. |
[13] | N. Berndt, M. S. Horger, S. Bulik and H. G. Holzhütter, A multiscale modelling approach to assess the impact of metabolic zonation and microperfusion on the hepatic carbohydrate metabolism, PLoS Comput. Biol., 14 (2018), e1006005. |
[14] | S. N. Bhatia, M. Toner, B. D. Foy, A. Rotem, K. M. O'Neil, R. G. Tompkins and M. L. Yarmush, Zonal liver cell heterogeneity: effects of oxygen on metabolic functions of hepatocytes, Cell. Eng., 1 (1996), 125-135. |
[15] | P. N. Black, C. Ahowesso, D. Montefusco, N. Saini and C. C. DiRusso, Fatty acid transport proteins: targeting FATP2 as a gatekeeper involved in the transport of exogenous fatty acids, Medchemcomm, 7 (2016), 612-622. |
[16] | P. N. Black, A. Sandoval, E. Arias-Barrau and C. C. DiRusso, Targeting the fatty acid transport proteins (FATP) to understand the mechanisms linking fatty acid transport to metabolism, Immunol. Endocr. Metab. Agents Med. Chem., 9 (2009), 11-17. |
[17] | E. M. Brunt, Pathology of nonalcoholic fatty liver disease, Nat. Rev. Gastroenterol. Hepatol., 7 (2010), 195-203. |
[18] | E. M. Brunt, Pathology of fatty liver disease, Mod. Pathol., 20 (2007), S40-S48. |
[19] | E. M. Brunt, C. G. Janney, A. M. Di Bisceglie, B. A. Neuschwander-Tetri and B. R. Bacon, Nonalcoholic steatohepatitis: a proposal for grading and staging the histological lesions, Am. J. Gastroenterol., 94 (1999), 2467-2474. |
[20] | X. Buqué, A. Cano, M. E. Miquilena-Colina, C. García-Monzón, B. Ochoa and P. Aspichueta, High insulin levels are required for FAT/CD36 plasma membrane translocation and enhanced fatty acid uptake in obese Zucker rat hepatocytes, Am. J. Physiol. Endocrinol. Metab., 303 (2012), E504-E514. |
[21] | D. Calvetti, A. Kuceyeski and E. Somersalo, Sampling-based analysis of a spatially distributed model for liver metabolism at steady state, Multiscale Model. Simul., 7 (2008), 407-431. |
[22] | G. D. Cartee, Mechanisms for greater insulin-stimulated glucose uptake in normal and insulin-resistant skeletal muscle after acute exercise, Am. J. Physiol. Endocrinol. Metab., 309 (2015), E949-E959. |
[23] | N. Chalasani, L. Wilson, D. E. Kleiner, O. W. Cummings, E. M. Brunt and A. Ünalp, Relationship of steatosis grade and zonal location to histological features of steatohepatitis in adult patients with non-alcoholic fatty liver disease, J. Hepatol., 48 (2008), 829-834. |
[24] | E. Chalhoub, L. Xie, V. Balasubramanian, J. Kim and J. Belovich, A distributed model of carbohydrate transport and metabolism in the liver during rest and high-intensity exercise, Ann. Biomed. Eng., 35 (2007), 474-491. |
[25] | M. Colletti, C. Cicchini, A. Conigliaro, L. Santangelo, T. Alonzi, E. Pasquini, M. Tripodi and L. Amicone, Convergence of Wnt signaling on the HNF4α-driven transcription in controlling liver zonation, Gastroenterology, 137 (2009), 660-672. |
[26] | S. W. Coppack, R. M. Fisher, G. F. Gibbons, S. M. Humphreys, M. J. McDonough, J. L. Potts and K. N. Frayn, Postprandial substrate deposition in human forearm and adipose tissues in vivo, Clin. Sci. (Lond.), 79 (1990), 339-348. |
[27] | M. E. Daly, C. Vale, M. Walker, A. Littlefield, K. G. Alberti and J. C. Mathers, Acute effects of insulin sensitivity and diurnal metabolic profiles of a high-sucrose compared with a high-starch diet, Am. J. Clin. Nutr., 67 (1998), 1186-1196. |
[28] | R. A. DeFronzo and E. Ferrannini, Influence of plasma glucose and insulin concentration on plasma glucose clearance in man, Diabetes, 31 (1982), 683-688. |
[29] | A. Deussen and J. B. Bassingthwaighte, Modeling [15O] oxygen tracer data for estimating oxygen consumption, Am. J. Physiol., 270 (1996), H1115-H1130. |
[30] | Gerda de Vries, Thomas Hillen, Mark Lewis, Johannes Muller and Birgitt Schonfisch, A Course in Mathematical Biology: Quantitative Modelling with Mathematical and Computational Methods, SIAM, Philadephia, 2006. |
[31] | G. Dimitriadis, P. Mitrou, V. Lambadiari, E. Boutati, E. Maratou, E. Koukkou, M. Tzanela, N. Thalassinos and S. A. Raptis, Glucose and lipid fluxes in the adipose tissue after meal ingestion in hyperthyroidism, J. Clin. Endocrinol. Metab., 91 (2006), 1112-1118. |
[32] | H. Doege, R. A. Baillie, A. M. Ortegon, B. Tsang, Q. Wu, S. Punreddy, D. Hirsch, N. Watson, R. Gimeno and A. Stahl, Targeted deletion of FATP5 reveals multiple functions in liver metabolism: alterations in hepatic lipid homeostasis, Gastroenterology, 130 (2006), 1245-1258. |
[33] | H. Doege, D. Grimm, A. Falcon, B. Tsang, T. A. Storm, H. Xu, A. M. Ortegon, M. Kazantzis, M. A. Kay and A. Stahl, Silencing of hepatic fatty acid transporter protein 5in vivo reverses diet-induced non-alcoholic fatty liver disease and improves hyperglycemia, J. Biol. Chem., 283 (2008), 22186-22192. |
[34] | B. Erdogmus, A. Tamer, R. Buyukkaya, B. Yazici, A. Buyukkaya, E. Korkut, A. Alcelik and U. Korkmaz, Portal vein hemodynamics in patients with non-alcoholic fatty liver disease, Tohoku J. Exp. Med., 215 (2008), 89-93. |
[35] | A. Falcon, H. Doege, A. Fluitt, B. Tsang, N. Watson, M. A. Kay and A. Stahl, FATP2 is a hepatic fatty acid transporter and peroxisomal very long-chain acyl-CoA synthetase, Am. J. Physiol. Endocrinol. Metab., 299 (2010), E384-E393. |
[36] | S. Fogler, Elements of Chemical Reaction Engineering, 3rd edition, Prentice Hall, New Jersey, 2001. |
[37] | E. L. Forker and B. Luxon, Hepatic transport kinetics and plasma disappearance curves: distributed modeling vs. conventional approach, Am. J. Physiol., 235 (1978), E648-E660. |
[38] | R. L. Fournier, Basic Transport Phenomena in Biomedical Engineering, Taylor & Francis, New York, 1998. |
[39] | X. Fu, J. P. Sluka, S. G. Clendenon1, K. W. Dunn, Z. Wang, J. E. Klaunig and J. A. Glazier, Modeling of Xenobiotic Transport and Metabolism in Virtual Hepatic Lobule Models, PLOS ONE, 13 (2018), e0198060. |
[40] | R. Gebhardt, Metabolic zonation of the liver: regulation and implications for liver function, Pharmacol. Ther., 53 (1992), 275-354. |
[41] | Z. Gong, E. Tas, S. Yakar and R. Muzumdar, Hepatic lipid metabolism and non-alcoholic fatty liver disease in aging, Mol. Cell. Endocrinol., 455 (2017), 115-130. |
[42] | M. R. Gray and Y. K. Tam, The series-compartment model for hepatic elimination, Drug Metab. Dispos., 15 (1987), 27-31. |
[43] | J. T. Haas, S. Francque and B. Staels, Pathophysiology and mechanisms of nonalcoholic fatty liver disease. Annu. Rev. Physiol., 78 (2016), 181-205. |
[44] | K. C. Hames, A. Vella, B. J. Kemp and M. D. Jensen, Free fatty acid uptake in humans with CD36 deficiency. Diabetes, 63 (2014), 3606-3614. |
[45] | T. Hardy, F. Oakley, Q. M. Anstee and C. P. Day, Nonalcoholic fatty liver disease: pathogenesis and disease spectrum. Annu. Rev. Pathol. Mech. Dis., 11 (2016), 451-496. |
[46] | B. S. Hijmans, A. Grefhorst, M. H. Oosterveer and A. K. Groen, Zonation of glucose and fatty acid metabolism in the liver: mechanism and metabolic consequences, Biochimie, 96 (2014), 121-129. |
[47] | W. Hundsdorfer and J. G. Verwer, Numerical Solution of Time-Dependent Advection-Diffusion-Reaction Equations, Springer, Berlin, Germany, 2003. |
[48] | K. Jungermann, Metabolic zonation of liver parenchyma: significance for the regulation of glycogen metabolism, gluconeogenesis, and glycolysis, Diabetes Metab. Rev., 3 (1987), 269-293. |
[49] | K. Jungermann and N. Katz, Functional hepatocellular heterogeneity, Hepatology, 2 (1982), 385-395. |
[50] | K. Jungermann and N. Katz, Functional specialization of different hepatocyte populations, Physiol. Rev., 69 (1989), 708-764. |
[51] | K. Jungermann and T. Kietzmann, Oxygen: modulator of metabolic zonation and disease of the liver, Hepatology, 31 (2000), 255-260. |
[52] | K. Jungermann and T. Kietzmann, Role of oxygen in the zonation of carbohydrate metabolism and gene expression in liver, Kidney Int., 51 (1997), 402-412. |
[53] | K. Jungermann and T. Kietzmann, Zonation of parenchymal and nonparenchymal metabolism in liver, Annu. Rev. Nutr., 16 (1996), 179-203. |
[54] | K. Jungermann and D. Sasse, Heterogeneity of liver parenchymal cells, Trends Biochem. Sci., 3 (1978), P198-P202. |
[55] | N. R. Katz, Metabolic heterogeneity of hepatocytes across the liver acinus, J. Nutr., 122 (1992), 843-849. |
[56] | T. Kietzmann, Metabolic zonation of the liver: the oxygen gradient revisited, Redox Biol.,11 (2017), 622-630. |
[57] | T. Kietzmann, E.Y. Dimova, D. Flúgel and J. G. Scharf, Oxygen: modulator of physiological and pathophysiological processes in the liver, Z. Gastroenterol., 44 (2006), 67-76. |
[58] | T. Kietzmann and K. Jungermann, Modulation by oxygen of zonal gene expression in liver studied in primary rat hepatocyte cultures, Cell. Biol. Toxicol., 13 (1997), 243-255. |
[59] | M. Kot, Elements of Mathematical Biology, Cambridge University Press, Cambridge, UK, 2001. |
[60] | Y. Li, C. C. Chow, A. B. Courville, A. D. Sumner and V. Periwal, Modeling glucose and free fatty acid kinetics in glucose and meal tolerance test, Theor. Biol. Med. Model., 13 (2016). |
[61] | D. Magalotti, G. Marchesini, S. Ramilli, A. Berzigotti, G. Bianchi and M. Zoli, Splanchnic haemodynamics in non-alcoholic fatty liver disease: effect of a dietary/pharmacological treatment: a pilot study, Dig. Liver Dis., 36 (2004), 406-411. |
[62] | D. G. Mashek, Hepatic fatty acid trafficking: multiple forks in the road, Adv. Nut., 4 (2013), 697-710. |
[63] | M. E. Miquilena-Colina, E. Lima-Cabello, S. Sánchez-Campos, M. V. García-Mediavilla, M. Fernández-Bermejo, T. Lozano-Rodríguez, J. Vargas-Castrillón, X. Buqe, B. Ochoa, P. Aspichueta, J. González-Gallego and C. García-Monzón, Hepatic fatty acid translocase CD36 upregulation is associated with insulin resistance, hyperinsulinaemia and increased steatosis in non-alcoholic steatohepatitis and chronic hepatitis C, Gut, 60 (2011), 1394-1402. |
[64] | H. Mitsuyoshi, K. Yasui, Y. Harano, M. Endo, K. Tsuji, M. Minami, Y. Itoh, T. Okanoue and T. Yoshikawa, Analysis of hepatic genes involved in the metabolism of fatty acids and iron in nonalcoholic fatty liver disease, Hepatol. Res., 39 (2009), 366-373. |
[65] | A. Mohammadi, M. Ghasemi-rad, H. Zahedi, G. Toldi and T. Alinia, Effect of severity of steatosis as assessed ultrasonographically on hepatic vascular indices in non-alcoholic fatty liver disease, Med. Ultrason., 13 (2011), 200-206. |
[66] | E. P. Newberry, Y. Xie, S. Kennedy, X. Han, K. K. Buhman, J. Luo, R. W. Gross and N. O. Davidson, Decreased hepatic triglyceride accumulation and altered fatty acid uptake in mice with deletion of the liver fatty acid-binding protein gene, J. Biol. Chem., 278 (2003), 51664-51672. |
[67] | H. Ohno, Y. Naito, H. Nakajima and M. Tomita, Construction of a biological tissue model based on a single-cell model: a computer simulation at metabolic heterogeneity in the liver lobule, Artif. Life, 14 (2008), 3-28. |
[68] | A. Okubo, Difffusion and Ecological Problems: Mathematical Models, Springer-Verlag, Berlin, Heidelberg, New York, 1980. |
[69] | K. S. Pang, M. Weiss and P. Macheras, Advanced pharmacokinetic models based on organ clearance, circulatory, and fractal concepts, AAPS J., 9 (2007), E268-E283. |
[70] | S. Park, S. H. J. Kim, G. E. P. Ropella, M. S. Roberts and C. A. Hunt, Tracing multiscale mechanisms of drug disposition in normal and diseased livers, J. Pharmacol. Exp. Ther., 334 (2010), 124-136. |
[71] | V. Periwal, C. C. Chow, R. N. Bergman, M. Rick, G. L. Vega and A. E. Sumner, Evaluation of quantitative models of the effect of insulin on lipolysis and glucose disposal, Am. J. Physiol. Regul. Integr. Comp. Physiol., 295 (2008), R1089-R1096. |
[72] | I. Probst, P. Schwartz and K. Jungermann, Induction in primary culture of `gluconeogenic' and `glycolytic' hepatocytes resembling periportal and perivenous cells, Eur. J. Biochem., 126 (1982), 271-278. |
[73] | G. Rajaraman, M. S. Roberts, D. Hung, G. Q. Wang and F. J. Burczynski, Membrane binding proteins are the major determinants for the hepatocellular transmembrane flux of long-chain fatty acids bound to albumin, Pharm. Res., 22 (2005), 1793-1804. |
[74] | V Rezania, D. Coombe and J. A. Tuszynski, A physiologically-based flow network model for hepatic drug elimination III: 2D/3D DLA lobule models, Theor. Biol. Med. Model., 13 (2016). |
[75] | T. Ricken, U. Dahmen and O. Dirsch, A biphasic model for sinusoidal liver perfusion remodeling after outflow obstruction. Biomech. Model. Mechanobiol., 9 (2010), 435-450. |
[76] | T. Ricken, D. Werner, H. G. Holzhütter, M. König, U. Dahmen and O. Dirsch, Modeling function-perfusion behavior in liver lobules including tissue, blood, glucose, lactate and glycogen by use of a coupled two-scale PDE-ODE approach, Biomech. Model. Mechanobiol., 14 (2015), 515-536. |
[77] | T. Ricken, N. Waschinsky and D. Werner, Simulation of steatosis zonation in liver lobule-a continuummechanical bi-scale, tri-phasic, multi-component approach, inLecture Notes in Applied and Computational Mechanics, Vol. 84, (eds. P. Wriggers and T. Lenarz), Springer, (2018), 15-33. |
[78] | M. S. Roberts and M. Rowland, Correlation between in-vitro microsomal enzyme activity and whole organ hepatic elimination kinetics: analysis with a dispersion model, J. Pharm. Pharmacol., 38 (1986), 177-181. |
[79] | M. S. Roberts and M. Rowland, Hepatic elimination-dispersion model, J. Pharm. Sci., 74 (1985), 585-587. |
[80] | S. I. Rubinow, Introduction to Mathematical Biology, John Wiley & Sons, New York, 1975. |
[81] | J. Schleicher, U. Dahmen, R. Guthke and S. Schuster, Zonation of hepatic fat accumulation: insights from mathematical modelling of nutrient gradients and fatty acid uptake, J. R. Soc. Interface, 14 (2017), 20170443. |
[82] | S. Sheikh-Bahaei, J. J. Maher and C. A. Hunt, Computational experiments reveal plausible mechanisms for changing patterns of hepatic zonation of xenobiotic clearance and hepatotoxicity, J. Theor. Biol., 265 (2010), 718-733. |
[83] | J. Shi and K. V. Kandror, Study of glucose uptake in adipose cells, Methods Mol. Biol., 456 (2008), 307-315. |
[84] | J. P. Sluka, X. Fu, Maciej Swat, J. M. Belmonte, A. Cosmanescu, S. G. Clendeno, J. F. Wambaugh, and J. A. Glazier, A liver-centric multiscale modeling framework for xenobiotics, PLOS ONE (11), 2016: DOI:10.1371/journal.pone.0162428. |
[85] | M. Soresi, L. Giannitrapani, D. Noto, A. Terranova, M. E. Campagna, A. B. Cefal, A. Giammanco and G. Montalto, Effects of steatosis on hepatic hemodynamics in patients with metabolic syndrome, Ultrasound Med. Biol., 41 (2015), 1545-1552. |
[86] | A. W. Thorburn, B. Gumbiner, F. Bulacan, P. Wallace and R. R. Henry, Intracellular glucose oxidation and glycogen synthase activity are reduced in non-insulin-dependent (type II) diabetes independent of impaired glucose uptake, J. Clin. Invest., 85 (1990), 522-529. |
[87] | C. Torre, C. Perret and S. Colnot, Molecular determinants of liver zonation, Prog. Mol. Biol. Transl. Sci., 97 (2010), 127-150. |
[88] | V. van Ginneken, E. de Vries, E. Verheij and J. van der Greef, Potential biomarkers for ``fatty liver'' (hepatic steatosis) and hepatocellular carcinoma (HCC) and an explanation of their pathogenesis, Gastroenterol. Liver Clin. Med., 1 (2017), 001. |
[89] | D. Wölfle and K. Jungermann, Long-term effects of physiological oxygen concentrations on glycolysis and gluconeogenesis in hepatocyte cultures, Eur. J. Biochem., 151, (1985), 299-303. |
[90] | Q. Wu, A. M. Ortegon, B. Tsang, H. Doege, K. R. Feingold and A. Stahl, FATP1 is an insulin-sensitive fatty acid transporter involved in diet-induced obesity, Mol. Cell. Biol., 26 (2006), 3455-3467. |
[91] | M. M. Yeh and E. M. Brunt, Pathological features of fatty liver disease, Gastroenterology, 147 (2014), 754-764. |
1. | Vahid Rezania, Dennis Coombe, Jack Tuszynski, Liver Bioreactor Design Issues of Fluid Flow and Zonation, Fibrosis, and Mechanics: A Computational Perspective, 2020, 11, 2079-4983, 13, 10.3390/jfb11010013 | |
2. | Yunjie Liao, Nathan A. Davies, I. David L. Bogle, Computational Modeling of Fructose Metabolism and Development in NAFLD, 2020, 8, 2296-4185, 10.3389/fbioe.2020.00762 | |
3. | Theodore R. Rieger, Richard J. Allen, Cynthia J. Musante, A Quantitative Systems Pharmacology Model of Liver Lipid Metabolism for Investigation of Non-Alcoholic Fatty Liver Disease, 2022, 13, 1663-9812, 10.3389/fphar.2022.910789 |
Eqns in H, HB, SB | Eqns in HB, SB | Eqns in H | |||
Glucose | | Insulin | | Glycogen | |
Triglycerides | | Glucagon | | Glucose-6-phosphate | |
Fatty acids | | Oxygen | | Glycerol-3-phosphate | |
Glycerol | | Acetyl-CoA | | ||
Lactate | | Inorganic phosphate | | ||
Guanosine Tri/Di-phosphate | | ||||
Uridine Tri/Di-phosphate | | ||||
Adenosine Tri/Di/Mono-phosphate | |
Molecule | Equation |
Glycogen | |
Glucose-6-phosphate | |
Glucose | |
Glycerol-3-phosphate | |
Lactate | |
Acetyl-CoA | |
Fatty acid | |
Triglyceride | |
Glycerol | |
Guanosine Tri-phosphate | |
Guanosine Di-phosphate | |
Uridine Tri-phosphate | |
Uridine Di-phosphate | |
Inorganic phosphate | |
Adenosine Tri-phosphate | |
Adenosine Di-phosphate | |
Adenosine Mono-phosphate |
Metabolic process | Abbrev. | Rate |
Glucokinase | | |
Glucose-6-phosphatase | | |
Glycogen Synthase | | |
Glycogen Phosphorylase | | |
Phosphoenolpyruvate Carboxykinase | | |
Fructose Biphosphatase | | |
Phosphofructokinase | | |
| ||
Pyruvate Kinase | |
Metabolic process | Abbrev. | Rate |
Pyruvate Oxidation | | |
| | |
Lipogenesis | | |
Lipolysis | | |
Triglyceride Synthesis | | |
Glycerol Kinase | |
Metabolic process | Abbrev. | Rate |
Adenosine Tri-phosphate Synthesis | | |
Uridine Diphosphate Kinase | | |
Guanosine Diphosphate Kinase | | |
Adenosine Kinase | | |
Additional Adenosine Tri-phosphate Use | | |
Control of Cellular Phosphate Levels | |
{Transport} | Rate |
Glycerol | |
Triglyceride | |
Fatty acid | |
Glucose (GLUT2) | |
Lactate | |
Molecule} | Equation |
Glycerol | |
Triglyceride | |
Fatty acid | |
Glucose | |
Lactate | |
Insulin | |
Glucagon | |
Oxygen | |
Mol. | AR Equation | ADR Equation |
| | |
| | |
| | |
| | |
| | |
| | |
| | |
| | |
Largest | Largest | Largest | |||
Variable | change | Variable | change | Variable | change |
| |||||
| |||||
| |||||
| |||||
| |||||
| |||||
| |||||
| |||||
Largest | Largest | Largest | |||
Variable | change | Variable | change | Variable | change |
| |||||
| |||||
| |||||
| |||||
| |||||
| |||||
| |||||
| |||||
Largest | Largest | Largest | |||
Variable | change | Variable | change | Variable | change |
Eqns in H, HB, SB | Eqns in HB, SB | Eqns in H | |||
Glucose | | Insulin | | Glycogen | |
Triglycerides | | Glucagon | | Glucose-6-phosphate | |
Fatty acids | | Oxygen | | Glycerol-3-phosphate | |
Glycerol | | Acetyl-CoA | | ||
Lactate | | Inorganic phosphate | | ||
Guanosine Tri/Di-phosphate | | ||||
Uridine Tri/Di-phosphate | | ||||
Adenosine Tri/Di/Mono-phosphate | |
Molecule | Equation |
Glycogen | |
Glucose-6-phosphate | |
Glucose | |
Glycerol-3-phosphate | |
Lactate | |
Acetyl-CoA | |
Fatty acid | |
Triglyceride | |
Glycerol | |
Guanosine Tri-phosphate | |
Guanosine Di-phosphate | |
Uridine Tri-phosphate | |
Uridine Di-phosphate | |
Inorganic phosphate | |
Adenosine Tri-phosphate | |
Adenosine Di-phosphate | |
Adenosine Mono-phosphate |
Metabolic process | Abbrev. | Rate |
Glucokinase | | |
Glucose-6-phosphatase | | |
Glycogen Synthase | | |
Glycogen Phosphorylase | | |
Phosphoenolpyruvate Carboxykinase | | |
Fructose Biphosphatase | | |
Phosphofructokinase | | |
| ||
Pyruvate Kinase | |
Metabolic process | Abbrev. | Rate |
Pyruvate Oxidation | | |
| | |
Lipogenesis | | |
Lipolysis | | |
Triglyceride Synthesis | | |
Glycerol Kinase | |
Metabolic process | Abbrev. | Rate |
Adenosine Tri-phosphate Synthesis | | |
Uridine Diphosphate Kinase | | |
Guanosine Diphosphate Kinase | | |
Adenosine Kinase | | |
Additional Adenosine Tri-phosphate Use | | |
Control of Cellular Phosphate Levels | |
{Transport} | Rate |
Glycerol | |
Triglyceride | |
Fatty acid | |
Glucose (GLUT2) | |
Lactate | |
Molecule} | Equation |
Glycerol | |
Triglyceride | |
Fatty acid | |
Glucose | |
Lactate | |
Insulin | |
Glucagon | |
Oxygen | |
Mol. | AR Equation | ADR Equation |
| | |
| | |
| | |
| | |
| | |
| | |
| | |
| | |
Largest | Largest | Largest | |||
Variable | change | Variable | change | Variable | change |
| |||||
| |||||
| |||||
| |||||
| |||||
| |||||
| |||||
| |||||
Largest | Largest | Largest | |||
Variable | change | Variable | change | Variable | change |
| |||||
| |||||
| |||||
| |||||
| |||||
| |||||
| |||||
| |||||
Largest | Largest | Largest | |||
Variable | change | Variable | change | Variable | change |