Stable and flexible system for glucose homeostasis

Heteroclinic switching between chimeras in a ring of six oscillator populations

In a network of coupled oscillators, a symmetry-broken dynamical state characterized by the coexistence of coherent and incoherent parts can spontaneously form. It is known as a chimera state. We study chimera states in a network consisting of six populations of identical Kuramoto-Sakaguchi phase oscillators. The populations are arranged in a ring, and oscillators belonging to one population are uniformly coupled to all oscillators within the same population and to those in the two neighboring populations. This topology supports the existence of different configurations of coherent and incoherent populations along the ring, but all of them are linearly unstable in most of the parameter space. Yet, chimera dynamics is observed from random initial conditions in a wide parameter range, characterized by one incoherent and five synchronized populations. These observable states are connected to the formation of a heteroclinic cycle between symmetric variants of saddle chimeras, which gives rise to a switching dynamics. We analyze the dynamical and spectral properties of the chimeras in the thermodynamic limit using the Ott-Antonsen ansatz and in finite-sized systems employing Watanabe-Strogatz reduction. For a heterogeneous frequency distribution, a small heterogeneity renders a heteroclinic switching dynamics asymptotically attracting. However, for a large heterogeneity, the heteroclinic orbit does not survive; instead, it is replaced by a variety of attracting chimera states.

Chaotic chimera attractors in a triangular network of identical oscillators

A prominent type of collective dynamics in networks of coupled oscillators is the coexistence of coherently and incoherently oscillating domains known as chimera states. Chimera states exhibit various macroscopic dynamics with different motions of the Kuramoto order parameter. Stationary, periodic and quasiperiodic chimeras are known to occur in two-population networks of identical phase oscillators. In a three-population network of identical Kuramoto-Sakaguchi phase oscillators, stationary and periodic symmetric chimeras were previously studied on a reduced manifold in which two populations behaved identically [Phys. Rev. E 82, 016216 (2010)1539-375510.1103/PhysRevE.82.016216]. In this paper, we study the full phase space dynamics of such three-population networks. We demonstrate the existence of macroscopic chaotic chimera attractors that exhibit aperiodic antiphase dynamics of the order parameters. We observe these chaotic chimera states in both finite-sized systems and the thermodynamic limit outside the Ott-Antonsen manifold. The chaotic chimera states coexist with a stable chimera solution on the Ott-Antonsen manifold that displays periodic antiphase oscillation of the two incoherent populations and with a symmetric stationary chimera solution, resulting in tristability of chimera states. Of these three coexisting chimera states, only the symmetric stationary chimera solution exists in the symmetry-reduced manifold.

Optimal Transport Flows for Distributed Production Networks

Network flows often exhibit a hierarchical treelike structure that can be attributed to the minimization of dissipation. The common feature of such systems is a single source and multiple sinks (or vice versa). In contrast, here we study networks with only a single source and sink. These systems can arise from secondary purposes of the networks, such as blood sugar regulation through insulin production. Minimization of dissipation in these systems leads to vascular shunting, a single vessel connecting the inlet and outlet. We show instead how optimizing the transport time yields network topologies that match those observed in the insulin-producing pancreatic islets. These are patterns of periphery-to-center and center-to-periphery flows. The obtained flow networks are broadly independent of how the flow velocity depends on the flow flux, but continuous and discontinuous phase transitions appear at extreme flux dependencies. Lastly, we show how constraints on flows can lead to buckling of the branches of the network, a feature that is also observed in pancreatic islets.

Tripartite cell networks for glucose homeostasis

