Modeling of hepatic elimination and organ distribution kinetics with the extended convection-dispersion model

The conventional convection-dispersion (also called axial 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. An extended form of the convection-dispersion model has been developed to adequately describe the outflow concentration-time profiles for vascular markers at both short and long times after bolus injections into perfused livers. The model, based on flux concentration and a convolution of catheters and large vessels, assumes that solute elimination in hepatocytes follows either fast distribution into or radial diffusion in hepatocytes. The model includes a secondary vascular compartment, postulated to be interconnecting sinusoids. Analysis of the mean hepatic transit time (MTT) and normalized variance (CV2) of solutes with extraction showed that the discrepancy between the predictions of MTT and CV2 for the extended and conventional models are essentially identical irrespective of the magnitude of rate constants representing permeability, volume, and clearance parameters, providing that there is significant hepatic extraction. In conclusion, the application of a newly developed extended convection-dispersion model has shown that the unweighted conventional convection-dispersion model can be used to describe the disposition of extracted solutes and, in particular, to estimate hepatic availability and clearance in booth experimental and clinical situations.

Commentary: Using the convection-dispersion model and transit time density functions in the analysis of organ distribution kinetics

The convection-dispersion model and its extended form have been used to describe solute disposition in organs and to predict hepatic availabilities. A range of empirical transit-time density functions has also been used for a similar purpose. The use of the dispersion model with mixed boundary conditions and transit-time density functions has been queried recently by Hisaka and Sugiyanaa in this journal. We suggest that, consistent with soil science and chemical engineering literature, the mixed boundary conditions are appropriate providing concentrations are defined in terms of flux to ensure continuity at the boundaries and mass balance. It is suggested that the use of the inverse Gaussian or other functions as empirical transit-time densities is independent of any boundary condition consideration. The mixed boundary condition solutions of the convection-dispersion model are the easiest to use when linear kinetics applies. In contrast, the closed conditions are easier to apply in a numerical analysis of nonlinear disposition of solutes in organs. We therefore argue that the use of hepatic elimination models should be based on pragmatic considerations, giving emphasis to using the simplest or easiest solution that will give a sufficiently accurate prediction of hepatic pharmacokinetics for a particular application. (C) 2000 Wiley-Liss Inc. and the American Pharmaceutical Association J Pharm Sci 89:1579-1586, 2000.

Distribution kinetics of solutes in the isolated in-situ perfused rat head using the multiple indicator dilution technique and a physiological two-barrier model

The purpose of this study was to determine the pharmacokinetics of [C-14]diclofenac, [C-14]salicylate and [H-3]clonidine using a single pass rat head perfusion preparation. The head was perfused with 3-[N-morpholino] propane-sulfonic acid-buffered Ringer's solution. Tc-99m-red blood cells and a drug were injected in a bolus into the internal carotid artery and collected from the posterior facial vein over 28 min. A two-barrier stochastic organ model was used to estimate the statistical moments of the solutes. Plasma, interstitial and cellular distribution volumes for the solutes ranged from 1.0 mL (diclofenac) to 1.6 mL (salicylate), 2.0 mL (diclofenac) to 4.2 mL (water) and 3.9 mL (salicylate) to 20.9 mL (diclofenac), respectively. A comparison of these volumes to water indicated some exclusion of the drugs from the interstitial space and salicylate from the cellular space. Permeability-surface area (PS) products calculated from plasma to interstitial fluid permeation clearances (CLPI) (range 0.02-0.40 mL s(-1)) and fractions of solute unbound in the perfusate were in the order: diclofenac>salicylate >clonidine>sucrose (from 41.8 to 0.10 mL s(-1)). The slow efflux of diclofenac, compared with clonidine and salicylate, may be related to its low average unbound fraction in the cells. This work accounts for the tail of disposition curves in describing pharmacokinetics in the head.

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.

Chlamydiaceae and Chlamydia-like organisms in the koala (Phascolarctos cinereus)--organ distribution and histopathological findings.

Chlamydial infections in koalas can cause life-threatening diseases leading to blindness and sterility. However, little is known about the systemic spread of chlamydiae in the inner organs of the koala, and data concerning related pathological organ lesions are limited. The aim of this study was to perform a thorough investigation of organs from 23 koalas and to correlate their histopathological lesions to molecular chlamydial detection. To reach this goal, 246 formalin-fixed and paraffin embedded organ samples from 23 koalas were investigated by histopathology, Chlamydiaceae real-time PCR and immunohistochemistry, ArrayTube Microarray for Chlamydiaceae species identification as well as Chlamydiales real-time PCR and sequencing. By PCR, two koalas were positive for Chlamydia pecorum whereas immunohistochemical labelling for Chlamydiaceae was detected in 10 tissues out of nine koalas. The majority of these (n=6) had positive labelling in the urogenital tract related to histopathological lesions such as cystitis, endometritis, pyelonephritis and prostatitis. Somehow unexpected was the positive labelling in the gastrointestinal tract including the cloaca as well as in lung and spleen indicating systemic spread of infection. Uncultured Chlamydiales were detected in several organs of seven koalas by PCR, and four of these suffered from plasmacytic enteritis of unknown aetiology. Whether the finding of Chlamydia-like organisms in the gastrointestinal tract is linked to plasmacytic enteritis is unclear and remains speculative. However, as recently shown in a mouse model, the gastrointestinal tract might play a role being the site for persistent chlamydial infections and being a source for reinfection of the genital tract.

Viral load, organ distribution, histopathological lesions, and cytokine mRNA expression in goats infected with a molecular clone of the caprine arthritis encephalitis virus

Caprine arthritis encephalitis virus (CAEV) is a lentivirus of goats that causes persistent infection characterized by the appearance of inflammatory lesions in various organs. To define the sites of persistence, 5 goats were infected with a molecular clone of CAEV, and the viral load was monitored by real-time-PCR and RT-PCR in different sites 8 years after infection. The lymph nodes proved to be an important virus reservoir, with moderate virus replication relative to what is reported for lentiviruses of primates. Mammary gland and milk cells were preferred sites of viral replication. The viral load varied significantly between animals, which points to an important role of the genetic background. We found a clear association between occurrence of histopathological lesions and viral load in specific sites. The mRNA expression analysis of several cytokines did not reveal differences between animals that could explain the considerable individual variations in viral load observed.

A cationic vaccine adjuvant based on a saturated quaternary ammonium lipid have different in vivo distribution kinetics and display a distinct CD4 T cell-inducing capacity compared to its unsaturated analog

Adjuvants are often composed of different constituents that can be divided into two groups based on their primary activity: the delivery system which carries and presents the vaccine antigen to antigen-presenting cells, and the immunostimulator that activates and modulates the ensuing immune response. Herein, we have investigated the importance of the delivery system and in particular its physical characteristics by comparing the delivery properties of two lipids which differ only in the degree of saturation of the acyl chains, rendering the liposomes either rigid (DDA, dimethyldioctadecylammonium) or highly fluid (DODA, dimethyldioleoylammonium) at physiological temperature. We show that these delivery systems are remarkably different in their ability to prime a Th1-directed immune response with the rigid DDA-based liposomes inducing a response more than 100 times higher compared to that obtained with the fluid DODA-based liposomes. Upon injection with a vaccine antigen, DDA-based liposomes form a vaccine depot that results in a continuous attraction of antigen-presenting cells that engulf a high amount of adjuvant and are subsequently efficiently activated as measured by an elevated expression of the co-stimulatory molecules CD40 and CD86. In contrast, the fluid DODA-based liposomes are more rapidly removed from the site of injection resulting in a lower up-regulation of co-stimulatory CD40 and CD86 molecules on adjuvant-positive antigen-presenting cells. Additionally, the vaccine antigen is readily dissociated from the DODA-based liposomes leading to a population of antigen-presenting cells that are antigen-positive but adjuvant-negative and consequently are not activated. These studies demonstrate the importance of studying in vivo characteristics of the vaccine components and furthermore show that physicochemical properties of the delivery system have a major impact on the vaccine-induced immune response. © 2012 Elsevier B.V. All rights reserved.

Structure-hepatic disposition relationships for cationic drugs in isolated perfused rat livers: Transmembrane exchange and cytoplasmic binding process

This work studied the structure-hepatic disposition relationships for cationic drugs of varying lipophilicity using a single-pass, in situ rat liver preparation. The lipophilicity among the cationic drugs studied in this work is in the following order: diltiazem. propranolol. labetalol. prazosin. antipyrine. atenolol. Parameters characterizing the hepatic distribution and elimination kinetics of the drugs were estimated using the multiple indicator dilution method. The kinetic model used to describe drug transport (the two-phase stochastic model) integrated cytoplasmic binding kinetics and belongs to the class of barrier-limited and space-distributed liver models. Hepatic extraction ratio (E) (0.30-0.92) increased with lipophilicity. The intracellular binding rate constant (k(on)) and the equilibrium amount ratios characterizing the slowly and rapidly equilibrating binding sites (K-S and K-R) increase with the lipophilicity of drug (k(on) : 0.05-0.35 s(-1); K-S : 0.61-16.67; K-R : 0.36-0.95), whereas the intracellular unbinding rate constant (k(off)) decreases with the lipophilicity of drug (0.081-0.021 s(-1)). The partition ratio of influx (k(in)) and efflux rate constant (k(out)), k(in)/k(out), increases with increasing pK(a) value of the drug [from 1.72 for antipyrine (pK(a) = 1.45) to 9.76 for propranolol (pK(a) = 9.45)], the differences in k(in/kout) for the different drugs mainly arising from ion trapping in the mitochondria and lysosomes. The value of intrinsic elimination clearance (CLint), permeation clearance (CLpT), and permeability-surface area product (PS) all increase with the lipophilicity of drug [CLint (ml . min(-1) . g(-1) of liver): 10.08-67.41; CLpT (ml . min(-1) . g(-1) of liver): 10.80-5.35; PS (ml . min(-1) . g(-1) of liver): 14.59-90.54]. It is concluded that cationic drug kinetics in the liver can be modeled using models that integrate the presence of cytoplasmic binding, a hepatocyte barrier, and a vascular transit density function.

Cationic drug pharmacokinetics in diseased livers determined by fibrosis index, hepatic protein content, microsomal activity, and nature of drug

The disposition kinetics of six cationic drugs in perfused diseased and normal rat livers were determined by multiple indicator dilution and related to the drug physicochemical properties and liver histopathology. A carbon tetrachloride (CCl4)induced acute hepatocellular injury model had a higher fibrosis index (FI), determined by computer-assisted image analysis, than did an alcohol-induced chronic hepatocellular injury model. The alcohol-treated group had the highest hepatic alpha(1)- acid glycoprotein, microsomal protein (MP), and cytochrome P450 (P450) concentrations. Various pharmacokinetic parameters could be related to the octanol-water partition coefficient (log P-app) of the drug as a surrogate for plasma membrane partition coefficient and affinity for MP or P450, the dependence being lower in the CCl4-treated group and higher in the alcohol-treated group relative to controls. Stepwise regression analysis showed that hepatic extraction ratio, permeability-surface area product, tissue-binding constant, intrinsic clearance, partition ratio of influx (k(in)) and efflux rate constant (k(out)), and k(in)/k(out) were related to physicochemical properties of drug (log P-app or pK(a)) and liver histopathology (FI, MP, or P450). In addition, hepatocyte organelle ion trapping of cationic drugs was evident in all groups. It is concluded that fibrosis-inducing hepatic disease effects on cationic drug disposition in the liver may be predicted from drug properties and liver histopathology.

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.

Catheter effects in organ perfusion experiments

In a typical isolated organ perfusion experiment, a substance is injected upstream of an organ and then collected at some distance downstream. To reach the organ from the injection site, and then from the organ to the collector, a solute passes through catheters, usually tubes with circular cross-sections. Catheters cause distortion to the concentration-time profile of the perfusion. In this paper, we analyse catheter distribution kinetics from a mathematical point of view, develop the function most suitable for modeling this distribution and successfully apply this function to experimental data. (C) 2002 Academic Press.

Kinetic analysis of vascular marker distribution in perfused rat livers after regeneration following partial hepatectomy

Background/Aims: Liver clearance models are based on information (or assumptions) on solute distribution kinetics within the microvasculatory system, The aim was to study albumin distribution kinetics in regenerated livers and in livers of normal adult rats, Methods: A novel mathematical model was used to evaluate the distribution space and the transit time dispersion of albumin in livers following regeneration after a two-thirds hepatectomy compared to livers of normal adult rats. Outflow curves of albumin measured after bolus injection in single-pass perfused rat livers were analyzed by correcting for the influence of catheters and fitting a long-tailed function to the data. Results: The curves were well described by the proposed model. The distribution volume and the transit time dispersion of albumin observed in the partial hepatectomy group were not significantly different from livers of normal adult rats. Conclusions: These findings suggest that the distribution space and the transit time dispersion of albumin (CV2) is relatively constant irrespective of the presence of rapid and extensive repair. This invariance of CV2 implies, as a first approximation, a similar degree of intrasinusoidal mixing, The finding that a sum of two (instead of one) inverse Gaussian densities is an appropriate empirical function to describe the outflow curve of vascular indicators has consequences for an improved prediction of hepatic solute extraction.

Drug targeting using OX7-immunoliposomes: correlation between Thy1.1 antigen expression and tissue distribution in the rat

OX7 monoclonal antibody F((ab')2) fragments directed against Thy1.1 antigen can be used for drug targeting by coupling to the surface of drug-loaded liposomes. Such OX7-conjugated immunoliposomes (OX7-IL) were used recently for drug delivery to rat glomerular mesangial cells, which are characterized by a high level of Thy1.1 antigen expression. In the present study, the relationship between OX7-IL tissue distribution and target Thy1.1 antigen localization in different organs in rat was investigated. Western blot and immunohistofluorescence analysis revealed a very high Thy1.1 expression in brain cortex and striatum, thymus and renal glomeruli. Moderate Thy1.1 levels were observed in the collecting ducts of kidney, lung tissue and spleen. Thy1.1 was not detected in liver and heart. There was a poor correlation between Thy1.1 expression levels and organ distribution of fluorescence- or (14)C-labeled OX7-IL. The highest overall organ density of OX7-IL was observed in the spleen, followed by lung, liver and kidney. Heart and brain remained negative. With respect to intra-organ distribution, a localized and distinct signal was observed in renal glomerular mesangial cells only. As a consequence, acute pharmacological (i.e. toxic) effects of doxorubicin-loaded OX7-IL were limited to renal glomeruli. The competition with unbound OX7 monoclonal antibody F((ab')2) fragments demonstrated that the observed tissue distribution and acute pharmacological effects of OX7-IL were mediated specifically by the conjugated OX7 antibody. It is concluded that both the high target antigen density and the absence of endothelial barriers are needed to allow for tissue-specific accumulation and pharmacological effects of OX7-IL. The liposomal drug delivery strategy used is therefore specific toward renal glomeruli and can be expected to reduce the risk of unwanted side effects in other tissues.

Metabolite mean transit times in the liver as predicted by various models of hepatic elimination

Predicted area under curve (AUC), mean transit time (MTT) and normalized variance (CV2) data have been compared for parent compound and generated metabolite following an impulse input into the liver, Models studied were the well-stirred (tank) model, tube model, a distributed tube model, dispersion model (Danckwerts and mixed boundary conditions) and tanks-in-series model. It is well known that discrimination between models for a parent solute is greatest when the parent solute is highly extracted by the liver. With the metabolite, greatest model differences for MTT and CV2 occur when parent solute is poorly extracted. In all cases the predictions of the distributed tube, dispersion, and tasks-in-series models are between the predictions of the rank and tube models. The dispersion model with mixed boundary conditions yields identical predictions to those for the distributed tube model (assuming an inverse gaussian distribution of tube transit times). The dispersion model with Danckwerts boundary conditions and the tanks-in series models give similar predictions to the dispersion (mixed boundary conditions) and the distributed tube. The normalized variance for parent compound is dependent upon hepatocyte permeability only within a distinct range of permeability values. This range is similar for each model but the order of magnitude predicted for normalized variance is model dependent. Only for a one-compartment system is the MIT for generated metabolite equal to the sum of MTTs for the parent compound and preformed metabolite administered as parent.