Am J Physiol Endocrinol Metab 290: E169-E176, 2006.
First published September 6, 2005; doi:10.1152/ajpendo.00473.2004
0193-1849/06 $8.00
A minimal model of insulin secretion and kinetics to assess hepatic insulin extraction
Gianna Toffolo,1
Marco Campioni,1
Rita Basu,2
Robert A. Rizza,2 and
Claudio Cobelli1
1Department of Information Engineering, University of Padua, Padua, Italy; and 2Department of Internal Medicine, Division of Endocrinology, Diabetes, Metabolism and Nutrition, Mayo Clinic and Foundation, Rochester, Minnesota
Submitted 5 October 2004
; accepted in final form 30 August 2005
 |
ABSTRACT
|
|---|
The liver is the principal site of insulin degradation, and assessing its ability to extract insulin is important to understand several pathological states. Noninvasive quantification of hepatic extraction (HE) in an individual requires comparing the profiles of insulin secretion (ISR) and posthepatic insulin delivery rate (IDR). To do this, we propose here the combined use of the classical C-peptide minimal model with a new minimal model of insulin delivery and kinetics. The models were identified on insulin-modified intravenous glucose tolerance test (IM-IVGTT) data of 20 healthy subjects. C-peptide kinetics were fixed to standard population values, whereas insulin kinetics were assessed in each individual, along with IDR parameters, thanks to the presence of insulin decay data observed after exogenous insulin administration. From the two models, profiles of ISR and IDR were predicted, and ISR and IDR indexes of
-cell responsivity to glucose in the basal state, as well as during first- and second-phase secretion, were estimated. HE profile, obtained by comparing ISR and IDR profiles, showed a rapid suppression immediately after the glucose administration. HE indexes, obtained by comparing ISR and IDR indexes, indicated that the liver is able to extract 70 ± 9% of insulin passing through it in the basal state and 54 ± 14% during IM-IVGTT. In conclusion, insulin secretion, kinetics, and hepatic extraction can be reliably assessed during an IM-IVGTT by using insulin and C-peptide minimal models.
hepatic extraction; intravenous glucose tolerance test; parameter estimation
THE LIVER PLAYS A MAJOR ROLE in determining circulating insulin levels because it is responsible for the elimination of the prevalent fraction of insulin. A quantitative assessment of hepatic insulin extraction (HE) in an individual, in the basal state and/or during a glucose perturbation, is an essential step in deriving an overall parametric picture of the glucose regulatory system. Because a direct measurement requires invasive protocols, with catheters placed in the artery and hepatic vein, indirect approaches based on mathematical models are essential to infer HE from plasma measurements (8). Several approaches have been proposed that combine the assessment of insulin secretion (ISR) from C-peptide data and insulin delivery rate (IDR) in plasma after its first pass through the liver from insulin data to provide HE noninvasively in humans during a glucose perturbation, either intravenously or orally (6, 7, 22, 24). Linking ISR and IDR to plasma measurements of C-peptide and insulin requires a description of C-peptide and insulin kinetics. Different compartmental structures were assumed for insulin (from 1 to 3 compartments) and C-peptide kinetics (1 or 2 compartments), as well as different functional descriptions for insulin secretion/delivery during glucose perturbation. Because all these approaches were aimed at simultaneously estimating kinetics and secretion, assessment of secretion could be biased because of undesired compensations that may arise in parameter estimation of kinetics and secretion processes (20). Moreover, HE was often assumed to be constant during the experiment, even if evidence exists that this process varies in time during a glucose/insulin perturbation (6, 15), e.g., due to saturable receptor-mediated mechanisms (12) and/or dependence of HE from blood flow (22).
The purpose here is to assess HE during an intravenous glucose load in humans by using the minimal-modeling approach. This approach has already been applied to C-peptide data, measured during an intravenous glucose tolerance test (IVGTT) (20) to develop the minimal model of C-peptide secretion and kinetics and then extended to an insulin-modified (IM)-IVGTT; i.e., a bolus or a 5-min square wave of insulin is infused 20 min after the glucose bolus (21). Standard population values were used for C-peptide kinetics (23) to avoid compensation between kinetic and secretion parameters, thus reliably measuring insulin secretion and indexes of
-cell responsivity to glucose. Unfortunately, the extension of this approach to assess insulin delivery rate and kinetics from insulin data is not straightforward, mainly because standard values are not available for insulin kinetics. To solve this problem, we took advantage of the unique features of IM-IVGTT: these data are more informative than standard IVGTT because the decay of insulin concentration observed after the exogenous insulin administration facilitates a reliable estimation of insulin kinetics in each individual, a prerequisite to avoiding errors in modeling secretion. The combined use of insulin and C-peptide minimal models was then explored as a tool to assess HE.
Database.
The data originated from a previous study (2), to which we refer for details on experimental protocols and analytical methods. Briefly, 20 normal subjects (16 males and 4 females, mean age 44 ± 5 yr, BMI 26.03 ± 0.68 kg/m2) underwent an IM-IVGTT that consisted of an intravenous injection of glucose (0.3 g/kg body wt) followed by a square wave (from 20 to 25 min) of insulin infusion (0.02 U/kg body wt). Blood was sampled at 0, 2, 4, 6, 8, 10, 15, 20, 22, 25, 26, 28, 31, 35, 45, 60, 75, 90, 120, 180, and 240 min after the glucose injection, and glucose, C-peptide, and insulin concentrations were measured.
 |
