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.

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.

On the validity of the dispersion model of hepatic drug elimination when intravascular transit time densities are long-tailed

The dispersion model with mixed boundary conditions uses a single parameter, the dispersion number, to describe the hepatic elimination of xenobiotics and endogenous substances. An implicit a priori assumption of the model is that the transit time density of intravascular indicators is approximated by an inverse Gaussian distribution. This approximation is limited in that the model poorly describes the tail part of the hepatic outflow curves of vascular indicators. A sum of two inverse Gaussian functions is proposed as ail alternative, more flexible empirical model for transit time densities of vascular references. This model suggests that a more accurate description of the tail portion of vascular reference curves yields an elimination rate constant (or intrinsic clearance) which is 40% less than predicted by the dispersion model with mixed boundary conditions. The results emphasize the need to accurately describe outflow curves in using them as a basis for determining pharmacokinetic parameters using hepatic elimination models. (C) 1997 Society for Mathematical Biology.

Relative dispersions of intra-albumin transit times across rat and elasmobranch perfused livers, and implications for intra- and inter-species scaling of hepatic clearance using microsomal data

It is recognized that vascular dispersion in the liver is a determinant of high first-pass extraction of solutes by that organ. Such dispersion is also required for translation of in-vitro microsomal activity into in-vivo predictions of hepatic extraction for any solute. We therefore investigated the relative dispersion of albumin transit times (CV2) in the livers of adult and weanling rats and in elasmobranch livers. The mean and normalized variance of the hepatic transit time distribution of albumin was estimated using parametric non-linear regression (with a correction for catheter influence) after an impulse (bolus) input of labelled albumin into a single-pass liver perfusion. The mean +/- s.e. of CV2 for albumin determined in each of the liver groups were 0.85 +/- 0.20 (n = 12), 1.48 +/- 0.33 (n = 7) and 0.90 +/- 0.18 (n = 4) for the livers of adult and weanling rats and elasmobranch livers, respectively. These CV2 are comparable with that reported previously for the dog and suggest that the CV2 Of the liver is of a similar order of magnitude irrespective of the age and morphological development of the species. It might, therefore, be justified, in the absence of other information, to predict the hepatic clearances and availabilities of highly extracted solutes by scaling within and between species livers using hepatic elimination models such as the dispersion model with a CV2 of approximately unity.

Interconnected-tubes model of hepatic elimination

The distributed-tubes model of hepatic elimination is extended to include intermixing between sinusoids, resulting in the formulation of a new, interconnected-tubes model. The new model is analysed for the simple case of two interconnected tubes, where an exact solution is obtained. For the case of many strongly-interconnected tubes, it is shown that a zeroth-order approximation leads to the convection-dispersion model. As a consequence the dispersion number is expressed, for the first time, in terms of its main physiological determinants: heterogeneity of flow and density of interconnections between sinusoids. The analysis of multiple indicator dilution data from a perfused liver preparation using the simplest version of the model yields the estimate 10.3 for the average number of interconnections. The problem of boundary conditions for the dispersion model is considered from the viewpoint that the dispersion-convection equation is a zeroth-order approximation to the equations for the interconnected-tubes model. (C) 1997 Academic Press Limited.

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.

Neonatal hepatic drug elimination

After the transition from in utero to newborn life, the neonate becomes solely reliant upon its own drug clearance processes to metabolise xenobiotics. Whilst most studies of neonatal hepatic drug elimination have focussed upon in vitro expression and activities of drug-metabolising enzymes, the rapid physiological changes in the early neonatal period of life also need to be considered. There are dramatic changes in neonatal liver blood how and hepatic oxygenation due to the loss of the umbilical blood supply, the increasing portal vein blood flow, and the gradual closure of the ductus venosus shunt during the first week of life. These changes which may well affect the capacity of neonatal hepatic drug metabolism. The hepatic expression of cytochromes P450 1A2, 2C, 2D6, 2E1 and 3A4 develop at different rates in the postnatal period, whilst 3A7 expression diminishes. Hepatic glucuronidation in the human neonate is relatively immature at birth, which contrasts with the considerably more mature neonatal hepatic sulfation activity. Limited in vivo studies show that the human neonate can significantly metabolise xenobiotics but clearance is considerably less compared with the older infant and adult. The neonatal population included in pharmacological studies is highly heterogeneous with respect to age, body weight, ductus venosus closure and disease processes, making it difficult to interpret data arising from human neonatal studies. Studies in the perfused foetal and neonatal sheep liver have demonstrated how the oxidative and conjugative hepatic elimination of drugs by the intact organ is significantly increased during the first week of life, highlighting that future studies will need to consider the profound physiological changes that may influence neonatal hepatic drug elimination shortly after birth.

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.

