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.

Hepatic disposition and metabolite kinetics of a homologous series of diflunisal esters

The hepatic disposition and metabolite kinetics of a homologous series of diflunisal O-acyl esters (acetyl, butanoyl, pentanoyl, anti hexanoyl) were determined using a single-pass perfused in situ rat liver preparation. The experiments were conducted using 2% BSA Krebs-Henseleit buffer (pH 7.4), and perfusions were performed at 30 mL/min in each liver. O-Acyl esters of diflunisal and pregenerated diflunisal were injected separately into the portal vein. The venous outflow samples containing the esters and metabolite diflunisal were analyzed by high performance liquid chromatography (HPLC). The normalized outflow concentration-time profiles for each parent ester and the formed metabolite, diflunisal, were analyzed using statistical moments analysis and the two-compartment dispersion model. Data (presented as mean +/- standard error for triplicate experiments) was compared using ANOVA repeated measures, significance level P < 0.05. The hepatic availability (AUC'), the fraction of the injected dose recovered in the outflowing perfusate, for O-acetyldiflunisal (C2D = 0.21 +/- 0.03) was significantly lower than the other esters (0.34-0.38). However, R-N/f(u), the removal efficiency number R-N divided by the unbound fraction in perfusate f(u), which represents the removal efficiency of unbound ester by the liver, was significantly higher for the most lipophilic ester (O-hexanoyldiflunisal, C6D = 16.50 +/- 0.22) compared to the other members of the series (9.57 to 11.17). The most lipophilic ester, C6D, had the largest permeability surface area (PS) product (94.52 +/- 38.20 mt min-l g-l liver) and tissue distribution value VT (35.62 +/- 11.33 mL g(-1) liver) in this series. The MTT of these O-acyl esters of diflunisal were not significantly different from one another. However, the metabolite diflunisal MTTs tended to increase with the increase in the parent ester lipophilicity (11.41 +/- 2.19 s for C2D to 38.63 +/- 9.81 s for C6D). The two-compartment dispersion model equations adequately described the outflow profiles for the parent esters and the metabolite diflunisal formed from the O-acyl esters of diflunisal in the liver.

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.

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.

Physiologically based synthetic models of hepatic disposition

Current Physiologically based pharmacokinetic (PBPK) models are inductive. We present an additional, different approach that is based on the synthetic rather than the inductive approach to modeling and simulation. It relies on object-oriented programming A model of the referent system in its experimental context is synthesized by assembling objects that represent components such as molecules, cells, aspects of tissue architecture, catheters, etc. The single pass perfused rat liver has been well described in evaluating hepatic drug pharmacokinetics (PK) and is the system on which we focus. In silico experiments begin with administration of objects representing actual compounds. Data are collected in a manner analogous to that in the referent PK experiments. The synthetic modeling method allows for recognition and representation of discrete event and discrete time processes, as well as heterogeneity in organization, function, and spatial effects. An application is developed for sucrose and antipyrine, administered separately and together PBPK modeling has made extensive progress in characterizing abstracted PK properties but this has also been its limitation. Now, other important questions and possible extensions emerge. How are these PK properties and the observed behaviors generated? The inherent heuristic limitations of traditional models have hindered getting meaningful, detailed answers to such questions. Synthetic models of the type described here are specifically intended to help answer such questions. Analogous to wet-lab experimental models, they retain their applicability even when broken apart into sub-components. Having and applying this new class of models along with traditional PK modeling methods is expected to increase the productivity of pharmaceutical research at all levels that make use of modeling and simulation.

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.

Fatty acid binding protein is a major determinant of hepatic pharmacokinetics of palmitate and its metabolites

Disposition kinetics of [H-3] palmitate and its low-molecular-weight metabolites in perfused rat livers were studied using the multiple-indicator dilution technique, a selective assay for [H-3] palmitate and its low-molecular-weight metabolites, and several physiologically based pharmacokinetic models. The level of liver fatty acid binding protein (L-FABP), other intrahepatic binding proteins (microsomal protein, albumin, and glutathione S-transferase) and the outflow profiles of [H-3] palmitate and metabolites were measured in four experimentalgroups of rats: 1) males; 2) clofibrate-treated males; 3) females; and 4) pregnant females. A slow-diffusion/bound model was found to better describe the hepatic disposition of unchanged [H-3] palmitate than other pharmacokinetic models. The L-FABP levels followed the order: pregnant female > clofibrate-treated male > female > male. Levels of other intrahepatic proteins did not differ significantly. The hepatic extraction ratio and mean transit time for unchanged palmitate, as well as the production of low-molecular-weight metabolites of palmitate and their retention in the liver, increased with increasing L-FABP levels. Palmitate metabolic clearance, permeability-surface area product, retention of palmitate by the liver, and cytoplasmic diffusion constant for unchanged [H-3] palmitate also increased with increasing L-FABP levels. It is concluded that the variability in hepatic pharmacokinetics of unchanged [H-3] palmitate and its low-molecular-weight metabolites in perfused rat livers is related to levels of L-FABP and not those of other intrahepatic proteins.

