 Research
 Open Access
 Published:
The necessity of identifying the basal glucose setpoint in the IVGTT for patients with Type 2 Diabetes
BioMedical Engineering OnLine volume 14, Article number: 18 (2015)
Abstract
Background
The modelbased dynamic insulin sensitivity and secretion test (DISST) uses fasting glucose (G _{ 0 }) as the basal glucose (G _{ B }) concentration when assessing insulin sensitivity (SI). However, this model was developed in a healthy, normoglycaemic cohort. We sought to determine the suitability the DISST model has for individuals with established type 2 diabetes (T2D).
Methods
14 participants with established T2D were recruited to take part in a dietary intervention study. Insulinmodified intravenous glucose tolerance tests (IMIVGTT) were undertaken at week 0, 12 and 24 and were used with DISST model to identify G _{ B }. A total of 36 tests were conducted across 12 participants throughout the study. Measured G _{ 0 } and identified G _{ B } values were compared using a KolmogorovSmirnov (KS) and signed rank (RS) test for the cohort.
Results
There were significant differences between the G _{ 0 } and identified G _{ B } values in this cohort (p_{rs} and p_{ks} < 0.0001), although both values were well correlated (R = 0.70). The residual plot demonstrates that the modified model captures the behaviour of the participants more accurately than the original model.
Conclusions
This analysis has shown that G _{ B } is an important variable for modelling the glycaemic behaviour in T2D. These findings suggest that the original DISST model, while appropriate for normoglycaemic cohorts, needs to model basal glucose level as a variable for assessing individuals with established T2D.
Background
Type 2 diabetes mellitus (T2DM) is a metabolic disease that affects the body’s ability to regulate glucose concentrations [13]. T2DM is characterized by fasting and postprandial hyperglycaemia [4] and causes comorbidities with significant personal and economic cost [5]. The hyperglycaemia is attributed to a combination of impaired insulin utilization (insulin resistance) and a limited ability to compensate with insulin production (net insulin deficiency).
The ability to quantify both insulin sensitivity and insulin secretion is essential to improving the understanding of the complex physiology underlying type 2 diabetes. Mathematical models of glycaemic dynamics have been coupled with clinical data to identify key aspects in the pathogenesis of type 2 diabetes. The dynamic insulin sensitivity and secretion test (DISST) incorporates a clinical protocol similar to the insulinmodified intravenous glucose tolerance test (IMIVGTT) [6,7] and measures the participant glucose, insulin and Cpeptide responses. The DISST data modelling and data fitting methods were customized to the clinical protocol and allow a robust measurement of insulin sensitivity (SI) that avoids the problems encountered with IVGTT assessment in insulin resistant patients [811]. The DISST SI value is highly correlated to the euglycaemic hyperinsulinaemic clamp (EIC) SI value (R = 0.81), which is widely regarded as the reference method [12].
However, the DISST model uses the participant’s measured fasting glucose concentration (G _{ 0 }) as their basal glucose concentration (G _{ B }). The G _{ B } term in the DISST model effectively has the role of determining the setpoint towards which the modelled glucose response tends, where this choice matches assumptions in all other modelbased tests [6,9,1315]. However, the DISST model was developed in a relatively healthy, normoglycemic cohort [16,17]. In contrast, studies show that G _{ 0 } levels and insulin concentrations are slightly higher in the morning than their overnight “basal” levels, especially for participants with diabetes [1821]. Therefore, the evidence suggests that G _{ B } and G _{ 0 } should be treated as separate entities for individuals with established diabetes as the levels are determined by relative insufficiencies in SI, endogenous insulin secretion (U _{ N }) and rates of gluconeogenesis [2224].
We sought to determine whether a novel modelling approach that identifies G _{ B } as a variable in individuals with established type 2 diabetes more accurately reflects glucose dynamics.
Methods
Participants
Fourteen individuals with established type 2 diabetes mellitus were recruited from the Wellington region of New Zealand to take part in an AtkinsBased low carbohydrate dietary intervention study. Recruited participants were aged between 30 and 65 with a BMI range of 34 to 46 kg · m^{−2} at baseline. Participants were excluded if they had major physiological or psychological illness at the time of testing. Pregnant or lactating females were also excluded. Two participants discontinued the intervention, the first citing personal reasons, and the second left the study due to a renal stone. Twelve participants each underwent three IMIVGTT over a 24 week period resulted in a total of 36 tests for the study. Participants had their age and BMI recorded (median [IQR]; 47.5 [42.5, 54.5] and 40.40 [37.48, 43.48], respectively). Full demographic details and results of the intervention study have been previously described [25]. Ethics approval for this study was provided by the New Zealand Ministry of Health, Central Regional Ethics Committee.
Clinical procedure
The clinical protocol utilised in this study was similar to the protocol defined by Ward et al. [7]. A 0.2 g · kg^{−1} glucose bolus was administered at t = 1 minute and then an infusion of insulin that was intended to replicate the insulinaemic response of a normoglycaemic individual was administered. An insulin infusion was started at t = 2 minutes at a rate of 3.5 mU · kg^{−1} · min^{−1} and was reduced to 0.5 mU · kg^{−1} · min^{−1} at t = 7 minutes. Further reductions occurred at t = 17 minutes, to 0.25 mU · kg^{−1} · min^{−1}, and at t = 50 minutes, to 0.1 mU · kg^{−1} · min^{−1}. The infusion of mU · kg^{−1} · min^{−1} was maintained for the remainder of the procedure. Venous blood samples were taken into fluoride oxalate tubes at times: t = −10, −5, −1, 0, 2, 3, 4, 5, 6, 8, 10, 12.5, 15, 20, 25, 30, 35, 40, 50, 60, 70, 80, 90, 100, 120, 140, 160, 180, 210, 240, 270 and 300 minutes. Blood samples were assayed for glucose and insulin concentration using standard commercial assays at an accredited laboratory.
Model
Dynamic Insulin Sensitivity and Secretion Test (DISST) Model
This analysis used the DISST models of interstitial insulin kinetics and glucose dynamics [16,17]:
where equation nomenclature is shown in Table 1.
V _{ Q }, n _{ I } and n _{ C } are defined apriori based on anatomical functions [17,2628] while DISST model sets p _{ G } as a constant at 0.004 min^{−1} [17].
Parameter identification
Interstitial insulin (Q) was simulated via integrating factors and a linear interpolation of I.
where Q _{ 0 } is determined assuming a steady state at t = −10 minutes:
The DISST model typically sets G _{ B } as equal to G _{ 0 } [16,17]. Hence, G _{ 0 } acts as a surrogate basal glucose concentration level. However, for individuals with elevated fasting glucose, this assumption may not be accurate [1821], and should be tested.
In this analysis, G _{ B } was identified in concert with SI and V _{ G }. The typical approach used with the DISST model identifies only SI and V _{ G }. Thus, the outcomes of the threeparameter (G _{ B }, SI, V _{ G }) model were compared to the outputs of the typical twoparameter model (SI, V _{ G }).
A Gauss Newton parameter identification method was used to identify the participantspecific models, with objective function defined:
where x_{ i } = [G _{ Bi }, SI _{ i }, V _{ Gi }]^{T} and i is the iteration number. The Jacobian matrix (J) and the residual matrix (ψ) are defined:
where n is the number of measured samples, G(x _{ i } ,t _{ 1 }) is the modelled glucose concentration at t = t _{ 1 } given x _{i}, G _{ S }(t _{ 1 }) is the measured glucose level at t = t _{ 1 }.
The Jacobian was numerically evaluated using perturbations of [δG _{ B }, δSI, δV _{ G }] = [10^{− 3}, 10^{8}, 10^{− 3}]. These perturbation values were 0.1% of the order of magnitude of the expected parameter values. Glucose samples between t = 1 and t = 10 minutes were disregarded by the identification methods, as this period is heavily influenced by mixing kinetics that are not captured by the whole body model of glucose metabolism [17,29].
Identifying G _{ B } in concert with SI and V _{ G } can cause identified parameter trade off in some cases [8]. The value of V _{ G } was thus limited to physiologically measured bounds from other studies [16,17,30,31]. In particular, V _{ G } was limited to the range of [0.12Bw, 0.25Bw] where bodyweight (Bw) is measured in kg and the coefficients have units of l · kg^{−1}, which is a standard estimation approach linking volume to an easily measured value. Similarly, G _{ B } was limited to a minimum of 3 mmol · L^{−1}.
Analysis
Model residuals and interpretation of population trends were used to assess the performance of the G _{ B } identified  DISST model. The pvalues are defined with signed ranksum (p_{rs}) and Kolmogrov Smirnov test (p_{ks}) to assess median and variability. All analysis was undertaken using MATLAB (R2013b, Mathworks, Inc., Natick, MA, USA).
Results
Figures 1 and 2 shows the individual relationships between fasting glucose, (G _{ 0 }), identified basal glucose (G _{ BID }) and identified SI from normal DISST and G _{ B }identified DISST model across all participants and tests. Note the bias about the 1:1 line indicating that on average, the identified, modelbased basal set point for glucose was significantly lower than the fasting rate for this cohort with diabetes. Figure 1 shows there were significant differences between the G _{ 0 } and G _{ BID } values in this cohort (Signedranksum: p_{rs} < 0.0001, Kolmogorov Smirnov: p_{ks} < 0.0001). In general, G _{ 0 } was higher than the G _{ BID } value, with only 4 exceptions over 36 results (11.1%). Although there was a significant difference in the levels of G _{ BID } and G _{ 0 }, they were relatively well correlated (R = 0.70), indicating moderately consistent bias in the relationship between values.
Figure 2 shows the effect that identifying G _{ BID } has on identified insulin sensitivity values. There is a reasonably strong correlation between the SI values between normal DISST and G _{ B }identified DISST model (R = 0.83). The bias indicates that by identifying basal glucose, the model captures consistently lower SI values for those with established T2DM.
Figure 3 shows the fitted glucose profiles and measured glucose data from 3 different participants. It also shows that the identified G _{ BID } levels are well below than G _{ 0 } as depicted in Figure 1. Figure 4 illustrates the residual errors of both the typical DISST model and the proposed three parameter identified model that identifies basal glucose. Note that the glucose samples taken within 10 minutes of glucose injection were ignored due to unmodelled mixing effects.
Discussion
This is the first study to demonstrate that glucose excursions are more accurately modelled using basal glucose as the variable in the DISST model, rather than fasting glucose for individuals with established type 2 diabetes. The typical approach employed when using the DISST model defines the fasting glucose (G _{ 0 }) as the basal glucose concentration (G _{ B }). Hence, the glucose response defined by the original model typically tends towards the measured basal value. However, this analysis has shown that this assumption is not valid for a cohort with established diabetes. This discrepancy in the assumption is evidenced by the significant distinction in the values of G _{ 0 } and G _{ B } that reduces as the participant glycaemic control improved across the time points of the dietary intervention study [25]. Figure 1 shows that while most participants had elevated fasting glucose levels, the identified basal level was often much closer to a lower value seen in healthy subjects. In particular, 14 of the 36 G _{ B } values identified were in the normal reference range of 4–5.6 mmol∙L^{−1} [32] while only two of 36 G _{ 0 } values were in that range.
However, there were some participants for whom G _{ BID } remained very high throughout the intervention. Of the three participants that exhibited G _{ BID } values greater than 9 mmol · L^{−1}, two were first diagnosed 10 years prior to the trial. In contrast, the mean duration of diabetes for the whole cohort was 4.4 years (SD = 1.0 year). This outcome indicates a possible mechanism of dysfunction in type 2 diabetes that develops during the course of the disease, and matches growing dysfunction over time in these individuals. However, this study lacks the numbers required for conclusive proof of this trend.
Figure 2 shows the effects of insulin sensitivity (SI) values have from identified G _{ B } values are significantly lower than G _{ 0 } values particularly for these type 2 diabetes participants. Hypothetically, if models to set G _{ B } to be equal to G _{ 0 }, SI will account for low glucose level rather than G _{ B } and is thus modelled as a higher SI value. A recent study shows that the type 2 diabetes subjects have SI values in the magnitude of 2–4 × 10^{−4} L · mU^{−1} · min^{−1} [16]. Although, there is not enough evidence to prove the agreeable range of SI value for type 2 diabetes participants, it is understandable that lower SI value contributes to the pathogenesis of type 2 diabetes as SI is inversely proportional to insulin resistance (IR) [3].
Figure 3 shows the blood glucose profiles of three participants as modelled by the typical twoparameter (SI, V _{ G }) and threeparameter (G _{ B }, SI, V _{ G }) DISST model. While the typical twoparameter DISST model fails to fully capture the responses of the participants, the amended threeparameter model captures the behaviours more closely. This outcome is confirmed by the residual plots in Figure 4 that indicate a much smaller, yet consistent trend about the measured data. This change implies the modified model captures previously unmodelled effects or poor apriori estimates in the interstitial insulin kinetic model.
The original DISST model was developed [17] and validated [12] in relatively normoglucose tolerant cohorts. In these cohorts, the incidence of impaired fasting glucose was relatively low, and thus, the assumption of G _{ 0 } equals to G _{ B } was well founded. However, the glycaemic behaviour of the cohort used in this analysis showed that this assumption was most likely to be invalid. In particular, the lower glucose levels achieved in the later part of the test would be falsely attributed to insulin sensitivity rather than a G _{ B } value that was lower than G _{ 0 }. The significantly biased residuals showed in Figure 4 show that the typical DISST model cannot capture all dynamics of this cohort without identifying G _{ B } directly.
Overall, these results indicate that the original twoparameter approach, while appropriate for normoglycaemic cohorts [16,17], is less suitable for individuals with established type 2 diabetes. Furthermore, it is most likely that the outcomes of this study would be applicable to prediabetic individuals that have elevated blood glucose. However, this assertion remains to be determined.
Neither DISST approach accurately captures the peak value of the measured blood glucose data. This particular result was due to the disregarded glucose data within 10 minutes of glucose injection. This data was rejected due to the unmodelled effects of intravascular mixing [29]. A second compartment to model local/global mixing kinetics could have been added. However, this addition was deemed unnecessary, as such compartments do not add value to the DISST modelled outcomes [17]. The approach used was intended to avoid overfitting to mixing effects or fitting the simple DISST model of global glucose dynamics to the local glucose mixing data.
Although, this analysis was done in a small cohort, the outcomes are significant as it has shown that G _{ B } is an important variable when modelling the glycaemic behaviour in type 2 diabetes. It also showed that G _{ B } can be quite different to the typically assumed G _{ 0 } value, and that it may also have some diagnostic value. These findings suggest that the G _{ B } value should be treated as a variable in DISST model for this cohort. Further validation in a much larger cohort will provide a solid foundation for these findings.
Conclusions
This analysis has shown the presence of a dysfunction in the basal (setpoint) glucose in individuals with type 2 diabetes. The magnitude of the dysfunction has been shown to be linked to insulin sensitivity and the degree of fasting glucose. This analysis suggests that the basal glucose is a more appropriate variable for individuals with type 2 diabetes, as using the fasting glucose measurement as the basal setpoint was shown to be a poor assumption for this cohort  although this requires confirmation in a larger study with a clamp as the reference.
Abbreviations
 .:

