A compartmental model of hepatic disposition kinetics: 1. Model development and application to linear kinetics

The conventional convection-dispersion model is widely used to interrelate hepatic availability (F) and clearance (Cl) with the morphology and physiology of the liver and to predict effects such as changes in liver blood flow on F and Cl. The extension of this model to include nonlinear kinetics and zonal heterogeneity of the liver is not straightforward and requires numerical solution of partial differential equation, which is not available in standard nonlinear regression analysis software. In this paper, we describe an alternative compartmental model representation of hepatic disposition (including elimination). The model allows the use of standard software for data analysis and accurately describes the outflow concentration-time profile for a vascular marker after bolus injection into the liver. In an evaluation of a number of different compartmental models, the most accurate model required eight vascular compartments, two of them with back mixing. In addition, the model includes two adjacent secondary vascular compartments to describe the tail section of the concentration-time profile for a reference marker. The model has the added flexibility of being easy to modify to model various enzyme distributions and nonlinear elimination. Model predictions of F, MTT, CV2, and concentration-time profile as well as parameter estimates for experimental data of an eliminated solute (palmitate) are comparable to those for the extended convection-dispersion model.

A two compartmental model of the respiratory system /

"UIUCDCS-R-73-616"

Mode locking in a spatially extended neuron model: active soma and compartmental tree

Understanding the mode-locked response of excitable systems to periodic forcing has important applications in neuroscience. For example it is known that spatially extended place cells in the hippocampus are driven by the theta rhythm to generate a code conveying information about spatial location. Thus it is important to explore the role of neuronal dendrites in generating the response to periodic current injection. In this paper we pursue this using a compartmental model, with linear dynamics for each compartment, coupled to an active soma model that generates action potentials. By working with the piece-wise linear McKean model for the soma we show how the response of the whole neuron model (soma and dendrites) can be written in closed form. We exploit this to construct a stroboscopic map describing the response of the spatially extended model to periodic forcing. A linear stability analysis of this map, together with a careful treatment of the non-differentiability of the soma model, allows us to construct the Arnol'd tongue structure for 1:q states (one action potential for q cycles of forcing). Importantly we show how the presence of quasi-active membrane in the dendrites can influence the shape of tongues. Direct numerical simulations confirm our theory and further indicate that resonant dendritic membrane can enlarge the windows in parameter space for chaotic behavior. These simulations also show that the spatially extended neuron model responds differently to global as opposed to point forcing. In the former case spatio-temporal patterns of activity within an Arnol'd tongue are standing waves, whilst in the latter they are traveling waves.

Malaria transmission model for different levels of acquired immunity and temperature-dependent parameters (vector)

OBJECTIVE: Describe the overall transmission of malaria through a compartmental model, considering the human host and mosquito vector. METHODS: A mathematical model was developed based on the following parameters: human host immunity, assuming the existence of acquired immunity and immunological memory, which boosts the protective response upon reinfection; mosquito vector, taking into account that the average period of development from egg to adult mosquito and the extrinsic incubation period of parasites (transformation of infected but non-infectious mosquitoes into infectious mosquitoes) are dependent on the ambient temperature. RESULTS: The steady state equilibrium values obtained with the model allowed the calculation of the basic reproduction ratio in terms of the model's parameters. CONCLUSIONS: The model allowed the calculation of the basic reproduction ratio, one of the most important epidemiological variables.

Steady-state brain glucose transport kinetics re-evaluated with a four-state conformational model.

Glucose supply from blood to brain occurs through facilitative transporter proteins. A near linear relation between brain and plasma glucose has been experimentally determined and described by a reversible model of enzyme kinetics. A conformational four-state exchange model accounting for trans-acceleration and asymmetry of the carrier was included in a recently developed multi-compartmental model of glucose transport. Based on this model, we demonstrate that brain glucose (G(brain)) as function of plasma glucose (G(plasma)) can be described by a single analytical equation namely comprising three kinetic compartments: blood, endothelial cells and brain. Transport was described by four parameters: apparent half saturation constant K(t), apparent maximum rate constant T(max), glucose consumption rate CMR(glc), and the iso-inhibition constant K(ii) that suggests G(brain) as inhibitor of the isomerisation of the unloaded carrier. Previous published data, where G(brain) was quantified as a function of plasma glucose by either biochemical methods or NMR spectroscopy, were used to determine the aforementioned kinetic parameters. Glucose transport was characterized by K(t) ranging from 1.5 to 3.5 mM, T(max)/CMR(glc) from 4.6 to 5.6, and K(ii) from 51 to 149 mM. It was noteworthy that K(t) was on the order of a few mM, as previously determined from the reversible model. The conformational four-state exchange model of glucose transport into the brain includes both efflux and transport inhibition by G(brain), predicting that G(brain) eventually approaches a maximum concentration. However, since K(ii) largely exceeds G(plasma), iso-inhibition is unlikely to be of substantial importance for plasma glucose below 25 mM. As a consequence, the reversible model can account for most experimental observations under euglycaemia and moderate cases of hypo- and hyperglycaemia.

