Effects of aerobic or resistance exercise training on cardiovascular autonomic function of subjects with type 2 diabetes: A pilot study.

BACKGROUND AND AIMS: Both aerobic (AER) and resistance (RES) training improve metabolic control in patients with type 2 diabetes (T2DM). However, information on the effects of these training modalities on cardiovascular autonomic control is limited. Our aim was to compare the effects of AER and RES training on cardiovascular autonomic function in these subjects. METHODS AND RESULTS: Cardiovascular autonomic control was assessed by Power Spectral Analysis (PSA) of Heart Rate Variability (HRV) and baroreceptors function indexes in 30 subjects with T2DM, randomly assigned to aerobic or resistance training for 4 months. In particular, PSA of HRV measured the Low Frequency (LF) and High Frequency (HF) bands of RR variations, expression of prevalent sympathetic and parasympathetic drive, respectively. Furthermore, we measured the correlation occurring between systolic blood pressure and heart rate during a standardized Valsalva maneuver using two indexes, b2 and b4, considered an expression of baroreceptor sensitivity and peripheral vasoactive adaptations during predominant sympathetic and parasympathetic drive, respectively. After training, the LF/HF ratio, which summarizes the sympatho-vagal balance in HRV control, was similarly decreased in the AER and RES groups. After AER, b2 and b4 significantly improved. After RES, changes of b2 were of borderline significance, whereas changes of b4 did not reach statistical significance. However, comparison of changes in baroreceptor sensitivity indexes between groups did not show statistically significant differences. CONCLUSION: Both aerobic and resistance training improve several indices of the autonomic control of the cardiovascular system in patients with T2DM. Although these improvements seem to occur to a similar extent in both training modalities, some differences cannot be ruled out. CLINICAL TRIAL REGISTRATION NUMBER: NCT01182948, clinicaltrials.gov.

Pharmacokinetics and pharmacodynamics of protamine zinc recombinant human insulin in healthy dogs.

Protamine zinc insulins are generally considered to be long acting, with slow absorption from subcutaneous tissue. Protamine zinc recombinant human insulin (PZIR) may be useful to treat diabetic dogs. The purpose of this study was to describe the pharmacokinetics and pharmacodynamics of PZIR in dogs. PZIR was administered subcutaneously to 10 healthy Beagles using an incomplete crossover design, at doses of 0.3 or 0.5 U/kg (each n=5), 0.8 U/kg (n=10), or 0.8 U/kg at three separate sites (n=6). Insulin and glucose concentrations were measured over 24 h. The shapes of insulin and glucose curves were variable among dogs, and the relationship between insulin dose, concentration, and glucose-lowering effect was nonlinear. For single-site 0.8 U/kg, median (range) onset of action was 3.5 h (0.5-10 h), time to glucose nadir was 14 h (5 to >24 h), and duration of action was >24 h (16 to >24 h). Mathematical model predictions of times to 50% and 90% insulin absorption, and fraction of insulin absorbed in 24 h, were not significantly different among protocols. Results confirm the tendency toward a late onset and long duration of action for PZIR in dogs. This insulin may be an alternative treatment option for diabetic dogs.

Fatty acid turnover, substrate oxidation, and heat production in lean and obese cats during the euglycemic hyperinsulinemic clamp.

Simultaneous application of the euglycemic hyperinsulinemic clamp (EHC) and indirect calorimetry was used to examine heat production, fat, and glucose metabolism in lean and obese adult neutered male and female cats. The results show that in lean insulin-sensitive cats glucose oxidation predominated during fasting, whereas lipid oxidation became more prominent in obese cats. Insulin infusion during the EHC in lean cats and obese male cats led to a large increase in glucose oxidation, glycogenesis, and lipogenesis. It also led to an increase in glucose oxidation and glycogenesis in obese female cats but it was significantly less compared to lean cats and obese males. This indicates that obese females show greater metabolic inflexibility. In obese cats of either gender, insulin caused greater suppression of non-esterified fatty acids compared to lean cats suggesting that obese cats show greater fatty acid clearance than lean cats. The heat production per metabolic size was lower in obese cats than lean cats. This would perpetuate obesity unless food intake is decreased. The higher glucose oxidation rate in obese neutered male cats suggests that they are able to replete their glycogen and lipid stores at a faster rate than females in response to insulin and it implies that they gain weight more rapidly. Further studies are needed to investigate if the different response to insulin of male cats is involved in their higher risk to develop diabetes.

Assessment and mathematical modeling of glucose turnover and insulin sensitivity in lean and obese cats.

