Next Article in Journal
Designing the Architecture of a Convolutional Neural Network Automatically for Diabetic Retinopathy Diagnosis
Next Article in Special Issue
Global Stability of a MERS-CoV Infection Model with CTL Immune Response and Intracellular Delay
Previous Article in Journal
NARX Deep Convolutional Fuzzy System for Modelling Nonlinear Dynamic Processes
Previous Article in Special Issue
Global Stability of a Reaction–Diffusion Malaria/COVID-19 Coinfection Dynamics Model
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Mathematical Modelling of Combined Intervention Strategies for the Management and Control of Plasma Glucose of a Diabetes Mellitus Patient: A System Dynamic Modelling Approach

by
Vincent O. Omwenga
1,*,
Vaishnav Madhumati
2,
Kumar Vinay
3,
Sathyanarayan Srikanta
2 and
Navakanta Bhat
2
1
School of Computing and Engineering Science, Strathmore University, Nairobi 59857-00200, Kenya
2
Samatvam Endocrinology Diabetes Centre, Bangalore 560041, India
3
Center for Nano Science and Engineering, Indian Institute of Science, Bangalore 560012, India
*
Author to whom correspondence should be addressed.
Mathematics 2023, 11(2), 306; https://doi.org/10.3390/math11020306
Submission received: 23 October 2022 / Revised: 26 December 2022 / Accepted: 30 December 2022 / Published: 6 January 2023

Abstract

:
With the rapid increase of diabetes mellitus cases in the world, management and control of the disease has become a complex and highly dynamic process. This challenge requires a multifaceted approach to manage and control the complications associated with the hyperglycaemia or hypoglycaemia conditions. This paper presents a mathematical model for determining the influence of combined intervention strategies in the management and control for the plasma glucose of the type II diabetes. System dynamics (SD) techniques were used in modelling the sub-compartments of biological systems of an Identifiable Patient (IP). The system dynamic model developed gave an illustration on how typical glucose-insulin dynamics occur at different intervention strategies involving varying amounts of carbohydrates taken, intensity of physical exercises, stress levels and the amount of exogenous insulin administered. The model was conceptualized within a semi-closed loop system representing the patient ecosystem by extending the Bergman Minimal Model. Stochastic differential equations (SDE) were used to capture the non-linear, continuous time varying interactions of the measurements associated with plasma glucose-insulin dynamics. The estimated results from the model showed combined intervention strategies of reduced amounts of carbohydrates intake, reduced stress levels and varying moderately high-to-low exercise intensity at a constant unit of exogenous insulin produced good plasma glucose levels control.

1. Introduction

