Geometric slow-fast analysis of a hybrid pituitary cell model with stochastic ion channel dynamics

To obtain explicit understanding of the behavior of dynamical systems, geometrical methods and slow-fast analysis have proved to be highly useful. Such methods are standard for smooth dynamical systems, and increasingly used for continuous, non-smooth dynamical systems. However, they are much less used for random dynamical systems, in particular for hybrid models with discrete, random dynamics. Indeed, the analysis of such systems has typically been done by studying the corresponding deterministic system and considering how noise perturbs the deterministic geometrical structures. Here we propose a geometrical method that works directly with the hybrid system. We illustrate our approach through an application to a hybrid pituitary cell model in which the stochastic dynamics of very few active large-conductance potassium (BK) channels is coupled to a deterministic model of the other ion channels and calcium dynamics. To employ our geometric approach, we exploit the slow-fast structure of the model. The random fast subsystem is analyzed by considering discrete phase planes, corresponding to the discrete number of open BK channels, and stochastic events correspond to jumps between these planes. The evolution within each plane can be understood from nullclines and limit cycles, and the overall dynamics, e.g., whether the model produces a spike or a burst, is determined by the location at which the system jumps from one plane to another. Our approach is generally applicable to other scenarios to study discrete random dynamical systems defined by hybrid stochastic-deterministic models.Comment: 15 pages, 8 figure

Dapagliflozin stimulates glucagon secretion at high glucose: experiments and mathematical simulations of human A-cells.

Glucagon is one of the main regulators of blood glucose levels and dysfunctional stimulus secretion coupling in pancreatic A-cells is believed to be an important factor during development of diabetes. However, regulation of glucagon secretion is poorly understood. Recently it has been shown that Na(+)/glucose co-transporter (SGLT) inhibitors used for the treatment of diabetes increase glucagon levels in man. Here, we show experimentally that the SGLT2 inhibitor dapagliflozin increases glucagon secretion at high glucose levels both in human and mouse islets, but has little effect at low glucose concentrations. Because glucagon secretion is regulated by electrical activity we developed a mathematical model of A-cell electrical activity based on published data from human A-cells. With operating SGLT2, simulated glucose application leads to cell depolarization and inactivation of the voltage-gated ion channels carrying the action potential, and hence to reduce action potential height. According to our model, inhibition of SGLT2 reduces glucose-induced depolarization via electrical mechanisms. We suggest that blocking SGLTs partly relieves glucose suppression of glucagon secretion by allowing full-scale action potentials to develop. Based on our simulations we propose that SGLT2 is a glucose sensor and actively contributes to regulation of glucagon levels in humans which has clinical implications

Recent advances in mathematical modeling and statistical analysis of exocytosis in endocrine cells

open5noMost endocrine cells secrete hormones as a result of Ca(2+)-regulated exocytosis, i.e., fusion of the membranes of hormone-containing secretory granules with the cell membrane, which allows the hormone molecules to escape to the extracellular space. As in neurons, electrical activity and cell depolarization open voltage-sensitive Ca(2+) channels, and the resulting Ca(2+) influx elevate the intracellular Ca(2+) concentration, which in turn causes exocytosis. Whereas the main molecular components involved in exocytosis are increasingly well understood, quantitative understanding of the dynamical aspects of exocytosis is still lacking. Due to the nontrivial spatiotemporal Ca(2+) dynamics, which depends on the particular pattern of electrical activity as well as Ca(2+) channel kinetics, exocytosis is dependent on the spatial arrangement of Ca(2+) channels and secretory granules. For example, the creation of local Ca(2+) microdomains, where the Ca(2+) concentration reaches tens of µM, are believed to be important for triggering exocytosis. Spatiotemporal simulations of buffered Ca(2+) diffusion have provided important insight into the interplay between electrical activity, Ca(2+) channel kinetics, and the location of granules and Ca(2+) channels. By confronting simulations with statistical time-to-event (or survival) regression analysis of single granule exocytosis monitored with TIRF microscopy, a direct connection between location and rate of exocytosis can be obtained at the local, single-granule level. To get insight into whole-cell secretion, simplifications of the full spatiotemporal dynamics have shown to be highly helpful. Here, we provide an overview of recent approaches and results for quantitative analysis of Ca(2+) regulated exocytosis of hormone-containing granules.openPedersen, Morten Gram; Tagliavini, Alessia; Cortese, Giuliana; Riz, Michela; Montefusco, FrancescoPedersen, MORTEN GRAM; Tagliavini, Alessia; Cortese, Giuliana; Riz, Michela; Montefusco, Francesc

