Validity of triple and dualtracer techniques to estimate glucose appearance.  
Jump to Full Text  
MedLine Citation:

PMID: 22454288 Owner: NLM Status: MEDLINE 
Abstract/OtherAbstract:

The tripletracer (TT) dilution technique has been proposed to be the gold standard method to measure postprandial glucose appearance. However, validation against an independent standard has been missing. We addressed this issue and also validated the simpler dualtracer (DT) technique. Sixteen young subjects with type 1 diabetes (age 19.5 ± 3.8 yr, BMI 23.4 ± 1.5 kg/m(2), HbA(1c) 8.7 ± 1.7%, diabetes duration 9.0 ± 6.9 yr, total daily insulin 0.9 ± 0.2 U·kg(1)·day(1), mean ± SD) received a variable intravenous 20% dextrose infusion enriched with [U(13)C]glucose over 8 h to achieve postprandialresembling glucose excursions while intravenous insulin was administered to achieve postprandialresembling levels of plasma insulin. Primed [6,6(2)H(2)]glucose was infused in a manner that mimicked the expected endogenous glucose production and [U(13)C; 1,2,3,4,5,6,6(2)H(7)]glucose was infused in a manner that mimicked the expected glucose appearance from a standard meal. Plasma glucose enrichment was measured by gas chromatographymass spectrometry. The intravenous dextrose infusion served as an independent standard and was reconstructed using the TT and DT techniques with the twocompartment Radziuk/Mari model and an advanced stochastic computational method. The difference between the infused and reconstructed dextrose profile was similar for the two methods (root mean square error 6.6 ± 1.9 vs. 8.0 ± 3.5 μmol·kg(1)·min(1), TT vs. DT, P = NS, paired ttest). The TT technique was more accurate in recovering the overall dextrose infusion (100 ± 9 and 92 ± 12%; P = 0.02). The root mean square error associated with the mean dextrose infusion profile was 2.5 and 3.3 μmol·kg(1)·min(1) for the TT and DT techniques, respectively. We conclude that the TT and DT techniques combined with the advanced computational method can measure accurately exogenous glucose appearance. The TT technique tends to outperform slightly the DT technique, but the latter benefits from reduced experimental and computational complexity. 
Authors:

A Haidar; D Elleri; J M Allen; J Harris; K Kumareswaran; M Nodale; C L Acerini; M E Wilinska; N Jackson; A M Umpleby; M L Evans; D B Dunger; R Hovorka 
Related Documents
:

8321068  Diabetes and diet. we are still learning. 16469808  Dietgenotype interactions in the development of the obese, insulinresistant phenotype... 20225188  Pancreatic fat is negatively associated with insulin secretion in individuals with impa... 8181258  Are the nutritional recommendations for insulindependent diabetic patients being achie... 10202168  Clinical features of 52 neonates with hyperinsulinism. 12365818  Insulin analogues. 
Publication Detail:

Type: Journal Article; Research Support, N.I.H., Extramural; Research Support, NonU.S. Gov't Date: 20120327 
Journal Detail:

Title: American journal of physiology. Endocrinology and metabolism Volume: 302 ISSN: 15221555 ISO Abbreviation: Am. J. Physiol. Endocrinol. Metab. Publication Date: 2012 Jun 
Date Detail:

Created Date: 20120618 Completed Date: 20120827 Revised Date: 20130626 
Medline Journal Info:

Nlm Unique ID: 100901226 Medline TA: Am J Physiol Endocrinol Metab Country: United States 
Other Details:

Languages: eng Pagination: E1493501 Citation Subset: IM 
Affiliation:

Metabolic Research Laboratories, Institute of Metabolic Science, University of Cambridge, United Kingdom. 
Export Citation:

APA/MLA Format Download EndNote Download BibTex 
MeSH Terms  
Descriptor/Qualifier:

Adolescent Algorithms Area Under Curve Blood Glucose / metabolism Carbon Radioisotopes / chemistry, diagnostic use Data Interpretation, Statistical Deuterium / chemistry, diagnostic use Female Gas ChromatographyMass Spectrometry Glucose / metabolism*, pharmacology Hemoglobin A, Glycosylated / analysis Humans Infusions, Intravenous Insulin / blood Intestinal Absorption Isotope Labeling LeastSquares Analysis Male Radioactive Tracers* Radioisotope Dilution Technique* Reproducibility of Results Stochastic Processes Young Adult 
Grant Support  
ID/Acronym/Agency:

1R01 DK085621/DK/NIDDK NIH HHS 
Chemical  
Reg. No./Substance:

0/Blood Glucose; 0/Carbon Radioisotopes; 0/Hemoglobin A, Glycosylated; 0/Insulin; 0/Radioactive Tracers; 50997/Glucose; 7782390/Deuterium 
Comments/Corrections 
Full Text  
Journal Information Journal ID (nlmta): Am J Physiol Endocrinol Metab Journal ID (isoabbrev): Am. J. Physiol. Endocrinol. Metab Journal ID (hwp): ajpendo Journal ID (pmc): ajpendo Journal ID (publisherid): AJPENDO ISSN: 01931849 ISSN: 15221555 Publisher: American Physiological Society, Bethesda, MD 
Article Information Copyright © 2012 the American Physiological Society openaccess: Received Day: 23 Month: 11 Year: 2011 Accepted Day: 21 Month: 3 Year: 2012 Print publication date: Day: 15 Month: 6 Year: 2012 Electronic publication date: Day: 27 Month: 3 Year: 2012 pmcrelease publication date: Day: 15 Month: 6 Year: 2013 Volume: 302 Issue: 12 First Page: E1493 Last Page: E1501 PubMed Id: 22454288 ID: 3378162 Publisher Id: E005812011 DOI: 10.1152/ajpendo.00581.2011 
Validity of triple and dualtracer techniques to estimate glucose appearance  
A. Haidar^{1}^{2}  
D. Elleri^{1}^{3}  
J. M. Allen^{1}^{3}  
J. Harris^{1}  
K. Kumareswaran^{1}  
M. Nodale^{1}  
C. L. Acerini^{3}  
M. E. Wilinska^{1}  
N. Jackson^{4}  
A. M. Umpleby^{4}  
M. L. Evans^{1}  
D. B. Dunger^{1}^{3}  
R. Hovorka^{1}^{3}  
^{1}Metabolic Research Laboratories, Institute of Metabolic Science, University of Cambridge, Cambridge, United Kingdom; 

