The effects of beta-cell mass and function, intercellular coupling, and islet synchrony on [Formula: see text] dynamics.

Type 2 diabetes (T2D) is a challenging metabolic disorder characterized by a substantial loss of [Formula: see text]-cell mass and alteration of [Formula: see text]-cell function in the islets of Langerhans, disrupting insulin secretion and glucose homeostasis. The mechanisms for deficiency in [Formula: see text]-cell mass and function during the hyperglycemia development and T2D pathogenesis are complex. To study the relative contribution of [Formula: see text]-cell mass to [Formula: see text]-cell function in T2D, we make use of a comprehensive electrophysiological model of human [Formula: see text]-cell clusters. We find that defect in [Formula: see text]-cell mass causes a functional decline in single [Formula: see text]-cell, impairment in intra-islet synchrony, and changes in the form of oscillatory patterns of membrane potential and intracellular [Formula: see text] concentration, which can lead to changes in insulin secretion dynamics and in insulin levels. The model demonstrates a good correspondence between suppression of synchronizing electrical activity and published experimental measurements. We then compare the role of gap junction-mediated electrical coupling with both [Formula: see text]-cell synchronization and metabolic coupling in the behavior of [Formula: see text] concentration dynamics within human islets. Our results indicate that inter-[Formula: see text]-cellular electrical coupling depicts a more important factor in shaping the physiological regulation of islet function and in human T2D. We further predict that varying the whole-cell conductance of delayed rectifier [Formula: see text] channels modifies oscillatory activity patterns of [Formula: see text]-cell population lacking intercellular coupling, which significantly affect [Formula: see text] concentration and insulin secretion.

Introduction

The human pancreatic -cells placed in islets of Langerhans organize a complex functional network[1-4], and ensure blood glucose homeostasis[5] through a pulsatile, well-regulated insulin secretion[6-8]. Like their rodent counterparts[9-11], human -cells respond to increasing plasma glucose concentrations with acceleration of metabolism and elevation in ATP levels, which in turn inhibits ATP-sensitive potassium () channels, leading to membrane depolarization, activation of voltage-dependent calcium channels (VDCCs), a rise in the cytosolic calcium concentration (), and triggering exocytosis of insulin granules[12-15].

Both murine and human dissociated -cells are intrinsically and functionally heterogeneous with respect to biophysical characteristics, metabolic responses, electrical dynamics, and insulin release[16-19]. Within rodent islets, gap junctions consisting of the Connexin36 (Cx36) protein form a strong intercellular electrical coupling between heterogeneous -cells, which is potentially important for coordination of oscillations in the -cell intracellular and insulin secretion across the islet, enhancing pulsatility of insulin secretion and regulating glucose homeostasis[20-24]. Similarly, in human pancreatic islets Cx36 gap junctions help to overcome the heterogeneity of individual -cells, and give rise to bursting behavior in synchrony upon glucose stimulation[24-26]. Additionally, Cx36 gap junction electrical coupling mediates a marked suppression of spontaneous elevations at basal glucose, generates a characteristic sigmoidal secretory response to increasing glucose, enhances the peak amplitude of first-phase insulin release, and coordinates the pulsatile second-phase insulin release[27-29], which in turn are important for glucose homeostasis. These observations illustrate a strong link between glucose-stimulated insulin secretion and gap junction function. Recent investigations have indicated the presence of heterogeneity in gap junctional conductance between the heterogeneous -cells[22,30,31] that influences the complex functional organization of -cells[32-34] and the spatiotemporal characteristics of waves under the islet mathematical models[35-38].

Besides electrical communications, metabolic communications are mediated by gap junction channels. These intercellular channels permit cell-to-cell diffusion of specific signaling ions and some glycolytic intermediates, which strongly affects the pattern of oscillations and insulin secretory profiles[39]. Intact islets display electrical behaviour consisting of so-called slow bursting pattern with a period of min, corresponding to the frequency of metabolism, , and insulin oscillations[40,41]. It is believed that this rhythmic islet activity is due to slow oscillatory dynamics of underlying glycolysis and metabolism observed in mouse and human -cells[42,43]. In particular, glycolytic oscillations are proposed to be the key player in the overall islet activity at stimulatory glucose concentrations.