METHODS
|
|---|
C-peptide model.
The minimal model of C-peptide secretion and kinetics that we described earlier (21) was used to assess ISR. C-peptide kinetics are described by the linear two-compartment model (Fig. 1) that was originally proposed by Eaton et al. (7)
 | (1) |
where CP1 and CP2 (pmol/l) are C-peptide concentrations (above basal) in the accessible and peripheral compartment respectively, k01, k21, and k12 (min1) are transfer rate parameters, and isr (pmol·l1·min1) is secretion above basal entering the accessible compartment, normalized by the volume of distribution of compartment 1.

View larger version (12K):
[in this window]
[in a new window]
|
Fig. 1. Minimal model of C-peptide secretion and kinetics. CP1 and CP2, C-peptide concentrations above basal in the accessible and peripheral compartments, respectively; k01, k21, and k12, trasfer rate parameters; isr, insulin secretion above basal; Y, provision of new insulin to the -cells controlled by glucose concentration; X0, amount of insulin secreted as an impulse during first-phase secretion in response to rapid increase of glucose concentration after glucose bolus.
|
|
It was assumed that isr was proportional to the amount X (pmol/l) of insulin in the
-cells, resulting from the balance between isr and provision Y (pmol·l1·min1) of new insulin to the
-cells.
 | (2) |
 | (3) |
X0 is responsible for first-phase secretion, whereas the slower second phase derives from provision Y, which is controlled by glucose according to the following equation:
 | (4) |
i.e., in response to an elevated glucose level, Y and thus isr tends with a time constant of 1/
(min), toward a steady-state value that is linearly related via parameter
(min1) to glucose concentration G (mmol/l) above the threshold value h.
From model parameters, the ISR profile (pmol/min) was reconstructed as
 | (5) |
where VC is the distribution volume of the accessible compartment, and C1b is the end-test C-peptide concentration.
Indexes of responsivity to glucose of first- and second-phase and basal secretion were also calculated. First-phase responsivity to glucose,
1 (dimensionless), is defined as the ratio between the incremental amount of insulin secreted during the first phase and the maximum increment of plasma glucose concentration
G (mmol/l).
 | (6) |
Second-phase responsivity to glucose,
2 (min1) equals parameter
, which describes the stimulatory effect of glucose concentration on provision
 | (7) |
Finally, basal responsivity to glucose,
b (min1), is defined as
 | (8) |
where Gb is the end-test glucose concentration.
Insulin model.
To assess IDR from plasma insulin concentration, a linear single-compartment model was used (Fig. 2) :
 | (9) |
where I (pmol/l) is insulin concentration above basal (end-test) value, n (min1) is the rate constant of insulin disappearance, idr (pmol·l1·min1) is insulin delivery rate above basal normalized to the volume of distribution, U (pmol/min) is the exogenous insulin input, different from zero in the 2025-min interval, and VI (l) is the insulin distribution volume.

View larger version (14K):
[in this window]
[in a new window]
|
Fig. 2. Minimal model of insulin delivery rate and kinetics. I, plasma insulin concentration above basal, accessible to measurement; n, rate constant of insulin disappearance; U, exogenous insulin input; idr, posthepatic insulin delivery rate above basal. Yidr and X0idr are its parameters, related to provision of new insulin to -cells and first-phase secretion, respectively.
|
|
For idr, a description similar to that proposed to describe insulin secretion from C-peptide data was adopted
 | (10) |
where secretory parameters and variables are labeled "idr" to indicate that they refer to insulin delivery rate and then account for secretion and hepatic extraction.
From model parameters, the IDR profile (pmol/min) was reconstructed as
 | (11) |
