Search

Abstract Non-invasive breath tests can serve as valuable diagnostic tools in medicine as they can determine particular enzymatic and metabolic functions in vivo. However, methodological pitfalls have limited the actual clinical application of those tests till today. A major challenge of non-invasive breath tests has remained the provision of individually reliable test results. To overcome these limitations, a better understanding of breath kinetics during non-invasive breaths tests is essential. This analysis compares the breath recovery of a 13C-methacetin breath test with the actual serum kinetics of the substrate. It is shown, that breath and serum kinetics of the same test are significantly different over a period of 60 minutes. The recovery of the tracer 13CO2 in breath seems to be significantly delayed due to intermediate storage in the bicarbonate pool. This has to be taken into account for the application of non-invasive breath test protocols. Otherwise, breath tests might display bicarbonate kinetics despite the metabolic capacity of the particular target enzyme.

Cells store triacylglycerol (TAG) within lipid droplets (LDs). A dynamic model describing complete LD formation at the endoplasmic reticulum (ER) membrane does not yet exist. A biochemical-biophysical model of LD synthesis is proposed. It describes the time-dependent accumulation of TAG in the ER membrane as the formation of a potential LD (pLD) bounded by spherical caps of the inner and outer monolayers of the membrane. The expansion rate of the pLD depends on the TAG supply, the elastic properties of the ER membrane, and the recruitment of phospholipids (PLs) to the cap-covering monolayers. Model simulations provided the following insights: (a) Marginal differences in the surface tension of the cap monolayers are sufficient to fully drive the expansion of the pLD towards the cytosol or lumen. (b) Selective reduction of PL supply to the luminal monolayer ensures stable formation of cytosolic LDs, irrespective of variations in the elasto-mechanical properties of the ER membrane. (c) The rate of TAG supply to the cytosolic monolayer has a major effect on the size and maturation time of LDs but has no significant effect on the TAG export per individual LD. The recruitment of additional PLs to the cap monolayers of pLDs critically controls the budding direction, size, and maturation time of LDs. The ability of cells to acquire additional LD initiation sites appears to be key to coping with acutely high levels of potentially toxic free fatty acids., (© 2024 The Author(s). The FEBS Journal published by John Wiley & Sons Ltd on behalf of Federation of European Biochemical Societies.)

Non-alcoholic fatty liver disease (NAFLD) is the most common chronic liver disease in children and is associated with overweight and insulin resistance (IR). Almost nothing is known about in vivo alterations of liver metabolism in NAFLD, especially in the early stages of non-alcoholic steatohepatitis (NASH). Here, we used a complex mathematical model of liver metabolism to quantify the central hepatic metabolic functions of 71 children with biopsy-proven NAFLD. For each patient, a personalized model variant was generated based on enzyme abundances determined by mass spectroscopy. Our analysis revealed statistically significant alterations in the hepatic carbohydrate, lipid, and ammonia metabolism, which increased with the degree of obesity and severity of NAFLD. Histologic features of NASH and IR displayed opposing associations with changes in carbohydrate and lipid metabolism but synergistically decreased urea synthesis in favor of the increased release of glutamine, a driver of liver fibrosis. Taken together, our study reveals already significant alterations in the NASH liver of pediatric patients, which, however, are differently modulated by the simultaneous presence of IR.

Background: Many heart diseases can result in reduced pumping capacity of the heart muscle. A mismatch between ATP demand and ATP production of cardiomyocytes is one of the possible causes. Assessment of the relation between myocardial ATP production (MV ATP ) and cardiac workload is important for better understanding disease development and choice of nutritional or pharmacologic treatment strategies. Because there is no method for measuring MV ATP in vivo, the use of physiology-based metabolic models in conjunction with protein abundance data is an attractive approach., Method: We developed a comprehensive kinetic model of cardiac energy metabolism (CARDIOKIN1) that recapitulates numerous experimental findings on cardiac metabolism obtained with isolated cardiomyocytes, perfused animal hearts, and in vivo studies with humans. We used the model to assess the energy status of the left ventricle of healthy participants and patients with aortic stenosis and mitral valve insufficiency. Maximal enzyme activities were individually scaled by means of protein abundances in left ventricle tissue samples. The energy status of the left ventricle was quantified by the ATP consumption at rest (MV ATP [rest]), at maximal workload (MV ATP [max]), and by the myocardial ATP production reserve, representing the span between MV ATP (rest) and MV ATP (max)., Results: Compared with controls, in both groups of patients, MV ATP (rest) was increased and MV ATP (max) was decreased, resulting in a decreased myocardial ATP production reserve, although all patients had preserved ejection fraction. The variance of the energetic status was high, ranging from decreased to normal values. In both patient groups, the energetic status was tightly associated with mechanic energy demand. A decrease of MV ATP (max) was associated with a decrease of the cardiac output, indicating that cardiac functionality and energetic performance of the ventricle are closely coupled., Conclusions: Our analysis suggests that the ATP-producing capacity of the left ventricle of patients with valvular dysfunction is generally diminished and correlates positively with mechanical energy demand and cardiac output. However, large differences exist in the energetic state of the myocardium even in patients with similar clinical or image-based markers of hypertrophy and pump function. Registration: URL: https://www.clinicaltrials.gov; Unique identifiers: NCT03172338 and NCT04068740.

Deep anaesthesia may impair neuronal, vascular and mitochondrial function facilitating neurological complications, such as delirium and stroke. On the other hand, deep anaesthesia is performed for neuroprotection in critical brain diseases such as status epilepticus or traumatic brain injury. Since the commonly used anaesthetic propofol causes mitochondrial dysfunction, we investigated the impact of the alternative anaesthetic isoflurane on neuro-metabolism. In deeply anaesthetised Wistar rats (burst suppression pattern), we measured increased cortical tissue oxygen pressure (p ti O 2 ), a ∼35% drop in regional cerebral blood flow (rCBF) and burst-associated neurovascular responses. In vitro , 3% isoflurane blocked synaptic transmission and impaired network oscillations, thereby decreasing the cerebral metabolic rate of oxygen (CMRO 2 ). Concerning mitochondrial function, isoflurane induced a reductive shift in flavin adenine dinucleotide (FAD) and decreased stimulus-induced FAD transients as Ca 2+ influx was reduced by ∼50%. Computer simulations based on experimental results predicted no direct effects of isoflurane on mitochondrial complexes or ATP-synthesis. We found that isoflurane-induced burst suppression is related to decreased ATP consumption due to inhibition of synaptic activity while neurovascular coupling and mitochondrial function remain intact. The neurometabolic profile of isoflurane thus appears to be superior to that of propofol which has been shown to impair the mitochondrial respiratory chain.