Cx36 knockout mice display a reduction in the synchrony of oscillations and disruption of calcium wave propagation resulting in impaired pulsatile patterns of insulin release and glucose intolerance[28,37,44]. Importantly, these changes in activity and insulin secretion dynamics have also been observed in the prediabetes stage and patients with type 2 diabetes[45-49]. This, along with mouse models of prediabetes which have demonstrated Cx36 disruption characteristic at this stage[50], suggest that changes in Cx36-mediated coupling may be a key determining factor in islet dysfunction and development of T2D[21,51]. Previous reports have shown that the expression of Cx36 protein is decreased in islets exposed to chronic hyperglycemia[52]. Altered Cx36 gap junction function makes the pancreatic islets more sensitive to -cell damage and lower efficiency in insulin secretion[28,53], which again suggests a potential role for decreased coupling in T2D characterized by progressive -cell death.

Plasma insulin levels depend on the absolute number of insulin secreting -cells (i.e., -cell mass) and the functional status of each of these cells (i.e., -cell function). Hereby, deficiency in either dynamics of -cell mass or function, or both, results in insulin insufficiency and the onset of hyperglycemia. Current research in diabetes reveals that besides a substantial decline in -cell mass, a significant defect in -cell function is evident in T2D patients[54-57]. Furthermore, other studies indicate that at the time of T2D diagnosis diabetic islets seem to have lost of their -cells[56-58], which is tightly correlated with islet dysfunction including a reduction in the amplitude of first-phase insulin secretion and impairment in the secretory pulses during second-phase insulin secretion[43,48]. Clinically, the most common therapeutic approaches for T2D aim to regenerate -cell mass or to preserve -cell function. Addressing the latter needs a deep understanding of the contribution and kinetics of -cell mass and function in T2D etiology and pathogenesis. Thus, in the present study, we examined whether -cell defects are intrinsically functional or whether a reduction in -cell mass is linked to -cell dysfunction. We first constructed a multicellular computational model of heterogeneous and heterogeneously coupled -cells, and analyzed the effects on functional behavior of single -cell caused by changes in -cell population of human islets. We obtained a nice agreement between theoretical results and experimental data regarding the disruption in normal oscillatory patterns of insulin secretion after -cell loss. We then compared the behavior of electrical activity and dynamics after reduction in gap junction coupling. By combining the effect of varying coupling strengths and glucose stimulations, we investigated how levels altered with loss of and -cell mass. Finally, we predicted the impact of delayed rectifier (Kv) channels on the electrical behavior of uncoupled -cell population.

Computational methods

Model of -cell

The Hodgkin–Huxley type model for human -cells has been developed by Pedersen[59], who carefully described the electrophysiological properties of ion channels in human -cells, and then Riz et al. included dynamics in the model[60]. For this study, we prefer to use such a formulation because it provides a firm explanation for the human -cell dynamics, confirmed by experimental investigations[61,62]. The model is composed of an electrical component and a glycolytic component[63]. It includes membrane potential activity, cytosol and submembrane dynamics of , and glucose metabolism. The glycolytic oscillatory component drives slow bursting patterns[59,60,64], underlying slow oscillations and insulin release pulses[40,41].

Briefly, the membrane potential of a single -cell i follows:where denotes the transmembrane current conducted by channel type . Full equations and parameters of the model can be found in the Supplementary Material.

Network of -cells

There exists evidence that the mean -cell number for each human pancreatic islet is [65-67]. For lattice structure of islet, we model cubic network including -cells such that each central cell is surrounded by 6 neighbors (the Von Neumann neighborhoods of cellular automata theory). These cells are coupled with adjacent -cells through both electrical and metabolic connections. The equation for membrane potential of the ith -cell surrounded by j neighboring cells is modified to simulate electrical gap junction coupling in the modeled islet:where refers to the electrical coupling conductance between cells i and j, and is all adjacent cells of cell i.

As in another modeling study[68], to account for the metabolic coupling among -cells we consider the diffusion of glucose-6-phosphate (G6P) between cells, which is assumed to be in rapid equilibrium with fructose-6-phosphate (F6P). The equation for the total concentrations of G6P and F6P in the ith -cell surrounded by j neighboring cells is defined by:where is the glucokinase reaction rate of cell i, refers to the phospho-fructokinase reaction rate, and the parameter describes the metabolic coupling strength computed by analyzing data of diffusion of glycolytic metabolites among the islet -cells.