where Ib is the end-test insulin concentration.
Indexes of responsivity to glucose of first- and second-phase and basal insulin delivery rate were also calculated by using formulas similar to Eqs. 68:
 | (12) |
 | (13) |
 | (14) |
where Gb is the end-test glucose concentration.
Hepatic insulin extraction.
From ISR and IDR profiles, the HE profile was reconstructed as
 | (15) |
From ISR and IDR indexes, two HE indexes, i.e., during IVGTT and in the basal state, were estimated as
 | (16) |
 | (17) |
where T is the time at which glucose, insulin, and C-peptide concentrations reach their end-test values after the IVGTT perturbation.
Model identification.
Identification of the two models was performed simultaneously in each subject by using SAAM II software (1), with the constraint that the threshold, h, assumes the same value in both models. Secretion parameters m,
,
, and X0 were estimated from plasma C-peptide, and glucose data with C-peptide kinetic parameters were fixed to standard population values, following the method proposed by Van Cauter et al. (23), whereas insulin kinetic (n and VI) and delivery rate parameters (mIDR,
IDR,
IDR, X0IDR) were estimated from plasma insulin and glucose data.
Errors in insulin and C-peptide measurements were assumed to be uncorrelated, Gaussian, zero mean with a variance linked to insulin (I), and C-peptide (C) measurements according to the models, Var(I) = 6 + 0.0055 x I2 and Var(C) = 2000 + 0.001 x C2. The coefficients of variations (CVs) of insulin measurements then ranged from
13% at basal level down to 8% (16) and that of C-peptide from 11% down to 4%.
Statistical analysis.
Data are presented as means ± SD. Pearson's correlation coefficient was used to evaluate univariate correlation.
 |
RESULTS
|
|---|
Average glucose, insulin, and C-peptide concentrations are shown in Fig. 3. A steady-state condition was reached at the end of the experiment at levels (4.58 ± 0.24 mmol/l, 18.3 ± 6.3 pmol/l, and 403 ± 114 pmol/l, respectively) somewhat lower than the pretest levels (5 ± 0.33 mmol/l, 25.3 ± 6.5 pmol/l, and 471 ± 121 pmol/l, respectively).

View larger version (10K):
[in this window]
[in a new window]
|
Fig. 3. Plasma glucose (A), insulin (B), and C-peptide (C) concentration during an IM-IVGTT (mean values in 20 normal subjects).
|
|
The ability of C-peptide and insulin minimal models to describe the data is shown in Fig. 4, where average weighted residuals are shown.

View larger version (22K):
[in this window]
[in a new window]
|
Fig. 4. Weighted residuals (means ± SD in 20 normal subjects) of C-peptide (A) and insulin (B) minimal models.
|
|
Estimated secretion parameters m,
,
, X0, and h and indexes
b,
1, and
2 are shown in Table 1. Estimated insulin kinetic parameters n and VI, delivery rate parameters mIDR,
IDR,
IDR, and X0IDR, and indexes
bIDR,
1IDR, and
2IDR in all individuals are summarized in Table 2.
Parameter h was, on average, 4.6 ± 0.3 mmol/l, not significantly different from the end-test value of glucose concentration.
The minimal model-derived average HE profile is shown in Fig. 5. Starting from a pretest value of 64 ± 12%, HE is rapidly suppressed after the glucose bolus and then slowly reaches, within 120 min, a steady-state value HEb of 70 ± 9%. Individual values of HEb and HE indexes are shown in Table 3 and indicate that
50% of secreted insulin is extracted by the liver during an IM-IVGTT.

View larger version (8K):
[in this window]
[in a new window]
|
Fig. 5. Hepatic extraction during an insulin-modified intravenous glucose tolerance test (mean in 20 normal subjects).
|
|
 |