Structural changes of soft tissues on the cellular level can be characterized by histopathology, but not longitudinally in the same tissue. Alterations of cellular structures and tissue matrix are associated with changes in biophysical properties which can be monitored longitudinally by quantitative diffusion-weighted imaging (DWI) and magnetic resonance elastography (MRE). In this work, DWI and MRE examinations were performed in a 0.5-Tesla compact scanner to investigate longitudinal changes in water diffusivity, stiffness and viscosity of ex-vivo rat livers for up to 20 h post-mortem (pm). The effect of blood on biophysical parameters was examined in 13 non-perfused livers (containing blood, NPLs) and 14 perfused livers (blood washed out, PLs). Changes in cell shape, cell packing and cell wall integrity were characterized histologically. In all acquisitions, NPLs presented with higher shear-wave speed (c) , higher shear-wave penetration rate (a) and smaller apparent-diffusion-coefficients (ADCs) than PL. Time-resolved analysis revealed three distinct phases: (i) an initial phase (up to 2 h pm) with markedly increased c and a and reduced ADCs; (ii) an extended phase with relatively stable values; and (iii) a degradation phase characterized by significant increases in a (10 h pm in NPLs and PLs) and ADCs (10 h pm in NPLs, 13 h pm in PLs). Histology revealed changes in cell shape and packing along with decreased cell wall integrity, indicating tissue degradation in NPLs and PLs 10 h pm. Taken together, our results demonstrate that the biophysical properties of fresh liver tissue rapidly change within 2 h pm, which seems to be an effect of both cytotoxic edema and vascular blood content. Several hours later, disruption of cell walls resulted in higher water diffusivity and wave penetration. These results reveal the individual contributions of vascular components and cellular integrity to liver elastography and provide a biophysical, imaging-based fingerprint of liver tissue degradation., Competing Interests: 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., (Copyright © 2021 Garczyńska, Tzschätzsch, Assili, Kühl, Häckel, Schellenberger, Berndt, Holzhütter, Braun, Sack and Guo.)

Proteasome-catalyzed peptide splicing (PCPS) of cancer-driving antigens could generate attractive neoepitopes to be targeted by T cell receptor (TCR)-based adoptive T cell therapy. Based on a spliced peptide prediction algorithm, TCRs were generated against putative KRAS G12V - and RAC2 P29L -derived neo-splicetopes with high HLA-A*02:01 binding affinity. TCRs generated in mice with a diverse human TCR repertoire specifically recognized the respective target peptides with high efficacy. However, we failed to detect any neo-splicetope-specific T cell response when testing the in vivo neo-splicetope generation and obtained no experimental evidence that the putative KRAS G12V - and RAC2 P29L -derived neo-splicetopes were naturally processed and presented. Furthermore, only the putative RAC2 P29L -derived neo-splicetopes was generated by in vitro PCPS. The experiments pose severe questions on the notion that available algorithms or the in vitro PCPS reaction reliably simulate in vivo splicing and argue against the general applicability of an algorithm-driven 'reverse immunology' pipeline for the identification of cancer-specific neo-splicetopes., Competing Interests: GW, CB, LI, GP, VS, AG, HH, TB, PK No competing interests declared, (© 2021, Willimsky et al.)

Metabolic reprogramming is a characteristic feature of cancer cells, but there is no unique metabolic program for all tumors. Genetic and gene expression studies have revealed heterogeneous inter- and intratumor patterns of metabolic enzymes and membrane transporters. The functional implications of this heterogeneity remain often elusive. Here, we applied a systems biology approach to gain a comprehensive and quantitative picture of metabolic changes in individual hepatocellular carcinoma (HCC). We used protein intensity profiles determined by mass spectrometry in samples of 10 human HCCs and the adjacent noncancerous tissue to calibrate Hepatokin1, a complex mathematical model of liver metabolism. We computed the 24-h profile of 18 metabolic functions related to carbohydrate, lipid, and nitrogen metabolism. There was a general tendency among the tumors toward downregulated glucose uptake and glucose release albeit with large intertumor variability. This finding calls into question that the Warburg effect dictates the metabolic phenotype of HCC. All tumors comprised elevated β-oxidation rates. Urea synthesis was found to be consistently downregulated but without compromising the tumor's capacity for ammonia detoxification owing to increased glutamine synthesis. The largest intertumor heterogeneity was found for the uptake and release of lactate and the size of the cellular glycogen content. In line with the observed metabolic heterogeneity, the individual HCCs differed largely in their vulnerability against pharmacological treatment with metformin. Taken together, our approach provided a comprehensive and quantitative characterization of HCC metabolism that may pave the way for a computational a priori assessment of pharmacological therapies targeting metabolic processes of HCC., (© 2020 The Authors. The FEBS Journal published by John Wiley & Sons Ltd on behalf of Federation of European Biochemical Societies.)

Nonalcoholic fatty liver disease (NAFLD) is the most common chronic liver disease in children and adolescents. About 30% of patients with NAFLD progress to the more severe condition of nonalcoholic steatohepatitis (NASH), which is typically diagnosed using liver biopsy. Liver stiffness (LS) quantified by elastography is a promising imaging marker for the noninvasive assessment of NAFLD and NASH in pediatric patients. However, the link between LS and specific histopathologic features used for clinical staging of NAFLD is not well defined. Furthermore, LS data reported in the literature can vary greatly due to the use of different measurement techniques. Uniquely, time-harmonic elastography (THE) based on ultrasound and magnetic resonance elastography (MRE) use the same mechanical stimulation, allowing us to compare LS in biopsy-proven NAFLD previously determined by THE and MRE in 67 and 50 adolescents, respectively. In the present work, we analyzed the influence of seven distinct histopathologic features on LS, including septal infiltration, bridging fibrosis, pericellular fibrosis, hepatocellular ballooning, portal inflammation, lobular inflammation, and steatosis. LS was highly correlated with periportal and lobular fibrosis as well as hepatocellular ballooning while no independent association was found for inflammation and steatosis. Based on this analysis, we propose a composite elastography score (CES) which includes the four key histopathologic features identified as mechanically relevant. Interestingly, CES-relevant histopathologic features were associated with zonal distribution patterns of pediatric NAFLD. Mechano-structural changes associated with NAFLD progression can be histopathologically staged using the CES, which is easily determined noninvasively based on LS measured by time-harmonic elastography., Competing Interests: Declaration of Competing Interest The authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper., (Copyright © 2021 Acta Materialia Inc. Published by Elsevier Ltd. All rights reserved.)

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

Background and Aims: Zone-dependent differences in expression of metabolic enzymes along the portocentral axis of the acinus are a long-known feature of liver metabolism. A prominent example is the preferential localization of the enzyme, glutamine synthetase, in pericentral hepatocytes, where it converts potentially toxic ammonia to the valuable amino acid, glutamine. However, with the exception of a few key regulatory enzymes, a comprehensive and quantitative assessment of zonal differences in the abundance of metabolic enzymes and, much more important, an estimation of the associated functional differences between portal and central hepatocytes is missing thus far., Approach and Results: We addressed this problem by establishing a method for the separation of periportal and pericentral hepatocytes that yields sufficiently pure fractions of both cell populations. Quantitative shotgun proteomics identified hundreds of differentially expressed enzymes in the two cell populations. We used zone-specific proteomics data for scaling of the maximal activities to generate portal and central instantiations of a comprehensive kinetic model of central hepatic metabolism (Hepatokin1)., Conclusions: The model simulations revealed significant portal-to-central differences in almost all metabolic pathways involving carbohydrates, fatty acids, amino acids, and detoxification., (© 2020 The Authors. Hepatology published by Wiley Periodicals LLC on behalf of American Association for the Study of Liver Diseases.)

