Original Research ARTICLE
Dynamic Metabolic Zonation of the Hepatic Glucose Metabolism Is Accomplished by Sinusoidal Plasma Gradients of Nutrients and Hormones
- 1Computational Biochemistry Group, Institute of Biochemistry, Charite—University Medicine Berlin, Berlin, Germany
- 2Institute for Computational and Imaging Science in Cardiovascular Medicine, Charite—University Medicine Berlin, Berlin, Germany
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.
The tightly controlled switch between hepatic uptake and release of glucose keeps the plasma glucose concentrations within a range between 4 and 10 mM despite largely varying carbohydrate intake and utilization. This homeostatic function of the liver with respect to plasma glucose is achieved by several enzyme-regulatory mechanisms acting on different time scales. On the short term, hormone-dependent reversible enzyme phosphorylation and changes of reaction rates elicited by concentration changes of reaction substrates/products and allosteric modulators enable a metabolic response within seconds or minutes. Recurrent activation of these fast regulatory modes is typically accompanied by slow changes in the abundance of metabolic enzymes on a time scale of hours to days (Hopgood et al., 1973; Weinberg and Utter, 1980). Both the fast and slow mode of enzyme regulation are important for the regulation of the glucose exchange flux between hepatocytes and blood plasma (Bulik et al., 2016), Owing to concentration gradients of oxygen, metabolites, hormones, and morphogens along the hepatic capillaries (sinusoids) the expression of metabolic enzymes may differ in various zones of the liver acinus. For example, the oxygen pressure decreases by 50% along the porto-central axis of the acinus (Jungermann and Kietzmann, 2000). This goes in line with the number and structure of mitochondria (Schmucker et al., 1978) and glycolytic capacities in the periportal and pericentral zone (Braeuning et al., 2006). Hepatocytes close to the portal pole (zone 1) experiencing the highest concentration of oxygen pressure are predestined for strong ATP-demanding anabolic pathways like gluconeogenesis and urea synthesis. Contrary, hepatocytes close to the venous pole of the acinus (zone 3) experience the lowest oxygen concentrations and thus possess a high glycolytic capacity, a typical feature of cells working under conditions of permanent oxygen deprivation. The heterogeneous allocation of gluconeogenetic and glycolytic capacities to different hepatocytes along the porto-central axis may even result in a situation where a certain fraction of glucose produced by periportal cells is used to fuel the glycolysis of pericentral cells (Berndt et al., 2018).
Several blood-born factors have been identified as regulators of zone-dependent gene expression of metabolic enzymes. Oxygen, glucose, the hormones glucagon and insulin, the morphogens Wnt and hedgehog and the growth factor HGF belong to the best studied factors. The various factors appear to act in a hierarchical fashion whereby the gradients of morphogens and growth factors create a basic expression pattern that is further modulated by nutrition-related factors such as oxygen, glucose, fatty acids and the hormones insulin and glucagon. In this work, we will focus on the latter group of modifiers, i.e., we restrict our model to the metabolic response of the liver to nutritional challenges and oxygen availability.
Metabolic adaptation of hepatocytes to varying oxygen pressures is mainly controlled by hypoxia-inducible transcription factors (HIFs), heterodimeric complexes consisting of a constitutively expressed β-subunit and an oxygen-sensitive α-subunit. In the liver, HIF-1α regulates primarily glycolytic genes whereas HIF-2α is known to primarily regulate genes involved in cell proliferation and iron metabolism (Ramakrishnan and Shah, 2017). In line with the falling oxygen pressure along the porto-central axis, HIFαs were found with higher levels in the less aerobic pericentral zone (Kietzmann et al., 2001). Besides oxygen, the pancreatic hormones insulin and glucagon are important drivers of zone-dependent differences in enzyme activities. The regulatory role of these hormones is 2-fold. They control the cellular cAMP level in an antagonistic manner and thus exert opposite effects on the reversible phosphorylation of key regulatory enzymes of glycolysis and gluconeogenesis as PFK2, PK, and PEPCK. The hepatic clearance of the two hormones by endocytic uptake into hepatocytes creates a concentration gradient along the porto-central axis which entails zone-dependent differences in the phosphorylation level of interconvertible enzymes. With respect to gene expression of metabolic enzymes, insulin and glucagon also control the efficiency of several transcription factors as ChREBP, SREBP-1c, CREB, and Foxo (Han et al., 2016). Both actions of glucagon and insulin are tightly interrelated and function in part through the same mechanisms. For example, the cAMP-activated protein kinase A (PKA) is responsible for phosphorylation of interconvertible enzymes such as FBPFK2 and PK, as well as for the phosphorylation of the transcription factors ChREBP and CREP (Uyeda and Repa, 2006). cAMP is produced by glucagon-induced activation of the adenylate cyclase and degraded by insulin-stimulated cAMP phosphodiesterase. Consequently, the protein level of key regulatory enzymes reflects the integral hormone levels over longer time periods.
In this work we included dynamic changes in the abundance of metabolic enzymes into our previously developed multi-scale tissue model of hepatic glucose metabolism (Berndt et al., 2018). The rates of protein synthesis and degradation were modeled by phenomenological rate equations which were parameterized by using experimentally determined protein levels at varying concentrations of oxygen, glucose, insulin, and glucagon. The central aims of our work were (i) to provide a proof of principle for integrating in a self-consistent manner the temporal gene expression of enzymes into kinetic models of cellular metabolism, (ii) to lend further support to the concept of post-differentiation patterning according to which metabolic zonation is driven by gradients of oxygen, nutrients and hormones in the capillary blood and (iii) to present a modeling approach that obviates the requirement to measure the cellular abundance of metabolic enzymes (e.g., by quantitative proteomics) in different physical states of the liver, a procedure burdened with many problems as, for example, invasive tissue sampling and protein quantification in cells separated from different zones.
The model combines a mathematical model of the sinusoidal tissue unit (STU) (Berndt et al., 2018) with a kinetic model of the protein turnover of key regulatory enzymes.
Compartment Model of Metabolite and Hormone Transport in the Sinusoidal Tissue Unit (STU)
Structurally, the STU is defined by a single sinusoid, the adjacent space of Disse and a monolayer of hepatocytes flanking the space of Disse (see Figure 1A). Functionally, the model describes the exchange of oxygen, metabolites and hormones between the sinusoidal blood, the space of Disse and the hepatocytes and the glucose metabolism within hepatocytes. The transport of metabolite and hormones within the STU is driven by diffusion and directional transport along the flow of water and blood. Lateral blood flow in the vessel is described by Hagen-Poiseuille law for fluid flow through a cylinder, water flow in the space of Disse is described by Hagen-Poiseuille law for fluid flow in a hollow cylinder. Exchange of water between the vessel and the space of Disse is driven by hydrostatic and oncotic pressure difference between the blood vessel and the space of Disse.
Figure 1. Schematic model representation (A) sinusoidal blood flow model describing blood flow, nutrient and hormone distribution within the sinusoids. The model encompasses a single sinusoid, the adjacent space of Disse and the surrounding layer of hepatocytes. It is described by morphological parameters (blood vessel radius, thickness of the space of Disse, hepatocyte thickness, hepatocyte number, sinusoid length, degree of fenestration), and systemic parameters (central and portal vein hydrostatic pressure, plasma and lymph oncotic pressure, diffusion coefficients); (B) model of carbohydrate metabolism encompassing glycolysis, glyconeogenesis and glycogen synthesis and utilization. These pathways comprise the enzymes Glucokinase (GK), Glucose-6-phosphate isomerase (GPI), Phosphofructokinase 1 (PFK1), Aldolase (ALD), Triosephosphate isomerase (TPI), Glyceraldehydephosphate dehydrogenase (GAPDH), Phosphoglycerate kinase (PGK), Phosphoglycerate mutase (PGM), Enolase (EN), Pyruvate kinase (PK) Lactate dehydrogenase (LDH), Glucose-6-phosphate phosphatase (G6P), Phosphofructokinase 2 (PFK2), Fructose-2,6-bisphosphatase (FBP2), Fructose-1,6-bisphosphatase (FBP1), Phosphoenolpyruvate carboxykinase (PEPCK), Pyruvate carboxylase (PC), Nucleoside-diphosphate kinase (NDK), Malate dehydrogenase (MDH), Pyrophosphatase (PPASE), Glucose-1-phosphate isomerase (P1PI), Glycuronosyltransferase (UGT), Glycogen phosphorylase (GP), Glycogen synthase (GS) and transporters (ER <-> cytosol: Glucose-6-phosphate transporter (Glc6PT), Glucose transporter (GlcT); mitochondrion <-> cytosol: Pyruvate transporter (PYRT), Phosphoenolpyruvate transporter (PEPT), Malate transporter (MALT); extern <-> cytosol: Glucose transporter 2 (GLUT2), Lactate transporter (LACT). Enzymes that are phosphorylated or dephosphorylated in response to insulin (Ins) and glucagon (Glu) stimulus are marked by a yellow P, allosteric modification of enzymes is marked by a red A. The model contains the metabolites: glucose (Glc), glucose-6-phosphate (Glc6P), fructose-6-phosphate (Fru6P), fructose-1,6-bisphosphate (Fru16P2), glyceraldehydephosphate (GraP), dihydroxyacetonephosphate (DHAP), 1,3-bisphosphoglycerate (13P2G), 3-phosphoglycerate (3PG), 2-phosphoglycerate (2PG), phosphoenolpyruvate (PEP), pyruvate (Pyr), lactate (Lac), malate (Mal), oxaloacetate (OA), glucose-1-phosphate (Glc1P), UDP-glucose (UDP-glc), glycogen, fructose-2,6-bisphosphate (Fru26P2). The cofactors NADH, NAD, ATP, ADP, UTP and UDP are not treated as dynamic variables. All physiological metabolites produced or consumed in the hepatocyte during glycolysis and gluconeogenesis are comprised into lactate. Figures adapted from Berndt et al. (2018).
Conceptually, the STU was divided into NH zones, where NH is the number of hepatocytes along the porto-central axis. Each zone is made up of the sinusoid volume, the space of Disse and the hepatocyte (Figure 1A). Within one zone, the concentration of metabolites and hormones is given by a single value. The mathematical description of the STU model and a complete list of parameters used can be found in Berndt et al. (2018).
Kinetic Model of Hepatocyte Glucose Metabolism
The reaction scheme for the glucose metabolism of a single hepatocyte is depicted in Figure 1B. It consists of the pathways for glycolysis, glyconeogenesis, glycogen synthesis and degradation. The time-dependent variation of metabolite concentrations is given by first-order differential equations. The liver specific enzymatic rate laws take into account substrate regulation, allosteric regulation and hormonal regulation by hormone-dependent reversible phosphorylation (Bulik et al., 2016).
Kinetic Model of Hormonal Signaling
The pancreatic hormones glucagon and insulin are released into the portal vein in response to the plasma glucose concentration and are partially cleared during their passage through the liver. Hence, there is a difference between their plasma concentrations determined in peripheral blood samples and effective intra-hepatic concentrations. This difference was taken into account by setting the concentration values of insulin and glucagon in the periportal blood to the 2-fold of their plasma values (Balks and Jungermann, 1984). The rate of intra-hepatic hormone clearance via receptor binding and subsequent endocytosis was put proportional to the binding and signaling strength of the hormone. We used empirical transfer functions to describe the relationship between glucose and hormone concentrations in the plasma and the relationship between and the phosphorylation state of interconvertible enzymes as described in Bulik et al. (2016). A detailed description of the functions and their construction can be found in Berndt et al. (2018).
Kinetic Model of Protein Turnover
The temporal change of the protein level PENZ of a metabolic enzyme (ENZ) is given by the difference between the rates of protein synthesis and protein degradation :
The right-hand side of equation (1), vENZ(E), represents the turnover rate of the enzyme protein. It is controlled by modulators affecting either the synthesis or the degradation or both. Note that the enzyme level PENZ scales linearly with the maximal rate of the enzyme. The rate equations of protein synthesis and degradation both depend on the momentary concentration of at least one of the four modulators Ei, i = 1 (insulin), i = 2 (glucagon), i = 3 (glucose), i = 4 (oxygen) considered in the model. The general structure of the rate equation for the protein synthesis of enzyme ENZ reads.
where is a constant determining the basal synthesis rate and fi is a nonlinear function of the i-th modulator. The “+” sign holds if Ei is an activator (inductor) of protein synthesis, the (–) sign holds If Ei is an inhibitor (repressor). If Ei has not been reported so far to exert an effect on the protein synthesis of enzyme ENZ it holds and Numerical values for the rate constants and were fixed in such a manner that for a normal 24 h plasma profile (see below) the zone- and time averaged protein levels coincided with the stationary protein levels as reported in Bulik et al. (2016). Numerical values of all other kinetic parameters were obtained by adjusting the rate equations to experimentally determined protein levels at varying concentrations of the four possible modulators (see Supplement 1). Table 1 depicts the rate equations for the synthesis and degradation of those enzymes of hepatic glucose metabolism possessing in the model timely variable protein levels.
In order to quantify the sensitivity of the turnover rate vENZ(E) of a protein against small changes of an modulator E, we used the sensitivity (elasticity) coefficient as defined in metabolic control analysis:
Figure 2 depicts the sensitivity coefficients for the turnover rates of the nine enzyme proteins with variable expression level as function of the four modulators oxygen (I), glucagon (II), insulin (III) and glucose (IV). Except for the sensitivities of the PEPCK and G6PP turnover with respect to glucagon and glucose, respectively, the extremum of all other sensitivity characteristics lies within the reported physiological range of the related modulators (green-shaded areas in Figure 2). The sensitivity of G6PP turnover with respect to glucose becomes important in the diabetic case, where glucose levels can exceed 20 mM (see below).
Figure 2. Sensitivity coefficients of protein turnover rates defined in equation (3) as function of modulator concentrations.I Sensitivity of protein turnover rates with respect to oxygen of (A) GK, (B) FBPFK2, (C) PK, (D) PEPCKII Sensitivity of protein turnover rates with respect to (E) G6PP, (F) FBP1, (G) FBPFK2, (H) PK, (I) PC, (J) PEPCKIII Sensitivity of protein turnover rates with respect to (K) GK, (L) PFK1, (M) PK, (N) PEPCKIV Sensitivity of protein turnover rates with respect to (O) GlcT, (P) G6PPThe green-shades areas indicate the reported physiological concentration range of the respective modulator.
Dynamic Metabolic Zonation in a Well-Fed State of the Rat
First, we used the model to simulate the temporal variation of enzyme abundances, metabolite concentrations and fluxes within the various zones along the porto-central axis of the STU. The simulation was initiated with identical abundance of enzymes along the sinusoid which we set to the stationary mean protein abundance used in Bulik et al. (2016). We used as model input the diurnal glucose profile reported for the healthy normal liver of a fed rat (La Fleur et al., 1999) and carried out the numerical integration of the model over several 24 h cycles until there was no change in the 24 h enzyme and metabolite profiles.
Even with identical enzyme abundances across all hepatocytes, there occurs a progressive decline of hormone plasma levels from the portal to the central pole due to the ongoing hormone uptake by hepatocytes in each zone. Moreover, oxygen uptake in one zone diminishes the available oxygen pressure seen by the cells in the adjacent zone toward the pericentral pole. As oxygen is not part of the model, we assumed a linear decrease in oxygen partial pressure from 90 mmHG in the periportal zone to 35 mm HG in the pericentral zone (Jungermann and Kietzmann, 2000; Allen et al., 2005). These initial gradients of hormones, glucose and oxygen feed back to the level of metabolic enzymes so that finally zone-dependent patterns of both enzyme abundances and metabolic variables (see Figures 3, 4) are generated.
Figure 3. Diurnal variations in the plasma levels of glucose (A), insulin (B), glucagon (C), cellular glycogen (D) and the glucose exchange flux (E) in different zones along the porto-central axis. The different curves refer to different spatial positions of hepatocytes, counted from periportal (red curve) to percentral (green curve). The bold blue line refers to the means values of the shown variable. Note that the red curves (= most portal cell) for the hormones and glucose are identical with their plasma profiles.
Figure 4. Diurnal variations in the relative abundance of glycolytic and gluconeogenetic enzymes within hepatocytes along the porto-central axis of a well-fed rat. The curves illustrate the relative deviation of protein abundances from the overall spatial and 24h- mean (= 24 h mean of the bold red curve). From the top to the bottom, the different curves refer to the spatial position of the hepatocyte counted from periportal to percentral. The bold red line refers to the average protein abundance across all cells. (A) GlcT, (B) G6PP, (C) GK, (D) PFK1, (E) FBPPFK1, (F) PFK2, (G) FBPPFK2, (H) PK, (I) PC, (J) PEPCK.
Figures 3A–C depicts the timely variation of glucose, insulin, glucagon in various sinusoidal compartments. Intriguingly, the highest glucose concentrations in the very portal zone (see red curve in Figure 3A) are paralleled by the lowest glucose concentrations in the very central zone (green curve). This seemingly paradoxical situation is due to the fact that the high level of insulin and low level of glucagon strongly increase the glucose uptake capacity of hepatocytes such that the otherwise strong zone-dependent differences in the glucose exchange flux (see Figure 3E) almost disappear. The simulation also reveals large zone-dependent differences in the cellular dynamics of glycogen (Figure 3D). The variation of the glycogen content in portal cells is much more pronounced than in central cells.
The time-dependent variation in the protein levels of key glycolytic and gluconeogenetic enzymes in different zones are depicted in Figure 4. The uniform overall shape of the curves reflects essentially the daily variation of the plasma glucose level. Generally, the daily fluctuations of enzyme levels around their 24 h mean hardly exceed 10%. Thus, as long as the liver is repeatedly confronted with the same 24 h plasma profile of metabolites and hormones, timely variations of protein levels should have only a marginal impact on the hepatic control of the plasma glucose level.
In contrast to the modest time-dependent variations of protein levels, the computed zone-dependent differences of enzyme levels display a large scatter. The maximal differences between the enzyme endowment of hepatocytes closest to the portal and central pole lie between 0.1 [e.g., glucose transporter (glcT) and phosphofructokinase 1 (pfk1)] and 4.5 [phosphoenolpyruvate kinase (pepck)]. For the validation of these computational predictions, we calculated the 24 h-average protein levels of the first (most portal) and last (most central) hepatocyte and compared the ratio of the computed average protein levels with experimental data (see first columns for each enzyme in Figure 5). We further compared the ratio of 24 h- and zone-averaged mean protein levels between a fed and fasted rat and a diabetic and normal rat.
Figure 5. Ratio of enzyme levels in hepatocytes. Red circles indicate the computed ratio of protein levels. Green circles and gray bars indicate the mean value and the range of variability of various experimentally determined ratios. First columns (pc): Portal-to-central ratios (= 24h-averaged protein levels between hepatocyte #1 and hepatocyte # 25) for the well-fed state of the rat. Experimental data were taken from (Sharma et al., 1964; Katz et al., 1977a,b; Trus et al., 1980; Jungermann and Katz, 1982, 1989; Probst et al., 1982; Teutsch and Lowry, 1982; Miethke et al., 1985; Quistorff, 1985; Wölfle and Jungermann, 1985; Morselt et al., 1987; Chen and Katz, 1988; Wals et al., 1988; Evans et al., 1989; Frederiks et al., 1991; Jones and Titheradge, 1996; Minchenko et al., 2003). Second columns (ff): Fed-to-fasted ratios of 24h- and zone-averaged mean enzyme levels. Experimental data taken from (Dipietro and Weinhouse, 1960; Wimhurst and Manchester, 1970b; Ballard and Hopgood, 1973; Bock et al., 1973; Zalitis and Pitot, 1979; Cladaras and Cottam, 1980; Neely et al., 1981; Bahnak and Gold, 1982; Van Schaftingen and Hers, 1983; Donofrio et al., 1984; Colosia et al., 1988; Crepin et al., 1988; Thorens et al., 1990; Giffin et al., 1993; Gannon and Nuttall, 1997). Third columns (nd): Normal-to-diabetic ratios of 24h- and zone-averaged mean enzyme levels. Experimental data taken from Dipietro and Weinhouse, 1960; Salas et al., 1963; Exton and Park, 1965; Wimhurst and Manchester, 1970a; Chang and Schneider, 1971; Dunaway et al., 1978; Weinberg and Utter, 1980; Neely et al., 1981; Donofrio et al., 1984; Miethke et al., 1985; Colosia et al., 1988; Crepin et al., 1988; Thorens et al., 1990; Miralpeix et al., 1992; Slieker et al., 1992; Gannon and Nuttall, 1997; Raju et al., 1999; Manna and Jain, 2012.
This analysis provided a good concordance between theoretical and experimental results. The only exception is the pyruvate carboxylase, a key regulatory enzyme of gluconeogenesis, were portal to central gradients could not be univocally explained by the reported oxygen dependency. Oxygen dependency accounts only for about 35% percent of the observed zonation (see Table 1 and Supplement 1).
Dynamic Metabolic Zonation of the Liver During Adaptation to Fasting
Next, we studied how the zonation of metabolic enzymes is affected if the liver has to cope with a fundamentally different nutritional regime. To this end, we simulated the zone-dependent dynamic changes of protein levels and metabolites during the transition from a fed state of the rat to a fasting state. The simulation started with the stable 24 h zonated enzyme profile that is established if the liver experiences recurrently the same plasma profile of a fed rat (see above). At time t = 24 h, the plasma profile of the fed rate was replaced by plasma profile of a fasted rat (La Fleur et al., 1999). The latter was only available for a time range of 24 h of fasting. After about 16 h of fasting, stable values of plasma metabolites were reached. Therefore, we used for time points t > 48 h (i.e., >24 h of fasting) for the model input a plasma profile that was composed of 6 repetitions of the last part (16–24 h) of the 24 h fasting plasma profile. As shown in Figure 6, the fed-to-fasting transition evokes a significant rise in the abundance of key gluconeogenetic enzymes (GlcT, FBP1, PC, PEPCK) and drop in the abundance of key glycolytic enzymes (GK, PFK1, PFK2, FBP2, PK) in all zones. For two enzymes, the GlcT and the PEPCK, the zone-dependent protein differences become more pronounced compared to the fed state. By contrast, for the GK, FBP1, GSPP, and PK the zone-dependent protein differences became smaller.
Figure 6. Diurnal variations in the relative abundance of glycolytic and gluconeogenetic enzymes within hepatocytes along the porto-central axis during the transition from a fed state (t = 0–24 h) to a fasted state (t > 24 h). The curves illustrate the relative deviation of protein abundances from the overall spatial and 24 h mean (= 24 h mean of the bold red curve). From the top to the bottom, the different curves refer to the spatial position of the hepatocyte counted from periportal to percentral. The bold red line refers to the average protein abundance across all cells. The vertical dotted line indicates onset of the starvation period. (A) GlcT, (B) G6PP, (C) GK, (D) PFK1, (E) FBPPFK1, (F) PFK2, (G) FBPPFK2, (H) PK, (I) PC, (J) PEPCK.
The computed changes of enzyme profiles toward a more gluconeogenetic phenotype are accompanied by significant alterations of the intra-sinusoidal glucose gradient and the zone-dependent differences in the glucose exchange rate (Figure 7). Compared with the fed state, the porto-venous glucose difference becomes much larger in the fasted state (Figure 7A). The same holds for the glucose exchange rate (Figure 7D).
Figure 7. Diurnal variations in the plasma levels of glucose (A), insulin (B), glucagon (C), cellular glycogen (D) and the glucose exchange flux (E) in different zones along the porto-central axis during the transition from a well-fed state (t = 0–24 h) to a fasted state (t > 24 h) of the rate. The different curves refer to different spatial positions of hepatocytes, counted from periportal (red curve) to percentral (green curve). The bold blue line refers to the means values of the shown variables. Note that the red curves (= most portal cell) for the hormones and glucose are identical with their plasma profiles.
In agreement with experimental data (Babcock and Cardell, 1974), the glycogen stores are almost depleted after about 1 day of fasting when the levels of insulin and glucagon have adopted a new stable temporal profile. Notably, also for the fasted state, the computed average protein abundance ratios are in good agreement with experimental data which further substantiates the reliability of the model (see Figure 5).
Figure 8 illustrates the importance of dynamic zonation for the adaptation of the porto-venous glucose difference (AVGD) to a specific nutritional regime. Regulation of interconvertible enzymes by hormone-dependent phosphorylation alone, i.e., at fixed protein levels of the fed state (blue line), would result in an AVGD of about 3.5 mM for the typical range of portal glucose concentrations in the fasted state (red shaded area). Dynamic adaption of protein levels enlarges the AVGD to about 7 mM (red line) thus rendering the liver to a strong glucose producer in the fasted state.
Figure 8. Hepatic porto-venous glucose difference (AVGD) for the well-fed and fasted nutritional state. The daily variations of plasma glucose in the fasted, fed and well-fed state are indicated by the red-, green- and blue-shaded areas. The solid lines represent the AVGD if in the fasted state (red), fed state (green) and the well-fed state (blue). Crosses depict experimentally measured AVGD (Huang and Veech, 1988).
Dynamic Metabolic Zonation of the Liver in Diabetes Type II (“Diabetic Liver”)
Late diabetes type 2 is characterized by long persistence of high postprandial plasma glucose levels, reduced insulin levels (hypo-insulinemia) and elevated glucagon levels (hyper-glucagonemia). It is mainly the shift in the insulin/glucagon ratio that renders the liver to a glucose producer which on top of the insulin-resistant muscle and adipose tissue contributes to high plasma glucose levels. We tested whether our model can also correctly describe this metabolic abnormality and the observed changes of protein abundances in different zones. To this end, we used the glucose-hormone transfer function constructed for the diabetic case (see Bulik et al., 2016) to calculate the phosphorylation state of interconvertible enzymes and confronted the model over 4 days repeatedly with a 24 h glucose profile of a diabetic rat until an almost stable 24 h pattern of protein abundances had established (Figure 9). As shown in Figure 9D, all hepatocytes work permanently as glucose producers, i.e., at all time points of the day the intra-sinusoidal glucose concentration increases along the portal-central axis (Figure 9A). The glycogen reserves are drastically diminished and the different time courses of glycogen emptying and filling between portal and central regions (see Figure 2D) are completely abolished. Taken together, the zonation of glucose metabolism in the diabetic liver bears a strong resemblance with that of the fasted liver. It is important to note that even in this pathophysiological case the computed average portal-to-central protein abundance ratios are in good agreement with experimental data (see Figure 5).
Figure 9. Diurnal variations in the plasma levels of glucose (A), insulin (B), glucagon (C), cellular glycogen (D), and the glucose exchange flux (E) in different zones along the porto-central axis of a diabetic rat.
Metabolic Zonation of Hepatic Glucose Metabolism Is Driven by Concentration Gradients of Hormones and Metabolites
In this work we used a mathematical model to study the dynamic zonation of the hepatic glucose metabolism. To this end we extended our previously published multi-scale tissue model of the hepatic carbohydrate metabolism (Berndt et al., 2018) by rendering the protein levels of key regulatory enzymes of glycolysis and gluconeogenesis as dynamic model variables which are controlled by timely variable synthesis and degradation in dependence from the concentration of the four modulators glucose, insulin, glucagon and oxygen. The model correctly replicates experimentally determined protein levels in different zones of the liver acinus as well as the adaptation of the liver to a well-fed, fasted and diabetic state. From this we draw four main conclusions. (1) Zonation of the hepatic glucose metabolism is a necessary consequence of the fact that the expression of key regulatory enzymes is controlled by modulators that display a porto-central concentration gradient along the sinusoid. (2) Mechanisms controlling the adaptation of enzyme abundances to varying external conditions necessarily lead to the zonation of hepatic carbohydrate metabolism. (3) The four modulators considered in the model are sufficient to describe the dynamic zonation of the glucose metabolism of the a normal liver. (4) The use of phenomenological transfer functions which directly relate the protein turnover to known modulators of gene expression appears a promising modeling strategy to include variable protein levels in kinetic models in view of the fact that in a foreseeable future explicit kinetic modeling of complex gene-regulatory network is out of reach.
The Proposed Multi-Scale Model Encompasses All Levels of Metabolic Regulation
An important feature of the cellular metabolic network of the liver is the ability to adapt its functional output to varying external conditions such as changes in nutrient supply and varying hormone levels. These adaptive mechanisms operate at two different time scales. The short term adaptation occurs within seconds or minutes and is brought about by activity changes in the present metabolic enzymes by substrate availability, allosteric regulation and reversible phosphorylation. The second adaptive mechanism operates within hours or days and is brought about by changes in the enzyme abundances. It is already known for a long time that the total protein content of liver enzymes may largely vary owing to enhanced protein degradation during fasting (providing glucogenic amino acids as substrate for gluconeogenesis) mediated by the hormone glucagon or enhanced protein synthesis by the hormone insulin (Hopgood et al., 1980). However, such general changes of the protein content do not tell anything about the changes in the abundance of individual enzymes, whose expression by insulin and glucagon differs profoundly. Therefore, it was necessary to develop empirical rate laws for the synthesis and degradation of individual enzymes. This phenomenological approach was chosen since currently biochemical information is insufficient to establish molecular-resolved kinetic models which include, for example, the interaction of transcription factors among each other and with specific DNA promotor regions, the processing of mRNA and the regulation of mRNA translation by micro RNAs and RNA-binding proteins. For example, PEPCK, probably the best-studied gluconeogenetic enzyme, is regulated by at least a dozen transcription factors with partially unknown interactions (Yang et al., 2009). Even if it was possible to explicitly model the mRNA transcription for individual enzymes, there is still a big gap in understanding post-transcriptional regulation and the processes of post-translational modification.
Metabolic Response of the Liver to Varying Nutritional Regimes
Our simulations suggest that in the presence of a constant daily nutritional regime the diurnal variation of enzyme abundances should be fairly moderate in the range of 10–20% around the mean. This is a lot less than daily variations in the abundance of the key regulatory enzyme of cholesterol synthesis, ßHMG-CoA reductase (Kirkpatrick et al., 1980), exhibiting a pronounced circadian rhythm or some enzymes of the amino acid metabolism as, for example, tyrosine transaminase the activity of which is almost four times as great several hours after nightfall as it is in the morning (Wurtman, 1974). In contrast, much larger changes of glycolytic and gluconeogenetic enzyme levels are elicited by a switch from well-fed to fasting conditions and vice versa. This metabolic adaptation occurs within a time span of several days (see Figure 6) and enables vertebrates to maintain the plasma glucose level in the absence of food. The slow change of enzyme concentration profiles implies that the capability of the fasted liver to clear a sudden excess of plasma glucose is diminished (impaired glucose tolerance) as the capacity of glycolytic enzymes and enzymes of the glycogen pathway a downregulated (Bulik et al., 2016). Intriguingly, there appears to be a striking similarity in the adaptive response of the liver to fasting conditions and diabetes type 2. In our modeling approach, this is mainly due to the fact that in both physiological settings the strong effect of insulin on the expression of glycolytic and gluconeogenetic enzymes is diminished whereas the effect of glucagon is more pronounced. As demonstrated in a previous model-based simulation study (König and Holzhütter, 2012), exposing the permanently glucose-releasing “starved” liver of the diabetic patient to a rigorous insulin treatment at persistently elevated glucagon level may increase the risk of severe hypoglycaemia. Thus, owing to the equally strong impact of both insulin and glucagon on the expression and phosphorylation state of liver enzymes, reconverting the “starved” liver of the diabetic patient into the normal metabolic phenotype requires the normalization of the plasma levels of both hormones.
Main Limitations of the Model and Outlook for Future Model Extensions
The used multiscale tissue model comprises a number of simplifications of the true anatomical structure of the liver which may impact on the simulated intra-sinusoidal concentration gradients of hormones and metabolites. For example, the blood flow rate within the pericentral zone of the sinusoid may vary if a sinusoid spreads out, forms anastomoses or merges with another sinusoid (Rappaport et al., 1954), anatomic peculiarities of liver parenchyma that are not yet considered in the STU model. Also, the number of periportal hepatocytes is higher (~2–3-fold) compared to pericentral hepatocytes. Regarding the concentration gradient of oxygen, hormones and metabolites in the sinusoid it may be of relevance that the terminal branches of the hepatic artery rarely join with the portal vein already before the blood enters the sinusoids, as presumed in our model. In the vast majority the merger of arterial blood with blood from the portal vein, occurs a few cells downstream within the sinusoid (Ekataksin and Kaneda, 1999) resulting in a local increase of the concentration of oxygen concentration at this site. Despite these limitations, the model correctly describes glucose exchange rates, gradients and indicator dilution curves for a structurally normal liver (Berndt et al., 2018). Hence, the functional implications of the above limitations and other neglected aspects of the real topology of liver tissue remain unclear. Therefore, in future work we aim to embed our metabolic cell model in a 3-D reconstructions of a mouse lobule (Hoehme et al., 2017).
The rate laws presented in this paper are effective transfer functions describing directly the relation between modulators (nutrients and hormones) and the turnover rate of a protein. Usage of effective transfer function raises the question which properties of the underlying regulatory network have to be captured. Obviously, not all known modulators of protein synthesis and degradation have been considered in the model. Ground-breaking experiments pointed initially to the oxygen gradient as the most important driving force of metabolic zonation (Jungermann and Kietzmann, 1997, 2000). Later experiments with isolated hepatocytes incubated with varying concentrations of insulin or glucagon (Probst et al., 1982) revealed an important role of these hormones for the establishment of liver zonation. Meanwhile a lot more potential modulators of metabolic zonation have been described in the literature, among them Wnt/β-catenin pathway (Torre et al., 2010; Vasilj et al., 2012), MAPK/ERK pathway (Zeller et al., 2013), Hnf4-alpha (Colletti et al., 2009), or thyroid hormones (Weinberg and Utter, 1979). However, as demonstrated in this study, the dynamic zonation of the glucose metabolism can be well described in different physiological settings by taken into account only the four modulators oxygen, glucose, glucagon and insulin. The central role of oxygen, glucose, glucagon and insulin for the dynamic zonation of the glucose metabolism does not exclude that morphogens and growth factors may have an important role in the zonation of other metabolic subsystems (Gebhardt and Matz-Soja, 2014). For example, the expression of the glutamine synthetase is restricted to last few hepatocytes close to the venous pole. Complementary, the urea cycle enzyme carbamoylphosphate synthetase I (CPS I) is present in the periportal, intermediate, and the first few layers of the perivenous zone. It has been clearly demonstrated that Wnt/ß-catenin signaling pathway plays a central role in the stable maintenance of these peculiar zonation profiles (Burke et al., 2009).
In summary, we propose a self-consistent model of liver carbohydrate metabolism that consistently takes into account variable gene expression of metabolic enzymes, regulation of metabolic pathways, exchange of metabolites and hormones between the blood and hepatocytes and microperfusion of the liver. Once the input of hormones and nutrients to the periportal region the liver acinus is known, the model allows to compute the metabolic phenotype of individual hepatocytes along the porto-central axis. The local hormone and metabolite concentrations determine the phosphorylation state of the interconvertible enzymes, hormonal clearance rates and expression level of metabolic enzymes. The metabolic phenotype in turn determines the functional output (here: glucose exchange rate) of each hepatocyte and this way the venous glucose output of the acinus. Integration across a representative set of acini yields finally the total glucose output of the liver.
NB developed the concept, implemented the model, carried out the simulation, wrote the manuscript. H-GH developed the concept, advised the implementation of the model, wrote the manuscript.
NB was funded by the German Systems Biology Program “LiSyM,” grant no. 31 L0057, sponsored by the German Federal Ministry of Education and Research (BMBF). The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.
Conflict of Interest Statement
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fphys.2018.01786/full#supplementary-material
Argaud, D., Zhang, Q., Pan, W., Maitra, S., Pilkis, S. J., and Lange, A. J. (1996). Regulation of rat liver glucose-6-phosphatase gene expression in different nutritional and hormonal states: gene structure and 5'-flanking sequence. Diabetes 45, 1563–1571. doi: 10.2337/diab.45.11.1563
Berndt, N., Horger, M. S., Bulik, S., and Holzhütter, H. G. (2018). A multiscale modelling approach to assess the impact of metabolic zonation and microperfusion on the hepatic carbohydrate metabolism. PLoS Comput. Biol. 14:e1006005. doi: 10.1371/journal.pcbi.1006005
Braeuning, A., Ittrich, C., Köhle, C., Hailfinger, S., Bonin, M., Buchmann, A., et al. (2006). Differential gene expression in periportal and perivenous mouse hepatocytes. FEBS J. 273, 5051–5061. doi: 10.1111/j.1742-4658.2006.05503.x
Bulik, S., Holzhütter, H. G., and Berndt, N. (2016). The relative importance of kinetic mechanisms and variable enzyme abundances for the regulation of hepatic glucose metabolism-insights from mathematical modeling. BMC Biol. 14:15. doi: 10.1186/s12915-016-0237-6
Burke, Z. D., Reed, K. R., Phesse, T. J., Sansom, O. J., Clarke, A. R., and Tosh, D. (2009). Liver zonation occurs through a beta-catenin-dependent, c-Myc-independent mechanism. Gastroenterology 136, 2316–2324 e2311–2313. doi: 10.1053/j.gastro.2009.02.063
Christ, B., Nath, A., Bastian, H., and Jungermann, K. (1988). Regulation of the expression of the phosphoenolpyruvate carboxykinase gene in cultured rat hepatocytes by glucagon and insulin. Eur. J. Biochem. 178, 373–379. doi: 10.1111/j.1432-1033.1988.tb14460.x
Colletti, M., Cicchini, C., Conigliaro, A., Santangelo, L., Alonzi, T., Pasquini, E., et al. (2009). Convergence of Wnt signaling on the HNF4 alpha-driven transcription in controlling liver zonation. Gastroenterology 137, 660–672. doi: 10.1053/j.gastro.2009.05.038
Colosia, A. D., Marker, A. J., Lange, A. J., el-Maghrabi, M. R., Granner, D. K., Tauler, A., et al. (1988). Induction of rat liver 6-phosphofructo-2-kinase/fructose-2,6-bisphosphatase mRNA by refeeding and insulin. J. Biol. Chem. 263, 18669–18677.
Crepin, K. M., Darville, M. I., Hue, L., and Rousseau, G. G. (1988). Starvation or diabetes decreases the content but not the mRNA of 6-phosphofructo-2-kinase in rat liver. Febs Lett. 227, 136–140. doi: 10.1016/0014-5793(88)80884-6
Dice, J. F., and Goldberg, A. L. (1975). Statistical-analysis of relationship between degradative rates and molecular-weights of proteins. Arch. Biochem. Biophys. 170, 213–219. doi: 10.1016/0003-9861(75)90112-5
Donofrio, J. C., Thompson, R. S., Reinhart, G. D., and Veneziale, C. M. (1984). Quantification of liver and kidney phosphofructokinase by radioimmunoassay in fed, starved and alloxan-diabetic rats. Biochem. J. 224, 541–547. doi: 10.1042/bj2240541
Dunaway, G. A., Leung, G. L., Thrasher, J. R., and Cooper, M. D. (1978). Turnover of hepatic phosphofructokinase in normal and diabetic rats-role of insulin and peptide stabilizing factor. J. Biol. Chem. 253, 7460–7463.
Dunaway, G. A., and Weber, G. (1974). Effects of hormonal and nutritional changes on rates of synthesis and degradation of hepatic phosphofructokinase isozymes. Arch. Biochem. Biophys. 162, 629–637. doi: 10.1016/0003-9861(74)90225-2
Frederiks, W. M., Marx, F., and van Noorden, C. J. (1991). Homogeneous distribution of phosphofructokinase in the rat liver acinus: a quantitative histochemical study. Hepatology 14, 634–639. doi: 10.1002/hep.1840140410
Gabbay, R. A., Sutherland, C., Gnudi, L., Kahn, B. B., O'Brien, R. M., Granner, D. K., et al. (1996). Insulin regulation of phosphoenolpyruvate carboxykinase gene expression does not require activation of the Ras/mitogen-activated protein kinase signaling pathway. J. Biol. Chem. 271, 1890–1897. doi: 10.1074/jbc.271.4.1890
Giffin, B. F., Drake, R. L., Morris, R. E., and Cardell, R. R. (1993). Hepatic lobular patterns of phosphoenolpyruvate carboxykinase, glycogen synthase, and glycogen phosphorylase in fasted and fed rats. J. Histochem. Cytochem. 41, 1849–1862. doi: 10.1177/41.12.8245433
Hoehme, S., Friebel, A., Hammad, S., Drasdo, D., and Hengstler, J. G. (2017). Creation of three-dimensional liver tissue models from experimental images for systems medicine. Methods Mol. Biol. 1506, 319–362. doi: 10.1007/978-1-4939-6506-9_22
Hopgood, M. F., Ballard, F. J., Reshef, L., and Hanson, R. W. (1973). Synthesis and degradation of phosphoenolpyruvate carboxylase in rat-liver and adipose-tissue-changes during a starvation-refeeding cycle. Biochem. J. 134, 445–453. doi: 10.1042/bj1340445
Hopgood, M. F., Clark, M. G., and Ballard, F. J. (1980). Protein-degradation in hepatocyte monolayers-effects of glucagon, adenosine 3'-5'-cyclic monophosphate and insulin. Biochem. J. 186, 71–79. doi: 10.1042/bj1860071
Hopkirk, T. J., and Bloxham, D. P. (1979). Studies on the biosynthesis of hepatic pyruvate-kinase and its correlation with enhanced hepatic lipogenesis in meal-trained rats. Biochem. J. 182, 383–397. doi: 10.1042/bj1820383
Hopkirk, T. J., and Bloxham, D. P. (1980). Biosynthesis of rat-liver pyruvate-kinase-measurement of enzyme lifetime and the rate of synthesis at weaning. Biochem. J. 192, 507–516. doi: 10.1042/bj1920507
Iynedjian, P. B., Jotterand, D., Nouspikel, T., Asfari, M., and Pilot, P. R. (1989). Transcriptional induction of glucokinase gene by insulin in cultured liver cells and its repression by the glucagon-cAMP system. J. Biol. Chem. 264, 21824–21829.
Jones, C. G., and Titheradge, M. A. (1996). Measurement of metabolic fluxes through pyruvate kinase, phosphoenolpyruvate carboxykinase, pyruvate dehydrogenase, and pyruvate carboxylase in hepatocytes of different acinar origin. Arch. Biochem. Biophys. 326, 202–206. doi: 10.1006/abbi.1996.0066
Katz, N., Teutsch, H. F., Jungermann, K., and Sasse, D. (1977a). Heterogeneous reciprocal localization of fructose-1,6-bisphosphatase and of glucokinase in microdissected periportal and perivenous rat-liver tissue. FEBS Lett. 83, 272–276. doi: 10.1016/0014-5793(77)81021-1
Katz, N., Teutsch, H. F., Sasse, D., and Jungermann, K. (1977b). Heterogeneous distribution of glucose-6-phosphatase in microdissected periportal and perivenous rat-liver tissue. FEBS Lett. 76, 226–230. doi: 10.1016/0014-5793(77)80157-9
Kietzmann, T., Cornesse, Y., Brechtel, K., Modaressi, S., and Jungermann, K. (2001). Perivenous expression of the mRNA of the three hypoxia-inducible factor alpha-subunits, HIF1alpha, HIF2alpha and HIF3alpha, in rat liver. Biochem. J. 354, 531–537. doi: 10.1042/bj3540531
Kietzmann, T., Roth, U., Freimann, S., and Jungermann, K. (1997). Arterial oxygen partial pressures reduce the insulin-dependent induction of the perivenously located glucokinase in rat hepatocyte cultures: mimicry of arterial oxygen pressures by H2O2. Biochem. J. 321(Pt 1),17–20.
Kirkpatrick, R. B., Robinson, S. F., and Killenberg, P. G. (1980). Diurnal-variation of rat-liver enzymes catalyzing bile-acid conjugation and sulfation. Biochim. Biophys. Acta. 620, 627–630. doi: 10.1016/0005-2760(80)90154-X
König, M., and Holzhütter, H. G. (2012). Kinetic modeling of human hepatic glucose metabolism in type 2 diabetes mellitus predicts higher risk of hypoglycemic events in rigorous insulin therapy. J. Biol. Chem. 287, 36978–36989. doi: 10.1074/jbc.M112.382069
La Fleur, S. E., Kalsbeek, A., Wortel, J., and Buijs, R. M. (1999). A suprachiasmatic nucleus generated rhythm in basal glucose concentrations. J. Neuroendocrinol. 11, 643–652. doi: 10.1046/j.1365-2826.1999.00373.x
Leskes, A., Siekevitz, P., and Palade, G. E. (1971). Differentiation of endoplasmic reticulum in hepatocytes : II. glucose-6-phosphatase in rough microsomes. J. Cell Biol. 49, 288–302. doi: 10.1083/jcb.49.2.288
Manna, P., and Jain, S. K. (2012). Decreased hepatic phosphatidylinositol-3,4,5-triphosphate (PIP3) levels and impaired glucose homeostasis in type 1 and type 2 diabetic rats. Cell. Physiol. Biochem. 30, 1363–1370. doi: 10.1159/000343325
Massillon, D. (2001). Regulation of the glucose-6-phosphatase gene by glucose occurs by transcriptional and post-transcriptional mechanisms-Differential effect of glucose and xylitol. J. Biol. Chem. 276, 4055–4062. doi: 10.1074/jbc.M007939200
Miethke, H., Wittig, B., Nath, A., Zierz, S., and Jungermann, K. (1985). Metabolic zonation in liver of diabetic rats. Zonal distribution of phosphoenolpyruvate carboxykinase, pyruvate kinase, glucose-6-phosphatase and succinate dehydrogenase. Biol. Chem. Hoppe. Seyler. 366, 493–501. doi: 10.1515/bchm3.1985.366.1.493
Minchenko, O., Opentanova, I., and Caro, J. (2003). Hypoxic regulation of the 6-phosphofructo-2-kinase/fructose-2,6-bisphosphatase gene family (PFKFB-1-4) expression in vivo. FEBS Lett. 554, 264–270. doi: 10.1016/S0014-5793(03)01179-7
Miralpeix, M., Carballo, E., Bartrons, R., Crepin, K., Hue, L., and Rousseau, G. G. (1992). Oral-administration of vanadate to diabetic rats restores liver 6-phosphofructo-2-kinase content and messenger-Rna. Diabetologia 35, 243–248. doi: 10.1007/BF00400924
Morselt, A. F., Frederiks, W. M., Copius Peereboom-Stegeman, J. H., and van Veen, H. A. (1987). Mechanism of damage to liver cells after chronic exposure to low doses of cadmium chloride. Arch. Toxicol. Suppl. 11, 213–215. doi: 10.1007/978-3-642-72558-6_34
Nauck, M., Wölfle, D., Katz, N., and Jungermann, K. (1981). Modulation of the glucagon-dependent induction of phosphoenolpyruvate carboxykinase and tyrosine aminotransferase by arterial and venous oxygen concentrations in hepatocyte cultures. Eur. J. Biochem. 119, 657–661. doi: 10.1111/j.1432-1033.1981.tb05658.x
Neely, P., El-Maghrabi, M. R., Pilkis, S. J., and Claus, T. H. (1981). Effect of diabetes, insulin, starvation, and refeeding on the level of rat hepatic fructose 2,6-bisphosphate. Diabetes 30, 1062–1064. doi: 10.2337/diab.30.12.1062
Noguchi, T., Inoue, H., and Tanaka, T. (1985). Transcriptional and post-transcriptional regulation of L-type pyruvate-kinase in diabetic rat-liver by insulin and dietary fructose. J. Biol. Chem. 260, 4393–4397.
Postic, C., Burcelin, R., Rencurel, F., Pegorier, J. P., Loizeau, M., Girard, J., et al. (1993). Evidence for a transient inhibitory effect of insulin on GLUT2 expression in the liver: studies in vivo and in vitro. Biochem. J. 293(Pt 1), 119–124.
Probst, I., Schwartz, P., and Jungermann, K. (1982). Induction in primary culture of gluconeogenic and glycolytic hepatocytes resembling periportal and perivenous cells. Eur. J. Biochem. 126, 271–278. doi: 10.1111/j.1432-1033.1982.tb06775.x
Quistorff, B. (1985). Gluconeogenesis in periportal and perivenous hepatocytes of rat liver, isolated by a new high-yield digitonin/collagenase perfusion technique. Biochem. J. 229, 221–226. doi: 10.1042/bj2290221
Rappaport, A. M., Borowy, Z. J., Lougheed, W. M., and Lotto, W. N. (1954). Subdivision of hexagonal liver lobules into a structural and functional unit - role in hepatic physiology and pathology. Anat. Rec. 119, 11–33. doi: 10.1002/ar.1091190103
Schmucker, D. L., Mooney, J. S., and Jones, A. L. (1978). Stereological analysis of hepatic fine structure in the Fischer 344 rat. Influence of sublobular location and animal age. J. Cell Biol. 78, 319–337. doi: 10.1083/jcb.78.2.319
Sibrowski, W., Staegemann, U., and Seitz, H. J. (1982). Accelerated turnover of hepatic glucokinase in starved and streptozotocin-diabetic rat. Eur. J. Biochem. 127, 571–574. doi: 10.1111/j.1432-1033.1982.tb06910.x
Slieker, L. J., Sundell, K. L., Heath, W. F., Osborne, H. E., Bue, J., Manetta, J., et al. (1992). Glucose transporter levels in tissues of spontaneously diabetic Zucker fa/fa rat (ZDF/drt) and viable yellow mouse (Avy/a). Diabetes 41, 187–193. doi: 10.2337/diab.41.2.187
Thorens, B., Flier, J. S., Lodish, H. F., and Kahn, B. B. (1990). Differential regulation of two glucose transporters in rat liver by fasting and refeeding and by diabetes and insulin treatment. Diabetes 39, 712–719. doi: 10.2337/diab.39.6.712
Uyeda, K., and Repa, J. J. (2006). Carbohydrate response element binding protein, ChREBP, a transcription factor coupling hepatic glucose utilization and lipid synthesis. Cell Metab. 4, 107–110. doi: 10.1016/j.cmet.2006.06.008
Van Schaftingen, E., and Hers, H. G. (1983). The role of fructose 2,6-bisphosphate in the long-term control of phosphofructokinase in rat liver. Biochem. Biophys. Res. Commun. 113, 548–554. doi: 10.1016/0006-291X(83)91760-6
Vasilj, A., Gentzel, M., Ueberham, E., Gebhardt, R., and Shevchenko, A. (2012). Tissue proteomics by one-dimensional gel electrophoresis combined with label-free protein quantification. J. Proteome Res. 11, 3680–3689. doi: 10.1021/pr300147z
Weinberg, M. B., and Utter, M. F. (1980). Effect of streptozotocin-induced diabetes mellitus on the turnover of rat liver pyruvate carboxylase and pyruvate dehydrogenase. Biochem. J. 188, 601–608. doi: 10.1042/bj1880601
Weinstein, S. P., O'Boyle, E., Fisher, M., and Haber, R. S. (1994). Regulation of GLUT2 glucose transporter expression in liver by thyroid hormone: evidence for hormonal regulation of the hepatic glucose transport system. Endocrinology 135, 649–654. doi: 10.1210/endo.135.2.8033812
Wimhurst, J. M., and Manchester, K. L. (1970a). A comparison of the effects of diabetes induced with either alloxan or streptozotocin and of starvation on the activities in rat liver of the key enzymes of gluconeogenesis. Biochem. J. 120, 95–103. doi: 10.1042/bj1200095
Wölfle, D., and Jungermann, K. (1985). Long-term effects of physiological oxygen concentrations on glycolysis and gluconeogenesis in hepatocyte cultures. Eur. J. Biochem. 151, 299–303. doi: 10.1111/j.1432-1033.1985.tb09100.x
Wurtman, R. J. (1974). Daily rhythms in tyrosine-transaminase and other hepatic enzymes that metabolize amino-acids-mechanisms and possible consequences. Life Sci. 15, 827–847. doi: 10.1016/0024-3205(74)90001-0
Yang, J., Reshef, L., Cassuto, H., Aleman, G., and Hanson, R. W. (2009). Aspects of the control of phosphoenolpyruvate carboxykinase gene transcription. J. Biol. Chem. 284, 27031–27035. doi: 10.1074/jbc.R109.040535
Keywords: metabolism, metabolic zonation, kinetic model, multiscale model, gene expression
Citation: Berndt N and Holzhütter H-G (2018) Dynamic Metabolic Zonation of the Hepatic Glucose Metabolism Is Accomplished by Sinusoidal Plasma Gradients of Nutrients and Hormones. Front. Physiol. 9:1786. doi: 10.3389/fphys.2018.01786
Received: 16 August 2018; Accepted: 28 November 2018;
Published: 12 December 2018.
Edited by:Steven Dooley, Universitätsmedizin Mannheim, Medizinische Fakultät Mannheim, Universität Heidelberg, Germany
Reviewed by:Adil Mardinoglu, Chalmers University of Technology, Sweden
Rolf Gebhardt, Leipzig University, Germany
Copyright © 2018 Berndt and Holzhütter. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Hermann-Georg Holzhütter, firstname.lastname@example.org