Next Article in Journal
Neutral and Pectic Heteropolysaccharides Isolated from Opuntia joconostle Mucilage: Composition, Molecular Dimensions and Prebiotic Potential
Next Article in Special Issue
Support Provided by Caregivers for Community-Dwelling Diabetic Hispanic Adults with Intellectual Disabilities and Comorbid Conditions
Previous Article in Journal
Is Hormone Replacement Therapy a Risk Factor or a Therapeutic Option for Alzheimer’s Disease?
Previous Article in Special Issue
Role of Impaired Glycolysis in Perturbations of Amino Acid Metabolism in Diabetes Mellitus
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

The Dynamics of Calcium Signaling in Beta Cells—A Discussion on the Comparison of Experimental and Modelling Data

1
Institute of Dynamics and Vibrations, Technische Universität Braunschweig, D-38106 Braunschweig, Germany
2
Institute of Pharmacology, Toxicology and Clinical Pharmacy, Technische Universität Braunschweig, D-38106 Braunschweig, Germany
*
Authors to whom correspondence should be addressed.
Int. J. Mol. Sci. 2023, 24(4), 3206; https://doi.org/10.3390/ijms24043206
Submission received: 23 December 2022 / Revised: 28 January 2023 / Accepted: 30 January 2023 / Published: 6 February 2023
(This article belongs to the Special Issue Diabetes Mellitus (DM) - Endocrine and Metabolic Disorders)

Abstract

:
The stimulus–secretion coupling of the pancreatic beta cell is particularly complex, as it integrates the availability of glucose and other nutrients with the neuronal and hormonal input to generate rates of insulin secretion that are appropriate for the entire organism. It is beyond dispute however, that the cytosolic Ca2+ concentration plays a particularly prominent role in this process, as it not only triggers the fusion of insulin granules with the plasma membrane, but also regulates the metabolism of nutrient secretagogues and affects the function of ion channels and transporters. In order to obtain a better understanding of the interdependence of these processes and, ultimately, of the entire beta cell as a working system, models have been developed based on a set of nonlinear ordinary differential equations, and were tested and parametrized on a limited set of experiments. In the present investigation, we have used a recently published version of the beta cell model to test its ability to describe further measurements from our own experimentation and from the literature. The sensitivity of the parameters is quantified and discussed; furthermore, the possible influence of the measuring technique is taken into account. The model proved to be powerful in correctly describing the depolarization pattern in response to glucose and the reaction of the cytosolic Ca2+ concentration to stepwise increases of the extracellular K+ concentration. Additionally, the membrane potential during a KATP channel block combined with a high extracellular K+ concentration could be reproduced. In some cases, however, a slight change of a single parameter led to an abrupt change in the cellular response, such as the generation of a Ca2+ oscillation with high amplitude and high frequency. This raises the question as to whether the beta cell may be a partially unstable system or whether further developments in modeling are needed to achieve a generally valid description of the stimulus–secretion coupling of the beta cell.

Graphical Abstract

1. Introduction

1.1. Physiological Background

The pancreatic beta cell transduces the increased availability of glucose and a few other carbohydrates, amino acids and keto acids into stimulated rates of secretion. While a number of neuronal and hormonal signals are able to modify the secretory response, it is the metabolic breakdown of these nutrients which produces the indispensable signals for stimulated secretion [1,2]. The KATP channel at the plasma membrane was identified as the link between the energy metabolism and the electrical activity of the beta cell, since the increased ATP/ADP ratio in the cytosol leads to a decreased potassium efflux, thus enabling the depolarization of the plasma membrane [3,4]. The view that the depolarization-induced Ca2+ influx via voltage-dependent Ca2+ channels is the final common pathway of stimulus–secretion coupling, and that the pattern of the cytosolic Ca2+ concentration determines the pattern of insulin secretion [5], was qualified by the identification of the amplifying pathway, which affects the secretion by signals circumventing plasma membrane depolarization [6,7].
Nevertheless, it is beyond dispute that the cytosolic Ca2+ concentration is of central importance not only for the onset of stimulated insulin secretion, but also for its entire duration. It regulates not only the fusion of insulin granules with the plasma membrane, but also the granule itinerary [8]. Additionally, the cytosolic Ca2+ concentration mediates multiple feedback loops by influencing the mitochondria, the endoplasmic reticulum and Ca2+-dependent K+ channels at the plasma membrane [9,10,11]. In view of the resultant complexity of the interrelations, modelling is necessary to achieve a quantitative description of the dynamics of the Ca2+ concentration in the cytosol and in subcellular compartments, together with changes in the membrane potential and in the glucose metabolism. This paper aims to compare published models to describe these dynamics, with actual measurements from the literature as well as from our own work. Because they are essential for the discussion, the established procedures for measuring and modelling these processes are subsequently presented in detail.

1.2. Methods for the Determination of the Plasma Membrane Potential

While the measurement of the plasma membrane potential has also been performed with fluorescent indicators such as oxonol or carbocyanine dyes [12], more quantitative observations of electrical activity are obtained with electrophysiological methods (see Figure A1 in Appendix A). In the original version, a sharp microelectrode had to be inserted into the cell under study, and the voltage was measured with reference to another electrode in the surrounding solution. This voltage is then fed into a comparator and compared with a previously set target voltage (i.e., the “command potential”). The current required to maintain the desired voltage is injected via a second inserted electrode, and the course of the voltage between the cell and the surrounding medium is calculated from this measurement. In contrast to the voltage clamp, a current can also be preset to flow to the cell. This is achieved by dropping a predetermined potential across a known resistance. Since no voltage is specified, a current clamp measurement results [13].
The microelectrode technique reports the entirety of events at the plasma membrane [14,15], which is not a disadvantage for the measurement of the membrane potential, but precludes the analysis of ion currents at the single-channel level. This disadvantage was overcome by the development of the patch clamp technique [16]. In this method, the cell membrane is not perforated, but a smooth pipette is placed on the cell and a small separated area (or “patch”) is created by slight suction applied to the pipette. The close contact (or “gigaseal”) isolates the membrane patch from the rest of the cell and enables single-channel measurements [16]. The electrode for the measurement is located in the pipette and serves both purposes of voltage measurement with reference to ground and current injection.
Further developments have widened the application of this technique beyond single-channel measurements. If stronger suction is applied to the attached cell, the membrane patch is ruptured and a direct access to the cell interior is generated. This configuration is referred to as “whole cell recording” [17], which corresponds to the measurement with an inserted sharp electrode but has a much larger connection between the pipette solution and the cytosol. To preclude the dialysis of the cytosol, the “perforated patch” configuration was devised, wherein the electrical contact between the pipette solution and the cell interior is created by pore-forming antibiotics such as nystatin [18].

1.3. Methods to Study the Cytosolic Calcium Concentration

The early measurements of the intracellular Ca2+ concentration were performed with electrodes in glass pipettes that contained a Ca2+ ionophore solution in the tip [19]. This required the impaling of the cells under study and was a manually demanding technique. The decisive methodological advance was achieved by the synthesis of fluorescent indicators which change the fluorescence emission, depending on the coordinative binding of Ca2+. In order to load the cells with these indicators, the cells must be membrane permeant; this was achieved by masking the Ca2+ binding groups with acetoxymethyl esters. The ester bond can be cleaved by esterases in the cytosol, and the now polar indicator compound is trapped within the cell [20]. While the accumulation of the indicator is needed for a good signal-to-noise ratio, an unavoidable drawback is that the binding of Ca2+ can decrease its cytosolic concentration and thus affect the signal to be reported when the concentration of the indicator is very high.
More recently, the use of genetically engineered fluorescent Ca2+-binding peptides (GECIs) has gained importance, since they can report localized Ca2+ changes when fused to a suitable cellular protein [21]. Thus far, the modeling of beta cell function is based on observations made with conventional small-molecule indicators, which report the Ca2+ concentration in the bulk cytosol. The Ca2+ binding affinity plays a major role in the selection of the indicator. Indicators with a high affinity such as Fura 2 [20], arguably the most used conventional Ca2+ indicator, tend to underestimate peak values occurring during strong Ca2+ influx, and vice versa; low affinity indicators will miss more subtle changes. Finally, the optical properties influence the selection. The Fluo series of indicators [22,23] are best excited at about 490 nm, which makes them useful for laser excitation. In contrast to the ratiometric Fura 2, the Fluo indicators are excited by a single wavelength, which may require the normalization of the data when a mean value of several measurements is calculated.

1.4. Evolution of Models

