Interaction of fast and slow dynamics in endocrine control systems with an application to β-cell dynamics

Endocrine dynamics spans a wide range of time scales, from rapid responses to physiological challenges to with slow responses that adapt the system to the demands placed on it. We outline a non-linear averaging procedure to extract the slower dynamics in a way that accounts properly for the non-linear dynamics of the faster time scale and is applicable to a hierarchy of more than two time scales, although we restrict our discussion to two scales for the sake of clarity. The procedure is exact if the slow time scale is infinitely slow (the dimensionless ε-quantity is the period of the fast time scale fluctuation times an upper bound to the slow time scale rate of change). However, even for an imperfect separation of time scales we find that this construction provides an excellent approximation for the slow-time dynamics at considerably reduced computational cost. Besides the computation advantage, the averaged equation provided a qualitative insight into the interaction of the time scales. We demonstrate the procedure and its advantages by applying the theory to the model described by Tolić et al. [I.M. Tolić, E. Mosekilde, J. Sturis, Modeling the insulin-glucose feedback system: the significance of pulsatile insulin secretion, J. Theor. Biol. 207 (2000) 361-375.] for ultradian dynamics of the glucose-insulin homeostasis feedback system, extended to include β-cell dynamics. We find that the dynamics of the β-cell mass are dependent not only on the glycemic load (amount of glucose administered to the system), but also on the way this load is applied (i.e. three meals daily versus constant infusion), effects that are lost in the inappropriate methods used by the earlier authors. Furthermore, we find that the loss of the protection against apoptosis conferred by insulin that occurs at elevated levels of insulin has a functional role in keeping the β-cell mass in check without compromising regulatory function. We also find that replenishment of β-cells from a rapidly proliferating pool of cells, as opposed to the slow turn-over which characterises fully differentiated β-cells, is essential to the prevention of type 1 diabetes.

► Dynamics on the long time scale is derived from the dynamics at the short time scale. ► The method is applied to the homeostasis of blood glucose levels and beta cell mass. ► Not just the average daily intake of glucose, but also its pattern of intake over a daily cycle influences long-term dynamics. ► Accepted rates of beta cell proliferation cannot account for rapid regeneration ► Neogenesis of beta cells is found be essential to long-term glucose homeostasis.