Insulin sensitivity (SI) of glucose disposal can be quantified with the euglycemic hyperinsulinemic clamp (EHC) with tracer glucose infusion. True steady state is, however, difficult to achieve, and non-steady state analysis of EHC data is preferred. This analysis requires information on glucose kinetics that can be obtained from bolus injection of cold and tracer glucose. The aim of this study was to assess glucose kinetics in cats. Mathematical modeling and non-steady state analysis was applied to assess effects of obesity on glucose turnover, glycolysis/glycogen synthesis, SI, and inhibition of endogenous glucose production (EGP) in lean cats (L) and obese cats (O). D-[3-(3)H]-glucose kinetics and 3H-H2O production were analyzed in 4 L and 4 O with three-compartment modeling. Frequently sampled insulin-modified intravenous glucose tolerance tests (FSIGT) with minimal model analysis were performed in 5L and 3 O to assess glucose kinetics and SI. EHC was performed in 10 L and 10 O with primed-constant infusion of 3H-glucose. Data were analyzed with a modified minimal model segregating suppression of EGP by insulin using a non-linear mixed-effects population approach. FSIGT provided estimates of SI, glucose effectiveness SG, and distribution volume. EHC provided estimates of SI, SG, glycolysis, and suprabasal insulin concentration for 50% EGP inhibition. Obesity appears to affect glucose distribution but not utilization at basal insulin, and reduces SI estimated by FSIGT and EHC. Differences in SI between FSIGT and EHC depend on different descriptions of EGP inhibition by insulin. Finally, glucose disposal at basal insulin appears to occur entirely through glycolysis, whereas significant amounts of glucose are sequestrated from oxidation during EHC.

The Valsalva manoeuvre: physiology and clinical examples.

The Valsalva manoeuvre (VM), a forced expiratory effort against a closed airway, has a wide range of applications in several medical disciplines, including diagnosing heart problems or autonomic nervous system deficiencies. The changes of the intrathoracic and intra-abdominal pressure associated with the manoeuvre result in a complex cardiovascular response with a concomitant action of several regulatory mechanisms. As the main aim of the reflex mechanisms is to control the arterial blood pressure (BP), their action is based primarily on signals from baroreceptors, although they also reflect the activity of pulmonary stretch receptors and, to a lower degree, chemoreceptors, with different mechanisms acting either in synergism or in antagonism depending on the phase of the manoeuvre. A variety of abnormal responses to the VM can be seen in patients with different conditions. Based on the arterial BP and heart rate changes during and after the manoeuvre several dysfunctions can be hence diagnosed or confirmed. The nature of the cardiovascular response to the manoeuvre depends, however, not only on the shape of the cardiovascular system and the autonomic function of the given patient, but also on a number of technical factors related to the execution of the manoeuvre including the duration and level of strain, the body position or breathing pattern. This review of the literature provides a comprehensive analysis of the physiology and pathophysiology of the VM and an overview of its applications. A number of clinical examples of normal and abnormal haemodynamic response to the manoeuvre have been also provided.

Power spectral analysis of heart-rate variations improves assessment of diabetic cardiac autonomic neuropathy.

Power spectral analysis (PSA) of heart-rate variations has recently proved a useful tool in evaluating cardiovascular autonomic activity. It offers the possibility of examining both the functioning of parasympathetic and sympathetic pathways through breakdown into two frequency bands, and of their effects on heart-rate cyclic variability. We applied an autoregressive model for PSA to study overall autonomic tone in 20 male age-matched control subjects and 53 insulin-dependent (type I) diabetic subjects, subdivided into three groups of 20, 15, and 18, each group presenting different degrees of autonomic involvement. We found that: 1) power spectrum density (PSD) values at high-frequency bands (parasympathetic dependent) were similar in diabetic subjects without cardiac autonomic neuropathy (CAN) and in control subjects, but differed significantly from diabetic subjects with mild CAN and severe CAN, both standing and lying; 2) PSD values at low frequency (mainly sympathetic dependent) were similar, or slightly different, in diabetic subjects without CAN and in control subjects, but differed significantly from diabetic subjects with mild and severe CAN, both standing and lying; 3) as an expression of parasympathetic versus sympathetic coherence, correlations, both standing and lying, existed between PSD values at low- and high-frequency bands in control and diabetic subjects without CAN, but not in diabetic subjects with CAN; and 4) different degrees of correlation characterized the PSD values of high and low frequencies versus traditional cardiovascular test values in the diabetic subjects. The best correlation was between PSD low-frequency values and the lying-to-standing maneuver.(ABSTRACT TRUNCATED AT 250 WORDS)

Elevated islet amyloid pancreatic polypeptide and proinsulin in lean gestational diabetes.

