Gap-junction coupling and ATP-sensitive potassium channels in human β-cell clusters: Effects on emergent dynamics

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.

From Isles of Königsberg to Islets of Langerhans: Examining the Function of the Endocrine Pancreas Through Network Science

Islets of Langerhans are multicellular microorgans located in the pancreas that play a central role in whole-body energy homeostasis. Through secretion of insulin and other hormones they regulate postprandial storage and interprandial usage of energy-rich nutrients. In these clusters of hormone-secreting endocrine cells, intricate cell-cell communication is essential for proper function. Electrical coupling between the insulin-secreting beta cells through gap junctions composed of connexin36 is particularly important, as it provides the required, most important, basis for coordinated responses of the beta cell population. The increasing evidence that gap-junctional communication and its modulation are vital to well-regulated secretion of insulin has stimulated immense interest in how subpopulations of heterogeneous beta cells are functionally arranged throughout the islets and how they mediate intercellular signals. In the last decade, several novel techniques have been proposed to assess cooperation between cells in islets, including the prosperous combination of multicellular imaging and network science. In the present contribution, we review recent advances related to the application of complex network approaches to uncover the functional connectivity patterns among cells within the islets. We first provide an accessible introduction to the basic principles of network theory, enumerating the measures characterizing the intercellular interactions and quantifying the functional integration and segregation of a multicellular system. Then we describe methodological approaches to construct functional beta cell networks, point out possible pitfalls, and specify the functional implications of beta cell network examinations. We continue by highlighting the recent findings obtained through advanced multicellular imaging techniques supported by network-based analyses, giving special emphasis to the current developments in both mouse and human islets, as well as outlining challenges offered by the multilayer network formalism in exploring the collective activity of islet cell populations. Finally, we emphasize that the combination of these imaging techniques and network-based analyses does not only represent an innovative concept that can be used to describe and interpret the physiology of islets, but also provides fertile ground for delineating normal from pathological function and for quantifying the changes in islet communication networks associated with the development of diabetes mellitus.

Hubs, diversity, and synchronization in FitzHugh-Nagumo oscillator networks: Resonance effects and biophysical implications

Using the FitzHugh-Nagumo equations to represent the oscillatory electrical behavior of β-cells, we develop a coupled oscillator network model with cubic lattice topology, showing that the emergence of pacemakers or hubs in the system can be viewed as a natural consequence of oscillator population diversity. The optimal hub to nonhub ratio is determined by the position of the diversity-induced resonance maximum for a given set of FitzHugh-Nagumo equation parameters and is predicted by the model to be in a range that is fully consistent with experimental observations. The model also suggests that hubs in a β-cell network should have the ability to "switch on" and "off" their pacemaker function. As a consequence, their relative amount in the population can vary in order to ensure an optimal oscillatory performance of the network in response to environmental changes, such as variations of an external stimulus.

Biophysical modeling of β-cells networks: Realistic architectures and heterogeneity effects

The β-cells dynamics is the regulator of insulin secretion in the pancreas, and its investigation is a central aspect in designing effective treatment strategies for diabetes. Despite great efforts, much is still unknown about the complex organization of such endocrine cells and realistic mathematical modeling represents a useful tool to elucidate key aspects of glucose control in humans. In this contribution, we study the human β-cells collective behaviour, by modeling their electric and metabolic coupling in a cluster, of size and architecture similar to human islets of Langerhans. We focus on the effect of coupling on various dynamics regimes observed in the islets, that are spiking and bursting on multiple timescales. In particular, we test the effect of hubs, that are highly glucose-sensitive β-cells, on the overall network dynamics, observing different modulation depending on the timescale of the dynamics. By properly taking into account the role of cells heterogeneity, recently emerged, our model effectively describes the effect of hubs on the synchronization of the islet response and the correlation of β-cells activity.

Biophysical modeling of C. elegans neurons: Single ion currents and whole-cell dynamics of AWCon and RMD

C. elegans neuronal system constitutes the ideal framework for studying simple, yet realistic, neuronal activity, since the whole nervous system is fully characterized with respect to the exact number of neurons and the neuronal connections. Most recent efforts are devoted to investigate and clarify the signal processing and functional connectivity, which are at the basis of sensing mechanisms, signal transmission, and motor control. In this framework, a refined modelof whole neuron dynamics constitutes a key ingredient to describe the electrophysiological processes, both at thecellular and at the network scale. In this work, we present Hodgkin-Huxley-based models of ion channels dynamics black, built on data available both from C. elegans and from other organisms, expressing homologous channels. We combine these channel models to simulate the electrical activity oftwo among the most studied neurons in C. elegans, which display prototypical dynamics of neuronal activation, the chemosensory AWCON and the motor neuron RMD. Our model properly describes the regenerative responses of the two cells. We analyze in detail the role of ion currents, both in wild type and in in silico knockout neurons. Moreover, we specifically investigate the behavior of RMD, identifying a heterogeneous dynamical response which includes bistable regimes and sustained oscillations. We are able to assess the critical role of T-type calcium currents, carried by CCA-1 channels, and leakage currents in the regulation of RMD response. Overall, our results provide new insights in the activity of key C. elegans neurons. The developed mathematical framework constitute a basis for single-cell and neuronal networks analyses, opening new scenarios in the in silico modeling of C. elegans neuronal system.

Critical transitions in heterogeneous networks: Loss of low-degree nodes as an early warning signal

A large number of real networks show abrupt phase transition phenomena in response to environmental changes. In this case, cascading phenomena can induce drastic and discontinuous changes in the system state and lead to collapse. Although complex network theory has been used to investigate these drastic events, we are still unable to predict them effectively. We here analyze collapse phenomena by proposing a minimal two-state dynamic on a complex network and introducing the effect of local connectivities on the evolution of network nodes. We find that a heterogeneous system of interconnected components presents a mixed response to stress and can serve as a control indicator. In particular, before the critical transition point is reached a severe loss of low-degree nodes is observed, masked by the minimal failure of higher-degree nodes. Accordingly, we suggest that a significant reduction in less connected nodes can indicate impending global failure.