Understanding post-operative temperature drop in cardiac surgery: a mathematical model

A mathematical model is presented to understand heat transfer processes during the cooling and re-warming of patients during cardiac surgery. Our compartmental model is able to account for many of the qualitative features observed in the cooling of various regions of the body including the central core containing the majority of organs, the rectal region containing the intestines and the outer peripheral region of skin and muscle. In particular, we focus on the issue of afterdrop: a drop in core temperature following patient re-warming, which can lead to serious post-operative complications. Model results for a typical cooling and re-warming procedure during surgery are in qualitative agreement with experimental data in producing the afterdrop effect and the observed dynamical variation in temperature between the core, rectal and peripheral regions. The influence of heat transfer processes and the volume of each compartmental region on the afterdrop effect is discussed. We find that excess fat on the peripheral and rectal regions leads to an increase in the afterdrop effect. Our model predicts that, by allowing constant re-warming after the core temperature has been raised, the afterdrop effect will be reduced.

An insulin infusion advisory system for type 1 diabetes patients based on non-linear model predictive control methods

In this paper, an Insulin Infusion Advisory System (IIAS) for Type 1 diabetes patients, which use insulin pumps for the Continuous Subcutaneous Insulin Infusion (CSII) is presented. The purpose of the system is to estimate the appropriate insulin infusion rates. The system is based on a Non-Linear Model Predictive Controller (NMPC) which uses a hybrid model. The model comprises a Compartmental Model (CM), which simulates the absorption of the glucose to the blood due to meal intakes, and a Neural Network (NN), which simulates the glucose-insulin kinetics. The NN is a Recurrent NN (RNN) trained with the Real Time Recurrent Learning (RTRL) algorithm. The output of the model consists of short term glucose predictions and provides input to the NMPC, in order for the latter to estimate the optimum insulin infusion rates. For the development and the evaluation of the IIAS, data generated from a Mathematical Model (MM) of a Type 1 diabetes patient have been used. The proposed control strategy is evaluated at multiple meal disturbances, various noise levels and additional time delays. The results indicate that the implemented IIAS is capable of handling multiple meals, which correspond to realistic meal profiles, large noise levels and time delays.

Effects of informative censoring on the hazard function: Application to the proportional hazards regression model

The standard analyses of survival data involve the assumption that survival and censoring are independent. When censoring and survival are related, the phenomenon is known as informative censoring. This paper examines the effects of an informative censoring assumption on the hazard function and the estimated hazard ratio provided by the Cox model.^ The limiting factor in all analyses of informative censoring is the problem of non-identifiability. Non-identifiability implies that it is impossible to distinguish a situation in which censoring and death are independent from one in which there is dependence. However, it is possible that informative censoring occurs. Examination of the literature indicates how others have approached the problem and covers the relevant theoretical background.^ Three models are examined in detail. The first model uses conditionally independent marginal hazards to obtain the unconditional survival function and hazards. The second model is based on the Gumbel Type A method for combining independent marginal distributions into bivariate distributions using a dependency parameter. Finally, a formulation based on a compartmental model is presented and its results described. For the latter two approaches, the resulting hazard is used in the Cox model in a simulation study.^ The unconditional survival distribution formed from the first model involves dependency, but the crude hazard resulting from this unconditional distribution is identical to the marginal hazard, and inferences based on the hazard are valid. The hazard ratios formed from two distributions following the Gumbel Type A model are biased by a factor dependent on the amount of censoring in the two populations and the strength of the dependency of death and censoring in the two populations. The Cox model estimates this biased hazard ratio. In general, the hazard resulting from the compartmental model is not constant, even if the individual marginal hazards are constant, unless censoring is non-informative. The hazard ratio tends to a specific limit.^ Methods of evaluating situations in which informative censoring is present are described, and the relative utility of the three models examined is discussed. ^

Modelos de transferência de metal pesado na cana-de-açúcar adubada com composto de lixo urbano