Mathematical modelling of stimulus–secretion coupling in the pancreatic beta cell began in the early 1980s [24]. However, at that time, the knowledge was too sketchy to permit sufficiently complex models. Only with the further development of electrophysiological and microfluorometric techniques could the oscillation of the membrane potential and of the cytosolic Ca2+ concentration during glucose stimulation be adequately characterized [25,26]. The generation of these oscillations and their relevance for the kinetics of insulin secretion and for glucose tolerance have remained a matter of debate ever since.
Both, the metabolism of glucose and the feedback between cytosolic Ca2+ levels and ion channel activity are likely to be relevant [27,28]. For a complete picture, the interactions within the islet have to be taken into account [29]. Correspondingly, the development of models progressed by taking more of the above mechanisms into account. The models are based on coupled non-linear ordinary differential equations to describe the multiplicity of oscillations. Special attention was paid to the fact that the dynamics consist of superimposed oscillations with widely separated period lengths.
Focusing initially on the electrical phenomena, the “phantom burster model“ was generated [30] based on four state variables. In addition to the membrane potential and the potassium flux, two abstract variables and their time-related constants were included. However, this did not yet allow the description of compound bursting. The next evolutionary step was the “dual oscillator model” (DOM) [31], which was based on seven state variables. Instead of the abstract state variable of the initial model, the corresponding Ca2+ concentrations were considered, and three equations described the contribution of the metabolism. In particular, the activity of the glycolytic enzyme phosphofructokinase (PFK) was considered a generator of oscillations.
A further layer of complexity resulted from the inclusion of ATP binding to the K+ channels and ATP activation of the Ca2+ pumps of the endoplasmic reticulum (ER) which in turn promotes PFK activity [32]. Here, the activity of PFK was calculated according to the model of [33]. In contrast to previous models, the model presented by [34] contains a direct coupling between the changes of Ca2+ and those of ATP. This relation is established via the regulatory effect of Ca2+ on the mitochondrial enzyme pyruvate dehydrogenase. Due to the closer relationship between the two origins of oscillation, this model was named the “integrated oscillator model” (IOM). This model incorporated the oscillatory nature of glycolysis as measured in pancreatic beta cells [35].

2. Results and Discussion

2.1. Modeling of Membrane Currents

The ion currents through the plasma membrane play a central role in the modelling. In light of this, it will first be analyzed whether fundamental measurement results can be described with the reference model, and how well the associated parameters meet the experimentally obtained quantities. This appeared necessary, since very little information on the exact method of parameter determination or on the basis of the experiment is contained in the pertaining publications [36,37].
In beta cells, the pattern of action potential spiking results mainly from the currents flowing through the voltage-dependent potassium and Ca2+ channels, whereby the transmembrane voltage determines both the open probability of the channels and the driving force of the ion current. To test the descriptive power of the equations and the given set of parameters, the fluxes through potassium and Ca2+ channels in human beta cells published by [38] were chosen (see Figure 1).
The Ca2+ current (left, A, the total current, marked as ■, is to be considered here) is characterized by the fact that at negative voltages the Ca2+ channels are only rarely open ( m   in Table A2, Appendix B is very small), which is why the corresponding current is also small. The probability of the channel being open increases with increasing voltage ( m   becomes larger), but at the same time, the distance to the Ca2+ equilibrium potential V C a decreases, which reduces the driving force and, consequently, the current.
The outward current through the potassium channels (right, B, the current without the addition of Co2+, an inhibitor of Ca2+ influx, is marked as ●) is also influenced by the above factors. Since the potassium equilibrium potential is in the negative range, the current increases with increasing voltage. The shoulder is caused by the current through the large-conductance Ca2+-activated K channels (BK channels), reacting to the opening of voltage-dependent Ca2+ channels.
Figure 2 shows the current–voltage relations for the Ca2+ flux (left) and the K+ flux (right), as calculated according to the associated formulas in Table A2, Appendix B. The course of the Ca2+ current is shown for two different equilibrium Ca2+ voltages V C a . These can take on very different values due to the variability of the cytosolic Ca2+ concentration. If a V C a value of 25 mV is set as in the reference model [36], the formula describes the blue curve, which deviates strongly from the sum of the measured Ca2+ currents shown in Figure 1A. If, on the other hand, the V C a value is set to 55 mV, the resulting orange curve represents the measured values well, both qualitatively and quantitatively.
The right panel in Figure 2 depicts the curves of the individual potassium ion currents depending on the voltage. For simplification, the steady state is assumed here for I K (i.e.,   n = n ), and open channels are assumed for I K ( C a ) and I K ( A T P ) . The violet curve can be compared with the curves in Figure 1 on the right. The increase between −20 mV and 0 mV can also be seen very well here, and its causation by the opening of the I K channels can be clearly identified (blue curve). A remaining difference is a shoulder with a peak at 40 mV, which does not appear, since the Ca2+ influx through voltage-dependent Ca2+ channels and its effect on I K ( C a ) were not included in the calculation.
These comparisons between measurements and simulation confirm the potential of the existing equations to represent complex phenomena. On the other hand, it also becomes clear that moderate changes in the parameters can have marked consequences on the system’s behavior.

2.2. Inflence of Agents

In this chapter, the reaction of the membrane potential to a set of insulinotropic stimuli is examined. In addition to glucose, which acts as a nutrient secretagogue (see above), the purely depolarizing stimuli, KCl and tolbutamide, are considered. Increasing the extracellular K+ concentration raises the potassium equilibrium potential, and tolbutamide reduces the outwardly directed K+ current by closing KATP channels.

2.2.1. Variation of Glucose Concentration

To test how well the model describes the changes of action potential spiking in response to stimulatory glucose concentrations, the measurements published by [39] were chosen. The development of the action potential pattern is depicted in Figure 3. It shows oscillations on different time scales (seconds and hundredths of a second). The high-frequency oscillations represent bursts of action potentials that last a few seconds, separated from each other by electrically silent phases of variable length. At all stimulatory glucose concentrations, the membrane potential oscillates between ca. −45 mV and −10 mV during the action potential bursting, and falls to about −55 mV during the interspersed repolarization phases. The measurements show that raising the glucose concentration to higher values diminishes the distance between the bursting phases until no repolarization occurs at 20 mM glucose.
To compare this measurement with the simulation, it was purposive to introduce a factor A G C to the ADP-term in the equation describing the relation between ATP and ADP (cf. Equation (A7) in Table A1, Appendix B), since the extracellular glucose concentration is not an explicit component of the model proposed by [36].
d A D P d t = A T P A G C e x p [ ( 1 + 2.2 J P D H 0.05 + J P D H ) ( 1 C a 0.35 ) ] A D P τ A
Actually, in the most recent publication by this group [40] a separate parameter G was included in the equations, but the effect of A G C , as used in the present paper, is mathematically comparable. A high value of A G C accelerates the decrease of ADP with increasing extracellular glucose concentrations. With respect to the reference parameter set, the A G C value is set to 1.0 for a glucose concentration of 10 mM, to match the measured curve as closely as possible. Figure 4 shows the calculated courses of the membrane potential and allows a direct comparison with the measurement in Figure 3. The agreement between both is excellent in many respects, and it particularly concerns the following:
  • The quantitative values of the membrane potential (except for the small deviations in the maximum potential during the bursts).
  • The characteristics of the curves (alternation between high-frequency and low-frequency fractions).
  • The length and variability of the phases between the action potentials. In particular, the ability to reproduce this variability with a set of deterministic ordinary differential equations is very remarkable.
  • The number of action potentials in the considered time interval.
Zoom-ins of the first action potential in each case are depicted on right-hand panels in Figure 3 and Figure 4. The authors do not have raw data for the measurement from the literature [39] in Figure 3, but from the representation and time scale, it can be read that the duration of this first action potential is approximately 8 s, during which about 20–30 spikes occur. The distance between the spikes is particularly small at the beginning. This property can be compared with the zoom-in from the simulation in Figure 4. Here, the duration of the first action potential is approximately 6 s, in which there are exactly 20 spikes. From this, it can be concluded from the comparison between the measurement and the simulation that
  • the duration of the action potentials is of a comparable size.
  • the number of spikes in relation to the duration matches very well.
  • The distance between the spikes within the action potential increases over time. This is a well-known feature, see, e.g., [41].
The values for A G C that achieve the best possible agreement at higher glucose concentrations were sought for the model, and gave ( G 15 A G C = 1.25 ,   G 20 A G C = 1.5 ) . Thus, the relation between the extracellular glucose and the AGCvalue is not linear, which is easily recognizable from the fact that a 50% increase ( A G C = 1.5 ) results in a depolarization pattern that fits the experimental observations obtained with 20 mM glucose, but not those with 15 mM glucose. If one now compares simulation and measurement from the point of view of varying glucose supply, the following conclusions can be drawn:
  • In the simulation as well as in the measurement, the time between two bursts of action potentials shortens with increasing A G C (increasing glucose supply).
  • At G 20 or A G C = 1.5 , respectively, there are practically no more resting phases; action potentials are permanently present.
  • In the experimental registrations, the peak value of the action potentials decreases with increasing G and time; in the simulation, this value remains nearly constant (at about 20   mV ).
