A mathematical model of the mitochondrial NADH shuttles and anaplerosis in the pancreatic beta-cell

Kinetic modelling of β-cell metabolism reveals control points in the insulin-regulating pyruvate cycling pathways

Insulin, a key hormone in the regulation of glucose homoeostasis, is secreted by pancreatic β-cells in response to elevated glucose levels. Insulin is released in a biphasic manner in response to glucose metabolism in β-cells. The first phase of insulin secretion is triggered by an increase in the ATP:ADP ratio; the second phase occurs in response to both a rise in ATP:ADP and other key metabolic signals, including a rise in the NADPH:NADP+ ratio. Experimental evidence indicates that pyruvate-cycling pathways play an important role in the elevation of the NADPH:NADP+ ratio in response to glucose. The authors developed a kinetic model for the tricarboxylic acid cycle and pyruvate cycling pathways. The authors successfully validated the model against experimental observations and performed a sensitivity analysis to identify key regulatory interactions in the system. The model predicts that the dicarboxylate carrier and the pyruvate transporter are the most important regulators of pyruvate cycling and NADPH production. In contrast, the analysis showed that variation in the pyruvate carboxylase flux was compensated by a response in the activity of mitochondrial isocitrate dehydrogenase (ICDm ) resulting in minimal effect on overall pyruvate cycling flux. The model predictions suggest starting points for further experimental investigation, as well as potential drug targets for the treatment of type 2 diabetes.

A pathway model of glucose-stimulated insulin secretion in the pancreatic β-cell

The pancreas plays a critical role in maintaining glucose homeostasis through the secretion of hormones from the islets of Langerhans. Glucose-stimulated insulin secretion (GSIS) by the pancreatic β-cell is the main mechanism for reducing elevated plasma glucose. Here we present a systematic modeling workflow for the development of kinetic pathway models using the Systems Biology Markup Language (SBML). Steps include retrieval of information from databases, curation of experimental and clinical data for model calibration and validation, integration of heterogeneous data including absolute and relative measurements, unit normalization, data normalization, and model annotation. An important factor was the reproducibility and exchangeability of the model, which allowed the use of various existing tools. The workflow was applied to construct a novel data-driven kinetic model of GSIS in the pancreatic β-cell based on experimental and clinical data from 39 studies spanning 50 years of pancreatic, islet, and β-cell research in humans, rats, mice, and cell lines. The model consists of detailed glycolysis and phenomenological equations for insulin secretion coupled to cellular energy state, ATP dynamics and (ATP/ADP ratio). Key findings of our work are that in GSIS there is a glucose-dependent increase in almost all intermediates of glycolysis. This increase in glycolytic metabolites is accompanied by an increase in energy metabolites, especially ATP and NADH. One of the few decreasing metabolites is ADP, which, in combination with the increase in ATP, results in a large increase in ATP/ADP ratios in the β-cell with increasing glucose. Insulin secretion is dependent on ATP/ADP, resulting in glucose-stimulated insulin secretion. The observed glucose-dependent increase in glycolytic intermediates and the resulting change in ATP/ADP ratios and insulin secretion is a robust phenomenon observed across data sets, experimental systems and species. Model predictions of the glucose-dependent response of glycolytic intermediates and biphasic insulin secretion are in good agreement with experimental measurements. Our model predicts that factors affecting ATP consumption, ATP formation, hexokinase, phosphofructokinase, and ATP/ADP-dependent insulin secretion have a major effect on GSIS. In conclusion, we have developed and applied a systematic modeling workflow for pathway models that allowed us to gain insight into key mechanisms in GSIS in the pancreatic β-cell.

A Consensus Model of Glucose-Stimulated Insulin Secretion in the Pancreatic β -Cell