Disposition kinetics of propranolol isomers in the perfused rat liver

The aim of this study was to define the determinants of the linear hepatic disposition kinetics of propranolol optical isomers using a perfused rat liver. Monensin was used to abolish the lysosomal proton gradient to allow an estimation of propranolol ion trapping by hepatic acidic vesicles. In vitro studies were used for independent estimates of microsomal binding and intrinsic clearance. Hepatic extraction and mean transit time were determined from outflow-concentration profiles using a nonparametric method. Kinetic parameters were derived from a physiologically based pharmacokinetic model. Modeling showed an approximate 34-fold decrease in ion trapping following monensin treatment. The observed model-derived ion trapping was similar to estimated theoretical values. No differences in ion-trapping values was found between R(+)- and S(-)- propranolol. Hepatic propranolol extraction was sensitive to changes in liver perfusate flow, permeability-surface area product, and intrinsic clearance. Ion trapping, microsomal and nonspecific binding, and distribution of unbound propranolol accounted for 47.4, 47.1, and 5.5% of the sequestration of propranolol in the liver, respectively. It is concluded that the physiologically more active S(-)- propranolol differs from the R(+)- isomer in higher permeability-surface area product, intrinsic clearance, and intracellular binding site values.

Reduced hepatic extraction of palmitate in steatosis correlated to lower level of liver fatty acid binding protein

Nonalcoholic fatty liver disease is the most common of all liver diseases. The hepatic disposition [H-3]palmitate and its low-molecular-weight metabolites in perfused normal and steatotic rat liver were studied using the multiple indicator dilution technique and a physiologically based slow diffusion/bound pharmacokinetic model. The steatotic rat model was established by administration of 17alpha-ethynylestradiol to female Wistar rats. Serum biochemistry markers and histology of treated and normal animals were assessed and indicated the presence of steatosis in the treatment group. The steatotic group showed a significantly higher alanine aminotransferase-to-aspartate aminotransferase ratio, lower levels of liver fatty acid binding protein and cytochrome P-450, as well as microvesicular steatosis with an enlargement of sinusoidal space. Hepatic extraction for unchanged [H-3]palmitate and production of low-molecular-weight metabolites were found to be significantly decreased in steatotic animals. Pharmacokinetic analysis suggested that the reduced extraction and sequestration for palmitate and its metabolites was mainly attributed to a reduction in liver fatty acid binding protein in steatosis.

Hepatic pharmacokinetics of propranolol in rats with adjuvant-induced systemic inflammation

Systemic inflammation is known to affect drug disposition in the liver. This study sought to relate and quantitate changes in hepatic pharmacokinetics of propranolol with changes in hepatic architecture and physiology in adjuvant-treated rats. Transmission electron microscopy was used to assess morphological changes in mitochondria and lysosomes of adjuvant-treated rat livers. The disposition of propranolol was assessed in the perfused rat liver using the multiple indicator dilution technique. Hepatic extraction and mean transit time were determined from outflow-concentration profiles using a nonparametric method. Kinetic parameters were derived from a two-phase physiologically based organ pharmacokinetic model. Possible relationships were then explored between the changes in hepatic drug disposition and cytochrome P-450 activity and iron concentration. Adjuvant treatment induced the appearance of mitochondrial inclusions/tubularization and irregularly shaped lysosomes in rat livers. Livers from adjuvant-treated rats had (relative to normal) significantly higher alpha(1)-acid glycoprotein (orosomucoid) and iron tissue concentrations but lower cytochrome P-450 content. The hepatic extraction, metabolism, and ion trapping of propranolol were significantly impaired in adjuvant-treated rats and could be correlated with altered iron store and cytochrome P-450 activity. It is concluded that adjuvant-induced systemic inflammation alters hepatocellular morphology and biochemistry and consequently influences hepatic disposition of propranolol.

