Glucose modulates [Ca2+]i oscillations in pancreatic islets via ionic and glycolytic mechanisms

Pancreatic islets of Langerhans display complex intracellular calcium changes in response to glucose that include fast (seconds), slow ( approximately 5 min), and mixed fast/slow oscillations; the slow and mixed oscillations are likely responsible for the pulses of plasma insulin observed in vivo. To better understand the mechanisms underlying these diverse patterns, we systematically analyzed the effects of glucose on period, amplitude, and plateau fraction (the fraction of time spent in the active phase) of the various regimes of calcium oscillations. We found that in both fast and slow islets, increasing glucose had limited effects on amplitude and period, but increased plateau fraction. In some islets, however, glucose caused a major shift in the amplitude and period of oscillations, which we attribute to a conversion between ionic and glycolytic modes (i.e., regime change). Raising glucose increased the plateau fraction equally in fast, slow, and regime-changing islets. A mathematical model of the pancreatic islet consisting of an ionic subsystem interacting with a slower metabolic oscillatory subsystem can account for these complex islet calcium oscillations by modifying the relative contributions of oscillatory metabolism and oscillatory ionic mechanisms to electrical activity, with coupling occurring via K(ATP) channels.     

Intra- and inter-islet synchronization of metabolically driven insulin secretion

Insulin secretion from pancreatic beta-cells is pulsatile with a period of 5-10 min and is believed to be responsible for plasma insulin oscillations with similar frequency. To observe an overall oscillatory insulin profile it is necessary that the insulin secretion from individual beta-cells is synchronized within islets, and that the population of islets is also synchronized. We have recently developed a model in which pulsatile insulin secretion is produced as a result of calcium-driven electrical oscillations in combination with oscillations in glycolysis. We use this model to investigate possible mechanisms for intra-islet and inter-islet synchronization. We show that electrical coupling is sufficient to synchronize both electrical bursting activity and metabolic oscillations. We also demonstrate that islets can synchronize by mutually entraining each other by their effects on a simple model "liver," which responds to the level of insulin secretion by adjusting the blood glucose concentration in an appropriate way. Since all islets are exposed to the blood, the distributed islet-liver system can synchronize the individual islet insulin oscillations. Thus, we demonstrate how intra-islet and inter-islet synchronization of insulin oscillations may be achieved.     

Gap junction coupling and calcium waves in the pancreatic islet

The pancreatic islet is a highly coupled, multicellular system that exhibits complex spatiotemporal electrical activity in response to elevated glucose levels. The emergent properties of islets, which differ from those arising in isolated islet cells, are believed to arise in part by gap junctional coupling, but the mechanisms through which this coupling occurs are poorly understood. To uncover these mechanisms, we have used both high-speed imaging and theoretical modeling of the electrical activity in pancreatic islets under a reduction in the gap junction mediated electrical coupling. Utilizing islets from a gap junction protein connexin 36 knockout mouse model together with chemical inhibitors, we can modulate the electrical coupling in the islet in a precise manner and quantify this modulation by electrophysiology measurements. We find that after a reduction in electrical coupling, calcium waves are slowed as well as disrupted, and the number of cells showing synchronous calcium oscillations is reduced. This behavior can be reproduced by computational modeling of a heterogeneous population of β-cells with heterogeneous levels of electrical coupling. The resulting quantitative agreement between the data and analytical models of islet connectivity, using only a single free parameter, reveals the mechanistic underpinnings of the multicellular behavior of the islet.     

Ca2+ dynamics in a population of smooth muscle cells: modeling the recruitment and synchronization