Physiological conditions
 T2DM:

Type 2 diabetes mellitus
 .:

Insulin sensitivity tests
 EIC:

Euglycaemic hyperinsulinaemic clamp
 IVGTT:

Intravenous glucose tolerance tests
 IMIVGTT:

Insulinmodified IVGTT
 DISST:

Dynamic insulin sensitivity and secretion test
 .:

Blood assays
 G _{ 0 } :

Fasting glucose
 G _{ B } :

Basal glucose
 G _{ BID } :

Identified basal glucose
 .:

Common model parameters
 SI :

Insulin sensitivity (L · mU^{−1} · min^{−1})
 U _{ N } :

endogenous insulin secretion (mU · min^{−1})
 I :

Plasma insulin concentration (mU · L^{−1})
 G :

Plasma glucose concentration (mmol · L^{−1})
 G _{ B } :

Basal plasma glucose concentration (mmol · L^{−1})
 Q :

Interstitial insulin concentration (mU · L^{−1})
 Q _{ B } :

Basal interstitial insulin concentration (mU · L^{−1})
 V _{ Q } :

Interstitial insulin distribution volume (L)
 n _{ I } :

Plasmainterstitial diffusion rate (L · min^{−1})
 n _{ C } :

Interstitial insulin degradation rate (min^{−1})
 P _{ X } :