Recent research indicates that islet amyloid pancreatic polypeptide (IAPP) might have a regulatory effect on beta-cell insulin processing and secretion. To study such interaction in more detail, IAPP secretion and kinetics and the serum concentrations of proinsulin were assessed both before and after delivery in lean pregnant women with gestational diabetes mellitus (GDM patients) in comparison to those with normal glucose tolerance (NGT) and to nonpregnant healthy lean (control) and obese insulin-resistant women during oral glucose tolerance tests. Kinetic analysis of IAPP was performed with mathematical modeling, providing indirect estimates of its secretion and fractional clearance. Total insulin secretion per 180 min was elevated by 30% in GDM patients (35 +/- 3 pmol/l) versus control subjects (27 +/- 1 pmol/l) (P < 0.05), but increased even more (190-250%) in obese insulin-resistant women, compared with all other groups (68 +/- 7 pmol/l, P < 0.0005). Pregnancy induced a more marked fourfold increase in apparent total IAPP secretion rate (TIR) (GDM patients, 172 +/- 31 pmol x 1(-1) x 3 h(-1); NGT subjects, 166 +/- 31 pmol x 1(-1) x 3 h(-1); control subjects, 40 +/- 1 pmol 1(-1) x 3 h(-1)) and a twofold rise in its fractional clearance versus control subjects (P < 0.01), whereas in GDM patients a 30% increase of IAPP secretion and a decreased clearance was found, compared with obese insulin-resistant women (TIR, 112 +/- 14 pmol x 1(-1) x 3 h(-1)). The increase in IAPP secretion in both pregnant groups was much higher than that of the insulin groups, resulting in a marked change of the IAPP-insulin cosecretion factor when compared with lean or obese nonpregnant women (P < 0.0005). Associated serum proinsulin and the postprandial (total divided by 180 min) proinsulin-to-insulin ratio were greater in GDM patients versus NGT and control subjects (0.11 +/- 0.01 vs. 0.07 +/- 0.01 and 0.08 +/- 0.01 pmol/l, P < 0.05), while the fasting proinsulin-to-insulin ratio was only increased in GDM patients versus control subjects (0.22 +/- 0.03 vs. 0.13 +/- 0.01 pmol/l, P < 0.05). After delivery, total IAPP secretion (52.4 +/- 1.5 pmol/l) was completely normalized in the GDM group, as were the clearance rate and the IAPP-insulin cosecretion factor. Similarly, serum proinsulin concentrations returned to normal, whereas proinsulin-to-insulin ratios remained elevated. In conclusion, IAPP hypersecretion is characteristic for pregnancy and might partially decrease hyperinsulinemia in pregnancy by inhibiting insulin secretion. Increased proinsulin concentrations and a raised proinsulin-to-insulin ratio, which did not abate following delivery, are specific to GDM and might thus serve as its marker and potentially even identify subjects at high risk for the development of NIDDM.

Pronounced insulin resistance and inadequate beta-cell secretion characterize lean gestational diabetes during and after pregnancy.

OBJECTIVE: To evaluate beta-cell secretion and glucose metabolism in lean subjects with gestational diabetes mellitus (GDM) compared with that in subjects with normal pregnancy and obesity. RESEARCH DESIGN AND METHODS: Insulin secretion, insulin sensitivity (S1), and hepatic insulin extraction were assessed in pregnant women with GDM before and after delivery and in those with normal glucose tolerance (NGT) in comparison to healthy nonpregnant lean and obese women. Kinetic analysis of glucose, insulin, and C-peptide plasma concentrations during oral and intravenous glucose tolerance tests was performed by mathematical modeling. RESULTS: S1 was blunted in pregnant women with GDM by 84% and in those with NGT by 66% compared with lean nonpregnant women (P < 0.005 vs. healthy nonpregnant lean control subjects; P < 0.05, GDM vs. pregnant women with NGT), whereas glucose effectiveness was decreased by 33% in both pregnant groups (P < 0.05 vs. healthy nonpregnant lean control subjects). Insulin secretion was 30% higher (P < 0.05) in subjects with GDM than in pregnant women with NGT or in nonpregnant lean women, but decreased (P < 0.005) when compared with obese women with a comparable degree of insulin resistance. Fractional hepatic insulin extraction was similar in both pregnant groups, being lower (P < 0.0001) by 30% versus nonpregnant females. beta-cell sensitivity to glucose for insulin release was decreased in subjects with GDM versus pregnant women with NGT as well as nonpregnant women by 40-50% (P < 0.01). Twelve weeks after delivery, GDM returned to normal glucose tolerance, but S1 remained 50% lower than that in lean nonpregnant women, while beta-cell sensitivity to glucose did not change (P < 0.01 vs. healthy nonpregnant lean control subjects). CONCLUSIONS: Pregnancy is characterized by insulin resistance, diminished hepatic insulin extraction, and glucose effectiveness. Lean subjects with GDM are additionally characterized by having more pronounced insulin resistance and inadequate insulin secretion, which persist after delivery. Compared with other insulin-resistant prediabetic states like impaired glucose tolerance (IGT), defective insulin secretion seems to be a predominant defect in lean GDM subjects, indicating that it might represent a specific prediabetic condition.