Diabetes, commonly referred to as diabetes mellitus, is one of the rapidly growing health concerns in the world [1,2]. It is considered to be a metabolic disorder disease caused by either the pancreas not being able to produce insulin which leads to type I diabetes or by insulin resistance developing near the body tissues on the glucose-regulatory action of the insulin which leads to type II diabetes. The latter is the most common accounting for almost 80% of the reported cases world-wide. In both type I and type II, the plasma glucose levels are not within the normal range. The amount of plasma glucose is controlled by a glucose homeostasis mechanism regulated by two hormones, namely insulin and glucagon. Insulin is secreted by the β -cells while glucagon is secreted by the -cells found in the Langerhans of the human pancreas [2].
For the body to function optimally, it requires optimal plasma glucose levels. The optimal plasma glucose levels fall within a very tight range of between 70 mg/dL ((3.90 mmol/L) and 110 mg/dL (7.00 mmol/L) [3]. Maintaining this ranges is essential in the management and control of diabetes. It also aids in delaying the onset of long-term complications caused by hyperglycaemia (>120 mg/dL (7.00 mmol/L)) or hypoglycaemia (<70 mg/dL (3.90 mmol/L). In the case of hyperglycaemia, the glucose molecules diffuse into the β -cells leading to an increase of the ratio between adenosine triphosphate (ATP) and the adenosine diphosphate (ADP) inside the β -cells that will cause an increased release of insulin [4]. Insulin resistance and an inadequate secretion of insulin hormone among diabetics are considered to be causal factors of hyperglycaemia. Insulin resistance occurs when the β -cells in the liver, muscles and adipose cells become less sensitive to the insulin hormones [5]. Due to this, any consumption of food rich in carbohydrates causes elevated levels of plasma glucose. The interaction between the insulin released and the glucose molecules leads to usage and storage of excessive plasma glucose in the liver, adipocyte tissues and skeletal muscles [5].
On the other hand, glucagon resistance leads to hypoglycaemia. Glucagon occurs when there is a limited response to the glycogenolysis and gluconeogenesis actions by the glucagon receptors in the -cells. The risk factors associated with hypoglycaemia includes kidney diseases, prolonged fasting and excess endogenous insulin secretion [6]. Other causes are insulin overdose and drugs used to treat diabetes mellitus [7]. When hypoglycaemia occurs, the β -cells are considered to be inactive causing an elevation of ATP to ADP ratio inside the -cells which leads to the release of glucagon hormone that stimulates glucagon receptors in the kidney, liver and adipose tissues [4]. The role of the glucagon receptors is to maintain the plasma glucose homeostasis by regulating the glycogenolysis and gluconeogenesis processes [8]. A failure in this regulatory process leads to the plasma glucose levels that are outside the normal range leading to complications experienced by diabetic patients.
It is also noted that there is a link between stress and type II diabetes mellitus [9]. This could be explained by the fact that some reactions by human beings to a stressor include developing unhealthy lifestyles—for instance, eating disorders. Moreover, the physiological changes that occur due to stress may affect the neuroendocrine and immune systems [10]. In the neuroendocrine system, a stressor causes the hypothalamus in the hypothalamo–pituitary–adrenal (HPA) axis to secrete corticotropin-releasing factor, which prompts the release of adrenocorticotropic. The adrenocorticotropic causes the adrenal cortex to secrete glucocorticoid hormones—specifically, the cortisol. The cortisol is considered to be one of the factors influencing the glucose metabolism and other metabolism processes including the breakdown of lipids, carbohydrates and proteins. It also leads to a rise of the plasma glucose levels by stimulating the hepatic gluconeogenesis while inhibiting the insulin actions [11]. The glucose released by the hepatic gluconeogenesis process, if not used, remains in the bloodstream leading to a rise in plasma glucose.
Therefore, a constant monitoring of the plasma glucose levels using glucometers by diabetes patients is needed. This approach is commonly used to ensure corrective interventions if any are instituted promptly whenever a significant deviation from the acceptable plasma glucose ranges is noticed. Some of the interventions include oral medication and insulin therapy, which involves the subcutaneous insulin injection at prescribed dossal units at regular time interval. For the accepted target to be realized, a strict adherence to monitoring and the institution of prescribed interventions are necessary. However, the constant monitoring of the plasma may sometimes fail due to factors beyond the control of the patient leading to complications such as retinopathy, neuropathy, comatose and dizziness. To overcome these challenges, scientist have developed mathematical and computer simulation models for the prediction of plasma glucose dynamics and trajectories based on the characteristics of insulin secretion and glucose utilization. Most of these models are designed to simulate the glucose-insulin dynamics in vivo [12]. A number of these models are based on Ordinary Differential Equations [13,14], Delay Differential Equations [15], Partial Differential Equations [16] and Stochastic Differential Equations (SDE) [17]. The underlying mathematical foundations of the developed system dynamics model is the SDE. SDE has been adopted due to its ability to capture non-linear, continuous time measurements in a more reliable manner, and its robustness in using maximum likelihood estimation techniques [17]. SDE is also able to model the random noise associated with system variables due to inherent uncertainties through the Weiner process term which models process diffusion.
A better prediction of the plasma glucose level trajectories using mathematical models requires proper understanding of the causal factors. Furthermore, modelling and determination of the glucose-insulin actions are complex and require a complete understanding of the glycaemic profile of an individual including the glycaemic disturbances and patient behaviour [18]. Moreover, the typical glucose-insulin models contain sub-compartments including the glucose model, insulin deliver model, gastro-intestinal tract, muscle and adipose tissue and the liver model; but a complete human glucose-insulin is influenced by many more biological systems and processes. Some of the factors that influence plasma glucose and, hence, the biological systems include the amount of carbohydrates ingested, the level of physical exercises, the physiological state of the person (stress), the body weight (BW) of the person, the genetic disposition of the individual, the metabolic characteristics of the individual, their familial history, the gestational diabetes, their life style, their knowledge of the diabetes and the access/availability of health services, among other factors. Therefore, understanding individual information such as risk factors is key in the development of glucose-insulin simulators since they provide an effective means of identifying diabetes complications [19].
In this study, the focus was on modelling the combined effects of the amount of carbohydrates ingested, the intensity of the physical exercise, the stress levels and the insulin therapy as the intervention strategies on the management and control of glucose dynamics of a diabetes population stock. Through this study, a clear understanding on the contribution in the management of diabetes is demonstrated. The study is conceptualized on a diagnosed diabetic state of the population stock assuming a dynamic state of a diabetic patient. This is because diagnosis is vital in the management and control of the risks associated with extreme glucose levels, both hypoglycaemia and hyperglycaemia. Moreover, diagnosis will determine the management interventions to be put in place.

2. Materials and Methods

A System Dynamics Modelling (SDM) has been adopted to build the model for determining the effect of the diabetes management and control intervention strategies involving meals, physical exercise, stress and insulin administration. The SDM is considered to be reliable in modelling complex health systems [20]. Moreover, it provides a mechanism of understanding the dynamical changes that occur in a complex situation, such as health, with a focus on analysing the cause-and-effect relations [21]. It also offers a whole-system-thinking approach that provides an understanding on how the interplaying factors interact and influence each other.
In SDM, causal and loop diagrams (CLDs) are used to demonstrate the interaction between sub-systems. CLDs are also used in analysing the feedback effects among the sub-systems and to show the interconnection between the state variables following the Markovian-linked chain paths approach. Furthermore, a Reinforcing-loop type of CLD is used to show the directional effect (the stocks and flows) of the head variables. In this study, meals, physical exercises and insulin administration are modelled as substantive system states while stress is modelled as a pseudo-system state since it can occur co-jointly with the other sub-systems using CLD.
Since the management and control of diabetes is a time-dependent process that occurs in a continuous manner, the SDM model is developed to capture the variables’ changes within sub-systems over a time period using Stochastic Differential Equations (SDEs). A random noise seed is also introduced to handle the time-dependent perturbations that will accompany the variables. The stocks and flows are measured per time units [20,22]. Therefore, the general SDE that describes the glucose-insulin dynamics involving the stocks and flow rates (inflows and outflows) of plasma glucose and insulin in the blood system is represented by the following:
dx t = f x t ,   u t ,   t ,   θ dt + σ u t ,   t ,   θ d ω
y i = g x i , u i , t ,   θ + ε i
where x t is a continuous time-dependent system that describes the glucose-insulin dynamics, f x t ,   u t ,   t ,   θ is drift component, σ u t ,   t ,   θ d ω is the diffusion component, y i is an observation for the glucose-insulin dynamics in a state given u i inputs, θ is the system parameters and ε i is the random noise term.
Equation (1) captures the system perturbations as a result of the diffusion term of the model. In this study, we modelled the physiological perturbations that occurs as stochastic behaviours due to meals taken, physical exercises, stress levels and insulin injection as part of the system parameters. The error term in the system observation equation provides a measure for uncorrelated system errors that are likely to occur in the case of randomly occurring events.
The model is built on certain assumptions about the future based on the current intervening variables. For purposes of modelling in this study, the interventions are assumed to exhibit current characteristics in the near future. This assumption only offers a starting point for future predictions and policy analysis but can be varied. Other assumptions considered in this study include the following: insulin is infused exogenously into the system, and, therefore, the study considered an insulin-dependent patient; ingested food has a direct effect on the plasma glucose concentrations; physical exercise intensity influences the plasma glucose concentrations and decreases linearly after attaining a certain maximum value; and plasma glucose levels are influenced by the stress severity levels which is intertwined with other state variables. Applying these assumptions and borrowing from [12] on the Bergman Minimal Model, the glucose dynamics as a consequence of carbohydrates ingested, physical exercises intensity, stress level and exogenously infused insulin, the sub-systems interactions are mathematically presented as follows:
dI si dt = 1 ρ 1 I in I cl I si dt + σ I si d ω 1
dI pi dt = 1 ρ 2 ( 1 + s I si I pi I ir dt + σ I pi d ω 2
where I si is the subcutaneous insulin concentration mU / L , I pi is the plasma insulin concentration mU / L , I in is the exogenous insulin infusion rate mU / min to the subcutaneous, s is the stress levels, I ir is the rate of insulin removal from the circulatory system as a result of exercise-induced physiological changes mU / L / min , σ is the scaling parameter for the diffusion term, d ω i represents the Wiener process which has Gaussian distribution characteristics and I cl ,   ρ i   and   ρ 2 are other system parameters defined in Table 1.
Equations (3) and (4) describe the plasma insulin and subcutaneous insulin pharmacokinetics. Equation (4) is modified from the original Bergman Minimal Model to include the terms for measuring the effects of exercise and stress on insulin levels in the plasma as a result of exercise and/or stress-induced physiological changes. It is also recognized that physical exercise improves insulin sensitivity [23] and stress decreases plasma insulin [2].
Figure 1 depicts the diagrammatic representation of the flow of glucose and insulin with associated interactions that take place in the extended Bergman Minimal Model compartments.
From Figure 1, dynamics interactions between insulin and glucose within the various subsystems are represented mathematically as follows:
dI e dt = T dr S s . I pi I e dt + σ I e d ω 3
dG p dt = I e . G p + EG pr 1 + s + CA I 2 ρ m + w vol gd G pd G gp G u dt + σ G p d ω 4
where I e measures the effect of insulin on glucose transportation min 1 ; S s is the measure of insulin sensitivity 1 / mU × min ; G p is the amount of plasma glucose mmol / L ; EG pr is the rate of production of the endogenous glucose when the patient has stress or no stress mg / L × min ; CA I 2 is the amount of carbohydrates absorption delayed in a compartment (mg); s is the factor for the stress levels; w is the weight (kg) of the patient, G pd is the rate of hepatic glucose production due to meal, exercise and stress; G gp represents the rate of glycogenolysis decline due to pro-longed exercise; and G u is the rate of glucose uptake during exercise.   T dr , vol gd and ρ m are other system parameters described in Table 1.
The effect of meals on the plasma glucose is dependent on the amount of carbohydrates ingested, the glucose distribution volume and the constant associated with carbohydrates movement. This effect is, therefore, mathematically represented in a two-stage compartmental model [24]:
dGA I 1 dt = Ca b . Ca I vol gd GA I 1 ρ m dt + σ GI 1 d ω 6
dGA I 2 dt = GA I 1 ρ m GA I 2 ρ m dt + σ GI 2 d ω 7
where GA I 1 represents the amount of carbohydrates absorbed (mg), GA I 2 is the amount of carbohydrates absorption delayed in a compartment (mg), Ca b is the bioavailability associated with the carbohydrates ingested and Ca I is the amount of ingested carbohydrate per minute g / min .
According to a study by [24], the rates of hepatic glucose production G pd , glucose uptake G u and plasma insulin removal I r as a result of physical exercise is represented by the following:
dG pd dt = r 1 . E cmax r 2 . G pd dt + σ G pd d ω 8
dG u dt = r 3 . E cmax r 4 . G u dt + σ G u d ω 9
dI r dt = r 5 . E cmax r 6 . I r dt + σ I r d ω 12
dE cmax dt = 0.8 . E in 0.8 . E cmax dt + σ E cmax d ω 11
where E cmax is a measure of the maximum exercise capacity level. Other parameters are described in Table 1.
It is also noted that the rate of glycogenolysis ( G gl t ), which influences glucose-insulin dynamics as a result of physical exercise, starts to decline at a certain maximum value G glmax and is dependent on exercise intensity and duration of the exercise [25]. This can be represented mathematically as follows:
G glmax = E t in T gl E t in
where T gly is the duration of exercise performed at exercise intensity ( E t in ) before the start of decrease in hepatic glycogenolysis [7,25].
The exercise intensity is also expected to decline over time following a linear fit [25] represented by the following:
E t in = 87.471 T gl 1.152
Therefore, substituting Equation (14) in Equation (13) gives the maximum value upon which glycogenolysis is likely to start to decrease as follows:
G glmax = 1.1521 E t in 2 + 87.471   E t in  
We, therefore, model G gl t as follows:
dG gl dt = r g G glmax C t dt + σ G gl d ω 10
where r g and C t are parameters described in Table 1.
Table 1. Parameters of the SDE model.
Table 1. Parameters of the SDE model.
ParameterUnitDescriptionRangeReferences
I cl L / min Insulin clearance from the subcutaneous rate0.54–2.01[26]
ρ 1 min Time in insulin movement action due to an increase in subcutaneous40–131[26]
ρ 2 min Time in insulin movement action due to an increase in plasma insulin10–70[26]
T dr min 1 Time delay ratio in action of insulin due to an increase of insulin in plasma8.14 × 10−3–2.23 × 10−2[26]
ρ m min Maximum time for carbohydrate absorptionFixed at 107[24]
r 1 mL / kg × min Rate of oxygen consumption during glucose production1.3 × 10−3–1.9 × 10−3[27]
r 2 min 1 Glucose production rate4.41 × 10−2–6.7 × 10−2[27]
r 3 mL / kg × min Rate of oxygen consumption during glucose uptake1.5 × 10−3–2.4 × 10−3[27]
r 4 min 1 Glucose uptake rate3.55 × 10−2–6.17 × 10−2[27]
r 5 mU / mL × min Insulin clearance rate due to exercise1.0 × 10−3–1.5 × 10−3[27]
r 6 min 1 Rate of insulin removal from the circulatory system5.88 × 10−2–9.12 × 10−2[27]
r g mL / kg × min Rate of change of glycogenolysis 8.5 × 10−3–1.31 × 10−2[27]
C t min Time constant required for glycogenolysis to return to basal value1.86–10.14[27]
S S L / mU × min Insulin sensitivity 9.64 × 10−5–1.73 × 10−3[26]
EG pr mg / dl × min Rate of glucose production6.0 × 10−1–3.45[26]
Vol gd dl/kgMean volume of glucose distribution1.93–4.14[26]
Ca b DimensionlessCarbohydrates bioavailabilityFixed at 0.9[26]
Equations (3)–(16) represent the plasma glucose-insulin dynamical interactions within the various subsystems shown in Figure 1. The equations are used to estimate the observable stock flows from and into the subsystems to predict the plasma glucose levels. The observation equation measuring the aggregate observations of the stock flows with an error term to account for any variabilities due to uncertainties and system imperfections is given by the following:
y G p = G p + exp e G p ,   exp e G p ϵ N 0 ,   S G p
where e G p represents the errors measurements which are assumed to be identically independently distributed (iid), S G p denotes variance measurement of the noise associated with the plasma glucose measurements.
To validate the mathematical models of the stocks and flows of the plasma glucose and insulin dynamics in the blood system, a SDM was developed using the Vensim®PLE software version 7.2a, Ventana Systems Inc, Harvard, MA, USA, 2018. Figure 2 shows the SDM of the causal interaction structure of the system dynamics for the plasma glucose control and diabetes management strategies. It depicts the stock inflows and outflows interactions in a semi-closed system loop involving the meals, exercise, stress and exogenous insulin. In the SDM model, the stocks that represent subsystems of the glucose-insulin dynamic system are identified in boxes and the flows that represent the input parameters are identified by arrows (pipes).
Equation (6) was used to test the influence of exogenous insulin, meal, physical exercise and stress on the plasma glucose. The outcome of Equation (6) is identified as a stock of labelled plasma glucose levels in Figure 2. The input parameter components of the Equation (6), which included exogenous insulin, stress, meals and physical exercises, were estimated using the various model equations implemented in the simulated SDM as follows: the exogenous insulin input parameters I e . G p were computed using Equation (5) that gave the rate of insulin uptake dI e dt and the amount of plasma glucose G p recorded which is used to measure the amount of insulin to be administered; the stress input parameters EG pr 1 + s were estimated as a scalar factor for the reason that stress was implemented as a pseudo sub-system in the SDM due to its direct influence on the rate of production of exogenous glucose and the rate of plasma insulin change occurring at the same time [28]; the input parameters for meals CA I 2 ρ m were estimated using Equation (7) and Equation (8) whose output is captured as plasma glucose due to carbs in the SDM representation in Figure 2, while the input parameters of physical exercise w vol gd G pd G gp G u were computed using Equations (9)–(16) with the influence of physical exercise on plasma glucose captured as glucose due to exercise in Figure 2. The scaling parameter for the diffusion term σ j and the Wiener process factor d ω i were used to give Equation (6) the Gaussian distribution characteristics. The result of Equation (6) which formed the input parameter G p of Equation (17) was used to compute the predicted plasma glucose value y G p while taking into consideration the uncertainties and system imperfections estimated as exp e G p in the equation.
To compare the performance of the proposed model, different values of the input parameters in the various equations and other parameters were used to generate simulated results using the SDM in the Vensim®PLE software version. The computations on the various model equations were to give the feedback, accumulation of flows into stocks and the time delays of the SDM. The final results of the predicted plasma glucose values obtained through Equation (7) at the end of the interactions of the model equations were compared with those recorded by the CGM to validate the model. The outcome of the dynamic interactions of the feedback, accumulation of flows into stocks and the time delays in the model provided an indication of the current system state’s behaviour when current input parameter values are entered into the simulated model. Using the generated current state behaviours, future predictions can be made by varying the model parameter values.

3. Results and Discussions

Validation of the mathematical models of the plasma glucose-insulin dynamics, as depicted in Figure 1, followed a series of tests which included analysis on both the structural and behavioural patterns. The structural tests were conducted to ascertain the real-world dynamics of the system interactions and to check if the theoretical foundations of the sub-systems are correctly mapped into the models. The tests involved confirmation of the interaction effects of the various model parameters on the anticipated outcomes. It, therefore, required the verification of both the parameters and the structures that were used in the model by comparing them with the existing literature [1,29]. The behavioural validation sought to examine the characteristics of the behavioural patterns of the outputs from the model. The numerically simulated behavioural patterns were compared with the expected, real-world data’s behavioural patterns.
An analysis to evaluate the behavioural performance of the model was conducted by examining the Autocorrelation function (ACF) for the standardized residuals and the one-step predictions using the maximum likelihood ratio tests. ACF was used to examine the independence of observations. Figure 3 and Figure 4 illustrate the ACF for the plasma glucose residuals recorded in the first 24 h of the experiment for the 24-hour plasma glucose residuals and the 4-day consecutive plasma glucose residual measurements, respectively.
The two ACFs show that the y G p residuals are uncorrelated and are, therefore, considered independent at a 95% confidence level for the observations recorded on the type II diabetes patient. The one-step prediction of the plasma glucose outputs using the proposed model are presented in Figure 5 and Figure 6 for the first 24 h and the 4 consecutive days’ observations, respectively.
From Figure 5 and Figure 6, the predicted values are within the acceptable range in the first 24 h of the experiment and in the four (4) consecutive days of the experiment. However, the sensitivity of the model seems to improve with more data as shown in Figure 6.
After confirming the behavioural performance of the model, different simulation runs were conducted on a Vensim®PLE simulation environment to further analyse the predictive ability of the model to estimate the overall plasma glucose given different diabetes management strategies. Table 2 presents a selected sample of outputs from the base run (the initial conditions) and simulation runs using clinical data from an IP at varying intervals of diabetes management and control intervention strategies. In this study, the base run refers to measurements of the glucose dynamics as recoded by the Continuous Glucose Monitoring (CGM) device while the diabetes management and control innervations refer to strategies involving carbohydrates intake, physical exercise intensity, stress levels and exogenously administered insulin. The CGM was done using Abbott FreeStyle® Libre Pro device (also referred to as Flash Glucose Monitoring (FGM)/Ambulatory Glucose Monitoring (AGM) System), Abbot, Chicago, IL, USA, 2016. The CGM sensor measured interstitial (in between the cells) plasma glucose and provided patterns and trends data that assist in making clinical decisions. The components of Abbott FreeStyle® Libre Pro device include a sensor that is inserted just under the skin into the subcutaneous tissue and lasts for 14 days (and gives calculated blood glucose readings every 15 min). The World Health Organization Well-being (1998 Version) tool was used to assess the stress levels of the patients during the study period.
A 24 h base run observation of the plasma glucose of the type II IP using the CGM was conducted. The recorded data gave an average of 9.22 mmol/L as shown in Table 2. This was considered above the set targets for the patient but was used as the initial value for predicting the potential plasma glucose trends. Figure 7 shows the plasma glucose profile levels for the patient generated using the proposed model considering the initial conditions (base run values).
The base run results were used to calibrate and check on the performance of the mathematical model as implemented using the SDM in Figure 2 to estimate the plasma glucose levels at different input parameters. Using the base run input parameter values, different simulation runs involving different intervention strategies in the management and control of diabetes were conducted. Table 2 shows the experimental estimated values of the average plasma glucose levels at varying input parameter values and different diabetes management and control intervention strategies as generated by the model.
Table 2 gives the input parameter values for the different strategies for the management and control of diabetes. Using the input parameter values, different plasma glucose trends using the proposed mathematical model were generated. These trends are discussed in case 1 to case 5.
  • Case 1: Increasing carbohydrates and stress levels while holding the exercise at base run levels
The findings from the simulation runs analysis show the amount of carbohydrates taken and stress levels are closely associated with the plasma glucose levels’ dynamics. It was observed that a 33.33% decrease in the amount of carbohydrates ingested while maintaining the base run levels of stress and exercise intensity led to a decrease of glucose level by 11.71%. This outcome is supported by previous studies that showed the amount of meal and, therefore, portion of carbohydrates ingested will have corresponding impact on the plasma glucose [12,30]. Figure 8 shows the plasma glucose profile of the simulated results from the intervention strategy involving the reduction of the amount of carbohydrates consumed below the base run values.
  • Case 2: Lowering the stress level while maintaining carbs and exercise intensity at base run levels
Similarly, a reduction in the stress level by 33.33%, while holding other intervention strategies at base run levels, led to a decrease in the glucose level by 2.93%. This finding affirms previous studies on brain glucose homeostasis and associated effects on plasma glucose levels. It underscores previous knowledge on the impact of stress on the metabolic processes [28]. For instance, stress is associated with the energy mobilization trigger for the fight-or-flight responses. This, consequently, influences the stimulation of glucocorticoid hormones which leads to increased plasma glucose concentrations in the blood stream system. It is also noted that high stress levels will have an effect on pancreatic glucagon production as a result of the influence of high stress on catecholamine [28]. An increase in catecholamine could lead to an increased glycogenolysis and gluconeogenesis in the liver. During this time, the insulin levels are considered to drop thereby causing an increase in hepatic glucose production [28]. Figure 9 shows the simulated results profile of the plasma glucose levels when stress level is reduced while maintaining the amount of carbohydrates ingested and the intensity of the physical exercise at the base run levels.
  • Case 3: Lowering carbohydrates, while maintaining the stress level and the exercise intensity at the base run levels
A change of the exercise intensity from moderately high to moderately low and/or near zero while holding both the carbohydrates and stress levels at the base run levels led to a decrease in plasma glucose levels by 0.33%. The minimal decline in the plasma glucose due to physical exercise is supported by previous studies that notes an increased hepatic glucose release due to enhanced glycogenolysis during the physical exercise period. Glycogenolysis activities are aimed at regulating the plasma glucose homeostasis. The effect of the physical exercise is also influenced by the oxygen intake rate and the duration of the exercise upon which a certain intensity is maintained. This, therefore, means that physical exercise is dependent on other factors which needs further investigations. Figure 10 shows the result of lowering exercise intensity while maintaining the amount of carbohydrates and stress levels at the base run levels. From the graph, episodes of swing of plasma glucose from high levels to low levels are evident. This illustrates poor attainment of control of the plasma glucose.
  • Case 4: Lowering carbohydrates and the stress levels, while maintaining exercise at base run levels
A combined reduction of amounts of carbohydrates intake and stress levels each at 33.33% while maintaining the exercise intensity at the base run levels led to a 14.75% reduction in the plasma glucose levels. This strategy had a relatively lower variability in the recorded plasma glucose levels over the time intervals. This implies a moderately good glycaemia control is achieved. Figure 11 shows the trend in the plasma glucose with a reduced amount of carbohydrates and stress levels while maintaining high intensity exercise. The trend demonstrates a highly controlled plasma glucose profile. From the graph, the variability of the observed values is small, indicating a good control of the plasma glucose-insulin.
  • Case 5: Increasing carbohydrates and stress levels, while maintaining the exercise intensity at the base run level
An increase in both the amount of carbohydrates intake and the stress level by 66.67% and 233.33%, respectively, while holding the exercise intensity at the base run level resulted in an increase of 30.59% in the plasma glucose levels. This scenario also had the largest variability in observed plasma glucose values. This means that the glucose management and control could be difficult as the glucose levels swing from one extreme values to the next. This, therefore, means that the corrective intervention will be both to regulate carbohydrates intake and to manage stress levels. Figure 12 shows the simulated outcome of the increased amount of carbohydrates and stress level while holding the physical exercise at the base run.
From the simulated experimental runs, the best plasma glucose control levels for the type II IP patient are estimated to be attained when both the carbohydrates portion and stress levels are reduced below the base run values while maintaining moderately high to low physical exercises intensity as illustrated in Figure 11. This combination of interventions strategies is estimated to give an average of 7.86 mmol/L which will be within the set target for the IP patient. It also projects the achievement of the least variability in the recorded observations of 1.827. This indicates that the patient will have better control of the plasma glucose levels.

4. Conclusions

Whereas the consequences of insulin resistance or the damage of β -cells due to plasma glucose concentration is known, a complete understanding of the root cause of hyperglycaemia or hypoglycaemia in diabetes patients is vital for overall management and control of the plasma glucose. A numerical analysis of the relationship between plasma glucose-insulin dynamics and the most common diabetes management and control intervention strategies is modelled in this study.
The study tries to expand the knowledge about the effect of carbohydrates, exercise and stress as causes of plasma glucose changes in diabetes patient. It also offers an understanding on how to select the best combined intervention strategies for good management and control of the plasma glucose. The parameters for plasma glucose changes were formulated and simulated. The parameters were used to analyse the influence of carbohydrates, exercise intensity and stress level as factors that lead to plasma glucose changes. To understand the influence of each factor on plasma glucose, sensitivity analysis of the parameters was conducted. Even though plasma glucose was observed to be sensitive to the insulin uptake rate, the sensitivity indices of other parameters were also observed to cause different plasma glucose concentrations. The parameter combinations with sensitivity indices leading to expected plasma glucose targets for a patient were considered to be the best management and control strategies.
In this study, many input parameters representing different conditions are considered. Some of the input parameters’ values are, however, difficult to measure in practical terms. This may limit the applicability of the proposed model in real-life situations—for example, predicting hyperglycaemia or hypoglycaemia in a patient’s blood system. We, therefore, recommend more investigation on the relationships between the parameters themselves which are considered in each factor. Analysis on the performance of the proposed model was limited to a small set of data and to type II diabetes, more in vivo and in vitro studies are needed to examine the contribution of different parameters in each condition and to determine the effect of each input parameter value on the overall plasma glucose.

Author Contributions

All the authors individually contributed as follows: conceptualization, investigation, data curation and formal analysis, V.O.O.; methodology, data curation and investigation, V.M. and S.S.; validation, K.V.; writing—review and editing, V.O.O.; methodology, investigation guiding and reviewing, N.B. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Institutional Review Board Statement

The study was conducted according to the guidelines of the Declaration of Helsinki, and approved by the Ethics Committee of Samatvam Endocrinology Diabetes Center, Bangalore, India.

Informed Consent Statement

Informed consent was obtained from all subjects involved in the study. Written informed consent has been obtained from the patient(s) to publish this paper.

Data Availability Statement

The data presented in this study are available on request from the corresponding author. The data are not publicly available due to privacy and ethical restrictions on the use of the data.

Acknowledgments

The authors acknowledge Strathmore University, Kenya; Department of Science & Technology; the Government of India and FICCI for the postdoctoral support under the C.V. Raman International Fellowship for African Researchers—2016; Centre for Nano Science and Engineering (CeNSE) at the Indian Institute of Science and the Samatvam Endocrinology Diabetes Centre, Bangaluru, India for their support.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Jones, A.P.; Homer, J.B.; Murphy, D.L.; Essien, J.D.K.; Milstein, B.; Seville, D.A. Understanding Diabetes Population Dynamics through Simulation Modelling and Experimentation. Am. J. Public Health 2006, 96, 488–494. [Google Scholar] [CrossRef] [PubMed]
  2. Mwita, P.S.; Shaban, N.; Mbalawata, I.S.; Mayige, M. Mathematical modelling of root causes of hyperglycemia and hypoglycemia in a diabetes mellitus patient. Sci. Afr. 2021, 14, e01042. [Google Scholar] [CrossRef]
  3. Felig, P.; Wahren, J. Fuel Homeostasis in Exercise. N. Engl. J. Med. 1975, 293, 1078–1084. [Google Scholar] [CrossRef] [PubMed]
  4. Gaisano, H.; MacDonald, P.E.; Vranic, M. Glucagon secretion and signalling in the development of diabetes. Front. Physiol. 2012, 3, 349. [Google Scholar] [CrossRef] [Green Version]
  5. Petersen, M.C.; Shulman, G.I. Mechanisms of insulin action and insulin resistance. Physiol. Rev. 2018, 98, 2133–2223. [Google Scholar] [CrossRef] [Green Version]
  6. Cryer, P.E.; Davis, S.N.; Shamoon, H. Hypoglycemia in diabetes. Diabetes Care 2003, 26, 1902–1912. [Google Scholar] [CrossRef] [Green Version]
  7. Pruett, E.D. Glucose and Insulin during prolonged work stress in men living on different diets. J. Appl. Physiol. 1970, 28, 199–208. [Google Scholar] [CrossRef]
  8. Hughes, D.S.; Narendran, P. Alpha cell function in type 1 diabetes. Br. J. Diabetes 2014, 14, 45–51. [Google Scholar] [CrossRef] [Green Version]
  9. Li, J.; Jarczok, M.N.; Loerbroks, A.; Schöllgen, I.; Siegrist, J.; Bosch, J.A.; Wilson, M.G.; Mauss, D.; Fischer, J.E. Work stress is associated with diabetes and prediabetes: Cross-sectional results from the MIPH Industrial Cohort Studies. Int. J. Behav. Med. 2013, 20, 495–503. [Google Scholar] [CrossRef] [Green Version]
  10. Falco, G.; Pirro, P.S.; Castellano, E.; Anfossi, M.; Borretta, G.; Gianotti, L. The Relationship between Stress and Diabetes Mellitus. J. Neurol. Psychol. 2015, 3, 7. [Google Scholar]
  11. Bergmann, N.; Gyntelberg, F.; Faber, J. The appraisal of chronic stress and the development of the metabolic syndrome: A systematic review of prospective cohort studies. Endocr. Connect. 2014, 3, R55–R80. [Google Scholar] [CrossRef]
  12. Adriana, A.; Holger, V.; Mohamed, D. Glucose-Insulin Minimal Model Augmented: Identification and Validation. In Proceedings of the International Conference on Modelling and Simulation, 25–28 November 2015. [Google Scholar]
  13. Bergman, R.; Phillips, L.; Cobelli, C. Physiologic evaluation of factors controlling glucose tolerance in man. J. Clin. Investig. 1981, 68, 1456–1467. [Google Scholar] [CrossRef] [Green Version]
  14. Parker, R.; Doyle, F.; Peppas, N. A Model-based algorithm for blood glucose control in type I diabetic patients. IEEE Trans. Biomed. Eng. 1999, 46, 148–157. [Google Scholar] [CrossRef]
  15. Bennett, D.; Gourley, S. Asymptotic properties of a delay differential equation model for the interaction of glucose with plasma and interstitial insulin. Appl. Math. Comput. 2004, 151, 189–207. [Google Scholar] [CrossRef]
  16. Keener, J. Diffusion induced oscillatory insulin secretion. Bull. Math. Biol. 2001, 63, 625–641. [Google Scholar] [CrossRef] [Green Version]
  17. De Vries, G.; Sherman, A. Channel sharing in pancreatic-cells revisited: Enhancement of emergent bursting by noise. J. Theor. Biol. 2000, 207, 513–530. [Google Scholar] [CrossRef] [Green Version]
  18. Mujahid, O.; Contreras, I.; Beneyto, A.; Conget, I.; Giménez, M.; Vehi, J. Conditional Synthesis of Blood Glucose Profiles for T1D Patients Using Deep Generative Models. Mathematics 2022, 10, 3741. [Google Scholar] [CrossRef]
  19. Alfian, G.; Syafrudin, M.; Fitriyani, N.L.; Anshari, M.; Stasa, P.; Svub, J.; Rhee, J. Deep Neural Network for Predicting Diabetic Retinopathy from Risk Factors. Mathematics 2020, 8, 1620. [Google Scholar] [CrossRef]
  20. Norhaslinda, Z.A.; Mustafa, M.; Tengku, H.; Brian, D.; Azizul, B. System Dynamics Modelling and its implications for Childhood Obesity Prevention: Evidence from Improving the Consumption of Portion Size and Meal Frequency. Appl. Math. Sci. 2014, 8, 3283–3296. [Google Scholar] [CrossRef] [Green Version]
  21. Forrester, J.W. Industrial Dynamics; MIT Press: Cambridge, MA, USA, 1961. [Google Scholar]
  22. Sterman, J.D. Business Dynamics: Systems Thinking and Modelling for a Complex World; Irwan McGraw-Hill: Boston, MA, USA, 2000. [Google Scholar]
  23. Derouich, M.; Boutayeb, A. The effect of physical exercise on the dynamics of glucose and insulin. J. Biomech. 2002, 35, 911–917. [Google Scholar] [CrossRef]
  24. Ahlborg, G.; Felig, P.; Hagenfeldt, L.; Hendler, R.; Wahren, J. Substrate turnover during prolonged exercise in man. J. Clin. Investig. 1974, 53, 1080–1090. [Google Scholar] [CrossRef] [PubMed]
  25. Wasserman, D.; Williams, P.; Lacy, D.; Goldstein, R.; Cherrington, A. Exercise-induced fall in insulin and hepatic carbohydrate metabolism during muscular work. Am. J. Physiol. 1989, 256, E500–E509. [Google Scholar] [CrossRef] [PubMed]
  26. Duun-Henriksen, A.; Schmidt, S.; Roge, R.; Moller, J.; Norgaard, K.; Jorgensen, J.; Madsen, H. Model identification using Stochastic Differential Equation grey-box models in diabetes. J. Diabetes Sci. Technol. 2013, 7, 431–440. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  27. Anirban, R.; Robert, S. Dynamic Modelling of Exercise effects on plasma glucose and insulin levels. J. Diabetes Sci. Technol. 2007, 1, 338–346. [Google Scholar] [CrossRef] [Green Version]
  28. Lu, G.; Hidenori, K. A Mathematical model of brain glucose homeostasis. Theor. Biol. Med. Model. 2009, 6, 1–24. [Google Scholar] [CrossRef] [Green Version]
  29. Abdel-Hamid, T.K. Exercise and diet in obesity treatment: An integrative system dynamics perspective. Med. Sci. Sport. Exerc. 2003, 35, 400–413. [Google Scholar] [CrossRef] [Green Version]
  30. Fisher, M. A semi-closed-loop algorithm for the control of blood glucose levels in diabetics. IEEE Trans. Biomed. Eng. 1991, 38, 57–61. [Google Scholar] [CrossRef]
Figure 1. Modified Bergman Minimal Model.
Figure 1. Modified Bergman Minimal Model.
Mathematics 11 00306 g001
Figure 2. Model of the semi-closed causal loop diagram of the interaction between meal, physical exercise and stress in glucose-insulin dynamics.
Figure 2. Model of the semi-closed causal loop diagram of the interaction between meal, physical exercise and stress in glucose-insulin dynamics.
Mathematics 11 00306 g002
Figure 3. ACF for first 24 hours’ plasma glucose residuals from the SDE model.
Figure 3. ACF for first 24 hours’ plasma glucose residuals from the SDE model.
Mathematics 11 00306 g003
Figure 4. ACF for four consecutive days’ plasma glucose residuals from the SDE model.
Figure 4. ACF for four consecutive days’ plasma glucose residuals from the SDE model.
Mathematics 11 00306 g004
Figure 5. CGM observed and one-step predicted values trend for plasma glucose for 24 h on the first day of the experiment.
Figure 5. CGM observed and one-step predicted values trend for plasma glucose for 24 h on the first day of the experiment.
Mathematics 11 00306 g005
Figure 6. CGM observed and one-step predicted values trend for plasma glucose for four consecutive days.
Figure 6. CGM observed and one-step predicted values trend for plasma glucose for four consecutive days.
Mathematics 11 00306 g006
Figure 7. Plasma glucose levels of a type II diabetic patient using base run data from the CGM.
Figure 7. Plasma glucose levels of a type II diabetic patient using base run data from the CGM.
Mathematics 11 00306 g007
Figure 8. Plasma glucose levels using identified parameters for the extended Bergman Minimal Model for a type II diabetic patient with a 33.33% increase in amount of ingested carbohydrates.
Figure 8. Plasma glucose levels using identified parameters for the extended Bergman Minimal Model for a type II diabetic patient with a 33.33% increase in amount of ingested carbohydrates.
Mathematics 11 00306 g008
Figure 9. Plasma glucose levels using identified parameters for the extended Bergman Minimal Model for a type II diabetic patient with a 33.33% reduction in stress levels.
Figure 9. Plasma glucose levels using identified parameters for the extended Bergman Minimal Model for a type II diabetic patient with a 33.33% reduction in stress levels.
Mathematics 11 00306 g009
Figure 10. Plasma glucose levels using identified parameters for the extended Bergman Minimal Model for a type II diabetic patient with reduction in the physical exercise intensity.
Figure 10. Plasma glucose levels using identified parameters for the extended Bergman Minimal Model for a type II diabetic patient with reduction in the physical exercise intensity.
Mathematics 11 00306 g010
Figure 11. Plasma glucose levels using identified parameters for the extended Bergman Minimal Model for a type II diabetic patient with reduction of 33.33% in both the amount of carbohydrates taken and stress levels while holding physical exercise intensity at the base run level.
Figure 11. Plasma glucose levels using identified parameters for the extended Bergman Minimal Model for a type II diabetic patient with reduction of 33.33% in both the amount of carbohydrates taken and stress levels while holding physical exercise intensity at the base run level.
Mathematics 11 00306 g011
Figure 12. Plasma glucose levels using identified parameters for the extended Bergman Minimal Model for a type II diabetic patient with an increase of 66.67% in the amount of carbohydrates is taken with an increase in stress levels by 233.33% while holding physical exercise intensity at the base run level.
Figure 12. Plasma glucose levels using identified parameters for the extended Bergman Minimal Model for a type II diabetic patient with an increase of 66.67% in the amount of carbohydrates is taken with an increase in stress levels by 233.33% while holding physical exercise intensity at the base run level.
Mathematics 11 00306 g012
Table 2. Selected simulated plasma glucose output at different intervention strategies.
Table 2. Selected simulated plasma glucose output at different intervention strategies.
Experimental and Modelling StrategiesInput Variables and MeasurementsAverage Plasma Glucose (mmol/L)Standard Deviation
Base Run: Scenario 1
No alteration of routine
Average carbs (grams) (120)9.221.875
Average   exercise   intensity   ( min 1 ) (Moderately high to low)
Stress level (0.15)
Scenario 2:
Increase both carbohydrates and stress levels while holding the exercise at base run levels
Average carbs (grams) (200)12.042.016
Average   exercise   intensity   ( min 1 ) (Moderately high to low)
Stress level (0.5)
Scenario 3:
Lower the stress level, maintain both carbs and exercise intensity at base run levels
Average carbs (grams) (120)8.951.871
Average   exercise   intensity   ( min 1 ) (Moderately high to low)
Stress level (0.1)
Scenario 4:
Lower the carbohydrates, maintain the stress level and the exercise intensity at the base run levels
Average carbs (grams) (80)8.141.831
Average   exercise   intensity   ( min 1 ) (Moderately high to low)
Stress level (0.15)
Scenario 5:
Lower the carbohydrates and the stress levels, maintain exercise at base run levels
Average carbs (grams) (80)7.861.827
Average   exercise   intensity   ( min 1 ) (Moderately high to low)
Stress level (0.1)
Scenario 6:
Lower the exercise intensity, maintain the carbohydrates and stress at the base run level
Average carbs (grams) (120)9.191.878
Average   exercise   intensity   ( min 1 ) (Moderately low to zero)
Stress level (0.15)
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Omwenga, V.O.; Madhumati, V.; Vinay, K.; Srikanta, S.; Bhat, N. Mathematical Modelling of Combined Intervention Strategies for the Management and Control of Plasma Glucose of a Diabetes Mellitus Patient: A System Dynamic Modelling Approach. Mathematics 2023, 11, 306. https://doi.org/10.3390/math11020306

AMA Style

Omwenga VO, Madhumati V, Vinay K, Srikanta S, Bhat N. Mathematical Modelling of Combined Intervention Strategies for the Management and Control of Plasma Glucose of a Diabetes Mellitus Patient: A System Dynamic Modelling Approach. Mathematics. 2023; 11(2):306. https://doi.org/10.3390/math11020306

Chicago/Turabian Style

Omwenga, Vincent O., Vaishnav Madhumati, Kumar Vinay, Sathyanarayan Srikanta, and Navakanta Bhat. 2023. "Mathematical Modelling of Combined Intervention Strategies for the Management and Control of Plasma Glucose of a Diabetes Mellitus Patient: A System Dynamic Modelling Approach" Mathematics 11, no. 2: 306. https://doi.org/10.3390/math11020306

APA Style

Omwenga, V. O., Madhumati, V., Vinay, K., Srikanta, S., & Bhat, N. (2023). Mathematical Modelling of Combined Intervention Strategies for the Management and Control of Plasma Glucose of a Diabetes Mellitus Patient: A System Dynamic Modelling Approach. Mathematics, 11(2), 306. https://doi.org/10.3390/math11020306

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop