Sage Journals: Discover world-class research

Abstract

Background:

Current mathematical models of postprandial glucose metabolism in people with normal and impaired glucose tolerance rely on insulin measurements and are therefore not applicable in clinical practice. This research aims to develop a model that only requires glucose data for parameter estimation while also providing useful information on insulin sensitivity, insulin dynamics and the meal-related glucose appearance (GA).

Methods:

The proposed glucose-only model (GOM) is based on the oral minimal model (OMM) of glucose dynamics and substitutes the insulin dynamics with a novel function dependant on glucose levels and GA. A Bayesian method and glucose data from 22 subjects with normal glucose tolerance are utilised for parameter estimation. To validate the results of the GOM, a comparison to the results of the OMM, obtained by using glucose and insulin data from the same subjects is carried out.

Results:

The proposed GOM describes the glucose dynamics with comparable precision to the OMM with an RMSE of 5.1 ± 2.3 mg/dL and 5.3 ± 2.4 mg/dL, respectively and contains a parameter that is significantly correlated to the insulin sensitivity estimated by the OMM (r = 0.7) Furthermore, the dynamic properties of the time profiles of GA and insulin dynamics inferred by the GOM show high similarity to the corresponding results of the OMM.

Conclusions:

The proposed GOM can be used to extract useful physiological information on glucose metabolism in subjects with normal glucose tolerance. The model can be further developed for clinical applications to patients with impaired glucose tolerance under the use of continuous glucose monitoring data.

Keywords

glucose-only model insulin sensitivity glucose appearance insulin dynamics

Introduction

Mathematical models are a powerful tool to describe and assess the body’s response to food intake in people with normal glucose tolerance as well as prediabetes and type 2 diabetes mellitus (T2DM). These models typically utilise glucose and insulin data after an oral glucose intake for parameter estimation. They have contributed significantly to the understanding of the metabolic processes responsible for the loss of glycaemic control.^1-4 Despite this success, the application of any of the proposed models in clinical practice, that is, for the diagnosis or treatment of individuals with impaired glucose tolerance, has yet to be seen. This lack of clinical application can mainly be attributed to the high cost, unreliability and dependence on venous access of insulin measurements, prohibiting widespread clinical or ambulatory insulin data collection.^5,6 This paper thus aims to develop a glucose-only model (GOM) that describes postprandial glucose dynamics and provides physiological information while only relying on glucose data for parameter estimation.

Excluding a vast number of GOMs for type 1 diabetes mellitus,⁷ where information on exogenous insulin administration can be used during model identification, a comparatively small number of GOMs applied to subjects with normal and impaired glucose tolerance has been published. A subgroup of these GOMs is based on the description of a harmonic oscillator with an impulse input. While this significantly limits their physiological interpretation, these GOMs have been shown to contain parameters that are dependent on glucose tolerance.^8-11 Other GOMs are based on physiological principles, but can only roughly approximate the postprandial glucose dynamics and have been applied to a very limited number of subjects.^12,13 The main weakness of all mentioned GOMs, however, is that their results have not been validated against the results of a model known to provide accurate physiological information. Specifically, this pertains to insulin sensitivity, insulin dynamics and the meal-related appearance of glucose (GA). To overcome this weakness, this work will develop a new GOM based on and validated by the results of the oral minimal model (OMM) of glucose dynamics, identified from glucose and insulin data.¹⁴ The OMM has been validated by gold-standard reference methods in the past and provides an estimation of insulin sensitivity and GA.^15-17 By identifying the novel GOM and the OMM from data of the same subjects, it is possible to validate and compare both models, particularly with respect to the GOM’s ability to provide physiological information on insulin sensitivity, insulin dynamics and GA.

Methods

Data Description

The dataset used in this work was collected by Ahmed et al.¹⁸ and Nuttall et al.¹⁹ and is publically accessible.²⁰ It contains plasma glucose and insulin profiles from subjects with normal glucose tolerance (NGT) collected over 12 hours in a single day, where subjects consumed three identical meals four hours apart. Blood samples were collected at the same time in each subject after meal consumption at 0, 2, 5, 10, 20, 30, 40, 50, 60 min, then every 15 min up to 120 min and then every 30 min up to 240 min. One additional fasting sample was collected before breakfast, that is, at −15 min.

In this work glucose and insulin profiles from 22 subjects consuming two different meal types of standard (STAND) and high carbohydrate (HCHO) macronutrient composition are used, leading to a total of 66 recorded responses. The average glucose and insulin profiles are shown in Figure 1. The absolute amount of macronutrients provided was scaled according to the body weight and female subjects received 12.5% fewer calories per body weight. Details on the subjects and consumed meals are provided in Table 1.

Figure 1.

