Insulin and C-peptide secretion and kinetics in humans: direct and model-based measurements during OGTT

Andrea Tura1, Bernhard Ludvik2, John J. Nolan3, Giovanni Pacini1, and Karl Thomaseth1

1 Institute of Systems Science and Biomedical Engineering, Italian National Research Council, 35127 Padua, Italy; 2 Division of Endocrinology and Metabolism, Department of Medicine 3, University of Vienna Medical School, A-1090 Vienna, Austria; and 3 Department of Endocrinology, St. James's Hospital, Trinity College, Dublin 8, Ireland


    ABSTRACT
TOP
ABSTRACT
INTRODUCTION
METHODS
RESULTS
DISCUSSION
REFERENCES

To directly evaluate prehepatic secretion of pancreatic hormones during a 3-h oral glucose tolerance test (OGTT), we measured insulin and C-peptide in six healthy control, six obese, and six type 2 diabetic subjects in the femoral artery and hepatic vein by means of the hepatic catheterization technique. Hypersecretion in obesity was confirmed (309 ± 66 nmol in obese vs. 117 ± 22 in control and 79 ± 13 in diabetic subjects, P <=  0.01), whereas early phase secretion was impaired in diabetes. We also measured hepatic insulin extraction (higher in diabetic than in control subjects, P = 0.03) and insulin clearance. The measured data were also used to validate a previously proposed mathematical model, developed to quantify prehepatic secretion, hepatic insulin extraction, and insulin clearance during OGTT, when C-peptide and insulin concentrations are systemically measured. We found good correspondence between experimental data and model estimates for prehepatic insulin secretion (P > 0.3, r2 = 0.93), whereas estimation of hepatic insulin extraction and insulin clearance needs further investigation for improvement.

oral glucose tolerance test; insulin clearance; pancreatic hormones; hepatic catheterization; mathematical modeling


    INTRODUCTION
TOP
ABSTRACT
INTRODUCTION
METHODS
RESULTS
DISCUSSION
REFERENCES

GLUCOSE HOMEOSTASIS in the postprandial state is regulated by the balanced interplay among the absorption of glucose from the gut, splanchnic glucose uptake, and secretion and effectiveness of pancreatic beta -cell hormones that regulate the uptake and production of glucose by target tissues (7). Knowledge of the secretion, kinetics, and clearance of pancreatic hormones and their possible interrelations during a physiological test is quite limited, because the direct assessment of portal levels of various hormones is not feasible in human subjects. On the other hand, measurements of peripheral levels do not always reflect the prehepatic concentration of the hormones and do not elucidate potential interactions of the liver on glucose metabolism.

Thus one of the aims of this study was to directly evaluate endogenous secretion of C-peptide and insulin and their kinetics during an oral glucose tolerance test (OGTT). This is a physiological test involving the normal route of glucose intake. We applied the hepatic catheterization technique, which allows the direct assessment of the transsplanchnic balance of pancreatic hormones.

A second aim was the validation, against the experimental data obtained from the hepatic catheterization technique, of a mathematical model (23) that yields quantitative information on prehepatic beta -cell secretion during the OGTT when C-peptide and insulin concentrations are systemically measured.

Glossary


BCS(t)   C-peptide and insulin secretion rate (pmol/min) from measurements (Eqs. 2-5 and 7-8)
ClearanceC-Pep   Systemic C-peptide clearance (l/min) from measurements (Eq. 7)
ClearanceIns   Systemic insulin clearance (l/min) from measurements (Eq. 8)
CP(t)   C-peptide plasma concentration (pmol/l) predicted by the model (Eq. 9)
CPa(t)   C-peptide concentration (pmol/l) measured in the artery (Eqs. 2 and 7)
CPv(t)   C-peptide concentration (pmol/l) measured in the hepatic vein (Eq. 2)
CPS(t)   C-peptide and insulin secretion rate (pmol · l-1 · min-1) from model estimations (Eqs. 9 and 10)
F   Posthepatic insulin fractional appearance (dimensionless) from model estimations (Eq. 10)
HBF(t)   Hepatic blood flow (l/min) (Eqs. 2-6)
HCL(t)   Hepatic insulin clearance (l/min) from measurements (Eq. 6)
HE(t)   Hepatic insulin degradation (pmol/min) from measurements (Eqs. 3 and 4)
HIFC(t)   Hepatic insulin fractional extraction (dimensionless) from measurements (Eqs. 4-6 and 8)
I(t)   Insulin plasma concentration (pmol/l) predicted by the model (Eq. 10)
Ia(t)   Insulin concentration (pmol/l) measured in the artery (Eqs. 3-5 and 8)
Iv(t)   Insulin concentration (pmol/l) measured in the hepatic vein (Eqs. 3 and 5)
k01   Systemic C-peptide fractional clearance (min-1). It is a fixed parameter in model estimations (Eq. 9)
n   Systemic insulin fractional clearance (min-1) from model estimations (Eq. 10)


    METHODS
