Projects

You are here:

Modeling of Cellular Processes in the Liver

Physiology-based large-scale kinetic model of liver metabolism

Reaction scheme of the metabolic sub-model.

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).

Publication: Berndt N, Bulik S, Wallach I, Wünsch T, König M, Stockmann M, Meierhofer D, Holzhütter HG. HEPATOKIN1 is a biochemistry-based model of liver metabolism for applications in medicine and pharmacology. Nat Commun. 2018 Jun 19;9(1):2386.

Project funding: Systems Biology Programs "Virtual Liver" (grant no. 0315741) and "LiSyM" (grant no. 31L0057 and 031L0058), as well as the e:Bio (Module I) project "HepatomaSys" (grant no.0316172A), all sponsored by the German Federal Ministry of Education and Research (BMBF).

Cooperation partners:

 

 

Multi-scale modeling of liver tissue

Schematic model representation. (A) Model of carbohydrate metabolism describing glycolysis, gluconeogenesis as well as glycogen synthesis and utilization. (B) Sinusoidal unit describing blood flow, nutrient and hormone distribution within the sinusoids.

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 (‘metabolic zonation’) have been reported. Understanding how this variability may affect the regional metabolic capacity of the liver is important for interpretation of functional liver tests and planning of pharmacological and surgical interventions. Here, we treat the liver as an ensemble of a large number (more than a million) of sinusoidal tissue units (STUs), each composed of a single sinusoid surrounded by the space of Disse and a monolayer of hepatocytes. We develop spatio-temporal kinetic models of the STU and calculate the total metabolic output of the liver (arterio-venous glucose difference) by integration across the metabolic output of a sufficiently large number of representative STUs differing in their anatomical structure (thickness and length of the sinusoid, number and size of hepatocytes etc.). Application of the model to the hepatic glucose metabolism provided the following insights: (i) At portal glucose concentrations between 6 to 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 showing the metabolic activities in various liver regions in a time-resolved manner. The model suggests that variations of tissue parameters are equally important as variations of enzyme activities for the control of the arterio-venous glucose difference.

Publications:

Project funding: Systems Biology Programs "Virtual Liver" (grant no. 0315741) and "LiSyM" (grant no. 31L0057 and 031L0058), as well as the e:Bio (Module I) project "HepatomaSys" (grant no.0316172A), all sponsored by the German Federal Ministry of Education and Research (BMBF)

Cooperation partners:

Multilayer control of cellular metabolism: hierarchical or democratic?

Schematic representation of the model of rat hepatocyte carbohydrate metabolism.

Adaptation of 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. 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.

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 – depending on the (patho)physiological condition – control of the hepatic glucose metabolism is mainly exerted by specific enzymes, which are differently controlled by alterations in enzyme abundance, reversible phosphorylation, and allosteric effects.

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.

Publication: Bulik S, Holzhütter HG, Berndt N. The relative importance of kinetic mechanisms and variable enzyme abundances for the regulation of hepatic glucose metabolism - insights from mathematical modeling. BMC Biology, 2016. 14:15

Project funding: Systems Biology Programs "Virtual Liver" (grant no. 0315741) and "LiSyM" (grant no. 31L0057), as well as the e:Bio (Module I) project "HepatomaSys" (grant no.0316172A), all sponsored by the German Federal Ministry of Education and Research (BMBF).

Integration of metabolism and signaling

The regulation of key reaction steps in mutually opposing pathways (e.g., glycolysis and gluconeogenesis, lipid synthesis and lipolysis) by hormone-dependent reversible enzyme phosphorylation represents an important regulatory principle to control the direction of the net flux [1]. The signaling part of the HEPATOKIN1 model [2] comprises the insulin and glucagon dependent regulation of key regulatory enzymes by reversible phosphorylation. The rate laws for These enzymes take into account that the phosphorylated and de-phosphorylated states of the enzyme possess differing maximal activities and kinetic properties. So far, we have used phenomenological mathematical functions to relate the enzyme’s phosphorylation state to the plasma concentrations of glucose [2,3].
To take better into account the mutual influence of the insulin, glucagon and epinephrine signaling pathways under normal and pathophysiological conditions such as diabetes type 2, we plan in a future project to set up kinetic models, which describe the dynamic state of individual constituents (receptors, kinases, phosphatases) by ordinary differential equations. These models will include different cellular compartments (cell membrane, cytosol, mitochondria, and endoplasmic reticulum). The aim of this integrated metabolic-signaling model is to predict the metabolic effects elicited by agonists and antagonists of the insulin, glucagon and epinephrine receptors.