Development of 3D in vitro models for prediction of hepatic metabolism and toxicity

The development of human cell models that recapitulate hepatic functionality allows the study of metabolic pathways involved in toxicity and disease. The increased biological relevance, cost-effectiveness and high-throughput of cell models can contribute to increase the efficiency of drug development in the pharmaceutical industry. Recapitulation of liver functionality in vitro requires the development of advanced culture strategies to mimic in vivo complexity, such as 3D culture, co-cultures or biomaterials. However, complex 3D models are typically associated with poor robustness, limited scalability and compatibility with screening methods. In this work, several strategies were used to develop highly functional and reproducible spheroid-based in vitro models of human hepatocytes and HepaRG cells using stirred culture systems. In chapter 2, the isolation of human hepatocytes from resected liver tissue was implemented and a liver tissue perfusion method was optimized towards the improvement of hepatocyte isolation and aggregation efficiency, resulting in an isolation protocol compatible with 3D culture. In chapter 3, human hepatocytes were co-cultivated with mesenchymal stem cells (MSC) and the phenotype of both cell types was characterized, showing that MSC acquire a supportive stromal function and hepatocytes retain differentiated hepatic functions, stability of drug metabolism enzymes and higher viability in co-cultures. In chapter 4, a 3D alginate microencapsulation strategy for the differentiation of HepaRG cells was evaluated and compared with the standard 2D DMSO-dependent differentiation, yielding higher differentiation efficiency, comparable levels of drug metabolism activity and significantly improved biosynthetic activity. The work developed in this thesis provides novel strategies for 3D culture of human hepatic cell models, which are reproducible, scalable and compatible with screening platforms. The phenotypic and functional characterization of the in vitro systems performed contributes to the state of the art of human hepatic cell models and can be applied to the improvement of pre-clinical drug development efficiency of the process, model disease and ultimately, development of cell-based therapeutic strategies for liver failure.

Hepatic structure-pharmacokinetic relationships: The hepatic disposition and metabolite kinetics of a homologous series of O-acyl derivatives of salicylic acid