TOP
ABSTRACT
INTRODUCTION
METHODS
RESULTS
DISCUSSION
REFERENCES

Subjects

Six male lean, nondiabetic [age 44.7 ± 4.4 yr, body mass index (BMI) 26.5 ± 0.9 kg/m2], six male obese, nondiabetic (44.7 ± 2.4 yr, BMI 35.3 ± 1.2 kg/m2), and six diabetic subjects (5 male, 1 female, 50.5 ± 5.3 yr, BMI 30.8 ± 1.1 kg/m2, diabetes duration 7.7 ± 2.7 yr) participated in the study. All subjects were admitted 3 days before the respective study to the San Diego Veterans Affairs Medical Center's Special Diagnostic and Treatment Unit and consumed a weight maintenance diet containing 55% carbohydrate, 30% fat, and 15% protein. None of the nondiabetic subjects had a positive family history for diabetes or was taking any medication known to affect glucose metabolism. The purpose, nature, and potential risks of the study were explained in detail to all subjects before their written consent was obtained. The study protocol was reviewed and approved by the Human Subjects Committee of the University of California San Diego. All studies were performed at 8.00 AM after a 10- to 12-h overnight fast.

Experimental Procedure

Under local anesthesia with 2% lidocaine, the femoral artery was punctured with an 18-gauge needle, and a 5-French Teflon catheter was introduced and positioned fluoroscopically at the level of the inferior end of the sacroiliac joint. The femoral vein was similarly punctured, and a 6.5-French polyethylene catheter was advanced under fluoroscopic control via the inferior vena cava into the right-sided hepatic vein in an area of adequate blood flow.

Hepatic blood flow was estimated by a primed continuous infusion of indocyanine green (19). The dye infusion was started via an antecubital vein 75 min before glucose ingestion and continued throughout the study. Blood was sampled simultaneously from the artery and the hepatic vein at 10-min intervals starting 45 min after the beginning of green dye infusion. At time 0, the subjects ingested 300 ml of a 75-g glucose solution over 5 min. Arterial and hepatic venous blood was sampled at 15-min intervals to determine the concentrations of glucose, C-peptide, insulin, and indocyanine green for 3 h after glucose ingestion. Hepatic plasma flow was calculated by dividing the green dye infusion rate by arteriohepatic venous dye concentration difference. Hepatic blood flow was estimated by dividing hepatic plasma flow by (1 - hematocrit).

Glucose was measured with a YSI automated glucose analyzer (Yellow Springs Instrument, Yellow Springs, OH) and C-peptide as described in Ref. 5. Insulin was assayed by double-antibody radioimmunoassay (4). Indocyanine green was analyzed by spectrophotometer after precipitation with sodium deoxycholate (8). The measurement errors, expressed as interassay coefficient of variation, were 1.5% for glucose, 5% for insulin, and 10% for C-peptide.

Direct-Measurement Data Analysis

Whole body kinetics were described with a circulatory model that includes the main processes involving the liver. In particular, this organ is represented as a compartment with inputs from the portal vein (pancreatic secretion) and the hepatic artery, and output in the hepatic vein. As a process of substrate disappearance from the liver, degradation in the hepatocytes is considered (Fig. 1). As the overall system can be assumed to be in a quasisteady state, given the slow dynamics of the OGTT, the mass flux of peptide across the liver was described as the steady-state equation
outflow<IT>=</IT>inflow<IT>+</IT>secretion<IT>−</IT>extraction (1)
which can be applied to both C-peptide and insulin.


View larger version (10K):
[in this window]
[in a new window]
 
Fig. 1.   Schematic representation of the physiological system under study. This model is used for calculating pancreatic secretion and extraction by the liver of insulin. The star denotes where the direct measurements are obtained.

For C-peptide, it is known that only a negligible proportion is degraded in the liver (20); therefore, Eq. 1 becomes
CPv(<IT>t</IT>) HBF(<IT>t</IT>)<IT>=</IT>CPa(<IT>t</IT>) HBF(<IT>t</IT>)<IT>+</IT>BCS(<IT>t</IT>) (2)
where CPv(t) and CPa(t) are C-peptide concentrations (pmol/l) in the hepatic vein and in the artery, respectively, HBF(t) is the measured hepatic blood flow (l/min), and BCS(t) is beta -cell C-peptide secretion rate (pmol/min). The only unknown is BCS(t), which can thus be calculated. Given the equimolar release of C-peptide and insulin, BCS(t) also represents beta -cell insulin secretion. The integral between 0 and 180 min of BCS(t) gives the total amount of insulin secretion (nmol). The ratio of total amount of insulin secretion to area under the curve (AUC) of glucose concentration provides an index of beta -cell sensitivity to glucose stimulation.