Many experimental studies have shown that arterial smooth muscle cells respond with cytosolic calcium rises to vasoconstrictor stimulation. A low vasoconstrictor concentration gives rise to asynchronous spikes in the calcium concentration in a few cells (asynchronous flashing). With a greater vasoconstrictor concentration, the number of smooth muscle cells responding in this way increases (recruitment) and calcium oscillations may appear. These oscillations may eventually synchronize and generate arterial contraction and vasomotion. We show that these phenomena of recruitment and synchronization naturally emerge from a model of a population of smooth muscle cells coupled through their gap junctions. The effects of electrical, calcium, and inositol 1,4,5-trisphosphate coupling are studied. A weak calcium coupling is crucial to obtain a synchronization of calcium oscillations and the minimal required calcium permeability is deduced. Moreover, we note that an electrical coupling can generate oscillations, but also has a desynchronizing effect. Inositol 1,4,5-trisphosphate diffusion does not play an important role to achieve synchronization. Our model is validated by published in vitro experiments obtained on rat mesenteric arterial segments.     

Model of intercellular calcium oscillations in hepatocytes: synchronization of heterogeneous cells.

Hepatocytes respond with repetitive cytosolic calcium spikes to stimulation by vasopressin and noradrenalin. In the intact liver, calcium oscillations occur in a synchronized fashion as periodic waves across whole liver lobules, but the mechanism of intercellular coupling remains unclear. Recently, it has been shown that individual hepatocytes can have very different intrinsic oscillation frequencies but become phase-locked when coupled by gap junctions. We investigate the gap junction hypothesis for intercellular synchronization by means of a mathematical model. It is shown that junctional calcium fluxes are effective in synchronizing calcium oscillations in coupled hepatocytes. An experimentally testable estimate is given for the junctional coupling coefficient required; it mainly depends on the degree of heterogeneity between cells. Intercellular synchronization by junctional calcium diffusion may occur also in other cell types exhibiting calcium-activated calcium release through InsP(3) receptors, if the gap junctional coupling is strong enough and the InsP(3) receptors are sufficiently sensitized by InsP(3).     

Long lasting synchronization of calcium oscillations by cholinergic stimulation in isolated pancreatic islets

Individual mouse pancreatic islets exhibit oscillations in [Ca²⁺]_i and insulin secretion in response to glucose in vitro, but how the oscillations of a million islets are coordinated within the human pancreas in vivo is unclear. Islet to islet synchronization is necessary, however, for the pancreas to produce regular pulses of insulin. To determine whether neurohormone release within the pancreas might play a role in coordinating islet activity, [Ca²⁺]_i changes in 4-6 isolated mouse islets were simultaneously monitored before and after a transient pulse of a putative synchronizing agent. The degree of synchronicity was quantified using a novel analytical approach that yields a parameter that we call the “Synchronization Index”. Individual islets exhibited [Ca²⁺]_i oscillations with periods of 3-6 min, but were not synchronized under control conditions. However, raising islet [Ca²⁺]_i with a brief application of the cholinergic agonist carbachol (25 μM) or elevated KCl in glucose-containing saline rapidly synchronized islet [Ca²⁺]_i oscillations for ≥30 min, long after the synchronizing agent was removed. In contrast, the adrenergic agonists clonidine or norepinephrine, and the K_ATP channel inhibitor tolbutamide, failed to synchronize islets. Partial synchronization was observed, however, with the K_ATP channel opener diazoxide. The synchronizing action of carbachol depended on the glucose concentration used, suggesting that glucose metabolism was necessary for synchronization to occur. To understand how transiently perturbing islet [Ca²⁺]_i produced sustained synchronization, we used a mathematical model of islet oscillations in which complex oscillatory behavior results from the interaction between a fast electrical subsystem and a slower metabolic oscillator. Transient synchronization simulated by the model was mediated by resetting of the islet oscillators to a similar initial phase followed by transient “ringing” behavior, during which the model islets oscillated with a similar frequency. These results suggest that neurohormone release from intrapancreatic neurons could help synchronize islets in situ. Defects in this coordinating mechanism could contribute to the disrupted insulin secretion observed in Type 2 diabetes.     

Calcium oscillations in a triplet of pancreatic acinar cells