DISCUSSION
|
|---|
The minimal-modeling approach, originally proposed to assess ISR during an IVGTT on the basis of C-peptide and glucose data, has been applied here to quantify IDR from insulin data, thus making it possible to noninvasively assess HE. The adopted model uses a functional description of the control of glucose on IDR similar to that proposed by us (20) for ISR. The model also includes a description of insulin kinetics to link secretion to insulin measurements taken in plasma. To reliably assess insulin kinetics, a prerequisite to avoid undesired compensations between parameters describing secretion and kinetics (20), the insulin minimal model was identified on IM-IVGTT data. This protocol associates the glucose bolus with a short insulin infusion 20 min after the bolus and was originally proposed to improve the identification of the glucose minimal model (9, 25). It also offers definite advantages for the identification of the insulin minimal model, because insulin decay data, measured after the insulin infusion, allowed a reliable estimation of insulin disappearance rate constant and distribution volume and, thus, of insulin clearance rate in each individual. Secretory parameters were estimated as well, thus allowing us to reconstruct the posthepatic secretion profile as well as derive indexes of first- and second-phase posthepatic secretion.
The combined use of insulin and C-peptide minimal models identified in the same subject provided quantitative information on the ability of the liver to extract insulin during any passage through it, because it enabled reconstruction of HE profile and calculation of two indexes that were able to quantify this process both in the basal state and during the glucose perturbation.
Insulin kinetics.
A variety of models have been used in the literature to describe insulin kinetics, with the number of compartments ranging from one to three. Moreover, it has been suggested that insulin clearance may vary during a glucose and/or insulin perturbation. This is invariably reported in the case of sustained perturbations such as, e.g., 3-h hyperglycemic glucose clamp in humans (19), 1-h insulin infusion in dogs (15), and 42-h low-rate glucose infusion in humans (4). In the case of transient perturbations, such as those induced by an IM-IVGTT, it is generally acknowledged that insulin kinetics can be considered linear until insulin levels of 100150 µU/ml (i.e., 600900 pmol/l) are reached (5, 10, 12, 15). In the present study, a linear single-compartment model was assumed to describe insulin kinetics. The linearity assumption appears reasonable because insulin concentration remains within the linearity range for the largest part of the experiment. The monocompartmental assumption renders the whole model (secretion plus kinetics) uniquely a priori identifiable, which allowed us to estimate kinetic parameters with acceptable precision. Moreover, it provided a reasonable fit, particularly of data collected immediately after the exogenous insulin administration, and reliable estimates of kinetic parameters; e.g., insulin plasma clearance rate (PCR = n·VI) was 1.7 ± 0.36 (means ± SD) l/min, very similar to values reported by some authors (4, 17, 19). These findings confirm that the linear single-compartment model provides a reasonable description of insulin kinetics in situations like the present one. They also support the validity of the minimal-model description of IDR, because we showed (20) that a reliable description of kinetic processes avoids errors in modeling the secretory processes.
IDR.
The model adopts a functional description of the control of glucose on IDR similar to that used in the C-peptide minimal model (20). Parameters and indexes estimated from insulin data (Table 2) assumed values that were different from their counterparts estimated from C-peptide data (Table 1), because they also accounted for HE. Only the threshold of glucose concentration that is able to stimulate ISR and IDR is forced to be the same. The estimated value is very close to the constant level reached by glucose at the end of the experiment, when insulin and C-peptide concentrations are almost constant. These concentrations were adopted as the basal concentrations to allow the two models to reach a steady-state condition at the end of the experiment. The basal responsivity to glucose 
was 0.74 ± 0.26·109·min1, representing 15% as an average of its pancreatic counterpart. By multiplying responsivities, which are related to secretion per unit of distribution volume, by insulin and C-peptide distribution volumes 9.17 and 4.5 l, respectively, basal responsivity accounted for 30% of pancreatic responsivity and correlated well with it (R = 0.53).
First- and second-phase responsivity to glucose were 57 ± 28·109·min1 and 2.2 ± 0.9·109·min1, accounting for 39 and 30%, respectively (77 and 59%, respectively, after volume correction), of their pancreatic counterparts.
HE.
HE was assessed by combining the information derived from the two minimal models. The HE profile showed a rapid suppression immediately after the glucose bolus administration, followed by a rapid increase that was attenuated at 20 min, when insulin was administered (Fig. 5). Then, it returned to steady state at the end of the experiment at a value (70 ± 9%) different from the pretest value (64 ± 12%), which was expected because insulin and C-peptide data are different in the two situations. These results indicate that HE varies in time during a transient glucose perturbation. They are consistent with previous observations, obtained, e.g., in dogs by using invasive methods (13) or in humans by using a modeling approach similar to the one adopted here (6). However, in evaluating the model-predicted pattern of HE, it is important to consider that its ingredients, the model-predicted ISR and IDR, are in fact the rate of appearance in plasma of C-peptide molecules resulting from secretion and transit through the liver and of newly secreted insulin molecules resulting from secretion, HE, and transit through the liver, respectively. However, the liver not only degrades a fraction of insulin molecules but also introduces a delay; e.g., for insulin, a transit time equal to 43 s in humans (8) and 2 min in dogs (7) was reported. Should the transit time of C-peptide be different than that of insulin, the HE profile shown in Fig. 5 accounts not only for the time course of HE but also for differences in the way the liver handles insulin and C-peptide molecules. By contrast, these differences do not affect HE indexes: HEb is calculated in steady-state conditions, i.e., in a situation where inward and outward hepatic fluxes are equal, whereas HE (Eq. 10) is defined from IDR and ISR integrals that are not affected by the liver delay. Our results (Table 3) indicate that the liver is able to extract 70 ± 9% of the insulin passing through it in the end-test basal state and 54 ± 14% during IM-IVGTT. A possible explanation is that liver saturation mechanisms are activated at insulin concentration levels observed during IM-IVGTT. In addition, one can also call for an increase in hepatic blood flow, induced by IM-IVGTT, which has been shown to be associated with a reduced HE (22). Our values are in good agreement with previous findings in the steady state, as reviewed by Ferrannini and Cobelli (8), and during an IVGTT perturbation (6).
Regarding HE estimation, our results indicate that HEb and HE can be estimated with very good precision, i.e., with CVs respectively equal to 3 and 7%. CVs were calculated by propagating the precision of all model parameters included in HEb and HE formulas (Eqs. 15 and 16). However, glucose data, as well as basal insulin and C-peptide levels, are assumed to be noise free, and C-peptide distribution volume VC was fixed to a population value and adjusted according to the body surface area, following the method proposed by Van Cauter et al. (23). Because VC plays a crucial role in HEb and HE calculations, we evaluated through a sensitivity analysis that the average percentage sensitivities of HEb and HE to percentage errors in VC are 0.45 ± 0.23 and 1.02 ± 0.87, respectively. The order of magnitude of errors in VC can be evaluated as 20% from a study on seven subjects for whom reference values were available (11). From the sensitivity analysis, they reflected a 10 and 20% increase in HEb and HE precision, thus moving from 3 and 7% (Table 3) to still acceptable values of 13 and 27%, respectively.
HE and insulin clearance.
HE is a component of insulin PCR, accounting for
58% of the total flux (from the data reviewed in Ref. 8). Our data in fact indicate a significant correlation of PCR with both HEb (R = 0.48, P = 0.03) and HE (R = 0.60, P = 0.005).
The reduction predicted by the two minimal models in HE during IM-IVGTT is apparently in contrast with the linearityassumption of insulin kinetics. However, the 23% reduction in HE estimated during an IM-IVGTT with respect to basal (from 70 to 54%; Table 3), although reflecting in a 23% reduction of the hepatic component of insulin clearance (from 0.58 PCR to 0.45 PCR), accounts only for a lower percentage decrease of insulin clearance (to 0.45 PCR + 0.42 PCR = 0.87 PCR) if nonhepatic insulin clearance is assumed to be linear during IM-IVGTT (thus remaining constant at 0.42 PCR). Because a 13% reduction of insulin clearance is hard to detect from our data, we can speculate that insulin clearance is likely to vary in time during an IM-IVGTT due to variations in HE, but these variations are modest and difficult to detect from plasma insulin data. Only a constant value of PCR can be estimated, which represents the overall clearance value during IM-IVGTT. To check the reliability of these assumptions, we have developed an alternative model that links insulin clearance to HE. The model, detailed in the APPENDIX, describes insulin clearance as the sum of two processes: hepatic clearance, varying in time with a pattern proportional to HE profile, and peripheral clearance, assumed to be constant. Numerical identification was successfully achieved only by fixing a priori the ratio between the two processes at basal state. This is a critical issue that is likely to vary among subjects/conditions. When a plausible value was used, e.g., at basal state, hepatic clearance accounted for 60% of the total clearance, the new model provided a fit similar to that of the original model (data not shown). However, values of basal insulin clearance were plausible in only 17 of our 20 subjects (in 3 subjects, implausibly high values of basal clearance of 5.8, 4.2, and 3.4 l/min were estimated, associated with elevated values of secretion indexes). In the remaining 17 subjects, the new model provided an average insulin clearance that was very similar and well correlated to that estimated with the original model (1.7 ± 0.5 vs. 1.6 ± 0.3 l/min, R = 0.74, P = 0.001). Overall, HE indexes were also very similar and well correlated to those derived from the original model (HEb = 69 ± 9 vs. 72 ± 7%, R = 0.71, P = 0.001; HE = 57 ± 12 vs. 58 ± 10%, R = 0.72, P < 0.001), thus indicating that HE indexes of the original model are robust with respect to variations in insulin clearance observed during an IVGTT.
HE and C-peptide/insulin molar ratios.
Finally, let us relate our HE indexes to C-peptide/insulin molar ratios, which are often used as indicators of hepatic insulin extraction. In Polonsky and Rubenstein (14), the use of the ratio between incremental areas under plasma C-peptide and insulin concentration curves as a reflection of HE was discussed, and it was pointed out that the marked difference in the plasma half-lives of insulin and C-peptide complicates the interpretation of changes in this ratio. By integrating Eq. 9, it is possible to link the incremental area under insulin concentration to the total amount of delivered insulin and to insulin clearance rate (nVI)
 | (18) |