Equation 1 applied to insulin is
Iv(<IT>t</IT>) HBF(<IT>t</IT>)<IT>=</IT>Ia(<IT>t</IT>) HBF(<IT>t</IT>)<IT>+</IT>BCS(<IT>t</IT>)<IT>−</IT>HE(<IT>t</IT>) (3)
where Iv(t) and Ia(t) are insulin concentrations (pmol/l) in the hepatic vein and in the artery, respectively, and HE(t) is the hepatic insulin degradation (pmol/min) and is expressed as a fraction of the amount of the hormone entering the liver, i.e.
HE(<IT>t</IT>)<IT>=</IT>HIFC(<IT>t</IT>)[Ia(<IT>t</IT>)HBF(<IT>t</IT>)<IT>+</IT>BCS(<IT>t</IT>)] (4)
where HIFC(t) is the hepatic insulin fractional extraction (dimensionless), calculated by substituting Eq. 4 in Eq. 3
HIFC(<IT>t</IT>)<IT>=</IT>1<IT>−</IT>{Iv(<IT>t</IT>) HBF(<IT>t</IT>)<IT>/</IT>[Ia(<IT>t</IT>) HBF(<IT>t</IT>)<IT>+</IT>BCS(<IT>t</IT>)]} (5)
Hepatic insulin clearance HCL(t) (l/min) can be computed as
HCL(<IT>t</IT>)<IT>=</IT>HIFC(<IT>t</IT>) HBF(<IT>t</IT>) (6)
Systemic C-peptide clearance (l/min) was calculated as the ratio of the time integral of secretion rate to that of C-peptide concentration in the artery (which is equal to mixed venous blood concentration)
clearance<SUB>C-Pep</SUB><IT>=</IT><FR><NU><LIM><OP>∫</OP><LL>0</LL><UL>180</UL></LIM> BCS(<IT>t</IT>)d<IT>t</IT></NU><DE><LIM><OP>∫</OP><LL>0</LL><UL>180</UL></LIM> CPa(<IT>t</IT>)d<IT>t</IT></DE></FR> (7)
Similarly, systemic insulin clearance (l/min), which does not include first-pass hepatic degradation, was calculated as the ratio of the time integral of secretion rate multiplied by hepatic insulin fractional delivery, 1 - HIFC(t), to that of arterial insulin concentration
clearance<SUB>Ins</SUB><IT>=</IT><FR><NU><LIM><OP>∫</OP><LL>0</LL><UL>180</UL></LIM> [1<IT>−</IT>HIFC(<IT>t</IT>)]BCS(<IT>t</IT>)d<IT>t</IT></NU><DE><LIM><OP>∫</OP><LL>0</LL><UL>180</UL></LIM> Ia(<IT>t</IT>)d<IT>t</IT></DE></FR> (8)

Model-Based Data Analysis

The aforementioned circulatory model was used for analyzing experimental data. Here, we briefly recall the model of insulin secretion and kinetics from OGTT that we want to validate against the measurements directly obtained with the hepatic catheterization protocol. A detailed description of the model with all the assumptions and hypotheses has been reported in Ref. 23.

For C-peptide, the following mathematical description was adopted
dCP(<IT>t</IT>)<IT>/</IT>d<IT>t=</IT>−<IT>k</IT><SUB>01</SUB>CP(<IT>t</IT>)<IT>+</IT>CPS(<IT>t</IT>) (9)
where CP(t) is the measured plasma C-peptide concentration (pmol/l), k01 is the disappearance constant, which represents the systemic C-peptide fractional clearance (min-1), and CPS(t) is the C-peptide secretion rate estimated by the model (pmol · l-1 · min-1); it also represents insulin secretion, since C-peptide is released equimolarly with insulin. The initial condition is provided by the basal C-peptide level, measured immediately before the glucose load. A considerable amount of insulin is extracted by the liver, and only a fraction F of CPS(t) constitutes the posthepatic appearance of the hormone in the peripheral circulation.

Insulin kinetics are described by
d I(<IT>t</IT>)<IT>/</IT>d<IT>t=</IT>−nI(<IT>t</IT>)<IT>+</IT>F<IT>×</IT>CPS(<IT>t</IT>) (10)
where I(t) is the measured plasma insulin concentration (pmol/l), n is the systemic insulin fractional clearance (min-1), and F × CPS(t) is the posthepatic insulin delivery; (1-F) represents the hepatic insulin fractional extraction (dimensionless). The initial condition is given by the basal insulin value. In this study, CP(t) and I(t) were fitted to CPa(t) and Ia(t) from the hepatic catheterization experiment, respectively. Equations 9 and 10 represent a monocompartmental description for C-peptide and insulin kinetics, which is an approach already adopted in other kinetic models (1, 2, 26).

