Modeling the spatiotemporal distribution of Ca 2+ during action potential firing in human pancreatic β -cells

Model of Calcium Dynamics Regulating [Formula: see text], ATP and Insulin Production in a Pancreatic [Formula: see text]-Cell

The calcium signals regulate the production and secretion of many signaling molecules like inositol trisphosphate ([Formula: see text]) and adenosine triphosphate (ATP) in various cells including pancreatic [Formula: see text]-cells. The calcium signaling mechanisms regulating [Formula: see text], ATP and insulin responsible for various functions of [Formula: see text]-cells are still not well understood. Any disturbance in these mechanisms can alter the functions of [Formula: see text]-cells leading to diabetes and metabolic disorders. Therefore, a mathematical model is proposed by incorporating the reaction-diffusion equation for calcium dynamics and a system of first-order differential equations for [Formula: see text], ATP-production and insulin secretion with initial and boundary conditions. The model incorporates the temporal dependence of [Formula: see text]-production and degradation, ATP production and insulin secretion on calcium dynamics in a [Formula: see text]-cell. The piecewise linear finite element method has been used for the spatial dimension and the Crank-Nicolson scheme for the temporal dimension to obtain numerical results. The effect of changes in source influxes and buffers on calcium dynamics and production of [Formula: see text], ATP and insulin levels in a [Formula: see text]-cell has been analyzed. It is concluded that the dysfunction of source influx and buffers can cause significant variations in calcium levels and dysregulation of [Formula: see text], ATP and insulin production, which can lead to various metabolic disorders, diabetes, obesity, etc. The proposed model provides crucial information about the changes in mechanisms of calcium dynamics causing proportionate disturbances in [Formula: see text], ATP and insulin levels in pancreatic cells, which can be helpful for devising protocols for diagnosis and treatment of various metabolic diseases.

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

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 Ca²⁺ 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 Ca²⁺ 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 Ca²⁺ 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.

Glucose-induced [Ca2+]i oscillations in β cells are composed of trains of spikes within a subplasmalemmal microdomain

Electrical activity and oscillations of cytosolic Ca²⁺ concentrations ([Ca²⁺]_i) that trigger insulin release in response to glucose are key functions of pancreatic β cells. Although oscillatory Ca²⁺ signals have been intensively studied in β cells, their lower frequency did not match that of electrical activity. In addition, the measured peak [Ca²⁺]_i did not reach levels that are typically required by synaptotagmins to elicit the release of insulin-containing vesicles in live-cell experiments. We therefore sought to resolve the Ca²⁺ dynamics in the subplasmalemmal microdomain that is critical for triggering fast exocytosis. Applying total internal reflection fluorescence (TIRF) microscopy in insulin-producing INS-1E and primary mouse β cells, we resolved extraordinary fast trains of Ca²⁺ spiking (frequency > 3 s^-1) in response to glucose exposure. Using a low-affinity [Ca²⁺]_i indicator dye, we provide experimental evidence that Ca²⁺ spikes reach low micromolar apparent concentrations in the vicinity of the plasma membrane_. Analysis of Ca²⁺ spikes evoked by repeated depolarization for 10 ms closely matched the Ca²⁺ dynamics observed upon glucose application. To our knowledge, this is the first study that experimentally demonstrates Ca²⁺ spikes in β cells with velocities that resemble those of bursting or continuously appearing trains of action potentials (APs) in non-patched cells.

Mata A, Godínez-Fernández JR. Reconstructing human pancreatic islet architectures using computational optimization

We outline a general methodology based on computational optimization and experimental data to reconstruct human pancreatic islet architectures. By using the nuclei coordinates of islet cells obtained through DAPI staining, cell types identified by immunostaining, and cell size distributions estimated from capacitance measurements, reconstructed islets composed of non-overlapping spherical cells were obtained through an iterative optimization procedure. In all cases, the reconstructed architectures included >99% of the experimental identified cells, each of them having a radius within the experimentally reported ranges. Given the wide use of mathematical modeling for the study of pancreatic cells, and recently, of cell-cell interactions within the pancreatic islets, the methodology here proposed, also capable of identifying cell-to-cell contacts, is aimed to provide with a framework for modeling and analyzing experimentally-based pancreatic islet architectures.

Electrophysiological models of the human pancreatic δ-cell: From single channels to the firing of action potentials

Minimal mathematical models were developed to describe the electrophysiological properties of human δ-cells. Markov models of single channels were first developed based on the analysis of electrophysiological data. Monte Carlo simulations of voltage-clamp experiments were performed in an iteratively optimization procedure to estimate the number of channels required to reproduce the main characteristics of the macroscopic currents recorded experimentally. A membrane model of the firing of action potentials was then developed based on the kinetic schemes of single channels and the number of channels estimated. We showed that macroscopic currents of human δ-cells can be reproduced by minimal models of single channels when the appropriate number of channels is considered. In addition, our simulations suggest that human δ-cells are capable of generating action potentials through the interaction of the ionic currents involved. Finally, we determined the relative contribution of the currents underlying the firing of action potentials in human pancreatic δ-cells, which allowed us to propose a qualitative model of an action potential in terms of the underlying ionic currents.