Controlling the excess and shortage of energy is a fundamental task for living organisms. Diabetes is a representative metabolic disease caused by the malfunction of energy homeostasis. The islets of Langerhans in the pancreas release long-range messengers, hormones, into the blood to regulate the homeostasis of the primary energy fuel, glucose. The hormone and glucose levels in the blood show rhythmic oscillations with a characteristic period of 5-10 min, and the functional roles of the oscillations are not clear. Each islet has [Formula: see text] and [Formula: see text] cells that secrete glucagon and insulin, respectively. These two counter-regulatory hormones appear sufficient to increase and decrease glucose levels. However, pancreatic islets have a third cell type, [Formula: see text] cells, which secrete somatostatin. The three cell populations have a unique spatial organization in islets, and they interact to perturb their hormone secretions. The mini-organs of islets are scattered throughout the exocrine pancreas. Considering that the human pancreas contains approximately a million islets, the coordination of hormone secretion from the multiple sources of islets and cells within the islets should have a significant effect on human physiology. In this review, we introduce the hierarchical organization of tripartite cell networks, and recent biophysical modeling to systematically understand the oscillations and interactions of [Formula: see text], [Formula: see text], and [Formula: see text] cells. Furthermore, we discuss the functional roles and clinical implications of hormonal oscillations and their phase coordination for the diagnosis of type II diabetes.

Lentiviral Mediated Gene Silencing in Human Pseudoislet Prepared in Low Attachment Plates

Various genetic tools are available to modulate genes in pancreatic islets of rodents to dissect function of islet genes for diabetes research. However, the data obtained from rodent islets are often not fully reproduced in or applicable to human islets due to well-known differences in islet structure and function between the species. Currently, techniques that are available to manipulate gene expression of human islets are very limited. Introduction of transgene into intact islets by adenovirus, plasmid, and oligonucleotides often suffers from low efficiency and high toxicity. Low efficiency is especially problematic in gene downregulation studies in intact islets, which require high efficiency. It has been known that enzymatically-dispersed islet cells reaggregate in culture forming spheroids termed pseudoislets. Size-controlled reaggregation of human islet cells creates pseudoislets that maintain dynamic first phase insulin secretion after prolonged culture and provide a window to efficiently introduce lentiviral short hairpin RNA (shRNA) with low toxicity. Here, a detailed protocol for the creation of human pseudoislets after lentiviral transduction using two commercially available multiwell plates is described. The protocol can be easily performed and allows for efficient downregulation of genes and assessment of dynamism of insulin secretion using human islet cells. Thus, human pseudoislets with lentiviral mediated gene modulation provide a powerful and versatile model to assess gene function within human islet cells.

Rhythmic synchronization and hybrid collective states of globally coupled oscillators

Macroscopic rhythms are often signatures of healthy functioning in living organisms, but they are still poorly understood on their microscopic bases. Globally interacting oscillators with heterogeneous couplings are here considered. Thorough theoretical and numerical analyses indicate the presence of multiple phase transitions between different collective states, with regions of bi-stability. Novel coherent phases are unveiled, and evidence is given of the spontaneous emergence of macroscopic rhythms where oscillators’ phases are always found to be self-organized as in Bellerophon states, i.e. in multiple clusters with quantized values of their average frequencies. Due to their rather unconditional appearance, the circumstance is paved that the Bellerophon states grasp the microscopic essentials behind collective rhythms in more general systems of interacting oscillators.

Design principles of the paradoxical feedback between pancreatic alpha and beta cells

Mammalian glucose homeostasis is controlled by the antagonistic hormones insulin and glucagon, secreted by pancreatic beta and alpha cells respectively. These two cell types are adjacently located in the islets of Langerhans and affect each others’ secretions in a paradoxical manner: while insulin inhibits glucagon secretion from alpha cells, glucagon seems to stimulate insulin secretion from beta cells. Here we ask what are the design principles of this negative feedback loop. We systematically simulate the dynamics of all possible islet inter-cellular connectivity patterns and analyze different performance criteria. We find that the observed circuit dampens overshoots of blood glucose levels after reversion of glucose drops. This feature is related to the temporal delay in the rise of insulin concentrations in peripheral tissues, compared to the immediate hormone action on the liver. In addition, we find that the circuit facilitates coordinate secretion of both hormones in response to protein meals. Our study highlights the advantages of a paradoxical paracrine feedback loop in maintaining metabolic homeostasis.

Finite-size scaling in the system of coupled oscillators with heterogeneity in coupling strength