Power spectral analysis of heart rate in hyperthyroidism.

The aim of the present study was to evaluate the impact of hyperthyroidism on the cardiovascular system by separately analyzing the sympathetic and parasympathetic influences on heart rate. Heart rate variability was evaluated by autoregressive power spectral analysis. This method allows a reliable quantification of the low frequency (LF) and high frequency (HF) components of the heart rate power spectral density; these are considered to be under mainly sympathetic and pure parasympathetic control, respectively. In 10 newly diagnosed untreated hyperthyroid patients with Graves' disease, we analyzed power spectral density of heart rate cyclic variations at rest, while lying, and while standing. In addition, heart rate variations during deep breathing, lying and standing, and Valsalva's maneuver were analyzed. The results were compared to those obtained from 10 age-, sex-, and body mass index-matched control subjects. In 8 hyperthyroid patients, the same evaluation was repeated after the induction of stable euthyroidism by methimazole. Heart rate power spectral analysis showed a sharp reduction of HF components in hyperthyroid subjects compared to controls [lying, 13.3 +/- 4.1 vs. 32.0 +/- 5.6 normalized units (NU; P < 0.01); standing, 6.0 +/- 2.7 vs. 15.0 +/- 4.0 NU (P < 0.01); mean +/- SEM]. On the other hand components were comparable in the 2 groups (lying, 64.0 +/- 6.9 vs. 62.0 +/- 6.5 NU; standing, 77.0 +/- 6.5 vs. 78.0 +/- 5.4 NU). Hence, the LF/HF ratio, which is considered an index of sympathovagal balance, was increased in hyperthyroid subjects while both lying (11.3 +/- 4.5 vs. 3.5 +/- 1.1; P < 0.05) and standing (54.0 +/- 12.6 vs. 9.8 +/- 2.6; P < 0.02). This parameter was positively correlated with both T3 (r = 0.61; P < 0.05) and free T4 (r = 0.63; P < 0.05) serum levels. Among traditional cardiovascular autonomic tests, the reflex response of heart rate during lying to standing was significantly lower in hyperthyroid patients than in controls (1.12 +/- 0.03 vs. 1.31 +/- 0.04; P < 0.002). No statistically significant difference in reflex responses between the two groups was found in deep breathing or Valsalva's maneuver. In the 8 patients reexamined after methimazole treatment, we observed complete normalization of altered cardiovascular parameters, with slight predominance of the vagal component compared with controls. These results suggest that thyroid hormone excess may determine reduced parasympathetic activity and, thus, a relative hypersympathetic tone.

Estimation of occupational exposure to 2,3,7,8-tetrachlorodibenzo-p-dioxin using a minimal physiologic toxicokinetic model.

In this study we investigated estimation of occupational exposure to 2,3,7,8-tetrachlorodibenzo-p-dioxin (TCDD) based on a minimal physiologic toxicokinetic model in humans. Our purpose was to obtain a mathematical tool for dose-response studies based on human data. We first simplified an existing model of TCDD kinetics in humans and estimated its parameters (i.e., liver elimination and background input of TCDD) using repeated measures of serum dioxin taken in Vietnam veterans (Ranch Hand data and data from an unexposed reference group). We carried out computer simulation and estimation of the model parameters both under a nonlinear weighted least-squares model (naive pooled data approach) and under a nonlinear mixed-effects model. The best parameter estimates were obtained with log-transformed data under a mixed-effects model: liver elimination parameter kf = 0.022 days-1 (95% confidence interval [CI] = 0.020, 0.024), and background input rate input = 0.1251 pg/kg/day (95% CI = 0.071, 0.179). The dioxin kinetic model and its estimated parameters were then used to provide dose estimates for a cohort of workers with exposure to TCDD at chemical plants in the United States. First, the model was used to estimate the rate of occupational intake of TCDD in a subset of the cohort consisting of 253 subjects for whom one measure of serum TCDD was available. A model of change in body-mass index over time was also identified for this subsample. The occupational exposure rate was estimated by linear regression using the above values of kinetic parameters and assuming an initial condition for serum TCDD of 7 ppt, i.e., the average level found in unexposed workers. The estimate of the occupational exposure parameter was 232.7 pg/kg/day (95% CI 192, 273). This value can be applied to the full cohort to obtain for each cohort member the time course of serum dioxin concentration from which exposure indices can be derived. Sensitivity coefficients to model parameters (background input, kf, occupational exposure, and the assumed TCDD concentration at hire) allow for a convenient recalculation of the serum TCDD curve and of the derived exposure indices for different assumed values of the model parameters.