As a further characteristic, the ratio of ATP/ADP can be specified for these systems at each point in time. The corresponding curve is shown in Figure 5.
This figure shows that the model modification described above couples the ATP/ADP ratio to the increasing extracellular glucose concentration, a relation that has repeatedly been demonstrated experimentally [42,43]. Primarily, this feature affects the mean value. During the electrical activity, a decrease in the ATP/ADP ratio does occur over time (best seen in the zoom-in to the right), but this is orders of magnitude smaller than that due to a decrease in glucose concentration.
In order to test the robustness of the model, the values of one of the parameters were changed. In beta cells, the influx of Ca2+ forms the upstroke of the action potential [44]. At a given open probability of the Ca2+ channels, the flux is proportional to the difference between the actual membrane potential, which is typically in the negative range, and the Ca2+ equilibrium potential ( V C a ), which is in the positive range. Consequently, the larger the value of V C a is, the larger the flux will be. V C a was only moderately changed; it decreased from 25 mV to 20 mV and increased from 25 mV to 30 mV.
In Figure 6, left panel shows the curves of V and C a for the first interval (G10, A G C = 1 ). The blue curve can be directly compared with the left curve in Figure 4. When V C a is reduced to 20 mV, no action potentials occur. The membrane potential and the cytosolic Ca2+ concentration change only marginally over time, and then reach a steady state. In the case of an increased V C a value (30 mV, Figure 6, right panel), the behavior is virtually the opposite; continuous large-amplitude action potentials are accompanied by a high-frequency oscillation of the cytosolic Ca2+ concentration.
This extreme reaction to a moderate change of one model parameter is an indication of how much the informative value of this model depends on the right choice of the model parameters. Since it cannot be assumed that a real beta cell has a similar sensitivity to the Ca2+ equilibrium potential, a further evolution of the model seems necessary to improve its robustness.

2.2.2. Variation of K+ Concentration

Raising the extracellular K+ concentration from a physiological value of about 5 mM to 30 or 40 mM, or even higher, is a widely used experimental procedure to elicit a Ca2+ influx via voltage-dependent Ca2+ channels. The resultant increase of the cytosolic Ca2+ is sufficient to stimulate the secretion of insulin [45]. To test the ability of the model to correctly describe a complex pattern of the cytosolic Ca2+ concentration in response to changes in the extracellular K+ concentration, we chose the measurement published by [39] and depicted in Figure 7. Here, the potassium concentration is increased stepwise from 7.5 mM to 30 mM, interrupted by 2.5   min intervals at the physiological value of 4.8 mM. The authors of the present paper have added the blue scale on the right of the figure to further quantify the results, in order to be able to compare individual quantities with the simulations.
Up to the increase in the potassium concentrations to 12   mM , a nearly continuous (linear over time, see solid blue line) increase of the cytosolic Ca2+ concentration predominates until a self-reinforcing effect occurs due to progressively stronger depolarisation, which causes the Ca2+ concentration to rise rapidly before it decreases again rapidly with a reduction in the potassium concentration. The difference between the temporary maxima (right before the potassium concentration is reduced again) and the temporary minima (immediately before the potassium concentration is increased again), increases with increasing potassium concentrations. These are clearly readable, especially in the last three oscillations (cf. Figure 7), and are listed in Table 1:
For the model, the influence of the changing potassium concentration is taken into account via the potassium equilibrium potential V K . For this purpose, the value of V K is changed abruptly for the respective time intervals, using the Nernst equation as a basis, and assuming a temperature of T = 310   K and an intracellular potassium concentration of K i n t r a = 130   mM :
V K = R T F l n ( K e x t r a K i n t r a ) = 8.314   kg   m 2 s 2   mol   K 310 K 96485.332 C mol l n ( K e x t r a 130   mM ) = 26.7 × 10 3 V   l n ( K e x t r a 130   mM )
This formula generated the curve for V K depending on the extracellular K+ concentration, as shown in Figure 8. The data tips represent the values used for the given K+ concentrations in the model.
With this procedure and the parameters set from the reference [36], the response of the cytosolic Ca2+ concentration and of the plasma membrane potential was simulated (Figure 9). With respect to the Ca2+ values (orange curve in Figure 9), certain similarities within the measurements can be seen, such as the increasing rise of the peaks for increasing K+ concentrations; however, clear differences are also visible, such as the Ca2+ increases at the low K+ concentrations. In particular, the bursts of action potentials (blue curve in Figure 9), occurring even before the K+ concentration was raised, are a cause for concern. Since the potassium channels remain open during K+ depolarization, no action potentials are generated by the beta cells [46]. In the given example, the open state of the KATP channels, which make up the major part of the K+ conductance, was ensured by the presence of diazoxide.
The consideration of this influence was not a major purpose of the reference model [36] and was therefore not taken into account. In the model, the permeability of the KATP channels is essentially determined by the parameter o , compare Table A2 in Appendix B. In the simulation described above, this value fluctuates slightly around a value of about 0.006. In order to adequately incorporate the influence of diazoxide, the parameter was subsequently set to the constant value of o = 0.03 , i.e., about five times as large as the reference value. This means that the KATP channels are significantly more permeable than in the calculated reference configuration. Keeping all other parameters unchanged, the curves shown in Figure 10 are obtained.
This change of parameter results in the absence of action potentials, and the rise of the cytosolic Ca2+ concentration shows greater similarity with the measurement. The quantitative comparison between the simulation results and the measurement is presented in Table 2:
It turns out that this relatively small change in system parameters matches the simulation much better to the measured behavior than the original set of parameters.
The remaining difference between the simulated Ca2+ response and the measured Ca2+ response is largely caused by the lacking linear increase of the minimal values between the K+ depolarizations. This was addressed by decreasing the terms describing the exchange of Ca2+ between the ER and the cytosol and between the mitochondria and the cytosol (see Equation (A4) in Table A1, Appendix B). While a decrease to 50% resulted in a very modest increase, a decrease to 10% generated a clearly visible increase, see Figure 11.
In conclusion, the model as such is very able to reproduce kinetic responses to a mechanistically simple stimulation with considerable robustness.

2.2.3. Effects of Pharmacological Block of KATP Channels by Tolbutamide

In contrast to the depolarization maneuver by high extracellular K+ concentration, the block of KATP channels produces a depolarization of the beta cell with action potential spiking. A virtually complete closure of the channels can be achieved by 500 μM of the sulfonylurea, tolbutamide. Interestingly, the amplitude of the action potentials diminishes to virtual non-existence when the extracellular K+ concentration is raised from 15 mM to 40 mM [47], see Figure 12.
The addition of 15 mM KCl leads to a moderate depolarization of the membrane potential (comparable to the behavior in Section 2.2.2). Only after the addition of tolbutamide and the resulting blockade of the KATP channels does a rapid depolarization of the membrane potential occur, as the depolarizing Ca2+ influx can no longer be balanced by a repolarizing K+ efflux. The additional depolarization opens the potassium channels through which the current I K now leads to potassium efflux, which decreases the membrane potential before Ca2+ influx again predominates and the membrane potential rises again. With the subsequent increase in the extracellular K+ concentration, the K+ equilibrium potential V K is increased, whereby all repolarizing K+ currents are diminished, inhibiting the generation of action potentials.
The simulation accounts for the addition of tolbutamide in the way that, according to Hill’s equation, the parameter g K ( A T P ) is modified, which considers the permeability of the ATP-dependent potassium channels. Specifically, in this model extension, the tolbutamide concentration T determines the tolbutamide-dependent permeability according to Equation (3):
g K ( A T P ) , T o l b = g K ( A T P ) , ( 1 T n T 0.5 n + T n )
with the half-effect concentration of T 0.5 n = 1.90   µ M [48] and the HILL exponent n = 1 , it can be obtained for the case T = 500   µ M :   g K ( A T P ) , T o l b = 2.69 × 10 3 g K ( A T P ) , which corresponds to a very low residual permeability.
With the slight modification of u m = 15   mV and s m = 6   mV (without these there would be no bursts) the behavior depicted in Figure 13a can be generated. It bears considerable resemblance to the experimental measurement shown in Figure 12. There are no bursts in the presence of 15 mM KCl, but the membrane potential is already elevated; however, the value of about −57 mV is still significantly below that of the measurement. With the addition of tolbutamide, the action potentials become visible. With a further increase of the KCl concentration to 22.5 mM, the mean membrane potential continues to rise, but the amplitude of the action potentials increases too, which does not conform with the measurement. With the increase to 30 mM KCl, the mean membrane potential continues to increase, but the amplitude of the oscillation decreases, which is closer to the characteristic of the measurement. Finally, at 40 mM KCl, the bursting ends, as in the measurement.
For these calculations, the equilibrium potentials according to Equation (2) have been used, which exclusively considers potassium ions. While the potassium conductance is by far dominating at the resting state, other conductances through voltage-dependent channels play a role during prolonged depolarization. To account for this influence phenomenologically and to get first ideas on this effect, another simulation was performed, wherein the potential calculated according to Equation (2) was artificially increased by 10 mV. The curves resulting from this modification are shown in Figure 13b. The simulated membrane potential curve now closely resembles the one in the measurement. The membrane potential before the addition of tolbutamide is −45 mV, and the mean membrane potential during action potential spiking agrees better with the measurement. Furthermore, the marked decrease in action potential amplitude is now also evident when the K+ concentration is increased from 15 mM to 22.5 mM.

2.3. Influence of the Measurement Itself