Exogenous glucose input rate (mmol · min^{−1})
 p _{ G } :

Noninsulin mediated glucose disposal rate (min^{−1})
 V _{ G } :

Glucose distribution volume (L)
References
 1.
DeFronzo RA, Ferrannini E. Insulin resistance. A multifaceted syndrome responsible for NIDDM, obesity, hypertension, dyslipidemia, and atherosclerotic cardiovascular disease. Diabetes Care. 1991;14:173–94.
 2.
Martin BC, Warram JH, Krolewski AS, Bergman RN, Soeldner JS, Kahn CR. Role of glucose and insulin resistance in development of type 2 diabetes mellitus: results of a 25year followup study. Lancet. 1992;340:925–9.
 3.
Ferrannini E. Insulin resistance is central to the burden of diabetes. Diabetes Metab Rev. 1997;13:81–6.
 4.
DeFronzo RA, Ferrannini E, Koivisto V. New concepts in the pathogenesis and treatment of noninsulindependent diabetes mellitus. Am J Med. 1983;74:52–81.
 5.
Santaguida PL, Balion C, Hunt D, Morrison K, Gerstein H, Raina P, et al. Diagnosis, prognosis, and treatment of impaired glucose tolerance and impaired fasting glucose. Evid Rep Technol Assess. 2005;128:1–11.
 6.