Modeling Doxorubicin Pharmacokinetics in Multiple Myeloma Suggests Mechanism of Drug Resistance

Objective: Multiple myeloma (MM) is a plasma cell malignancy often treated with chemotherapy drugs. Among these, doxorubicin (DOXO) is commonly employed, sometimes in combined-drug therapies, but it has to be optimally administered in order to maximize its efficacy and reduce possible side effects. To support DOXO studies and treatment optimization, here we propose an experimental/modeling approach to establish a model describing DOXO pharmacokinetics (PK) in MM cells. Methods: A series of in vitro experiments were performed in MM1R and MOLP-2 cells. DOXO was administered at two dosages (200 nM, 450 nM) at [Formula: see text]=0 and removed at [Formula: see text]=3 hrs. Intracellular DOXO concentration was measured via fluorescence microscopy during both drug uptake ([Formula: see text]=0-3 hrs) and release phases ([Formula: see text]=3-8 hrs). Four PK candidate models were identified, and were compared and selected based on their ability to describe DOXO data and numerical parameter identification. Results: The most parsimonious model consists of three compartments describing DOXO distribution between the extracellular space, the cell cytoplasm and the nucleus, and defines the intracellular DOXO efflux rate through a Hill function, simulating a threshold/saturation drug resistance mechanism. This model predicted DOXO data well in all the experiments and provided precise parameter estimates (mean ± standard deviation coefficient of variation: 15.8±12.2%). Conclusions: A reliable PK model describing DOXO uptake and release in MM cells has been successfully developed. Significance: The proposed PK model, once integrated with DOXO pharmacodynamics, has the potential of allowing the study and the optimization of DOXO treatment strategies in MM

Complex Patterns of Metabolic and Ca²⁺ Entrainment in Pancreatic Islets by Oscillatory Glucose

Glucose-stimulated insulin secretion is pulsatile and driven by intrinsic oscillations in metabolism, electrical activity, and Ca(2+) in pancreatic islets. Periodic variations in glucose can entrain islet Ca(2+) and insulin secretion, possibly promoting interislet synchronization. Here, we used fluorescence microscopy to demonstrate that glucose oscillations can induce distinct 1:1 and 1:2 entrainment of oscillations (one and two oscillations for each period of exogenous stimulus, respectively) in islet Ca(2+), NAD(P)H, and mitochondrial membrane potential. To our knowledge, this is the first demonstration of metabolic entrainment in islets, and we found that entrainment of metabolic oscillations requires voltage-gated Ca(2+) influx. We identified diverse patterns of 1:2 entrainment and showed that islet synchronization during entrainment involves adjustments of both oscillatory phase and period. All experimental findings could be recapitulated by our recently developed mathematical model, and simulations suggested that interislet variability in 1:2 entrainment patterns reflects differences in their glucose sensitivity. Finally, our simulations and recordings showed that a heterogeneous group of islets synchronized during 1:2 entrainment, resulting in a clear oscillatory response from the collective. In summary, we demonstrate that oscillatory glucose can induce complex modes of entrainment of metabolically driven oscillations in islets, and provide additional support for the notion that entrainment promotes interislet synchrony in the pancreas

Complex Patterns of Metabolic and Ca²⁺ Entrainment in Pancreatic Islets by Oscillatory Glucose

Glucose-stimulated insulin secretion is pulsatile and driven by intrinsic oscillations in metabolism, electrical activity, and Ca(2+) in pancreatic islets. Periodic variations in glucose can entrain islet Ca(2+) and insulin secretion, possibly promoting interislet synchronization. Here, we used fluorescence microscopy to demonstrate that glucose oscillations can induce distinct 1:1 and 1:2 entrainment of oscillations (one and two oscillations for each period of exogenous stimulus, respectively) in islet Ca(2+), NAD(P)H, and mitochondrial membrane potential. To our knowledge, this is the first demonstration of metabolic entrainment in islets, and we found that entrainment of metabolic oscillations requires voltage-gated Ca(2+) influx. We identified diverse patterns of 1:2 entrainment and showed that islet synchronization during entrainment involves adjustments of both oscillatory phase and period. All experimental findings could be recapitulated by our recently developed mathematical model, and simulations suggested that interislet variability in 1:2 entrainment patterns reflects differences in their glucose sensitivity. Finally, our simulations and recordings showed that a heterogeneous group of islets synchronized during 1:2 entrainment, resulting in a clear oscillatory response from the collective. In summary, we demonstrate that oscillatory glucose can induce complex modes of entrainment of metabolically driven oscillations in islets, and provide additional support for the notion that entrainment promotes interislet synchrony in the pancreas