We consider a mean-field model of coupled phase oscillators with random heterogeneity in the coupling strength. The system that we investigate here is a minimal model that contains randomness in diverse values of the coupling strength, and it is found to return to the original Kuramoto model [Y. Kuramoto, Prog. Theor. Phys. Suppl. 79, 223 (1984)10.1143/PTPS.79.223] when the coupling heterogeneity disappears. According to one recent paper [H. Hong, H. Chaté, L.-H. Tang, and H. Park, Phys. Rev. E 92, 022122 (2015)10.1103/PhysRevE.92.022122], when the natural frequency of the oscillator in the system is "deterministically" chosen, with no randomness in it, the system is found to exhibit the finite-size scaling exponent ν[over ¯]=5/4. Also, the critical exponent for the dynamic fluctuation of the order parameter is found to be given by γ=1/4, which is different from the critical exponents for the Kuramoto model with the natural frequencies randomly chosen. Originally, the unusual finite-size scaling behavior of the Kuramoto model was reported by Hong et al. [H. Hong, H. Chaté, H. Park, and L.-H. Tang, Phys. Rev. Lett. 99, 184101 (2007)10.1103/PhysRevLett.99.184101], where the scaling behavior is found to be characterized by the unusual exponent ν[over ¯]=5/2. On the other hand, if the randomness in the natural frequency is removed, it is found that the finite-size scaling behavior is characterized by a different exponent, ν[over ¯]=5/4 [H. Hong, H. Chaté, L.-H. Tang, and H. Park, Phys. Rev. E 92, 022122 (2015)10.1103/PhysRevE.92.022122]. Those findings brought about our curiosity and led us to explore the effects of the randomness on the finite-size scaling behavior. In this paper, we pay particular attention to investigating the finite-size scaling and dynamic fluctuation when the randomness in the coupling strength is considered.

Local complexity predicts global synchronization of hierarchically networked oscillators

We study the global synchronization of hierarchically-organized Stuart-Landau oscillators, where each subsystem consists of three oscillators with activity-dependent couplings. We considered all possible coupling signs between the three oscillators, and found that they can generate different numbers of phase attractors depending on the network motif. Here, the subsystems are coupled through mean activities of total oscillators. Under weak inter-subsystem couplings, we demonstrate that the synchronization between subsystems is highly correlated with the number of attractors in uncoupled subsystems. Among the network motifs, perfect anti-symmetric ones are unique to generate both single and multiple attractors depending on the activities of oscillators. The flexible local complexity can make global synchronization controllable.

Synchronization and Bellerophon states in conformist and contrarian oscillators

The study of synchronization in generalized Kuramoto models has witnessed an intense boost in the last decade. Several collective states were discovered, such as partially synchronized, chimera, π or traveling wave states. We here consider two populations of globally coupled conformist and contrarian oscillators (with different, randomly distributed frequencies), and explore the effects of a frequency-dependent distribution of the couplings on the collective behaviour of the system. By means of linear stability analysis and mean-field theory, a series of exact solutions is extracted describing the critical points for synchronization, as well as all the emerging stationary coherent states. In particular, a novel non-stationary state, here named as Bellerophon state, is identified which is essentially different from all other coherent states previously reported in the Literature. A robust verification of the rigorous predictions is supported by extensive numerical simulations.

Paracrine regulation of glucagon secretion: the β/α/δ model

The regulation of glucagon secretion in the pancreatic α-cell is not well understood. It has been proposed that glucose suppresses glucagon secretion either directly through an intrinsic mechanism within the α-cell or indirectly through an extrinsic mechanism. Previously, we described a mathematical model for isolated pancreatic α-cells and used it to investigate possible intrinsic mechanisms of regulating glucagon secretion. We demonstrated that glucose can suppress glucagon secretion through both ATP-dependent potassium channels (K_ATP) and a store-operated current (SOC). We have now developed an islet model that combines previously published mathematical models of α- and β-cells with a new model of δ-cells and use it to explore the effects of insulin and somatostatin on glucagon secretion. We show that the model can reproduce experimental observations that the inhibitory effect of glucose remains even when paracrine modulators are no longer acting on the α-cell. We demonstrate how paracrine interactions can either synchronize α- and δ-cells to produce pulsatile oscillations in glucagon and somatostatin secretion or fail to do so. The model can also account for the paradoxical observation that glucagon can be out of phase with insulin, whereas α-cell calcium is in phase with insulin. We conclude that both paracrine interactions and the α-cell's intrinsic mechanisms are needed to explain the response of glucagon secretion to glucose.