Numerical methods

All equations of the model are written and implemented in a Python algorithm, and the forth-order Runge–Kutta numerical scheme is used for solving the ODE systems, both electrical and metabolic components with a time-step of 0.02 and 0.05 ms, respectively.

In the present study, biological heterogeneity of the human -cells is introduced by some crucial parameters that control the electrical behavior of the modeled -cells. Specifically, cellular heterogeneity is represented in the conductances of the gap junction channels and the delayed rectifier channels , and the maximal reaction rate of glucokinase enzyme , which has strong effects on the slow oscillation frequency, and is also connected to heterogeneity glucose sensitivity. The values of parameters (and similarly ) and select from normal distributions with mean value equal to original value of the parameters and standard deviation is set to and of the mean, respectively.

Results

Does -cell loss of mass primarily cause functional -cell defect in type 2 diabetes?

To investigate the contribution of -cell mass and function to the insufficient insulin release and progression of T2D, we eliminated the modeled -cells in islet network randomly, in order to more closely mimic the clinical conditions. We observed that with removing the simulated cells, not only the summed activity was decreased (Fig. 1A) but surprisingly, of single -cell was reduced (Fig. 1B). These changes in intracellular levels of an active -cell were caused by substantial changes in the shape of electrical activity, in addition to impairment of coordinated electrical behavior in -cell islet. When of the islet cells were lost, the pattern of oscillatory membrane potential, which correlated with the pattern of oscillations, was noteworthy different from the deletion of -cell mass (Fig. 1E (I and II)). In fact, the peak level of was significantly lower in than loss due to widely varying patterns of the membrane potential oscillations, which occurred at a major decline in -cell mass. Therefore, changes in -cell mass caused the functional alterations in single -cell, leading to changes in insulin concentration and secretion dynamics. This result can support to the hypothesis that deficit in -cell mass induces various abnormalities in single -cell function observed in patients with type 2 diabetes.

Figure 1

Loss of - cell mass causes defect of -cell function. (A) Linear decrease in the summed activity during reduction of -cell mass. The number of cells reduced is represented as the of all cells in the modeled islet (1000 cells). (B) Mean activity during reduction of -cell mass. Note that deficiency in -cell mass has the strong effect on level. (C, D) Dependence of the synchronization of and membrane potential oscillations as a function of the number of cells reduced ( of islet). The intra-islet synchronization is disrupted after loss of -cell mass. (E) Representative simulated membrane potential and time-courses across the islet lacking -cells (subpanel (I)) and lacking -cells (subpanel (II)). A massive decline in -cell mass causes changes in the shape of electrical bursting. Data in (A–D) represent the mean ± S.E.M. over 10 modeled islets.

Additionally, we surprisingly noticed that -cell electrical activity across the islet, specially oscillations of -cell , completely synchronized before loss of the cells, and then the intra-islet synchronization dropped with loss of (Fig. 1C,D), resembling experimental recordings which disruption in pulsatile insulin secretion linked to impairment in coordination of oscillations can be found after reduction of the -cell mass in the pathogenesis of type 2 diabetes[56-58]. On the other hand, the summed activity was partly linearly decreased with deletion of -cells (Fig. 1A). Therefore, considering that the phase transition occurred in the intra-islet synchrony and that this behavior was not observed in the summed response, it seems that the lack of insulin pulsatility patterns is a more important factor in type 2 diabetes. These simulations confirm previous studies, which show that loss of -cell mass is a critical point[56,58].

Intercellular coupling is more effective than islet synchronization on concentrations

We simulated slow oscillations with a period of min, and compared the influence of gap junction and synchronization on levels of single human -cell. To examine whether intercellular coupling or islet synchronization is more effective on the calcium dynamics, we remove randomly gap junction coupling between -cells in the islet lattice (Fig. 2). With reduction in -cell coupling, the average cytosolic calcium concentration began to fall, initially (Fig. 2A) while after deletion of gap junction connections, the synchrony of electrical dynamics across the islet started to reduce (Fig. 2B,C). In other words, as less than connections among islet -cells were lost, of single -cell was gradually decreased whereas the whole-islet activity remained synchronous. This result suggests that gap junctional conductance of -cell is a more important factor than islet synchrony of electrical patterns for intracellular dynamics and levels of insulin secretion.

Figure 2

Gap junction coupling is more effective on activity. (A) Mean response across the modeled islet as a function of fraction of gap junctions removed. Note that level of decreases faster after loss of connections. (B, C) The synchrony of and membrane potential oscillations depends on fraction of gap junctions removed. Data represent the mean ± S.E.M. over 10 modeled islets.

Moreover, Fig. 2A displayed the dual behavior of the islet -cells: before deletion of gap junctions, the slope of average was less than after deletion of intercellular connections. In fact, level of single -cell decreased substantially faster after a greater () loss of cells compared with the loss of . These results indicate that a massive decline in the islet gap junction coupling, resulting in disruption of synchrony in human -cell islet, plays a key role in more fast decrease in the average , poorly coordinated calcium dynamics, and plausibly, impairment of pulsatile insulin release and development of diabetes.

The interplay of gap junction electrical coupling and metabolic coupling in activity

Loppini et al. considered metabolic coupling, in addition to electrical coupling, resulting in diffusion of glycolytic metabolites such as G6P and F6P molecules among the islet -cells[68]. We analyzed the role of electrical and metabolic coupling strength on the -cell calcium activity by graded changes in the gap junction conductance and the metabolic coupling strength.

In this model, we found that the mean behavior of for single -cell changed only slightly when fixing the parameter and elevating the value of (Fig. 3A). In other words, for any fixed the level of intracellular concentration did not depend strongly on . It should be noted that -cell was noticeably low at , because purely electrical coupling did not synchronize metabolic oscillations, giving rise to out-of-phase slow bursting and small amplitude oscillations. This predicted that to increase the concentration, there had to be metabolic diffusion between islet -cells, even a very small value.

Figure 3

Electrical coupling is more effective on activity. (A) The interplay between electrical coupling and metabolic coupling strength in the mean for all -cells in a human islet model. Note that gap junction electrical coupling is more effective on intracellular activity. (B) Representative simulated membrane potential and time-courses corresponding to three points (I), (II), and (III) indicated in (A). Changes in electrical behavior cause varying levels of concentration. (C, D) The interplay between electrical coupling and metabolic coupling strength in the synchronization of and membrane potential oscillations. Data in (A), (C), and (D) represent the mean ± S.E.M. over 10 modeled islets.

At fixed , the mean of single -cell significantly changed with varying value of such that the concentration of an active -cell was maximum in certain ranges of , and was then reduced (Fig. 3A). These results demonstrate that electrical coupling has a greater effect compare with metabolic coupling on the -cell activity in human islets.

Figure 3A displayed that changes in -cell was bimodal. Three points were indicated in the plot at which the average level was significantly different. These differences in occurred with changes in multicellular electrical behavior (Fig. 3B). In fact, different oscillations of the -cell membrane potential, imposing oscillation patterns, caused varying cytosolic calcium concentration. Furthermore, Fig. 3C and D depicted that electrical dynamics, especially intracellular activity, in the islet were almost fully synchronous whereas behavior was bimodal. Therefore, as we already mentioned, the islet gap junction coupling is mostly important for synchronization of oscillatory activity.

The interplay of gap junction electrical coupling and glucose stimulation in activity

We examined the impact of islet gap junction coupling, and elevated glucose concentrations on the dynamics of intracellular . At low levels of glucose , despite increasing intercellular coupling strength and nearly complete synchrony of electrical activity, -cells showed globally quiescent behavior, and were inexcitable within the islet (Fig. 4 A and B and C). In other words, the coupling strength between human -cells had no considerable affect on concentration of single -cell in basal levels of glucose. At glucose, the islet -cells sharply transferred from inactive to active states, and glucose acutely elevated . In fact, a phase transition between globally active behavior and global quiescence of -cells was observed, as cellular excitability approached a critical threshold. For high levels of glucose , the slope of elevation significantly decreased, and intracellular activity comparatively achieved saturation. Indeed, in the specific range , the -cell dynamics and concentrations were highly sensitive to increasing glucose level and then was only slightly altered, leading to near-independent insulin behavior of the glucose gradient. Also, at high levels of glucose, we observed bimodal behavior in cytosolic calcium concentrations, similar to Fig. 3, caused by varying values of and very different oscillation patterns. The critical point did not depend on the value of gap junction conductance, whereas the strength of coupling could affect the behaviour above the threshold point such that in small values of , intracellular calcium concentration was saturated at higher glucose levels. Fig. 4B and C showed that interestingly, change in glucose levels had no significant effect on global synchronization of electrical activity across the islet.