Individualized time courses of C-peptide and insulin secretion rate CPS(t) were estimated by use of the approach proposed previously (23), which is based on a parametric mathematical representation of CPS(t) in terms of continuous piecewise polynomials (splines). These are made up of a series of quadratic polynomials that are joined at specific knot points (3), i.e., 0, 0, 15, 30, 60, 90, 120, 150, 180, 180, and 180. This sequence is slightly different from that used in Ref. 23 to maintain close resemblance to the sampling schedule of the hepatic catheterization experiments, i.e., 0, 15, 30, 45, 60, 75, 90, 105, 120, 135, 150, 165, and 180 min. Double knots at t = 0 were introduced to enable the description of a rapid increase of insulin secretion at the beginning of the experiment; triple knots at t = 180 min were used to describe the sustained insulin release at the end of the observation interval (23).

Nonlinear weighted least squares were used for estimating the unknown parameters in the spline representation of the C-peptide and insulin secretion rate CPS(t) and the fraction F characterizing insulin appearance. The adopted weights were the inverse of the variance of the measurement errors. Total amount of insulin secretion (nmol/l) was estimated by integrating CPS(t) over 180 min. In the original model (23), insulin fractional clearance n was a fixed parameter assumed equal to that calculated from other studies and not estimated. Here, n was estimated as a model parameter to increase the possibility of using the method in subjects who can exhibit variations at the level of systemic insulin clearance. The C-peptide fractional clearance k01 was still maintained constant (=0.062 min-1), in accord with previous studies (13, 24).

Comparison Between Model Estimates and Direct-Measurement Data

For each set of data, the individualized reconstructed time course CPS(t) was compared with BCS(t) directly measured from the hepatic catheterization experiments. This comparison was also performed for the total amount of insulin secretion in the 180-min interval. It must be noticed that BCS(t) is expressed in picomoles per minute, whereas CPS(t) is expressed per unit volume. Thus, to compare the two variables with the same units, it was necessary to divide BCS(t) by the individual insulin distribution volume, calculated on the assumption of a distribution volume per unit of body weight of 78 ml/kg (6) (6.0 ± 0.2, 7.9 ± 0.3, and 6.9 ± 0.2 l for control, obese, and diabetic subjects, respectively). The measured clearances of insulin and C-peptide, calculated according to Eqs. 7 and 8, were compared with the corresponding model parameters n and k01, respectively. For this purpose, measured clearances were divided by the distribution volume to obtain fractional clearances (min-1). To compare measured and model estimated hepatic insulin fractional extraction, the time average of the measured extraction over the 180-min interval was computed.

Calculations and Statistical Analysis

All numerical calculations and parameter estimations were performed using MATLAB (The Mathworks), and numerical simulations were performed using the PANSYM software (22).

Results are presented as means ± SE unless otherwise designated. Nonparametric tests were used for statistical comparisons. In particular, comparisons between different groups were performed by the Mann-Whitney U-test, whereas those between model and experimental results were done by the Wilcoxon signed-rank test. The relationship between model and experimental results was also investigated by linear regression analysis. Regression was also used to investigate the relationship between measured hepatic blood flow and hepatic insulin fractional extraction and hepatic insulin clearance.


    RESULTS
TOP
ABSTRACT
INTRODUCTION
METHODS
RESULTS
DISCUSSION
REFERENCES

Direct Measurements from Hepatic Catheterization

OGTT. The time courses of the concentrations in the hepatic artery of the measured compounds after the administration of the oral glucose load are presented in Fig. 2. The patterns of C-peptide and insulin were qualitatively similar, with the typical hyperinsulinemia in obese and hypoinsulinemia in diabetic subjects. Glucose levels were not different between control and obese subjects, whereas diabetic subjects showed a marked hyperglycemia. Similar patterns, but higher values, were observed for the concentrations in the hepatic vein (not shown). The basal values and the AUCs for C-peptide, insulin, and glucose are reported in Table 1.


View larger version (20K):
[in this window]
[in a new window]
 
Fig. 2.   Measured concentrations (means ± SE) in the artery of glucose (A), insulin (B), and C-peptide (C) during the oral glucose tolerance test (OGTT). , Control; open circle , obese; black-triangle, diabetic subjects.


                              
View this table:
[in this window]
[in a new window]
 
Table 1.   Measured values and calculated parameters from the hepatic catheterization experiments

Insulin and C-peptide secretion. beta -Cell C-peptide and insulin secretion time course BCS(t) is shown in Fig. 3. The total amount of hormone released by the beta -cell (in 180 min) was higher in obese than in control and diabetic subjects (Table 1). When the total amounts of hormone released in control and diabetic subjects were compared, no difference was found; however, when only the first 60 min were considered, the two secretory patterns were different (41 ± 11 nmol in control and 17 ± 4 in diabetic subjects, P = 0.037), confirming the lack of early-phase release typical of diabetic patients. Moreover, when the total amount of insulin secretion was normalized to the distribution volume of each subject, it was found to be lower in diabetic than in control subjects in the 180-min interval also (Table 2). The beta -cell sensitivity index was 0.08 ± 0.01 nmol/(mmol · l-1 · min-1) in control, 0.19 ± 0.04 in obese and 0.03 ± 0.005 in diabetic subjects, different in each class with respect to the others (P < 0.03).