Bergman RN, Ider YZ, Bowden CR, Cobelli C. Quantitative estimation of insulin sensitivity. Am J Physiol. 1979;236:E667–77.
 7.
Ward GM, Walters JM, Barton J, Alford FP, Boston RC. Physiologic modeling of the intravenous glucose tolerance test in type 2 diabetes: a new approach to the insulin compartment. Metab Clin Exp. 2001;50:512–9.
 8.
Docherty P, Chase JG, Lotz T, Desaive T. A graphical method for practical and informative identifiability analyses of physiological models: A case study of insulin kinetics and sensitivity. Biomed Eng Online. 2011;10:39.
 9.
Caumo A, Vicini P, Zachwieja JJ, Avogaro A, Yarasheski K, Bier DM, et al. Undermodeling affects minimal model indexes: insights from a twocompartment model. Am J Physiol. 1999;276:E1171–93.
 10.
Cobelli C, Pacini G, Toffolo G, Sacca L. Estimation of insulin sensitivity and glucose clearance from minimal model  New insights from labeled ivgtt. Am J Physiol. 1986;250:E591–8.
 11.
Krudys KM, Kahn SE, Vicini P. Population approaches to estimate minimal model indexes of insulin sensitivity and glucose effectiveness using full and reduced sampling schedules. Am J Physiol Endocrinol Metab. 2006;291:E716–23.
 12.