Design Principles of Pancreatic Islets: Glucose-Dependent Coordination of Hormone Pulses

Pancreatic islets are functional units involved in glucose homeostasis. The multicellular system comprises three main cell types; β and α cells reciprocally decrease and increase blood glucose by producing insulin and glucagon pulses, while the role of δ cells is less clear. Although their spatial organization and the paracrine/autocrine interactions between them have been extensively studied, the functional implications of the design principles are still lacking. In this study, we formulated a mathematical model that integrates the pulsatility of hormone secretion and the interactions and organization of islet cells and examined the effects of different cellular compositions and organizations in mouse and human islets. A common feature of both species was that islet cells produced synchronous hormone pulses under low- and high-glucose conditions, while they produced asynchronous hormone pulses under normal glucose conditions. However, the synchronous coordination of insulin and glucagon pulses at low glucose was more pronounced in human islets that had more α cells. When β cells were selectively removed to mimic diabetic conditions, the anti-synchronicity of insulin and glucagon pulses was deteriorated at high glucose, but it could be partially recovered when the re-aggregation of remaining cells was considered. Finally, the third cell type, δ cells, which introduced additional complexity in the multicellular system, prevented the excessive synchronization of hormone pulses. Our computational study suggests that controllable synchronization is a design principle of pancreatic islets.

Traveling wave in a three-dimensional array of conformist and contrarian oscillators

We consider a system of conformist and contrarian oscillators coupled locally in a three-dimensional cubic lattice and explore collective behavior of the system. The conformist oscillators attractively interact with the neighbor oscillators and therefore tend to be aligned with the neighbors' phase. The contrarian oscillators interact repulsively with the neighbors and therefore tend to be out of phase with them. In this paper, we investigate whether many peculiar dynamics that have been observed in the mean-field system with global coupling can emerge even with local coupling. In particular, we pay attention to the possibility that a traveling wave may arise. We find that the traveling wave occurs due to coupling asymmetry and not by global coupling; this observation confirms that the global coupling is not essential to the occurrence of a traveling wave in the system. The traveling wave can be a mechanism for the coherent rhythm generation of the circadian clock or of hormone secretion in biological systems under local coupling.

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.

A conserved rule for pancreatic islet organization

Morphogenesis, spontaneous formation of organism structure, is essential for life. In the pancreas, endocrine , , and cells are clustered to form islets of Langerhans, the critical micro-organ for glucose homeostasis. The spatial organization of endocrine cells in islets looks different between species. Based on the three-dimensional positions of individual cells in islets, we computationally inferred the relative attractions between cell types, and found that the attractions between homotypic cells were slightly, but significantly, stronger than the attractions between heterotypic cells commonly in mouse, pig, and human islets. The difference between cell attraction and cell attraction was minimal in human islets, maximizing the plasticity of islet structures. Our result suggests that although the cellular composition and attractions of pancreatic endocrine cells are quantitatively different between species, the physical mechanism of islet morphogenesis may be evolutionarily conserved.

Periodic synchronization and chimera in conformist and contrarian oscillators

We consider a system of phase oscillators that couple with both attractive and repulsive interaction under a pinning force and explore collective behavior of the system. The oscillators can be divided into two subpopulations of "conformist" oscillators with attractive interaction and "contrarian" ones with repulsive interaction. We find that the interplay between the pinning force and the opposite relationship of the conformist and contrarian oscillators induce peculiar dynamic states: periodic synchronization, breathing chimera, and fully pinned state depending on the fraction of the conformists. Using the Watanabe-Strogatz transformation, we reduce the dynamics into a low-dimensional one and find that the above dynamic states are generated from the reduced dynamics.