We use a mathematical model of calcium dynamics in pancreatic acinar cells to investigate calcium oscillations in a ring of three coupled cells. A connected group of cells is modeled in two different ways: 1), as coupled point oscillators, each oscillator being described by a spatially homogeneous model; and 2), as spatially distributed cells coupled along their common boundaries by gap-junctional diffusion of inositol trisphosphate and/or calcium. We show that, although the point-oscillator model gives a reasonably accurate general picture, the behavior of the spatially distributed cells cannot always be predicted from the simpler analysis; spatially distributed diffusion and cell geometry both play important roles in determining behavior. In particular, oscillations in which two cells are in synchrony, with the third phase-locked but not synchronous, appears to be more dominant in the spatially distributed model than in the point-oscillator model. In both types of model, intercellular coupling leads to a variety of synchronous, phase-locked, or asynchronous behaviors. For some parameter values there are multiple, simultaneous stable types of oscillation. We predict 1), that intercellular calcium diffusion is necessary and sufficient to coordinate the responses in neighboring cells; 2), that the function of intercellular inositol trisphosphate diffusion is to smooth out any concentration differences between the cells, thus making it easier for the diffusion of calcium to synchronize the oscillations; 3), that groups of coupled cells will tend to respond in a clumped manner, with groups of synchronized cells, rather than with regular phase-locked periodic intercellular waves; and 4), that enzyme secretion is maximized by the presence of a pacemaker cell in each cluster which drives the other cells at a frequency greater than their intrinsic frequency.     

Diffusion of extracellular K+ can synchronize bursting oscillations in a model islet of Langerhans.

Electrical bursting oscillations of mammalian pancreatic beta-cells are synchronous among cells within an islet. While electrical coupling among cells via gap junctions has been demonstrated, its extent and topology are unclear. The beta-cells also share an extracellular compartment in which oscillations of K+ concentration have been measured (Perez-Armendariz and Atwater, 1985). These oscillations (1-2 mM) are synchronous with the burst pattern, and apparently are caused by the oscillating voltage-dependent membrane currents: Extracellular K+ concentration (Ke) rises during the depolarized active (spiking) phase and falls during the hyperpolarized silent phase. Because raising Ke depolarizes the cell membrane by increasing the potassium reversal potential (VK), any cell in the active phase should recruit nonspiking cells into the active phase. The opposite is predicted for the silent phase. This positive feedback system might couple the cells' electrical activity and synchronize bursting. We have explored this possibility using a theoretical model for bursting of beta-cells (Sherman et al., 1988) and K+ diffusion in the extracellular space of an islet. Computer simulations demonstrate that the bursts synchronize very quickly (within one burst) without gap junctional coupling among the cells. The shape and amplitude of computed Ke oscillations resemble those seen in experiments for certain parameter ranges. The model cells synchronize with exterior cells leading, though incorporating heterogeneous cell properties can allow interior cells to lead. The model islet can also be forced to oscillate at both faster and slower frequencies using periodic pulses of higher K+ in the medium surrounding the islet. Phase plane analysis was used to understand the synchronization mechanism. The results of our model suggest that diffusion of extracellular K+ may contribute to coupling and synchronization of electrical oscillations in beta-cells within an islet.     

Mitochondrial metabolism reveals a functional architecture in intact islets of Langerhans from normal and diabetic Psammomys obesus

The cells within the intact islet of Langerhans function as a metabolic syncytium, secreting insulin in a coordinated and oscillatory manner in response to external fuel. With increased glucose, the oscillatory amplitude is enhanced, leading to the hypothesis that cells within the islet are secreting with greater synchronization. Consequently, non-insulin-dependent diabetes mellitus (NIDDM; type 2 diabetes)-induced irregularities in insulin secretion oscillations may be attributed to decreased intercellular coordination. The purpose of the present study was to determine whether the degree of metabolic coordination within the intact islet was enhanced by increased glucose and compromised by NIDDM. Experiments were performed with isolated islets from normal and diabetic Psammomys obesus. Using confocal microscopy and the mitochondrial potentiometric dye rhodamine 123, we measured mitochondrial membrane potential oscillations in individual cells within intact islets. When mitochondrial membrane potential was averaged from all the cells in a single islet, the resultant waveform demonstrated clear sinusoidal oscillations. Cells within islets were heterogeneous in terms of cellular synchronicity (similarity in phase and period), sinusoidal regularity, and frequency of oscillation. Cells within normal islets oscillated with greater synchronicity compared with cells within diabetic islets. The range of oscillatory frequencies was unchanged by glucose or diabetes. Cells within diabetic (but not normal) islets increased oscillatory regularity in response to glucose. These data support the hypothesis that glucose enhances metabolic coupling in normal islets and that the dampening of oscillatory insulin secretion in NIDDM may result from disrupted metabolic coupling.     

