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.

The blood volumes of the primary and secondary circulatory system in the Atlantic cod Gadus morhua L., using plasma bound Evans Blue and compartmental analysis

The volume of the primary (PCS) and secondary (SCS) circulatory system in the Atlantic cod Gadus morhua was determined using a modified dye dilution technique. Cod (N=10) were chronically cannulated in the second afferent branchial artery with PE-50 tubing. Evans Blue dye was bound to harvested fish plasma at a concentration of 1 mg dye ml(-1) plasma, and injected at a concentration of 1 mg kg(-1) body mass. Serial sampling from the cannula produced a dye dilution curve, which could be described by a double exponential decay equation. Curve analysis enabled the calculation of the primary circulatory and total distribution volume. The difference between these volumes is assumed to be the volume of the SCS. From the dilution curve, it was also possible to calculate flow rates between and within the systems. The results of these experiments suggest a plasma volume in the PCS of 3.42+/-0.89 ml 100 g(-1) body mass, and in the SCS of 1.68+/-0.35 ml 100 g(-1) body mass (mean +/- S.D.) or approximately 50% that of the PCS. Flow rates to the SCS were calculated as 2.7% of the resting cardiac output. There was an allometric relationship between body mass and blood volumes. Increasing condition factor showed a tendency towards smaller blood volumes of the PCS, expressed as percentage body mass, but this was not evident for the volume of the SCS.

Sex differences in urinary levels of several biological indicators of exposure: a simulation study using a compartmental based toxicokinetic model

Toxicokinetic modeling is a useful tool to describe or predict the behavior of a chemical agent in the human or animal organism. A general model based on four compartments was developed in a previous study in order to quantify the effect of human variability on a wide range of biological exposure indicators. The aim of this study was to adapt this existing general toxicokinetic model to three organic solvents, which were methyl ethyl ketone, 1-methoxy-2-propanol and 1,1,1,-trichloroethane, and to take into account sex differences. We assessed in a previous human volunteer study the impact of sex on different biomarkers of exposure corresponding to the three organic solvents mentioned above. Results from that study suggested that not only physiological differences between men and women but also differences due to sex hormones levels could influence the toxicokinetics of the solvents. In fact the use of hormonal contraceptive had an effect on the urinary levels of several biomarkers, suggesting that exogenous sex hormones could influence CYP2E1 enzyme activity. These experimental data were used to calibrate the toxicokinetic models developed in this study. Our results showed that it was possible to use an existing general toxicokinetic model for other compounds. In fact, most of the simulation results showed good agreement with the experimental data obtained for the studied solvents, with a percentage of model predictions that lies within the 95% confidence interval varying from 44.4 to 90%. Results pointed out that for same exposure conditions, men and women can show important differences in urinary levels of biological indicators of exposure. Moreover, when running the models by simulating industrial working conditions, these differences could even be more pronounced. In conclusion, a general and simple toxicokinetic model, adapted for three well known organic solvents, allowed us to show that metabolic parameters can have an important impact on the urinary levels of the corresponding biomarkers. These observations give evidence of an interindividual variablity, an aspect that should have its place in the approaches for setting limits of occupational exposure.

Changes in cerebral compartmental compliances during mild hypocapnia in patients with traumatic brain injury.