View larger version (22K):
[in this window]
[in a new window]
 
Fig. 3.   Pancreatic secretion rate of insulin and C-peptide (means ± SE) during OGTT calculated from measured data. , Control; open circle , obese; black-triangle, diabetic subjects.


                              
View this table:
[in this window]
[in a new window]
 
Table 2.   Statistical comparison between model-estimated parameters and the corresponding parameters obtained from the hepatic catheterization experiments

Clearances and hepatic extraction. Time average over the 180-min interval of the hepatic insulin fractional extraction (expressed as percentage of the insulin amount entering the liver) was found to be higher in diabetic than in control but not in obese subjects (Table 1), whereas hepatic insulin clearance was not different in the three groups (Table 1). A significant negative linear relationship (regression coefficient = -4.46 ± 1.05 10-4, P < 0.0001, r2 = 0.37) was found between hepatic blood flow and hepatic insulin fractional extraction in obese subjects (Fig. 4). Conversely, no significant relationship between the two variables was found in control (P = 0.07) and diabetic subjects (P = 0.2). Regression between hepatic blood flow and hepatic insulin clearance was not significant in any class (P = 0.9 in control, P = 0.16 in obese, and P = 0.10 in diabetic subjects). The time average over the 180-min interval of the hepatic blood flow is reported in Table 1.


View larger version (19K):
[in this window]
[in a new window]
 
Fig. 4.   Regression plot of measured time samples of hepatic blood flow vs. those of hepatic insulin fractional extraction for the obese subjects.

Systemic C-peptide and insulin clearances are reported in Table 1. That of C-peptide was higher in obese than in control subjects, but when it was normalized to the distribution volume to obtain the C-peptide fractional clearance, no difference was observed (Table 2). Insulin clearance was not different in the three groups (Table 1). Fractional insulin clearance of obese subjects was lower than in diabetic but not lower than in control subjects (Table 2).

Model-Derived Estimates

Model fit of C-peptide and insulin data. Figure 5 shows the model fit for the mean C-peptide and insulin concentration data along with the pattern of residuals. In Table 2, the model-estimated parameters are shown for the three groups. Their precision, assessed by the coefficient of variation, i.e., the ratio of the diagonal element of the covariance matrix to the parameter value, was 47% for F, 35% for n, and 3% for CPS(t).


View larger version (35K):
[in this window]
[in a new window]
 
Fig. 5.   Insulin (left) and C-peptide (right) concentration values for the control (A), obese (B) and diabetic (C) group directly measured (; means ± SE), and reconstructed by the model (solid curve, representing the average model fit of the data). For every fit, the time course of the residuals is shown (insets; means ± SE).

Comparison between model-derived and measured variables. Model-reconstructed CPS(t) was compared with the measured insulin secretion BCS(t) normalized to the distribution volume, and the two patterns were similar in every group (Fig. 6). A good correlation was found between measured and model-estimated mean values at any sampling time points (Fig. 7). Similarly, the normalized total amount of insulin secretion (Table 2) was not different from the corresponding measured quantity (Table 2), and the two were highly correlated (Fig. 8, A). Moreover, according to experimental results, model estimates confirmed the difference in the normalized total amount of insulin secretion in each class with respect to the others (Table 2).


View larger version (21K):
[in this window]
[in a new window]
 
Fig. 6.   Normalized pancreatic secretion rate directly measured () in the 3 groups (means ± SE) compared with the same variable obtained by the model [mean (solid curves) and confidence interval at 95% (dashed curves)], for control (A), obese (B), and diabetic (C) subjects, respectively.



View larger version (22K):
[in this window]
[in a new window]
 
Fig. 7.   Plot of mean time samples of normalized pancreatic secretion of insulin and C-peptide from measured data vs. model estimates (, control; open circle , obese; black-triangle, diabetic subjects). Unit slope line (dashed line) and regression line at zero intercept (solid line) are reported (regression coefficient = 0.87 ± 0.03, P < 0.0001, r2 = 0.97).



View larger version (21K):
[in this window]
[in a new window]
 
Fig. 8.   Plot of individual values from measured data vs. model estimates for total amount of insulin secretion (A), hepatic insulin extraction (B), and systemic insulin fractional clearance (C) (, control; open circle , obese; black-triangle, diabetic subjects). Unit slope line (dashed line) and regression line at zero intercept (solid line) are reported (regression coefficient = 0.79 ± 0.05, P < 0.0001, r2 = 0.93 for secretion, 1.21 ± 0.09, P < 0.0001, r2 = 0.92 for extraction, 1.06 ± 0.25, P = 0.0005, r2 = 0.52 for clearance, respectively).

Although the correlation among individual values was not excellent (Fig. 8, B), estimated hepatic insulin fractional extraction was not different from that directly measured (Table 2). Despite the fact that individual values were found to be poorly correlated (Fig. 8, C), estimated insulin fractional clearance was not different from the corresponding measured quantity (Table 2). The value assumed for C-peptide fractional clearance was not different from the corresponding measured values (Table 2).


    DISCUSSION
