e Young Cha

Ionic mechanisms and Ca2+ dynamics underlying the glucose response of pancreatic β cells: a simulation study.

To clarify the mechanisms underlying the pancreatic β-cell response to varying glucose concentrations ([G]), electrophysiological findings were integrated into a mathematical cell model. The Ca2+ dynamics of the endoplasmic reticulum (ER) were also improved. The model was validated by demonstrating quiescent potential, burst–interburst electrical events accompanied by Ca2+ 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 Ca2+ 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 Ca2+- or Na+-dependent currents, which were generated by the plasma membrane Ca2+ pump, Na+/K+ pump, Na+/Ca2+ exchanger, and TRPM channel. Accumulation and release of Ca2+ 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.

Time-dependent changes in membrane excitability during glucose-induced bursting activity in pancreatic β cells.

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.

A Novel Method to Quantify Contribution of Channels and Transporters to Membrane Potential Dynamics

The action potential, once triggered in ventricular or atrial myocytes, automatically proceeds on its time course or is generated spontaneously in sinoatrial node pacemaker cells. It is induced by complex interactions among such cellular components as ion channels, transporters, intracellular ion concentrations, and signaling molecules. We have developed what is, to our knowledge, a new method using a mathematical model to quantify the contribution of each cellular component to the automatic time courses of the action potential. In this method, an equilibrium value, which the membrane potential is approaching at a given moment, is calculated along the time course of the membrane potential. The calculation itself is based on the time-varying conductance and the reversal potentials of individual ion channels and electrogenic ion transporters. Since the equilibrium potential moves in advance of the membrane potential change, we refer to it as the lead potential, V(L). The contribution of an individual current was successfully quantified by comparing dV(L)/dt before and after fixing the time-dependent change of a component of interest, such as the variations in the open probability of a channel or the turnover rate of an ion transporter. In addition to the action potential, the lead-potential analysis should also be applicable in all types of membrane excitation in many different kinds of cells.

Analyzing electrical activities of pancreatic β cells using mathematical models.

Bursts of repetitive action potentials are closely related to the regulation of glucose-induced insulin secretion in pancreatic β cells. Mathematical studies with simple β-cell models have established the central principle that the burst-interburst events are generated by the interaction between fast membrane excitation and slow cytosolic components. Recently, a number of detailed models have been developed to simulate more realistic β cell activity based on expanded findings on biophysical characteristics of cellular components. However, their complex structures hinder our intuitive understanding of the underlying mechanisms, and it is becoming more difficult to dissect the role of a specific component out of the complex network. We have recently developed a new detailed model by incorporating most of ion channels and transporters recorded experimentally (the Cha-Noma model), yet the model satisfies the charge conservation law and reversible responses to physiological stimuli. Here, we review the mechanisms underlying bursting activity by applying mathematical analysis tools to representative simple and detailed models. These analyses include time-based simulation, bifurcation analysis and lead potential analysis. In addition, we introduce a new steady-state I-V (ssI-V) curve analysis. We also discuss differences in electrical signals recorded from isolated single cells or from cells maintaining electrical connections within multi-cell preparations. Towards this end, we perform simulations with our detailed pancreatic β-cell model.

Characterization of the cardiac Na+/K+ pump by development of a comprehensive and mechanistic model

A large amount of experimental data on the characteristics of the cardiac Na(+)/K(+) pump have been accumulated, but it remains difficult to predict the quantitative contribution of the pump in an intact cell because most measurements have been made under non-physiological conditions. To extrapolate the experimental findings to intact cells, we have developed a comprehensive Na(+)/K(+) pump model based on the thermodynamic framework (Smith and Crampin, 2004) of the Post-Albers reaction cycle combined with access channel mechanisms. The new model explains a variety of experimental results for the Na(+)/K(+) pump current (I(NaK)), including the dependency on the concentrations of Na(+) and K(+), the membrane potential and the free energy of ATP hydrolysis. The model demonstrates that both the apparent affinity and the slope of the substrate-I(NaK) relationship measured experimentally are affected by the composition of ions in the extra- and intracellular solutions, indirectly through alteration in the probability distribution of individual enzyme intermediates. By considering the voltage dependence in the Na(+)- and K(+)-binding steps, the experimental voltage-I(NaK) relationship could be reconstructed with application of experimental ionic compositions in the model, and the view of voltage-dependent K(+) binding was supported. Re-evaluation of charge movements accompanying Na(+) and K(+) translocations gave a reasonable number for the site density of the Na(+)/K(+) pump on the membrane. The new model is relevant for simulation of cellular functions under various interventions, such as depression of energy metabolism.

A Model of Na+/H+ Exchanger and Its Central Role in Regulation of pH and Na+ in Cardiac Myocytes