where the integral of the exogenous input U(t) is the administered insulin dose.
Similarly, by integrating Eq. 1, one has a relationship between the incremental area under C-peptide concentration, the total amount of secreted hormone, and C-peptide clearance rate (k01VC):
 | (19) |
During an IM-IVGTT, the ratio between incremental C-peptide and insulin areas (Eqs. 18 and 19) represents only a semiquantitative indicator of HE, because it also depends on the total amount of infused insulin and insulin and C-peptide clearance rates. Conversely, the minimal-model index HE (Eq. 16) represents a correct measure of HE in terms of the incremental areas of pancreatic and posthepatic secretion.
In conclusion, insulin secretion, kinetics, and hepatic extraction can be assessed by using the insulin and C-peptide minimal models on IM-IVGTT data. These models, used in conjunction with the minimal model of glucose kinetics (3), thus provide indexes such as insulin sensitivity, glucose effectiveness, first- and second-phase responsivity to glucose of insulin secretion, insulin clearance, and hepatic extraction, which represent a rich parametric picture of the glucose-regulatory system that is usable for studying pathophysiological states such as impaired glucose tolerance, obesity, and diabetes.
 |
GRANTS
|
|---|
This study was partially supported by National Institutes of Health Grants DK-29953 and RR-00585 and by Ministero dell'Università e della Ricerca Scientifica, Italy
 |