Publications:

  1. König M., Bulik S. and Holzhütter HG. Quantifying the Contribution of the Liver to the Homeostasis of Plasma Glucose: A Detailed Kinetic Model of Hepatic Glucose Metabolism. PLoS Comput Biol. 2012 Jun;8(6):e1002577. Epub 2012 Jun 21. 
  2. Berndt N, Bulik S, Wallach I, Wünsch T, König M, Stockmann M, Meierhofer D, Holzhütter HG. HEPATOKIN1 is a biochemistry-based model of liver metabolism for applications in medicine and pharmacology. Nat Commun, 2018. 9(1): p. 2386.
  3. Bulik S, Holzhütter HG, Berndt N. The relative importance of kinetic mechanisms and variable enzyme abundances for the regulation of hepatic glucose metabolism - insights from mathematical modeling. BMC Biology, 2016. 14:15

Project funding:
Systems Biology Programs "Virtual Liver" (grant no. 0315741) and "LiSyM" (grant no. 31L0057 and 031L0058), as well as the e:Bio (Module I) project "HepatomaSys" (grant no.0316172A), all sponsored by the German Federal Ministry of Education and Research (BMBF).

Hepatic lipid droplet metabolism

Schematic representation of the processes included into the LD model.

The liver responds to elevated plasma concentrations of free fatty acids (FFAs) with enhanced uptake and esterification of FFAs to triacylglycerol (TAG). This may result in massive hepatic TAG accumulation called fatty liver (steatosis hepatis), the first stage on the route towards more serious liver diseases, such as cirrhosis, fibrosis or hepatocellular carcinoma. In hepatocytes, the poor water-soluble TAG is packed in lipid droplets (LDs) serving as transient cellular deposit or lipoproteins transporting TAG and cholesterol esters to extra-hepatic tissues. The dynamics of these ‘organelles’ is controlled by a variety of regulatory surface proteins (RSPs). Knockdown or overexpression of RSPs 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 TAG metabolism and the main molecular processes governing the dynamics of LDs [1]. We are using the model to investigate LD size distributions in human hepatocytes under physiological and pathological conditions such as steatosis, fibrosis, cirrhosis or hepatocellular carcinoma [2].

Publications:

  1. Wallstab C, Eleftheriadou D, Schulz T, Damm G, Seehofer D, Borlak J, Holzhütter HG, Berndt N. A unifying mathematical model of lipid droplet metabolism reveals key molecular players in the development of hepatic steatosis. FEBS J, 2017. 284(19): p. 3245-3261.
  2. Berndt N, Eckstein J, Heucke N, Gajowski R, Stockmann M, Meierhofer D, Holzhütter HG. Characterization of Lipid and Lipid Droplet Metabolism in Human HCC. Cells 2019, 8(5), 512.

Project funding:Graduate school "Computational Systems Biology" (GRK 1722)sponsored by the DFG (German Research Foundation) and the Systems Biology Program "LiSyM"sponsored by the German Federal Ministry of Education and Research (BMBF) (grant no. 31L0057 and 031L0058) and the Max Planck Society.

Cooperation partners:

Membrane domain formation and lipid secretion into the bile

Liquid ordered (blue, red) and liquid disordered (dark blue, green) membrane domains with raft proteins (white) and non-raft proteins (black) of different sizes (4/3/2 nm radius).
Simulated (shadowed line) and experimental values (dots) of bile salt dependent secretion of cholesterol (left) an phospholipid (right) into the bile for wild type (blue) and Abcb4(+/-) knock out mice.