Mean and standard deviation (shaded areas) of the glucose and insulin profiles above basal levels for the two meal types of (a) standard (STAND) and (b) high carbohydrate (HCHO) composition utilised in this paper. The basal level is calculated for each subject individually as the average of the −15, 0, 2 and 5 min measurement points. The vertical dashed lines indicate the time of meal consumption.

Table 1.

Details on the Subject Characteristics and Different Meal Types Containing Standard (STAND) and High Carbohydrate (HCHO) Mixtures of Macronutrient Content.

	STAND	HCHO
Number of subjects (females)	12 (5)	10 (4)
Age	23 ± 1	25 ± 3
Body weight males (females) [kg]	76 ± 5 (59 ± 1)	77 ± 4 (59 ± 5)
Meal composition [% CHO/Fat/Protein]	40/49/11	63/27/10
CHO per meal (females) [g/kg body weight]	1.2 (1.1)	2 (1.8)
Calories per meal (females) [kcal/kg body weight]	13 (11)	13 (11)

The meal composition is given in percentage of calories contained in the respective macronutrient content. Data are given as mean ± standard error.

Model Formulation

The GOM proposed in this work is based on the following generalised formulation of the OMM¹⁴:

$\begin{array}{l} \frac{d G (t)}{d t} = - G (t) X (t) - p_{1} [G (t) - G_{b}] \\ + \frac{R a_{P L} (t) + R a p (t)}{V} G (0) = G_{0}, \end{array}$ (1)

$\frac{d X (t)}{d t} = - p_{2} (X (t) - S_{I} [I (t) - I_{b}]), X (0) = X_{0},$ (2)

The glucose concentration, its basal level and initial condition are represented by $G$ , $G_{b}$ and $G_{0}$ (mg/dL), respectively. Parameters $p_{1}$ (min⁻¹) and $V$ (dL/kg) represent the glucose effectiveness and distribution volume of glucose relative to body weight, respectively. The state $X$ (min⁻¹) and its initial condition $X_{0}$ represent the insulin action in a remote compartment with parameter $p_{2}$ (min⁻¹) governing its decay dynamics and $S_{I}$ (min⁻¹ per mU/L) representing insulin sensitivity. The insulin concentration $I$ (mU/L) and its basal level $I_{b}$ are considered to be known, error-free inputs. The input function $R a_{P L}$ (mg/kg/min) describes the meal-related, posthepatic GA and is described by a piecewise linear function with seven breakpoints at adjustable heights and a fixed area under the curve (AUC), calculated based on the carbohydrate content of the meal (Figure 2). The function $R a p$ represents the persisting GA originating from a previously consumed meal. The measurement process of the glucose levels is considered to be affected by an additive, normally distributed error with zero mean and a known standard deviation. The unknown parameters to be estimated from glucose and insulin data are $p_{1}$ , $p_{2}$ , $S_{I}$ and the adjustable heights of GA function $R a_{P L}$ . The details of the model and parameter estimation procedure have been described previously.¹⁴

Figure 2.

Example of the piecewise linear GA function $R a_{P L}$ used in the oral minimal model (OMM) and the log-normally based GA function $R a_{L N}$ used in the glucose-only model (GOM) with associated 95% confidence intervals (shaded area). The confidence intervals overlap for 75% (OL value) of the response duration of 240 min. The dashed lines indicate the two components of $R a_{L N}$ .

To formulate a GOM based on the OMM, it is necessary to remove the measured insulin levels, that is, $I$ and $I_{b}$ as a known input. For that, the following model is proposed:

$\begin{array}{l} \frac{d G (t)}{d t} = - G (t) X (t) - p_{1} [G (t) - G_{b}] \\ + \frac{R a_{L N} (t) + R a p (t)}{V} G (0) = G_{0}, \end{array}$ (3)

$\frac{d X (t)}{d t} = - p_{2} [X (t) - S_{G} Z (t)], X (0) = X_{0},$ (4)

$Z (t) = \frac{G (t) - G_{b}}{\underset{Z_{P O S}}{\underset{︸}{1 + \exp [- α (G (t) - G_{b})]}}} + β \frac{R a_{L N} (t)}{V}$ (5)

where the GA function $R a_{L N}$ is defined as follows:

$\begin{matrix} R a_{L N} (t) = A (1 - R_{H}) f_{L N} (t, T_{1}, W_{1}) \\ + A R_{H} f_{L N} (t, T_{2}, W_{2}), \end{matrix}$ (6)