McAuley KA, Berkeley JE, Docherty PD, Lotz TF, Morenga LAT, Shaw GM, et al. The dynamic insulin sensitivity and secretion testa novel measure of insulin sensitivity. Metab Clin Exp. 2011;60:1748–56.
 13.
Bergman RN, Phillips LS, Cobelli C. Physiologic evaluation of factors controlling glucose tolerance in man: measurement of insulin sensitivity and betacell glucose sensitivity from the response to intravenous glucose. J Clin Invest. 1981;68:1456–67.
 14.
Bergman RN, Prager R, Volund A, Olefsky JM. Equivalence of the insulin sensitivity index in Man derived by the minimal model method and the euglycemic glucose clamp. J Clin Investig. 1987;79:790–800.
 15.
Boston RC, Stefanovski D, Moate PJ, Sumner AE, Watanabe RM, Bergman RN. MINMOD Millennium: a computer program to calculate glucose effectiveness and insulin sensitivity from the frequently sampled intravenous glucose tolerance test. Diabetes Technol The. 2003;5:1003–15.
 16.
Lotz T. High resolution clinical modelbased assessment of insulin sensitivity. PhD Thesis. University of Canterbury, Christchurch, New Zealand, Department of Mechanical Engineering; 2007.
 17.
Lotz TF, Chase JG, McAuley KA, Shaw GM, Docherty PD, Berkeley JE, et al. Design and clinical pilot testing of the modelbased dynamic insulin sensitivity and secretion test (DISST). J Diabetes Sci Technol. 2010;4:1408–23.
 18.