Non-alcoholic fatty liver disease (NAFLD) is a significant health burden in obese children for which there is currently no specific therapy. Preclinical studies indicate that epoxyeicosanoids, a class of bioactive lipid mediators that are generated by cytochrome P450 (CYP) epoxygenases and inactivated by the soluble epoxide hydrolase (sEH), play a protective role in NAFLD. We performed a comprehensive lipidomics analysis using liver tissue and blood samples of 40 children with NAFLD. Proteomics was performed to determine CYP epoxygenase and sEH expressions. Hepatic epoxyeicosanoids significantly increased with higher grades of steatosis, while their precursor PUFAs were unaltered. Concomitantly, total CYP epoxygenase activity increased while protein level and activity of sEH decreased. In contrast, hepatic epoxyeicosanoids showed a strong decreasing trend with higher stages of fibrosis, accompanied by a decrease of CYP epoxygenase activity and protein expression. These findings suggest that the CYP epoxygenase/sEH pathway represents a potential pharmacologic target for the treatment of NAFLD., (Copyright © 2020 Elsevier Ltd. All rights reserved.)

During pregnancy, the body's hyperestrogenic state alters hepatic metabolism and synthesis. While biochemical changes related to liver function during normal pregnancy are well understood, pregnancy-associated alterations in biophysical properties of the liver remain elusive. In this study, we investigated 26 ex vivo fresh liver specimens harvested from pregnant and non-pregnant rats by diffusion-weighted imaging (DWI) and magnetic resonance elastography (MRE) in a 0.5-Tesla compact magnetic resonance imaging (MRI) scanner. Water diffusivity and viscoelastic parameters were compared with histological data and blood markers. We found livers from pregnant rats to have (i) significantly enlarged hepatocytes (26 ± 15%, p < 0.001), (ii) increased liver stiffness (12 ± 15%, p = 0.012), (iii) decreased viscosity (-23 ± 14%, p < 0.001), and (iv) increased water diffusivity (12 ± 11%, p < 0.001). In conclusion, increased stiffness and reduced viscosity of the liver during pregnancy are mainly attributable to hepatocyte enlargement. Hypertrophy of liver cells imposes fewer restrictions on intracellular water mobility, resulting in a higher hepatic water diffusion coefficient. Collectively, MRE and DWI have the potential to inform on structural liver changes associated with pregnancy in a clinical context., (Copyright © 2020 Garczyńska, Tzschätzsch, Kühl, Morr, Lilaj, Häckel, Schellenberger, Berndt, Holzhütter, Braun, Sack and Guo.)

The principle of dynamic liver function breath tests is founded on the administration of a 13 C-labeled drug and subsequent monitoring of 13 CO 2 in the breath, quantified as time series delta over natural baseline 13 CO 2 (DOB) liberated from the drug during hepatic CYP-dependent detoxification. One confounding factor limiting the diagnostic value of such tests is that only a fraction of the liberated 13 CO 2 is immediately exhaled, while another fraction is taken up by body compartments from which it returns with delay to the plasma. The aims of this study were to establish a novel variant of the methacetin-based breath test LiMAx that allows to estimate and to eliminate the confounding effect of systemic 13 CO 2 distribution on the DOB curve and thus enables a more reliable assessment of the hepatic detoxification capacity compared with the conventional LiMAx test. We designed a new test variant (named "2DOB") consisting of two consecutive phases. Phase 1 is initiated by the intravenous administration of 13 C-bicarbonate. Phase 2 starts about 30 min later with the intravenous administration of the 13 C-labelled test drug. Using compartment modelling, the resulting 2-phasic DOB curve yields the rate constants for the irreversible elimination and the reversible exchange of plasma 13 CO 2 with body compartments (phase 1) and for the detoxification and exchange of the drug with body compartments (phase 2). We carried out the 2DOB test with the test drug 13 C-methacetin in 16 subjects with chronic liver pathologies and 22 normal subjects, who also underwent the conventional LiMAx test. Individual differences in the systemic CO 2 kinetics can lead to deviations up to a factor of 2 in the maximum of DOB curves (coefficient of variation CV ≈ 0.2) which, in particular, may hamper the discrimination between subjects with normal or mildly impaired detoxification capacities. The novel test revealed that a significant portion of the drug is not immediately metabolized, but transiently taken up into a storage compartment. Intriguingly, not only the hepatic detoxification rate but also the storage capacity of the drug, turned out to be indicative for a normal liver function. We thus used both parameters to define a scoring function which yielded an excellent disease classification (AUC = 0.95) and a high correlation with the MELD score (R Spearman  = 0.92). The novel test variant 2DOB promises a significant improvement in the assessment of impaired hepatic detoxification capacity. The suitability of the test for the reliable characterization of the natural history of chronic liver diseases (fatty liver-fibrosis-cirrhosis) has to be assessed in further studies.

Background: Metabolic alterations can serve as targets for diagnosis and cancer therapy. Due to the highly complex regulation of cellular metabolism, definite identification of metabolic pathway alterations remains challenging and requires sophisticated experimentation., Methods: We applied a comprehensive kinetic model of the central carbon metabolism (CCM) to characterise metabolic reprogramming in murine liver cancer., Results: We show that relative differences of protein abundances of metabolic enzymes obtained by mass spectrometry can be used to assess their maximal velocity values. Model simulations predicted tumour-specific alterations of various components of the CCM, a selected number of which were subsequently verified by in vitro and in vivo experiments. Furthermore, we demonstrate the ability of the kinetic model to identify metabolic pathways whose inhibition results in selective tumour cell killing., Conclusions: Our systems biology approach establishes that combining cellular experimentation with computer simulations of physiology-based metabolic models enables a comprehensive understanding of deregulated energetics in cancer. We propose that modelling proteomics data from human HCC with our approach will enable an individualised metabolic profiling of tumours and predictions of the efficacy of drug therapies targeting specific metabolic pathways.

In developing neurons, phosphoinositide 3-kinases (PI3Ks) control axon growth and branching by positively regulating PI3K/PI(3,4,5)P 3 , but how neurons are able to generate sufficient PI(3,4,5)P 3 in the presence of high levels of the antagonizing phosphatase PTEN is difficult to reconcile. We find that normal axon morphogenesis involves homeostasis of elongation and branch growth controlled by accumulation of PI(3,4,5)P 3 through PTEN inhibition. We identify a plasma membrane-localized protein-protein interaction of PTEN with plasticity-related gene 2 (PRG2). PRG2 stabilizes membrane PI(3,4,5)P 3 by inhibiting PTEN and localizes in nanoclusters along axon membranes when neurons initiate their complex branching behavior. We demonstrate that PRG2 is both sufficient and necessary to account for the ability of neurons to generate axon filopodia and branches in dependence on PI3K/PI(3,4,5)P 3 and PTEN. Our data indicate that PRG2 is part of a neuronal growth program that induces collateral branch growth in axons by conferring local inhibition of PTEN., (Copyright © 2019 The Author(s). Published by Elsevier Inc. All rights reserved.)