Power spectral analysis of heart rate in hypothyroidism.

OBJECTIVE: The aim of the present study was to evaluate the impact of hypothyroidism on the autonomic regulation of the cardiovascular system by analysing separately sympathetic and parasympathetic influences on the heart. DESIGN: In seven newly diagnosed untreated hypothyroid patients we analysed power spectral density of heart rate cyclic variations at rest, while lying, and while standing. The same protocol was repeated after the induction of stable euthyroidism by levothyroxine (L-T(4)) treatment. The results were also compared with those obtained from seven age-, sex- and body mass index-matched control subjects. METHODS: Heart rate variability was evaluated by autoregressive power spectral analysis (PSA). This method allows reliable quantification of low frequency (LF) and high frequency (HF) components of the heart rate power spectral density. These are considered to be under mainly sympathetic and purely parasympathetic control respectively. In addition, heart rate variations during deep breathing, lying to standing, and Valsalva's manoeuvre were assessed. RESULTS: PSA showed a sharp reduction in the HF (parasympathetic) component in hypothyroid subjects compared with controls (lying, 29.4+/-5.4 vs 47.7+/-6.3 normalized units (NU) (means +/- s.e.m.), P<0.05; standing, 14.0+/-3.5 vs 32.1+/-3.6NU, P<0.005). Conversely, the LF (mainly sympathetic) component was higher in hypothyroid subjects than in controls (lying, 61.6+/-6.4 vs 45.4+/-6.7 NU; standing, 71.7+/-8.0 vs 53.1+/-5.6NU), this difference being significant in the standing position. Hence, the LF/HF ratio, which is considered an index of sympathovagal balance, was increased in hypothyroid subjects while both lying (2.75+/-0.6 vs 1.16+/-0.3; P<0.05) and standing (10.0+/-3.7 vs 1.85+/-0.3; P<0. 02). Total heart rate variability, expressed as total power spectral density, was lower in hypothyroid patients than in control subjects, this difference being significant in the lying position (574+/-126 vs 2302+/-994ms(2), P<0.05). In patients re-examined after L-T(4) therapy, complete normalization of cardiovascular parameters was observed (LF/HF ratio, lying, 1.26+/-0.4; standing, 2.56+/-0.8, both P<0.01 vs baseline values). The response to conventional cardiovascular autonomic tests was not significantly different between hypothyroid patients and healthy controls, and did not change in patients after therapy. CONCLUSIONS: These results suggest that, contrary to the clinical picture, thyroid hormone deficiency is associated with an increased sympathetic influence on the autonomic cardiovascular system. The changes in sympathetic function could be explained by a secondary adaptation to an altered cardiovascular responsiveness.

Uncertainty in estimating exposure using a toxicokinetic model. The example of 2,3,7,8-tetrachlorodibenzo-p-dioxin.

This paper deals with sources of uncertainty in the use of a minimal physiological toxicokinetic model to obtain dose estimates for a dose-response analysis of cancer in an occupational cohort. Toxicokinetic models make it possible to construct exposure parameters that are more closely related to the individual dose than traditional measures of exposures to toxic agents. However, the process introduces a wide array of sources of uncertainty. Selecting a model structure to describe the kinetics of a toxic agent implies necessarily making simplifications and assumptions that influence the range of applicability of the model. Once a model has been selected, the value of certain model parameters (constants) must be assigned, for example, from anthropometric data. The question then arises of how sensitive the model predictions are to variations in the values of these constants. Other model parameters, typically those describing the kinetics of the agent, are next estimated from actual data. There may be limitations in the data concerning, for example, sparseness (too few observations per subject) or missing values. The methods used for parameter estimation carry their own set of assumptions that need to be appropriate to the situation at hand. In summary, the dioxin example is used to characterize the sources of uncertainty at different levels, such as model structure, methods and data used for parameter estimation, estimation of occupational exposure, and imputation of missing values in exposure indices derived from the kinetic model.

Beta-cell activity and hepatic insulin extraction following dexamethasone administration in healthy subjects.

Glucocorticoids induce an increase of hepatic glucose production and peripheral resistance to insulin action. It is further assumed that dexamethasone administration in humans causes insulin hypersecretion, although inferences on beta-cell activity have been made in absolute terms and mostly from observations of systemic insulin concentration. In fact, the role of hepatic insulin extraction in humans treated long-term with glucocorticoids has not been investigated. The aim of the present study was to factor out quantitatively the main components of the insulin pathway that are responsible for the peripheral hypersecretion observed after steroids. Frequently sampled intravenous (FSIGT) and oral (OGTT) glucose tolerance tests were performed in healthy subjects before and after 5 days of oral dexamethasone administration (4 mg/d). Insulin sensitivity, beta-cell secretion, and hepatic insulin extraction were estimated by means of mathematical modeling. After steroids, insulin sensitivity decreased from 6.00 +/- 1.29 to 4.23 +/- 1.04 min-1/(microU/mL) (P < .04). Basal beta-cell secretion increased from 45 +/- 7 to 104 +/- 26 pmol/L . min-1 (P < .004) during the FSIGT and from 40 +/- 6 to 88 +/- 21 (P < .05) during the OGTT; total insulin release increased from 19 +/- 5 to 36 +/- 7 nmol/L in 180 minutes (P < .005) and from 33 +/- 5 to 50 +/- 10 (P < .02), respectively, FSIGT data also showed that first-phase beta-cell sensitivity increased from 236 +/- 39 to 309 +/- 33 pmol/L . min-1/(mg/dL) (P < .04), and second-phase from 631 +/- 154 to 1,103 +/ 196 10(4) pmol/L . min-2/(mg/dL) (P < .03). Posthepatic insulin delivery increased only insignificantly during the FSIGT (from 3.4 +/- 0.6 to 4.5 +/- 0.5 nmol/L, P = .073) due to an augmented hepatic insulin extraction from 73.0% +/- 7.2% to 83.0% +/- 3.5% (P < .05). During the OGTT, posthepatic insulin delivery increased after treatment from 6.6 +/- 1.2 to 11.4 +/- 2.5 nmol/L (P < .035) due to an increase, although slight, of hepatic insulin extraction from 77.4% +/- 1.9% to 79.3% +/- 3.3% (P = .319). The increased overall beta-cell activity during both tests was observed also by analyzing OGTT profiles of islet amyloid polypeptide (IAPP), the secretion of which was higher after steroids (basal, 0.081 +/- 0.012 v 0.272 +/- 0.082 pmol/L/min, P < .02; total, 35 +/- 8 v 116 +/- 48 mpmol/L in 3 hours, P < .05). In conclusion, after dexamethasone administration, peripheral hyperinsulinemia due to marked prehepatic beta-cell insulin hypersecretion is partially ameliorated by a concomitant increase of hepatic insulin clearance, which is more evident during a FSIGT. Model-derived secretion parameters from the OGTT and FSIGT produced comparable results, indicating that both tests, when properly analyzed, are feasible tools to evaluate insulin secretion.

Zinc supplementation improves glucose disposal in patients with cirrhosis.

Zinc deficiency is common in cirrhosis, and was proved to affect nitrogen metabolism. In experimental animals, zinc status may also affect glucose disposal, and acute zinc supplementation improves glucose tolerance in healthy subjects. This study was aimed at measuring the effects of long-term oral zinc supplements on glucose tolerance in cirrhosis. The time courses of glucose, insulin, and C-peptide in response to an intravenous (i.v.) glucose load were analyzed by the minimal-model technique before and after long-term oral zinc supplements (200 mg three times per day for 60 days) in 10 subjects with advanced cirrhosis and impaired glucose tolerance or diabetes. The test was performed using a simplified procedure, based on 20 blood samples collected within 4 hours from the glucose load. Normal values were obtained in 25 age-matched healthy subjects. Zinc levels were low to normal or reduced before treatment, and were normalized by oral zinc. Glucose disappearance improved by greater than 30% in response to treatment. There were no changes in pancreatic insulin secretion and systemic delivery, or in the hepatic extraction of insulin. Insulin sensitivity (SI), which was reduced by 80% before treatment, did not change. Glucose effectiveness (SG) was nearly halved in cirrhosis before treatment (0.013 [SD 0.007] min(-1) v. 0.028 [SD 0.009] in controls; P < .001), and increased to 0.017 (SD 0.009) after zinc (P < .05 v. baseline). The return to normal of plasma zinc levels after long-term zinc treatment in advanced cirrhosis improves glucose tolerance via an increase of the effects of glucose per se on glucose metabolism. Poor zinc status may contribute to the impaired glucose tolerance and diabetes of cirrhosis.