^{2}Centre for Intelligent Machines, McGill University, Montreal, Quebec, Canada; 

^{3}Department of Paediatrics, University of Cambridge, Cambridge, United Kingdom; and 

^{4}Postgraduate Medical School, University of Surrey, Guilford, United Kingdom 

Correspondence: Address for reprint requests and other correspondence: R. Hovorka, Univ. of Cambridge Metabolic Research Laboratories, and NIHR Cambridge Biomedical Research Centre, Inst. of Metabolic Science, Box 289, Addenbrooke's Hospital, Hills Road, Cambridge CB2 0QQ, UK (email: rh347@cam.ac.uk). 
ASSESSMENTS of postprandial glucose turnover facilitate understanding of carbohydrate physiology in health and diseases (^{7}, ^{22}, ^{30}). In the postprandial state, endogenous glucose production (EGP) and glucose disposal (R_{d}) vary with time together with the appearance of glucose originating from the meal (R_{a meal}). Direct measurements of postprandial turnover are generally not possible, and indirect assessments using the tracer dilution methodology are hindered by the nonsteady state of the glucoregulatory system (^{1}) and illconditioning (small measurement errors propagate into large estimation errors) of the computational problem (^{12}).
The dualtracer (DT) technique is the original method to assess postprandial glucose fluxes (^{29}) and has been adopted in a range of studies (^{19}, ^{21}, ^{22}, ^{30}). The DT technique uses one glucose tracer infused intravenously to measure EGP while a second glucose tracer mixed with the meal allows oral glucose appearance to be determined (^{35}). Infusing the intravenous tracer in a manner that mimics the expected EGP pattern minimizes the change over time in the specific activity or the tracertotracee ratio (TTR) and improves accuracy of EGP estimates (^{3}, ^{18}).
The tripletracer (TT) technique utilizes a third tracer, which is infused in a manner that mimics the expected glucose appearance from the meal (^{1}). The TT technique has been proposed to be the gold standard to measure postprandial glucose appearance (^{1}). The TT technique provides a modelindependent estimate of R_{a meal} if TTR of the meal tracer over the mealmimicking tracer does not change over time. However, for practical reasons, variations in the TTR are inevitable, and errors of unknown magnitude arise resulting from model misspecification and ill conditioning.
Validation of the TT and DT methods against an independent standard is part of our research to measure gut absorption rates of meals containing complex carbohydrates in subjects with type 1 diabetes (T1D) (^{6}). In the present study, we assessed the accuracy of the TT and DT techniques by contrasting the administered rate of intravenous dextrose against the reconstructed rate. The delivery pattern of intravenous dextrose mimicked a mealderived glucose appearance serving as an independent standard. We used an advanced computational approach based on population stochastic modeling to reduce the effect of ill conditioning.
Sixteen young subjects with T1D (age 19.5 ± 3.8 yr, BMI 23.4 ± 1.5 kg/m^{2}, HbA_{1c} 8.7 ± 1.7%, diabetes duration 9.0 ± 6.9 yr, total daily insulin 0.9 ± 0.2 U·kg^{−1}·day^{−1}, mean ± SD) participated in the study performed at the Wellcome Trust Clinical Research Facility, Addenbrooke's Hospital, Cambridge, UK. The study was approved by the local Ethics Committee and subjects provided written informed consent.
Figure 1 shows the study design. Subjects arrived at the clinical facility after breakfast. Long intermediateacting insulin was withdrawn 12 h before admission and replaced with shortacting insulin. From 1000 to 1800 during the stabilization period, subjects fasted and received variable intravenous insulin (Actrapid; Novo Nordisk, Bagsvaerd, Denmark) titrated every 15 min to achieve plasma glucose at 6.0 mmol/l.
A primed constant continuous intravenous infusion of [6,6^{2}H_{2}]glucose (6 mg/kg prime, 0.06 mg·kg^{−1}·min^{−1} continuous infusion; Cambridge Isotope Laboratories, Andover, MA) was given for 150 min from 1530. From 1800, the infusion rate was altered in a manner that mimicked the expected postprandiallike suppression of the EGP [6,6^{2}H_{2}]glucose was infused at a rate of 70, 60, 50, 35, 45, 50, 55, 65, and 70% (100% represents the rate of 0.06 mg·kg^{−1}·min^{−1}) at 0, 10, 20, 30, 220, 240, 280, 300, and 400 min after 1800.
From 1800 until 0200, a variable intravenous insulin infusion was given to achieve plasma insulin concentrations comparable to that resulting from a subcutaneous bolus of rapidacting insulin analog. The total amount of variable insulin corresponded to an insulin bolus associated with a meal containing 120 g of carbohydrates. The calculations utilized individualized insulintocarbohydrate ratios. Starting at 1800, a variable intravenous 20% dextrose infusion (the independent standard mimicking appearance of exogenous glucose from a standard meal) enriched with [U^{13}C]glucose (to represent the tracer that would normally be added to the meal; 13 mg/g, Cambridge Isotope Laboratories) was infused over 8 h to replicate postmeal plasma glucose excursions. Postprandial glucose patterns were taken from a previous study in young subjects with T1D who consumed a standard meal containing 120 g of carbohydrates (^{6}). An adaptive modelpredictive controller (gMPC, version 1.0.2; University of Cambridge, Cambridge, UK) was used to adjust the dextrose infusion utilizing information about subject's total daily insulin dose, the predetermined intravenous insulin delivery, and the desirable plasma glucose profile.
Starting at 1800, [U^{13}C; 1,2,3,4,5,6,6^{2}H_{7}]glucose (mealmimicking tracer; Cambridge Isotope Laboratories) was infused in a manner that mimicked the expected glucose appearance from a standard meal; [U^{13}C; 1,2,3,4,5,6,6^{2}H_{7}]glucose was infused at 20, 40, 55, 70, 90, 100, 90, 70, 50, 40, 30, 20, 10, 5, and 0% of 0.1 mg·kg^{−1}·min^{−1} at 0, 10, 20, 30, 40, 50, 90, 120, 150, 180, 195, 225, 270, 330, and 405 min after 1800.
From 1000, venous blood samples were taken every 10–15 min for the determination of plasma glucose. Samples were taken at 1500, 1515, and 1530 to determine background glucose enrichments. From 1700, venous blood samples were taken every 10–30 min for the determination of plasma insulin, [U^{13}C]glucose, [6,6^{2}H_{2}]glucose, and [U^{13}C; 1,2,3,4,5,6,6^{2}H_{7}]glucose. The samples were immediately centrifuged and separated. Plasma glucose was measured immediately while other samples were stored at −80°C until assayed.
Plasma glucose was measured using a YSI2300 STAT Plus Analyser (YSI, Lynchford House, Farnborough, UK). Plasma insulin was measured using an immunochemiluminometric assay (Invitron, Monmouth, UK) with an intraassay coefficient of variation of 4.7% and an interassay coefficient of variation of 7.2–8.1%. Plasma samples were derivatized to form the methoximetrimethylsilyl derivative for gas chromatography massspectrometry analysis (Agilent 5975C inert XL EI/CI MSD; Agilent Technologies, Wokingham, UK) to measure ions m/z 319.2 (M + 0), 321.2 (M + 2), 322.2 (M + 3), 323.2 (M + 4), 324.2 (M + 5), and 327.2 (M + 8) (^{26}).
TTRs were calculated using a variation of a method described previously (^{12}, ^{14}). In brief, ions m/z M+0, M+2, M+4, and M+8 were used to calculate TTRs of [6,6^{2}H_{2}]glucose, [U^{13}C]glucose, and [U^{13}C; 1,2,3,4,5,6,6^{2}H_{7}]glucose corrected for recycled glucose using ions m/z M+0, M+3, and M+5. The calculations accounted for spectra overlap (^{23}). It was assumed that [U^{13}C; 1,2,3,4,5,6,6^{2}H_{7}]glucose is recycled equally into glucose molecules with m/z of M+4 and M+5 and that [U^{13}C]glucose is recycled equally into glucose molecules with m/z of M+2 and M+3 (^{15}, ^{16}). A system of algebraic equations was solved analytically, and the solution was reduced and validated using Mathematica (Wolfram Research, Champaign, IL). The endogenous glucose concentration was calculated using a modelindependent method as described previously (^{1}).
A stochastic modeling approach was used to reconstruct the exogenous dextrose infusions (R_{a exo}) and to calculate EGP and R_{d}. The DT method utilized the plasma concentration of unlabeled glucose, [6,6^{2}H_{2}]glucose, and [U^{13}C]glucose. The TT methods utilized additionally the plasma concentration of [U^{13}C; 1,2,3,4,5,6,6^{2}H_{7}]glucose. A hierarchical Bayes model was created implementing the Radziuk/Mari twocompartment model of glucose kinetics (^{18}, ^{22}). The method used Bayesian inference with a regularizing prior that assumes smooth glucose fluxes (^{12}), with individualized smoothness levels drawn from a population distribution (^{27}). The method has been described previously (^{9}), and additional details are given in the appendixes.
Values are presented as means ± SD or median (interquartile range) as appropriate. The difference between the actual and reconstructed dextrose infusion rates was determined using the root mean square error (RMSE).
RMSE=∑i=1N(Ra exo,i−Ra exo,i¯)2N 
Plasma glucose and plasma insulin concentrations are shown in Fig. 2. The plasma glucose concentration before the start of dextrose infusion averaged 5.7 ± 1.1 mmol/l and peaked at 11.3 ± 3.7 mmol/l after 120 min, and then decreased slowly to around 8.0 mmol/l by the end of the experiment. The concentration of plasma insulin was 96 (57–197) pmol/l before the start of the dextrose infusion and peaked at 879 (726–986) pmol/l after 70 min. The plasma insulin concentration returned to basal level by 360 min.
The TTRs are shown in Fig. 3. Due to the infusion pattern of [6,6^{2}H_{2}]glucose and [U^{13}C; 1,2,3,4,5,6,6^{2}H_{7}]glucose, the changes in the ratios [6,6^{2}H_{2}]glucose over endogenous glucose and [U^{13}C; 1,2,3,4,5,6,6^{2}H_{7}]glucose over [U^{13}C]glucose were reduced, minimizing the modeldependent error when estimating EGP and R_{a exo} by the TT method.
The ratio [6,6^{2}H_{2}]glucose over [U^{13}C]glucose is shown in Fig. 3. The [6,6^{2}H_{2}]glucose was infused in a manner that mimicked the expected EGP, and this led to marked changes in this particular ratio, potentially leading to a modeldependent error in the estimation of R_{a exo} by the DT method.
Figure 4 shows a sample dextrose infusion profile and its reconstructed profiles by the TT and DT techniques.
The difference between individual actual and individual reconstructed dextrose infusion rates assessed by RMSE is shown in Table 1. The RMSEs associated with the TT and DT techniques were comparable (P = NS). The DT technique slightly underestimated the overall glucose appearance compared with the TT method (P = 0.02). The TT method provided 100% recovery.
The mean dextrose infusion and its reconstructed patterns are shown in Fig. 5. The calculations indicate nearly identical mean profiles obtained with the two methods, documenting the ability of the two methods to reliably estimate exogenous glucose appearance. The RMSE associated with the mean dextrose infusion profile is shown in Table 1.
Estimation errors over time are shown in Fig. 5. The DT technique tended to underestimate the early portion of the dextrose infusion rate. This could be ameliorated by higher sampling frequency. The estimation error of the DT method tended to be less than zero, in concordance with the observation that the overall dextrose recovery was slightly underestimated.
Figure 6 documents a nearly identical EGP and R_{d} obtained with the two methods. The plot of differences shown in Fig. 6 suggests that, on an individual basis, EGP estimates by the two methods were nearly identical as documented by the low variability of the differences at each time point. Differences in individual estimates of R_{d} occurred. The TT technique provided slightly higher estimates of R_{d} at 30 and 180 min.
Table 2 shows partial AUCs associated with the EGP and R_{d} profiles. AUCs of EGP were statistically lower with the DT technique, but the difference lacked clinical significance.
As expected, the fractional clearance rate obtained by the two methods with the use of [6,6^{2}H_{2}]glucose were virtually superimposable (Fig. 7). The two methods are fitting the same [6,6^{2}H_{2}]glucose concentration. Minor differences in the fractional glucose clearance resulted from the way the parameter estimation process was set up utilizing simultaneous fit to all glucose species (see Appendix A). We also confirm that the TT and DT methods were equivalent when estimating EGP. The time course of the factional clearance rate calculated with [U^{13}C; 1,2,3,4,5,6,6^{2}H_{7}]glucose by the TT method was slightly higher during the initial study period. This likely reflects a high relative measurement error associated with low concentrations of [U^{13}C; 1,2,3,4,5,6,6^{2}H_{7}]glucose at the start of the study. This overestimation is linked to a higher [U^{13}C]glucose turnover calculated by the TT method contributing to a slightly higher overall R_{d} (Fig. 6).
The present study suggests that the TT and DT techniques provide accurate estimates of exogenous glucose appearance. The mean dextrose infusion was accurately reconstructed. The total amount of delivered dextrose was well recovered by the TT technique, whereas the DT technique underestimated the amount by 8%.
We adopted a unique validation method mimicking postprandial conditions and utilizing an independent standard represented by the intravenous dextrose infusion. A concomitant intravenous insulin infusion achieved typical postprandial insulin excursions in subjects with T1D. We focused on T1D because accurate estimation of gut absorption is required for optimum timing and amount of prandial insulin (^{2}, ^{25}). Optimal prandial insulin dosing is important to minimize postprandial glucose excursions, which are important determinants of overall glucose control in T1D (^{24}). Examining T1D avoids the need to consider endogenous insulin secretion, facilitating a simpler characterization of the relationship between plasma insulin and glucose excursions. In combination with mathematical modeling, estimates of gut absorption can be used to explore, in silico (^{34}), alternative insulin dosing strategies including closedloop systems (^{5}, ^{10}, ^{11}, ^{13}, ^{33}). Given the common underlying computational principles and experimental procedures, it is reasonable to assume applicability of present results to other populations.
On theoretical grounds, the model misspecification error is the only source for differences between the DT and TT techniques. If the Radziuk/Mari model provided an accurate description of glucose kinetics, the two methods would return identical results. In the early portion of the experiment, when dextrose delivery was rapidly increasing, the DT method underestimated the dextrose infusion. The overall amount of dextrose was slightly underestimated by the DT method, in concordance with results obtained by others (^{1}, ^{32}); otherwise, results provided by the two techniques were comparable. The present study thus suggests that the twocompartment Radziuk/Mari model provides a reasonable approximation of glucose kinetics during meallike experimental scenarios. This is supported by the observation that the fractional turnover rates were comparable apart from the initial part of the study. We argue that the DT technique is the preferred approach, given the reduced cost and reduced experimental and analytic complexity, unless highly accurate measurements of the early postprandial period are required.
The DT approach is equally suitable to estimate EGP as is the TT approach. From a physiological viewpoint, the results were comparable between the two techniques despite the statistical significance in partial AUCs. Indirectly, the present study also validates estimates of EGP. Given the comparability of the computational approach to estimate EGP and dextrose infusion, it is reasonable to assume validity for estimates of EGP, given that the dextrose infusion was accurately reconstructed. The main difference is that EGP originated endogenously but the dextrose infusion exogenously, but this should not affect the validity of the calculations.
The infusion pattern of EGPmimicking tracer was designed to minimize changes over time of the TTR G_{EM}/G_{E} but led to considerable fluctuations in the ratio G_{EM}/G_{M} (Fig. 3), possibly accentuating model misspecification errors by the DT method. Despite glucose fluxes calculated by the DT and TT methods being similar, they indicate that the model misspecification error may have been mitigated by the use of the novel computational method. Further studies are warranted to determine whether the amount of exogenous glucose would be less underestimated by the DT method if the infusion of the EGPmimicking tracer was kept constant, as is the case in many (^{21}, ^{30}), but not all (^{20}), DT studies. We suggest that the infusion pattern of EGPmimicking tracer should be designed to match the research question. Studies focusing on measuring EGP may adopt a pattern similar to that adopted in the present study.
The low initial concentrations of the mealmimicking tracer and the meal tracer at the start of the experiment caused inaccurate estimation of the initial portion of the fractional clearance rate k_{01,MM} (Fig. 7). However, unlike the traditional calculation approach using an analytic solution and the fractional clearance rate (^{1}, ^{18}), our approach avoids the use of fractional clearance rate, shown here for illustrative purposes, and utilizes model fit to measured concentrations. This renders the calculations more robust and less prone to ill conditioning. As the measurement errors associated with the meal tracer and the mealmimicking tracer at the early portion of the experiment are large compared with actual measurements, the initial part of the dextrose infusion estimate will be less dependent on the measured concentrations and will more reflect the assumption of “smooth glucose fluxes”. This is a drawback of the TT study design; the DT approach is not affected by similar considerations.
The novel computational approach may have increased the accuracy of the calculations. The approach assumed that glucose fluxes were smooth to limit the effect of ill conditioning (^{12}, ^{31}). The extent of smoothing was shared among individuals by assuming a population distribution from which individual smoothing estimates were drawn. This “pulled” individual smoothing factors closer to their mean value (^{27}) and avoided aberrations and oscillations, which can be present when data are processed separately. Additionally, this novel approach used an appropriate model of the measurement error. The traditional method (^{1}), although computationally simpler, transforms data, which leads to transformations of the measurement error. This is of particular concern at the initial part of the study when concentrations of meal tracer are low and the associated relative measurement errors are high. This was addressed previously by excluding the initial portion of the data from calculations (^{1}), but all data can be utilized with the present approach.
The computational approach estimates all glucose fluxes simultaneously by fitting measurements of all glucose species in all subjects within a single parameter estimation run (Appendix A). This facilitates interaction between fits to different glucose species. This also explains why the DT and TT methods provide slightly different estimates of endogenous glucose production and fractional clearance rate k01_{EM} while fitting identical EGPmimicking tracer concentrations. The TT method will additionally be fitting the mealmimicking tracer concentration, and this will propagate slightly into estimates of EGP and the fractional clearance rate.
In conclusion, our results suggest that tripletracer and dualtracer techniques combined with advanced computational methods can measure accurately and reliably glucose appearance during postprandiallike conditions. The tripletracer technique tends to outperform slightly the dualtracer technique, but the latter benefits from reduced experimental and analytic complexity.
This work was supported by Juvenile Diabetes Research Foundation (nos. 2220061113, 2220071801, 222009801, 222009802), Diabetes UK (BDA07/0003549), the National Institute of Diabetes and Digestive and Kidney Diseases (1R01 DK085621), Medical Research Council Centre for Obesity and Related metabolic Diseases, and the National Institute for Health Research Cambridge Biomedical Research Centre.
No conflicts of interest, financial or otherwise, are declared by the author(s).
Author contributions: A.H., M.N., N.J., M.U., and R.H. analyzed data; A.H., D.E., and R.H. interpreted results of experiments; A.H. and R.H. prepared figures; A.H. and R.H. drafted manuscript; A.H., D.E., M.W., M.L.E., D.D., and R.H. edited and revised manuscript; A.H., D.E., J.A., J.H., K.K., M.N., C.A., M.W., N.J., M.U., M.L.E., D.D., and R.H. approved final version of manuscript; D.E., C.A., M.L.E., D.D., and R.H. conception and design of research; D.E., J.A., J.H., K.K., M.N., M.W., and R.H. performed experiments.
Angie Watts provided laboratory support. The Diabetes Research Network Laboratory Wales (Dr. Steve Luzio) measured plasma insulin. Tomas Hovorka devised TTR calculations and validated the reduced formulae. We acknowledge support by the staff at the Addenbrooke's Wellcome Trust Clinical Research Facility. We are grateful to study volunteers for their participation.
EGP  Endogenous glucose production 
R_{a exo}  Appearance of unlabeled glucose due to exogenous infusion of dextrose 
R_{d}  Disposal of unlabeled glucose 
R_{dE}  Disposal of unlabeled glucose originating endogenously 
R_{dM}  Disposal of unlabeled glucose originating exogenously 
E  Glucose originating endogenously 
EM  Meal tracer (glucose tracer included with the dextrose) 
MM  Mealmimicking tracer 
T  Total glucose (native and all tracers) 
The computational method utilizes the Bayes theory (^{4}) to calculate the posterior density of unknown fluxes and parameters f(uy,v) defined as:
f(uy,v)=f(y,vu).f(u)/f(y,v) 
The prior and the likelihood functions with the TT method adopting the Radziuk/Mari model are described below. The prior and the likelihood functions for the DT method can be obtained using a similar approach. In what follows, R_{a exo} represents the dextrose infusion, R_{dE} represents the disposal of endogenous glucose, and R_{dM} represents the disposal of glucose originating from the dextrose infusion.
Let n be the number of subjects and m be the number of samples taken in each subject. The prior distribution of the unknown glucose flux f, f = EGP, R_{a exo}, R_{dE}, and R_{dM}, is defined as random walk:
(A1)
fj(i+1)=fj(i)+RWj(i) 
RWj(i)~N(0,εji) 
εji=εjti+1−ti 
exp(εj)~N(ε,τ) 
ε~Γ(0.001,0.001) 
τ~Γ(0.001,0.001) 
Ra exoj(0)=0 
EGPj(0)~N(EGP,τ2) 
τ2~Γ(0.001,0.001) 
EGP~N(10,1−10) 
RdEj(0)=EGPj(0) 
RdMj(0)=0 
QM1(0)=QM2(0)=QMM1(0)=QMM2(0)=0 
QEM1(0)=uEM(0)⋅QE1(0)/RdE(0) 
QEM2(0)=QEM1(0)⋅k21/k12 
QE2(0)=QE1(0)⋅k21/k12 
QE1(0)~N(QE1,τ3) 
QE1~N(600,1−10) 
τ3~Γ(0.001,0.001) 
The time courses of glucose masses of different species as calculated by the Radziuk/Mari model are related to the measurements of tracer and total plasma glucose concentrations G^EM(i), G^M(i), G^MM(i), and G^T(i) as follows:
G^EM(i)=QEM1(i)/V+eEM(i) 
G^M(i)=QM1(i)/V+eM(i) 
G^MM(i)=QMM1(i)/V+eMM(i) 
G^T(i)=QT1(i)/V+eT(i) 
QT1(i)=QEM1(i)+QE1(i)+QMM1(i)+QM1(i)⋅(1+1/TTRdex) 
The method estimates individual parameters such as vertices of random walks RW(i) for each of the four timevarying profiles EGP, R_{a exo}, R_{dE}, and R_{dM}, one “smoothness” parameter ε per random walk (all together four smoothness parameters), EGP at time 0 EGP(0) and the mass of glucose in the accessible compartment at time 0 Q_{E1}(0), and the population parameters EGP,Q_{E1},ε,, τ, τ_{2}, and τ_{3}. In total there are 147 parameters per individual and 6 population parameters. Individual parameters from all subjects and population parameters are estimated simultaneously during a single parameter estimation run (in total, 2,358 parameters) using all available measurements of the total and labeled glucose concentrations.
We use the Markov chain Monte Carlo (MCMC) method (^{8}) to obtain the posterior distributions of the fluxes. The MCMC was implemented using WinBUGS version 1.4 (^{17}, ^{28}), extended by the WBDiff package version 1.9.4 (MRC Biostatistics Unit, Cambridge, UK). The latter is used to solve numerically the differential equations defining the Radziuk/Mari model. WinBUGS is a public domain package facilitating Bayesian estimation using MCMC methods. An outline of MCMC principles is provided below.
For the purposes of parameter estimation, measurement errors were assumed to be normally distributed with zero mean. The measurement errors associated with the total glucose and [6,6^{2}H_{2}]glucose were assumed to be multiplicative with coefficients of variation (CV) of 2 and 5%, respectively. The measurement errors associated with [U^{13}C]glucose and [U^{13}C; 1,2,3,4,5,6,6^{2}H_{7}]glucose were assumed to be multiplicative with 5% CV if tracer concentrations were greater than 0.015 mmol/l and 0.01 mmol/l, respectively, and additive with standard deviations of 0.0075 and 0.005 mmol/l, respectively, otherwise. The additive measurement errors were determined empirically and expressed that, at low tracer concentrations, occurring at the early portion of the experiment after the start of tracer infusions, the measurement instrumentation has precision independent of the measured value and that instrumentation detection limits apply.
The kinetics of glucose species S, S = E, EM, M, and MM are described by the following differential equations.
dQ1,S(t)dt=uS(t)−(k01,S(t)+k21)Q1,S(t)+k12Q2,S(t) 
dQ2,S(t)dt=k21Q1,S(t)−k12Q2,S(t) 
GS(t)=Q1,S(t)V 
The fractional clearance k_{01,S}(t) is defined as:
A2
k01,S(t)={RdM(t)/Q1,MifS=M,MMRdE(t)/Q1,EifS=E,EM 
The total unlabeled glucose disposal R_{d}(t) is equal to summation of R_{dM}(t) and R_{dE}(t). The glucose appearance rate is defined as follows:
uS(t)={Ra exo(t)if S=MEGP(t)ifS=EuEM(t)if S=EMuMM(t)if S=MM 
MCMC methods are a class of algorithms that allow sampling from complex intractable posterior Bayesian distributions. The MCMC methods create a Markov chain that, under fairly general conditions, converges to the posterior distribution (^{8}). The MetropolisHasting algorithm is a generally applicable schema to generate the Markov chain, although faster approaches such as Gibbs sampler exist (^{8}). Basic steps of the MetropolisHasting algorithm applied to present computational problem are given below.
Let the vectorY=[Y(0),...,Y(m−1)],Y(i)=[G^T(i)G^EM(i)G^M(i)G^MM(i)] 
1. Define a Markov chain θ^{(i)}, i = 1,2,..., as θ^{(i)} = [ U^{i}P^{i}].
2. Initialize the Markov chain θ^{(1)} by setting suitable initial values for each of the stochastic variables, e.g., glucose fluxes, smoothing parameters, initial individual EGP values, population means, etc.
3. Propose a new state of the Markov chain θ̄^{(i+1)} that depends only on the previous state θ^{(i)}.
4. Solve numerically the differential equations of the Radziuk/Mari model with accompanying initial conditions and individual glucose fluxes as defined in the proposed state θ̄^{(i+1)}.
5. The proposed state θ̄^{(i+1)} is accepted, i.e., θ^{(i+1)} = θ̄^{(i+1)}, if α, drawn from the standard uniform distribution, satisfies
α<P(θ¯(i+1)Y)f(θ(i+1)θ¯(i+1))P(θ(i+1)Y)f(θ¯(i+1)θ(i+1))=P(Yθ¯(i+1))P(θ¯(i+1))f(θ(i+1)θ¯(i+1))P(Yθ(i+1))P(θ(i+1))f(θ¯(i+1)θ(i+1)) 
Go to Step 3.
After a sufficient number of samples, the Markov chain converges to the joint posterior distribution of the stochastic parameters, and all subsequent samples are considered as posterior realizations of the distribution. Medians of posterior distributions are then used to infer point estimates of unknown parameters. In the present study, the Fisherian approach was utilized to calculate population mean (or median as appropriate) of glucose fluxes across the individuals.
REFERENCES
1..  Basu R,Di Camillo B,Toffolo G,Basu A,Shah P,Vella A,Rizza R,Cobelli C. Use of a novel tripletracer approach to assess postprandial glucose metabolism. Am J Physiol Endocrinol Metab284: E55–E69, Year: 200312485809 
2..  Chase HP,Saib SZ,MacKenzie T,Hansen MM,Garg SK. Postprandial glucose excursions following four methods of bolus insulin administration in subjects with type 1 diabetes. Diabetic Med19: 317–321, Year: 200211943004 
3..  Cobelli C,Mari A,Ferrannini E. Nonsteady state: error analysis of Steele's model and developments for glucose kinetics. Am J Physiol Endocrinol Metab252: E679–E689, Year: 1987 
4..  Congdon P. Bayesian Statistical Modelling. Chichester, UK; Hoboken, NJ: Wiley & Sons, Year: 2006 
5..  ElKhatib FH,Russell SJ,Nathan DM,Sutherlin RG,Damiano ER. A Bihormonal closedloop artificial pancreas for type 1 diabetes. Sci Transl Med2: 27ra27, Year: 2010 
6..  Elleri D,Allen JM,Harris J,Haidar A,Kumareswaran K,Leelarthna L,Nodale M,Wilinska ME,Weston J,Acerini CL,Jackson N,Umpleby AM,Evans ML,Dunger DB,Hovorka R. Glucose appearance of large evening meals with low and high glycaemic loads in type 1 diabetes. Diabetologia54, Suppl 1: S54, Year: 2011 
7..  Ferrannini E,Bjorkman O,Reichard GA Jr,Pilo A,Olsson M,Wahren J,DeFronzo RA. The disposal of an oral glucose load in healthy subjects A quantitative study. Diabetes34: 580–588, Year: 19853891471 
8..  Gilks WR,Richardson S,Spiegelhalter DJ. Markov Chain Monte Carlo in Practice. Boca Raton, FL: Chapman & Hall, Year: 1998 
9..  Haidar A,Potocka E,Boulet B,Umpleby AM,Hovorka R. Estimating postprandial glucose fluxes using hierarchical Bayes modelling. Computer Methods Programs Biomed. Year: 2012Month: 2 Day: 22 [Epub ahead of print] PMID: 22364961. 
10..  Hovorka R. Closedloop insulin delivery: from bench to clinical practice. Nat Rev Endocrinol7: 385–395, Year: 201121343892 
11..  Hovorka R,Allen JM,Elleri D,Chassin LJ,Harris J,Xing D,Kollman C,Hovorka T,Larsen AM,Nodale M,De Palma A,Wilinska ME,Acerini CL,Dunger DB. Manual closedloop insulin delivery in children and adolescents with type 1 diabetes: a phase 2 randomised crossover trial. Lancet375: 743–751, Year: 201020138357 
12..  Hovorka R,Jayatillake H,Rogatsky E,Tomuta V,Hovorka T,Stein DT. Calculating glucose fluxes during meal tolerance test: a new computational approach. Am J Physiol Endocrinol Metab293: E610–E619, Year: 200717519281 
13..  Hovorka R,Kumareswaran K,Harris J,Allen JM,Elleri D,Xing D,Kollman C,Nodale M,Murphy HR,Dunger DB,Amiel SA,Heller SR,Wilinska ME,Evans ML. Overnight closed loop insulin delivery (artificial pancreas) in adults with type 1 diabetes: crossover randomised controlled studies. BMJ342: d1855, Year: 201121493665 
14..  Hovorka R,ShojaeeMoradie F,Carroll PV,Chassin LJ,Gowrie IJ,Jackson NC,Tudor RS,Umpleby AM,Jones RH. Partitioning glucose distribution/transport, disposal, and endogenous production during IVGTT. Am J Physiol Endocrinol Metab282: E992–E1007, Year: 200211934663 
15..  Kalderon B,Korman SH,Gutman A,Lapidot A. Glucose recycling and production in glycogenosis type I and III: stable isotope technique study. Am J Physiol Endocrinol Metab257: E346–E353, Year: 1989 
16..  Katz J,Lee WN,Wals PA,Bergner EA. Studies of glycogen synthesis and the Krebs cycle by mass isotopomer analysis with [U^{13}C]glucose in rats. J Biol Chem264: 12994–13004, Year: 19892753898 
17..  Lunn DJ,Thomas A,Best N,Spiegelhalter D. WinBUGS—a Bayesian modelling framework: concepts, structure, and extensibility. Stat Comput10: 325–337, Year: 2000 
18..  Mari A. Estimation of the rate of appearance in the nonsteady state with a twocompartment model. Am J Physiol Endocrinol Metab263: E400–E415, Year: 1992 
19..  Pennant ME,Bluck LJ,Marcovecchio ML,Salgin B,Hovorka R,Dunger DB. Insulin administration and rate of glucose appearance in people with type 1 diabetes. Diabetes Care31: 2183–2187, Year: 200818650373 
20..  Potocka E,Hovorka R,Baughman RA,Umpleby M,Diaz MLM,Chen R,Boss AH,Richardson PC. AFRESA (TM) suppresses endogenous glucose production earlier than a rapidacting analog (Lispro) and inhaled exubera. Diabetes58, Suppl 1: A61, Year: 2009 
21..  RabasaLhoret R,Burelle Y,Ducros F,Bourque J,Lavoie C,Massicotte D,Peronnet F,Chiasson JL. Use of an alphaglucosidase inhibitor to maintain glucose homoeostasis during postprandial exercise in intensively treated Type 1 diabetic subjects. Diabetes Med18: 739–744, Year: 2001 
22..  Radziuk J,McDonald TJ,Rubenstein D,Dupre J. Initial splanchnic extraction of ingested glucose in normal man. Metabolism27: 657–669, Year: 1978651653 
23..  Rosenblatt J,Chinkes D,Wolfe M,Wolfe RR. Stable isotope tracer analysis by GCMS, including quantification of isotopomer effects. Am J Physiol Endocrinol Metab263: E584–E596, Year: 1992 
24..  Scavone G,Manto A,Pitocco D,Gagliardi L,Caputo S,Mancini L,Zaccardi F,Ghirlanda G. Effect of carbohydrate counting and medical nutritional therapy on glycaemic control in Type 1 diabetic subjects: a pilot study. Diabetic Med27: 477–479, Year: 201020536522 
25..  Schernthaner G,Wein W,Shnawa N,Bates PC,Birkett Preprandial vs MA. postprandial insulin lisproa comparative crossover trial in patients with Type 1 diabetes. Diabetic Med21: 279–284, Year: 200415008840 
26..  ShojaeeMoradie F,Jackson NC,Jones RH,Mallet AI,Hovorka R,Umpleby AM. Quantitative measurement of 3OmethylDglucose by gas chromatography mass spectrometry as a measure of glucose transport in vivo. J Mass Spectrom31: 961–966, Year: 19968831149 
27..  Spiegelhalter D,Thomas A. Graphical modeling for complex stochastic systems: the BUGS project. IEEE Intelligent Syst Appl13: 14–15, Year: 1998 
28..  Spiegelhalter D,Thomas A,Best N,Lunn D. WinBugs User Manual. Cambridge, UK: Medical Research Council Biostatistics Unit, Year: 2003 
29..  Steele R,Bjerknes C,Rathgeb I,Altszuler N. Glucose uptake and production during the oral glucose tolerance test. Diabetes17: 415–421, Year: 19684875169 
30..  Thorburn A,Litchfield A,Fabris S,Proietto J. Abnormal transient rise in hepatic glucose production after oral glucose in noninsulindependent diabetic subjects. Diabetes Res Clin Pract28: 127–135, Year: 19957587920 
31..  Tikhonov AN,Goncharsky A,Stepanov VV,Yagola AG. Numerical Methods for the Solution of IllPosed Problems. Dordrecht; Boston, MA: Kluwer Academic, Year: 1995 
32..  Toffolo G,Basu R,Dalla Man C,Rizza R,Cobelli C. Assessment of postprandial glucose metabolism: conventional dual vs. tripletracer method. Am J Physiol Endocrinol Metab291: E800–E806, Year: 200616720627 
33..  Weinzimer SA,Steil GM,Swan KL,Dziura J,Kurtz N,Tamborlane WV. Fully automated closedloop insulin delivery versus semiautomated hybrid control in pediatric patients with type 1 diabetes using an artificial pancreas. Diabetes Care31: 934–939, Year: 200818252903 
34..  Wilinska ME,Chassin LJ,Acerini CL,Allen JM,Dunger DB,Hovorka R. Simulation environment to evaluate closedloop insulin delivery systems in type 1 diabetes. J Diabet Sci Technol4: 132–144, Year: 2010 
35..  Wolfe R. Radioactive and Stable Isotope Tracers in Biomedicine: Principles and Practice of Kinetic Analysis. New York: WileyLiss, Year: 1992 
Figures
Tables
RSME of individual and mean dextrose profiles and percentage of dextrose recovery as calculated by tripletracer and dualtracer methods
Triple Tracer  Dual Tracer  P Value  

RMSE associated with individual dextrose profiles, μmol•kg^{−1}•min^{−1}  6.6 ± 1.9  8.0 ± 3.5  NS 
RMSE associated with the mean dextrose profile, μmol•kg^{−1}•min^{−1}  2.5  3.3  
Percentage of recovered dextrose infusion^{*}  100 ± 9  92 ± 12  0.02 
Values are means ± SD; n = 16. RMSE, root mean square error; NS, not significant.
*Difference between actual and estimated total dextrose delivery.
AUCs for the first 60 and 480 min of EGP and R_{d} profiles as estimated by dualtracer and tripletracer methods
Triple Tracer  Dual Tracer  P Value  

AUC_{0–480} R_{d}, mmol/kg per 480 min  13.0 ± 2.9  12.6 ± 2.9  NS 
AUC_{0–480} EGP, mmol/kg per 480 min  1.9 ± 0.8  1.8 ± 0.8  <0.001 
AUC_{0–60} R_{d}, mmol/kg per 60 min  1.5 ± 0.5  1.3 ± 0.2  NS 
AUC_{0–60} EGP, mmol/kg per 60 min  0.4 ± 0.2  0.4 ± 0.2  <0.001 
Values are means ± SD; n = 16. AUC area under the curve; EGP, endogenous glucose production.
Article Categories:
Keywords: meal tolerance test, tripletracer technique, dualtracer technique, glucose flux, glucose tracer, Radziuk/Mari model. 
Previous Document: Prolonged infusion of amino acids increases leucine oxidation in fetal sheep.
Next Document: Effects of methyldeficient diets on methionine and homocysteine metabolism in the pregnant rat.