Human hepatocellular carcinoma (HCC) is the most common type of primary liver cancer in adults and the most common cause of death in people with cirrhosis. While previous metabolic studies of HCC have mainly focused on the glucose metabolism (Warburg effect), less attention has been paid to tumor-specific features of the lipid metabolism. Here, we applied a computational approach to analyze major pathways of fatty acid utilization in individual HCC. To this end, we used protein intensity profiles of eleven human HCCs to parameterize tumor-specific kinetic models of cellular lipid metabolism including formation, enlargement, and degradation of lipid droplets (LDs). Our analysis reveals significant inter-tumor differences in the lipid metabolism. The majority of HCCs show a reduced uptake of fatty acids and decreased rate of β-oxidation, however, some HCCs display a completely different metabolic phenotype characterized by high rates of β-oxidation. Despite reduced fatty acid uptake in the majority of HCCs, the content of triacylglycerol is significantly enlarged compared to the tumor-adjacent tissue. This is due to tumor-specific expression profiles of regulatory proteins decorating the surface of LDs and controlling their turnover. Our simulations suggest that HCCs characterized by a very high content of triglycerides comprise regulatory peculiarities that render them susceptible to selective drug targeting without affecting healthy tissue., Competing Interests: The authors declare no conflict of interest.

Cortical information processing comprises various activity states emerging from timed synaptic excitation and inhibition. However, the underlying energy metabolism is widely unknown. We determined the cerebral metabolic rate of oxygen (CMRO 2 ) along a tissue depth of <0.3 mm in the hippocampal CA3 region during various network activities, including gamma oscillations and sharp wave-ripples that occur during wakefulness and sleep. These physiological states associate with sensory perception and memory formation, and critically depend on perisomatic GABA inhibition. Moreover, we modelled vascular oxygen delivery based on quantitative microvasculature analysis. (1) Local CMRO 2 was highest during gamma oscillations (3.4 mM/min), medium during sharp wave-ripples, asynchronous activity and isoflurane application (2.0-1.6 mM/min), and lowest during tetrodotoxin application (1.4 mM/min). (2) Energy expenditure of axonal and synaptic signaling accounted for >50% during gamma oscillations. (3) CMRO 2 positively correlated with number and synchronisation of activated synapses, and neural multi-unit activity. (4) The median capillary distance was 44 µm. (5) The vascular oxygen partial pressure of 33 mmHg was needed to sustain oxidative phosphorylation during gamma oscillations. We conclude that gamma oscillations featuring high energetics require a hemodynamic response to match oxygen consumption of respiring mitochondria, and that perisomatic inhibition significantly contributes to the brain energy budget.

Objectives: Today, nonalcoholic fatty liver disease (NAFLD) is the most common chronic liver disease in children and adults alike. Yet, the noninvasive evaluation of disease severity remains a diagnostic challenge. In this study, we apply multifrequency magnetic resonance elastography (mMRE) for the quantification of liver steatosis and fibrosis in adolescents with NAFLD., Methods: Fifty adolescents (age range, 10-17 years; mean BMI, 33.9 kg/m; range, 21.4-42.1 kg/m) with biopsy-proven NAFLD were included in this prospective study. Multifrequency magnetic resonance elastography was performed using external multifrequency vibrations of 30 to 60 Hz and tomoelastography postprocessing, resulting in penetration rate (a) and shear wave speed (c). Hepatic fat fraction was determined using Dixon method. The diagnostic accuracy of mMRE in grading liver steatosis and staging liver fibrosis was assessed by receiver operating characteristic curve analysis., Results: Multifrequency magnetic resonance elastography parameters c and a were independently sensitive to fibrosis and steatosis, respectively, providing area under the receiver operating characteristic values of 0.79 (95% confidence interval [CI], 0.66-0.92), 0.91 (95% CI, 0.83-0.99), and 0.90 (95% CI, 0.80-0.99) for the detection of any (≥F1), moderate (≥F2), and advanced (≥F3) fibrosis, and 0.87 (95% CI, 0.76-0.97) and 0.87 (95% CI, 0.77-0.96) for the detection of moderate (≥S2) and severe (S3) steatosis., Conclusions: One mMRE measurement provides 2 independent parameters with very good diagnostic accuracy in detecting moderate and advanced fibrosis as well as moderate and severe steatosis in pediatric NAFLD.

Background and Aims: Non-alcoholic fatty liver disease (NAFLD) is the most common chronic liver disease in children and adolescents today. In comparison with adult disease, paediatric NAFLD may show a periportal localization, which is associated with advanced fibrosis. This study aimed to assess the role of genetic risk variants for histological disease pattern and severity in childhood NAFLD., Methods: We studied 14 single nucleotide polymorphisms (SNP) in a cohort of 70 adolescents with biopsy-proven NAFLD. Genotype was compared to an adult control cohort (n = 200) and analysed in relation to histological disease severity and liver tissue proteomics., Results: Three of the 14 SNPs were significantly associated with paediatric NAFLD after FDR adjustment, rs738409 (PNPLA3, P = 2.80 × 10 -06 ), rs1044498 (ENPP1, P = 0.0091) and rs780094 (GCKR, P = 0.0281). The severity of steatosis was critically associated with rs738409 (OR=3.25; 95% CI: 1.72-6.52, FDR-adjusted P = 0.0070). The strongest variants associated with severity of fibrosis were rs1260326, rs780094 (both GCKR) and rs659366 (UCP2). PNPLA3 was associated with a portal pattern of steatosis, inflammation and fibrosis. Proteome profiling revealed decreasing levels of GCKR protein with increasing carriage of the rs1260326/rs780094 minor alleles and downregulation of the retinol pathway in rs738409 G/G carriers. Computational metabolic modelling highlighted functional relevance of PNPLA3, GCKR and UCP2 for NAFLD development., Conclusions: This study provides evidence for the role of PNPLA3 as a determinant of portal NAFLD localization and severity of portal fibrosis in children and adolescents, the risk variant being associated with an impaired hepatic retinol metabolism., (© 2018 John Wiley & Sons A/S. Published by John Wiley & Sons Ltd.)

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.

Propofol is the most frequently used intravenous anesthetic for induction and maintenance of anesthesia. Propofol acts first and formost as a GABA A -agonist, but effects on other neuronal receptors and voltage-gated ion channels have been described. Besides its direct effect on neurotransmission, propofol-dependent impairment of mitochondrial function in neurons has been suggested to be responsible for neurotoxicity and postoperative brain dysfunction. To clarify the potential neurotoxic effect in more detail, we investigated the effects of propofol on neuronal energy metabolism of hippocampal slices of the stratum pyramidale of area CA3 at different activity states. We combined oxygen-measurements, electrophysiology and flavin adenine dinucleotide (FAD)-imaging with computational modeling to uncover molecular targets in mitochondrial energy metabolism that are directly inhibited by propofol. We found that high concentrations of propofol (100 µM) significantly decrease population spikes, paired pulse ratio, the cerebral metabolic rate of oxygen consumption (CMRO 2 ), frequency and power of gamma oscillations and increase FAD-oxidation. Model-based simulation of mitochondrial FAD redox state at inhibition of different respiratory chain (RC) complexes and the pyruvate-dehydrogenase show that the alterations in FAD-autofluorescence during propofol administration can be explained with a strong direct inhibition of the complex II (cxII) of the RC. While this inhibition may not affect ATP availability under normal conditions, it may have an impact at high energy demand. Our data support the notion that propofol may lead to neurotoxicity and neuronal dysfunction by directly affecting the energy metabolism in neurons.