Holman RR, Turner RC. Diabetes: the quest for basal normoglycaemia. Lancet. 1977;1:469–74.
 19.
Holman RR, Turner RC. Basal normoglycemia attained with chlorpropamide in mild diabetes. Metab Clin Exp. 1978;27:539–47.
 20.
Holman RR, Turner RC. Maintenance of basal plasma glucose and insulin concentrations in maturityonset diabetes. Diabetes. 1979;28:227–30.
 21.
Holman RR, Turner RC. The basal plasma glucose: a simple relevant index of maturityonset diabetes. Clin Endocrinol. 1981;14:279–86.
 22.
Cahill Jr GF. The Banting Memorial Lecture 1971. Physiology of insulin in man. Diabetes. 1971;20:785–99.
 23.
Cahill Jr GF, Ashmore J, Renold AE, Hastings AB. Blood glucose and the liver. Am J Med. 1959;26:264–82.
 24.
Steele R, Bjerknes C, Rathgeb I, Altszuler N. Glucose uptake and production during the oral glucose tolerance test. Diabetes. 1968;17:415–21.
 25.
Krebs JD, Bell D, Hall R, ParryStrong A, Docherty PD, Clarke K, et al. Improvements in glucose metabolism and insulin sensitivity with a Lowcarbohydrate diet in obese patients with type 2 diabetes. J Am Coll Nutr. 2013;31:11–7.
 26.
Van Cauter E, Mestrez F, Sturis J, Polonsky KS. Estimation of insulin secretion rates from Cpeptide levels. Comparison of individual and standard kinetic parameters for Cpeptide clearance. Diabetes. 1992;41:368–77.
 27.
Docherty PD. Evaluation and Development of the Dynamic Insulin Sensitivity and Secretion Test for Numenrous Clinical Applications. Ph.D Thesis. University of Canterbury, New Zealand, Mechanical Engineering; 2011.
 28.
Barrett EJ, Eggleston EM, Inyard AC, Wang H, Li G, Chai W, et al. The vascular actions of insulin control its delivery to muscle and regulate the ratelimiting step in skeletal muscle insulin action. Diabetologia. 2009;52:752–64.
 29.
Edsberg B, Herly D, Hildebrandt P, Kuhl C. Insulin bolus given by sprinkler needle: effect on absorption and glycaemic response to a meal. Br Med J (Clin Res Ed). 1987;294:1373–6.
 30.
Defronzo RA, Tobin JD, Andres R. Glucose clamp technique  method for quantifying insulinsecretion and resistance. Am J Physiol. 1979;237:E214–23.
 31.
Ferrannini E, Mari A. How to measure insulin sensitivity. J Hypertens. 1998;16:895–906.
 32.
ADA. Summary of revisions for the 2010 clinical practice recommendations. Diabetes Care. 2010;33:S3.
Acknowledgements
Ministry of Higher Education, Malaysia and the University of Technology MARA, Penang, Malaysia are acknowledged for providing research funding.
Author information
Affiliations
Corresponding author
Additional information
Competing interests
The authors declare that they have no competing interests.
Authors’ contributions
NAO contributed to the concept design, analysis, and drafting of the manuscript. PDD contributed to the concept design, analysis and edited the manuscript. JGC contributed to the analysis and manuscript editing. JDK and DAB contributed to the manuscript editing and data presentation. None of the authors have any conflict of interest associated with the contents of this manuscript. All authors read and approved the final manuscript.
Rights and permissions
This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly credited. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.
About this article
Cite this article
Othman, N.A., Docherty, P.D., Krebs, J.D. et al. The necessity of identifying the basal glucose setpoint in the IVGTT for patients with Type 2 Diabetes. BioMed Eng OnLine 14, 18 (2015). https://doi.org/10.1186/s1293801500157
Received:
Accepted:
Published:
Keywords
 Basal glucose concentration
 DISST
 Insulin resistance
 Insulin sensitivity