APPENDIX
|
|---|
To include in the model a link between insulin clearance and HE, insulin clearance was described as the sum of two processes: hepatic clearance, varying in time with a pattern proportional to hepatic extraction; and peripheral clearance, which was assumed to be constant. Thus the equation describing insulin kinetics is
 | (A1) |
where nH HE(t) (min1) represents hepatic clearance per unit distribution volume, modulated by HE, and nP (min1) represents the contribution of peripheral clearance, which was assumed to be constant during the experiment.
For HE, Eq. 15 was used
 | (A2) |
where ISR is described by Eq. 5, whereas IDR now accounts for both hepatic and peripheral clearance:
 | (A3) |
The new insulin model consists of Eq. A2 coupled with Eq. 10 to describe idr. It was simultaneously identified with the C-peptide model, Eqs. 14, which provides a description for ISR appearing in Eq. A2. To achieve numerical identification, a constraint was introduced on the relative contribution of hepatic and peripheral clearance in the basal state; namely, it was assumed that, in all subjects, hepatic and peripheral clearance account for 60 and 40%, respectively, of the basal clearance:
 | (A4) |
 |
FOOTNOTES
|
|---|
Address for reprint requests and other correspondence: Claudio Cobelli, Dept. of Information Engineering, Via Gradenigo 6a, 35131, Padua, Italy (e-mail: cobelli{at}dei.unipd.it)
The costs of publication of this article were defrayed in part by the payment of page charges. The article must therefore be hereby marked "advertisement" in accordance with 18 U.S.C. Section 1734 solely to indicate this fact.
 |