The benefit of induced hyperventilation for intracranial pressure (ICP) control after severe traumatic brain injury (TBI) is controversial. In this study, we investigated the impact of early and sustained hyperventilation on compliances of the cerebral arteries and of the cerebrospinal (CSF) compartment during mild hyperventilation in severe TBI patients. We included 27 severe TBI patients (mean 39.5 ± 3.4 years, 6 women) in whom an increase in ventilation (20% increase in respiratory minute volume) was performed during 50 min as part of a standard clinical CO(2) reactivity test. Using a new mathematical model, cerebral arterial compliance (Ca) and CSF compartment compliance (Ci) were calculated based on the analysis of ICP, arterial blood pressure, and cerebral blood flow velocity waveforms. Hyperventilation initially induced a reduction in ICP (17.5 ± 6.6 vs. 13.9 ± 6.2 mmHg; p < 0.001), which correlated with an increase in Ci (r(2) = 0.213; p = 0.015). Concomitantly, the reduction in cerebral blood flow velocities (CBFV, 74.6 ± 27.0 vs. 62.9 ± 22.9 cm/sec; p < 0.001) marginally correlated with the reduction in Ca (r(2) = 0.209; p = 0.017). During sustained hyperventilation, ICP increased (13.9 ± 6.2 vs. 15.3 ± 6.4 mmHg; p < 0.001), which correlated with a reduction in Ci (r(2) = 0.297; p = 0.003), but no significant changes in Ca were found during that period. The early reduction in Ca persisted irrespective of the duration of hyperventilation, which may contribute to the lack of clinical benefit of hyperventilation after TBI. Further studies are needed to determine whether monitoring of arterial and CSF compartment compliances may detect and prevent an adverse ischemic event during hyperventilation.

Further analysis of open-respirometry systems: an a-compartmental mechanistic approach

A system is said to be "instantaneous" when for a given constant input an equilibrium output is obtained after a while. In the meantime, the output is changing from its initial value towards the equilibrium one. This is the transient period of the system and transients are important features of open-respirometry systems. During transients, one cannot compute the input amplitude directly from the output. The existing models (e.g., first or second order dynamics) cannot account for many of the features observed in real open-respirometry systems, such as time lag. Also, these models do not explain what should be expected when a system is speeded up or slowed down. The purpose of the present study was to develop a mechanistic approach to the dynamics of open-respirometry systems, employing basic thermodynamic concepts. It is demonstrated that all the main relevant features of the output dynamics are due to and can be adequately explained by a distribution of apparent velocities within the set of molecules travelling along the system. The importance of the rate at which the molecules leave the sensor is explored for the first time. The study approaches the difference in calibrating a system with a continuous input and with a "unit impulse": the former truly reveals the dynamics of the system while the latter represents the first derivative (in time) of the former and, thus, cannot adequately be employed in the apparent time-constant determination. Also, we demonstrate why the apparent order of the output changes with volume or flow.

A general framework for multi-compartmental analysis of drug chemotherapy dynamics in human immunodeficiency virus type-1 infected individuals

An optimal control strategy for the highly active antiretroviral therapy associated to the acquired immunodeficiency syndrome should be designed regarding a comprehensive analysis of the drug chemotherapy behavior in the host tissues, from major viral replication sites to viral sanctuary compartments. Such approach is critical in order to efficiently explore synergistic, competitive and prohibitive relationships among drugs and, hence, therapy costs and side-effect minimization. In this paper, a novel mathematical model for HIV-1 drug chemotherapy dynamics in distinct host anatomic compartments is proposed and theoretically evaluated on fifteen conventional anti-retroviral drugs. Rather than interdependence between drug type and its concentration profile in a host tissue, simulated results suggest that such profile is importantly correlated with the host tissue under consideration. Furthermore, the drug accumulative dynamics are drastically affected by low patient compliance with pharmacotherapy, even when a single dose lacks. (C) 2012 Elsevier Inc. All rights reserved.

Effects of Hypoxia on 13NH4+ Fluxes in Rice Roots1 : Kinetics and Compartmental Analysis

