To clarify the mechanisms underlying the pancreatic β-cell response to varying glucose concentrations ([G]), electrophysiological findings were integrated into a mathematical cell model. The Ca²⁺ dynamics of the endoplasmic reticulum (ER) were also improved. The model was validated by demonstrating quiescent potential, burst–interburst electrical events accompanied by Ca²⁺ transients, and continuous firing of action potentials over [G] ranges of 0–6, 7–18, and >19 mM, respectively. These responses to glucose were completely reversible. The action potential, input impedance, and Ca²⁺ transients were in good agreement with experimental measurements. The ionic mechanisms underlying the burst–interburst rhythm were investigated by lead potential analysis, which quantified the contributions of individual current components. This analysis demonstrated that slow potential changes during the interburst period were attributable to modifications of ion channels or transporters by intracellular ions and/or metabolites to different degrees depending on [G]. The predominant role of adenosine triphosphate–sensitive K⁺ current in switching on and off the repetitive firing of action potentials at 8 mM [G] was taken over at a higher [G] by Ca²⁺- or Na⁺-dependent currents, which were generated by the plasma membrane Ca²⁺ pump, Na⁺/K⁺ pump, Na⁺/Ca²⁺ exchanger, and TRPM channel. Accumulation and release of Ca²⁺ by the ER also had a strong influence on the slow electrical rhythm. We conclude that the present mathematical model is useful for quantifying the role of individual functional components in the whole cell responses based on experimental findings.

In our companion paper, the physiological functions of pancreatic β cells were analyzed with a new β-cell model by time-based integration of a set of differential equations that describe individual reaction steps or functional components based on experimental studies. In this study, we calculate steady-state solutions of these differential equations to obtain the limit cycles (LCs) as well as the equilibrium points (EPs) to make all of the time derivatives equal to zero. The sequential transitions from quiescence to burst–interburst oscillations and then to continuous firing with an increasing glucose concentration were defined objectively by the EPs or LCs for the whole set of equations. We also demonstrated that membrane excitability changed between the extremes of a single action potential mode and a stable firing mode during one cycle of bursting rhythm. Membrane excitability was determined by the EPs or LCs of the membrane subsystem, with the slow variables fixed at each time point. Details of the mode changes were expressed as functions of slowly changing variables, such as intracellular [ATP], [Ca2+], and [Na+]. In conclusion, using our model, we could suggest quantitatively the mutual interactions among multiple membrane and cytosolic factors occurring in pancreatic β cells.

Glucagon-like peptide-1 (GLP-1) elevates intracellular concentration of cAMP ([cAMP]) and facilitates glucose-dependent insulin secretion in pancreatic β-cells. There has been much evidence to suggest that multiple key players such as the GLP-1 receptor, G_sprotein, adenylate cyclase (AC), phosphodiesterase (PDE), and intracellular Ca²⁺concentration ([Ca²⁺]) are involved in the regulation of [cAMP]. However, because of complex interactions among these signaling factors, the kinetics of the reaction cascade as well as the activities of ACs and PDEs have not been determined in pancreatic β-cells. We have constructed a minimal mathematical model of GLP-1 receptor signal transduction based on experimental findings obtained mostly in β-cells and insulinoma cell lines. By fitting this theoretical reaction scheme to key experimental records of the GLP-1 response, the parameters determining individual reaction steps were estimated. The model reconstructed satisfactorily the dynamic changes in [cAMP] and predicted the activities of cAMP effectors, protein kinase A (PKA), and cAMP-regulated guanine nucleotide exchange factor [cAMP-GEF or exchange protein directly activated by cAMP (Epac)] during GLP-1 stimulation. The simulations also predicted the presence of two sequential desensitization steps of the GLP1 receptor that occur with fast and very slow reaction rates. The cross talk between glucose- and GLP-1-dependent signal cascades for cAMP synthesis was well reconstructed by integrating the direct regulation of AC and PDE by [Ca²⁺]. To examine robustness of the signaling system in controlling [cAMP], magnitudes of AC and PDE activities were compared in the presence or absence of GLP-1 and/or the PDE inhibitor IBMX.¹