$f_{L N} (t, T, W) = {\begin{matrix} 0 & i f t = 0 \\ \frac{1}{t \sqrt{π W}} \exp [\frac{- {(\log (\frac{t}{T}) - \frac{W}{2})}^{2}}{W}] & i f t > 0 . \end{matrix}$ (7)

The process of observing the glucose levels is considered to be identical to the OMM (details in section 1.2 of the supplementary information). Furthermore, the parameters $p_{1}$ , $p_{2}$ , $G_{b}$ , $G_{0}$ and $V$ as well as the variables $G$ and $R a p$ from expressions (3) to (4) keep the same interpretation as in the OMM. Instead of the piecewise linear GA function $R a_{P L}$ , the GOM features the fully differentiable function $R a_{L N}$ in expressions (6) to (7) that is based on two overlapping components defined by a modified structure of the log-normal distribution (Figure 2). The function has a total fixed AUC of $A$ (mg/kg) which is calculated from the carbohydrate content of the meals and has parameters $T_{1}$ (min), $T_{2}$ (min), $W_{1}$ and $W_{2}$ governing the peak times and general widths of the respective components. The parameter $R_{H}$ is restricted to the range (0,1) which ensures positivity of $R a_{L N}$ and determines the contributions of each component to the total AUC. The function $R a_{L N}$ has been suggested previously as a replacement for the piecewise linear function $R a_{P L}$ in the context of the OMM, where additional details on the function can be found.¹⁴

The main adaptation of the GOM (3) to (7) in comparison to the OMM (1) to (2) is the introduction of the variable $Z$ (mg/dL) in place of the insulin profile above baseline $[I - I_{b}]$ . This adaptation results in the fact that the state $X$ (min⁻¹) and its initial condition $X_{0}$ in expression (4) of the GOM no longer represent the active insulin. Instead, the state $X$ is interpreted as a general glucose-lowering effect and the parameter $S_{G}$ (min⁻¹ per mg/dL) replaces the insulin parameter $S_{I}$ in the OMM. It is thus expected that the parameter $S_{G}$ contains similar information as the insulin parameter $S_{I}$ .

The formulation of the variable $Z$ in expression (5) is based on the general similarity between glucose and insulin dynamics, especially during the initial rise of a meal response, as demonstrated in Figure 1. This similarity allows the assumption that the information contained in the insulin data can be partially recovered from the glucose data. A similar supposition is made in several models of insulin secretion,^21-24 where glucose levels are considered to be a known input and the primary driver of insulin secretion and therefore insulin levels. Despite the similarity between glucose and insulin dynamics, Figure 1 reveals two main differences that need to be considered by the GOM.

Firstly, it is far more prevalent for glucose levels to fall below the basal level $G_{b}$ than it is for insulin levels to fall below $I_{b}$ . This effect is incorporated by using the function $Z_{P O S}$ (mg/dL) in expression (5), shown in Figure 3. The function has an approximately linear relationship to its input $G (t)$ for $G (t) > G_{b}$ , but approaches zero for $G (t) < G_{b}$ , with the parameter $α$ (dL/mg) governing the shape of the transition (Figure 3). Secondly, the average glucose and insulin profiles above baseline in Figure 1 indicate that, after a simultaneous rise, insulin levels often remain elevated for longer and decay slower toward the baseline levels in comparison to glucose levels. This is especially prominent in all responses to the STAND meal and after breakfast in the HCHO meal. The GOM accommodates this behaviour by including the GA function $R a_{L N}$ in the description of $Z$ . This allows the variable $Z$ to remain elevated even when glucose levels have reached basal levels and the contribution of $Z_{P O S}$ vanishes. The parameter $β$ (min) acts as a unit conversion factor and adjusts the strength of the coupling between $R a_{L N}$ and the variable $Z$ . A comparable feature is also included in the previously mentioned models of insulin secretion,^22-24 where the rate of change of glucose levels, which is highly dependent on $R a_{L N}$ as indicated by expression (3), is thought to affect insulin secretion and therefore its levels.

Figure 3.

Example of the function $Z_{P O S}$ for varying values of the shape parameter $α$ . Overlayed are the simulatenously measured glucose and insulin samples above baseline to illustrate their relationship approximated by the function $Z_{P O S}$ .

Parameter Estimation

The dataset contains three consecutive meal responses from each subject that are considered separately during parameter estimation in the GOM, that is, one set of unknown parameters is estimated from every meal response. To incorporate the overlapping effects of consecutive meals, the parameter estimation procedure previously described for the OMM is utilised.¹⁴ The procedure adapts the initial conditions of the states, $G_{0}$ and $X_{0}$ , as well as the persisting GA $R a p$ based on the time of meal consumption, while keeping the basal level of glucose $G_{b}$ constant throughout the day (details in section 1.1 of the supplementary information).

The following parameters of the GOM (3) to (7) are considered for estimation: system parameters $p_{1}$ , $p_{2}$ , $S_{G}$ and $β$ , and parameters $T_{1}$ , $T_{2}$ , $W_{1}$ , $W_{2}$ and $R_{H}$ governing the log-normally based GA function $R a_{L N}$ . Using the observability rank criterion,^25,26 it can be shown that these parameters are structurally locally identifiable (details in section 1.3 of the supplementary information). The shape parameter $α$ of $Z_{P O S}$ is fixed to a value of 0.1 dL/mg as a stochastic sensitivity analysis revealed that it is practically unidentifiable, that is, it cannot be estimated with an acceptable level of precision (details in section 1.4 of the supplementary information). The value of 0.1 dL/mg is chosen as it approximates the relationship between glucose and insulin data suitably (Figure 3).

The parameter estimation is carried out using a variational Bayesian approach,^27-29 which has been used previously to identify low dimensional models including the OMM.^10,14,30-32 This approach provides a probabilistic treatment of unknown parameters which allows the estimation of parameter uncertainty and requires the specification of prior distributions over unknown parameters. All unknown parameters of the GOM are specified as log-normally distributed and characterised by their median and coefficient of variation (CV) since the parameters are only physiologically plausible when positive. The exception to that is the parameter $R_{H}$ which is restricted to the range (0,1). The details of the chosen prior distributions are provided in section 1.4 of the supplementary information. For the parameters $p_{1}$ and $p_{2}$ as well as the GA function parameters, the same prior distributions as in the OMM are used.¹⁴ For the newly introduced parameters $S_{G}$ and $β$ , a stochastic sensitivity analysis was carried out to ensure that the chosen prior distributions can capture the variabitly of response from the data (details in section 1.4 of the supplementary information). Here, it should be mentioned that due to the chosen GOM formulation, parameters $β$ and $S_{G}$ show significant covariance, which leads to poor estimation precision when both parameters have wide prior distributions. Due to the importance of the parameter $S_{G}$ for carrying information on insulin sensitivity, its prior CV was chosen to be 50%, while simultaneously using a narrow prior distribution (CV of 10 %) for the parameter $β$ .

Validation

The validity of the results produced by the GOM is assessed by comparing the corresponding results of the OMM obtained with the identical approach from the same dataset.¹⁴ In particular, the following aspects are compared between the OMM and the GOM:

Model fit as assessed through the time profile of residuals between the model-inferred and observed glucose levels, weighted by the measurement error and the root mean squared residuals (RMSE).

Information on insulin sensitivity as assessed through correlation and comparison between the parameter $S_{G}$ of the GOM and parameter $S_{I}$ of the OMM.

Agreement of the inferred time profiles of GA, that is, the piecewise-linear GA function $R a_{P L}$ of the OMM and the log-normally based GA function $R a_{L N}$ . To quantify this agreement for every response, the confidence interval (CI) associated with the individual GA profiles is inferred from the posterior distributions of the unknown parameters of $R a_{P L}$ and $R a_{L N}$ . Subsequently, the time during which the 95 % CIs of $R a_{P L}$ and $R a_{L N}$ overlap (OL) is calculated and expressed as the share of the total response time of 240 min (Figure 2). High OL percentages thus indicate high similarity between $R a_{P L}$ and $R a_{L N}$ .

Agreement of the time course of insulin dynamics represented by $S_{I} [I (t) - I_{b}]$ in expression (2) for the OMM, abbreviated as $Y_{O M M}$ , and inferred by $S_{G} Z (t)$ in expression (4) by the GOM, abbreviated as $Y_{G O M}$ . As these quantities enter the description of the state $X$ at the same position in both OMM and GOM, $Y_{G O M}$ could carry information on insulin dynamics. Analogous to the GA profiles, the agreement of $Y_{O M M}$ and $Y_{G O M}$ is quantified as the share of time during with the 95 % CIs overlap. To calculate the CI for $Y_{O M M}$ , the reported insulin assay CV of 13 % is used.¹⁸

Results

The parameter estimation of the GOM was carried out in all 66 recorded responses (22 subjects with 3 responses each), and the individual results are provided in section 2.1 of the supplementary information. The time profile of weighted residuals between model-inferred and observed glucose levels is displayed in Figure 4. These results demonstrate that the model is capable of describing the glucose data well, as all average weighted residuals are contained within the −1/+1 range. Additionally, it is demonstrated that in comparison to the OMM results, the GOM shows a smaller error in the first 30 min of the responses. The RMSE values of the GOM are statistically equivalent to the RMSE values of the OMM, that is, 5.1 ± 2.3 mg/dL for the GOM and 5.3 ± 2.4 mg/dL for the GOM (P = .73).

Figure 4.

Mean and standard deviation of weighted residuals between the model-inferred and observed glucose levels for the oral minimal model (OMM) and the glucose-only model (GOM) identified on the same dataset.

The comparison between parameter $S_{G}$ of the GOM and parameter $S_{I}$ of the OMM are displayed in Figure 5. Firstly, the parameter $S_{G}$ can be estimated with good precision as indicated by a small posterior CV of 9.0 ± 2.5 % which is comparable to the posterior CV of $S_{I}$ (6.8 ± 5.1 %). The values of $S_{G}$ and $S_{I}$ are significantly correlated with r = 0.59 and r = 0.7 for the STAND and HCHO meals, respectively (Figure 5a). The posterior results of all parameters are provided in section 2.2 of the supplementary material.

Figure 5.

Results of the parameter $S_{G}$ from glucose-only model (GOM) and the corresponding insulin sensitivity parameter $S_{I}$ from oral minimal model (OMM). The results are compared between the meals of standard (STAND) and high carbohydrate (HCHO) composition. Plot (a) shows a correlation and linear regression analysis with the horizontal and vertical lines for each data point indicating the one-sigma range of the posterior log-normal parameter distribution. Plots (b, c) show boxplots of $S_{I}$ and $S_{G}$ separated by meal type and subject sex.

GA profiles from the OMM ( $R a_{P L}$ ) and the GOM ( $R a_{L N}$ ) are presented in Figure 6. The average profiles in plots (a) and (b) display similar dynamic properties, that is, the shoulder of GA in the STAND meal and secondary peak in the HCHO meal are correctly inferred by the GOM. There is, however, a larger difference in the first 30 min of the response due to the different mathematical formulations of the GA functions. The distribution of OL values in Figure 6c indicates no difference between meal types (P = .63) and shows that a majority of OL values lie above 65 %, indicating a good agreement between inferred GA profiles.

Figure 6.

Comparison between the glucose appearance (GA) profiles estimated by the piecewise linear function $R a_{P L}$ in the oral minimal model (OMM) and the log-normally based function $R a_{L N}$ in the glucose-only model (GOM). Plot (a) displays to the results from the meal of standard composition (STAND) and plot (b) to the meal of high carbohydrate composition (HCHO). The results are given as mean (solid line) and standard deviation (shaded area) of all responses. Plot (c) gives boxplots of the share of time during which the 95% confidence intervals (CIs) of $R a_{P L}$ and $R a_{L N}$ overlap (OL values).

Analogous to the GA profiles, the time courses of $Y_{O M M}$ and $Y_{G O M}$ are compared in Figure 7. The average profiles in plots (a) and (b) show similar dynamic properties, e.g. a shoulder after the initial rise in the case of the HCHO meal. This agreement is confirmed by the distribution of OL values in Figure 7c. Of note is that the OL values of the HCHO meal are increased in comparison to the STAND meal (P = .08), which could be connected to the increased correlation between $S_{I}$ and $S_{G}$ estimates in the HCHO meal.

Figure 7.

Comparison between the estimated profiles of $Y_{O M M} = S_{I} [I (t) - I_{b}]$ in the oral minimal model (OMM) and $Y_{G O M} = S_{G} Z (t)$ in the glucose-only model (GOM). Plot (a) shows to the results from the meal of standard composition (STAND) and plot (b) to the meal of high carbohydrate composition (HCHO). The results are given as mean (solid line) and standard deviation (shaded area) of all responses. Plot (c) gives boxplots of the share of time during which the 95% confidences intervals (CIs) of $Y_{O M M}$ and $Y_{G O M}$ overlap (OL values).

Discussion

A glucose-based model to describe postprandial glucose responses from different meals in subjects with NGT is presented. This new GOM has been formulated and validated using the physiology-based OMM. Analysing the weighted residuals (Figure 4) and RMSE, it can be concluded that the GOM can describe the glucose data equally well and possess sufficient flexibility to account for the large variability in the responses (Figure 1).

The ability of the GOM to provide information on insulin sensitivity is indicated by a significant correlation between the parameters $S_{G}$ and $S_{I}$ (Figure 5a). Especially for the HCHO meal, the correlation coefficient of 0.7 is comparable to commonly used surrogate indices of insulin sensitivity such as HOMA-IR and the Matsuda index, where correlation values of 0.65 and 0.73, respectively, against clamping results have been reported.³³ Further evidence for the informative value of the parameter $S_{G}$ is given by the fact that it displayes the same differences between male and female subjects for the different meal types, as the parameter $S_{I}$ (Figure 5b and 5c). In contrast, the interpretability of $S_{G}$ values across meal types is weakened by the fact that there is a significant difference between meal types (P = .04) that is not observed in $S_{I}$ values (P = .45).

There are two inherent limitations in the approach to using only glucose data to assess insulin sensitivity. Firstly, it could be rarely the case that, the dynamic properties of glucose and insulin levels, e.g. the timing and existence of peaks, can exhibit very little similarity, thus violating one of the modeling assumptions. The second limitation stems from the fact that absolute levels of insulin are not always correlated to absolute glucose levels, even when the dynamical properties of both signals are identical. This means that two subjects could have quantitatively similar glucose profiles but exhibit vastly different absolute insulin levels and thus also have different insulin sensitivities. Detecting this difference using glucose data alone is thus an inherent limitation.

In terms of GA, the results show that average profiles inferred by the GOM and OMM show very similar dynamic properties, with a larger difference in the first 30 min of the response (Figure 6). As the weighted residuals of the GOM are closer to zero in that same period (Figure 4), a more realistic estimation of GA with the log-normally based function $R a_{L N}$ and the GOM during this period is indicated. A very similar observation was made when $R a_{L N}$ was used in conjunction with the OMM.¹⁴

Similar to GA, the GOM’s ability to infer information on insulin dynamics is demonstrated by the similarity of average profiles of $Y_{O M M}$ and $Y_{G O M}$ (Figure 7). To assess the agreement between the individual results of GA and insulin dynamics, the OL value was introduced. The results of both GA and insulin dynamics indicate a satisfactory agreement but also exhibit high variability between responses (see Figures 6c and 7c as well as the individual results in the supplementary information). This variability in OL values reflects the variability in overall glucose and insulin responses (Figure 1). The interpretability of GOM results is thus less reliable on an individual level.

A general weakness of the current study is the use of a dataset that only contains subjects with NGT. To assess the model’s applicability in patients with prediabetes and T2DM, further validation and adaptation with appropriate datasets, e.g. from Peter et al.³⁴ is required.

While the dataset used in this research contained glucose data from blood sampling collected in a controlled clinical setting, it would also be possible to identify the proposed GOM from more easily obtainable, ambulatory datasets. For instance, glucose profiles recorded with continuous glucose monitoring (CGM) at home, where meals are typically consumed at irregular intervals and contain varying amounts of carbohydrates, could be used. An application of the GOM to these types of datasets is in part possible as the GOM features the differentiable input function $R a_{L N}$ that is independent of the considered response duration and easily adaptable to meals with greatly varying carbohydrate content.

Conclusion

This paper, for the first time, proposed a glucose-based model for the successful extraction of useful physiological information on glucose metabolism in subjects with NGT, thereby overcoming the weaknesses of existing GOM approaches.^8-13 The model’s independence from insulin measurements and exclusive use of easily accessible data enable further developments and its potential application in research and clinical practice to a large number of subjects. In particular, the proposed model could allow a more sophisticated physiological interpretation of CGM profiles collected under ambulatory conditions. It could thus support the design of personalised dietary interventions in prediabetes and T2DM or examine the glycaemic derangement in gestational diabetes mellitus.

Supplemental Material

sj-docx-1-dst-10.1177_19322968211026978 – Supplemental material for A Glucose-Only Model to Extract Physiological Information from Postprandial Glucose Profiles in Subjects with Normal Glucose Tolerance

Supplemental material, sj-docx-1-dst-10.1177_19322968211026978 for A Glucose-Only Model to Extract Physiological Information from Postprandial Glucose Profiles in Subjects with Normal Glucose Tolerance by Manuel M. Eichenlaub, Natasha A. Khovanova, Mary C. Gannon, Frank Q. Nuttall and John G. Hattersley in Journal of Diabetes Science and Technology

Footnotes

Abbreviations

AUC,area under the curve;CGM,continuous glucose monitoring;CI,confidence interval;CV,coefficient of variation;GA,meal-related glucose appearance;GOM,glucose-only model;HCHO,meal of high carbohydrate macronutrient composition;NGT,normal glucose tolerance;OL,overlap;OMM,oral minimal model;RMSE,root mean square error;STAND,meal of standard macronutrient composition;T2DM,type 2 diabetes mellitus.

Declaration of Conflicting Interests

The author(s) declared the following potential conflicts of interest with respect to the research,authorship,and/or publication of this article: A preliminary version of this work was presented as a poster at the 2019 Diabetes Technology Meeting and has been published as an abstract in the Journal of Diabetes Science and Technology 14 (2) entitled “Modelling of Glucose Dynamics and Estimation of Insulin Sensitivity from Glucose Data Only”.

Funding

The author(s) disclosed receipt of the following financial support for the research,authorship,and/or publication of this article: The study was jointly funded by the School of Engineering at the University of Warwick and the CRF Human Metabolic Research Unit at the University Hospitals Coventry and Warwickshire NHS Trust. Additional support came from the EPSRC UK grant EP/T013648/1.

ORCID iDs

Manuel M. Eichenlaub

Natasha A. Khovanova

Supplemental Material

Supplemental material for this article is available online.

References

Cobelli

Dalla Man

Pedersen

Bertoldo

Toffolo

Advancing our understanding of the glucose system via modeling: a perspective. IEEE Trans Biomed Eng. 2014;61(5):1577-1592. doi:10.1109/TBME.2014.2310514

Palumbo

Ditlevsen

Bertuzzi

De Gaetano

Mathematical modeling of the glucose-insulin system: a review. Math Biosci. 2013;244(2):69-81. doi:10.1016/j.mbs.2013.05.006

Ajmera

Swat

Laibe

Le Novere

Chelliah

The impact of mathematical modeling on the understanding of diabetes and related complications. CPT Pharmacometrics Syst Pharmacol. 2013;2(7):54. doi:10.1038/psp.2013.30

Rathee

SN.

ODE models for the management of diabetes: a review. Int J Diabetes Dev Ctries. 2017;37(1):4-15. doi:10.1007/s13410-016-0475-8

Staten

Stern

Miller

Steffes

Campbell

Insulin Standardization Workgroup. Insulin assay standardization: leading to measures of insulin sensitivity and secretion for practical clinical care. Diabetes Care. 2010;33(1):205-206. doi:10.2337/dc09-1206

Hoerber

Achenbach

Schleicher

Peter

Harmonization of immunoassays for biomarkers in diabetes mellitus. Biotechnol Adv. Published online February 23, 2019. doi:10.1016/j.biotechadv.2019.02.015

Oviedo

Vehi

Calm

Armengol

A review of personalized blood glucose prediction strategies for T1DM patients. Int J Numer Method Biomed Eng. 2017;33(6):e2833. doi:10.1002/cnm.2833

Ackerman

Rosevear

McGuckin

WF.

A mathematical model of the glucose-tolerance test. Phys Med Biol. 1964;9(2):203.

Khovanova

Zhang

Holt

. Generalised stochastic model for characterisation of subcutaneous glucose time series. In: IEEE-EMBS International Conference on Biomedical and Health Informatics (BHI), Valencia, Spain; 2014:484-487. doi:10.1109/BHI.2014.6864408

10.

Zhang

Holt

Khovanova

A data driven nonlinear stochastic model for blood glucose dynamics. Comput Methods Programs Biomed. 2016;125:18-25. doi:10.1016/j.cmpb.2015.10.021

11.

Vargas

Moreles

Peña

Monroy

Alavez

Estimation and SVM classification of glucose-insulin model parameters from OGTT data: a comparison with the ADA criteria. Int J Diabetes Dev Ctries. Published online September 12, 2020. doi:10.1007/s13410-020-00851-2

12.

Chen

C-L

Tsai

H-W.

Modeling the physiological glucose insulin system on normal and diabetic subjects. Comput Methods Programs Biomed. 2010;97(2):130-140. doi:10.1016/j.cmpb.2009.06.005

13.

Goel

Parkhi

Barua

Shah

Ghaskadbi

A minimal model approach for analyzing continuous glucose monitoring in type 2 diabetes. Front Physiol. 2018;9:673. doi:10.3389/fphys.2018.00673

14.

Eichenlaub

Hattersley

Gannon

Nuttall

Khovanova

NA.

Bayesian parameter estimation in the oral minimal model of glucose dynamics from non-fasting conditions using a new function of glucose appearance. Comput Methods Programs Biomed. 2021;200:105911. doi:10.1016/j.cmpb.2020.105911

15.

Dalla Man

Caumo

Cobelli

. The oral glucose minimal model: estimation of insulin sensitivity from a meal test. IEEE Trans Biomed Eng. 2002;49(5):419-429. doi:10.1109/10.995680

16.

Dalla Man

Caumo

Basu

Rizza

Toffolo

Cobelli

. Minimal model estimation of glucose absorption and insulin sensitivity from oral test: validation with a tracer method. Am J Physiol Endocrinol Metab. 2004;287(4):E637-E643. doi:10.1152/ajpendo.00319.2003

17.

Dalla Man

Yarasheski

Caumo

, et al Insulin sensitivity by oral glucose minimal models: validation against clamp. Am J Physiol Endocrinol Metab. 2005;289(6):E954-E959. doi:10.1152/ajpendo.00076.2005

18.

Ahmed

Gannon

Nuttall

FQ.

Postprandial plasma glucose, insulin, glucagon and triglyceride responses to a standard diet in normal subjects. Diabetologia. 1976;12(1):61-67.

19.

Nuttall

Gannon

Wald

Ahmed

Plasma glucose and insulin profiles in normal subjects ingesting diets of varying carbohydrate, fat, and protein content. J Am Coll Nutr. 1985;4(4):437-450.

20.

Eichenlaub

Hattersley

Gannon

Nuttall

Khovanova

. Data for Bayesian parameter estimation in the oral minimal model of glucose dynamics from non-fasting conditions using a new function of glucose appearance. Published January 5, 2021. Accessed March 4, 2021. https://wrap.warwick.ac.uk/146758/

21.

Hovorka

Chassin

Luzio

Playle

Owens

DR.

Pancreatic beta-cell responsiveness during meal tolerance test: model assessment in normal subjects and subjects with newly diagnosed noninsulin-dependent diabetes mellitus. J Clin Endocrinol Metab. 1998;83(3):744-750. doi:10.1210/jcem.83.3.4646

22.

Ruan

Thabit

Wilinska

Hovorka

Modelling endogenous insulin concentration in type 2 diabetes during closed-loop insulin delivery. Biomed Eng Online. 2015;14:19. doi:10.1186/s12938-015-0009-5

23.

Mari

Schmitz

Gastaldelli

Oestergaard

Nyholm

Ferrannini

Meal and oral glucose tests for assessment of β-cell function: modeling analysis in normal subjects. Am J Physiol Endocrinol Metab. 2002;283(6):E1159-E1166. doi:10.1152/ajpendo.00093.2002

24.

Breda

Cavaghan

Toffolo

Polonsky

Cobelli

Oral glucose tolerance test minimal model indexes of beta-cell function and insulin sensitivity. Diabetes. 2001;50(1):150-158.

25.

Villaverde

Barreiro

Papachristodoulou

Structural identifiability of dynamic systems biology models. PLOS Comput Biol. 2016;12(10):e1005153. doi:10.1371/journal.pcbi.1005153

26.

Villaverde

Evans

Chappell

Banga

JR.

Input-dependent structural identifiability of nonlinear systems. IEEE Control Syst Lett. 2019;3(2):272-277. doi:10.1109/LCSYS.2018.2868608

27.

Daunizeau

Friston

Kiebel

SJ.

Variational Bayesian identification and prediction of stochastic nonlinear dynamic causal models. Physica D. 2009;238(21):2089-2118. doi:10.1016/j.physd.2009.08.002

28.

Daunizeau

Adam

Rigoux

VBA: a probabilistic treatment of nonlinear models for neurobiological and behavioural data. PLoS Comput Biol. 2014;10(1):e1003441. doi:10.1371/journal.pcbi.1003441

29.

Eichenlaub

. On the relationship between a Gamma distributed precision parameter and the associated standard deviation in the context of Bayesian parameter inference. ArXiv210106289 Cs Stat. Published January 15, 2021. Accessed February 12, 2021. http://arxiv.org/abs/2101.06289

30.

Zhang

Briggs

Lowe

, et al A new data-driven model for post-transplant antibody dynamics in high risk kidney transplantation. Math Biosci. 2017;284:3-11. doi:10.1016/j.mbs.2016.04.008

31.

Eichenlaub

Hattersley

Khovanova

. A minimal model approach for the description of postprandial glucose responses from glucose sensor data in diabetes mellitus. In: Conference Proceedings: IEEE Engineering in Medicine and Biology Society; 2019:265-268. doi:10.1109/EMBC.2019.8857195

32.

Eichenlaub

Khovanova

Hattersley

. A model describing the multiphasic dynamics of mixed meal glucose responses in healthy subjects. In: IFMBE Proceedings. Springer Singapore; 2019:577-581.

33.

Bonadonna

Boselli

Dei Cas

Trombetta

. Methods to assess in vivo insulin sensitivity and insulin secretion. In: Bonora

DeFronzo

, eds. Diabetes Epidemiology, Genetics, Pathogenesis, Diagnosis, Prevention, and Treatment. Springer International Publishing; 2018:317-367.

34.

Peter

Dunseath

Luzio

Chudleigh

Roy Choudhury

Owens

DR.

Daytime variability of postprandial glucose tolerance and pancreatic B-cell function using 12-h profiles in persons with Type 2 diabetes. Diabet Med. 2010;27(3):266-273. doi:10.1111/j.1464-5491.2010.02949.x

Supplementary Material

Please find the following supplemental material available below.

For Open Access articles published under a Creative Commons License, all supplemental material carries the same license as the article it is associated with.

For non-Open Access articles published, all supplemental material carries a non-exclusive license, and permission requests for re-use of supplemental material or any part of supplemental material shall be sent directly to the copyright owner as specified in the copyright notice associated with the article.

2.46 MB

0.00 MB