Purpose To measure in vivo liver stiffness by using US time-harmonic elastography in a cohort of pediatric patients who were overweight to extremely obese with nonalcoholic fatty liver disease (NAFLD) and to evaluate the diagnostic value of time-harmonic elastography for differentiating stages of fibrosis associated with progressive disease. Materials and Methods In this prospective study, 67 consecutive adolescents (age range, 10-17 years; mean body mass index, 34.7 kg/m 2 ; range, 21.4-50.4 kg/m 2 ) with biopsy-proven NAFLD were enrolled. Liver stiffness was measured by using time-harmonic elastography based on externally induced continuous vibrations of 30 Hz to 60 Hz frequency and real-time B-mode-guided wave profile analysis covering tissue depths of up to 14 cm. The diagnostic accuracy of time-harmonic elastography in staging liver fibrosis was assessed with area under the receiver operating characteristic curve (AUC) analysis. Liver stiffness cutoffs for the differentiation of fibrosis stages were identified based on the highest Youden index. Results Time-harmonic elastography was feasible in all patients (0% failure rate), including 70% (n = 47) of individuals with extreme obesity (body mass index above the 99.5th percentile). AUC analysis for the detection of any fibrosis (≥ stage F1), moderate fibrosis (≥ stage F2), and advanced fibrosis (≥ stage F3) was 0.88 (95% confidence interval [CI]: 0.80, 0.96), 0.99 (95% CI: 0.98, 1.00), and 0.88 (95% CI: 0.80, 0.96), respectively. The best liver stiffness cutoffs were 1.52 m/sec for at least stage F1, 1.62 m/sec for at least stage F2, and 1.64 m/sec for at least stage F3. Conclusion US time-harmonic elastography allows accurate detection of moderate fibrosis even in pediatric patients with extreme obesity. Larger clinical trials are warranted to confirm the accuracy of US time-harmonic elastography., (© RSNA, 2018 Online supplemental material is available for this article.)

The epidemic increase of non-alcoholic fatty liver diseases (NAFLD) requires a deeper understanding of the regulatory circuits controlling the response of liver metabolism to nutritional challenges, medical drugs, and genetic enzyme variants. As in vivo studies of human liver metabolism are encumbered with serious ethical and technical issues, we developed a comprehensive biochemistry-based kinetic model of the central liver metabolism including the regulation of enzyme activities by their reactants, allosteric effectors, and hormone-dependent phosphorylation. The utility of the model for basic research and applications in medicine and pharmacology is illustrated by simulating diurnal variations of the metabolic state of the liver at various perturbations caused by nutritional challenges (alcohol), drugs (valproate), and inherited enzyme disorders (galactosemia). Using proteomics data to scale maximal enzyme activities, the model is used to highlight differences in the metabolic functions of normal hepatocytes and malignant liver cells (adenoma and hepatocellular carcinoma).

Alternative models explaining the biliary lipid secretion at the canalicular membrane of hepatocytes exist: successive lipid extraction by preformed bile salt micelles, or budding of membrane fragments with formation of mixed micelles. To test the feasibility of the latter mechanism, we developed a mathematical model that describes the formation of lipid microdomains in the canalicular membrane. Bile salt monomers intercalate into the external hemileaflet of the canalicular membrane, to form a rim to liquid disordered domain patches that then pinch off to form nanometer-scale mixed micelles. Model simulations perfectly recapitulate the measured dependence of bile salt-dependent biliary lipid extraction rates upon modulation of the membrane cholesterol (lack or overexpression of the cholesterol transporter Abcg5-Abcg8) and phosphatidylcholine (lack of Mdr2, also known as Abcb4) content. The model reveals a strong dependence of the biliary secretion rate on the protein density of the membrane. Taken together, the proposed model is consistent with crucial experimental findings in the field and provides a consistent explanation of the central molecular processes involved in bile formation., Competing Interests: Competing interestsThe authors declare no competing or financial interests., (© 2018. Published by The Company of Biologists Ltd.)

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

Renal oncocytomas are rare benign tumors of the kidney and characterized by a deficient complex I (CI) enzyme activity of the oxidative phosphorylation (OXPHOS) system caused by mitochondrial DNA (mtDNA) mutations. Yet, little is known about the underlying molecular mechanisms and alterations of metabolic pathways in this tumor. We compared renal oncocytomas with adjacent matched normal kidney tissues on a global scale by multi-omics approaches, including whole exome sequencing (WES), proteomics, metabolomics, and metabolic pathway simulation. The abundance of proteins localized to mitochondria increased more than 2-fold, the only exception was a strong decrease in the abundance for CI subunits that revealed several pathogenic heteroplasmic mtDNA mutations by WES. We also observed renal oncocytomas to dysregulate main metabolic pathways, shunting away from gluconeogenesis and lipid metabolism. Nevertheless, the abundance of energy carrier molecules such as NAD + , NADH, NADP, ATP, and ADP were significantly higher in renal oncocytomas. Finally, a substantial 5000-fold increase of the reactive oxygen species scavenger glutathione can be regarded as a new hallmark of renal oncocytoma. Our findings demonstrate that renal oncocytomas undergo a metabolic switch to eliminate ATP consuming processes to ensure a sufficient energy supply for the tumor., Competing Interests: CONFLICTS OF INTEREST No competing interests are declared.

The liver responds to elevated plasma concentrations of free fatty acids (FFAs) with an enhanced uptake of FFAs and their esterification to triacylglycerol (TAG). On the long term, this may result in massive hepatic TAG accumulation called steatosis hepatitis. In hepatocytes, the poor water-soluble TAG is packed in specialized organelles: Lipid droplets (LDs) serving as transient cellular deposit and lipoproteins (LPs) transporting TAG and cholesterol esters to extra-hepatic tissues. The dynamics of these organelles is controlled by a variety of regulatory surface proteins (RSPs). Assembly and export of VLDLs are mainly regulated by the microsomal transfer protein (MTP) and apoprotein B100. Formation and lipolysis of LDs are regulated by several RSPs. The best studied regulators belong to the PAT (Perilipin/Adipophilin/TIP47) and CIDE families. Knockdown or overexpression of SRPs may significantly affect the total number and size distribution of LDs. Intriguingly, a large cell-to-cell heterogeneity with respect to the number and size of LDs has been found in various cell types including hepatocytes. These findings suggest that the extent of cellular lipid accumulation is determined not only by the imbalance between lipid supply and utilization but also by variations in the expression of RSPs and metabolic enzymes. To better understand the relative regulatory impact of individual processes involved in the cellular TAG turnover, we developed a comprehensive kinetic model encompassing the pathways of the fatty acid and triglyceride metabolism and the main molecular processes governing the dynamics of LDs. The model was parametrized such that a large number of experimental in vitro and in vivo findings are correctly recapitulated. A control analysis of the model revealed that variations in the activity of FFA uptake, diacylglycerol acyltransferase (DGAT) 2, and adipose triglyceride lipase (ATGL) have the strongest influence on the cellular TAG level. We used the model to simulate LD size distributions in human hepatoma cells and hepatocytes exposed to a challenge with FFAs. A random fold change by a factor of about two in the activity of RSPs was sufficient to reproduce the large diversity of droplet size distributions observed in individual cells. Under the premise that the same extent of variability of RSPs holds for the intact organ, our model predicts variations in the TAG content of individual hepatocytes by a factor of about 3-6 depending on the nutritional regime. Taken together, our modeling approach integrates numerous experimental findings on individual processes in the cellular TAG metabolism and LD dynamics metabolism to a consistent state-of-the-art dynamic network model that can be used to study how changes in the external conditions or systemic parameters will affect the TAG content of hepatocytes., (© 2017 Federation of European Biochemical Societies.)