TOP
ABSTRACT
INTRODUCTION
METHODS
RESULTS
DISCUSSION
REFERENCES

The present study addresses two main issues: the analysis of experimental data from hepatic catheterization and the validation of a mathematical model of insulin secretion and kinetics.

By use of hepatic catheterization, direct measurements were obtained for secretion and kinetics of insulin and C-peptide during an OGTT. This test was chosen because it is simple to perform, its use is widespread, and it is considered a "physiological" dynamic test. This study confirms insulin hypersecretion in obese compared with lean subjects and impairment of insulin secretion in type 2 diabetic patients. These patients also showed increased arterial glucose levels, as expected, mirrored by their decreased beta -cell sensitivity to glucose.

In many studies, it is assumed that the level of insulinemia depends almost entirely on insulin secretion; however, the role of metabolic clearance should also be taken into account. Hepatic catheterization allowed the direct assessment of hepatic extraction and clearance as well as systemic insulin and C-peptide clearances. Because no difference in absolute or fractional insulin clearance was observed in obese and diabetic compared with control subjects, the hyperinsulinemia typical of obesity and the hypoinsulinemia of diabetes seem to be due only to enhanced or reduced secretion, respectively. A minor role, however, is played in diabetic subjects by hepatic extraction, which was found to be slightly higher in this group, contributing, therefore, to further lowering of systemic insulin. C-peptide clearance was higher in obese subjects, probably due to the larger distribution volume in these subjects. In fact, when C-peptide clearance was normalized to the distribution volume, no difference was observed between the three groups.

The second aim of this study was the validation of a mathematical model (23) developed to assess prehepatic C-peptide and insulin secretion and hepatic insulin extraction from systemic measurements during an OGTT. Mathematical modeling plays a key role in the estimation of these metabolic parameters that can hardly be measured directly in human subjects, and certainly not in the clinical routine. For that validation, we compared model estimates of insulin secretion time courses and parameters with the corresponding variables directly measured with the hepatic catheterization experiment.

As shown in Fig. 6, the average time courses of prehepatic secretion that the model reconstructs were similar to those measured directly, despite measured patterns that were quite variable. From a quantitative point of view, the values of the total amount of released insulin were very similar. To compare the prehepatic secretion from the experiments with the prehepatic secretion estimated by the model, it was necessary to normalize the former to the insulin distribution volume to have both variables in the same units. The same procedure was used to compare clearances. The choice of the insulin distribution volume is therefore a critical aspect, but there is still a lack of knowledge about this physiological variable. Because we were unable to measure the distribution volume from our experiments, we had to choose a value reported in the literature (6), obtained from a noncompartmental analysis of pork insulin bolus during a hypoglycemic clamp. The same distribution volume per unit of body weight was used for control, obese, and diabetic subjects. In fact, to the best of our knowledge, only a few studies have been carried out on possible differences in the insulin distribution volume between different classes of subjects. McGuire et al. (14) found a reduced distribution volume in obese with respect to nonobese subjects (~30%), and in diabetic with respect to nondiabetic subjects (10-20%). However, the absolute value of the estimate varied by as much as a factor of five, depending on the insulin kinetics model used. Opposite results with regard to diabetic compared with nondiabetic subjects (increase of ~40%) were found by Navalesi et al. (15). Thus, because of the evident uncertainty in the estimation of this parameter, we neglected possible variations that could occur between subjects with different characteristics and pathologies. It is worth noting that we used the distribution volume only to allow comparison between results expressed in different units. If different values were selected, the comparisons would be affected in terms of absolute values, but not in terms of correlation between experimental and model based findings.

Hepatic insulin extraction in each class of subjects as estimated by the model was in accord with the measured mean value, although individual correlation was not excellent. It is worth noting, however, that hepatic insulin extraction from experimental data is likely to be less accurate than the measurement of insulin secretion, due to error propagation in its calculation (see Eqs. 2-5). This is reflected in uncertainty in the comparison between model-based and experimental findings. Another possible reason for the low correlation may be the small number of subjects. However, because the experimental protocol is complex and invasive in humans, it has been possible to carry out this investigation only in a few subjects per class.

One of the assumptions of the model is a constant hepatic extraction for each subject. However, experimental data revealed a significant relationship between hepatic blood flow and hepatic extraction in obese subjects, in accord with conclusions from a previous study (24). Thus these results support the hypothesis that hepatic extraction is not constant in these subjects and that it depends on hepatic blood flow, which is not constant during an OGTT. Consequently, at least for obese subjects, the model estimation of a constant parameter has to be considered an approximation. On the other hand, it is not possible to further increase the number of parameters estimated by the model and maintain an acceptable accuracy. Moreover, both hepatic blood flow and hepatic extraction from experimental data do not actually exhibit marked changes with respect to their time average (standard error with respect to time average equal to 1.8, 2.4, and 2.8% in control, obese, and diabetic subjects, respectively, for hepatic blood flow and 7.5, 4.1, and 4.0% for hepatic extraction). Thus the assumption of a constant model parameter for representing hepatic extraction can be considered a reasonable choice, although it may partially contribute to the observed differences between model estimates and experimental results.