Figure 4

(A–C) The interplay between gap junction coupling strength and glucose concentration in the mean response, and in the synchronization of and membrane potential oscillations across the intact islet. (D–F) As in panels (A–C), but across the islet lacking -cells. Note that after loss of -cell mass the level of intracellular is more slowly saturated, and the islet synchronization is more sensitive to increasing glucose level. Data represent the mean ± S.E.M. over 10 modeled islets.

Additionally, we eliminated the islet -cells, and then analyzed the role of and glucose levels on activity (Fig. 4D–F). As above, intracellular concentrations of single -cell and the release of insulin exhibited a steep sigmoidal secretory response to increasing glucose levels. The transition between quiet state and avalanche occurred at a position equivalent to glucose, similar to intact islet, while at glucose, the slope of elevation was less in islet lacking -cells compared with the intact islet, i.e., the level of intracellular was more slowly saturated. During high glucose concentrations, in islets lacking -cells was relatively unchanged from the intact islet, although in small values of there were noticeable differences (Fig. 4D). Furthermore, when the level of glucose concentrations was elevated, islet synchrony of electrical patterns had complex behavior. In fact, the synchronization of islet lacking -cells was less than intact islet and more sensitive to increasing glucose level (Fig. 4E,F).

Kv-channels increase levels in the absence of gap junction coupling

We predict that Kv-channels modify the electrical patterns of -cells lacking gap junctions and that changes in -cell burst behavior cause quantitatively alterations in the intracellular concentration. To test this prediction, we first considered the simulated behavior of interconnected -cells, which had Kv-conductances picked from a normal distribution (Fig. 5A), and then determined how different values affected the slow electrical burst patterns. When the modeled -cells were isolated from the islet, not only the stochastic and heterogeneous behavior occurred but the shape of the oscillatory pattern of the membrane potential and changed, which for the amplitude of oscillations was much lower compared to the coupled scenario (Fig. 5B). In continuation, we decreased slightly Kv-channel conductance values of dissociated -cells, and compared activity patterns of these cells with the coupled -cell behavior. Surprisingly, the electrical patterns, specifically activity and the metabolic oscillations resulting from a variation of conductance nicely resembled the global islet behavior mediated through gap junctional communications and synchronizing dynamics, as summarized in Fig. 5A and C. These findings suggest that modifying forms of electrical activity in human -cells to enhance insulin levels can arise out of another regulatory mechanism without coupling and synchrony across the islet, i.e., small changes in the expression of Kv-channels lead to unexpected changes in the integrated behavior of all -cells. In sum, according to this view changes in oscillatory behavior of -cells involving different mechanisms play a physiologically important role in the level of and insulin secretion.

Figure 5

Effect of Kv-channels on patterns. (A) Membrane potential, cytosolic calcium, and FBP time series for coupled human -cells with . (B) As in panel (A), but for uncoupled -cells. (C) Model simulation of uncoupled -cells as in panel (B), but with .

Discussion

Type 2 diabetes mellitus (T2DM) is a complex metabolic disorder triggered by a progressive decline in -cell mass and an increase in insulin resistance. In T2D, there are fundamental changes in insulin secretion dynamics, which can lead to the development of insulin insufficiency. The total amount of released insulin is predominantly the product of the morphological mass of -cells and insulin output of each of these cells. Thus, inadequate levels of insulin are importantly the result of deficiency in either -cell mass or function, or both, leading to hyperglycemia and diabetes. Currently, diabetes research shows that the actual cause of developing T2D seems to be strongly correlated with pancreatic -cells; however it has long been assumed that insulin resistance is the major risk component. In patients with type 2 diabetes, apoptotic -cell death can be induced by various etiological factors, such as exposure to chronic hyperglycemia and carbohydrate metabolites[69], oxidative stres[70,71], pro-inflammatory mediators (e.g., cytokines)[72], proinsulin misfolding[73,74], human islet amyloid polypeptide misfolding[75], as well as endoplasmic reticulum stress[76], which these sequential events initiate apoptosis, increase -cell workload and stress, culminate in exhaustion, and finally -cell death[77-81]. Therefore, creating a successful treatment for T2D will need to specifically include targeting insulin resistance, regenerating -cell mass, and restoring appropriate insulin release by recovery and increase of -cell function. Despite these findings, the precise relationship between morphological and functional -cell in progressing impairments in insulin secretion is a thorny issue, and dynamics of -cell dysfunction and -cell loss in human T2D phenotype is still under debate[82-84].