Excitation wave propagation as a possible mechanism for signal transmission in pancreatic islets of Langerhans

In response to glucose application, beta-cells forming pancreatic islets of Langerhans start bursting oscillations of the membrane potential and intracellular calcium concentration, inducing insulin secretion by the cells. Until recently, it has been assumed that the bursting activity of beta-cells in a single islet of Langerhans is synchronized across the whole islet due to coupling between the cells. However, time delays of several seconds in the activity of distant cells are usually observed in the islets of Langerhans, indicating that electrical/calcium wave propagation through the islets can occur. This work presents both experimental and theoretical evidence for wave propagation in the islets of Langerhans. Experiments with Fura-2 fluorescence monitoring of spatiotemporal calcium dynamics in the islets have clearly shown such wave propagation. Furthermore, numerical simulations of the model describing a cluster of electrically coupled beta-cells have supported our view that the experimentally observed calcium waves are due to electric pulses propagating through the cluster. This point of view is also supported by independent experimental results. Based on the model equations, an approximate analytical expression for the wave velocity is introduced, indicating which parameters can alter the velocity. We point to the possible role of the observed waves as signals controlling the insulin secretion inside the islets of Langerhans, in particular, in the regions that cannot be reached by any external stimuli such as high glucose concentration outside the islets.     

Anti-phase calcium oscillations in astrocytes via inositol (1, 4, 5)-trisphosphate regeneration

Observations in cultured mouse astrocytes suggest anti-phase synchronization of cytosolic calcium concentrations in nearest neighbor cells that are coupled through gap junctions. A mathematical model is used to investigate physiologic conditions under which diffusion of the second messenger inositol (1, 4, 5)-trisphosphate (IP(3)) through gap junctions can facilitate synchronized anti-phase Ca(2+) oscillations. Our model predicts anti-phase oscillations in both cytosolic calcium and IP(3) concentrations if (a) the gap junction permeability is within a window of values and (b) IP(3) is regenerated in the astrocytes via, e.g. phospholipase C(delta). This result sheds new light on the current dispute on the mechanism of intercellular calcium signaling. It provides indirect evidence for a partially regenerative mechanism as the model excludes anti-phase synchrony in the absence of IP(3) regeneration.     

Size distribution of mouse Langerhans islets

Pancreatic beta-cells are clustered in islets of Langerhans, which are typically a few hundred micrometers in a variety of mammals. In this study, we propose a theoretical model for the growth of pancreatic islets and derive the islet size distribution, based on two recent observations: First, the neogenesis of new islets becomes negligible after some developmental stage. Second, islets grow via a random process, where any cell in an islet proliferates with the same rate regardless of the present size of the islet. Our model predicts either log-normal or Weibull distributions of the islet sizes, depending on whether cells in an islet proliferate coherently or independently. To confirm this, we also measure the islet size by selectively staining islets, which are exposed from exocrine tissues in mice after enzymatic treatment. Indeed revealed are skewed distributions with the peak size of approximately 100 cells, which fit well to the theoretically derived ones. Interestingly, most islets turned out to be bigger than the expected minimal size (approximately 10 or so cells) necessary for stable synchronization of beta-cells through electrical gap-junction coupling. The collaborative behavior among cells is known to facilitate synchronized insulin secretion and tends to saturate beyond the critical (saturation) size of approximately 100 cells. We further probe how the islets change as normal mice grow from young (6 weeks) to adult (5 months) stages. It is found that islets may not grow too large to maintain appropriate ratios between cells of different types. Our results implicate that growing of mouse islets may be regulated by several physical constraints such as the minimal size required for stable cell-to-cell coupling and the upper limit to keep the ratios between cell types. Within the lower and upper limits the observed size distributions of islets can be faithfully regenerated by assuming random and uncoordinated proliferation of each beta-cell at appropriate rates.     

