Evaluating the Importance of Glucagon in the Insulin-Glucose Regulatory System: A Mechanistic Modeling Approach

Evaluating the Importance of Glucagon in the Insulin-Glucose Regulatory System: A Mechanistic Modeling Approach | bioRxiv /* */ /* */ <!-- <!-- /*! * yepnope1.5.4 * (c) WTFPL, GPLv2 */ (function(a,b,c){function d(a){return"[object Function]"==o.call(a)}function e(a){return"string"==typeof a}function f(){}function g(a){return!a||"loaded"==a||"complete"==a||"uninitialized"==a}function h(){var a=p.shift();q=1,a?a.t?m(function(){("c"==a.t?B.injectCss:B.injectJs)(a.s,0,a.a,a.x,a.e,1)},0):(a(),h()):q=0}function i(a,c,d,e,f,i,j){function k(b){if(!o&&g(l.readyState)&&(u.r=o=1,!q&&h(),l.onload=l.onreadystatechange=null,b)){"img"!=a&&m(function(){t.removeChild(l)},50);for(var d in y[c])y[c].hasOwnProperty(d)&&y[c][d].onload()}}var j=j||B.errorTimeout,l=b.createElement(a),o=0,r=0,u={t:d,s:c,e:f,a:i,x:j};1===y[c]&&(r=1,y[c]=[]),"object"==a?l.data=c:(l.src=c,l.type=a),l.width=l.height="0",l.onerror=l.onload=l.onreadystatechange=function(){k.call(this,r)},p.splice(e,0,u),"img"!=a&&(r||2===y[c]?(t.insertBefore(l,s?null:n),m(k,j)):y[c].push(l))}function j(a,b,c,d,f){return q=0,b=b||"j",e(a)?i("c"==b?v:u,a,b,this.i++,c,d,f):(p.splice(this.i++,0,a),1==p.length&&h()),this}function k(){var a=B;return a.loader={load:j,i:0},a}var l=b.documentElement,m=a.setTimeout,n=b.getElementsByTagName("script")[0],o={}.toString,p=[],q=0,r="MozAppearance"in l.style,s=r&&!!b.createRange().compareNode,t=s?l:n.parentNode,l=a.opera&&"[object Opera]"==o.call(a.opera),l=!!b.attachEvent&&!l,u=r?"object":l?"script":"img",v=l?"script":u,w=Array.isArray||function(a){return"[object Array]"==o.call(a)},x=[],y={},z={timeout:function(a,b){return b.length&&(a.timeout=b[0]),a}},A,B;B=function(a){function b(a){var a=a.split("!"),b=x.length,c=a.pop(),d=a.length,c={url:c,origUrl:c,prefixes:a},e,f,g;for(f=0;f<d;f++)g=a[f].split("="),(e=z[g.shift()])&&(c=e(c,g));for(f=0;f<b;f++)c=x[f](c);return c}function g(a,e,f,g,h){var i=b(a),j=i.autoCallback;i.url.split(".").pop().split("?").shift(),i.bypass||(e&&(e=d(e)?e:e[a]||e[g]||e[a.split("/").pop().split("?")[0]]),i.instead?i.instead(a,e,f,g,h):(y[i.url]?i.noexec=!0:y[i.url]=1,f.load(i.url,i.forceCSS||!i.forceJS&&"css"==i.url.split(".").pop().split("?").shift()?"c":c,i.noexec,i.attrs,i.timeout),(d(e)||d(j))&&f.load(function(){k(),e&&e(i.origUrl,h,g),j&&j(i.origUrl,h,g),y[i.url]=2})))}function h(a,b){function c(a,c){if(a){if(e(a))c||(j=function(){var a=[].slice.call(arguments);k.apply(this,a),l()}),g(a,j,b,0,h);else if(Object(a)===a)for(n in m=function(){var b=0,c;for(c in a)a.hasOwnProperty(c)&&b++;return b}(),a)a.hasOwnProperty(n)&&(!c&&!--m&&(d(j)?j=function(){var a=[].slice.call(arguments);k.apply(this,a),l()}:j[n]=function(a){return function(){var b=[].slice.call(arguments);a&&a.apply(this,b),l()}}(k[n])),g(a[n],j,b,n,h))}else!c&&l()}var h=!!a.test,i=a.load||a.both,j=a.callback||f,k=j,l=a.complete||f,m,n;c(h?a.yep:a.nope,!!i),i&&c(i)}var i,j,l=this.yepnope.loader;if(e(a))g(a,0,l,0);else if(w(a))for(i=0;i (function(w,d,s,l,i){w[l]=w[l]||[];w[l].push({'gtm.start':new Date().getTime(),event:'gtm.js'});var f=d.getElementsByTagName(s)[0];var j=d.createElement(s);var dl=l!='dataLayer'?'&l='+l:'';j.src='//www.googletagmanager.com/gtm.js?id='+i+dl;j.type='text/javascript';j.async=true;f.parentNode.insertBefore(j,f);})(window,document,'script','dataLayer','GTM-M677548'); Skip to main content Home About Submit ALERTS / RSS Search for this keyword Advanced Search New Results Evaluating the Importance of Glucagon in the Insulin-Glucose Regulatory System: A Mechanistic Modeling Approach Mackenzie Dalton , Emmanuel Asante-Asamani , View ORCID Profile James Greene doi: https://doi.org/10.1101/2024.07.31.606029 Mackenzie Dalton a Mathematics Department, Clarkson University , Potsdam, New York, 13699 Find this author on Google Scholar Find this author on PubMed Search for this author on this site Emmanuel Asante-Asamani a Mathematics Department, Clarkson University , Potsdam, New York, 13699 Find this author on Google Scholar Find this author on PubMed Search for this author on this site James Greene a Mathematics Department, Clarkson University , Potsdam, New York, 13699 Find this author on Google Scholar Find this author on PubMed Search for this author on this site ORCID record for James Greene For correspondence: jgreene{at}clarkson.edu Abstract Full Text Info/History Metrics Preview PDF Abstract The dynamics of insulin and glucose are tightly regulated. The pancreatic islets of Langerhans contain both beta and alpha cells which produce insulin and glucagon, respectively. Insulin is the only hormone in the body that lowers blood glucose levels by acting like a key for glucose to enter cells. Without insulin, cells cannot utilize glucose, their primary source of energy. In contrast, glucagon functions as a hormone which elevates blood glucose levels by promoting the breakdown of glycogen in the liver. Maintaining blood glucose within a safe range is vital since both excessively high and low levels can be life-threatening (hyperglycemia and hypoglycemia, respectively), and these two hormones work together to achieve this balance. In this work we aim to underscore the significance of glucagon in the insulin-glucose regulatory system. We construct a three-compartment mechanistic model that includes insulin, glucose, and glucagon, which is then validated by fitting to publicly available from an intravenous glucose tolerance test (IVGTT). After model validation, we investigate how removing glucose feedback from insulin secretion, as seen in insulin-dependent diabetes, disrupts the regulation of glucose and glucagon. To do this, we simulate the model (a) when insulin secretion is reduced to mimic an insufficient dose of insulin, (b) when the peak of insulin action is delayed mimicking a dosing delay of insulin, and (c) when both occur simultaneously. Lastly, we test different half-lives of insulin to evaluate how an increased half-life of manufactured insulin may further disrupt the system. We find that when insulin secretion is decreased, glucagon still responds to high glucose levels by decreasing glucagon production. This suggests that in cases of type 2 diabetes, where glucagon secretion is elevated despite high levels of glucose, a lack of insulin response may not be the sole cause for glucagon dysfunction. We also find that delaying insulin secretion increases the risk of a hypoglycemic event through a suppression of glucagon production. Initially, the spike in glucose causes glucagon secretion to be reduced; this is then followed by the delay in insulin peak which then continues to suppress glucagon despite blood glucose levels falling, leading to a lack of response by glucagon and a subsequent hypoglycemic event. Furthermore, we find that a higher half-life of insulin causes it to remain longer in the blood stream, inhibiting glucagon’s response to severely low glucose levels (glucose levels less than 3.9 mmol/L). This sheds light on why patients taking exogenous insulin, which has a longer half-life than endogenous insulin, may have difficulty recovering from hypoglycemic events. Hence, our model suggests that keeping the half-life of exogenous insulin below 10 minutes and administering it immediately after meals could help reduce the risk of hypoglycemic events in patients with type 1 or insulin dependent diabetes. Overall, we highlight how a disruption in the feedback between insulin and glucose not only alters blood glucose levels, but also glucagon response, which may lead to further disruption of the system. 1. Introduction Blood glucose levels are a crucial component of homeostatic regulation. Therefore, maintaining healthy blood glucose levels is critical to health. The regulation of glucose levels in the body is accomplished through a complex interplay of hormones, primarily insulin and glucagon, secreted by the pancreas. These hormones work together to ensure that blood glucose levels remain within a narrow and well-defined range. The significance of maintaining healthy glucose levels cannot be overstated, as levels outside a healthy range can have severe consequences. Elevated blood glucose levels, known as hyper-glycemia, can lead to a range of complications such as cardiovascular diseases, kidney dysfunction, and nerve damage [ 22 , 32 ]. In extreme cases, uncontrolled hyperglycemia can lead to diabetic ketoacidosis (DKA) which is life-threatening and is the leading cause of death in children with type 1 diabetes [ 47 ]. Similarly, low blood glucose levels, or hypoglycemia, also pose significant risks. Glucose is the primary energy source for the brain and other vital organs. When blood glucose levels drop too low, it can result in symptoms such as dizziness, confusion, and, in severe cases, loss of consciousness. In situations where the body cannot access an adequate supply of glucose, the consequences can be fatal [ 2 ]. The pancreatic islets of Langerhans are cell clusters found in the pancreas that respond to glucose levels by secreting either a blood glucose lowering hormone, insulin, or a blood glucose raising hormone, glucagon [ 34 ]. Alpha ( α ) cells, beta ( β ) cells, and delta ( δ ) cells are the primary cells composing the pancreatic islets of Langerhans ( Figure 1 ), each of which respond to both increases and decreases in blood sugar levels differently [ 34 ]. For simplicity, we do not include δ cell dynamics in our model since glucagon and insulin are sufficient to capture data utilized in this work (see Section 4.1 ). Download figure Open in new tab Figure 1: Pancreatic islet of Langerhans and the hormones secreted by α cells, β cells, and δ cells. [ 45 ]. To raise blood glucose levels, the pancreatic α cells secrete the hormone glucagon, also known as anti-insulin, to raise blood glucose levels. Glucagon works by triggering the liver to convert stored glucose (glycogen) into glucose, which is then released into the bloodstream; this process is known as glycogenolysis [ 23 ]. Moreover, glucagon can also raise blood glucose levels by preventing the liver from taking in and storing glucose so that more glucose remains in the blood stream [ 23 ]. To lower blood glucose levels, the pancreatic β cells secrete the hormone insulin to lower blood glucose levels. Insulin works by acting as a “key” that allows glucose to enter into a cell and subsequently lowers blood glucose levels. Without insulin, glucose cannot enter the cell [ 10 ]. Moreover, insulin blocks the production of glucagon from α cells, which leads to less production of glucose by the body [ 3 ]. In 1979 Richard Bergman and colleagues introduced a model of the insulin-glucose regulatory system [ 6 ] to help understand data from an intravenous glucose tolerance test (IVGTT). The test involves an injection of glucose into the blood stream followed by monitoring of blood glucose and insulin levels for a period of time. The test can help in the diagnosis of type 2 diabetes by calculating an individual’s insulin sensitivity - a measurement of the ability of insulin to stimulate the uptake of glucose. This model was coined the minimal model , as it was the simplest quantitative biological framework that was able to successfully capture insulin-glucose longitudinal data [ 12 ]. However, while this model is still utilized as of writing, it does possess limitations. Two of the main criticisms of the minimal model are that (a) insulin production is dependent on time, making the system non-autonomous and (b) a non-observable remote insulin compartment is introduced that aims to model the insulin-dependent uptake of glucose [ 33 ]. Beyond the minimal model, there have been a variety of other models that aim to describe insulin-glucose dynamics [ 41 , 26 , 33 , 11 ]. We note that a number of these models aim to remove the time dependence and remote insulin compartment from the model; however, this often leads to highly complex and overparameterized (i.e. nonidentifiable) systems and the appeal of the minimal model being minimal while still capturing the dynamics of the system should not be be overlooked. Moreover, in alignment with the growing biological research on glucagon, additional models have been introduced to explore the relationship between insulin, glucose, and glucagon [ 41 , 16 , 39 , 48 ]. However, many of these models are derived from the minimal model, which results in the aforementioned criticisms. Efforts to understand the role of α cells in secreting glucagon have recently started to gain traction, as the role of glucagon in maintaining glucose homeostasis is more well understood [ 27 , 17 ]. Some sources suggest that α cells are primarily regulated by neighboring β cells, while others contend that α cells possess intrinsic mechanisms for detecting and responding to glucose [ 25 , 8 ]. Moreover, research has suggested that glucagon may play a role in the development of both type 1 and type 2 diabetes [ 28 , 21 ]. More generally, there is also research investigating the role of glucagon in regulating blood glucose levels, as well as possible mechanisms that may cause glucagon dysfunction in both type 1 and type 2 diabetes [ 25 ]. In this work, we develop a model of the insulin-glucose regulatory system that is self-regulatory and includes a direct interaction between insulin and glucose (i.e.there is no remote insulin compartment). Furthermore, we add glucagon as a third compartment of our model. In fact, we note that the addition of glucagon is necessary to fit to time course data of insulin and glucose, highlighting an important role for glucagon in insulin-glucose dynamics. The model is presented in Section 2 , and results of model fitting can be seen in Section 4.1 ; Section 3 contains methods related to the (publicly available) experimental data incorporated into the mathematical model, as well as parameter estimation techniques. After model validation, we investigate in Section 4.3 how removing glucose feedback from insulin secretion, as is seen in insulin-dependent diabetes where exogenous insulin is administered, disrupts the regulation of glucose and glucagon. To do this, we simulate the model (a) when insulin secretion is reduced to mimic an insufficient dose of insulin, (b) when the peak of insulin action is delayed mimicking a dosing delay of insulin, and (c) when both occur simultaneously. Lastly, we varied the half-life of insulin to evaluate how an increased half-life of manufactured insulin (when compared with endogenous insulin) may further disrupt the system. We found that when insulin secretion is decreased, glucagon still responds to high glucose levels by decreasing glucagon production ( Figure 8a ). Next, we found that delaying insulin secretion increases the risk of a hypoglycemic event through a suppression of glucagon production ( Figure 8b ). Lastly, we found that the higher the half-life of insulin, the longer it remains in the blood stream, inhibiting glucagon’s response to severely low glucose levels (glucose levels less than 3.9 mmol/L; Figure 9 ). Overall, we highlight how a disruption in the feedback between insulin and glucose does not only alter blood glucose levels but also glucagon response, leading to further disruption of the system. 2. Model We develop a three compartment model of insulin and glucose dynamics that includes insulin, glucose, and glucagon. Mass action kinetics are assumed for most interactions, unless otherwise specified. The major biological features of the model are the following: Glucose stimulates the production of insulin [ 14 ]. The production of insulin was assumed to saturate as a function of glucose [ 12 ]. That is, the maximum insulin production rate ( k 1 ) is independent of glucose. Insulin regulates the insulin dependent uptake of glucose ( δ 2 I ) which is the primary mechanism by which glucose enters cells [ 13 ]. A small amount of glucose can enter cells through an insulin independent uptake of glucose ( δ 3 ) [ 13 ]. The production of glucagon is inhibited by both high levels of insulin ( p I ) and glucose ( z ( G − h ) + )[ 38 ]. The inhibition by glucose is modeled as z ( G − h ) + due to the assumption that when G is below a baseline value h , glucose does not inhibit glucagon at all, as observed in the response of glucagon to low glucose levels. Thus, when G < h glucose no longer inhibits glucagon production. For x ∈ ℝ, recall that x + denotes the positive part of x : The production of glucose in the body (that is, not due to a glucose injection) is regulated by glucagon ( k 3 G 𝓁 ) [ 18 ]. This occurs via glucagon signaling to the liver to break down stored glucose and releasing it into the blood stream [ 23 ]. Insulin and glucagon both decay naturally ( δ 1 and δ 4 , respectively) [ 37 ]. Under these assumptions, the model schematic of the three compartment insulin-glucose-glucagon model is shown in Figure 2 . A system of three ordinary differential equations incorporating the above assumptions is given below: Download figure Open in new tab Figure 2: Three compartment model schematic representing insulin-glucose-glucagon model (1) - (3) . It can be shown that the above system possess a unique positive steady state (see Section 7.1 in the Supplementary Material). We note that we have also developed a two compartment model of insulin and glucose alone to assess the necessity of glucagon in explaining insulin-glucose dynamics; see Section 7.2 in the Supplemental Material. Results suggest that three compartments are necessary to fully capture the utilized experimental data (see Figure 11 ). To simulate a glucose injection from the IVGTT, we incorporate an external input u ( t ) to equation (2) of the form Here, is the applied (constant) dosage which was taken to match the peak glucose value of the experimental data (see Table 1 for value), t i is start time of the dosage, and T is the length of infusion. In this case, system (1) - (3) takes the form View this table: View inline View popup Download powerpoint Table 1: Table of initial estimate ranges, best fit parameter values, and units for the insulin-glucose-glucagon model 1 - 3 . View this table: View inline View popup Download powerpoint Table 2: Table of initial estimate ranges, best fit parameter values, and units for the insulin-glucose model (32) - (33) . 3. Methods 3.1. Experimental Data We utilized data collected from a study conducted by Manell et al. [ 29 ], where ten pigs were subjected to an Intravenous Glucose Tolerance Tests (IVGTT). These pigs underwent measurements of glucose, insulin, and glucagon, which were taken both before glucose infusion and at 2, 5, 10, 20, 30, 45, 60, 90, 120, and 180 minutes post infusion. Of the ten pigs, the final glucagon measurement was missing for three of them. We fit the model (4) - (6) to the median (with respect to the ten pigs) temporal data. We used the median to avoid outliers skewing the estimated parameter values, particularly with respect to glucagon, where the measurements exhibited a large variance (see Figure 3 ). Furthermore, we computed a 95% confidence interval with respect to the median via bootstrapping [ 1 ]. Specifically, we generated 10,000 new data sets, with each set consisting of 10 pigs. We then calculated the median of each new data set to generate 10,000 medians which gave a distribution of the medians at each time point. From this data we constructed a sampling distribution for the median, which was used to calculate the 95% confidence interval. The longitudinal pig data, together with the median and estimated confidence intervals are provided in Figure 3 . We highlight that while glucose and insulin exhibited minimal variation over time, the levels of glucagon displayed substantial variability, particularly during the 30 minutes following the IVGTT. Download figure Open in new tab Figure 3: Experimental Data. Plot of the pig data (insulin, glucose, and glucagon over a 180 minute period), medians, and estimated 95% confidence intervals for an IVGTT. Left is insulin, middle is glucose and right is glucagon. 3.2. Parameter Estimation To estimate rate parameters in 1-3, we fit to the median of the pig data presented in Section 3.1 . The median of the pig data is measured at times t i , i = 1, 2,…, n ; thus the experimental data consists of tuples of the form . For example, G d, i denotes the median glucose level at time t i . We utilize a normalized least-squares inspired cost functional which penalizes deviations from all three compartments (insulin, glucose, and glucagon). Specifically, we define where Here I ( t i ; θ ) denotes the solution of the insulin compartment of the ODE system 1-3 evaluated at time t i , as a function of the model parameters θ ; similar statements holds for G ( t i , θ ) and G 𝓁 ( t i , θ ). The parameter vector is then estimated to minimize C , with respect to non-negativity constraints, utilizing fmincon in M at-lab ’s nonlinear optimization toolbox. We note that when calibrating the minimal model to the same data set, glucagon is not explicitly considered (see equations 12 - 14 ), and thus the following analogous cost functional is utilized: Multi-start optimization was performed, and the initial parameter ranges utilized are provided in Table 1 . We note that while the initial conditions for the optimization algorithm are contained in this range, the best fit parameter value may lie outside of it. To ensure that a global minimum was achieved, we repeated the optimization 1000 times while randomizing initial estimates between the ranges seen in Table 1 . We used an n -dimensional Sobol set transformed to the hypercube [0, γ ] ρ where ρ is the length of the parameter vector θ and γ is the upper bound for each parameter value. During optimization, all parameter values were constrained to remain non-negative. After the multi-start optimization was performed, we constructed tornado plots to ensure that we located the global minimum of the cost function. The tornado plots were generated by varying one parameter at a time within 10% of its computed optimal value, and subsequently calculating the cost functional value. If a change in parameter value led to a decrease in the cost functional, the parameter value was changed accordingly, and the process was repeated until all local changes in parameter values led to only increases in the cost 5. 4. Results and Discussion 4.1. Model fitting We begin by estimating parameters in the insulin-glucose-glucagon model ( (1) - (3) ) by fitting the model to the data presented in Section 3.1 , as discussed in Section 3.2 ; the code to generate all figures can be found on github. The optimal parameter set is defined to be the one that minimizes the cost functional defined in equation (7) . The results are provided in Figure 4 . We see that our model is able to capture the median of the data well. The cost functional at the optimal parameters is C = 0.96976, with an adjusted R 2 value of 0.9283. Download figure Open in new tab Figure 4: Model Fits. The best fit for the insulin-glucose-glucagon model 1-3. Here, the constant dosage was used with initial time t i = 0 and dose length T = 1. The initial conditions for insulin, glucose and glucagon were taken to be the initial values of the data. The cost functions evaluated at the optimal set of parameters are C = 0.96976 (insulin-glucose-glucagon cost) and (insulin-glucose cost). See Table 1 for the corresponding estimated parameters. Left is insulin, middle is glucose and right is glucagon. The estimated parameter values are provided in Table 1 . Figure 5 also shows tornado plots when individually varying each best fit parameter within 5%, 10%, 15%, and 20% of its original value. Note that varying each parameter does not lead to a decrease in the cost function C , providing evidence that we are indeed at a global minimum. Download figure Open in new tab Figure 5: Local Minimum Validation. Tornado plots for the change in the cost function ( C ) when varying each parameter by the given percentage of its optimal value. The vertical line represents the minimum residual. Parameter values in (a) are changed by 5%, (b) are changed 10%, (c) are changed by 15% and (d) are changed by 20%. For comparison, we also fit the minimal model developed by Bergman [ 4 ] to the same data set. The minimal model consists of three compartments: insulin ( I ), glucose ( G ) and remote insulin ( X ). Here, insulin does not directly inhibit glucose; rather insulin enters a remote compartment ( X ) and X acts to inhibit glucose. Insulin production is a function of both time and the degree to which glucose exceeds a baseline threshold value h . Glucose production is constant with respect to a baseline level of glucose, which is dependent on the specific data set. More precisely, G b is taken to be the initial value of glucose prior to the start of the IVGTT. Note that G b is not estimated from the data. The production of remote insulin depends on I , and glucose uptake occurs through both an insulin dependent and independent pathway. Lastly, both insulin and remote insulin have natural decay rates n and P 2 , respectively. The differential equations describing these dynamics are provided below; for more details on the model formulation, see [ 4 ]. Figure 6 shows the result of the best fit of the minimal model (12) - (14) to the same data set utilizing the procedure outlined in Section 3.2 compared with the best fit for the insulin-glucose-glucagon model. Recall that when estimating parameters in the minimal model, we consider the cost function (11) . The cost for the minimal model is 1.5 whereas cost for the insulin-glucose-glucagon model is 0.78. To adjust for the insulin-glucose-glucagon model possessing a larger number of parameters, we also calculated the adjusted R 2 value for each model; the adjusted R 2 for the minimal model is 0.77, while that of the insulin-glucose-glucagon model is 0.92. The best fit parameter values and tornado plots for the minimal model can be found in the Supplemental Material (see Section 7.3 , Table 3 , Figure 14 ). View this table: View inline View popup Download powerpoint Table 3: Table of initial estimate ranges, best fit parameter values, and units for the minimal model (12) - (14) . Download figure Open in new tab Figure 6: Minimal Model Fits. The best fit for the minimal model (12) - (14) together with the best fit of the three compartment model (1) - (3) . The cost function at the given parameters for the minimal model fit is , with a corresponding R 2 value of 0.7569. Left is insulin and right is glucose. 4.2. Insulin Sensitivity and Glucose Effectiveness One of the main applications of the minimal model is its ability to provide a patient’s glucose effectiveness and insulin sensitivity to both aid in diagnosing a patient with type 2 diabetes, as well as understanding the severity of the disease [ 43 , 15 ].We now explore how to calculate both insulin sensitivity and glucose effectiveness using our three compartment insulin, glucose and glucagon model (1) - (3) . We utilize Bergman’s definitions to calculate both glucose effectiveness and insulin sensitivity. Specifically, Bergman defines glucose effectiveness ( E ) as the quantitative enhancement of glucose disappearance due to an increase in the plasma glucose concentrations [ 5 ]. In other words, glucose effectiveness is the ability of glucose to stimulate its own uptake and is defined quantitatively as follows: Here Ġ is the rate of change of the plasma glucose concentration G , i.e. effectiveness for model (1) - (3) is given by Ġ = dG / dt . Thus, glucose effectiveness for model (1) - (3) is given by Bergman also defines the insulin sensitivity ( S ) as the quantitative influence of insulin to increase the enhancement of glucose disappearance [ 5 ]. That is, insulin sensitivity is the ability of insulin to stimulate glucose disappearance and is given by the following formula: Here the subscript SS indicates that the quantities are measured at steady state. Using the above definition, insulin sensitivity for model (1) - (3) is The insulin sensitivity factor is easy to interpret as δ 2 is the parameter that controls the impact of insulin on glucose clearance. A high δ 2 value would result in an increased clearance of glucose and indicate sensitivity to insulin, whereas a low δ 2 would result in less clearance of glucose and indicate insulin resistance [ 24 ]. Our insulin sensitivity factor estimated by fitting the median of the pig data is 0.17. In [ 5 ], Bergman estimates insulin sensitivity factors between (using the conversion factor 1 µU = 0.0347 µg ; [ 20 ]). Hence, our insulin sensitivity factor appears to be biologically reasonable. 4.3. Removing Glucose Feedback on Insulin Diabetes is a disease where the body either becomes resistant to insulin (type 2 diabetes) or there is a low/absent production of insulin due to an autoimmune attack on the insulin producing β cells (type 1 diabetes). In the latter case, and occasionally in the former as well, the patient requires insulin injections in order to regulate blood glucose levels which results in the loss of direct feedback of glucose on insulin production. This loss of communication has a large impact on the insulin-glucose regulatory system; indeed, the impact is significant to the degree so that glucagon no longer functions properly despite α cells remaining functional [ 7 ]. Understanding the disruption in glucagon dynamics caused by insulin disruption may be critical to better understanding the role of glucagon both before and after the diagnosis of diabetes. In this section, we investigate this loss of glucose feedback on both glucose and glucagon dynamics using a simplified system derived from (1) - (3) . A schematic of the considered system with no glucose feedback on insulin is provided in Figure 7 . We are thus analyzing a system where the communication between insulin and glucose is disrupted, so that insulin can be considered an external input, which operates independently of glucose ( G ) and glucagon ( G 𝓁 ). For clarity, we denote this external insulin concentration with the symbol Î = Î ( t ). This leads to the following two-dimensional system with external input Î : Download figure Open in new tab Figure 7: Schematic of the insulin-glucose-glucagon system with glucose feedback on insulin removed. The red X denotes that the interaction is not considered in system (19) - (20) , and insulin is considered as an external input to the system ( Î ). Download figure Open in new tab Figure 8: Dynamics in the Absence of Glucose Feedback. Dynamics of insulin, glucose, and glucagon (19) - (20) when altering the insulin profile from the model fit ( Figure 4 ). Recall that in this system, there is no feedback of glucose on insulin. The insulin profile in (a) is reduced, in (b) is delayed and in (c) is both reduced and delayed. Download figure Open in new tab Figure 9: Varying Insulin Half-lives. Insulin, glucose, and glucagon dynamics in system (19) - (20) when varying the half-life of insulin to 10, 60, and 120 minutes. Insulin dynamics are given by (21) . Left is insulin, middle is glucose and right is glucagon. To investigate how the removal of glucose feedback on insulin production impacts glucose and glucagon dynamics, we utilize the fit insulin profile from Figure 4 as a baseline, and vary it in different ways to simulate simplified scenarios that occur in diabetes (specifically in type 1 diabetes). We note that in the absence of any perturbation to the insulin profile Î , the solution to system (19) - (20) produces the dynamics observed in Figure 4 precisely (i.e. to the solution of the original model system (19) - (20) ). To simulate different scenarios occurring in diabetes, we alter the insulin input in two different ways and then combine these two alterations and investigate the impact of these alterations on glucose and glucagon dynamics. First, we simulate the dynamics of glucose and glucagon when decreasing insulin secretion; specifically, we take the insulin profile from Figure 4 and reduce the total (integral) insulin concentration by 50% ( Figure 8a , top left), which provides a new insulin profile Î . A reduction in insulin secretion can be seen in both type 1 and type 2 diabetes. In type 1 diabetes, a reduction in insulin secretion due to loss in beta cells mass may arise during disease progression. In type 2 diabetes, a reduction in insulin secretion may occur due to a decline in the function of beta cells [ 9 ]. As opposed to the case of a simple reduction in insulin, we also investigate the dynamics of glucose and glucagon when the peak insulin concentration is delayed. Here, we again take the curve associated with the insulin profile in Figure 4 , and delay the insulin dynamics by 20 minutes to produce the external signal Î (i.e. we shift insulin dynamics by 20 minutes, Figure 8b top right). This again has applications in both type 1 and type 2 diabetes. In type 1 diabetes, patients may inject insulin after a glucose spike is already observed resulting in a delayed reaction of insulin to a glucose spike, while in type 2 diabetes, it has been shown that insulin secretion may be delayed [ 42 ]. Lastly, we combine both the reduction in insulin secretion and the time delay to create a perturbed insulin profile Î ( Figure 8c bottom), both of which have application in type 1 and type 2 diabetes for the aforementioned reasons. The results of these three scenarios are shown in Figure 8 . Figure 8a shows the response of system (19) - (20) when the total external insulin concentration is reduced by 50%. As expected, we see that a reduction in insulin secretion leads to higher glucose levels. That is, when compared with the original insulin profile (around 20 mins), glucose levels remain higher over time; however, we note that the glucose spike is relatively constant with respect to both models (approximately 20 mmol/L), meaning that less insulin does not increase the glucose peak. On the other hand, we observe that glucagon has an increased minimum value compared to the original insulin profile, which is most likely due to the reduction in insulin peak. This is because the reduction in insulin leads to less suppression of glucagon production. Interestingly though, the maximum of glucagon is greatly reduced despite lower levels of insulin. Therefore, this reduction in glucagon production is caused by the higher blood glucose levels. This is particularly interesting because while it has been documented that glucagon secretion depends on both insulin and glucose, the role each plays is not fully understood. The results appearing in Figure 8a thus suggest that high glucose levels can greatly reduce glucagon production in the near absence of insulin. This demonstrates that the irregular patterns of glucagon secretion in type 2 diabetes (in particular, over secretion) may not be explained solely by the dysfunction of β cells, supporting research that has suggested the α cells (the glucagon-secreting cells) are also dysfunctional in type 2 diabetes [ 19 , 31 , 46 ]. Figure 8b shows the simulation results when delaying insulin secretion, which is relevant in type 1 diabetes in the case when insulin is administered after, as opposed to during, a meal. Here, the insulin profile Î is obtained by a time delay of 20 minutes compared to figure Figure 4 . In this case, we now observe a significant hypoglycemic event taking place between 40 and 80 minutes post-glucose injection. Interestingly, in this simulation the maximum insulin concentration was not reduced, yet glucose levels significantly decreased when compared to the original insulin profile. This is most likely due to the dynamics of glucagon in the system. We see that glucagon initially decreases due to the rise in glucose levels. This reduction in glucagon secretion then reduces the production of glucose. Furthermore, glucagon remains suppressed due to increased insulin levels. That is, the insulin delay prevents glucagon from responding to lower glucose levels leading to a significant hypoglycemic event ( ≈ 3 mmol/L). Then, when insulin begins to fall and glucagon begins to respond to lower glucose levels, we observe a larger peak in glucagon which causes glucose levels to significantly overshoot its steady state value. The responses of both glucose and glucagon to the time shift in insulin suggest that a small delay in insulin action can lead to an increased risk of hypoglycemia, an observation that is important for patients administering insulin [ 44 ]. Furthermore, these results suggest that glucagon plays an important role regulating glucose levels by reducing the production of glucose. Thus, glucagon has dual roles in the regulation of blood sugar: in both helping lower glucose levels, as well as its primary role in raising blood glucose. Figure 8c shows simulation results when combining both a reduction and a delay of the insulin profile. Note here that glucose levels return close to baseline values without a hypoglycemic event although they remain elevated for a longer period of time. Moreover, the decrease in glucose levels is delayed when compared with the original simulation. Glucagon dynamics closely resemble the original glucagon dynamics, with the exception of the time delay and a reduction in the peak. These results suggest that when insulin is administered later, less insulin may be needed to avoid hypoglycemic events, which is most likely due to the proper regulation of glucagon. In the feedback system (1) - (3) , insulin and glucagon interact synergistically to regulate glucose levels. However, when feedback on insulin is removed, a time delay in insulin results in disruption to glucagon dynamics, and thus leads to impaired regulation of glucose. In order to optimize glucagon’s regulatory role in this case, our simulation suggest it would be best to time insulin administration so that the insulin peak and glucagon minimum occur at similar times, which seems to occur directly following a glucose spike. In fact, delays in insulin administration leading to low glucose levels is a fairly well documented phenomenon, particularly in type 1 diabetes [ 40 , 35 ] As a final application, we consider the response of the system to doses of insulin with different half-lives. The aim of this is to more closely resemble patients with diabetes, as insulins that are administered have a longer half-life compared to endogenous insulin [ 36 ]. We consider a simple exponentially decaying external insulin profile Î , with decay rate k and administered (constant) dosage U . The dose was assumed to be given at t = 0.01 minutes after the glucose injection, so that Recall that the decay rate is related to the half-life τ via The dose U was taken to be the maximum of the median insulin data minus the initial value of the median insulin data to closely resemble the endogenous insulin response upon glucose injection. That is, We considered three possible half-lives: 10 minutes to closely replicate endogenous insulin, 60 minutes to replicate a majority of meal time FDA approved insulins, and 120 minutes to investigate the effects of longer-lasting insulin. Note here that insulin was taken to be an intravenous dose rather than a sub-cutaneous dose for simplicity (i.e. pharmocokinetics are ignored). The results are shown in Figure 9 . We observe that as the half-life of insulin increases, the response of glucagon to low blood glucose levels becomes significantly more impaired, which leads to longer and more severe hypoglycemia. The failure of glucagon to respond to low blood glucose levels in type 1 diabetes is well documented [ 30 ]. Our simulations suggest that a main cause of the lack of response of glucagon is due to the clearance of insulin. In the original simulations, we see that insulin peaks and then decreases rapidly allowing for glucagon to respond to lower blood glucose levels. However, in the case of exogenous insulin with longer half-lives, the reduction in insulin clearance means that glucagon is suppressed for longer and therefore cannot respond appropriately to low blood glucose levels as demonstrated in Figure 9 . This also suggests that blood glucose levels may be more challenging to regulate in patients administering insulin, since if glucagon is not responding appropriately, a patient may no longer rely on the inhibition of glucagon to regulate glucose levels. This is especially clear in Figure 4 , where we see that when glucose undershoots its steady state value, glucagon is able to quickly reduce levels back to steady state. This demonstrates the degree to which the body relies on the anti-insulin properties of glucagon to regulate glucose levels, and highlights the importance of developing insulins that have half-lives closely resembling endogenous insulin in order to avoid impairing the function of glucagon. Conclusions and Future Work In this work we developed a three-compartment model of insulin, glucose and glucagon dynamics and validated our model by fitting it to the median of publicly available IVGTT pig data. We demonstrated that a two-compartment model of insulin and glucose was not sufficient to model the IVGTT pig data (see Supplementary Material Section 7.2 ), and thus highlighted the importance of glucagon in the insulinglucose regulatory system. We then removed the feedback of glucose on insulin production (see Figure 7 ), which allowed us to investigate the impact the loss of communication between insulin and the system has on glucose and, perhaps more interestingly, glucagon dynamics. This is a simplified phenomena of what occurs in both type 1 and type 2 diabetes. Overall, our model highlights the importance of glucagon in insulin-glucose dynamics. Without glucagon, we were unable to capture insulin-glucose dynamics well (Supplemental Material Figure 11 ). Furthermore, when glucose feedback on insulin production was removed, changes in insulin from the endogenous response leads to impairment of glucagon response ( Figures 8 and 9 ); inhibition of glucagon subsequently leads to prolonged and severe hypoglycemia ( Figures 8b and 9 ). Moreover, Figure 4 demon-strates the importance of the synergy between insulin and glucagon to optimally regulate glucose. Specifically, we observe that an undershoot of the glucose steady state value is rapidly counteracted via glucagon; note that this is not observed in Figure 9 , where feedback was removed. Moreover, upon the removal of feedback, it is still possible for glucagon to regulate the system, but only if the insulin profile is adjusted (see Figure 8(c) ); this highlights the important role of communication between glucagon and insulin to ensure that both are responding properly to glucose levels. However, for patients with diabetes, achieving insulin dynamics via external mechanisms similar to endogenous insulin secretion can be nearly impossible, and thus dysfunction of glucagon is often observed. Future work will continue to focus on the deregulation of insulin. We intend to use the proposed model (1) - (3) to understand how the progression of type 1 diabetes impacts this system by incorporating both the immune system and β cell loss. We will then combine this model to the insulin-glucose regulatory system via the insulin compartment (i.e. decreased β cell mass implies reduced insulin secretion) to understand how insulin, glucose and glucagon become increasingly deregulated throughout the progression of the disease. 6. Funding sources This research did not receive any specific grant from funding agencies in the public, commercial, or not-for-profit sectors. 7. Supplementary Material 7.1. Uniqueness of steady state of insulin-glucose-glucagon model In this section, we show that system (1) - (3) has a unique positive steady state. Let ( I ∗ , G ∗ , G 𝓁,∗ ) denote a steady state of system (1) - (3) . Thus, ( I ∗ , G ∗ , G 𝓁,∗ ) satisfy the following system of equations: Equation (24) implies that where Solving (25) for G 𝓁,∗ , and using I ∗ = f 1 ( G ∗ ) yields Define Thus, equation (26) can be written entirely in terms of G ∗ : Since f 3 ( G ∗ ) is increasing and f 3 (0) = 0, the first term on the left-hand side of (31) is decreasing (as a function of G ∗ ) and takes the value k 4 / δ 4 at G ∗ = 0. The second term, f 2 ( G ∗ ) is increasing and satisfies f 2 (0) = 0. Thus, equation (31) has precisely one positive solution G ∗ . As G 𝓁,∗ and I ∗ are determined uniquely from f 2 and f 1 , respectively, and are also strictly positive (if G ∗ is), the claim is shown: system (1) - (3) has a unique positive steady state. 7.2. Two Compartment Model Here we propose a two compartment model of insulin-glucose dynamics; see Section 2 for a discussion of specific expression used in the below system of differential equations. A schematic of the modeled interactions is provided in Figure 10 . Note that in system (32) - (33) we do not incorporate the role of glucagon on stimulating glucose production, and thus there is no inhibition of glucose or insulin in this system (see equation (3) ). Download figure Open in new tab Figure 10: Two compartment insulin-glucose model schematic Download figure Open in new tab Figure 11: The best fit for the insulin-glucose model 32 - 33. The cost function evaluated at the optimal set of parameters is and the adjusted R 2 value is 0.7754. See Table 2 for the corresponding estimated parameters. The results of fitting the insulin and glucose pig data to system (32) - (33) are provided in Figure 11 . Recall that in this case, the cost function that we minimize is given by (11) . For details on the external data and calibration methods, see Section 3 . Tornado plots to ensure minimality are provided in Figure 12 , and the optimal set of parameter values are shown in Table 2 . Download figure Open in new tab Figure 12: Tornado plots for the change in cost function when varying each parameter by the given percentage of its optimal value. The vertical line represents the minimum residual. 7.3. Estimating parameters in the minimal model We provide the best fit solution ( Figure 13 ), estimated parameters ( Table 3 ), and tornado plots ( Figure 14 ) for the minimal model (12) - (14) . For details on the external data and calibration methods, see Section 3 . A comparison between the insulin-glucose-glucagon model and the minmimal model is provided in Figure 6 in the main text of the article. Download figure Open in new tab Figure 13: The best fit for the minimal model (12) - (14) . The cost function evaluated at the optimal set of parameters is . See Table 3 for the corresponding estimated parameters. Download figure Open in new tab Figure 14: Tornado plots for the change in the cost function when varying each parameter by the given percentage of its optimal value. The vertical line represents the minimum residual. References [1]. ↵ 11.2.1 - Bootstrapping Methods | STAT 500. en . URL: https://online.stat.psu.edu/stat500/lesson/11/11.2/11.2.1(visited on 02/10/2024). [2]. ↵ Stephanie A. Amiel . “ The consequences of hypoglycaemia ”. In: Diabetologia 64 . 5 ( 2021 ), pp. 963 – 970 . ISSN: 0012-186X . DOI: 10.1007/s00125-020-05366-3 . URL: https://www.ncbi.nlm.nih.gov/pmc/articles/PMC8012317/ (visited on 02/20/2024). OpenUrl CrossRef PubMed [3]. ↵ Dr Marc Barton . The Discovery of Insulin. en-GB . June 2018 . URL: https://www.pastmedicalhistory.co.uk/the-discovery-of-insulin/ (visited on 01/17/2024). [4]. ↵ RN Bergman , LS Phillips , and C Cobelli . “ Physiologic evaluation of factors controlling glucose tolerance in man: measurement of insulin sensitivity and beta-cell glucose sensitivity from the response to intravenous glucose .” In: Journal of Clinical Investigation 68 . 6 ( Dec . 1981 ), pp. 1456 – 1467 . ISSN: 0021-9738 . URL: https://www.ncbi.nlm.nih.gov/pmc/articles/PMC370948/ (visited on 03/20/2023). OpenUrl [5]. ↵ R. N. Bergman et al. “ Quantitative estimation of insulin sensitivity ”. eng. In: The American Journal of Physiology 236 . 6 ( June 1979 ), E667 – 677 . ISSN: 0002-9513 . DOI: 10.1152/ajpendo.1979.236.6.E667 . OpenUrl CrossRef [6]. ↵ Richard N. Bergman . “ Origins and History of the Minimal Model of Glucose Regulation ”. In: Frontiers in Endocrinology 11 ( 2021 ). ISSN: 1664-2392 . URL : doi: 10.3389/fendo.2020.583016 (visited on 02/20/2024). OpenUrl CrossRef [7]. ↵ Mads Bisgaard Bengtsen and Niels Møller . “ Mini-review: Glucagon responses in type 1 diabetes - a matter of complexity ”. eng. In: Physiological Reports 9 . 16 ( Aug . 2021 ), e15009 . ISSN: 2051-817X . DOI: 10.14814/phy2.15009 . OpenUrl CrossRef [8]. ↵ Linford Briant et al. “ Glucagon secretion from pancreatic -cells ”. In: Upsala Journal of Medical Sciences 121 . 2 ( May 2016 ), pp. 113 – 119 . ISSN: 0300-9734 . DOI: 10.3109/03009734.2016.1156789 . URL: https://www.ncbi.nlm.nih.gov/pmc/articles/PMC4900066/ (visited on 02/10/2024). OpenUrl CrossRef PubMed [9]. ↵ James Cantley and Frances M. Ashcroft . “ Q&A: insulin secretion and type 2 diabetes: why do -cells fail? ” In: BMC Biology 13 ( May 2015 ), p. 33 . ISSN: 1741-7007 . DOI: 10.1186/s12915-015-0140-6 . URL: https://www.ncbi.nlm.nih.gov/pmc/articles/PMC4435650/ (visited on 03/07/2024). OpenUrl CrossRef [10]. ↵ CDC . The Insulin Resistance–Diabetes Connection. en-us . June 2022 . URL: https://www.cdc.gov/diabetes/basics/insulin-resistance.html (visited on 02/20/2024). [11]. ↵ Sebastián Contreras et al. “ A Novel Synthetic Model of the Glucose-Insulin System for Patient-Wise Inference of Physiological Parameters From Small-Size OGTT Data ”. English. In: Frontiers in Bio-engineering and Biotechnology 8 ( Mar . 2020 ). Publisher: Frontiers . ISSN: 2296-4185 . DOI: 10.3389/fbioe.2020.00195 . URL: https://www.frontiersin.org/articles/10.3389/fbioe.2020.00195 (visited on 03/10/2024). OpenUrl CrossRef [12]. ↵ Paul D Docherty et al. “ The Identification of Insulin Saturation Effects During the Dynamic Insulin Sensitivity Test ”. In: The Open Medical Informatics Journal 4 ( July 2010 ), pp. 141 – 148 . ISSN: 1874-4311 . DOI: 10.2174/1874431101004010141 . URL: https://www.ncbi.nlm.nih.gov/pmc/articles/PMC3096059/ (visited on 03/10/2024). OpenUrl CrossRef [13]. ↵ P. Ebeling , H. A. Koistinen , and V. A. Koivisto . “ Insulin-independent glucose transport regulates insulin sensitivity ”. eng. In: FEBS letters 436 . 3 ( Oct . 1998 ), pp. 301 – 303 . ISSN: 0014-5793 . DOI: 10.1016/s0014-5793(98)01149-1 . OpenUrl CrossRef [14]. ↵ Zhuo Fu , Elizabeth R. Gilbert , and Dongmin Liu . “ Regulation of Insulin Synthesis and Secretion and Pancreatic Beta-Cell Dysfunction in Diabetes ”. In: Current diabetes reviews 9 . 1 ( Jan . 2013 ), pp. 25 – 53 . ISSN: 1573-3998 . URL: https://www.ncbi.nlm.nih.gov/pmc/articles/PMC3934755/ (visited on 02/10/2024). OpenUrl [15]. ↵ Ana Gabriela Gallardo-Hernández et al. “ Minimally-Invasive and Efficient Method to Accurately Fit the Bergman Minimal Model to Diabetes Type 2 ”. In: Cellular and Molecular Bioengineering 15 . 3 ( Feb . 2022 ), pp. 267 – 279 . ISSN: 1865-5025 . DOI: 10.1007/s12195-022-00719-x . URL: https://www.ncbi.nlm.nih.gov/pmc/articles/PMC9124285/ (visited on 03/20/2024). OpenUrl CrossRef [16]. ↵ Immacolata Garzilli and Shalev Itzkovitz . “ Design principles of the paradoxical feedback between pancreatic alpha and beta cells ”. en. In: Scientific Reports 8 . 1 ( July 2018 ). Publisher: Nature Publishing Group , p. 10694 . ISSN: 2045-2322 . DOI: 10.1038/s41598-018-29084-4 . URL: https://www.nature.com/articles/s41598-018-29084-4 (visited on 04/05/2024). OpenUrl CrossRef [17]. ↵ C. Gauthier , M. Vranic , and G. Hetenyi . “ Importance of glucagon in regulatory rather than emergency responses to hypoglycemia ”. In: American Journal of Physiology-Endocrinology and Metabolism 238 . 2 ( Feb . 1980 ). Publisher: American Physiological Society , E131 – E140 . ISSN: 0193-1849 . DOI: 10.1152/ajpendo.1980.238.2.E131 . URL: https://journals.physiology.org/doi/abs/10.1152/ajpendo.1980.238.2.e131 (visited on 02/10/2024). OpenUrl CrossRef [18]. ↵ Glucagon - an overview | ScienceDirect Topics . URL: https://www.sciencedirect.com/topics/medicine-and-dentistry/glucagon (visited on 02/10/2024). [19]. ↵ Norikiyo Honzawa , Kei Fujimoto , and Tadahiro Kitamura . “ Cell Autonomous Dysfunction and Insulin Resistance in Pancreatic Cells ”. In: International Journal of Molecular Sciences 20 . 15 ( July 2019 ), p. 3699 . ISSN: 1422-0067 . DOI: 10.3390/ijms20153699 . URL: https://www.ncbi.nlm.nih.gov/pmc/articles/PMC6695724/ (visited on 03/20/2024). OpenUrl CrossRef [20]. ↵ How can I convert from international units (IU) to milligrams or micrograms? | Socratic . URL: https://socratic.org/questions/how-can-i-convert-from-international-units-iu-to-milligrams-or-micrograms (visited on 06/14/2024). [21]. ↵ David Simon Hughes and Parth Narendran . “ Alpha cell function in type 1 diabetes ”. en. In: British Journal of Diabetes 14 . 2 ( June 2014 ). Number: 2, pp. 45 – 51 . ISSN: 2397-6241 . DOI: 10.15277/bjdvd.2014.014 . URL: https://www.bjd-abcd.com/index.php/bjd/article/view/12 (visited on 06/14/2024). OpenUrl CrossRef [22]. ↵ Hyperglycemia in diabetes-Hyperglycemia in diabetes - Symptoms & causes. en . URL: https://www.mayoclinic.org/diseases-conditions/hyperglycemia/symptoms-causes/syc-20373631 (visited on 03/18/2024). [23]. ↵ Insulin and glucagon: How they regulate blood sugar levels. en . Mar . 2019 . URL: https://www.medicalnewstoday.com/articles/316427 (visited on 01/17/2024). [24]. ↵ Insulin and Insulin Resistance: The Ultimate Guide . URL: https://www.healthline.com/nutrition/insulin-and-insulin-resistance (visited on 03/20/2024). [25]. ↵ David A. Jacobson , Barton L. Wicksteed , and Louis H. Philipson . “ The α-Cell Conundrum: ATP-Sensitive K+ Channels and Glucose Sensing ”. In: Diabetes 58 . 2 ( Feb . 2009 ), pp. 304 – 306 . ISSN: 0012-1797 . DOI: 10.2337/db08-1618 . URL: https://www.ncbi.nlm.nih.gov/pmc/articles/PMC2628601/ (visited on 02/10/2024). OpenUrl FREE Full Text [26]. ↵ Ross A. Kelly et al. “ Modelling the effects of glucagon during glucose tolerance testing ”. In: Theoretical Biology and Medical Modelling 16 . 1 ( Dec . 2019 ), p. 21 . ISSN: 1742-4682 . DOI: 10.1186/s12976-019-0115-3 . URL: 10.1186/s12976-019-0115-3(visited on 02/20/2024). OpenUrl CrossRef [27]. ↵ J. E. Liljenquist et al. “Evidence for an important role of glucagon in the regulation of hepatic glucose production in normal man.” en . In: The Journal of Clinical Investigation 59 . 2 ( Feb . 1977 ). Publisher: American Society for Clinical Investigation , pp. 369 – 374 . ISSN: 0021-9738 . DOI: 10.1172/JCI108649 . URL: https://www.jci.org/articles/view/108649 (visited on 02/10/2024). OpenUrl CrossRef PubMed Web of Science [28]. ↵ Asger Lund et al. “ Glucagon and type 2 diabetes: the return of the alpha cell ”. eng. In: Current Diabetes Reports 14 . 12 ( Dec . 2014 ), p. 555 . ISSN: 1539-0829 . DOI: 10.1007/s11892-014-0555-4 . OpenUrl CrossRef PubMed [29]. ↵ Elin Manell et al. “ Establishment of a Refined Oral Glucose Tolerance Test in Pigs, and Assessment of Insulin, Glucagon and Glucagon-Like Peptide-1 Responses ”. en. In: PLOS ONE 11 . 2 ( Feb . 2016 ). Publisher: Public Library of Science , e0148896 . ISSN: 1932-6203 . DOI: 10.1371/journal.pone.0148896 . URL: https://journals.plos.org/plosone/article?id=10.1371/journal.pone.0148896 (visited on 08/25/2023). OpenUrl CrossRef [30]. ↵ Mini-review: Glucagon responses in type 1 diabetes - a matter of complexity - PubMed . URL: https://pubmed.ncbi.nlm.nih.gov/34405569/ (visited on 03/20/2024). [31]. ↵ Jun Sung Moon and Kyu Chang Won . “ Pancreatic -Cell Dysfunction in Type 2 Diabetes: Old Kids on the Block ”. In: Diabetes & Metabolism Journal 39 . 1 ( Feb . 2015 ), pp. 1 – 9 . ISSN: 2233-6079 . DOI: 10.4093/dmj.2015.39.1.1 . URL: https://www.ncbi.nlm.nih.gov/pmc/articles/PMC4342530/ (visited on 03/20/2024). OpenUrl CrossRef PubMed [32]. ↵ MIchelle Mouri and Madhu Badireddy . “ Hyperglycemia ”. eng. In: StatPearls. Treasure Island (FL) : StatPearls Publishing , 2024 . URL: http://www.ncbi.nlm.nih.gov/books/NBK430900/ (visited on 03/18/2024). [33]. ↵ Pasquale Palumbo et al. “ Mathematical modeling of the glucose–insulin system: A review ”. In: Mathematical Biosciences 244 . 2 ( Aug . 2013 ), pp. 69 – 81 . ISSN: 0025-5564 . DOI: 10.1016/j.mbs.2013.05.006 . URL: https://www.sciencedirect.com/science/article/pii/S0025556413001302 (visited on 03/10/2024). OpenUrl CrossRef PubMed [34]. ↵ Pancreas—Islets of Langerhans | SEER Training . URL: https://training.seer.cancer.gov/anatomy/endocrine/glands/pancreas.html (visited on 02/20/2024). [35]. ↵ Allison Petznick . “ Insulin Management of Type 2 Diabetes Mellitus ”. en-US. In: American Family Physician 84 . 2 ( July 2011 ), pp. 183 – 190 . URL: https://www.aafp.org/pubs/afp/issues/2011/0715/p183.html (visited on 02/15/2024). OpenUrl [36]. ↵ Pharmacokinetics common insulin preps - UpToDate . URL: https://www.uptodate.com/contents/image?imageKey=ENDO/130874 (visited on 03/20/2024). [37]. ↵ Plasma Half Life - an overview | ScienceDirect Topics . URL: https://www.sciencedirect.com/topics/biochemistry-genetics-and-molecular-biology/plasma-half-life (visited on 02/10/2024). [38]. ↵ Kenneth R. Feingold et al. Iben Rix et al. “ Glucagon Physiology ”. eng. In: Endotext . Ed. by Kenneth R. Feingold et al. South Dartmouth (MA ): MDText.com, Inc ., 2000 . URL: http://www.ncbi.nlm.nih.gov/books/NBK279127/ (visited on 02/10/2024). [39]. ↵ Subramaniam Sangeetha et al. “ Modeling and Control of the Glucose-Insulin-Glucagon System in Type I Diabetis Mellitus ”. en. In: ChemBioEng Reviews 7 . 3 ( 2020 ). _eprint: https://onlinelibrary.wiley.com/doi/pdf/10.100 pp. 89 – 100 . ISSN: 2196-9744 . DOI: 10.1002/cben.201900015 . URL: https://onlinelibrary.wiley.com/doi/abs/10.1002/cben.201900015 (visited on 04/05/2024). OpenUrl CrossRef [40]. ↵ Amar Singh . Missing Insulin Injections - What To Do & Late Injections. en-GB . Jan . 2019 . URL: https://www.diabetes.co.uk/insulin/missing-insulin-injection.html (visited on 02/15/2024). [41]. ↵ Vijaya Subramanian et al. “ A glucose-insulin-glucagon coupled model of the isoglycemic intravenous glucose infusion experiment ”. In: Frontiers in Physiology 13 ( Sept . 2022 ), p. 911616 . ISSN: 1664-042X . DOI: 10.3389/fphys.2022.911616 . URL: https://www.ncbi.nlm.nih.gov/pmc/articles/PMC9485803/ (visited on 02/20/2024). OpenUrl CrossRef [42]. ↵ Yun Sun et al. “ Delayed insulin secretion response during an OGTT is associated with an increased risk for incidence of diabetes in NGT subjects ”. eng. In: Journal of Diabetes and Its Complications 30 . 8 ( 2016 ), pp. 1537 – 1543 . ISSN: 1873-460X . DOI: 10.1016/j.jdiacomp.2016.07.029 . OpenUrl CrossRef [43]. ↵ Toward Physiological Understanding of Glucose Tolerance: Minimal-Model Approach | Diabetes | American Diabetes Association . URL: https://diabetesjournals.org/diabetes/article/38/12/1512/9005/Toward-Physiological-Understanding-of-Glucose (visited on 03/20/2024). [44]. ↵ Paula M. Trief et al. “ Incorrect Insulin Administration: A Problem That Warrants Attention ”. In: Clinical Diabetes: A Publication of the American Diabetes Association 34 . 1 ( Jan . 2016 ), pp. 25 – 33 . ISSN: 0891-8929 . DOI: 10.2337/diaclin.34.1.25 . URL: https://www.ncbi.nlm.nih.gov/pmc/articles/PMC4714726/ (visited on 03/20/2024). OpenUrl Abstract / FREE Full Text [45]. ↵ Robert D. Utiger . Pancreas | anatomy | Britannica . URL: https://www.britannica.com/science/pancreas (visited on 02/17/2023). [46]. ↵ Senthil K. Venugopal , Parvathy Sankar , and Ishwarlal Jialal . “ Physiology, Glucagon ”. eng. In: Stat-Pearls . Treasure Island (FL ): StatPearls Publishing , 2024 . URL: http://www.ncbi.nlm.nih.gov/books/NBK537082/ (visited on 03/20/2024). [47]. ↵ Diana R Wasag et al. “ Excess all-cause mortality before age 30 in childhood onset type 1 diabetes: data from the Brecon Group Cohort in Wales ”. en. In: Archives of Disease in Childhood 103 . 1 ( Jan . 2018 ), pp. 44 – 48 . ISSN: 0003-9888, 1468-2044 . DOI: 10.1136/archdischild-2016-312581 . URL: https://adc.bmj.com/lookup/doi/10.1136/archdischild-2016-312581 (visited on 02/20/2024). OpenUrl Abstract / FREE Full Text [48]. ↵ Nan Wu and Guohua An . “ A Quantitative Systems Pharmacology Model of the Incretin Hormones GIP and GLP1, Glucagon, Glucose, Insulin, and the Small Molecule DPP-4 Inhibitor, Linagliptin ”. In: Journal of Pharmaceutical Sciences 113 . 1 ( Jan . 2024 ), pp. 278 – 289 . ISSN: 0022-3549 . DOI: 10.1016/j.xphs.2023.09.006 . URL: https://www.sciencedirect.com/science/article/pii/S0022354923003696 (visited on 04/05/2024). OpenUrl CrossRef View the discussion thread. Back to top Previous Next Posted August 01, 2024. Download PDF Email Thank you for your interest in spreading the word about bioRxiv. NOTE: Your email address is requested solely to identify you as the sender of this article. Your Email * Your Name * Send To * Enter multiple addresses on separate lines or separate them with commas. You are going to email the following Evaluating the Importance of Glucagon in the Insulin-Glucose Regulatory System: A Mechanistic Modeling Approach Message Subject (Your Name) has forwarded a page to you from bioRxiv Message Body (Your Name) thought you would like to see this page from the bioRxiv website. Your Personal Message CAPTCHA This question is for testing whether or not you are a human visitor and to prevent automated spam submissions. Share Evaluating the Importance of Glucagon in the Insulin-Glucose Regulatory System: A Mechanistic Modeling Approach Mackenzie Dalton , Emmanuel Asante-Asamani , James Greene bioRxiv 2024.07.31.606029; doi: https://doi.org/10.1101/2024.07.31.606029 Share This Article: Copy Citation Tools Evaluating the Importance of Glucagon in the Insulin-Glucose Regulatory System: A Mechanistic Modeling Approach Mackenzie Dalton , Emmanuel Asante-Asamani , James Greene bioRxiv 2024.07.31.606029; doi: https://doi.org/10.1101/2024.07.31.606029 Citation Manager Formats BibTeX Bookends EasyBib EndNote (tagged) EndNote 8 (xml) Medlars Mendeley Papers RefWorks Tagged Ref Manager RIS Zotero Tweet Widget Facebook Like Google Plus One Subject Area Systems Biology Subject Areas All Articles Animal Behavior and Cognition (7644) Biochemistry (17726) Bioengineering (13916) Bioinformatics (42033) Biophysics (21486) Cancer Biology (18635) Cell Biology (25549) Clinical Trials (138) Developmental Biology (13397) Ecology (19940) Epidemiology (2067) Evolutionary Biology (24361) Genetics (15620) Genomics (22541) Immunology (17763) Microbiology (40468) Molecular Biology (17207) Neuroscience (88739) Paleontology (667) Pathology (2842) Pharmacology and Toxicology (4834) Physiology (7659) Plant Biology (15175) Scientific Communication and Education (2047) Synthetic Biology (4304) Systems Biology (9834) Zoology (2272)