To get a more detailed insight into the distinct role of human -cell mass and function as cause for type 2 diabetes, we made use of multicellular computational approach of interconnected -cells, based on the theoretical model of Riz et al.[60]. To incorporate known particularities, -cell heterogeneity was characterized by picking parameter values from normal distributions, and -cell mass reduction was established in the islet network by removing simulated cells, randomly. We found that contribution of impaired -cell function to insulin inadequacy, correlated with irregular activity dynamics, could be deduced from -cell mass deficit in type 2 diabetes (Fig. 1). Besides alterations in the summed activity after removal of insulin secreting cells (Fig. 1A), surprisingly the human islet showed decreased levels in the average of an active -cell (Fig. 1B), caused by reduced concentration uptake and asynchronous oscillations (Fig. 1E (I and II)). Our results indicate that in addition to disruption in -cell population electrical activity, changes in the form of oscillatory patterns of membrane potential and intracellular concentration evoke functional deteriorations in -cells, such as insulin secretory dysfunction, which are caused by changes in -cell mass. However, determining the specific role of -cell mass on -cell function is beyond the scope of this manuscript. Therefore, future studies will be needed to precisely examine whether a decrease in -cell mass primarily is the leading factor in -cell dysfunction.

Pulsatile insulin secretion is tightly associated with the synchronous oscillations in electrical activity, particular oscillatory dynamics of within the islet. Insulin pulsatility leads to complex rhythms in blood glucose concentration[85,86], enhances hepatic insulin action and post-receptor signaling[87], protects against insulin resistance and shows a greater efficiency than continuous levels of secretion[43]. Previous experimental investigations have revealed that a decline in first-phase insulin response and an impairment in regular rhythmic secretion patterns can be found at the onset of T2D and thereafter[43,48], which occurred in deficit of -cell mass[56-58]. Our computational modeling confirms these experimental observations, as of -cells are lost, the synchronous membrane potential and excursions, underlying the peak amplitude of early-phase insulin release and the coordination of insulin pulsatility fashion, began to be disrupted within the islet cellular network (Fig. 1C,D).

Cx36 gap junctions have potential roles in dynamics and physiological function of the islet[21,26], and a reduction in -cell coupling has been suggested to occur in type 2 diabetes[49,51,52,88]. Inter--cellular coupling has also been implicated in protecting -cells against a variety of cytotoxic factors, regulating -cell differentiation and maturation, and supporting islet development and fitness[89-91]. Chronic hyperglycemia characterizing the onset of diabetes leads to impaired gap junctional communication[52], which makes the islets more sensitive to -cell death[28,53]. Likewise, it has been demonstrated that a number of chronic insults, including glucotoxicity[92], lipotoxicity[93,94], and pro-inflammatory cytokines and oxidative stress[95] target the expression of Cx36 transcript and inhibit gap junction functionality. Therefore, human -cell islets exposed to such insults exhibit a lack in intra-islet synchrony of oscillations, a suppression and limited propagation of calcium waves, a disruption in plasma insulin pulsatility, and glucose intolerance[96-98]. A significant decrease in the peak level of first-phase insulin release and loss of pulsatile second-phase secretion is mainly the secretory defects observed in human T2D, resulting in disrupted glucose homeostasis[48]. For the most part, the peak elevation of first-phase secretion and the second-phase pulses are dependent on the coordinated pulsatility of individual islets and eventually abolish as the development of diabetes. Furthermore, cell-cell interactions via gap junctional channels are a prerequisite for the oscillatory patterns of electrical activity within the islet and regulate the dynamics of insulin secretion.