Urea as a marker of adequacy in hemodialysis: lesson from in vivo urea dynamics monitoring.

BACKGROUND: "Dialysis dose," a concept developed by Sargent and Gotch based on urea kinetic modeling, is a useful and recognized tool that is used to quantitate and optimize a dialysis-efficacy program. However, it has been shown that oversimplification of the "dialysis adequacy" concept to the Kt/V index might lead to dramatic underdialysis and subsequent deleterious consequences on morbidity and mortality of dialysis patients. With this perspective, the determination of Kt/V must be very cautious and rely on accurate measurement of postdialysis urea concentration and its use integrated as a tool in a quality-assurance process. METHODS: In this study, we analyzed urea dynamics by means of a blood side (ultrafiltrate) continuous online urea monitoring system interfaced with a two-pool model hosted in a microcomputer. The study was designed to provide instantaneous dialysis performances (body and dialyzer clearances, dialyzer mass transfer coefficient) and to determine the in vivo functional permeability characteristics of the patient [intercompartment urea mass transfer coefficient (Kc)]. Thirteen end-stage renal disease patients (age 54 +/- 16 years; 12 male and 1 female) were studied during nine consecutive dialysis sessions (3 weeks). RESULTS: Urea kinetics obtained from the urea monitoring system fitted closely the urea kinetic modeling prediction, confirming the validity of the double-pool model structure. Effective in vivo urea mass transfer coefficient averaged 912 +/- 235 mL/min/1.73 m2, a value close to those reported with more invasive methods. Large variations ranging from 363 to 1249 mL/min were observed among patients, confirming very large interindividual patient permeability differences. Interestingly, the urea mass transfer coefficient was inversely correlated with the postdialysis rebound values. Intraindividual variations were also noted as a function of time denoting functional changes in urea mass transfer coefficient values. The urea distribution volume was 38.1 +/- 7, 8 L (53 +/- 8% body weight). V1 referring to the extracellular volume and V2 to the intracellular volume were 9 +/- 2 L (13 +/- 2% body weight) and 29.2 +/- 6.6 L (41 +/- 1.3% body wt), respectively. The extracellular/intracellular volume ratio was 0.31 (approximately one third) and was not as usually defined by the paradigm 1/2 ratio. CONCLUSION: Online double-pool urea kinetic modeling gave a new insight in urea kinetic modeling approach. Urea dynamics fit perfectly a double-compartment model structure. Accessible extracellular volume to hemodialysis is smaller than expected. The in vivo urea mass transfer coefficient must be considered as an individual and variable characteristic of ESRD patients that should be taken into consideration when prescribing the hemodialysis schedule.

Analysis of errors due to assigned model constants in pharmacokinetic calculations based on functional optimization.

Mathematical models in pharmacokinetic studies frequently include, in addition to model parameters that are estimated from experimental data, model constants that are directly measured or derived from standardized, anthropometric formulae. Although these assigned constants are usually considered as error-free, uncertainties in their values affect the estimation precision of the remaining model parameters as well as the reliability of other derived quantities such as predictions of drug concentration levels or prescriptions of drug dosage. The present paper addresses the problem of quantifying the effects of errors in assigned model constants on the outcome of pharmacokinetic calculations. The error analysis is performed starting from the general assumption that the optimization of some objective function, such as encountered in parameter estimation, is involved. The proposed approach is based on the linearization of theoretical conditions that must be satisfied by the optimal solution of the problems. It requires quite complex mathematical calculations which can be however performed automatically by specific software based on symbolic computer algebra. Examples of applications are reported showing the usefulness of this kind of error analysis, in particular for detecting potentially critical error sources in the modeling of physiological systems and for suggesting improvements to experimental protocols for increasing the confidence in diagnostic and therapeutic decisions.

Use of a toxicokinetic model in the analysis of cancer mortality in relation to the estimated absorbed dose of dioxin (2,3,7,8-tetrachlorodibenzo-p-dioxin, TCDD).