The pancreas plays a critical role in maintaining glucose homeostasis through the secretion of hormones from the islets of Langerhans. Glucose-stimulated insulin secretion (GSIS) by the pancreatic β -cell is the main mechanism for reducing elevated plasma glucose. Here we present a systematic modeling workflow for the development of kinetic pathway models using the Systems Biology Markup Language (SBML). Steps include retrieval of information from databases, curation of experimental and clinical data for model calibration and validation, integration of heterogeneous data including absolute and relative measurements, unit normalization, data normalization, and model annotation. An important factor was the reproducibility and exchangeability of the model, which allowed the use of various existing tools. The workflow was applied to construct the first consensus model of GSIS in the pancreatic β -cell based on experimental and clinical data from 39 studies spanning 50 years of pancreatic, islet, and β -cell research in humans, rats, mice, and cell lines. The model consists of detailed glycolysis and equations for insulin secretion coupled to cellular energy state (ATP/ADP ratio). Key findings of our work are that in GSIS there is a glucose-dependent increase in almost all intermediates of glycolysis. This increase in glycolytic metabolites is accompanied by an increase in energy metabolites, especially ATP and NADH. One of the few decreasing metabolites is ADP, which, in combination with the increase in ATP, results in a large increase in ATP/ADP ratios in the β -cell with increasing glucose. Insulin secretion is dependent on ATP/ADP, resulting in glucose-stimulated insulin secretion. The observed glucose-dependent increase in glycolytic intermediates and the resulting change in ATP/ADP ratios and insulin secretion is a robust phenomenon observed across data sets, experimental systems and species. Model predictions of the glucose-dependent response of glycolytic intermediates and insulin secretion are in good agreement with experimental measurements. Our model predicts that factors affecting ATP consumption, ATP formation, hexokinase, phosphofructokinase, and ATP/ADP-dependent insulin secretion have a major effect on GSIS. In conclusion, we have developed and applied a systematic modeling workflow for pathway models that allowed us to gain insight into key mechanisms in GSIS in the pancreatic β -cell.

Pharmacological inhibition of thioredoxin reductase increases insulin secretion and diminishes beta cell viability

Apparently, both a decrease in beta cell function and in beta cell mass contribute to the progressive worsening of type 2 diabetes. So, it is of particular interest to define factors which are relevant for the regulation of insulin secretion and at the same time for the maintenance of beta cell mass. The NADPH-thioredoxin system has a candidate role for such a dual function. Here, we have characterized the effects of a highly specific inhibitor of thioredoxin reductase, AM12, on the viability and function of insulin-secreting MIN6 cells and isolated NMRI mouse islets. Viability was checked by MTT testing and the fluorescent live-dead assay. Apoptosis was assessed by annexin V assay. Insulin secretion of perifused islets was measured by ELISA. The cytosolic Ca²⁺ concentration was measured by the Fura technique. Acute exposure of perifused pancreatic islets to 5 μM AM12 was without significant effect on insulin secretion. Islets cultured for 24 h in 0.5 or 5 μM AM12 showed unchanged basal secretion during perifusion, but the response to 30 mM glucose was significantly enhanced by 5 μM. Twenty-four-hour exposure to 5 μM AM12 proved to be without effect on the viability of MIN6 cells, whereas longer exposure was clearly toxic. Islets were more susceptible, showing initial signs of apoptosis after 24-h exposure to 5 μM AM12. The activity of the NADPH-thioredoxin system is indispensable for beta cell viability but may have a limiting effect on glucose-induced insulin secretion.

Simple Kinetic Models in Molecular Chronobiology

Circadian rhythms are constituted by a complex dynamical system with intertwined feedback loops, molecular switches, and self-sustained oscillations. Mathematical modeling supports understanding available heterogeneous kinetic data, highlights basic mechanisms, and can guide experimental research. Here, we introduce the basic steps from a biological question to simple models providing insight into gene-regulatory mechanisms. We illustrate the general approach by three examples: modeling decay processes, clock-controlled genes, and self-sustained oscillations.

Mathematical model of metabolism and electrophysiology of amino acid and glucose stimulated insulin secretion: in vitro validation using a β-cell line