The development of drugs that can inactivate disease-causing cells (e.g. cancer cells or parasites) without causing collateral damage to healthy or to host cells is complicated by the fact that many proteins are very similar between organisms. Nevertheless, due to subtle, quantitative differences between the biochemical reaction networks of target cell and host, a drug can limit the flux of the same essential process in one organism more than in another. We identified precise criteria for this 'network-based' drug selectivity, which can serve as an alternative or additive to structural differences. We combined computational and experimental approaches to compare energy metabolism in the causative agent of sleeping sickness, Trypanosoma brucei, with that of human erythrocytes, and identified glucose transport and glyceraldehyde-3-phosphate dehydrogenase as the most selective antiparasitic targets. Computational predictions were validated experimentally in a novel parasite-erythrocytes co-culture system. Glucose-transport inhibitors killed trypanosomes without killing erythrocytes, neurons or liver cells.

Background: Adaptation of the cellular metabolism to varying external conditions is brought about by regulated changes in the activity of enzymes and transporters. Hormone-dependent reversible enzyme phosphorylation and concentration changes of reactants and allosteric effectors are the major types of rapid kinetic enzyme regulation, whereas on longer time scales changes in protein abundance may also become operative. Here, we used a comprehensive mathematical model of the hepatic glucose metabolism of rat hepatocytes to decipher the relative importance of different regulatory modes and their mutual interdependencies in the hepatic control of plasma glucose homeostasis., Results: Model simulations reveal significant differences in the capability of liver metabolism to counteract variations of plasma glucose in different physiological settings (starvation, ad libitum nutrient supply, diabetes). Changes in enzyme abundances adjust the metabolic output to the anticipated physiological demand but may turn into a regulatory disadvantage if sudden unexpected changes of the external conditions occur. Allosteric and hormonal control of enzyme activities allow the liver to assume a broad range of metabolic states and may even fully reverse flux changes resulting from changes of enzyme abundances alone. Metabolic control analysis reveals that control of the hepatic glucose metabolism is mainly exerted by enzymes alone, which are differently controlled by alterations in enzyme abundance, reversible phosphorylation, and allosteric effects., Conclusion: In hepatic glucose metabolism, regulation of enzyme activities by changes of reactants, allosteric effects, and reversible phosphorylation is equally important as changes in protein abundance of key regulatory enzymes.

Cellular metabolism basically consists of the conversion of chemical compounds taken up from the extracellular environment into energy (conserved in energy-rich bonds of organic phosphates) and a wide array of organic molecules serving as catalysts (enzymes), information carriers (nucleic acids), and building blocks for cellular structures such as membranes or ribosomes. Metabolic modeling aims at the construction of mathematical representations of the cellular metabolism that can be used to calculate the concentration of cellular molecules and the rates of their mutual chemical interconversion in response to varying external conditions as, for example, hormonal stimuli or supply of essential nutrients. Based on such calculations, it is possible to quantify complex cellular functions as cellular growth, detoxification of drugs and xenobiotic compounds or synthesis of exported molecules. Depending on the specific questions to metabolism addressed, the methodological expertise of the researcher, and available experimental information, different conceptual frameworks have been established, allowing the usage of computational methods to condense experimental information from various layers of organization into (self-) consistent models. Here, we briefly outline the main conceptual frameworks that are currently exploited in metabolism research.

Background: Disrupted bile secretion leads to liver damage characterized by inflammation, fibrosis, eventually cirrhosis, and hepatocellular cancer. As obstructive cholestasis often progresses insidiously, markers for the diagnosis and staging of the disease are urgently needed. To this end, we compiled a comprehensive data set of serum markers, histological parameters and transcript profiles at 8 time points of disease progression after bile duct ligation (BDL) in mice, aiming at identifying a set of parameters that could be used as robust biomarkers for transition of different disease progression phases., Results: Statistical analysis of the more than 6,000 data points revealed distinct temporal phases of disease. Time course correlation analysis of biochemical, histochemical and mRNA transcript parameters (=factors) defined 6 clusters for different phases of disease progression. The number of CTGF-positive cells provided the most reliable overall measure for disease progression at histological level, bilirubin at biochemical level, and metalloproteinase inhibitor 1 (Timp1) at transcript level. Prominent molecular events exhibited by strong transcript peaks are found for the transcriptional regulator Nr0b2 (Shp) and 1,25-dihydroxyvitamin D(3) 24-hydroxylase (Cyp24a1) at 6 h. Based on these clusters, we constructed a decision tree of factor combinations potentially useful as markers for different time intervals of disease progression. Best prediction for onset of disease is achieved by fibronectin (Fn1), for early disease phase by Cytochrome P450 1A2 (Cyp1a2), passage to perpetuation phase by collagen1α-1 (Col1a1), and transition to the progression phase by interleukin 17-a (Il17a), with early and late progression separated by Col1a1. Notably, these predictions remained stable even for randomly chosen small sub-sets of factors selected from the clusters., Conclusion: Our detailed time-resolved explorative study of liver homogenates following BDL revealed a well-coordinated response, resulting in disease phase dependent parameter modulations at morphological, biochemical, metabolic and gene expression levels. Interestingly, a small set of selected parameters can be used as diagnostic markers to predict disease stages in mice with cholestatic liver disease.

Proteasomal protein degradation is a key determinant of protein half-life and hence of cellular processes ranging from basic metabolism to a host of immunological processes. Despite its importance the mechanisms regulating proteasome activity are only incompletely understood. Here we use an iterative and tightly integrated experimental and modelling approach to develop, explore and validate mechanistic models of proteasomal peptide-hydrolysis dynamics. The 20S proteasome is a dynamic enzyme and its activity varies over time because of interactions between substrates and products and the proteolytic and regulatory sites; the locations of these sites and the interactions between them are predicted by the model, and experimentally supported. The analysis suggests that the rate-limiting step of hydrolysis is the transport of the substrates into the proteasome. The transport efficiency varies between human standard- and immuno-proteasomes thereby impinging upon total degradation rate and substrate cleavage-site usage.