We performed an analysis of All cancer and Lung cancer mortality in relation to estimated absorbed dose of dioxin (2,3,7,8-tetrachlorodibenzo-p-dioxin, TCDD) in the cohort of chemical workers at 12 US plants assembled by the US National Institute for Occupational Safety and Health (NIOSH) (n = 5172). Estimates of cumulative exposure to TCDD were based on a minimal physiologic toxicokinetic model (MPTK) that accounts for inter- and intra-individual variations in body mass index (BMI) over time. Population-level parameters related to liver elimination and background (input or concentration) of TCDD were estimated from separate data with repeated measures of serum TCDD (US Air Force Health Study). An occupational TCDD input parameter was estimated based on one-point-in-time TCDD data available for a subset (n = 253) of the NIOSH cohort. Model-based time-dependent cumulative dose estimates (area under the curve (AUC) of the lipid-adjusted serum TCDD concentration over time) were obtained for members of the full cohort with recorded body height and weight (n = 4049), as this information is required by the MPTK model to compute dose. Missing-value problems arose in the estimation of the occupational input parameter (n = 42) and in TCDD-dose calculation in the full cohort (n = 886) and they were handled with multiple imputation methods. Risk-regression analyses were based on Cox log-linear models including age at entry, year of entry and duration of employment as categorical covariates in addition to the logarithm of cumulative TCDD dose in ppt-years. Risk sets were stratified on birth cohort. Estimates of the unlagged exposure coefficient in these models were 0.1249 [95% confidence interval (CI) 0.0144, 0.2354] for All cancer and 0.2158 (95% CI 0.02376, 0.4078) for lung cancer. A 10-year lag produced an increase in the estimate for all cancer (0.1539, 95% CI 0.0387, 0.2691), whereas, the estimate for lung cancer was not affected much (0.2125, 95% CI 0.0138, 0.4112). At a dose level of 100 times the background the estimates obtained with a 10-year lag translate into a relative risk of 2.03 (95% CI 1.19-3.45) for all cancer and of 2.66 (95% CI 1.07-6.64) for lung cancer. Higher estimates of the exposure coefficients were obtained after imputation of missing values. This increase in risk seemed due to the inclusion of short-term workers, who may exhibit a higher mortality for reasons other than dioxin exposure.

Evaluation of the vagal-sympathetic interaction in diabetics with autonomic neuropathy through power spectrum density analysis of the heart rate. A critical revision of the natural history of diabetic autonomic neuropathy is possible.

In this paper we apply spectral analysis methods to heart rate variability to assess the autonomic nervous system activity in normal subjects and in patients affected by different degrees of diabetic autonomic neuropathy. The current opinion, based on different clinical tests, is that parasympathetic impairment occurs earlier in autonomic dysfunctions. However, the use of power spectrum density analysis based on a single parameter (heart rate) suggests a simultaneous involvement of parasympathetic and sympathetic pathways leading to the conclusion that perhaps the natural history of diabetic autonomic neuropathy should be substantially rewritten.

PANSYM: a symbolic equation generator for mathematical modelling, analysis and control of metabolic and pharmacokinetic systems.

Software is presented for automatic generation of first-order ordinary differential equations (ODE) that arise from lumped parameter representations of metabolic and pharmacokinetic systems. The definition of system structures is accomplished by fractional transfer rates between state variables, together with input/output equations and initial conditions of state variables. General non-linear mathematical expressions can be assigned to all structure definition items. The software parses and interprets the system definitions and generates symbolically the mathematical expression of the model's set of ODE. In addition, symbolic derivatives of state equations are determined with respect to model parameters, state variables and external inputs. These derivatives represent the constituents of systems of sensitivity-differential and adjoint-differential equations that arise in identification and optimal control problems. Finally, output routines generate source code that, once compiled and linked to simulation programs, allows efficient numerical integration of the system of ODE. This software has been developed in PROLOG on Macintosh computers and has been extensively used with the programming environment MATLAB. Possible applications of this software include model building, sensitivity analysis, identification, optimal experiment design and numerical solution of optimal control problems.

Design of optimal two-sample bolus injection tests for measuring low plasma clearance rates.

The bolus injection technique is widely used for assessing plasma clearance rates of substances and is based on multiexponential data analysis of multiple concentration measurements. This approach can be simplified to a monoexponential description with only two measurements if intercompartmental mixing occurs at a much faster rate than elimination, i.e. with low plasma clearance rates. In this context initial transients must be ignored because the first measurement time instant affects the accuracy of clearance estimates if mixing is still incomplete. Moreover, measurement noise affects estimation precision which can be optimized by suitable choice of the sampling schedule and of the injected test dose. The aim of this study is to design two-sample bolus injection tests for measuring low plasma clearance rates with preassigned precision and with minimum amount of injected test substance. This paper provides equations for evaluating whether plasma clearance is sufficiently low to allow a monoexponential description of the plasma disappearance curve and for choosing the first sampling time instant. Closed form equations are proposed for determining the optimal test dose and the second sampling time. Result are derived for a particular heteroscedastic measurement noise description, and the problem of robustness with respect to interpatient variability of kinetic parameters is addressed.