We developed a mathematical model of lateral diffusion of lipids and proteins in cellular membranes. The movement of lipids and proteins along the membrane surface is modeled as a movement on a triangular lattice, governed by nearest neighbor interactions. The lipids may switch between two alternative states of ordering energy resulting in different mobilities. Minimizing the ordering energies results in the formation of liquid ordered or liquid disordered phase domains. The model also includes proteins of two different species that have a high affinity for either one of the two phases. The lipid and protein mobilities were parameterized using experimental data from different model membranes. The influence of protein size and density on the formation of lipid domains can be studied.

Model simulations provided support for a budding mechanism of lipid transfer into the bile consisting in the bile salt dependent extraction of membrane patches from liquid disordered microdomains of the canalicular membrane. We applied the model to the canalicular membrane of hepatocytes to study how changes of the lipid composition and protein density may influence the size distribution of microdomains and efficiency of lipid extraction into the bile. Our simulations recapitulate the dependence of lipid secretion from the bile salt secretion measured in mouse models.

Publications:

Project funding:SFB 618 "Theoretical Biology: Robustness, Modularity and Evolutionary Design of Living Systems" (grant no. 5485271) and the Research Training Group (Graduiertenkolleg 1772) "Computational Systems Biology", both sponsored by the DFG (German Research Foundation) as well as the Systems Biology Program "LiSyM" (grant no. 31L0057 and 031L0058), sponsored by the German Federal Ministry of Education and Research (BMBF).