1 The hepatic disposition and metabolite kinetics of a homologous series of O-acyl (acetyl, propionyl, butanoyl, pentanoyl, hexanoyl and octanoyl) esters of salicylic acid (C2SA, C3SA, C4SA, C5SA, C6SA and C8SA, respectively) was determined using a single-pass, in-sills rat liver preparation. 2 The hepatic venous outflow profiles for the parent esters and the generated metabolite, salicylic acid (SA) were analysed by HPLC. Non-parametric moments analysis was used to determine the area under the curve (AUC'), mean transit time (MTT) and normalized variance (CV2) for the parent esters and generated SA. 3 Pregenerated SA ([C-14]-salicylic acid) was injected into each liver with the parent ester to determine its distribution characteristics. 4 The overall recovery of ester plus metabolite was 89% of the ester dose injected and independent of the ester carbon number, suggesting that ester extraction was due to hepatic metabolism to salicylic acid. 5 The metabolite AUC' value increased directly with the lipophilicity of the parent ester (from 0.12 for C2SA to 0.95 for C8SA). By contrast, the parent AUC' decreased with the lipophilicity (from 0.85 for C2SA to zero for C8SA). The metabolite MTT value also showed a trend to increase with the lipophilicity of the parent ester (from 15.72 s for C3SA to 61.97 s for C8SA). However, the parent MTT value shows no significant change across the series. 6 The two-compartment dispersion model was used to derive the kinetic parameters for parent ester, pregenerated SA and generated SA. Consequently, these parameters were used to estimate the values of AUG', MITT and CV2 for the parent ester and metabolite. The moments values obtained using the two-compartment dispersion model show similar trends to the corresponding moments values obtained from the outflow profiles using a non-parametric approach. 7 The more lipophilic aspirin analogues are more confined to the portal circulation after oral administration than aspirin due to their more extensive hepatic elimination avoiding systemic prostacyclin inhibition. Given that aspirin's selectivity as an anti-thrombotic agent has been postulated to be due to selective anti-platelet effects in the portal circulation, the more lipophilic and highly extracted analogues are potentially more selective anti-thrombotic agents than aspirin.

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.

Experimental methods for studying drug uptake in the head and brain

A number of techniques have been developed to study the disposition of drugs in the head and, in particular, the role of the blood-brain barrier (BBB) in drug uptake. The techniques can be divided into three groups: in-vitro, in-vivo and in-situ. The most suitable method depends on the purpose(s) and requirements of the particular study being conducted. In-vitro techniques involve the isolation of cerebral endothelial cells so that direct investigations of these cells can be carried out. The most recent preparations are able to maintain structural and functional characteristics of the BBB by simultaneously culturing endothelial cells with astrocytic cells,The main advantages of the in-vitro methods are the elimination of anaesthetics and surgery. In-vivo methods consist of a diverse range of techniques and include the traditional Brain Uptake Index and indicator diffusion methods, as well as microdialysis and positron emission tomography. In-vivo methods maintain the cells and vasculature of an organ in their normal physiological states and anatomical position within the animal. However, the shortcomings include renal acid hepatic elimination of solutes as well as the inability to control blood flow. In-situ techniques, including the perfused head, are more technically demanding. However, these models have the ability to vary the composition and flow rate of the artificial perfusate. This review is intended as a guide for selecting the most appropriate method for studying drug uptake in the brain.

Relative dispersion of intravascular transit times during isolated human limb perfusions for recurrent melanoma

Aims We have characterized the relative dispersion of vascular and extravascular markers in the limbs of three patients undergoing isolated limb perfusions with the cytotoxic melphalan for recurrent malignant melanoma both before and after melphalan dosing. Methods A bolus of injectate containing [Cr-51] labelled red blood cells, [C-14]-sucrose and [H-3]-water was injected into an iliac or femoral artery and outflow samples collected at 1 s intervals by a fraction collector. The radioactivity due to each isotype was analysed by either gamma [Cr-51] or beta [C-14 and H-3] counting. The moments of the outflow fraction-time profiles were estimated by a nonparametric (numerical integration) method and a parametric model (sum of two inverse Gaussian functions). Results The availability, mean transit time and normalised variance (CV2) obtained for labelled red blood cells, sucrose and water were similar before and after melphalan dosing and with the two methods of calculation but varied between the patients. Conclusions The vascular space is not well-stirred but characterized by a CV2 similar that reported previously for in situ rat hind limb and rat liver perfusions. A flow-limited blood-tissue exchange was observed for the permeating indicators. Administration of melphalan did not influence the distribution characteristics of the indicators.

Synergistic acceleration of thyroid hormone degradation by phenobarbital and the PPAR alpha agonist WY14643 in rat hepatocytes

Energy balance is maintained by controlling both energy intake and energy expenditure. Thyroid hormones play a crucial role in regulating energy expenditure. Their levels are adjusted by a tight feed back-control led regulation of thyroid hormone production/incretion and by their hepatic metabolism. Thyroid hormone degradation has previously been shown to be enhanced by treatment with phenobarbital or other antiepileptic drugs due to a CAR-dependent induction of phase 11 enzymes of xenobiotic metabolism. We have recently shown, that PPAR alpha agonists synergize with phenobarbital to induce another prototypical CAR target gene, CYP2B1. Therefore, it was tested whether a PPAR alpha agonist could enhance the phenobarbital-dependent acceleration of thyroid hormone elimination. In primary cultures of rat hepatocytes the apparent half-life of T3 was reduced after induction with a combination of phenobarbital and the PPARa agonist WY14643 to a larger extent than after induction with either Compound alone. The synergistic reduction of the half-life could be attributed to a synergistic induction of CAR and the CAR target genes that code for enzymes and transporters involved in the hepatic elimination of T3, such as OATP1A1, OATP1A3, UGT1A3 and UCT1A10. The PPAR alpha-dependent CAR induction and the subsequent induction of T3-eliminating enzymes might be of physiological significance for the fasting-incluced reduction in energy expenditure by fatty acids as natural PPARa ligands. The synergism of the PPAR alpha agonist WY14643 and phenobarbital in inducing thyroid hormone breakdown might serve as a paradigm for the synergistic disruption of endocrine control by other combinations of xenobiotics. (C) 2009 Elsevier Inc. All rights reserved.

Selenite and Selenate Effects on Mercury (Hg2+) Uptake and Distribution in Tilapia, Oreochromis niloticus L., Assessed by Chronic Bioassay

Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP)