The hyperbolic effect of density and strength of inter beta-cell coupling on islet bursting: a theoretical investigation

A primer on modelling pancreatic islets: from models of coupled β-cells to multicellular islet models

Pancreatic islets are mini-organs composed of hundreds or thousands of ɑ, β and δ-cells, which, respectively, secrete glucagon, insulin and somatostatin, key hormones for the regulation of blood glucose. In pancreatic islets, hormone secretion is tightly regulated by both internal and external mechanisms, including electrical communication and paracrine signaling between islet cells. Given its complexity, the experimental study of pancreatic islets has been complemented with computational modeling as a tool to gain a better understanding about how all the mechanisms involved at different levels of organization interact. In this review, we describe how multicellular models of pancreatic cells have evolved from the early models of electrically coupled β-cells to models in which experimentally derived architectures and both electrical and paracrine signals have been considered.

Multilevel synchronization of human β-cells networks

β-cells within the endocrine pancreas are fundamental for glucose, lipid and protein homeostasis. Gap junctions between cells constitute the primary coupling mechanism through which cells synchronize their electrical and metabolic activities. This evidence is still only partially investigated through models and numerical simulations. In this contribution, we explore the effect of combined electrical and metabolic coupling in β-cell clusters using a detailed biophysical model. We add heterogeneity and stochasticity to realistically reproduce β-cell dynamics and study networks mimicking arrangements of β-cells within human pancreatic islets. Model simulations are performed over different couplings and heterogeneities, analyzing emerging synchronization at the membrane potential, calcium, and metabolites levels. To describe network synchronization, we use the formalism of multiplex networks and investigate functional network properties and multiplex synchronization motifs over the structural, electrical, and metabolic layers. Our results show that metabolic coupling can support slow wave propagation in human islets, that combined electrical and metabolic synchronization is realized in small aggregates, and that metabolic long-range correlation is more pronounced with respect to the electrical one.

Integrated multiscale mathematical modeling of insulin secretion reveals the role of islet network integrity for proper oscillatory glucose-dose response

The integrated multiscale mathematical model we present in this paper is built on two of our previous ones: a model of electrical oscillation in β-cells connected to neighboring cells within a three-dimensional (3D) network, and a model of glucose-induced β-cell intracellular insulin granule trafficking and insulin secretion. In order to couple these two models, we assume that the rate at which primed and release-ready insulin granules fuse at the cell membrane increases with the intracellular calcium concentration. Moreover, by assuming that the fraction of free K_ATP-channels decreases with increasing glucose concentration, we take into account the effect of glucose dose on membrane potential and, indirectly via the effect on the potential, on intracellular calcium. Numerical analysis of our new model shows that a single step increase in glucose concentration yields the experimentally observed characteristic biphasic insulin release. We find that the biphasic response is typically oscillatory in nature for low and moderate glucose concentrations. The plateau fraction (the time that the β-cells spend in their active firing phase) increases with increasing glucose dose, as does the total insulin secretion. At high glucose concentrations, the oscillations tend to vanish due to a constantly elevated membrane potential of the β-cells. Our results also demonstrate how insulin secretion characteristics in various glucose protocols depend on the degree of β-cell loss, highlighting the potential impact from disease. In particular, both the secretory capacity (average insulin secretion rate per β-cell) and the oscillatory response diminish as the islet cell network becomes compromised.

Investigation of the role of ion channels in human pancreatic β-cell hubs: A mathematical modeling study

In many cellular networks, the structure of the network follows a scale-free organization, where a limited number of cells are strongly coupled to other cells. These cells are called hub cells and their critical roles are well accepted. Despite their importance, there have been only a few studies investigating the characteristic features of these cells. In this paper, a computational approach is proposed to study the possible role of different ion channels in distinguishing between the hub and non-hub cells. The results show that the P/Q-type and T-type calcium channels may have an especial role in the β-cell hubs because the high-level expressions of these channels make a pancreatic β-cell more potent to force other coupled cells to follow it. In addition, in order to consider the variation of the coupling strength with voltage, a novel mathematical model is proposed for the gap junction coupling between the pancreatic β-cells. The proposed approach is validated based on the data from the literature.

Role of topology in complex functional networks of beta cells

The activity of pancreatic β cells can be described by biological networks of coupled nonlinear oscillators that, via electrochemical synchronization, release insulin in response to augmented glucose levels. In this work, we analyze the emergent behavior of regular and percolated β-cells clusters through a stochastic mathematical model where "functional" networks arise. We show that the emergence and robustness of the synchronized dynamics depend both on intrinsic and extrinsic parameters. In particular, cellular noise level, glucose concentration, network spatial architecture, and cell-to-cell coupling strength are the key factors for the generation of a rhythmic and robust activity. Their role in the functional network topology associated with β-cells clusters is analyzed and discussed.

Development, growth and maintenance of β-cell mass: models are also part of the story

