Nikolaus Berndt

Oxygen consumption rates during three different neuronal activity states in the hippocampal CA3 network

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 μm 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.

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.

Reduced activity of brain α-ketoglutarate dehydrogenase complex (KGDHC) occurs in a number of neurodegenerative diseases like Parkinsons disease and Alzheimers disease. In order to quantify the relation between diminished KGDHC activity and the mitochondrial ATP generation, redox state, transmembrane potential, and generation of reactive oxygen species (ROS) by the respiratory chain (RC), we developed a detailed kinetic model. 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. By contrast, progressive inhibition of the enzyme aconitase had only little impact on these mitochondrial parameters. As KGDHC is susceptible to ROS-dependent inactivation, we also investigated the reduction state of those sites of the RC proposed to be involved in ROS production. The reduction state of all sites except one decreased with increasing degree of KGDHC inhibition suggesting an ROS-reducing effect of KGDHC inhibition. Our model underpins the important role of reduced KGDHC activity in the energetic breakdown of neuronal cells during development of neurodegenerative diseases.

A Conceptual Mathematical Model of the Dynamic Self-Organisation of Distinct Cellular Organelles

Formation, degradation and renewal of cellular organelles is a dynamic process based on permanent budding, fusion and inter-organelle traffic of vesicles. These processes include many regulatory proteins such as SNAREs, Rabs and coats. Given this complex machinery, a controversially debated issue is the definition of a minimal set of generic mechanisms necessary to enable the self-organization of organelles differing in number, size and chemical composition. We present a conceptual mathematical model of dynamic organelle formation based on interacting vesicles which carry different types of fusogenic proteins (FP) playing the role of characteristic marker proteins. Our simulations (ODEs) show that a de novo formation of non-identical organelles, each accumulating a different type of FP, requires a certain degree of disproportionation of FPs during budding. More importantly however, the fusion kinetics must indispensably exhibit positive cooperativity among these FPs, particularly for the formation of larger organelles. We compared different types of cooperativity: sequential alignment of corresponding FPs on opposite vesicle/organelles during fusion and pre-formation of FP-aggregates (equivalent, e.g., to SNARE clusters) prior to fusion described by Hill kinetics. This showed that the average organelle size in the system is much more sensitive to the disproportionation strength of FPs during budding if the vesicular transport system gets along with a fusion mechanism based on sequential alignments of FPs. Therefore, pre-formation of FP aggregates within the membranes prior to fusion introduce robustness with respect to organelle size. Our findings provide a plausible explanation for the evolution of a relatively large number of molecules to confer specificity on the fusion machinery compared to the relatively small number involved in the budding process. Moreover, we could speculate that a specific cooperativity which may be described by Hill kinetics (aggregates or Rab/SNARE complex formation) is suitable if maturation/identity switching of organelles play a role (bistability).

Physiology-Based Kinetic Modeling of Neuronal Energy Metabolism Unravels the Molecular Basis of NAD(P)H Fluorescence 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 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.

Regulation of Liver Metabolism by the Endosomal GTPase Rab5

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

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

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

The relative importance of kinetic mechanisms and variable enzyme abundances for the regulation of hepatic glucose metabolism – insights from mathematical modeling

BackgroundAdaptation 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.ResultsModel 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.ConclusionIn 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.

Metabolic gradients as key regulators in zonation of tumor energy metabolism: A tissue-scale model-based study

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.

The influence of the chloride currents on action potential firing and volume regulation of excitable cells studied by a kinetic model.

In excitable cells, the generation of an action potential (AP) is associated with transient changes of the intra- and extracellular concentrations of small ions such as Na(+), K(+) and Cl(-). If these changes cannot be fully reversed between successive APs cumulative changes of trans-membrane ion gradients will occur, impinging on the cell volume and the duration, amplitude and frequency of APs. Previous computational studies focused on effects associated with excitation-induced changes of potassium and sodium. Here we present a model based study on the influence of chloride on the fidelity of AP firing and cellular volume regulation during excitation. Our simulations show that depending on the magnitude of the basal chloride permeability two complementary types of responsiveness and volume variability exist: (i) At high chloride permeability (typical for muscle cells), large excitatory stimuli are required to elicit APs; repetitive stimuli of equal strength result in almost identical spike train patterns (Markovian behavior), however, long excitation may lead to after discharges due to an outward directed current of intracellular chloride ions which accumulate during excitation; cell volume changes are large. (ii) At low chloride permeability (e.g., neurons), small excitatory stimuli are sufficient to elicit APs, repetitive stimuli of equal strength produce spike trains with progressively changing amplitude, frequency and duration (short-term memory effects or non-Markovian behavior); cell volume changes are small. We hypothesize that variation of the basal chloride permeability could be an important mechanism of neuronal cells to adapt their responsiveness to external stimuli during learning and memory processes.

Local oxygen homeostasis during various neuronal network activity states in the mouse hippocampus

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