We integrated biological experimental data with mathematical modelling to gain insights into the role played by L-alanine in amino acid-stimulated insulin secretion (AASIS) and in D-glucose-stimulated insulin secretion (GSIS), details important to the understanding of complex β-cell metabolic coupling relationships. We present an ordinary differential equations (ODEs) based simplified kinetic model of core metabolic processes leading to ATP production (glycolysis, TCA cycle, L-alanine-specific reactions, respiratory chain, ATPase and proton leak) and Ca(2+) handling (essential channels and pumps in the plasma membrane) in pancreatic β-cells and relate these to insulin secretion. Experimental work was performed using a clonal rat insulin-secreting cell line (BRIN-BD11) to measure the consumption or production of a range of important biochemical parameters (D-glucose, L-alanine, ATP, insulin secretion) and Ca(2+) levels. These measurements were then used to validate the theoretical model and fine-tune the parameters. Mathematical modelling was used to predict L-lactate and L-glutamate concentrations following D-glucose and/or L-alanine challenge and Ca(2+) levels upon stimulation with a non metabolizable L-alanine analogue. Experimental data and mathematical model simulations combined suggest that L-alanine produces a potent insulinotropic effect via both a stimulatory impact on β-cell metabolism and as a direct result of the membrane depolarization due to Ca(2+) influx triggered by L-alanine/Na(+) co-transport. Our simulations indicate that both high intracellular ATP and Ca(2+) concentrations are required in order to develop full insulin secretory responses. The model confirmed that K(+) ATP channel independent mechanisms of stimulation of intracellular Ca(2+) levels, via generation of mitochondrial coupling messengers, are essential for promotion of the full and sustained insulin secretion response in β-cells.

Plasma membrane electron transport in pancreatic β-cells is mediated in part by NQO1

Plasma membrane electron transport (PMET), a cytosolic/plasma membrane analog of mitochondrial electron transport, is a ubiquitous system of cytosolic and plasma membrane oxidoreductases that oxidizes cytosolic NADH and NADPH and passes electrons to extracellular targets. While PMET has been shown to play an important role in a variety of cell types, no studies exist to evaluate its function in insulin-secreting cells. Here we demonstrate the presence of robust PMET activity in primary islets and clonal β-cells, as assessed by the reduction of the plasma membrane-impermeable dyes WST-1 and ferricyanide. Because the degree of metabolic function of β-cells (reflected by the level of insulin output) increases in a glucose-dependent manner between 4 and 10 mM glucose, PMET was evaluated under these conditions. PMET activity was present at 4 mM glucose and was further stimulated at 10 mM glucose. PMET activity at 10 mM glucose was inhibited by the application of the flavoprotein inhibitor diphenylene iodonium and various antioxidants. Overexpression of cytosolic NAD(P)H-quinone oxidoreductase (NQO1) increased PMET activity in the presence of 10 mM glucose while inhibition of NQO1 by its inhibitor dicoumarol abolished this activity. Mitochondrial inhibitors rotenone, antimycin A, and potassium cyanide elevated PMET activity. Regardless of glucose levels, PMET activity was greatly enhanced by the application of aminooxyacetate, an inhibitor of the malate-aspartate shuttle. We propose a model for the role of PMET as a regulator of glycolytic flux and an important component of the metabolic machinery in β-cells.

Explaining the enigmatic K(M) for oxygen in cytochrome c oxidase: a kinetic model

We present a mathematical model for the functioning of proton-pumping cytochrome c oxidase, consisting of cyclic conversions between 26 enzyme states. The model is based on the mechanism of oxygen reduction and linked proton translocation postulated by Wikström and Verkhovsky (2007). It enables the calculation of the steady-state turnover rates and enzyme-state populations as functions of the cytochrome c reduction state, oxygen concentration, membrane potential, and pH on either side of the inner mitochondrial membrane. We use the model to explain the enigmatic decrease in oxygen affinity of the enzyme that has been observed in mitochondria when the proton-motive force is increased. The importance of the 26 transitions in the mechanism of cytochrome oxidase for the functional properties of cytochrome oxidase is compared through Metabolic Control Analysis. The control of the K(M) value is distributed mainly between the steps in the mechanism that involve electrogenic proton movements, with both positive and negative contributions. Positive contributions derive from the same steps that control enzyme turnover rate in the model. Limitations and possible further applications of the model are discussed.