As mentioned in Section 1.3, the cytosolic Ca2+ concentration is usually measured by fluorescent indicators, whereas the plasma membrane potential is best measured by electrophysiological techniques. Thus, it is possible to combine the measurement of the membrane potential with the microfluorimetric measurement of the cytosolic Ca2+ concentration. Due to the fact that the occurring voltages and currents are very small, there are high demands on the measurement system itself, especially to exclude a feedback effect of the measurement on the system behavior.
An indication that there may be a feedback effect of the measurement on the system is shown in Figure 14. In this experiment, the K+ concentration was increased stepwise in the presence of tolbutamide, first to 15 mM, then to 40 mM. In the left panel, the plasma membrane potential was simultaneously measured by current clamping in the “perforated patch” configuration. In the right panel, the same configuration was generated, but no current clamp was applied.
It is obvious that the Fluo 4-fluorescence and thus the cytosolic Ca2+ concentration is significantly smaller with the simultaneous measurement of the membrane potential than that without. If the measurements are feedback free, there should be no difference between the measurements.
In the following, we attempt to derive an explanation for this effect, based on the model presented. For this purpose, a single measurement from the mean values in the left panel of Figure 14 will be used first (Figure 15, left panel).
The calculated curves for this experiment as obtained with the parameters from Table A3 (in Appendix B) are shown in Figure 15 on the right. The depolarization pattern is approximately well matched, and the voltage is in the correct range. However, no action potentials are generated by the addition of tolbutamide, whereas the transition from 15 mM to 40 mM KCl abolishes the preexistent action potentials, conforming with the measurements. For the cytosolic Ca2+ concentration, a moderate increase comparable to the measurement occurs with both the addition of tolbutamide and the subsequent KCl addition. In contrast to the measurement, however, the brief Ca2+ decrease after the addition of 40 mM KCl is missing, and the subsequent increase is much stronger in the model than in the measurement. This is still true when taking into consideration that the simulation gives cytosolic Ca2+ values, whereas the measurement shows that Fluo-4 fluorescence under-reports large increases due to its highly non-linear correlation. To illustrate this point, it is remarkable that the peak value of the Fluo-4 fluorescence in the presence of 15 mM KCl is virtually the same as the one in the presence of 40 mM KCl. Of note, in the control experiments, the peak value at 40 mM KCl was significantly higher than the one at 15 mM KCl (Figure 15, left panel).
Thus, the simulation can reproduce the features of this experiment in a first approximation. So, it was interesting to test whether it could be used to investigate the possible reasons as to why the current clamp mode affects simultaneously measured cytosolic Ca2+ concentrations.
In order to understand the possible mechanisms, detailed information about the electrical circuit of the measuring apparatus is necessary. A principle schematic is shown in Figure 16 [49]. In addition to the operational amplifiers and impedance converters, which provide the voltages that cause the specified current, the main electrical components are the resistance of the pipette and its capacitance. For the simulation, both the resistance and the capacitance of the pipette are to be understood as a series connection to the resistance and the capacitance of the cell membrane (compare Figure A1). Resistances connected in series are added. The resistance of the pipette is much smaller than the resistance of the cell membrane, and can therefore be neglected.
The situation is rather different with the capacitances. The capacitance of the pipette C p is in the range of 1–5 pF [49] and is thus in the order of magnitude of the capacitance of the cell (5.3 pF was assumed in the model, see Appendix B). A series connection of two capacitances is made by adding the reciprocal values.
The effective capacity is therefore
1 C r e s = 1 C m + 1 C p C r e s = C m C p C m + C p
Depending on whether one now assumes the minimum, maximum, or an average value of 3 pF for the pipette capacitance, this results in C r e s , m a x = 2.57   pF , C r e s , m i n = 0.84   pF and C r e s , m e a n = 1.92   pF .
In any case, the resulting capacitance is significantly below the membrane capacitance. In order to investigate this influence, tests were subsequently carried out on the influence of the “true” capacitance (see Equation (4)). The model showed instabilities at very small capacitances, which is why moderate changes in capacitance are considered below. In this context, the original value of 5.3 pF was set to the values 3.0 pF and 4.0 pF, respectively. The corresponding results are shown in Figure 17.
Comparing these curves with the simulation resulting from the use of the original capacitance value (Figure 15, right) reveals clear differences. The reduction of the capacitance results in a higher propensity of the system to produce action potentials, and at the same time, in a diminished increase of the cytosolic Ca2+ concentration. When no action potentials or action potentials with a small amplitude are present, the Ca2+ concentration approaches the same steady-state value as in the original simulations. This is also the reason that the response of the Ca2+ concentration to the transition from 15 mM to 40 mM KCl is nearly unchanged.
In conclusion, the measurement of the plasma membrane potential by current clamping appears to affect the simultaneous measurement of the cytosolic Ca2+ concentration because of the reduced effective capacitance of the system. This effect is the likely explanation for the observations depicted in Figure 14. Bearing this in mind, a suitable model should be able to subsequently calculate the effect of the measurement and to compensate for the feedback effect of the measurement system.

3. Materials and Methods

3.1. Chemicals

Collagenase NB8 (Serva, Heidelberg, Germany) and collagenase P from Roche (Sigma-Aldrich, Deisenhofen, Germany) were used for islet isolation. Tolbutamide was obtained from Sigma, and Fluo 4 from Thermo-Fisher Scientific (Schwerte, Germany). Cell culture medium RPMI 1640 (without glucose) was from Sigma, and fetal calf serum (FCS Gold ADD) was from Bio & Sell (Nürnberg-Feucht, Germany). All other reagents of analytical grade were from E. Merck (Darmstadt, Germany).

3.2. Tissue and Tissue Culture

Islets were isolated from the pancreas of NMRI mice (12–16 weeks old, fed ad libitum) by a collagenase digestion technique and hand-picked under a stereomicroscope. Single islet cells were obtained by incubation of the islets for 10 min in a Ca2+-free Krebs-Ringer medium and vortex mixing. Single cells were cultured in a humidified atmosphere of 95% air and 5% CO2 at 37 °C for 24 h in RPMI containing 10% FCS. The glucose concentration was 10 mM for the first 2 h, then 5 mM. Animal care is supervised by the regional authority (LAVES, Lower Saxony, Germany) and conforms to the current EU regulations.

3.3. Measurement of the Plasma Membrane Potential

The plasma membrane potential of 1-day- or 2-day-cultured single islet cells was measured using the patch clamp technique in the perforated patch configuration. Pipettes were pulled from borosilicate glass (2 mm o.d., 1.4 mm i.d., Hilgenberg, Germany) by a two-stage vertical puller (HEKA-Electronics, Lambrecht, Germany) and had resistances between 3 and 6 MΩ when filled with solution. The measurements were performed using an EPC 7 patch clamp amplifier (HEKA-Electronics) in the current clamp mode and pClamp 6.03 software (Axon Instruments, Foster City, CA, USA). The data were stored on a hard disk and analysed offline using GraphPad Prism5 software (GraphPad, LaJolla, CA, USA). A slow bath perfusion system was used, and all experiments were performed at room temperature (20–22 °C). The pipette solution contained the following (mM): 10 KCl, 10 NaCl, 70 K2SO4, 7 MgCl2, and 5 HEPES, pH 7.15, plus 125 µg/mL of the pore-forming agent nystatine (2.5% DMSO final concentration). The extracellular solution contained (mM) 140 NaCl, 5.6 KCl, 1.2 MgCl2, 2.6 CaCl2, 10 HEPES, and 1 glucose, pH 7.4.

3.4. Microfluorimetric Measurement of the Cytosolic Ca2+ Concentration

During the measurement of the plasma membrane potential, the cytosolic Ca2+ concentration of single islet cells was measured using Fluo 4 (loaded at a concentration of 1 µM for 20 min at 37 °C). Excitation was at 490 nm, the emission (>510 nm) was collected with a Leitz fluorite objective (40×, 1.3 N.A.) and registered by a photon-counting multiplier (Hamamatsu H6240-01) at 2 Hz under the control of LabView 7.1 (National Instruments, Munich, Germany).

3.5. Modelling

The advanced state of modeling, as represented by the “integrated oscillator model” (IOM), was chosen to test whether it has achieved general applicability to describe the changes observed in single beta cells (Figure 18) by independent experimentation. The number of eight state variables in this model (Table 3) permitted us to model various modes of bursting, which include fast bursting, slow bursting, and compound bursting [36,37]. The corresponding equations and parameters are listed in Table A1, Table A2 and Table A3 in Appendix B. All the following numerical studies have been carried out with the tool MATLAB R 2020b. Specifically, the present study aims to answer the following questions:
  • How robust is the model when tested considering novel experimental results?
  • How good is the response of the model to pharmacological agents?
  • Is the model helpful to understand the influence of the measuring technique?
The measurements, in which various agents are applied, are always designed in such a way that they start from a steady state. In order to be able to compare the measurement results with the simulation results, it is therefore very important to also start from the steady state in the simulations. This is particularly important because, in this highly non-linear system, there are dynamics on different time scales that can significantly influence the system. For this reason, a simulation was carried out for each simulation before the actual time integration, in which practically the agent-free state is characterised and in which a steady state is calculated over sufficiently long simulation times. This state forms the initial condition for the actual simulations. The numerical time integration has been implemented using the conventional MATLAB intern integrator ode45.

4. Conclusions and Outlook