The model presented in this study also allows estimation of systemic insulin clearance; that in the original version (23) was fixed to values obtained from other studies (10-12) or from the literature (18, 21). The data set of the present study, in fact, allowed estimation of at least one more parameter, maintaining an acceptable accuracy in the estimates. We chose to estimate insulin clearance, because it has been reported to exhibit changes in different conditions and to be quite variable, even in normal subjects (6). On the contrary, C-peptide clearance was left at a fixed value, as in the original version, because several studies (16-18, 25) have reported a virtually unchanged C-peptide clearance in many different pathophysiological conditions. Those authors' experience with the intravenous glucose test demonstrated that only in subjects with a clear impairment in renal function was C-peptide clearance significantly reduced (9).

Similarly to the hepatic extraction, the model estimated well in each class the mean values of the measured systemic insulin clearance, but individual values were poorly correlated. This could be due to the limits of the description of insulin kinetics of the model used, although in this case also problems of error propagation in the computation of this variable from experimental data (see Eq. 8), as well as the modest number of subjects, could have played a significant role. A larger number of subjects would help to clarify the reason for the observed poor correlation. C-peptide clearances from experimental results confirmed the adequacy of the choice of a fixed parameter, being similar in the three groups and not different from the value assumed in the model.

In conclusion, in the experimental part of this study, we provided direct confirmation of altered insulin secretion in obese and diabetic subjects. Information was obtained on clearance, which is less often investigated. Finally, it was shown that a previously introduced model of insulin secretion and kinetics during an OGTT in humans provides accurate estimates of prehepatic insulin secretion, as well as some information on hepatic insulin extraction and clearance, in a single subject. Thus, because the direct measure of these variables is not feasible in the clinical routine, the presented model-based approach is a useful tool for the assessment of insulin secretion and kinetics during an OGTT in different pathophysiological conditions, although changes in the model formulation are probably necessary to improve the estimation of insulin clearance and hepatic extraction.


    ACKNOWLEDGEMENTS

We are grateful to Dr. A. Mari of the Institute of Systems Science and Biomedical Engineering (LADSEB) for helpful comments and advice.


    FOOTNOTES

This study was supported in part by the "Progetto Strategico Metodi e Modelli Matematici nello Studio dei Fenomeni Biologici" [Decreto Consiglio Nazionale delle Ricerche (CNR) 82185-19/10/98]. Hepatic catheterization experiments were performed in San Diego (Veterans Affairs Hospital, University of California, San Diego) during a visit of J. Nolan and B. Ludvik, with the support of Dr. J. Olefsky. During the early stage of this study, G. Pacini and K. Thomaseth were supported in part by the CNR with a grant "Progetto Bilaterale" assigned to K. Thomaseth for cooperative projects between LADSEB and the Third Medical Clinic of the University of Vienna.