REFERENCES
|
|---|
- Barrett PH, Bell BM, Cobelli C, Golde H, Schumitzky A, Vicini P, and Foster DM. SAAM II: simulation, analysis and modeling software for tracer and pharmacokinetics studies. Metabolism 47: 484492, 1998.[CrossRef][ISI][Medline]
- Basu R, Breda E, Oberg AL, Powell CC, Dalla Man C, Basu A, Vittone JL, Klee GG, Arora P, Jensen MD, Toffolo G, Cobelli C, and Rizza RA. Mechanisms of the age-associated deterioration in glucose tolerance: contribution of alterations in insulin secretion, action, and clearance. Diabetes 52: 17381748, 2003.[Abstract/Free Full Text]
- Bergman RN, Ider YZ, Bowden CR, and Cobelli C. Quantitative estimation of insulin sensitivity. Am J Physiol Endocrinol Metab Gastrointest Physiol 236: E667E677, 1979.[Abstract/Free Full Text]
- Byrne MM, Sturis J, and Polonsky KS. Insulin secretion and clearance during low-dose graded glucose infusion. Am J Physiol Endocrinol Metab 268: E21E27, 1995.[Abstract/Free Full Text]
- Cobelli C, Mari A, and Ferrannini E. On linearity of insulin kinetics. Am J Physiol Endocrinol Metab 251: E247E250, 1986.[Free Full Text]
- Cobelli C and Pacini G. Insulin secretion and hepatic extraction in humans by minimal modeling of C-peptide and insulin data. Diabetes 37: 223231, 1988.[Abstract]
- Eaton RP, Allen RC, Schade DS, Erickson KM, and Standefer J. Prehepatic insulin production in man: kinetic analysis using peripheral connecting peptide behaviour. J Clin Endocrinol Metab 51: 520528, 1980.[Abstract]
- Ferrannini E and Cobelli C. The kinetics of insulin in man. II. Role of the liver. Diabetes Metab Rev 3: 365397, 1987.[Medline]
- Finegood DT, Hramiak IM, and Dupre J. A modified protocol for estimation of insulin sensitivity with the minimal model of glucose kinetics in patients with insulin dependent diabetes. J Clin Endocrinol Metab 70: 15381549, 1990.[Abstract]
- Hovorka R, Poeriw JK, Smith GD, Sucharda P, and Sonka J. Five-compartment model of insulin kinetics and its use to investigate action of chloroquine in NIDDM. Am J Physiol Endocrinol Metab 265: E162E175, 1993.[Abstract/Free Full Text]
- Magni P, Bellazzi R, Sparacino G, and Cobelli C. Bayesian identification of a population compartmental model of C-peptide kinetics. Ann Biomed Eng 28: 812823, 2000.[Medline]
- Morishima T, Bradshaw C, and Radziuk J. Measurement using tracers of steady-state turnover and metabolic clearance of insulin in dogs. Am J Physiol Endocrinol Metab 248: E203E208, 1985.[Abstract/Free Full Text]
- Morishima T, Pye S, Bradshaw C, and Radziuk J. Posthepatic rate of appearance of insulin: measurement and validation in the nonsteady state. Am J Physiol Endocrinol Metab 263: E772E779, 1992.[Abstract/Free Full Text]
- Polonsky KS and Rubenstein AH. C-peptide as a measure of the secretion and hepatic extraction of insulin. Pitfalls and limitations. Diabetes 33: 486494, 1984.[Abstract]
- Pye S, Wattarai T, Davies G, and Radziuk J. Comparison of the continuously calculated fractional splanchnic extraction of insulin with its fractional disappearance using a new double tracer technique. Metabolism 42: 145153, 1993.[Medline]
- Robbins DC, Andersen L, Bowsher R, Chance R, Dinesen B, Frank B, Gingerich R, Goldstein D, Widemeyer HM, Haffner S, Hales CN, Jarett L, Polonsky K, Porte D, Skyler J, Webb G, and Gallagher K Report of the American Diabetes Association's task force on standardization of the insulin assay. Diabetes 45: 242256, 1996.[Abstract]
- Shapiro ET, Tillil H, Rubenstein AH, and Polonsky KS. Peripheral insulin parallels changes in insulin secretion more closely than C-peptide after bolus intravenous glucose administration. J Clin Endocrinol Metab 67: 10941099, 1988.[Abstract]
- Tillil H, Shapiro ET, Miller MA, Karrison T, Frank BH, Galloway GA, Rubenstein AH, and Polonsky KS. Dose-dependent effects of oral and intravenous glucose on insulin secretion and clearance in normal humans. Am J Physiol Endocrinol Metab 254: E349E357, 1988.[Abstract/Free Full Text]
- Tillil H, Shapiro T, Rubenstein AH, Galloway JA, and Polonsky KS. Reduction of insulin clearance during hyperglycemic clamp. Dose-response study in normal humans. Diabetes 37: 13511357, 1988.[Abstract]
- Toffolo G, De Grandi F, and Cobelli C. Estimation of beta-cell sensitivity from IVGTT C-peptide data. Knowledge of the kinetics avoids errors in modeling the secretion. Diabetes 44: 845854, 1995.[Abstract]
- Toffolo G, Cefalu W, and Cobelli C. Beta-cell function during insulin modified IVGTT successfully assessed by the C-peptide minimal model. Metabolism 48: 11621166, 1999.[Medline]
- Tura A, Ludvik B, Nolan JJ, Pacini G, and Thomaseth K. Insulin and C-peptide secretion and kinetics in humans: direct and model-based measurements during OGTT. Am J Physiol Endocrinol Metab 281: E966E974, 2001.[Abstract/Free Full Text]
- Van Cauter E, Mestrez F, Sturie J, and Polonsky KS. Estimation of insulin secretion rates from C-peptide levels. Comparison of individual and standard kinetic parameters for C-peptide clearance. Diabetes 41: 368377, 1992.[Abstract]
- Volund A, Polonsky KS, and Bergman RN. Calculated pattern of intraportal insulin appearance without independent assessment of C-peptide kinetics. Diabetes 36: 11951202, 1987.[Abstract]
- Welch S, Gebhart SS, Bergman RN, and Phillips LS. Minimal modeling analysis of intravenous glucose tolerance test-derived insulin sensitivity in diabetic subjects. J Clin Endocrinol Metab 71: 15081518, 1990.[Abstract]
This article has been cited by other articles:

