×

zbMATH — the first resource for mathematics

Where to look and how to look: combining global sensitivity analysis with fast/slow analysis to study multi-timescale oscillations. (English) Zbl 1425.92046
Summary: Parameterized systems of nonlinear ordinary differential equations, the type of system that is often used in mathematical models for biological systems, can be of sufficient complexity that it can take years to appreciate the full range of behaviors that can be produced. Global sensitivity analysis is one tool that has been developed for determining which parameters have the largest impact on the behavior of the model. Thus, it provides the user with a tool to know where to look in parameter space for important changes in behavior. However, it says nothing about the underlying mechanism mediating a change in behavior. For this, other tools exist. If the system dynamics occur over multiple highly-separated time scales then one useful analysis tool is fast/slow geometric analysis, also known as geometric singular perturbation analysis. This is based on bifurcation analysis of a fast or slow subsystem, and can shed light on the influence that a parameter has on structures of either subsystem, and thus on the system dynamics. Hence, once one knows where to look in parameter space for interesting behavior, this technique describes how to look at the system to extract information about how parameter changes influence the behavior of the system. In this study, we combine the two techniques in the analysis of bursting behavior in a model of insulin-secreting pancreatic \(\beta\)-cells, with the goal of determining the key parameters setting the period of the bursting oscillations, and understanding why they are so influential. This can be viewed as a case study for combining mathematical techniques to build on the strengths of each and thereby achieve a better understanding of what most influences the range of model behaviors and how this influence is brought about.
MSC:
92C30 Physiology (general)
92C42 Systems biology, networks
34D15 Singular perturbations of ordinary differential equations
Software:
LSODE
PDF BibTeX XML Cite
Full Text: DOI
References:
[1] Bertram, R.; Previte, J.; Sherman, A.; Kinard, T. A.; Satin, L. S., The phantom Burster model for pancreatic β-cells, Biophys. J., 79, 2880-2892, (2000)
[2] Bertram, R.; Rubin, J., Multi-timescale systems and fast-slow analysis, Math. Biosci., 287, 105-121, (2017) · Zbl 1377.92035
[3] Bertram, R.; Sherman, A., A calcium-based phantom bursting model for pancreatic islets, Bull. Math. Biol., 66, 1313-1344, (2004) · Zbl 1334.92089
[4] Bertram, R.; Sherman, A., Negative calcium feedback: The road from Chay-Keizer, (Coombes, S.; Bressloff, P. C., Bursting: The Genesis of Rhythm in the Nervous System, (2005), World Scientific: World Scientific Singapore), 19-48 · Zbl 1404.92043
[5] Bertram, R.; Sherman, A.; Satin, L. S., Metabolic and electrical oscillations: partners in controlling pulsatile insulin secretion, Am. J. Physiol., 293, E890-E900, (2007)
[6] Bordin, S.; Boschero, A. C.; Carneiro, E. M.; Atwater, I., Ionic mechanisms involved in the regulation of insulin secretion by muscarinic agonists, J. Membrane Biol., 148, 177-184, (1995)
[7] Butera Jr., R. J.; Rinzel, J.; Smith, J. C., Models of respiratory rhythm generation in the pre-Bötzinger complex. i. bursting pacemaker neurons, J. Neurophysiol., 82, 382-397, (1999)
[8] Canavier, C. C.; Clark, J. W.; Byrne, J. H., Simulation of the bursting activity of neuron \(r_{15}\) in aplysia: role of ionic currents, calcium balance, and modulatory transmitters, J. Neurophysiol., 66, 2107-2124, (1991)
[9] Clewley, R.; Soto-Trevino, C.; Nadim, F., Dominant ionic mechanisms explored in spiking and bursting using local low-dimensional reductions of a biophysically realistic model neuron, J. Comput. Neurosci., 26, 75-90, (2009)
[10] Cook, D. L.; Crill, W. E.; Porte, D., Glucose and acetylcholine have different effects on the plateau pacemaker of pancreatic islet cells, Diabetes, 30, 558-561, (1981)
[11] Cook, D. L.; Porte, D.; Crill, W. E., Voltage dependence of rhythmic plateau potentials of pancreatic islet cells, Am. J. Physiol., 240, E290-E296, (1981)
[12] Deister, C. A.; Teagarden, M. A.; Wilson, C. J.; Paladini, C. A., An intrinsic neuronal oscillator underlies dopaminergic neuron bursting, J. Neurosci., 29, 15888-15897, (2009)
[13] Falke, L. C.; Gillis, K. D.; Pressel, D. M.; Misler, S., Perforated patch recording allows long-term monitoring of metabolite-induced electrical activity and voltage-dependent ca \({}^{2 +}\) currents in pancreatic islet β cells, FEBS Lett., 251, 167-172, (1989)
[14] Fletcher, P.; Bertram, R.; Tabak, J., From global to local: exploring the relationship between parameters and behaviors in models of electrical excitability, J. Comput. Neurosci., 40, 331-345, (2016)
[15] Henquin, J. C.; Meissner, H. P.; Schmeer, W., Cyclic variations of glucose-induced electrical activity in pancreatic β cells, Pflügers Arch., 393, 322-327, (1982)
[16] Iooss, B.; Lemaître, P., A review on global sensitivity analysis methods, (Dellino, G.; Meloni, C., Uncertainty Management in Simulation-Optimization of Complex Systems. Uncertainty Management in Simulation-Optimization of Complex Systems, Operations Research/Computer Science Interfaces Series, 59, (2015), Springer US), 101-122
[17] Jansen, M. J.W., Analysis of variance designs for model output, Comput. Phys. Commun., 117, 1-2, 35-43, (1999) · Zbl 1015.68218
[18] Kinard, T. A.; de Vries, G.; Sherman, A.; Satin, L. S., Modulation of the bursting properties of single mouse pancreatic β-cells by artificial conductances, Biophys. J., 76, 3, 1423-1435, (1999)
[19] Larsson, O.; Kindmark, H.; Bränström, R.; Fredholm, B.; Berggren, P.-O., Oscillations in \(k_{ ATP }\) channel activity promote oscillations in cytoplasmic free ca \({}^{2 +}\) concentration in the pancreatic β-cell, Proc. Natl. Acad. Sci. USA, 93, 5161-5165, (1996)
[20] Nunemaker, C. S.; Zhang, M.; Wasserman, D. H.; McGuinness, O. P.; Powers, A. C.; Bertram, R.; Sherman, A.; Satin, L. S., Individual mice can be destinguished by the period of their islet calcium oscillations: is there an intrinsic islet period that is imprinted in vivo?, Diabetes, 54, 3517-3522, (2005)
[21] Owen, A. B., Quasi-monte carlo sampling, Monte Carlo Ray Tracing: Siggraph, 1, 69-88, (2003)
[22] Owen, A. B., Better estimation of small Sobol sensitivity indices, ACM T. Model. Comput. S., 23, 2, 11, (2013) · Zbl 06860009
[23] Radhakrishnan, K.; Hindmarsh, A. C., Description and use of LSODE, the livermore solver for ordinary differential equations, Lawrence Livermore Natl. Labor. Rep., (1993)
[24] Rinzel, J.; Ermentrout, G. B., Analysis of neural excitability and oscillations, (Koch, C.; Segev, I., Methods in Neuronal Modeling: From Synapses to Networks, (1989), MIT Press: MIT Press Cambridge), 135-169
[25] Saltelli, A.; Annoni, P.; Azzini, I.; Campolongo, F.; Ratto, M.; Tarantola, S., Variance based sensitivity analysis of model output. design and estimator for the total sensitivity index, Comput. Phys. Commun., 181, 2, 259-270, (2010) · Zbl 1219.93116
[26] Smith, P. A.; Ashcroft, F. M.; Rorsman, P., Simultaneous recordings of glucose dependent electrical activity and ATP-regulated k \({}^+\)-currents in isolated mouse pancreatic β-cells, FEBS Lett., 261, 187-190, (1990)
[27] Sobol’, I. M., On the distribution of points in a cube and the approximate evaluation of integrals, Zhurnal Vychislitel’noi Matematiki i Matematicheskoi Fiziki, 7, 4, 784-802, (1967)
[28] Sobol’, I. M., Global sensitivity indices for nonlinear mathematical models and their monte carlo estimates, Math. Comput. Simul., 55, 271-280, (2001) · Zbl 1005.65004
This reference list is based on information provided by the publisher or from digital mathematics libraries. Its items are heuristically matched to zbMATH identifiers and may contain data conversion errors. It attempts to reflect the references listed in the original paper as accurately as possible without claiming the completeness or perfect precision of the matching.