In the present paper, the model of Marinelli et al., 2018 [36], used to describe the events of stimulus secretion-coupling within the pancreatic beta cell, was tested using published and unpublished measurements of the plasma membrane potential and the cytosolic Ca2+ concentration. We confirm that the model is well suited to represent the dynamics and interactions of membrane potential, Ca2+ transport, and the metabolism of glucose. Furthermore, the usefulness of simulating experimental results h shown by investigating the unexpected observation that measuring the membrane potential affects the simultaneous measurement of the cytosolic Ca2+ concentration.
However, it was also found that the model is very sensitive to slight variations in some parameters, and may have to be adjusted differently to correctly describe a wider range of observations obtained by other groups. A high degree of sensitivity is in particular due to the description of the permeability of the different ion channels. For example, in the case of the ATP-dependent potassium channels, there is a close interaction between tolbutamide and ADP or ATP that has not yet been taken into account in the model.
A similar example of yet-to-be-considered complexity is given by direct interactions between intracellular Ca2+ stores, such as those recently described by Klec et al. for a mitochondrially directed Ca2+ leak of the ER in beta cells [50,51]. So far, the interactions considered mathematically are based on the assumption that Caer and Cam are only linked via Ca, the free Ca2+ concentration in the bulk cytosol which equally affects processes at the plasma membrane. Therefore, this particular effect is not included in the model’s present stage of evolution, but can be considered in future works.
To achieve a generally valid description of the stimulus–secretion coupling within the beta cell, it seems desirable to include the effects of non-glucidic nutrient secretagogues such as alpha-ketoisocaproic acid, or the combination of glutamine with leucine. Since the effects of these stimuli include membrane depolarization and Ca2+ influx, their integration should be possible within the current conceptual framework. Another desirable feature to include in the model is the physiologically relevant endpoint, which is, after all, insulin secretion (see, e.g., the modeling of insulin secretion as a function of the actin cytoskeleton, [52]).
To achieve this goal, two further aspects have to be considered. First, Ca2+ microdomains in the vicinity of the Ca2+ channels are considered to be the immediate fusion signal, and may not be well represented in the bulk cytosolic Ca2+ concentration [53,54]. Based on recent work by Felix-Martinez et al. [55], spatially resolved calcium concentrations can also be calculated, and lead to very significant gradients within the cytosol (they are also correlated with insulin secretion models [56]). For future studies, this would also raise questions about the assumption of a homogenised consideration of the Ca2+ in the system of equations discussed here. Corresponding model extensions could be used as a possible approach to address this issue.
Second, the efficiency of the Ca2+ signal is modulated by the amplifying signals [6,7] which have thus far eluded unequivocal identification. So, the task of modeling the stimulus–secretion coupling of the pancreatic beta cell offers a large bandwidth of opportunities for further work.

Author Contributions

Conceptualization, M.M. and I.R.; methodology, M.M., J.W., N.S. and I.R.; software, M.M. and J.W.; validation, M.M., J.W. and I.R.; formal analysis, J.W.; investigation, M.M., J.W. and I.R.; resources, M.M. and I.R.; data curation, M.M., J.W., N.S. and I.R.; writing—original draft preparation, M.M. and I.R.; writing—review and editing, J.W. and N.S.; visualization, M.M. and J.W.; supervision, M.M. and I.R.; project administration, M.M. and I.R.; funding acquisition, I.R. All authors have read and agreed to the published version of the manuscript.

Funding

The authors’ research mentioned in this publication was supported by the Deutsche Forschungsgemeinschaft (Ru 368/5-1 to 5-4) and by grants from the Deutsche Diabetes Gesellschaft (to I.R.). We acknowledge support by the Open Access Publication Funds of Technische Universität Braunschweig.

Data Availability Statement

The data obtained in this study, including the MATLAB program, are available upon reasonable request.

Acknowledgments

Some of the measurements presented here were obtained during the thesis work of N.S. The authors thank Paul Thunich for restructuring the MATLAB code and Eric Lütgens for suggesting some of the study parameters carried out in this work.

Conflicts of Interest

The authors declare no conflict of interest.

Appendix A. Simplified Electrical Model for the Cell Membrane

The electrical events at the plasma membrane can be considered as taking place in the simplified circuit (see Figure A1) to allow further calculations. The standard model represents the membrane by its electrical resistance, capacity, and resting potential. In contrast to technical circuits, the charge is not transported by electrons flowing through a conductor, but by ions. These can have different charges, such as, e.g., calcium ions with their double positive charge [57].
Figure A1. Simplified circuit for the electrical system on a cell membrane during the patch-clamp process.
Figure A1. Simplified circuit for the electrical system on a cell membrane during the patch-clamp process.
Ijms 24 03206 g0a1
In this model, the total current flow I m is composed of the capacitive flow I C and the ionic current I i . is mainly determined by the specific membrane capacitance C m . This depends primarily on the cell type and is almost constant. The ion current depends on the resistance of the membrane R m , which is dependent on the transmembrane voltage and the Nernst potential E r , which is determined by the ion concentration on both sides of the membrane. For this circuit, this results in the following equation for the current I m :
I m = I C + I i = C m d V m d t + V m E r R m
The total ion current is composed of different ion currents, such as the flow of sodium, potassium, calcium, or chloride ions.

Appendix B. Governing Equations and Parameters

Table A1. The concrete equations as set up in [36].
Table A1. The concrete equations as set up in [36].
d V d t = 1 C m ( I C a + I K + I K ( C a ) + I K ( A T P ) ) (A2)
d n d t = n ( V ) n τ n (A3)
d C a d t = f C a ( J m e m J m J e r ) (A4)
d C a m d t = f C a σ m J m (A5)
d C a e r d t = f C a σ e r J e r (A6)
d A D P d t = A T P e x p ( 1 + 2.2 J P D H 0.05 + J P D H ) ( 1 C a 0.35 ) A D P τ A (A7)
d F B P d t = J P F K 1 2 J P D H (A8)
d F 6 P d t = 0.3 ( J G K J P F K ) (A9)
Table A2. Measures, additional formulae and meaning of terms related to the integrated oscillator model, taken from [36].
Table A2. Measures, additional formulae and meaning of terms related to the integrated oscillator model, taken from [36].
Measure and FormulaMeaning
I C a = g C a m ( V ) ( V V C a )
Current due to Calcium ion influx
m ( V ) = ( 1 + exp ( ( u m V ) / s m ) ) 1
Degree   of   opening   of   the   Ca 2 +   channels
I K = g K n ( V V K )
Current due to Potassium ion efflux
I K ( C a ) = g K ( C a ) q ( C a ) ( V V K )
Current due to Calcium-dependent potassium ion efflux
q ( C a ) = C a 2 / ( k d 2 + C a 2 )
Degree   of   opening   of   the   Calcium - dependent   K + channels
I K ( A T P ) = g K ( A T P ) o ( A D P , A T P ) ( V V K )
Current due to ATP-dependent potassium ion efflux
o ( A D P , A T P ) = 0.08 + 0.89 ( M g A D P k d d ) 2 + 0.16 ( M g A D P k d d ) ( 1 + M g A D P k d d ) 2 ( 1 + A T P 4 k t t + A D P 3 k t d )
ATP - dependent   degree   of   opening   of   K A T P channels
n ( V ) = ( 1 + exp ( ( u n V ) / s n ) ) 1
Steady-state potassium ion activity
J m e m = ( α I C a / V c y t + k P M C a C a )
Calcium influx into the cell
J e r = k S e r C a C a p l e a k ( C a e r C a )
Calcium influx into the ER
J m = k u n i C a k N a C a ( C a m C a )
Calcium influx into the mitochondria
J P F K = v P F K w 1110 + k P F K i , j , l { 0 , 1 } w i j 1 l i , j , k , l { 0 , 1 } w i j k l
PFK activity according to [33]
w i j k l = ( A M P / K 1 ) i ( F B P / K 2 ) j ( F 6 P / K 3 ) k ( A T P / K 4 ) l f 13 i k f 23 j k f 41 i l f 42 j l f 43 k l
Weight factor according to [33]
A M P = A D P 2 / A T P
Amount of adenosine monophosphate
J P D H = v P D H s ( C a m ) F B P
Pyravate dehydrogenase reaction rate
s ( C a m ) = C a m / ( K P D H + C a m )  
Feedback on pyravate dehydrogenase by calcium in the mitochondria
A T P = 1 2 ( A t o t + 4 A D P 2 + ( A t o t A D P ) 2 A D P )
Amount of ATP
Table A3. Parameter set for the reference model, taken from [36] and the reference therein.
Table A3. Parameter set for the reference model, taken from [36] and the reference therein.
C m = 5300   fF k t t = 1   µ M k P F K = 0.06
g C a = 1000   pS k t d = 26   µ M K 1 = 30   µ M
g K = 2700   pS f C a = 0.01 K 2 = 1   µ M
g K ( C a ) = 700   pS   ( 2.2.1 ) ,   50   pS (else) α = 5.18 × 10 18   µ mol / ( fA · ms ) K 3 = 5 × 10 4   µ M
g K ( A T P ) = 25000   pS V c y t = 1.15 × 10 12   l K 4 = 1 × 10 3   µ M
V C a = 25   mV k P M C A = 0.2   ms 1 f 13 = 0.02
V K = 75   mV k S E R C A = 0.4   ms 1 f 23 = 0.2
V K = 20   mV p l e a k = 2 × 10 4   ms 1 f 41 = 20
s m = 12   mV k u n i = 0.4   ms 1 f 42 = 20
u n = 16   mV k N a C a = 0.001   ms 1 f 43 = 20
s n = 5   mV σ m = 100 v P D H = 0.009   µ M / ms
τ n = 20   ms σ e r = 31 K P D H = 200   µ M
k d = 0.5   µ M J G K = 0.0015   µ M / ms ( 2.2.1 ) ,   J G K = 0.001   µ M / ms   ( else ) τ a = 3 × 10 5   ms
k d d = 17   µ M v P F K = 0.01   µ M / ms A t o t = 3 × 10 3   µ M