|
 |

|
 |
 
R. Basu, C. D. Man, M. Campioni, A. Basu, K. S. Nair, M. D. Jensen, S. Khosla, G. Klee, G. Toffolo, C. Cobelli, et al.
Effect of 2 Years of Testosterone Replacement on Insulin Secretion, Insulin Action, Glucose Effectiveness, Hepatic Insulin Clearance, and Postprandial Glucose Turnover in Elderly Men
Diabetes Care,
August 1, 2007;
30(8):
1972 - 1978.
[Abstract]
[Full Text]
[PDF]
|
 |
|

|
 |

|
 |
 
C. Cobelli, G. M. Toffolo, C. D. Man, M. Campioni, P. Denti, A. Caumo, P. Butler, and R. Rizza
Assessment of beta-cell function in humans, simultaneously with insulin sensitivity and hepatic extraction, from intravenous and oral glucose tests
Am J Physiol Endocrinol Metab,
July 1, 2007;
293(1):
E1 - E15.
[Abstract]
[Full Text]
[PDF]
|
 |
|

|
 |

|
 |
 
G. Bock, C. Dalla Man, M. Campioni, E. Chittilapilly, R. Basu, G. Toffolo, C. Cobelli, and R. Rizza
Effects of Nonglucose Nutrients on Insulin Secretion and Action in People With Pre-Diabetes
Diabetes,
April 1, 2007;
56(4):
1113 - 1119.
[Abstract]
[Full Text]
[PDF]
|
 |
|

|
 |

|
 |
 
R. Basu, C. Dalla Man, M. Campioni, A. Basu, K. S. Nair, M. D. Jensen, S. Khosla, G. Klee, G. Toffolo, C. Cobelli, et al.
Two Years of Treatment With Dehydroepiandrosterone Does Not Improve Insulin Secretion, Insulin Action, or Postprandial Glucose Turnover in Elderly Men or Women
Diabetes,
March 1, 2007;
56(3):
753 - 766.
[Abstract]
[Full Text]
[PDF]
|
 |
|

|
 |

|
 |
 
M. Campioni, G. Toffolo, L. T. Shuster, F. J. Service, R. A. Rizza, and C. Cobelli
Incretin effect potentiates beta-cell responsivity to glucose as well as to its rate of change: OGTT and matched intravenous study
Am J Physiol Endocrinol Metab,
January 1, 2007;
292(1):
E54 - E60.
[Abstract]
[Full Text]
[PDF]
|
 |
|

|
 |

|
 |
 
G. Bock, C. Dalla Man, M. Campioni, E. Chittilapilly, R. Basu, G. Toffolo, C. Cobelli, and R. Rizza
Pathogenesis of Pre-Diabetes: Mechanisms of Fasting and Postprandial Hyperglycemia in People With Impaired Fasting Glucose and/or Impaired Glucose Tolerance
Diabetes,
December 1, 2006;
55(12):
3536 - 3549.
[Abstract]
[Full Text]
[PDF]
|
 |
|

|
 |

|
 |
 
R. Basu, C. Dalla Man, M. Campioni, A. Basu, G. Klee, G. Toffolo, C. Cobelli, and R. A. Rizza
Effects of age and sex on postprandial glucose metabolism: differences in glucose turnover, insulin secretion, insulin action, and hepatic insulin extraction.
Diabetes,
July 1, 2006;
55(7):
2001 - 2014.
[Abstract]
[Full Text]
[PDF]
|
 |
|
Copyright © 2006 by the American Physiological Society.