Imaging of the cellular fluorescence of the reduced form of nicotinamide adenine dinucleotide (phosphate) (NAD(P)H) is one of the few metabolic readouts that enable noninvasive and time-resolved monitoring of the functional status of mitochondria in neuronal tissues. Stimulation-induced transient changes in NAD(P)H fluorescence intensity frequently display a biphasic characteristic that is influenced by various molecular processes, e.g., intracellular calcium dynamics, tricarboxylic acid cycle activity, the malate-aspartate shuttle, the glycerol-3-phosphate shuttle, oxygen supply or adenosine triphosphate (ATP) demand. To evaluate the relative impact of these processes, we developed and validated a detailed physiologic mathematical model of the energy metabolism of neuronal cells and used the model to simulate metabolic changes of single cells and tissue slices under different settings of stimulus-induced activity and varying nutritional supply of glucose, pyruvate or lactate. Notably, all experimentally determined NAD(P)H responses could be reproduced with one and the same generic cellular model. Our computations reveal that (1) cells with quite different metabolic status may generate almost identical NAD(P)H responses and (2) cells of the same type may quite differently contribute to aggregate NAD(P)H responses recorded in brain slices, depending on the spatial location within the tissue. Our computational approach reconciles different and sometimes even controversial experimental findings and improves our mechanistic understanding of the metabolic changes underlying live-cell NAD(P)H fluorescence transients.

Analyzing multiple single nucleotide polymorphisms (SNPs) is a promising approach to finding genetic effects beyond single-locus associations. We proposed the use of multilocus stepwise regression (MSR) to screen for allele combinations as a method to model joint effects, and compared the results with the often used genetic risk score (GRS), conventional stepwise selection, and the shrinkage method LASSO. In contrast to MSR, the GRS, conventional stepwise selection, and LASSO model each genotype by the risk allele doses. We reanalyzed 20 unlinked SNPs related to type 2 diabetes (T2D) in the EPIC-Potsdam case-cohort study (760 cases, 2193 noncases). No SNP-SNP interactions and no nonlinear effects were found. Two SNP combinations selected by MSR (Nagelkerke's R² = 0.050 and 0.048) included eight SNPs with mean allele combination frequency of 2%. GRS and stepwise selection selected nearly the same SNP combinations consisting of 12 and 13 SNPs (Nagelkerke's R² ranged from 0.020 to 0.029). LASSO showed similar results. The MSR method showed the best model fit measured by Nagelkerke's R² suggesting that further improvement may render this method a useful tool in genetic research. However, our comparison suggests that the GRS is a simple way to model genetic effects since it does not consider linkage, SNP-SNP interactions, and no non-linear effects., (© 2015 John Wiley & Sons Ltd/University College London.)

This study focuses on a two-scale, continuum multicomponent model for the description of blood perfusion and cell metabolism in the liver. The model accounts for a spatial and time depending hydro-diffusion-advection-reaction description. We consider a solid-phase (tissue) containing glycogen and a fluid-phase (blood) containing glucose as well as lactate. The five-component model is enhanced by a two-scale approach including a macroscale (sinusoidal level) and a microscale (cell level). The perfusion on the macroscale within the lobules is described by a homogenized multiphasic approach based on the theory of porous media (mixture theory combined with the concept of volume fraction). On macro level, we recall the basic mixture model, the governing equations as well as the constitutive framework including the solid (tissue) stress, blood pressure and solutes chemical potential. In view of the transport phenomena, we discuss the blood flow including transverse isotropic permeability, as well as the transport of solute concentrations including diffusion and advection. The continuum multicomponent model on the macroscale finally leads to a coupled system of partial differential equations (PDE). In contrast, the hepatic metabolism on the microscale (cell level) was modeled via a coupled system of ordinary differential equations (ODE). Again, we recall the constitutive relations for cell metabolism level. A finite element implementation of this framework is used to provide an illustrative example, describing the spatial and time-depending perfusion-metabolism processes in liver lobules that integrates perfusion and metabolism of the liver.

The liver maintains glucose and lipid homeostasis by adapting its metabolic activity to the energy needs of the organism. Communication between hepatocytes and extracellular environment via endocytosis is key to such homeostasis. Here, we addressed the question of whether endosomes are required for gluconeogenic gene expression. We took advantage of the loss of endosomes in the mouse liver upon Rab5 silencing. Strikingly, we found hepatomegaly and severe metabolic defects such as hypoglycemia, hypercholesterolemia, hyperlipidemia, and glycogen accumulation that phenocopied those found in von Gierke's disease, a glucose-6-phosphatase (G6Pase) deficiency. G6Pase deficiency alone can account for the reduction in hepatic glucose output and glycogen accumulation as determined by mathematical modeling. Interestingly, we uncovered functional alterations in the transcription factors, which regulate G6Pase expression. Our data highlight a requirement of Rab5 and the endosomal system for the regulation of gluconeogenic gene expression that has important implications for metabolic diseases., (Copyright © 2015 The Authors. Published by Elsevier Inc. All rights reserved.)

Temporal changes of gene expression are a well-known regulatory feature of all cells, which is commonly perceived as a strategy to adapt the proteome to varying external conditions. However, temporal (rhythmic and non-rhythmic) changes of gene expression are also observed under virtually constant external conditions. Here we hypothesize that such changes are a means to render the synthesis of the metabolic output more efficient than under conditions of constant gene activities. In order to substantiate this hypothesis, we used a flux-balance model of the cellular metabolism. The total time span spent on the production of a given set of target metabolites was split into a series of shorter time intervals (metabolic phases) during which only selected groups of metabolic genes are active. The related flux distributions were calculated under the constraint that genes can be either active or inactive whereby the amount of protein related to an active gene is only controlled by the number of active genes: the lower the number of active genes the more protein can be allocated to the enzymes carrying non-zero fluxes. This concept of a predominantly protein-limited efficiency of gene expression clearly differs from other concepts resting on the assumption of an optimal gene regulation capable of allocating to all enzymes and transporters just that fraction of protein necessary to prevent rate limitation. Applying this concept to a simplified metabolic network of the central carbon metabolism with glucose or lactate as alternative substrates, we demonstrate that switching between optimally chosen stationary flux modes comprising different sets of active genes allows producing a demanded amount of target metabolites in a significantly shorter time than by a single optimal flux mode at fixed gene activities. Our model-based findings suggest that temporal expression of metabolic genes can be advantageous even under conditions of constant external substrate supply.

The bile fluid contains various lipids that are secreted at the canalicular membrane of hepatocytes. As the secretion mechanism is still a matter of debate and a direct experimental observation of the secretion process is not possible so far, we used a mathematical model to simulate the extraction of the major bile lipids cholesterol, phosphatidylcholine and sphingomyelin from the outer leaflet of the canalicular membrane. Lipid diffusion was modeled as random movement on a triangular lattice governed by next-neighbor interaction energies. Phase separation in liquid-ordered and liquid-disordered domains was modeled by assigning two alternative ordering states to each lipid species and minimization of next-neighbor ordering energies. Parameterization of the model was performed such that experimentally determined diffusion rates and phases in ternary lipid mixtures of model membranes were correctly recapitulated. The model describes the spontaneous formation of nanodomains in the external leaflet of the canalicular membrane in a time window between 0.1 ms to 10 ms at varying lipid proportions. The extraction of lipid patches from the bile salt soluble nanodomain into the bile reproduced observed biliary phospholipid compositions for a physiological membrane composition. Comparing the outcome of model simulations with available experimental observations clearly favors the extraction of tiny membrane patches composed of about 100-400 lipids as the likely mechanism of biliary lipid secretion.