The research approaches recycling of urban waste compost (UWC) as an alternative fertilizer for sugarcane crop and as a social and environmental solution to the solids residuals growth in urban centers. A mathematical model was used in order to know the metal dynamics as decision support tool, aiming to establish of criteria and procedures for UWC's safe use, limited by the amount of heavy metal. A compartmental model was developed from experimental data in controlled conditions and partially checked with field data. This model described the heavy metal transference in the system soil-root-aerial portion of sugarcane plants and concluded that nickel was metal to be concern, since it takes approximately three years to be attenuated in the soil, reaching the aerial portions of the plant at high concentrations. Regarding factors such as clay content, oxide level and soil pH, it was observed that for soil with higher buffering capacity, the transfer of the majority of the metals was slower. This model may become an important tool for the attainment of laws regarding the UWC use, aiming to reduce environment contamination the waste accumulation and production costs.

Metformin pharmacokinetics in nondiabetic pregnant women with polycystic ovary syndrome

The use of metformin throughout gestation by pregnant women with polycystic ovary syndrome (PCOS) significantly reduces the number of first trimester spontaneous abortions and the rate of occurrence of gestational diabetes. The objective of this study was to investigate the pharmacokinetics and the placental transfer of metformin in pregnant women with PCOS. Eight pregnant women with PCOS taking 850 mg metformin every 12 h during the third trimester of pregnancy were evaluated. Maternal blood samples were collected at steady state during the dose interval (0-12 h). Maternal and umbilical cord blood samples were also obtained at delivery. Metformin plasma concentrations were analyzed by high-performance liquid chromatography, and pharmacokinetic parameters were determined using a non-compartmental model. Data are reported as median and minimum and maximum values. Metformin pharmacokinetic parameters were: t(A1/2), 3.8 (2.8-5.4) h; t(max), 2.0 (0.5-3.0) h; C(max), 1.4 (0.5-2.1) mg/L; C(mean), 0.5 (0.2-0.9) mg/L; AUC(0-12), 6.4 (1.1-9.2) mg h/L; Cl/f, 105 (60-274) L/h; Vd/f, 551 (385-1173) L; median fluctuation, 89 (79-95)%. Umbilical/maternal metformin plasma concentration ratios were 0.7 (0.4-1.3). Metformin oral clearance (Cl/f) had increased in our patients relative to nonpregnant healthy volunteers or diabetic patients. Therefore, lower plasma metformin concentrations were observed for nondiabetic pregnant women with PCOS. Future studies should be conducted to demonstrate the therapeutic efficacy of metformin during pregnancy. Caution is warranted as umbilical/maternal metformin plasma concentrations ratios of around 0.7 require metformin dosage adjustment.

Kinetics of propranolol uptake in muscle, skin, and fat of the isolated perfused rat hindlimb