Usurping the mitochondrial supremacy: extramitochondrial sources of reactive oxygen intermediates and their role in beta cell metabolism and insulin secretion

Insulin secretion from pancreatic beta cells is a process dependent on metabolism. While oxidative stress is a well-known inducer of beta cell toxicity and impairs insulin secretion, recent studies suggest that low levels of metabolically-derived reactive oxygen intermediates (ROI) also play a positive role in insulin secretion. Glucose metabolism is directly correlated with ROI production, particularly in beta cells in which glucose uptake is proportional to the extracellular concentration of glucose. Low levels of exogenously added ROI or quinones, which stimulate ROI production, positively affect insulin secretion, while antioxidants block insulin secretion, suggesting that ROI activate unidentified redox-sensitive signal transduction components within these cells. The mitochondria are one source of ROI: increased metabolic flux increases mitochondrial membrane potential resulting in electron leakage and adventitious ROI production. A second source of ROI are cytosolic and plasma membrane oxidoreductases which oxidize NAD(P)H and directly produce ROI through the reduction of molecular oxygen. The mechanism of ROI-mediated potentiation of insulin secretion remains an important topic for future study.

Glucose sensing in the pancreatic beta cell: a computational systems analysis

Background

Pancreatic beta-cells respond to rising blood glucose by increasing oxidative metabolism, leading to an increased ATP/ADP ratio in the cytoplasm. This leads to a closure of K_ATP channels, depolarization of the plasma membrane, influx of calcium and the eventual secretion of insulin. Such mechanism suggests that beta-cell metabolism should have a functional regulation specific to secretion, as opposed to coupling to contraction. The goal of this work is to uncover contributions of the cytoplasmic and mitochondrial processes in this secretory coupling mechanism using mathematical modeling in a systems biology approach.

Methods

We describe a mathematical model of beta-cell sensitivity to glucose. The cytoplasmic part of the model includes equations describing glucokinase, glycolysis, pyruvate reduction, NADH and ATP production and consumption. The mitochondrial part begins with production of NADH, which is regulated by pyruvate dehydrogenase. NADH is used in the electron transport chain to establish a proton motive force, driving the F₁F₀ ATPase. Redox shuttles and mitochondrial Ca²⁺ handling were also modeled.

Results

The model correctly predicts changes in the ATP/ADP ratio, Ca²⁺ and other metabolic parameters in response to changes in substrate delivery at steady-state and during cytoplasmic Ca²⁺ oscillations. Our analysis of the model simulations suggests that the mitochondrial membrane potential should be relatively lower in beta cells compared with other cell types to permit precise mitochondrial regulation of the cytoplasmic ATP/ADP ratio. This key difference may follow from a relative reduction in respiratory activity. The model demonstrates how activity of lactate dehydrogenase, uncoupling proteins and the redox shuttles can regulate beta-cell function in concert; that independent oscillations of cytoplasmic Ca²⁺ can lead to slow coupled metabolic oscillations; and that the relatively low production rate of reactive oxygen species in beta-cells under physiological conditions is a consequence of the relatively decreased mitochondrial membrane potential.

Conclusion

This comprehensive model predicts a special role for mitochondrial control mechanisms in insulin secretion and ROS generation in the beta cell. The model can be used for testing and generating control hypotheses and will help to provide a more complete understanding of beta-cell glucose-sensing central to the physiology and pathology of pancreatic β-cells.

Role of NADH/NAD+ transport activity and glycogen store on skeletal muscle energy metabolism during exercise: in silico studies