Disposition of venlafaxine enantiomers in rats with hepatic encephalopathy after chronic drug treatment.

Portacaval shunted (PCS) rats, a model of hepatic encephalopathy, and control animals were administered racemic venlafaxine for 14 days (10 mg/kg). The levels of the S- and R-enantiomers and the S/R-enantiomer ratios of venlafaxine and its metabolites were assessed by an enantiomer-selective chromatographic assay in serum, brain parenchyma, and brain dialysate of both groups. Higher levels of the S- and R-enantiomers of venlafaxine were found in serum and brain of PCS vs. normal rats (median values of S- and R-venlafaxine in serum: 290 and 201 nM in PCS; 97 and 66 nM in normal rats; median values of S- and R-venlafaxine in cortex: 956 and 939 nM in PCS; 357 and 318 nM in normal rats). Interestingly, similar S/R-venlafaxine ratios were observed in PCS and normal rats both in serum (S/R = 1.4) and brain compartments (S/R = l.0-1.1). These findings may have clinical relevance for the safety of venlafaxine in chronic hepatic encephalopathy.

Reduction of enantioselectivity in the kinetic disposition and metabolism of verapamil in rats exposed to toluene

Toluene and verapamil are subject to extensive oxidative metabolism mediated by CYP enzymes, and their interaction can be stereoselective. In the present study we investigated the influence of toluene inhalation on the enantioselective kinetic disposition of verapamil and its metabolite, norverapamil, in rats. Male Wistar rats (n = 6 per group) received a single dose of racemic verapamil (10 mg/kg) orally at the fifth day of nose-only toluene or air (control group) inhalation for 6 h/day (25, 50, and 100 ppm). Serial blood samples were collected from the tail up to 6 h after verapamil administration. The plasma concentrations of verapamil and norverapamil enantiomers were analyzed by LC-MS/MS by using a Chiralpak AD column. Toluene inhalation did not influence the kinetic disposition of verapamil or norverapamil enantiomers (p > 0.05, Kruskal-Wallis test) in rats. The pharmacokinetics of verapamil was enantioselective in the control group, with a higher plasma proportion of the S-verapamil (AUC 250.8 versus 120.4 ng.h.mL(-1); p <= 0.05, Wilcoxon test) and S-norverapamil (AUC 72.3 versus 52.3 ng.h.mL(-1); p <= 0.05, Wilcoxon test). Nose-only exposure to toluene at 25, 50, or 100 ppm resulted in a lack of enantioselectivity for both verapamil and norverapamil. The study demonstrates the importance of the application of enantioselective methods in studies on the interaction between solvents and chiral drugs.

Polyinosinic acid and polycationic liposomes attenuate the hepatic clearance of circulating plasmid DNA

DNA that enters the circulation is rapidly cleared both by tissue uptake and by DNase-mediated degradation. In this study, we have examined the uptake of linear plasmid DNA in an isolated perfused liver model and following intra-arterial administration to rats. We found that the DNA was rapidly taken up by the isolated perfused liver without degradation. The single-pass extraction ratio was 0.76 +/- 0.05, the mean transit time was 15.3 +/- 3.6 s, and the volume of distribution was 0.29 +/- 0.07 ml/g. Hepatic uptake was saturable and was inhibited by polyinosinic acid or polycationic liposomes but not by condensation of the DNA with polylysine. When the linear plasmid DNA was administered in vivo, plasma half-life was 3.1 +/- 0.2 min, volume of distribution was 670 +/- 85 ml/kg, and clearance was 32 +/- 4 min. Coadministration of cationic liposomes decreased the volume of distribution to 180 +/- 28 ml/kg as well as the half-life (2.6 +/- 0.2 min). By contrast, polyinosinic acid significantly increased the circulating half-life (7.7 +/- 0.5 min), decreased the volume of distribution (95 +/- 17 ml/kg), and partially inhibited DNA degradation. When administered along with the liposomes and the polyinosinic acid, the distribution of plasmid-derived radioactivity decreased in the liver and increased in most other peripheral tissues. This study shows that pharmacological manipulation of the uptake and degradation of DNA can alter its distribution and clearance in vivo. These results may be useful in optimizing gene delivery procedures for in vivo gene therapy.