Proteasomes are key proteases involved in a variety of processes ranging from the clearance of damaged proteins to the presentation of antigens to CD8+ T-lymphocytes. Which cleavage sites are used within the target proteins and how fast these proteins are degraded have a profound impact on immune system function and many cellular metabolic processes. The regulation of proteasome activity involves different mechanisms, such as the substitution of the catalytic subunits, the binding of regulatory complexes to proteasome gates and the proteasome conformational modifications triggered by the target protein itself. Mathematical models are invaluable in the analysis; and potentially allow us to predict the complex interactions of proteasome regulatory mechanisms and the final outcomes of the protein degradation rate and MHC class I epitope generation. The pioneering attempts that have been made to mathematically model proteasome activity, cleavage preference variation and their modification by one of the regulatory mechanisms are reviewed here.

Background: Accounts of evidence are vital to evaluate and reproduce scientific findings and integrate data on an informed basis. Currently, such accounts are often inadequate, unstandardized and inaccessible for computational knowledge engineering even though computational technologies, among them those of the semantic web, are ever more employed to represent, disseminate and integrate biomedical data and knowledge., Results: We present SEE (Semantic EvidencE), an RDF/OWL based approach for detailed representation of evidence in terms of the argumentative structure of the supporting background for claims even in complex settings. We derive design principles and identify minimal components for the representation of evidence. We specify the Reasoning and Discourse Ontology (RDO), an OWL representation of the model of scientific claims, their subjects, their provenance and their argumentative relations underlying the SEE approach. We demonstrate the application of SEE and illustrate its design patterns in a case study by providing an expressive account of the evidence for certain claims regarding the isolation of the enzyme glutamine synthetase., Conclusions: SEE is suited to provide coherent and computationally accessible representations of evidence-related information such as the materials, methods, assumptions, reasoning and information sources used to establish a scientific finding by adopting a consistently claim-based perspective on scientific results and their evidence. SEE allows for extensible evidence representations, in which the level of detail can be adjusted and which can be extended as needed. It supports representation of arbitrary many consecutive layers of interpretation and attribution and different evaluations of the same data. SEE and its underlying model could be a valuable component in a variety of use cases that require careful representation or examination of evidence for data presented on the semantic web or in other formats.

It is estimated that maintenance of the resting potential of neurons consumes between 15% (in gray matter) and 44% (in fully myelinated white matter) of the brain's total energy budget [1]. This poses the intriguing question why evolution has not strived to lower the permeability of passive ion channels to cut the high resting-state energy budget of the brain. Based on a conceptual mathematical model of neuronal ion currents and action potential (AP) firing we demonstrate that a neuron endowed with small leak currents and correspondingly low energy consumption by the Na(+)/K(+)-ATPase in the resting state may indeed recapitulate all features of normal AP firing. However, the activation and inactivation of such a "low-energy-cost neuron" turns out to be extremely sensitive to small fluctuation of Na(+) currents associated with Na(+)-dependent secondary-active transport that is indispensable for the metabolic integrity of the cell and neurotransmitter recycling. We provide evidence that sufficiently large leak currents function as important stabilizers of the membrane potential and thus are required to allow robust AP firing. Our simulations suggest that the energy demand of the Na(+)/K(+)-ATPase needed to counterbalance passive leak currents cannot be significantly dropped below observed values.

Steadily growing experimental evidence suggests that mitochondrial dysfunction plays a key role in the age-dependent impairment of nerve cells underlying several neurodegenerative diseases. In particular, the citric acid cycle enzyme complex α-ketoglutarate dehydrogenase (KGDHC) and respiratory chain complex I of the respiratory chain often show reduced activities in the dopaminergic neurons involved in Parkinson's disease, both giving rise to an impaired mitochondrial energy metabolism as demonstrated in a number of in vitro studies with cell lines as well as isolated mitochondria. To understand the metabolic regulation underlying these experimental findings we used a detailed kinetic model of mitochondrial energy metabolism. First, we investigated the effect of complex I inhibition on energy production and formation of reactive oxygen species (ROS). Next, we applied the model to a situation where both KGDHC and complex I exhibit reduced activities. These calculations reveal synergistic effects with respect to the energy metabolism but antagonistic effects with respect to ROS formation: the drop in the ATP production capacity is more pronounced than at inhibition of either enzyme complex alone. Interestingly, however, the reduction state of the ROS-generating sites of the impaired complex I becomes significantly lowered if additionally the activity of the KGDHC is reduced. We discuss the pathophysiological consequences of these intriguing findings., (© 2013 FEBS.)

Proteasomes are cellular proteases involved in the degradation of numerous cellular proteins. The 20S proteasome is a cylindrical 28-mer protein complex composed of two outer heptameric α-rings forming the entrance for the protein substrate and two inner heptameric β-rings carrying the catalytic sites. Numerous in vitro studies have provided evidence that the 20S proteasome may degrade peptides of various lengths and even unfolded full-length polypeptide chains. However, a direct demonstration that the 20S proteasome may also cleave surface-attached immobilized peptides is lacking so far. To this end, we used a model system by coupling peptides from different source proteins covalently to the surface of glass beads and applied nanoLC/MS analysis to monitor the generation of proteolytic fragments in the presence of the 20S proteasome. Detectable amounts of cleavage products occurred within a few minutes indicating a much higher cleavage rate than observed with the same substrates in solution. Our finding lends support to the idea that proteasomes may directly degrade segments of membrane-bound proteins protruding into the aqueous phase., (Copyright © 2013 European Peptide Society and John Wiley & Sons, Ltd.)

Characteristics of many tumor types are the reprogramming of metabolism and the occurrence of regional hypoxia. In this work, we investigated the hypothesis that metabolic reprogramming in combination with metabolic zonation of cellular energy metabolism are important factors in promotion of the growth capacity of solid tumors. A tissue-scale model of the two main ATP delivering pathways, glycolysis (GLY) and oxidative phosphorylation (OXP), was used to simulate the energy metabolism within solid tumors under various metabolic strategies. Remarkably, despite the high diversity in the usage of glucose, lactate and oxygen in various spatial regions, the tumor as a whole clearly displays the hallmark of the so-called Warburg effect, i.e. a high rate of glucose consumption and lactate production in the presence of sufficiently high levels of oxygen. Our simulations suggest that an increase in GLY capacity and concomitant decrease in OXP capacity from the periphery towards the center of the tumor improves the availability of oxygen to pericentral tumor cells. The found relationship between the regional oxygen level and the relative share of GLY and OXP capacities supports the view that metabolite availability functions as key regulator of tumor energy metabolism., (Copyright © 2013 WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim.)