Cooperation partner: Frank Lammert (Saarland University Medical Center and Saarland University Faculty of Medicine, Gastroenterology and Endocrinology

Non-invasive liver function diagnostics

Organization of the project in three interconnected subprojects.

The ultimate (overall) objective of this research project is the establishment of a completely novel systems medicine approach to liver diagnostics that combines multi-scale modeling of liver metabolism and tissue architecture with non-invasively generated data from liver function tests (breath tests) and magnetic resonance elastography (MRE) for the construction of a whole-organ spatially resolved 3D-liver function model (3D-LFM) depicting regional changes in central metabolic functions of the liver. This novel diagnostic tool will lend support to the clinician’s assessment of disease severity, prognosis of disease progression and functional characterization of the potential liver remnant prior to a partial hepatectomy. Besides this important medical application, the 3D-LFM will help to decipher the relative contribution of changes in blood perfusion, tissue architecture and activity of metabolic pathways in the regional metabolic performance of the liver.

Publications:

Project funding: Systems Biology Programs "Virtual Liver" (grant no. 0315741) and "LiSyM" (grant no. 31L0057 and 031L0058) sponsored by the German Federal Ministry of Education and Research (BMBF).

Cooperation partners:

Metabolic alterations in liver pathologies

Publications:

Liver metabolism of adolescents with non-alcoholic fatty liver disease (NAFLD)

Publications: 

Modeling of Cellular Processes in Neuronal Cells

Oxygen consumption rates in brain slices

Depth profiles of partial oxygen pressure (pO2) during three different activity states. (A) Representative sample traces of pO2 depth profiles in the absence of spiking (TTX, black trace), spontaneous network activity (SPON, dark gray trace), and cholinergically induced gamma oscillations (GAM, light gray trace). (B) Quantification of lowest pO2 values as determined during the three different activity states. (C) Quantification of pO2 values at five defined depths in slice cultures.

The brain is an organ with high metabolic rate. However, little is known about energy utilization during different activity states of neuronal networks. We addressed this issue in area CA3 of hippocampal slice cultures under well-defined recording conditions using a 20% O2 gas mixture. We combined recordings of local field potential and interstitial partial oxygen pressure (pO2) during three different activity states, namely fast network oscillations in the gamma-frequency band (30 to 100 Hz), spontaneous network activity and absence of spiking (action potentials). Oxygen consumption rates were determined by pO2 depth profiles with high spatial resolution and a mathematical model that considers convective transport, diffusion, and activity-dependent consumption of oxygen. We show that: (1) Relative oxygen consumption rate during cholinergic gamma oscillations was 2.2-fold and 5.3-fold higher compared with spontaneous activity and absence of spiking, respectively. (2) Gamma oscillations were associated with a similar large decrease in pO2 as observed previously with a 95% O2 gas mixture. (3) Sufficient oxygenation during fast network oscillations in vivo is ensured by the calculated critical radius of 30 to 40 mm around a capillary. We conclude that the structural and biophysical features of brain tissue permit variations in local oxygen consumption by a factor of about five [1].

Publications:

  1. Huchzermeyer C*, Berndt N*, Holzhütter HG*, Kann O*. Oxygen consumption rates during three different neuronal activity states in the hippocampal CA3 network.J Cereb Blood Flow Metab. 2013 Feb;33(2):263-71.

Project funding: Collaborative Research Center for "Theoretical Biology: Robustness, Modularity and Evolutionary Design of Living Systems"SFB 618(grant no. 5485271) sponsored by the DFG (German Research Foundation).

Cooperation partners:

How NAD(P)H fluorescence mirrors neuronal energy metabolism

(A) Reactions and transport processes included in the single-cell kinetic model. (B) Schematic representation of the slice model used to simulate spatial oxygen gradients within a brain slice. (C) Schematic representation of the tissue model used to simulate in vivo NADH transients.

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 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 [1]. 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. In a subsequent study, we investigated the energy metabolism underlying cortical information processing [2]. We concluded 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. In summary, our data show that energy expenditure is strongly dependent on the neuronal network activity state and may reach critical levels during higher brain functions.

Publications:

  1. Berndt N, Kann O, Holzhütter HG. Physiology-based kinetic modeling of neuronal energy metabolism unravels the molecular basis of temporal NAD(P)H fluorescence profiles. J Cereb Blood Flow Metab. 2015 Sep;35(9):1494-506.
  2. Schneider J*, Berndt N*, Papageorgiou IE, Maurer J, Bulik S, Both M, Draguhn A, Holzhütter HG, Kann O. Local oxygen homeostasis during various neuronal network activity states in the mouse hippocampus.J Cereb Blood Flow Metab. 2019 May;39(5):859-873.

Project funding: The projects were in part funded by the German Systems Biology Program “Virtual Liver” (grant no. 0315741) sponsored by the German Federal Ministry of Education and Research (BMBF) and by the German Research Foundation (DFG) within the Collaborative Research Center (SFB) 1134.

Cooperation partner: Oliver Kann (Medical Faculty of Heidelberg, Institute of Physiology and Pathophysiology)

Metabolic alterations in neurodegenerative diseases

Schematic of the mathematical model of mitochondrial energy metabolism.

Steadily growing experimental evidence suggests that mitochondrial dysfunction plays a key role in the age-dependent impairment of nerve cells underlying several neurodegenerative diseases. Especially, reduced activity of brain α-ketoglutarate dehydrogenase complex (KGDHC), reduced activity of complex I of the respiratory chain (RC) and increased reactive oxygen species (ROS) production occurs in a number of neurodegenerative diseases like Parkinson's disease and Alzheimer's disease. To understand the metabolic Regulation underlying these experimental findings we developed and applied a detailed kinetic model of mitochondrial energy metabolism. Model simulations revealed a threshold-like decline of the ATP production rate at about 60% inhibition of KGDHC accompanied by a significant increase of the mitochondrial Membrane potential. We also showed that the reduction state of those sites of the respiratory chain proposed to be involved in ROS production decreased with increasing degree of KGDHC inhibition suggesting a ROS-reducing effect of KGDHC inhibition [1].
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 [2].

Publications:

  1. Berndt N, Bulik S, Holzhütter HG. Kinetic Modeling of the Mitochondrial Energy Metabolism of Neuronal Cells: The Impact of Reduced α-Ketoglutarate Dehydrogenase Activities on ATP Production and Generation of Reactive Oxygen Species. Int J Cell Biol. 2012;2012:757594.
  2. Berndt N, Holzhütter HG, Bulik S. Implications of enzyme deficiencies on mitochondrial energy metabolism and reactive oxygen species formation of neurons involved in rotenone-induced Parkinson's disease: a model-based analysis. FEBS J. 2013 Oct;280(20):5080-93.

Project funding: The project was in part funded by the German Systems Biology Program “Virtual Liver” (grant no. 0315741) sponsored by the German Federal Ministry of Education and Research (BMBF).

Impact of anesthetics on cerebral energy metabolism during light and deep anesthesia

Illustration of the effects of propofol on neuronal functionality during and after anesthesia.

General anesthesia is a drug-induced, reversible state of unconsciousness, amnesia, analgesia and akinesia. The cortical electroencephalogram displays typical dose-dependent changes during anesthesia with characteristic stages of neuronal activity. Despite undisputable improvements in anesthesiology, major concerns related to the long-term effects of anesthetics on the central nervous system are rising. Specifically, deep anesthesia has been associated with postoperative delirium, long lasting postoperative cognitive dysfunction and increased mortality. The underlying role of anesthetics in these neurological complications remains unclear and needs urgent clarification.
Propofol is the most frequently used intravenous anesthetic for induction and maintenance of anesthesia acting primarily as a GABAA-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 (CMRO2), 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.
In a current study, we are investigating the effect of the gas anesthetics isoflurane in neuronal transmission and metabolism in anesthetized Wistar rats and in brain slices of the same species using the same methods as above.

Publication: Berndt N, Rösner J, Haq RU, Kann O, Kovács R, Holzhütter HG, Spies C, Liotta A. Possible neurotoxicity of the anesthetic propofol: evidence for the inhibition of complex II of the respiratory chain in area CA3 of rat hippocampal slices. Arch Toxicol. 2018 Oct;92(10):3191-3205.

Project funding: This work is in part funded by the DFG grant no. 650953 and 408355133 as well as the German Systems Biology Program "LiSyM" (grant no. 31L0057) sponsored by the German Federal Ministry of Education and Research (BMBF). Agustin Liotta is participant in the BIH Charité Clinician Scientist Program funded by the CharitéUniversitätsmedizin Berlin and the Berlin Institute of Health.

Cooperation partners:

Influence of vessel structure on neuronal energy metabolism in different species

 

The Neurovascular unit: neurovascular coupling in patients with traumatic brain injury

 

Spreading depolarization

 

Modeling of Cellular Processes in the Heart

Integrative model of cardiac metabolism

Reaction scheme of the metabolic sub-model.

The heart is energetically one of the most expensive organs. One third of the cellular volume of cardiac myocytes is occupied by mitochondria. Per gram tissue, the heart has the highest oxygen consumption rate and the ATP turnover during one day amounts to 20 times its own weight. This requires a robust and high rate of ATP production to maintain cardiac functionality. ATP is spent on electrophysiological processes of ion pumping as well as on mechanical work in its contractile apparatus. Perturbations in ATP-generating processes may therefore directly affect contractile function. The heart can rely on any energy source available like carbohydrates, amino acids, lipids and ketone bodies. Under normal conditions, oxidation of free fatty acids is the prevailing energy source contributing around 70% to ATP production rate, while the utilization of glucose and pyruvate oxidation become increasingly important during ischemia, hypoxia or increased work load. The use of different substrates is tightly regulated under physiological conditions and there is ample crosstalk between the different metabolic pathways.

Kinetic modeling of cardiac metabolism has a long tradition starting in the late 70’s, but all of the available models neglect crucial factors determining the energetic status of the heart, such as influence of alternating substrate supply, hormonal metabolic control or variable gene expression of key metabolic enzymes necessary for the understanding of metabolic alterations in heart disease. In this work, we developed a kinetic multi-pathway model for cardiomyocytes with hitherto unprecedented scope and level of detail. The model includes the regulation of enzyme activities by allosteric effectors, hormone-dependent reversible phosphorylation and variable protein abundances. For each enzyme, rate equations have been developed that take into account the enzyme’s kinetic and regulatory features as obtained in in vitro assays. One possible application of the model is its ability to reproduce the metabolic state of the healthy heart under different nutritional and load conditions.

Cooperation partner: Titus Kühne (Charité, Institute for Imaging Science and Computational Modelling in Cardiovascular Medicine)

Previous publications:

Cardiac metabolism of patients with heart failure

Mean substrate utilization and standard deviations for all controls (green) and patients with aortic stenosis (red) in response to alterations in substrate availability. (A) glucose utilization; (B) lactate utilization; (C) oxygen consumption rate; (D) fatty acid utilization; (E) reciprocal utilization of glucose and fatty acids in response to external glucose utilization.

For proper functionality, the heart relies on coordinated utilization of different energy providing substrates like glucose, fatty acids, glycogen, and lactate. Depending on substrate availability and energy demand, the heart needs to adapt its internal energy delivering pathways to ensure demand matching energy supply. In pathological situations like aortic stenosis (AS), the efficiency of cardiac muscle activity is disturbed and maladaptation might lead to metabolic alterations contributing to declined cardiac function. Experimental assessment of cardiac energy metabolism is not possible due to ethical and technical restrictions.

In this project, we present a detailed, comprehensive, biochemistry-based kinetic model of the central cardiac metabolism including the regulation of enzymes by kinetic allosteric and hormonal regulation. We show the ability of the model to investigate substrate utilization under different conditions. We use the model to investigate the alterations in cardiac energy metabolism in a cohort of patients with AS showing that a subgroup of patients displays metabolic alterations. The figure shows that mean glucose uptake rate (A) as well as mean oxygen consumption rate (C) are decreased in patients with AS. The same holds for lactate utilization (B) at physiological glucose concentrations below 10 mM and fatty acid utilization at low glucose levels below 5 mM (D). These changes amount to a decreased fatty acid vs. glucose uptake relation in AS patients (E). Importantly, for all simulated cases, the exchange rates for the various metabolites have much larger standard deviations for the patients with AS (red), indicating that this group is much less homogenous than the controls (green).

Cooperation partner: Titus Kühne (Charité, Institute for Imaging Science and Computational Modelling in Cardiovascular Medicine)

Cardiac metabolism of diabetic patients

Graphical work plan description of the DFG project “Mathematical modeling of the metabolic implications of the diabetic heart”.

Diabetes mellitus is an epidemically growing disease worldwide having an overall prevalence of 9.8% in Germany in 2015, with the vast majority of cases (9.5%) attributable to type2 diabetes mellitus. Heart failure, the most common cardiovascular disease associated with diabetes, is a clinical syndrome in which myocardial pump function is inadequate for maintaining and supporting an individual’s physiological requirements. Heart failure in a patient with diabetes may arise from myocardial damage resulting from an ischemic, thrombotic event. In many cases, however, heart failure cannot be attributed to any cardiovascular disease, such as hypertension or coronary artery disease.
Adaptive processes start often at the cellular level by changes in signaling and metabolic pathways, typically evolve to changes in the structural organization of the tissue as, for example, enhanced formation of extracellular matrix (fibrosis) and finally result in alterations of functional parameters such as the cardiac output. A major problem in the treatment of cardiovascular diseases consists in the poor predictability of the responses that are potentially elicited by a medical intervention, whether it is dietary, pharmacologically or surgically. In the worst case, treatment-induced adaptive changes can even exacerbate the pathological situation. A promising approach to overcome this dilemma consists in the use of mathematical models, which integrate existing knowledge on central molecular and physiological circuits operative at the cellular levels and provide reliable predictions of the heart functional capacity and performance in response to intervention.
The goal of this project is to systematically investigate the metabolic and functional changes associated with the diabetic heart. To this end, we will develop, test and verify a computational model of cardiac energy metabolism. The main objective is to understand the short-term and long-term metabolic adaptation of the cardiomyocyte and the functional metabolic changes arising from changes in metabolic enzyme abundance and signaling pathways in dependence of external Substrate supply, hormonal stimuli and internal demand.

Project funding: This project is funded by the German Research Foundation (DFG) (grant no. 422215721).

Cooperation partner: Tilman Grune (German Institute of Human Nutrition Potsdam-Rehbrücke (DIfE), Dept. of Molecular Toxicology)

Metabolic alteration of the heart in systemic inflammation during viral infection