Techniques of compartmental (efflux) and kinetic influx analyses with the radiotracer 13NH4+ were used to examine the adaptation to hypoxia (15, 35, and 50% O2 saturation) of root N uptake and metabolism in 3-week-old hydroponically grown rice (Oryza sativa L., cv IR72) seedlings. A time-dependence study of NH4+ influx into rice roots after onset of hypoxia (15% O2) revealed an initial increase in the first 1 to 2.5 h after treatment imposition, followed by a decline to less than 50% of influx in control plants by 4 d. Efflux analyses conducted 0, 1, 3, and 5 d after the treatment confirmed this adaptation pattern of NH4+ uptake. Half-lives for NH4+ exchange with subcellular compartments, cytoplasmic NH4+ concentrations, and efflux (as percentage of influx) were unaffected by hypoxia. However, significant differences were observed in the relative amounts of N allocated to NH4+ assimilation and the vacuole versus translocation to the shoot. Kinetic experiments conducted at 100, 50, 35, and 15% O2 saturation showed no significant change in the Km value for NH4+ uptake with varying O2 supply. However, Vmax was 42% higher than controls at 50% O2 saturation, unchanged at 35%, and 10% lower than controls at 15% O2. The significance of these flux adaptations is discussed.

A two compartmental model of the respiratory system /

"UIUCDCS-R-73-616"

Estimation of pharmacokinetic parameters from non-compartmental variables using Microsoft Excel((R))

This study was conducted to develop a method, termed 'back analysis (BA)', for converting non-compartmental variables to compartment model dependent pharmacokinetic parameters for both one- and two-compartment models. A Microsoft Excel((R)) spreadsheet was implemented with the use of Solver((R)) and visual basic functions. The performance of the BA method in estimating pharmacokinetic parameter values was evaluated by comparing the parameter values obtained to a standard modelling software program, NONMEM, using simulated data. The results show that the BA method was reasonably precise and provided low bias in estimating fixed and random effect parameters for both one- and two-compartment models. The pharmacokinetic parameters estimated from the BA method were similar to those of NONMEM estimation.

A compartmental CFD-PBM model of high shear wet granulation

The conventional, geometrically lumped description of the physical processes inside a high shear granulator is not reliable for process design and scale-up. In this study, a compartmental Population Balance Model (PBM) with spatial dependence is developed and validated in two lab-scale high shear granulation processes using a 1.9L MiPro granulator and 4L DIOSNA granulator. The compartmental structure is built using a heuristic approach based on computational fluid dynamics (CFD) analysis, which includes the overall flow pattern, velocity and solids concentration. The constant volume Monte Carlo approach is implemented to solve the multi-compartment population balance equations. Different spatial dependent mechanisms are included in the compartmental PBM to describe granule growth. It is concluded that for both cases (low and high liquid content), the adjustment of parameters (e.g. layering, coalescence and breakage rate) can provide a quantitative prediction of the granulation process.

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.

The Basic Reproduction Number Obtained From Jacobian And Next Generation Matrices - A Case Study Of Dengue Transmission Modelling.

The basic reproduction number is a key parameter in mathematical modelling of transmissible diseases. From the stability analysis of the disease free equilibrium, by applying Routh-Hurwitz criteria, a threshold is obtained, which is called the basic reproduction number. However, the application of spectral radius theory on the next generation matrix provides a different expression for the basic reproduction number, that is, the square root of the previously found formula. If the spectral radius of the next generation matrix is defined as the geometric mean of partial reproduction numbers, however the product of these partial numbers is the basic reproduction number, then both methods provide the same expression. In order to show this statement, dengue transmission modelling incorporating or not the transovarian transmission is considered as a case study. Also tuberculosis transmission and sexually transmitted infection modellings are taken as further examples.

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.

Rumor Propagation Model: An Equilibrium Study

Compartmental epidemiological models have been developed since the 1920s and successfully applied to study the propagation of infectious diseases. Besides, due to their structure, in the 1960s an interesting version of these models was developed to clarify some aspects of rumor propagation, considering that spreading an infectious disease or disseminating information is analogous phenomena. Here, in an analogy with the SIR (Susceptible-Infected-Removed) epidemiological model, the ISS (Ignorant-Spreader-Stifler) rumor spreading model is studied. By using concepts from the Dynamical Systems Theory, stability of equilibrium points is established, according to propagation parameters and initial conditions. Some numerical experiments are conducted in order to validate the model.