The tissue distribution kinetics of a highly bound solute, propranolol, was investigated in a heterogeneous organ, the isolated perfused limb, using the impulse-response technique and destructive sampling. The propranolol concentration in muscle, skin, and fat as well as in outflow perfusate was measured up to 30 min after injection. The resulting data were analysed assuming (1) vascular, muscle, skin and fat compartments as well mixed (compartmental model) and (2) using a distributed-in-space model which accounts for the noninstantaneous intravascular mixing and tissue distribution processes but consists only of a vascular and extravascular phase (two-phase model). The compartmental model adequately described propranolol concentration-time data in the three tissue compartments and the outflow concentration-time curve (except of the early mixing phase). In contrast, the two-phase model better described the outflow concentration-time curve but is limited in accounting only for the distribution kinetics in the dominant tissue, the muscle. The two-phase model well described the time course of propranolol concentration in muscle tissue, with parameter estimates similar to those obtained with the compartmental model. The results suggest, first that the uptake kinetics of propranolol into skin and fat cannot be analysed on the basis of outflow data alone and, second that the assumption of well-mixed compartments is a valid approximation from a practical point of view las, e.g., in physiological based pharmacokinetic modelling). The steady-state distribution volumes of skin and fat were only 16 and 4%, respectively, of that of muscle tissue (16.7 ml), with higher partition coefficient in fat (6.36) than in skin (2.64) and muscle (2.79. (C) 2000 Elsevier Science B.V. All rights reserved.

Evaluation of a new bioelectrical impedance instrument for the prediction of body cell mass independently of height or weight

The objective of the present study was to evaluate the performance of a new bioelectrical impedance instrument, the Soft Tissue Analyzer (STA), which predicts a subject's body composition. A cross-sectional population study in which the impedance of 205 healthy adult subjects was measured using the STA. Extracellular water (ECW) volume (as a percentage of total body water, TBW) and fat-free mass (FFM) were predicted by both the STA and a compartmental model, and compared according to correlation and limits of agreement analysis, with the equivalent data obtained by independent reference methods of measurement (TBW measured by D2O dilution, and FFM measured by dual-energy X-ray absorptiometry). There was a small (2.0 kg) but significant (P < 0.02) difference in mean FFM predicted by the STA, compared with the reference technique in the males, but not in the females (-0.4 kg) or in the combined group (0.8 kg). Both methods were highly correlated. Similarly, small but significant differences for predicted mean ECW volume were observed. The limits of agreement for FFM and ECW were -7.5-9.9 and -4.1-3.0 kg, respectively. Both FFM and ECW (as a percentage of TBW) are well predicted by the STA on a population basis, but the magnitude of the limits of agreement with reference methods may preclude its usefulness for predicting body composition in an individual. In addition, the theoretical basis of an impedance method that does not include a measure of conductor length requires further validation. (C) Elsevier Science Inc. 2000.

Methimazole and propylthiouracil equally cross the perfused human term placental lobule

Propylthiouracil (PTU) is widely believed to cross the placenta less freely than methimazole (MMI) and is therefore regarded as the preferred drug for treatment of hyperthyroidism in pregnancy. Clinical studies comparing the two drugs show, however, no differences in maternal or fetal thyroid function. We investigated transfer from the maternal to the fetal circuit in the isolated perfused term human placental lobule of low and high doses of PTU (4 mu g/mL and 40 mu g/mL) and MMI(1.5 mu g/mL and 15 mu g/mL) in protein-free perfusate and low doses of both drugs with addition of 40 g/L of bovine albumin. Both drugs readily crossed the placenta, reaching equilibrium in all experiments in about 2 h. Drug concentrations in the two circuits fitted a two compartmental model. Transfer kinetics for the two drugs were similar, nonsaturable, and unaffected by addition of albumin. Clearances (mL.min(-1).g(-1), means +/- SD) of PTU from maternal to fetal circuits were: 0.229 +/- 0.110, 0.216 +/- 0.065, and 0.170 +/- 0.032; and for transfer of MMI: 0.165 +/- 0.025, 0.232 +/- 0.153, and 0.174 +/- 0.009 (for low doses without, low doses with, and high doses without albumin, respectively). Clearances of PTU from fetal to maternal circuits were: 0.147 +/- 0.072, 0.109 +/- 0.014, and 0.116 +/- 0.028; and for transfer of MMI: 0.095 +/- 0.029, 0.122 +/- 0.088, and 0.12 +/- 0.005 (in the same experiments). There was no significant difference between drugs or drug doses and no effect of addition of albumin. We conclude that PTU and MMI have similar placental transfer kinetics.

Capillaries within compartments: Microvascular interpretation of dynamic positron emission tomography data

Measurement of exchange of substances between blood and tissue has been a long-lasting challenge to physiologists, and considerable theoretical and experimental accomplishments were achieved before the development of the positron emission tomography (PET). Today, when modeling data from modern PET scanners, little use is made of earlier microvascular research in the compartmental models, which have become the standard model by which the vast majority of dynamic PET data are analysed. However, modern PET scanners provide data with a sufficient temporal resolution and good counting statistics to allow estimation of parameters in models with more physiological realism. We explore the standard compartmental model and find that incorporation of blood flow leads to paradoxes, such as kinetic rate constants being time-dependent, and tracers being cleared from a capillary faster than they can be supplied by blood flow. The inability of the standard model to incorporate blood flow consequently raises a need for models that include more physiology, and we develop microvascular models which remove the inconsistencies. The microvascular models can be regarded as a revision of the input function. Whereas the standard model uses the organ inlet concentration as the concentration throughout the vascular compartment, we consider models that make use of spatial averaging of the concentrations in the capillary volume, which is what the PET scanner actually registers. The microvascular models are developed for both single- and multi-capillary systems and include effects of non-exchanging vessels. They are suitable for analysing dynamic PET data from any capillary bed using either intravascular or diffusible tracers, in terms of physiological parameters which include regional blood flow. (C) 2003 Elsevier Ltd. All rights reserved.

Assessing the effects of global warming and local social and economic conditions on the malaria transmission

OBJECTIVE: To show how a mathematical model can be used to describe and to understand the malaria transmission. METHODS: The effects on malaria transmission due to the impact of the global temperature changes and prevailing social and economic conditions in a community were assessed based on a previously presented compartmental model, which describes the overall transmission of malaria. RESULTS/CONCLUSIONS: The assessments were made from the scenarios produced by the model both in steady state and dynamic analyses. Depending on the risk level of malaria, the effects on malaria transmission can be predicted by the temperature ambient or local social and-economic conditions.