A model for glucose-induced wave propagation in pancreatic islets of Langerhans

A reaction-diffusion type model is constructed, describing the spatio-temporal dynamics of the basic intracellular variables assumed to be involved in the initiation of the insulin secretion process by beta -cells in the pancreatic islets of Langerhans. The model includes equations for the electric membrane potential of the cells, with respective kinetics for ionic currents, for concentrations of both free and stored intracellular calcium, and for the intra- and extracellular concentrations of glucose. An empirical expression connecting the equation for the intracellular glucose concentration to the electrical equation is introduced. The model reproduces the events observed in experiments in vitro upon external glucose application to the islets of Langerhans, such as usual bursting oscillations of the membrane potential and corresponding oscillations of the intracellular calcium concentration. It also allows simulation of electric wave propagation through the islet, initiated by the spatial gradient of glucose concentration within the islet. The gradient emerges due to glucose diffusing into the islets from the external medium, being high at the edges. The latter results show that glucose diffusion presents a means for wave initiation in the islets, which supports our previous assumption (Aslanidi et al., 2001).     

Calcium and glycolysis mediate multiple bursting modes in pancreatic islets

Pancreatic islets of Langerhans produce bursts of electrical activity when exposed to stimulatory glucose levels. These bursts often have a regular repeating pattern, with a period of 10-60 s. In some cases, however, the bursts are episodic, clustered into bursts of bursts, which we call compound bursting. Consistent with this are recordings of free Ca2+ concentration, oxygen consumption, mitochondrial membrane potential, and intraislet glucose levels that exhibit very slow oscillations, with faster oscillations superimposed. We describe a new mathematical model of the pancreatic beta-cell that can account for these multimodal patterns. The model includes the feedback of cytosolic Ca2+ onto ion channels that can account for bursting, and a metabolic subsystem that is capable of producing slow oscillations driven by oscillations in glycolysis. This slow rhythm is responsible for the slow mode of compound bursting in the model. We also show that it is possible for glycolytic oscillations alone to drive a very slow form of bursting, which we call "glycolytic bursting." Finally, the model predicts that there is bistability between stationary and oscillatory glycolysis for a range of parameter values. We provide experimental support for this model prediction. Overall, the model can account for a diversity of islet behaviors described in the literature over the past 20 years.     

Effects of gap junction to Ca(2+) and to IP(3) on the synchronization of intercellular calcium oscillations in hepatocytes

The frequency of free cytosolic calcium concentration ([Ca(2+)]) oscillations elicited by a given agonist concentration differs between individual hepatocytes. However, in multicellular systems of rat hepatocytes and even in the intact liver, [Ca(2+)] oscillations are synchronized and highly coordinated. In this paper, we have investigated theoretically the gap junction permeable to calcium and to IP(3) on intercellular synchronization by means of a mathematical model, respectively. It is shown that gap junction permeable to calcium and to IP(3) are effective on synchronizing calcium oscillations in coupled hepatocytes. Our theoretical results are similar either for the case of Ca(2+) acting as coordinating messenger or for the case of IP(3) as coordinating messenger. There exists an optimal coupling strength for a pair of connected hepatocytes. Appropriate coupling strength and IP(3) level can induce various harmonic locking of intercellular [Ca(2+)] oscillations. Furthermore, a phase diagram in two-dimensional parameter space of the coupling strength and IP(3) level (or the velocity of IP(3) synthesis) has been predicted, in which the synchronization region is similar to Arnol'd tongue.     

Emergent properties of electrically coupled smooth muscle cells