References

  1. Coore, H.G.; Randle, P.J. Regulation of insulin secretion studied with pieces of rabbit pancreas incubated in vitro. Biochem. J. 1964, 93, 66. [Google Scholar] [CrossRef]
  2. Malaisse, W.J. Insulin release: The fuel concept. Diabète Métabolisme 1983, 9, 313–320. [Google Scholar]
  3. Cook, D.L.; Hales, N. Intracellular ATP directly blocks K+ channels in pancreatic B-cells. Nature 1984, 311, 271–273. [Google Scholar] [CrossRef] [PubMed]
  4. Inagaki, N.; Gonoi, T.; Clement IV, J.P.; Namba, N.; Inazawa, J.; Gonzalez, G.; Aguilar-Bryan, L.; Seino, S.; Bryan, J. Reconstitution of I KATP: An inward rectifier subunit plus the sulfonylurea receptor. Science 1995, 270, 1166–1170. [Google Scholar] [CrossRef]
  5. Bergsten, P.; Grapengiesser, E.; Gylfe, E.; Tengholm, A.; Hellman, B. Synchronous oscillations of cytoplasmic Ca2+ and insulin release in glucose-stimulated pancreatic islets. J. Biol. Chem. 1994, 269, 8749–8753. [Google Scholar] [CrossRef]
  6. Henquin, J.-C. Triggering and amplifying pathways of regulation of insulin secretion by glucose. Diabetes 2000, 49, 1751–1760. [Google Scholar] [CrossRef] [PubMed]
  7. Rustenbeck, I.; Schulze, T.; Morsi, M.; Alshafei, M.; Panten, U. What is the metabolic amplification of insulin secretion and is it (still) relevant? Metabolites 2021, 11, 355. [Google Scholar] [CrossRef] [PubMed]
  8. Gaus, B.; Brüning, D.; Groß, S.; Müller, M.; Rustenbeck, I. The changing view of insulin granule mobility: From conveyor belt to signaling hub. Front. Endocrinol. 2022, 13, 983152. [Google Scholar] [CrossRef]
  9. Georgiadou, E.; Rutter, G.A. Control by Ca2+ of mitochondrial structure and function in pancreatic β-cells. Cell Calcium 2020, 91, 102282. [Google Scholar] [CrossRef]
  10. Hughes, E.; Lee, A.K.; Tse, A. Dominant role of sarcoendoplasmic reticulum Ca2+-ATPase pump in Ca2+ homeostasis and exocytosis in rat pancreatic β-cells. Endocrinology 2006, 147, 1396–1407. [Google Scholar] [CrossRef]
  11. Houamed, K.M.; Sweet, I.R.; Satin, L.S. BK channels mediate a novel ionic mechanism that regulates glucose-dependent electrical activity and insulin secretion in mouse pancreatic β-cells. J. Physiol. 2010, 588, 3511–3523. [Google Scholar] [CrossRef]
  12. Plášek, J.; Sigler, K. Slow fluorescent indicators of membrane potential: A survey of different approaches to probe response analysis. J. Photochem. Photobiol. B Biol. 1996, 33, 101–124. [Google Scholar] [CrossRef] [PubMed]
  13. Molleman, A. Patch Clamping: An Introductory Guide to Patch Clamp Electrophysiology; John Wiley & Sons: Hoboken, NJ, USA, 2003; ISBN 0470856513. [Google Scholar]
  14. Dean, P.M.; Matthews, E.K. Glucose-induced electrical activity in pancreatic islet cells. J. Physiol. 1970, 210, 255–264. [Google Scholar] [CrossRef] [PubMed]
  15. Henquin, J.-C.; Meissner, H.P.; Schmeer, W. Cyclic variations of glucose-induced electrical activity in pancreatic B cells. Pflüg. Arch. 1982, 393, 322–327. [Google Scholar] [CrossRef] [PubMed]
  16. Neher, E.; Sakmann, B. Single-channel currents recorded from membrane of denervated frog muscle fibres. Nature 1976, 260, 799–802. [Google Scholar] [CrossRef]
  17. Hamill, O.P.; Marty, A.; Neher, E.; Sakmann, B.; Sigworth, F.J. Improved patch-clamp techniques for high-resolution current recording from cells and cell-free membrane patches. Pflüg. Arch. 1981, 391, 85–100. [Google Scholar] [CrossRef]
  18. Horn, R.; Marty, A. Muscarinic activation of ionic currents measured by a new whole-cell recording method. J. Gen. Physiol. 1988, 92, 145–159. [Google Scholar] [CrossRef]
  19. Tsien, R.Y.; Rink, T.J. Neutral carrier ion-selective microelectrodes for measurement of intracellular free calcium. Biochim. Biophys. Acta (BBA)-Biomembr. 1980, 599, 623–638. [Google Scholar] [CrossRef]
  20. Grynkiewicz, G.; Poenie, M.; Tsien, R.Y. A new generation of Ca2+ indicators with greatly improved fluorescence properties. J. Biol. Chem. 1985, 260, 3440–3450. [Google Scholar] [CrossRef]
  21. McCombs, J.E.; Palmer, A.E. Measuring calcium dynamics in living cells with genetically encodable calcium indicators. Methods 2008, 46, 152–159. [Google Scholar] [CrossRef]
  22. Minta, A.; Kao, J.P.Y.; Tsien, R.Y. Fluorescent indicators for cytosolic calcium based on rhodamine and fluorescein chromophores. J. Biol. Chem. 1989, 264, 8171–8178. [Google Scholar] [CrossRef]
  23. Thomas, D.; Tovey, S.C.; Collins, T.J.; Bootman, M.D.; Berridge, M.J.; Lipp, P. A comparison of fluorescent Ca2+ indicator properties and their use in measuring elementary and global Ca2+ signals. Cell Calcium 2000, 28, 213–223. [Google Scholar] [CrossRef] [PubMed]
  24. Chay, T.R.; Keizer, J. Minimal model for membrane oscillations in the pancreatic beta-cell. Biophys. J. 1983, 42, 181–189. [Google Scholar] [CrossRef] [PubMed]
  25. Liu, Y.-J.; Tengholm, A.; Grapengiesser, E.; Hellman, B.; Gylfe, E. Origin of slow and fast oscillations of Ca2+ in mouse pancreatic islets. J. Physiol. 1998, 508, 471–481. [Google Scholar] [CrossRef]
  26. Rorsman, P.; Eliasson, L.; Kanno, T.; Zhang, Q.; Gopel, S. Electrophysiology of pancreatic β-cells in intact mouse islets of Langerhans. Prog. Biophys. Mol. Biol. 2011, 107, 224–235. [Google Scholar] [CrossRef]
  27. Kjems, L.L.; Ravier, M.A.; Jonas, J.-C.; Henquin, J.-C. Do oscillations of insulin secretion occur in the absence of cytoplasmic Ca2+ oscillations in β-cells? Diabetes 2002, 51, S177–S182. [Google Scholar] [CrossRef]
  28. Satin, L.S.; Butler, P.C.; Ha, J.; Sherman, A.S. Pulsatile insulin secretion, impaired glucose tolerance and type 2 diabetes. Mol. Asp. Med. 2015, 42, 61–77. [Google Scholar] [CrossRef]
  29. Zmazek, J.; Klemen, M.S.; Markovič, R.; Dolenšek, J.; Marhl, M.; Stožer, A.; Gosak, M. Assessing different temporal scales of calcium dynamics in networks of beta cell populations. Front. Physiol. 2021, 12, 612233. [Google Scholar] [CrossRef]
  30. Bertram, R.; Previte, J.; Sherman, A.; Kinard, T.A.; Satin, L.S. The phantom burster model for pancreatic β-cells. Biophys. J. 2000, 79, 2880–2892. [Google Scholar] [CrossRef] [PubMed]
  31. Bertram, R.; Sherman, A. Filtering of calcium transients by the endoplasmic reticulum in pancreatic β-cells. Biophys. J. 2004, 87, 3775–3785. [Google Scholar] [CrossRef] [PubMed]
  32. Bertram, R.; Sherman, A.; Satin, L.S. Metabolic and electrical oscillations: Partners in controlling pulsatile insulin secretion. Am. J. Physiol. -Endocrinol. Metab. 2007, 293, E890–E900. [Google Scholar] [CrossRef]
  33. Smolen, P. A model for glycolytic oscillations based on skeletal muscle phosphofructokinase kinetics. J. Theor. Biol. 1995, 174, 137–148. [Google Scholar] [CrossRef]
  34. McKenna, J.P.; Ha, J.; Merrins, M.J.; Satin, L.S.; Sherman, A.; Bertram, R. Ca2+ effects on ATP production and consumption have regulatory roles on oscillatory islet activity. Biophys. J. 2016, 110, 733–742. [Google Scholar] [CrossRef] [PubMed]
  35. Merrins, M.J.; van Dyke, A.R.; Mapp, A.K.; Rizzo, M.A.; Satin, L.S. Direct Measurements of Oscillatory Glycolysis in Pancreatic Islet β-Cells Using Novel Fluorescence Resonance Energy Transfer (FRET) Biosensors for Pyruvate Kinase M2 Activity. J. Biol. Chem. 2013, 288, 33312–33322. [Google Scholar] [CrossRef]
  36. Marinelli, I.; Vo, T.; Gerardo-Giorda, L.; Bertram, R. Transitions between bursting modes in the integrated oscillator model for pancreatic β-cells. J. Theor. Biol. 2018, 454, 310–319. [Google Scholar] [CrossRef] [PubMed]
  37. McKenna, J.P.; Bertram, R. Fast-slow analysis of the Integrated Oscillator Model for pancreatic β-cells. J. Theor. Biol. 2018, 457, 152–162. [Google Scholar] [CrossRef] [PubMed]
  38. Braun, M.; Ramracheya, R.; Bengtsson, M.; Zhang, Q.; Karanauskaite, J.; Partridge, C.; Johnson, P.R.; Rorsman, P. Voltage-gated ion channels in human pancreatic β-cells: Electrophysiological characterization and role in insulin secretion. Diabetes 2008, 57, 1618–1628. [Google Scholar] [CrossRef]
  39. Gilon, P.; Chae, H.-Y.; Rutter, G.A.; Ravier, M.A. Calcium signaling in pancreatic β-cells in health and in Type 2 diabetes. Cell Calcium 2014, 56, 340–361. [Google Scholar] [CrossRef]
  40. Marinelli, I.; Thompson, B.M.; Parekh, V.S.; Fletcher, P.A.; Gerardo-Giorda, L.; Sherman, A.S.; Satin, L.S.; Bertram, R. Oscillations in K (ATP) conductance drive slow calcium oscillations in pancreatic β-cells. Biophys. J. 2022, 121, 1449–1464. [Google Scholar] [CrossRef]
  41. Sanchez-Andres, J.V.; Gomis, A.; Valdeolmillos, M. The electrical activity of mouse pancreatic beta-cells recorded in vivo shows glucose-dependent oscillations. J. Physiol. 1995, 486, 223–228. [Google Scholar] [CrossRef]
  42. Detimary, P.; Jonas, J.-C.; Henquin, J.-C. Possible links between glucose-induced changes in the energy state of pancreatic B cells and insulin release. Unmasking by decreasing a stable pool of adenine nucleotides in mouse islets. J. Clin. Investig. 1995, 96, 1738–1745. [Google Scholar] [CrossRef] [PubMed]
  43. Brüning, D.; Morsi, M.; Früh, E.; Scherneck, S.; Rustenbeck, I. Metabolic Regulation of Hormone Secretion in Beta-Cells and Alpha-Cells of Female Mice: Fundamental Differences. Endocrinology 2022, 163, bqac125. [Google Scholar] [CrossRef]
  44. Rorsman, P.; Trube, G. Calcium and delayed potassium currents in mouse pancreatic beta-cells under voltage-clamp conditions. J. Physiol. 1986, 374, 531–550. [Google Scholar] [CrossRef] [PubMed]
  45. Hatlapatka, K.; Willenborg, M.; Rustenbeck, I. Plasma membrane depolarization as a determinant of the first phase of insulin secretion. Am. J. Physiol.-Endocrinol. Metab. 2009, 297, E315–E322. [Google Scholar] [CrossRef]
  46. Willenborg, M.; Hatlapatka, K.; Ghaly, H.; Belz, M.; Panten, U.; Rustenbeck, I. Studies of first phase insulin secretion using imposed plasma membrane depolarization. Front. Biosci.-Sch. 2011, 3, 662–679. [Google Scholar]
  47. Belz, M.; Willenborg, M.; Görgler, N.; Hamada, A.; Schumacher, K.; Rustenbeck, I. Insulinotropic effect of high potassium concentration beyond plasma membrane depolarization. Am. J. Physiol.-Endocrinol. Metab. 2014, 306, E697–E706. [Google Scholar] [CrossRef] [PubMed]
  48. Babenko, A.P.; Gonzalez, G.; Bryan, J. The tolbutamide site of SUR1 and a mechanism for its functional coupling to KATP channel closure. FEBS Lett. 1999, 459, 367–376. [Google Scholar] [CrossRef]
  49. Gold, R.S. The Axon Guide for Electrophysiology and Biophysics, Labaratory Techniques; Comp. Axon Instruments; Molecular Devices: Sunnyvale, CA, USA, 1993. [Google Scholar]
  50. Klec, C.; Madreiter-Sokolowski, C.T.; Stryeck, S.; Sachdev, V.; Duta-Mare, M.; Gottschalk, B.; Depaoli, M.R.; Rost, R.; Hay, J.; Waldeck-Weiermair, M. Glycogen synthase kinase 3 beta controls Presenilin-1-mediated endoplasmic reticulum Ca2+ leak directed to mitochondria in pancreatic islets and β-cells. Cell. Physiol. Biochem. Int. J. Exp. Cell. Physiol. Biochem. Pharmacol. 2019, 52, 57. [Google Scholar]
  51. Klec, C.; Madreiter-Sokolowski, C.T.; Ziomek, G.; Stryeck, S.; Sachdev, V.; Duta-Mare, M.; Gottschalk, B.; Depaoli, M.R.; Rost, R.; Hay, J. Presenilin-1 established ER-Ca2+ leak: A follow up on its importance for the initial insulin secretion in pancreatic islets and beta-cells upon elevated glucose. Cell. Physiol. Biochem. 2019, 53, 573–586. [Google Scholar]
  52. Müller, M.; Glombek, M.; Powitz, J.; Brüning, D.; Rustenbeck, I. A Cellular Automaton Model as a First Model-Based Assessment of Interacting Mechanisms for Insulin Granule Transport in Beta Cells. Cells 2020, 9, 1487. [Google Scholar] [CrossRef]
  53. Parekh, A.B. Ca2+ microdomains near plasma membrane Ca2+ channels: Impact on cell function. J. Physiol. 2008, 586, 3043–3054. [Google Scholar] [CrossRef] [PubMed]
  54. Hoppa, M.B.; Collins, S.; Ramracheya, R.; Hodson, L.; Amisten, S.; Zhang, Q.; Johnson, P.; Ashcroft, F.M.; Rorsman, P. Chronic palmitate exposure inhibits insulin secretion by dissociation of Ca2+ channels from secretory granules. Cell Metab. 2009, 10, 455–465. [Google Scholar] [CrossRef] [PubMed]
  55. Félix-Martínez, G.J.; Godínez-Fernández, J.R. Modeling the spatiotemporal distribution of Ca2+ during action potential firing in human pancreatic β-cells. Biomed. Phys. Eng. Express 2017, 3, 25020. [Google Scholar] [CrossRef]
  56. Félix-Martínez, G.J.; Godínez-Fernández, J.R. Effect of cell size on insulin secretion in human β-cells: A simulation study. Health Technol. 2020, 10, 379–383. [Google Scholar] [CrossRef]
  57. Johnston, D.; Wu, S.M.-S. Foundations of Cellular Neurophysiology; MIT Press: Cambridge, MA, USA, 1994; ISBN 0262293498. [Google Scholar]