A new kinetic model of the Na(+)/H(+) exchanger (NHE) was developed by fitting a variety of major experimental findings, such as ion-dependencies, forward/reverse mode, and the turnover rate. The role of NHE in ion homeostasis was examined by implementing the NHE model in a minimum cell model including intracellular pH buffer, Na(+)/K(+) pump, background H(+), and Na(+) fluxes. This minimum cell model was validated by reconstructing recovery of pH(i) from acidification, accompanying transient increase in [Na(+)](i) due to NHE activity. Based on this cell model, steady-state relationships among pH(i), [Na(+)](I), and [Ca(2+)](i) were quantitatively determined, and thereby the critical level of acidosis for cell survival was predicted. The acidification reported during partial blockade of the Na(+)/K(+) pump was not attributed to a dissipation of the Na(+) gradient across the membrane, but to an increase in indirect H(+) production. This NHE model, though not adapted to the dimeric behavioral aspects of NHE, can provide a strong clue to quantitative prediction of degree of acidification and accompanying disturbance of ion homeostasis under various pathophysiological conditions.

Quantitative Decomposition of Dynamics of Mathematical Cell Models: Method and Application to Ventricular Myocyte Models

Mathematical cell models are effective tools to understand cellular physiological functions precisely. For detailed analysis of model dynamics in order to investigate how much each component affects cellular behaviour, mathematical approaches are essential. This article presents a numerical analysis technique, which is applicable to any complicated cell model formulated as a system of ordinary differential equations, to quantitatively evaluate contributions of respective model components to the model dynamics in the intact situation. The present technique employs a novel mathematical index for decomposed dynamics with respect to each differential variable, along with a concept named instantaneous equilibrium point, which represents the trend of a model variable at some instant. This article also illustrates applications of the method to comprehensive myocardial cell models for analysing insights into the mechanisms of action potential generation and calcium transient. The analysis results exhibit quantitative contributions of individual channel gating mechanisms and ion exchanger activities to membrane repolarization and of calcium fluxes and buffers to raising and descending of the cytosolic calcium level. These analyses quantitatively explicate principle of the model, which leads to a better understanding of cellular dynamics.

Ionic Basis of the Pacemaker Activity of SA Node Revealed by the Lead Potential Analysis

Ionic mechanisms of spontaneous action potential in sinoatrial (SA) node pacemaker cells have been discussed for decades. Although a number of theoretical studies have proposed different mathematical models, no scientific consensus has been achieved yet, because of the complexity and variations in experimental findings used for developing models. Here, we introduce a theoretical method in simulation study, the lead potential analysis, which enabled us to isolate the contribution of individual currents from the secondary effect of modified channel activities. We compared three models, suggesting different ionic mechanisms (Himeno et al. model, Kurata et al. model, and Maltsev and Lakatta model), and contributions of Ca2+ through activation of I NaCa is estimated. Finally, the effect of catecholamine stimulation is discussed based on a SA node cell model with β1-adrenergic signaling cascade and mechanisms of the positive chronotropy are analyzed.

Steady-state solutions of cell volume in a cardiac myocyte model elaborated for membrane excitation, ion homeostasis and Ca2+ dynamics.

The cell volume continuously changes in response to varying physiological conditions, and mechanisms underlying volume regulation have been investigated in both experimental and theoretical studies. Here, general formulations concerning cell volume change are presented in the context of developing a comprehensive cell model which takes Ca(2+) dynamics into account. Explicit formulas for charge conservation and steady-state volumes of the cytosol and endoplasmic reticulum (ER) are derived in terms of membrane potential, amount of ions, Ca(2+)-bound buffer molecules, and initial cellular conditions. The formulations were applied to a ventricular myocyte model which has plasma-membrane Ca(2+) currents with dynamic gating mechanisms, Ca(2+)-buffering reactions with diffusive and non-diffusive buffer proteins, and Ca(2+) uptake into or release from the sarcoplasmic reticulum (SR) accompanied by compensatory cationic or anionic currents through the SR membrane. Time-dependent volume changes in cardiac myocytes induced by varying extracellular osmolarity or by action potential generation were successfully simulated by the novel formulations. Through application of bifurcation analysis, the existence and uniqueness of steady-state solutions of the cell volume were validated, and contributions of individual ion channels and transporters to the steady-state volume were systematically analyzed. The new formulas are consistent with previous fundamental theory derived from simple models of minimum compositions. The new formulations may be useful for examination of the relationship between cell function and volume change in other cell types.

Modeling the cardiac Na+/H+ exchanger based on major experimental findings

Na+-H+ exchanger (NHE) is the main acid extruder in cardiac myocytes. We review the experimental findings of iondependency of NHE activity, and the mathematical modeling developed so far. In spite of extensive investigation, many unsolved questions still remain. We consider that the precise description of NHE activity with mathematical models elucidates the roles of NHE in maintaining ionic homeostasis, especially under pathophysiological conditions.