Dynamics and control of the central carbon metabolism in hepatoma cells

Parental sex-dependent effects of either maternal or paternal eNOS deficiency on the offspring's phenotype without transmission of the parental eNOS deficiency to the offspring

Background: Preclinical animal studies and clinical studies indicate that both maternal as well as paternal genetic alterations/gene defects might affect the phenotype of the next-generation without transmissions of the affected gene. Currently, the question of whether the same genetic defect present in the mother or father leads to a similar phenotype in the offspring remains insufficiently elucidated. Methods: In this head-to-head study, we crossbred female and male mice with heterozygous endothelial eNOS knockout (eNOS+/-) with male and female wild-type (wt) mice, respectively. Subsequently, we compared the phenotype of the resulting wt offspring with that of wt offspring born to parents with no eNOS deficiency. Results: Wt female offspring of mothers with heterozygous eNOS showed elevated liver fat accumulation, while wt male offspring of fathers with heterozygous eNOS exhibited increased fasting insulin, heightened insulin levels after a glucose load, and elevated liver glycogen content. By quantitative mass-spectrometry it was shown that concentrations of six serum metabolites (lysoPhosphatidylcholine acyl C20:3, phosphatidylcholine diacyl C36:2, phosphatidylcholine diacyl C38:1, phosphatidylcholine acyl-alkyl C34:1, phosphatidylcholine acyl-alkyl C36:3, and phosphatidylcholine acyl-alkyl C42:5 (PC ae C42:5) as well as four liver carbon metabolites (fructose 6-phosphate, fructose 1,6-bisphosphate, glucose 6-phosphate and fumarate) were different between wt offspring with eNOS+/- mothers and wt offspring with eNOS+/- fathers. Importantly, fumarate was inversely correlated with the liver fat accumulation in female offspring with eNOS+/- mothers and increased liver glycogen in offspring of both sexes with eNOS+/- fathers. The qRT-PCR results revealed that the gene expression patterns were different between wt offspring with eNOS+/- mothers and those offspring with eNOS+/- fathers. Different gene expression patterns were correlated with different observed phenotypic changes in male/female offspring born to mothers or fathers with a heterozygous eNOS genotype. Conclusion: The identical parental genetic alteration (heterozygous eNOS deficiency), without being passed on to the offspring, results in distinct metabolic, liver phenotype, and gene expression pattern variations depending on whether the genetic alteration originated from the father or the mother.

Monitoring and modelling the dynamics of the cellular glycolysis pathway: A review and future perspectives

BACKGROUND

The dynamics of the cellular glycolysis pathway underpin cellular function and dysfunction, and therefore ultimately health, disease, diagnostic and therapeutic strategies. Evolving our understanding of this fundamental process and its dynamics remains critical.

SCOPE OF REVIEW

This paper reviews the medical relevance of glycolytic pathway in depth and explores the current state of the art for monitoring and modelling the dynamics of the process. The future perspectives of label free, vibrational microspectroscopic techniques to overcome the limitations of the current approaches are considered.

MAJOR CONCLUSIONS

Vibrational microspectroscopic techniques can potentially operate in the niche area of limitations of other omics technologies for non-destructive, real-time, in vivo label-free monitoring of glycolysis dynamics at a cellular and subcellular level.

CD147 a direct target of miR-146a supports energy metabolism and promotes tumor growth in ALK+ ALCL

We recently reported that miR-146a is differentially expressed in ALK+ and ALK− anaplastic large cell lymphoma (ALCL). In this study, the downstream targets of miR-146a in ALK+ ALCL were investigated by transcriptome analysis, identifying CD147 as potential target gene. Because CD147 is differentially expressed in ALK+ ALCL versus ALK− ALCL and normal T cells, this gene emerged as a strong candidate for the pathogenesis of this tumor. Here we demonstrate that CD147 is a direct target of miR-146 and contributes to the survival and proliferation of ALK+ ALCL cells in vitro and to the engraftment and tumor growth in vivo in an ALK+ ALCL-xenotransplant mouse model. CD147 knockdown in ALK+ ALCL cells resulted in loss of monocarboxylate transporter 1 (MCT1) expression, reduced glucose consumption and tumor growth retardation, as demonstrated by [¹⁸F]FDG-PET/MRI analysis. Investigation of metabolism in vitro and in vivo supported these findings, revealing reduced aerobic glycolysis and increased basal respiration in CD147 knockdown. In conclusion, our findings indicate that CD147 is of vital importance for ALK+ ALCL to maintain the high energy demand of rapid cell proliferation, promoting lactate export, and tumor growth. Furthermore, CD147 has the potential to serve as a novel therapeutic target in ALK+ ALCL, and warrants further investigation.

Integrating a dynamic central metabolism model of cancer cells with a hybrid 3D multiscale model for vascular hepatocellular carcinoma growth

We develop here a novel modelling approach with the aim of closing the conceptual gap between tumour-level metabolic processes and the metabolic processes occurring in individual cancer cells. In particular, the metabolism in hepatocellular carcinoma derived cell lines (HEPG2 cells) has been well characterized but implementations of multiscale models integrating this known metabolism have not been previously reported. We therefore extend a previously published multiscale model of vascular tumour growth, and integrate it with an experimentally verified network of central metabolism in HEPG2 cells. This resultant combined model links spatially heterogeneous vascular tumour growth with known metabolic networks within tumour cells and accounts for blood flow, angiogenesis, vascular remodelling and nutrient/growth factor transport within a growing tumour, as well as the movement of, and interactions between normal and cancer cells. Model simulations report for the first time, predictions of spatially resolved time courses of core metabolites in HEPG2 cells. These simulations can be performed at a sufficient scale to incorporate clinically relevant features of different tumour systems using reasonable computational resources. Our results predict larger than expected temporal and spatial heterogeneity in the intracellular concentrations of glucose, oxygen, lactate pyruvate, f16bp and Acetyl-CoA. The integrated multiscale model developed here provides an ideal quantitative framework in which to study the relationship between dosage, timing, and scheduling of anti-neoplastic agents and the physiological effects of tumour metabolism at the cellular level. Such models, therefore, have the potential to inform treatment decisions when drug response is dependent on the metabolic state of individual cancer cells.

Spatio-Temporal Multiscale Analysis of Western Diet-Fed Mice Reveals a Translationally Relevant Sequence of Events during NAFLD Progression

Mouse models of non-alcoholic fatty liver disease (NAFLD) are required to define therapeutic targets, but detailed time-resolved studies to establish a sequence of events are lacking. Here, we fed male C57Bl/6N mice a Western or standard diet over 48 weeks. Multiscale time-resolved characterization was performed using RNA-seq, histopathology, immunohistochemistry, intravital imaging, and blood chemistry; the results were compared to human disease. Acetaminophen toxicity and ammonia metabolism were additionally analyzed as functional readouts. We identified a sequence of eight key events: formation of lipid droplets; inflammatory foci; lipogranulomas; zonal reorganization; cell death and replacement proliferation; ductular reaction; fibrogenesis; and hepatocellular cancer. Functional changes included resistance to acetaminophen and altered nitrogen metabolism. The transcriptomic landscape was characterized by two large clusters of monotonously increasing or decreasing genes, and a smaller number of 'rest-and-jump genes' that initially remained unaltered but became differentially expressed only at week 12 or later. Approximately 30% of the genes altered in human NAFLD are also altered in the present mouse model and an increasing overlap with genes altered in human HCC occurred at weeks 30-48. In conclusion, the observed sequence of events recapitulates many features of human disease and offers a basis for the identification of therapeutic targets.

Sex-dependent dynamics of metabolism in primary mouse hepatocytes

The liver is one of the most sexually dimorphic organs. The hepatic metabolic pathways that are subject to sexual dimorphism include xenobiotic, amino acid and lipid metabolism. Non-alcoholic fatty liver disease and hepatocellular carcinoma are among diseases with sex-dependent prevalence, progression and outcome. Although male and female livers differ in their abilities to metabolize foreign compounds, including drugs, sex-dependent treatment and pharmacological dynamics are rarely applied in all relevant cases. Therefore, it is important to consider hepatic sexual dimorphism when developing new treatment strategies and to understand the underlying mechanisms in model systems. We isolated primary hepatocytes from male and female C57BL6/N mice and examined the sex-dependent transcriptome, proteome and extracellular metabolome parameters in the course of culturing them for 96 h. The sex-specific gene expression of the general xenobiotic pathway altered and the female-specific expression of Cyp2b13 and Cyp2b9 was significantly reduced during culture. Sex-dependent differences of several signaling pathways increased, including genes related to serotonin and melatonin degradation. Furthermore, the ratios of male and female gene expression were inversed for other pathways, such as amino acid degradation, beta-oxidation, androgen signaling and hepatic steatosis. Because the primary hepatocytes were cultivated without the influence of known regulators of sexual dimorphism, these results suggest currently unknown modulatory mechanisms of sexual dimorphism in vitro. The large sex-dependent differences in the regulation and dynamics of drug metabolism observed during cultivation can have an immense influence on the evaluation of pharmacodynamic processes when conducting initial preclinical trials to investigate potential new drugs.

Role of mitochondria in liver metabolic health and diseases

The liver is a major organ that coordinates the metabolic flexibility of the whole body, which is characterized by the ability to adapt dynamically in response to fluctuations in energy needs and supplies. In this context, hepatocyte mitochondria are key partners in fine-tuning metabolic flexibility. Here we review the metabolic and signalling pathways carried by mitochondria in the liver, the major pathways that regulate mitochondrial function and how they function in health and metabolic disorders associated to obesity, i.e. insulin resistance, non-alcoholic steatosis and steatohepatitis and hepatocellular carcinoma. Finally, strategies targeting mitochondria to counteract liver disorders are discussed.

Mutual Zonated Interactions of Wnt and Hh Signaling Are Orchestrating the Metabolism of the Adult Liver in Mice and Human

The Hedgehog (Hh) and Wnt/β-Catenin (Wnt) cascades are morphogen pathways whose pronounced influence on adult liver metabolism has been identified in recent years. How both pathways communicate and control liver metabolic functions are largely unknown. Detecting core components of Wnt and Hh signaling and mathematical modeling showed that both pathways in healthy liver act largely complementary to each other in the pericentral (Wnt) and the periportal zone (Hh) and communicate mainly by mutual repression. The Wnt/Hh module inversely controls the spatiotemporal operation of various liver metabolic pathways, as revealed by transcriptome, proteome, and metabolome analyses. Shifting the balance to Wnt (activation) or Hh (inhibition) causes pericentralization and periportalization of liver functions, respectively. Thus, homeostasis of the Wnt/Hh module is essential for maintaining proper liver metabolism and to avoid the development of certain metabolic diseases. With caution due to minor species-specific differences, these conclusions may hold for human liver as well.

Gene Suppression of Transketolase-Like Protein 1 (TKTL1) Sensitizes Glioma Cells to Hypoxia and Ionizing Radiation

In several tumor entities, transketolase-like protein 1 (TKTL1) has been suggested to promote the nonoxidative part of the pentose phosphate pathway (PPP) and thereby to contribute to a malignant phenotype. However, its role in glioma biology has only been sparsely documented. In the present in vitro study using LNT-229 glioma cells, we analyzed the impact of TKTL1 gene suppression on basic metabolic parameters and on survival following oxygen restriction and ionizing radiation. TKTL1 was induced by hypoxia and by hypoxia-inducible factor-1α (HIF-1α). Knockdown of TKTL1 via shRNA increased the cells’ demand for glucose, decreased flux through the PPP and promoted cell death under hypoxic conditions. Following irradiation, suppression of TKTL1 expression resulted in elevated levels of reactive oxygen species (ROS) and reduced clonogenic survival. In summary, our results indicate a role of TKTL1 in the adaptation of tumor cells to oxygen deprivation and in the acquisition of radioresistance. Further studies are necessary to examine whether strategies that antagonize TKTL1 function will be able to restore the sensitivity of glioma cells towards irradiation and antiangiogenic therapies in the more complex in vivo environment.

Mammalian target of rapamycin complex 1 activation sensitizes human glioma cells to hypoxia-induced cell death

Glioblastomas are characterized by fast uncontrolled growth leading to hypoxic areas and necrosis. Signalling from EGFR via mammalian target of rapamycin complex 1 (mTORC1) is a major driver of cell growth and proliferation and one of the most commonly altered signalling pathways in glioblastomas. Therefore, epidermal growth factor receptor and mTORC1 signalling are plausible therapeutic targets and clinical trials with inhibitors are in progress. However, we have previously shown that epidermal growth factor receptor and mTORC1 inhibition triggers metabolic changes leading to adverse effects under the conditions of the tumour microenvironment by protecting from hypoxia-induced cell death. We hypothesized that conversely mTORC1 activation sensitizes glioma cells to hypoxia-induced cell death. As a model for mTORC1 activation we used gene suppression of its physiological inhibitor TSC2 (TSC2sh). TSC2sh glioma cells showed increased sensitivity to hypoxia-induced cell death that was accompanied by an earlier ATP depletion and an increase in reactive oxygen species. There was no difference in extracellular glucose consumption but an altered intracellular metabolic profile with an increase of intermediates of the pentose phosphate pathway. Mechanistically, mTORC1 upregulated the first and rate limiting enzyme of the pentose phosphate pathway, G6PD. Furthermore, an increase in oxygen consumption in TSC2sh cells was detected. This appeared to be due to higher transcription rates of genes involved in mitochondrial respiratory function including PPARGC1A and PPARGC1B (also known as PGC-1α and -β). The finding that mTORC1 activation causes an increase in oxygen consumption and renders malignant glioma cells susceptible to hypoxia and nutrient deprivation could help identify glioblastoma patient cohorts more likely to benefit from hypoxia-inducing therapies such as the VEGFA-targeting antibody bevacizumab in future clinical evaluations.

Maternal eNOS deficiency determines a fatty liver phenotype of the offspring in a sex dependent manner

Maternal environmental factors can impact on the phenotype of the offspring via the induction of epigenetic adaptive mechanisms. The advanced fetal programming hypothesis proposes that maternal genetic variants may influence the offspring's phenotype indirectly via epigenetic modification, despite the absence of a primary genetic defect. To test this hypothesis, heterozygous female eNOS knockout mice and wild type mice were bred with male wild type mice. We then assessed the impact of maternal eNOS deficiency on the liver phenotype of wild type offspring. Birth weight of male wild type offspring born to female heterozygous eNOS knockout mice was reduced compared to offspring of wild type mice. Moreover, the offspring displayed a sex specific liver phenotype, with an increased liver weight, due to steatosis. This was accompanied by sex specific differences in expression and DNA methylation of distinct genes. Liver global DNA methylation was significantly enhanced in both male and female offspring. Also, hepatic parameters of carbohydrate metabolism were reduced in male and female offspring. In addition, male mice displayed reductions in various amino acids in the liver. Maternal genetic alterations, such as partial deletion of the eNOS gene, can affect liver metabolism of wild type offspring without transmission of the intrinsic defect. This occurs in a sex specific way, with more detrimental effects in females. This finding demonstrates that a maternal genetic defect can epigenetically alter the phenotype of the offspring, without inheritance of the defect itself. Importantly, these acquired epigenetic phenotypic changes can persist into adulthood.

RNAi in murine hepatocytes: the agony of choice--a study of the influence of lipid-based transfection reagents on hepatocyte metabolism

Primary hepatocyte cell cultures are widely used for studying hepatic diseases with alterations in hepatic glucose and lipid metabolism, such as diabetes and non-alcoholic fatty liver disease. Therefore, small interfering RNAs (siRNAs) provide a potent and specific tool to elucidate the signaling pathways and gene functions involved in these pathologies. Although RNA interference (RNAi) in vitro is frequently used in these investigations, the metabolic alterations elucidated by different siRNA delivery strategies have hardly been investigated in transfected hepatocytes. To elucidate the influence of the most commonly used lipid-based transfection reagents on cultured primary hepatocytes, we studied the cytotoxic effects and transfection efficiencies of INTERFERin(®), Lipofectamine(®)RNAiMAX, and HiPerFect(®). All of these transfection agents displayed low cytotoxicity (5.6-9.0 ± 1.3-3.4%), normal cell viability, and high transfection efficiency (fold change 0.08-0.13 ± 0.03-0.05), and they also favored the satisfactory down-regulation of target gene expression. However, when effects on the metabolome and lipidome were studied, considerable differences were observed among the transfection reagents. Cellular triacylglycerides levels were either up- or down-regulated [maximum fold change: INTERFERin(®) (48 h) 2.55 ± 0.34, HiPerFect(®) (24 h) 0.79 ± 0.08, Lipofectamine(®)RNAiMAX (48 h) 1.48 ± 0.21], and mRNA levels of genes associated with lipid metabolism were differentially affected. Likewise, metabolic functions such as amino acid utilization from were perturbed (alanine, arginine, glycine, ornithine, and pyruvate). In conclusion, these findings demonstrate that the choice of non-viral siRNA delivery agent is critical in hepatocytes. This should be remembered, especially if RNA silencing is used for studying hepatic lipid homeostasis and its regulation.

Glycolysis is governed by growth regime and simple enzyme regulation in adherent MDCK cells

Due to its vital importance in the supply of cellular pathways with energy and precursors, glycolysis has been studied for several decades regarding its capacity and regulation. For a systems-level understanding of the Madin-Darby canine kidney (MDCK) cell metabolism, we couple a segregated cell growth model published earlier with a structured model of glycolysis, which is based on relatively simple kinetics for enzymatic reactions of glycolysis, to explain the pathway dynamics under various cultivation conditions. The structured model takes into account in vitro enzyme activities, and links glycolysis with pentose phosphate pathway and glycogenesis. Using a single parameterization, metabolite pool dynamics during cell cultivation, glucose limitation and glucose pulse experiments can be consistently reproduced by considering the cultivation history of the cells. Growth phase-dependent glucose uptake together with cell-specific volume changes generate high intracellular metabolite pools and flux rates to satisfy the cellular demand during growth. Under glucose limitation, the coordinated control of glycolytic enzymes re-adjusts the glycolytic flux to prevent the depletion of glycolytic intermediates. Finally, the model's predictive power supports the design of more efficient bioprocesses.

Glycolysis generates biomass precursors and energy from sugars and is therefore a key element in the metabolism of mammalian cells. Changes in its activity greatly affect cellular function which is often recognized as metabolic disease but also as opportunity for the design of efficient bioprocesses. Metabolic research discovered that continuously growing mammalian cells often exhibit a high glycolytic activity but also delivered seemingly endless facets in the pathway operation. The latter call for a systems-level understanding regarding capacity and regulation for a broad range of cultivation conditions. In this work, we couple a cell growth model to a simple kinetic description of glycolysis to consistently explain intracellular metabolite pool dynamics of the Madin-Darby canine kidney cell line over a variety of experiments and time scales while considering the growth status and cultivation history of the cells. We argue that the many different dynamics in glycolysis result from an interplay between a growth-dependent sugar uptake together with simple intrinsic enzyme regulation.

KiMoSys: a web-based repository of experimental data for KInetic MOdels of biological SYStems

BACKGROUND

The kinetic modeling of biological systems is mainly composed of three steps that proceed iteratively: model building, simulation and analysis. In the first step, it is usually required to set initial metabolite concentrations, and to assign kinetic rate laws, along with estimating parameter values using kinetic data through optimization when these are not known. Although the rapid development of high-throughput methods has generated much omics data, experimentalists present only a summary of obtained results for publication, the experimental data files are not usually submitted to any public repository, or simply not available at all. In order to automatize as much as possible the steps of building kinetic models, there is a growing requirement in the systems biology community for easily exchanging data in combination with models, which represents the main motivation of KiMoSys development.

DESCRIPTION

KiMoSys is a user-friendly platform that includes a public data repository of published experimental data, containing concentration data of metabolites and enzymes and flux data. It was designed to ensure data management, storage and sharing for a wider systems biology community. This community repository offers a web-based interface and upload facility to turn available data into publicly accessible, centralized and structured-format data files. Moreover, it compiles and integrates available kinetic models associated with the data.KiMoSys also integrates some tools to facilitate the kinetic model construction process of large-scale metabolic networks, especially when the systems biologists perform computational research.

CONCLUSIONS

KiMoSys is a web-based system that integrates a public data and associated model(s) repository with computational tools, providing the systems biology community with a novel application facilitating data storage and sharing, thus supporting construction of ODE-based kinetic models and collaborative research projects.The web application implemented using Ruby on Rails framework is freely available for web access at http://kimosys.org, along with its full documentation.

Structural control of metabolic flux

Organisms have to continuously adapt to changing environmental conditions or undergo developmental transitions. To meet the accompanying change in metabolic demands, the molecular mechanisms of adaptation involve concerted interactions which ultimately induce a modification of the metabolic state, which is characterized by reaction fluxes and metabolite concentrations. These state transitions are the effect of simultaneously manipulating fluxes through several reactions. While metabolic control analysis has provided a powerful framework for elucidating the principles governing this orchestrated action to understand metabolic control, its applications are restricted by the limited availability of kinetic information. Here, we introduce structural metabolic control as a framework to examine individual reactions' potential to control metabolic functions, such as biomass production, based on structural modeling. The capability to carry out a metabolic function is determined using flux balance analysis (FBA). We examine structural metabolic control on the example of the central carbon metabolism of Escherichia coli by the recently introduced framework of functional centrality (FC). This framework is based on the Shapley value from cooperative game theory and FBA, and we demonstrate its superior ability to assign “share of control” to individual reactions with respect to metabolic functions and environmental conditions. A comparative analysis of various scenarios illustrates the usefulness of FC and its relations to other structural approaches pertaining to metabolic control. We propose a Monte Carlo algorithm to estimate FCs for large networks, based on the enumeration of elementary flux modes. We further give detailed biological interpretation of FCs for production of lactate and ATP under various respiratory conditions.

Insight into the functioning of metabolic control to meet changing demands is a first step in rational engineering of biological systems towards a desired behavior. Metabolic control analysis provides the means to examine the impact of change of reaction fluxes on a specific target flux based on kinetic modeling, but suffers from limitations of the kinetic approach. Here, we introduce and analyze structural metabolic control as a framework to overcome these limitations. We utilize functional centrality, a framework based on the Shapley value from cooperative game theory and flux balance analysis, to determine the contribution of individual reactions to the functions accomplished by a metabolic network. These contributions, in turn, depend on the control exerted on the remaining network. Functional centrality provides the mathematical means to gain further understanding of metabolic control. The potential applications range from facilitating strategies of rational strain design to drug target identification.

Mass isotopomer analysis of metabolically labeled nucleotide sugars and N- and O-glycans for tracing nucleotide sugar metabolisms

Nucleotide sugars are the donor substrates of various glycosyltransferases, and an important building block in N- and O-glycan biosynthesis. Their intercellular concentrations are regulated by cellular metabolic states including diseases such as cancer and diabetes. To investigate the fate of UDP-GlcNAc, we developed a tracing method for UDP-GlcNAc synthesis and use, and GlcNAc utilization using (13)C6-glucose and (13)C2-glucosamine, respectively, followed by the analysis of mass isotopomers using LC-MS. Metabolic labeling of cultured cells with (13)C6-glucose and the analysis of isotopomers of UDP-HexNAc (UDP-GlcNAc plus UDP-GalNAc) and CMP-NeuAc revealed the relative contributions of metabolic pathways leading to UDP-GlcNAc synthesis and use. In pancreatic insulinoma cells, the labeling efficiency of a (13)C6-glucose motif in CMP-NeuAc was lower compared with that in hepatoma cells. Using (13)C2-glucosamine, the diversity of the labeling efficiency was observed in each sugar residue of N- and O-glycans on the basis of isotopomer analysis. In the insulinoma cells, the low labeling efficiencies were found for sialic acids as well as tri- and tetra-sialo N-glycans, whereas asialo N-glycans were found to be abundant. Essentially no significant difference in secreted hyaluronic acids was found among hepatoma and insulinoma cell lines. This indicates that metabolic flows are responsible for the low sialylation in the insulinoma cells. Our strategy should be useful for systematically tracing each stage of cellular GlcNAc metabolism.

Linking metabolism and epigenetic regulation in development of hepatocellular carcinoma

Hepatocellular carcinoma (HCC) is the fifth most common form of cancer globally and is rarely curable once detected. The 5-year survival rate of patients diagnosed with late-stage HCC may be as low as 27%. HCC is a cancer largely driven by epigenetic changes that arise from exposure to exogenous environmental factors rather than coding sequence mutations. The liver is susceptible to effects from Hepatitis C and Hepatitis B viruses, exposure to aflatoxin and continuous excessive consumption of alcohol. The liver is a highly metabolic organ balancing many vital biochemical processes; exposure to any of the above environmental factors is associated with loss of liver function and is a major risk factor for the development of HCC. Emerging studies aim to examine the underlying metabolic processes that are abrogated in cancer and lead to the altered flux and availability of key metabolites important for epigenetic processes. Metabolites have been shown to act as substrates for many canonical epigenetic regulators. These enzymes are responsible for regulating histone modification, DNA methylation and micro RNA expression. By studying the impact of altered liver metabolism, we may better understand the long-term epigenetic processes, which lead to the development and progression of HCC.

Recent advances in 2D and 3D in vitro systems using primary hepatocytes, alternative hepatocyte sources and non-parenchymal liver cells and their use in investigating mechanisms of hepatotoxicity, cell signaling and ADME

This review encompasses the most important advances in liver functions and hepatotoxicity and analyzes which mechanisms can be studied in vitro. In a complex architecture of nested, zonated lobules, the liver consists of approximately 80 % hepatocytes and 20 % non-parenchymal cells, the latter being involved in a secondary phase that may dramatically aggravate the initial damage. Hepatotoxicity, as well as hepatic metabolism, is controlled by a set of nuclear receptors (including PXR, CAR, HNF-4α, FXR, LXR, SHP, VDR and PPAR) and signaling pathways. When isolating liver cells, some pathways are activated, e.g., the RAS/MEK/ERK pathway, whereas others are silenced (e.g. HNF-4α), resulting in up- and downregulation of hundreds of genes. An understanding of these changes is crucial for a correct interpretation of in vitro data. The possibilities and limitations of the most useful liver in vitro systems are summarized, including three-dimensional culture techniques, co-cultures with non-parenchymal cells, hepatospheres, precision cut liver slices and the isolated perfused liver. Also discussed is how closely hepatoma, stem cell and iPS cell-derived hepatocyte-like-cells resemble real hepatocytes. Finally, a summary is given of the state of the art of liver in vitro and mathematical modeling systems that are currently used in the pharmaceutical industry with an emphasis on drug metabolism, prediction of clearance, drug interaction, transporter studies and hepatotoxicity. One key message is that despite our enthusiasm for in vitro systems, we must never lose sight of the in vivo situation. Although hepatocytes have been isolated for decades, the hunt for relevant alternative systems has only just begun.

Metabolomic changes in gastrointestinal tissues after whole body radiation in a murine model

Exposure to ionizing radiation (IR) elicits a set of complex biological responses involving gene expression and protein turnover that ultimately manifest as dysregulation of metabolic processes representing the cellular phenotype. Although radiation biomarkers have been reported in urine and serum, they are not informative about IR mediated tissue or organ specific injury. In the present study we report IR induced metabolic changes in gastrointestinal (GI) tissue of CD2F1 mice using ultra-performance liquid chromatography (UPLC) coupled with electrospray time-of-flight mass spectrometry. Post-radiation GI injury is a critical determinant of survival after exposure to IR. Our results show a distinct dose and time dependent response to GI tissue injury.

The human hepatocyte cell lines IHH and HepaRG: models to study glucose, lipid and lipoprotein metabolism

Metabolic diseases reach epidemic proportions. A better knowledge of the associated alterations in the metabolic pathways in the liver is necessary. These studies need in vitro human cell models. Several human hepatoma models are used, but the response of many metabolic pathways to physiological stimuli is often lost. Here, we characterize two human hepatocyte cell lines, IHH and HepaRG, by analysing the expression and regulation of genes involved in glucose and lipid metabolism. Our results show that the glycolysis pathway is activated by glucose and insulin in both lines. Gluconeogenesis gene expression is induced by forskolin in IHH cells and inhibited by insulin in both cell lines. The lipogenic pathway is regulated by insulin in IHH cells. Finally, both cell lines secrete apolipoprotein B-containing lipoproteins, an effect promoted by increasing glucose concentrations. These two human cell lines are thus interesting models to study the regulation of glucose and lipid metabolism.

A whole-body model for glycogen regulation reveals a critical role for substrate cycling in maintaining blood glucose homeostasis

Timely, and sometimes rapid, metabolic adaptation to changes in food supply is critical for survival as an organism moves from the fasted to the fed state, and vice versa. These transitions necessitate major metabolic changes to maintain energy homeostasis as the source of blood glucose moves away from ingested carbohydrates, through hepatic glycogen stores, towards gluconeogenesis. The integration of hepatic glycogen regulation with extra-hepatic energetics is a key aspect of these adaptive mechanisms. Here we use computational modeling to explore hepatic glycogen regulation under fed and fasting conditions in the context of a whole-body model. The model was validated against previous experimental results concerning glycogen phosphorylase a (active) and glycogen synthase a dynamics. The model qualitatively reproduced physiological changes that occur during transition from the fed to the fasted state. Analysis of the model reveals a critical role for the inhibition of glycogen synthase phosphatase by glycogen phosphorylase a. This negative regulation leads to high levels of glycogen synthase activity during fasting conditions, which in turn increases substrate (futile) cycling, priming the system for a rapid response once an external source of glucose is restored. This work demonstrates that a mechanistic understanding of the design principles used by metabolic control circuits to maintain homeostasis can benefit from the incorporation of mathematical descriptions of these networks into “whole-body” contextual models that mimic in vivo conditions.

Homeostasis of blood glucose concentrations during circadian shifts in survival-related activities, sleep and food availability is crucial for the survival of mammals. This process depends upon glucose intake, short-term storage as glycogen, and gluconeogenesis. The integration of hepatic glycogen anabolic and catabolic dynamics with whole body energetics is critical for survival. In this paper we use computational modeling to investigate the potential survival advantage of substrate (futile) cycling of glycogen and glycogen precursors. Our simulations, combined with published experimental results of other researchers, indicate that as the body enters a state of fasting, the activity of enzymes involved in the synthesis of glycogen increases leading to increased substrate cycling. This increase in substrate cycling allows the system to respond more rapidly once new external sources of glucose become available. The whole-body computational model developed for this work allows the metabolic control circuitry to be studied under simulated in vivo conditions, providing functional insights that are not evident when individual modules of glycogen regulatory circuitry are examined in isolation.

Oncogenic stress induced by acute hyper-activation of Bcr-Abl leads to cell death upon induction of excessive aerobic glycolysis

In response to deregulated oncogene activation, mammalian cells activate disposal programs such as programmed cell death. To investigate the mechanisms behind this oncogenic stress response we used Bcr-Abl over-expressing cells cultivated in presence of imatinib. Imatinib deprivation led to rapid induction of Bcr-Abl activity and over-stimulation of PI3K/Akt-, Ras/MAPK-, and JAK/STAT pathways. This resulted in a delayed necrosis-like cell death starting not before 48 hours after imatinib withdrawal. Cell death was preceded by enhanced glycolysis, glutaminolysis, and amino acid metabolism leading to elevated ATP and protein levels. This enhanced metabolism could be linked to induction of cell death as inhibition of glycolysis or glutaminolysis was sufficient to sustain cell viability. Therefore, these data provide first evidence that metabolic changes induced by Bcr-Abl hyper-activation are important mediators of oncogenic stress-induced cell death.

During the first 30 hours after imatinib deprivation, Bcr-Abl hyper-activation did not affect proliferation but resulted in cellular swelling, vacuolization, and induction of eIF2α phosphorylation, CHOP expression, as well as alternative splicing of XPB, indicating endoplasmic reticulum stress response. Cell death was dependent on p38 and RIP1 signaling, whereas classical death effectors of ER stress, namely CHOP-BIM were antagonized by concomitant up-regulation of Bcl-xL.

Screening of 1,120 compounds for their potential effects on oncogenic stress-induced cell death uncovered that corticosteroids antagonize cell death upon Bcr-Abl hyper-activation by normalizing cellular metabolism. This protective effect is further demonstrated by the finding that corticosteroids rendered lymphocytes permissive to the transforming activity of Bcr-Abl. As corticosteroids are used together with imatinib for treatment of Bcr-Abl positive acute lymphoblastic leukemia these data could have important implications for the design of combination therapy protocols.

In conclusion, excessive induction of Warburg type metabolic alterations can cause cell death. Our data indicate that these metabolic changes are major mediators of oncogenic stress induced by Bcr-Abl.

Regulation and control of metabolic fluxes in microbes

After about ten years of research renaissance in metabolism, the present challenge is to understand how metabolic fluxes are controlled by a complex interplay of overlapping regulatory mechanisms. Reconstruction of various regulatory network topologies is steaming, illustrating that we underestimated the broad importance of post-translational modifications such as enzyme phosphorylation or acetylation for microbial metabolism. With the growing topological knowledge, the functional relevance of these regulatory events becomes an even more pressing need. A major knowledge gap resides in the regulatory network of protein-metabolite interactions, simply because we lacked pertinent methods for systematic analyses - but a start has now been made. Perhaps most dramatic was the conceptual shift in our perception of metabolism from an engine of cellular operation to a generator of input and feedback signals for regulatory circuits that govern many important decisions on cell proliferation, differentiation, death, and naturally metabolism.