Pancreatic β-cells in the islets of Langerhans play a crucial role in regulating glucose homeostasis in the circulation. Loss of β-cell mass or function due to environmental, genetic and immunological factors leads to the manifestation of diabetes mellitus. The mechanisms regulating the dynamics of pancreatic β-cell mass during normal development and diabetes progression are complex. To fully unravel such complexity, experimental and clinical approaches need to be combined with mathematical and computational models. In the natural sciences, mathematical and computational models have aided the identification of key mechanisms underlying the behavior of systems comprising multiple interacting components. A number of mathematical and computational models have been proposed to explain the development, growth and death of pancreatic β-cells. In this review, we discuss some of these models and how their predictions provide novel insight into the mechanisms controlling β-cell mass during normal development and diabetes progression. Lastly, we discuss a handful of the major open questions in the field.

Phase transitions in pancreatic islet cellular networks and implications for type-1 diabetes

In many aspects the onset of a chronic disease resembles a phase transition in a complex dynamic system: Quantitative changes accumulate largely unnoticed until a critical threshold is reached, which causes abrupt qualitative changes of the system. In this study we examine a special case, the onset of type-1 diabetes (T1D), a disease that results from loss of the insulin-producing pancreatic islet β cells. Within each islet, the β cells are electrically coupled to each other via gap-junctional channels. This intercellular coupling enables the β cells to synchronize their insulin release, thereby generating the multiscale temporal rhythms in blood insulin that are critical to maintaining blood glucose homeostasis. Using percolation theory we show how normal islet function is intrinsically linked to network connectivity. In particular, the critical amount of β-cell death at which the islet cellular network loses site percolation is consistent with laboratory and clinical observations of the threshold loss of β cells that causes islet functional failure. In addition, numerical simulations confirm that the islet cellular network needs to be percolated for β cells to synchronize. Furthermore, the interplay between site percolation and bond strength predicts the existence of a transient phase of islet functional recovery after onset of T1D and introduction of treatment, potentially explaining the honeymoon phenomenon. Based on these results, we hypothesize that the onset of T1D may be the result of a phase transition of the islet β-cell network.

A stochastic mathematical model to study the autoimmune progression towards type 1 diabetes

BACKGROUND

The integrity of the interactions and the 3D architecture among beta cell populations in pancreatic islets is critical for proper biosynthesis, storage and release of insulin. The aim of this study was to evaluate the effect on electrophysiological signalling of beta cells that is produced by progressive lymphocytic islet cell infiltration (insulitis), by modelling the disruption of pancreatic islet anatomy as a consequence of insulitis and altered glucose concentrations.

METHODS

On the basis of histopathological images of murine islets from non-obese diabetic mice, we simulated the electrophysiological dynamics of a 3D cluster of mouse beta cells via a stochastic model. Progressive damage was modelled at different glucose concentrations, representing the different glycaemic states in the autoimmune progression towards type 1 diabetes.

RESULTS

At 31% of dead beta cells (normoglycaemia) and 69% (hyperglycaemia), the system appeared to be biologically robust to maintain regular Ca(2+) ion oscillations guaranteeing an effective insulin release. Simulations at 84%, 94% and 98% grades (severe hyperglycemia) showed that intracellular calcium oscillations were absent. In such conditions, insulin pulsatility is not expected to occur.

CONCLUSIONS

Our results suggest that the islet tissue is biophysically robust enough to compensate for high rates of beta cell loss. These predictions can be experimentally tested in vitro by quantifying space and time electrophysiological dynamics of animal islets kept at different glucose gradients. The model indicates the necessity of maintaining glycaemia within the physiological range as soon as possible after diabetes onset to avoid a dramatic interruption of Ca(2+) pulsatility and the consequent drop of insulin release.

Mathematical modeling of insulin secretion and the role of glucose-dependent mobilization, docking, priming and fusion of insulin granules

In this paper we develop a new mathematical model of glucose-induced insulin secretion from pancreatic islet β-cells, and we use this model to investigate the rate limiting factors. We assume that insulin granules reside in different pools inside each β-cell, and that all β-cells respond homogeneously to glucose with the same recruitment thresholds. Consistent with recent experimental observations, our model also accounts for the fusion of newcomer granules that are not pre-docked at the plasma membrane. In response to a single step increase in glucose concentration, our model reproduces the characteristic biphasic insulin release observed in multiple experimental systems, including perfused pancreata and isolated islets of rodent or human origin. From our model analysis we note that first-phase insulin secretion depends on rapid depletion of the primed, release-ready granule pools, while the second phase relies on granule mobilization from the reserve. Moreover, newcomers have the potential to contribute significantly to the second phase. When the glucose protocol consists of multiple changes in sequence (a so-called glucose staircase), our model predicts insulin spikes of increasing height, as has been seen experimentally. This increase stems from the glucose-dependent increase in the fusion rate of insulin granules at the plasma membrane of single β-cells. In contrast, previous mathematical models reproduced the staircase experiment by assuming heterogeneous β-cell activation. In light of experimental data indicating limited heterogeneous activation for β-cells within intact islets, our findings suggest that a graded, dose-dependent cell response to glucose may contribute to insulin secretion patterns observed in multiple experiments, and thus regulate in vivo insulin release. In addition, the strength of insulin granule mobilization, priming and fusion are critical limiting factors in determining the total amount of insulin release.