Skeletal muscle can maintain ATP concentration constant during the transition from rest to exercise, whereas metabolic reaction rates may increase substantially. Among the key regulatory factors of skeletal muscle energy metabolism during exercise, the dynamics of cytosolic and mitochondrial NADH and NAD+ have not been characterized. To quantify these regulatory factors, we have developed a physiologically based computational model of skeletal muscle energy metabolism. This model integrates transport and reaction fluxes in distinct capillary, cytosolic, and mitochondrial domains and investigates the roles of mitochondrial NADH/NAD+ transport (shuttling) activity and muscle glycogen concentration (stores) during moderate intensity exercise (60% maximal O2 consumption). The underlying hypothesis is that the cytosolic redox state (NADH/NAD+) is much more sensitive to a metabolic disturbance in contracting skeletal muscle than the mitochondrial redox state. This hypothesis was tested by simulating the dynamic metabolic responses of skeletal muscle to exercise while altering the transport rate of reducing equivalents (NADH and NAD+) between cytosol and mitochondria and muscle glycogen stores. Simulations with optimal parameter estimates showed good agreement with the available experimental data from muscle biopsies in human subjects. Compared with these simulations, a 20% increase (or approximately 20% decrease) in mitochondrial NADH/NAD+ shuttling activity led to an approximately 70% decrease (or approximately 3-fold increase) in cytosolic redox state and an approximately 35% decrease (or approximately 25% increase) in muscle lactate level. Doubling (or halving) muscle glycogen concentration resulted in an approximately 50% increase (or approximately 35% decrease) in cytosolic redox state and an approximately 30% increase (or approximately 25% decrease) in muscle lactate concentration. In both cases, changes in mitochondrial redox state were minimal. In conclusion, the model simulations of exercise response are consistent with the hypothesis that mitochondrial NADH/NAD+ shuttling activity and muscle glycogen stores affect primarily the cytosolic redox state. Furthermore, muscle lactate production is regulated primarily by the cytosolic redox state.

Altered NAD(P)H production in neonatal rat islets resistant to H2O2

AIMS

We determined the involvement of NAD(P)H generation ability on the resistance of pancreatic islets B-cells to oxidative stress caused by culture exposition to H2O2.

MAIN METHODS

We cultured isolated neonatal Wistar rat islets for four days in medium containing 5.6 or 20 mM glucose, with or without H2O2 (200 microM), and analyzed several parameters associated with islet survival in different media. High glucose was used since it protects neonatal islets against the loss of GSIS.

KEY FINDINGS

While none of the culture conditions increased the rate of NAD(P)H content at 16.7 mM glucose, the islets resistant to H2O2 and those exposed to 20 mM glucose showed a greater use of the pentose phosphate pathway and increased ATP synthesis from glucose.

SIGNIFICANCE

Oxidative stress contributes to the loss of glucose-induced insulin secretion (GSIS) during the onset of diabetes mellitus. Although immature rat islets have reduced GSIS compared to mature islets, they adapt better to oxidative stress and are a good model for understanding the causes involved in the destruction or survival of islet cells. These data support the idea that GSIS and resistance against oxidative stress in immature islets rely on NADH shuttle activities, with little contribution of reduced equivalents from the tricarboxylic acid cycle (TCAC).

A kinetic core model of the glucose-stimulated insulin secretion network of pancreatic beta cells

The construction and characterization of a core kinetic model of the glucose-stimulated insulin secretion system (GSIS) in pancreatic beta cells is described. The model consists of 44 enzymatic reactions, 59 metabolic state variables, and 272 parameters. It integrates five subsystems: glycolysis, the TCA cycle, the respiratory chain, NADH shuttles, and the pyruvate cycle. It also takes into account compartmentalization of the reactions in the cytoplasm and mitochondrial matrix. The model shows expected behavior in its outputs, including the response of ATP production to starting glucose concentration and the induction of oscillations of metabolite concentrations in the glycolytic pathway and in ATP and ADP concentrations. Identification of choke points and parameter sensitivity analysis indicate that the glycolytic pathway, and to a lesser extent the TCA cycle, are critical to the proper behavior of the system, while parameters in other components such as the respiratory chain are less critical. Notably, however, sensitivity analysis identifies the first reactions of nonglycolytic pathways as being important for the behavior of the system. The model is robust to deletion of malic enzyme activity, which is absent in mouse pancreatic beta cells. The model represents a step toward the construction of a model with species-specific parameters that can be used to understand mouse models of diabetes and the relationship of these mouse models to the human disease state.