Our analysis regarding the interplay between loss of gap junction coupling and impairment of synchronous electrical activity patterns in inadequate levels of plasma insulin and progression of T2D suggest a essentially important impression for intercellular connections. In fact, firstly the level of concentration immediately decreased following reduction in -cell coupling, however the oscillations of -cell membrane potential and synchronized (Fig. 2). Secondly, the -cell population in islets showed poorly coordinated behavior, especially dynamics, after -cell coupling loss, although was harmonize as (Fig. 2 B and C). The data discussed above strongly supports that the role of gap junctional coupling in affecting the cytosolic calcium concentration and the amount of secreted insulin is more subtle and fundamental than synchronous oscillations of -cell activity.

Available studies suggest that mouse and human -cells show different electrical dynamics; mouse -cell population display islet-wide synchrony in response to glucose, whereas human islet synchrony of oscillations is constrained to localized subpopulations[66,99]. These differences probably relate to differences in mouse and human islet architectures: mouse islets have a large, highly connected -cell core, whereas human islets are composed of distinct clusters of gap junction coupled -cells[66,67,100]. Noteworthy, in the case of intra-islet synchronization, our results seem to be in contrast to less coordinated behavior in human islets. These findings highlight that -cells in human islets occur in distinct clusters separated by other cell types, notably -cells and vascular cells[66].

Using computer simulations, we proposed that human -cells exhibit great different behavior in dynamics caused by a substantial decrease in islet gap junction coupling, i.e., lack of -cell connections sped up reducing level in single -cell (Fig. 2A) and lost inter--cell synchronization (Fig. 2B,C), which is believed to lead to impairment of calcium waves and normal oscillatory insulin secretion[44]. In particular, our results show -cell functional insufficiency, such that there exist specific changes in the oscillation patterns of -cell electrical activity (not shown), which most likely be due to combined increased -cell apoptosis and workload, and finally result in functional exhaustion and persistent hyperglycemia.

Gap junctional coupling between -cells provide the intra-islet synchrony of glycolysis oscillations, which is a prerequisite for pulsatile insulin secretion[101]. Slow oscillations in metabolism of glucose-stimulated -cells coupled to electrical activity patterns by oscillations in ATP production and closure of channels are mediated by the positive feedback on the allosteric enzyme phosphofructokinase (PFK) via its product fructose-1,6-bisphosphate (FBP). A rise in substrate G6P, which is converted to F6P, from glucose consumption leads to FBP production, which increases the autocatalytic activation of the enzyme PFK with an eventual crash in the FBP level due to depletion of substrate G6P[102]. Gap junctional permeability leads to diffusion of G6P among islet -cells, which is considerably smaller than other glycolytic metabolites[39].

In order to investigate the effects of coupling electrically and metabolically on the oscillatory behavior of human pancreatic -cells, we used dynamical mathematical model incorporating glycolytic oscillations. Diffusive coupling of metabolites, specifically at low , made an important contribution to the normal function of -cells, however significance of this diffusion process is likely less compared to electrical conduction in regulating islet activity. Namely, metabolic communications, if did not exist, could have a negative effect on the levels (Fig. 3A) and a desynchronizing effect on the bursting activity as well as dynamics of system coupled electrically through gap junctions (Fig. 3C,D), nevertheless increasing electrical coupling strength further was able to overcome islet asynchronous. Addition of a very small degree of metabolic coupling could contribute to intracellular elevation and robust coordination of membrane potential and oscillations. On the other hand, our results revealed different modes of -cell behavior, strongly depending on the relative strength of electrical coupling. Incorporation of enhanced electrical coupling in addition to gap junctional diffusion persuaded great alterations in the intercellular calcium levels (Fig. 3A), and more importantly imposed the bimodal behavior on the single -cell achieved through widely transforming patterns of membrane potential and calcium oscillations (Fig. 3B).

-cells within islets of Langerhans respond to stimulatory glucose level by insulin secretion. Electrical coupling across the islets, principally via Cx36 gap junction channels, mediates oscillatory dynamics of membrane depolarization and to propagate waves and robust pulsatile insulin release under elevated glucose, as well as efficient suppression of spontaneous elevations under basal glucose[6,29], illustrating a strong link between glucose-stimulated insulin secretion and gap junction function. Upon a glucose gradient, a characteristic sigmoidal secretory response is observed in intact islets, indicating critical behavior that depends on physiological properties of gap junction conductance. Additionally, intact islets exhibit more insulin response to increasing glucose than dispersed -cells and that, at nonstimulatory glucose concentrations, insulin levels from dispersed -cells are significantly higher than from intact islets[28]. It is therefore important to consider cellular communication for regulating insulin secretory dynamics and ultimately glucose homeostasis. The disruption of gap junctional coupling results in reduced first-phase amplitude of insulin secretion, loss of coordinated oscillations leading to lack of pulsatile second-phase insulin release, and disrupted glucose homeostasis[28,29,53], similar to defects seen in human patients with type 2 diabetes[43,48,51]. More striking is the fact that islets lacking gap junctions have statistically normal insulin levels and insulin sensitivity, despite glucose intolerance due to altered dynamics of insulin secretion.

We simulated -cell behavior under varying the glucose concentration periodically from 0 to 14 mM to determine the relative role of glucose levels and gap junction activity in shaping the glucose-stimulated of an active -cell and its ability to insulin secretion. In clusters with -cell loss, responses after glucose stimulation were characterized by a sharp transition phase between quiescent and active behavior, and plateau phase that followed, similar to what was measured in intact islets (Fig. 4A,D). Most importantly, our results showed that gap junction coupling strength did not significantly impact the level of intracellular concentration at basal glucose and the position of activation threshold, while the plateau fraction of elevation was gap junctional dependent, later saturating at low coupling conductances. In addition, the synchrony of electrical dynamics across the islet seemed to be almost independent on glucose levels, in case of loss (Fig. 4B,C) and loss (Fig. 4E,F) with increasing the coupling strength.

In the islet, the shape of electrical behavior is highly dependent on the biophysical characterizations of ion channels expressed in a single -cell. Variability in the gating dynamics of specific channels between -cells leads to generate variable patterns of membrane potential and oscillations. Pedersen has demonstrated the capability of Kv-channels to change spiking behaviors to bursting patterns in human -cells[59]. Riz et al.[103] and Montefusco et al.[104] also investigated the contribution of channels in shaping -cell electrical activity and controlling insulin secretion.

Our simulation data revealed that human -cells lacking intercellular coupling exhibit similar electrical patterns to coupled -cells within the islet by smoothly changing the conductance Kv-channel gating, which significantly affected the level of intracellular concentrations. The multicellular behavior of the islet was analyzed, based on the absence of gap junctional connections, to quantitatively describe changes in -cell dynamics after small changes in the expression of delayed rectifying potassium channels. When the Kv-channel conductance was slowly reduced, the shape of electrical activity and oscillatory introduced by low amplitude excursions in Fig. 5B modified, which were similarly observed before disruptions to gap junction coupling in intact islet (Fig. 5A,C). It appears that the spatiotemporal organization of response are likely governed by two different mechanisms, characterized by the introduction of gap junction coupling and synchronizing dynamics, or Kv-channel properties, affecting the burst behavior of -cells and quantitatively the intercellular calcium events. Additionally, these data demonstrated that the form of membrane potential oscillations, correlated with concentration oscillations, is a necessary factor in -cell calcium elevation, in addition to inter--cellular communications and islet synchrony. As changes in the gating of -channels can yield an excess of large events in the patterns and activity of , and the pulse mass of insulin secretion, it will be necessary to uncover the underlying mechanisms of normal Kv-channel function for potential diabetes therapies, however, the exact reason for this remains still unclear.

Conclusion

Our knowledge about the differential contribution of human -cell mass and function in hyperglycemia development and T2D pathogenesis can provide key information for regenerating -cell mass or preserving -cell function. This study demonstrates that -cell mass reduction is an important factor in -cell dysfunction, impairment in intra-islet synchrony, and changes in the shape of electrical bursting, which cause changes in insulin secretion dynamics and insulin levels. The role of gap junction-mediated electrical coupling in affecting the behavior of intracellular dynamics is more significant compared with both metabolic coupling and synchronous oscillations of islet activity. Our results reveal that in human -cells lacking gap junctions modifying electrical patterns to enhance levels and the amount of secreted insulin can arise from changes in the expression of Kv-channels, pointing towards a prominent role of Kv-channels in T2D development and therapy.

Supplementary Information.