Address for reprint requests and other correspondence: K. Thomaseth, LADSEB-CNR, Corso Stati Uniti, 4, 35127 Padua, Italy (E-mail: Karl.Thomaseth{at}ladseb.pd.cnr.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.

Received 29 November 2000; accepted in final form 14 June 2001.


    REFERENCES
TOP
ABSTRACT
INTRODUCTION
METHODS
RESULTS
DISCUSSION
REFERENCES

1.   Bergman, RN, Phillips LS, and Cobelli C. Physiologic evaluation of factors controlling glucose tolerance in man: measurement of insulin sensitivity and beta -cell glucose sensitivity from the response to intravenous glucose. J Clin Invest 68: 1456-1467, 1981[ISI][Medline].

2.   Cobelli, C, and Pacini G. Insulin secretion and hepatic extraction in humans by minimal modeling of C-peptide and insulin kinetics. Diabetes 37: 223-231, 1988[Abstract].

3.   De Boor, C. A Practical Guide to Splines. New York: Springer Verlag, 1978.

4.   Desbuquois, B, and Aurbach GD. Use of polyethylene glycol to separate free and antibody-bound peptide hormones in radioimmunoassays. J Clin Endocrinol Metab 33: 732-738, 1971[ISI][Medline].

5.   Faber, OK, Binder C, Markussen J, Heding LG, Naithani VK, Kuzuya H, Blix P, Horwitz DL, and Rubenstein AH. Characterization of seven C-peptide antisera. Diabetes 27, Suppl1: 170-177, 1978[ISI][Medline].

6.   Ferrannini, E, and Cobelli C. The kinetics of insulin in man. Diabetes Metab Rev 3: 335-397, 1987[Medline].

7.   Fery, F, and Balasse EO. Glucose metabolism during the starved-to-fed transition in obese patients with NIDDM. Diabetes 43: 1418-1425, 1994[Abstract].

8.   Gasic, S, Kleinbloesem CH, Heinz G, and Waldhäusl W. Contribution of splanchnic and peripheral vascular tissues to the disposal of angiotensin-II and to regional conversion rates of angiotensin-I: a pilot study in humans. J Cardiovasc Pharmacol 17: 615-620, 1991[ISI][Medline].

9.   Kautzky-Willer, A, Pacini G, Barnas U, Ludvik B, Streli C, Graf H, and Prager R. Intravenous calcitriol normalizes insulin sensitivity in uremic patients. Kidney Int 47: 200-206, 1995[ISI][Medline].

10.   Kautzky-Willer, A, Pacini G, Ludvik B, Schernthaner G, and Prager R. beta -Cell hypersecretion and not reduced hepatic insulin extraction is the main cause of hyperinsulinemia in obese nondiabetic subjects. Metabolism 41: 1304-1312, 1992[ISI][Medline].

11.   Kautzky-Willer, A, Pacini G, Weissel M, Capek M, Ludvik B, and Prager R. Elevated hepatic insulin extraction in essential hypertension. Hypertension 21: 646-653, 1993[Abstract].

12.   Kautzky-Willer, A, Thomaseth K, Clodi M, Ludvik B, Waldhäusl W, Prager R, and Pacini G. beta -Cell activity and hepatic insulin extraction following dexamethasone administration in healthy subjects. Metabolism 45: 486-491, 1996[ISI][Medline].

13.   Kautzky-Willer, A, Thomaseth K, Ludvik B, Nowotny P, Rabensteiner D, Waldhäusl W, Pacini G, and Prager R. Elevated islet amyloid pancreatic polypeptide and proinsulin in lean gestational diabetes. Diabetes 46: 607-614, 1997[Abstract].

14.   McGuire, EA, Tobin JD, Berman M, and Andres R. Kinetics of native insulin in diabetic, obese, and aged man. Diabetes 28: 110-120, 1979[Abstract].

15.   Navalesi, R, Pilo A, and Ferrannini E. Kinetic analysis of plasma insulin disappearance in nonketotic diabetic patients and in normal subjects. A tracer study with 125I-insulin. J Clin Invest 61: 197-208, 1978[ISI][Medline].

16.   Pacini, G. Mathematical models of insulin secretion in physiological and clinical investigations. Comput Methods Programs Biomed 41: 269-285, 1994[ISI][Medline].

17.   Pacini, G, and Thomaseth K. beta -Cell secretion from modeling intravenous and oral glucose tests. In: The Minimal Model Approach and Determinants of Glucose Tolerance, edited by Bergman RN, and Lovejoy JC. Baton Rouge: LSU Press, 1997, p. 243-259.

18.   Polonsky, KS, Licinio-Paixao J, Given BD, Pugh W, Rue P, Galloway J, Karrison T, and Frank B. Use of biosynthetic human C-peptide in the measurement of insulin secretion rates in normal volunteers and type I diabetic patients. J Clin Invest 77: 98-105, 1986[ISI][Medline].

19.   Rowell, LB, Blackmon JR, and Bruce RA. Indocyanine green clearance and estimated hepatic blood flow during mild to maximal exercise in upright man. J Clin Invest 43: 1677-1690, 1964[ISI].

20.   Rubenstein, AH, Pottenger LA, Mako M, Getz GS, and Steiner DF. The metabolism of proinsulin and insulin by the liver. J Clin Invest 51: 912-921, 1971[ISI].

21.   Sherwin, RS, Kramer KJ, Tobin JD, Insel PA, Liljenquist JE, Berman M, and Andres R. A model of the kinetics of insulin in man. J Clin Invest 53: 1481-1492, 1974[ISI][Medline].

22.   Thomaseth, K. PANSYM: a tool for mathematical modelling, analysis and control of metabolic and pharmacokinetic systems. Comput Methods Programs Biomed 41: 99-112, 1994.

23.   Thomaseth, K, Kautzky-Willer A, Ludvik B, Prager R, and Pacini G. Integrated mathematical model to assess beta -cell activity during the oral glucose test. Am J Physiol Endocrinol Metab 270: E522-E531, 1996[Abstract/Free Full Text].

24.   Thomaseth, K, Pacini G, Clodi M, Kautzky-Willer A, Nolan JJ, Prager R, Olefsky JM, and Ludvik B. Amylin release during oral glucose tolerance test. Diabet Med 14: S29-S34, 1997[ISI][Medline].

25.   Van Cauter, E, Mestrez F, Sturis 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: 368-377, 1992[Abstract].

26.   Vølund, A, Polonsky KS, and Bergman RN. Calculated pattern of intraportal insulin appearance without independent assessment of C-peptide kinetics. Diabetes 36: 1195-1202, 1987[Abstract].


Am J Physiol Endocrinol Metab 281(5):E966-E974
0193-1849/01 $5.00 Copyright © 2001 the American Physiological Society