Asynchronous and synchronous calcium oscillations occur in a variety of cells. A well-established pathway for intercellular communication is provided by gap junctions which connect adjacent cells and can mediate electrical and chemical coupling. Several experimental studies report that cells presenting only a transient increase when freshly dispersed may oscillate when they are coupled. Such observations suggest that the role of gap junctions is not only to coordinate calcium oscillations of adjacent cells. Gap junctions may also be important to generate oscillations. Here we illustrate the emergent properties of electrically coupled smooth muscle cells using a model that we recently proposed. A bifurcation analysis in the case of two cells reveals that synchronous and asynchronous calcium oscillations can be induced by electrical coupling. In a larger population of smooth muscle cells, electrical coupling may result in the creation of groups of cells presenting synchronous calcium oscillations. The elements of one group may be distant from each other. Moreover, our results highlight a general mechanism by which gap junctional electrical coupling can give rise to out of phase calcium oscillations in smooth muscle cells that are non-oscillating when uncoupled. All these observations remain true in the case of non-identical cells, except that the solution corresponding to synchronous calcium oscillations disappears and that the formation of groups is sensitive to the degree of heterogeneity. The first two authors contributed equally to this work.     

Three roads to islet bursting: emergent oscillations in coupled phantom bursters

Glucose-induced membrane potential and Ca(2+) oscillations in isolated pancreatic beta-cells occur over a wide range of frequencies, from >6/min (fast) to <1/min (slow). However, cells within intact islets generally oscillate with periods of 10-60 s (medium). The phantom bursting concept addresses how beta-cells can generate such a wide range of frequencies. Here, we explore an updated phantom bursting model to determine how heterogeneity in a single parameter can explain both the broad frequency range observed in single cells and the rarity of medium oscillations. We then incorporate the single-cell model into an islet model with parameter heterogeneity. We show that strongly coupled islets must be composed of predominantly medium oscillating single cells or a mixture of fast and slow cells to robustly produce medium oscillations. Surprisingly, we find that this constraint does not hold for moderate coupling, and that robustly medium oscillating islets can arise from populations of single cells that are essentially all slow or all fast. Thus, with coupled phantom bursters, medium oscillating islets can be constructed out of cells that are either all fast, all slow, or a combination of the two.     

Multiple factors influence calcium synchronization in arterial vasomotion

The intercellular synchronization of spontaneous calcium (Ca(2+)) oscillations in individual smooth muscle cells is a prerequisite for vasomotion. A detailed mathematical model of Ca(2+) dynamics in rat mesenteric arteries shows that a number of synchronizing and desynchronizing pathways may be involved. In particular, Ca(2+)-dependent phospholipase C, the intercellular diffusion of inositol trisphosphate (IP(3), and to a lesser extent Ca(2+)), IP(3) receptors, diacylglycerol-activated nonselective cation channels, and Ca(2+)-activated chloride channels can contribute to synchronization, whereas large-conductance Ca(2+)-activated potassium channels have a desynchronizing effect. Depending on the contractile state and agonist concentrations, different pathways become predominant, and can be revealed by carefully inhibiting the oscillatory component of their total activity. The phase shift between the Ca(2+) and membrane potential oscillations can change, and thus electrical coupling through gap junctions can mediate either synchronization or desynchronization. The effect of the endothelium is highly variable because it can simultaneously enhance the intercellular coupling and affect multiple smooth muscle cell components. Here, we outline a system of increased complexity and propose potential synchronization mechanisms that need to be experimentally tested.     

Dynamics of spiking neurons connected by both inhibitory and electrical coupling

We study the dynamics of a pair of intrinsically oscillating leaky integrate-and-fire neurons (identical and noise-free) connected by combinations of electrical and inhibitory coupling. We use the theory of weakly coupled oscillators to examine how synchronization patterns are influenced by cellular properties (intrinsic frequency and the strength of spikes) and coupling parameters (speed of synapses and coupling strengths). We find that, when inhibitory synapses are fast and the electrotonic effect of the suprathreshold portion of the spike is large, increasing the strength of weak electrical coupling promotes synchrony. Conversely, when inhibitory synapses are slow and the electrotonic effect of the suprathreshold portion of the spike is small, increasing the strength of weak electrical coupling promotes antisynchrony (see Fig. 10). Furthermore, our results indicate that, given a fixed total coupling strength, either electrical coupling alone or inhibition alone is better at enhancing neural synchrony than a combination of electrical and inhibitory coupling. We also show that these results extend to moderate coupling strengths.