Figure 1. Measurement of current–voltage correlations for Ca2+-channels (A) and K+-channels (B); for (A): ■ sum of all currents, ◯: L-Type, ▲ P/Q-Type ▽: T-Type; for (B): ● current without further agents, ◯: current in the presence of Co2+, Reprinted with permission from [38], 2008, American Diabetes Association.
Figure 1. Measurement of current–voltage correlations for Ca2+-channels (A) and K+-channels (B); for (A): ■ sum of all currents, ◯: L-Type, ▲ P/Q-Type ▽: T-Type; for (B): ● current without further agents, ◯: current in the presence of Co2+, Reprinted with permission from [38], 2008, American Diabetes Association.
Ijms 24 03206 g001
Figure 2. Currents for varying voltages, calculated according to Table A2 (Appendix B), left: current due to calcium ion influx ( I C a ) for varying Ca2+ equilibrium voltages, right: current due to Potassium ion efflux, sum and separated into the components I K ,   I K ( C a ) and I K ( A T P ) .
Figure 2. Currents for varying voltages, calculated according to Table A2 (Appendix B), left: current due to calcium ion influx ( I C a ) for varying Ca2+ equilibrium voltages, right: current due to Potassium ion efflux, sum and separated into the components I K ,   I K ( C a ) and I K ( A T P ) .
Ijms 24 03206 g002
Figure 3. Electric activity of a beta cell within an intact mouse islet in response to glucose, Reprinted with permission from [39], 2014, Elsevier. G10 stands for raising the extracellular glucose concentration from 3 mM to 10 mM, and G15 and G20 for raising to 15 mM and to 20 mM, respectively.
Figure 3. Electric activity of a beta cell within an intact mouse islet in response to glucose, Reprinted with permission from [39], 2014, Elsevier. G10 stands for raising the extracellular glucose concentration from 3 mM to 10 mM, and G15 and G20 for raising to 15 mM and to 20 mM, respectively.
Ijms 24 03206 g003
Figure 4. Simulation of the electrical activity as carried out with the model by [36] with varying A G C -values.
Figure 4. Simulation of the electrical activity as carried out with the model by [36] with varying A G C -values.
Ijms 24 03206 g004
Figure 5. ATP-ADP ratio over time for the processes simulated in Figure 4.
Figure 5. ATP-ADP ratio over time for the processes simulated in Figure 4.
Ijms 24 03206 g005
Figure 6. Effect of a slight change in resting Ca2+ potential; simulated curves for the membrane potential (blue) and Ca2+ concentration (orange) for (with respect to the reference configuration) (a) a slight decrease in V C a (b) a slight increase in V C a .
Figure 6. Effect of a slight change in resting Ca2+ potential; simulated curves for the membrane potential (blue) and Ca2+ concentration (orange) for (with respect to the reference configuration) (a) a slight decrease in V C a (b) a slight increase in V C a .
Ijms 24 03206 g006
Figure 7. Dynamics of the cytosolic Ca2+ concentration under different given extracellular potassium concentrations in the presence of 10 mM glucose and 100 μM diazoxide. , Reprinted with permission from [39], 2014, Elsevier.
Figure 7. Dynamics of the cytosolic Ca2+ concentration under different given extracellular potassium concentrations in the presence of 10 mM glucose and 100 μM diazoxide. , Reprinted with permission from [39], 2014, Elsevier.
Ijms 24 03206 g007
Figure 8. Course of the K+ equilibrium potential over the external potassium concentration and data points for the values used in the model.
Figure 8. Course of the K+ equilibrium potential over the external potassium concentration and data points for the values used in the model.
Ijms 24 03206 g008
Figure 9. Dynamics of Ca2+ concentration and corresponding membrane potential under different given extracellular potassium concentrations; reference curve calculated with the model from [36].
Figure 9. Dynamics of Ca2+ concentration and corresponding membrane potential under different given extracellular potassium concentrations; reference curve calculated with the model from [36].
Ijms 24 03206 g009
Figure 10. Dynamics of Ca2+ concentration and corresponding membrane potential under different given extracellular potassium concentrations, calculated with a higher value of o ( o ) = 0.03 .
Figure 10. Dynamics of Ca2+ concentration and corresponding membrane potential under different given extracellular potassium concentrations, calculated with a higher value of o ( o ) = 0.03 .
Ijms 24 03206 g010
Figure 11. Dynamics of Ca2+ concentration and corresponding membrane potential under different given extracellular potassium concentrations, calculated with a higher value of o ( o ) = 0.03 and highly reduced (1/10 of the reference values) intracellular Ca2+ transport parameters ( p l e a k = 2 × 10 5   ms 1 and k N a C a = 1 × 10 4   ms 1 ).
Figure 11. Dynamics of Ca2+ concentration and corresponding membrane potential under different given extracellular potassium concentrations, calculated with a higher value of o ( o ) = 0.03 and highly reduced (1/10 of the reference values) intracellular Ca2+ transport parameters ( p l e a k = 2 × 10 5   ms 1 and k N a C a = 1 × 10 4   ms 1 ).
Ijms 24 03206 g011
Figure 12. Membrane potential of a pancreatic beta cell in the presence of 500 µM tolbutamide and stepwise increases of the potassium concentration (Reprinted with permission from [47], 2014, the American Physiological Society).
Figure 12. Membrane potential of a pancreatic beta cell in the presence of 500 µM tolbutamide and stepwise increases of the potassium concentration (Reprinted with permission from [47], 2014, the American Physiological Society).
Ijms 24 03206 g012
Figure 13. Simulated membrane potential and cytosolic Ca2+ concentration in the presence of 500 µM tolbutamide and stepwise increases of the potassium concentration, considering (a) the original equilibrium potential function (Equation (2)), and (b) a modified equilibrium potential function (10 mV larger than the original).
Figure 13. Simulated membrane potential and cytosolic Ca2+ concentration in the presence of 500 µM tolbutamide and stepwise increases of the potassium concentration, considering (a) the original equilibrium potential function (Equation (2)), and (b) a modified equilibrium potential function (10 mV larger than the original).
Ijms 24 03206 g013
Figure 14. Measurement of the cytosolic Ca2+ concentration of a pancreatic beta cell in the presence of 500 µM tolbutamide and stepwise increases of the potassium concentration; left: with simultaneous current clamp measurement; right: with the patch pipette sealed on the cell, but without applied voltage. Data are means ± SEM of four experiments.
Figure 14. Measurement of the cytosolic Ca2+ concentration of a pancreatic beta cell in the presence of 500 µM tolbutamide and stepwise increases of the potassium concentration; left: with simultaneous current clamp measurement; right: with the patch pipette sealed on the cell, but without applied voltage. Data are means ± SEM of four experiments.
Ijms 24 03206 g014
Figure 15. Simultaneous measurement of the membrane potential and the cytosolic Ca2+ concentration of a pancreatic beta cell in the presence of 500 µM tolbutamide and stepwise increases of the potassium concentration; (left): original measurement; (right): simulated values.
Figure 15. Simultaneous measurement of the membrane potential and the cytosolic Ca2+ concentration of a pancreatic beta cell in the presence of 500 µM tolbutamide and stepwise increases of the potassium concentration; (left): original measurement; (right): simulated values.
Ijms 24 03206 g015
Figure 16. Principle sketch for patch clamp measurements.
Figure 16. Principle sketch for patch clamp measurements.
Ijms 24 03206 g016
Figure 17. Membrane potential and Ca2+ concentration under the addition of 500 µM tolbutamide and varying potassium concentrations with reduced resulting capacities.
Figure 17. Membrane potential and Ca2+ concentration under the addition of 500 µM tolbutamide and varying potassium concentrations with reduced resulting capacities.
Ijms 24 03206 g017
Figure 18. Interactions considered mathematically in the integrated oscillator model—red: eight state variables; black: associated flows, transport and reaction kinetics; green: explicit mathematical interactions; violet: mathematical coupling between the state variables. The direction of the arrow indicates which state variable generates immediate input into the dynamics of the respective other state variables. The state variables interact with each other via currents, transports and reaction kinetics (each in blackletters). The associated equations, measures, and parameter values are listed in Table A1, Table A2 and Table A3 in Appendix B.
Figure 18. Interactions considered mathematically in the integrated oscillator model—red: eight state variables; black: associated flows, transport and reaction kinetics; green: explicit mathematical interactions; violet: mathematical coupling between the state variables. The direction of the arrow indicates which state variable generates immediate input into the dynamics of the respective other state variables. The state variables interact with each other via currents, transports and reaction kinetics (each in blackletters). The associated equations, measures, and parameter values are listed in Table A1, Table A2 and Table A3 in Appendix B.
Ijms 24 03206 g018
Table 1. Characteristic values for the Ca2+ concentration in arbitrary units, as derived from Figure 7.
Table 1. Characteristic values for the Ca2+ concentration in arbitrary units, as derived from Figure 7.
C a m a x   C a m i n D = C a m a x C a m i n D K 30 / D K , i
K300.300.120.181
K200.220.090.131.38
K150.140.080.063.0
Table 2. Characteristic values for the simulated Ca2+ concentration and comparison to the measurement in [39].
Table 2. Characteristic values for the simulated Ca2+ concentration and comparison to the measurement in [39].
Ref .   Simulation   o 0.006 Simulation   with   o = 0.03 Measurement
C a m a x C a m a x , K 30 / C a m a x , K , i C a m a x C a m a x , K 30 / C a m a x , K , i D K 30 / D K , i
K300.5010.4611
K200.461.090.331.391.38
K150.431.160.133.543.0
Table 3. Eight state variables of the integrated oscillator model and their meaning.
Table 3. Eight state variables of the integrated oscillator model and their meaning.
State VariableMeaning
V Membrane potential
n Potassium ion activity
C a Calcium ion concentration in the cytosol
C a m Calcium ion concentration in the mitochondria
C a e r Calcium ion concentration in the endoplasmic reticulum
F 6 P Amount of fructose-6-phosphate
F B P Amount of fructose-6-biphosphate
A D P Amount of ADP
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Müller, M.; Walkling, J.; Seemann, N.; Rustenbeck, I. The Dynamics of Calcium Signaling in Beta Cells—A Discussion on the Comparison of Experimental and Modelling Data. Int. J. Mol. Sci. 2023, 24, 3206. https://doi.org/10.3390/ijms24043206

AMA Style

Müller M, Walkling J, Seemann N, Rustenbeck I. The Dynamics of Calcium Signaling in Beta Cells—A Discussion on the Comparison of Experimental and Modelling Data. International Journal of Molecular Sciences. 2023; 24(4):3206. https://doi.org/10.3390/ijms24043206

Chicago/Turabian Style

Müller, Michael, Jonas Walkling, Nele Seemann, and Ingo Rustenbeck. 2023. "The Dynamics of Calcium Signaling in Beta Cells—A Discussion on the Comparison of Experimental and Modelling Data" International Journal of Molecular Sciences 24, no. 4: 3206. https://doi.org/10.3390/ijms24043206

APA Style

Müller, M., Walkling, J., Seemann, N., & Rustenbeck, I. (2023). The Dynamics of Calcium Signaling in Beta Cells—A Discussion on the Comparison of Experimental and Modelling Data. International Journal of Molecular Sciences, 24(4), 3206. https://doi.org/10.3390/ijms24043206

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop