In vitro and in vivo assessment of hepatic and extrahepatic metabolism of diazepam in the rat

Utility of Minimal Physiologically Based Pharmacokinetic Models for Assessing Fractional Distribution, Oral Absorption, and Series-Compartment Models of Hepatic Clearance

Minimal physiologically based pharmacokinetic (mPBPK) models are physiologically relevant, require less information than full PBPK models, and offer flexibility in pharmacokinetics (PK). The well-stirred hepatic model (WSM) is commonly used in PBPK, whereas the more plausible dispersion model (DM) poses computational complexities. The series-compartment model (SCM) mimics the DM but is easier to operate. This work implements the SCM and mPBPK models for assessing fractional tissue distribution, oral absorption, and hepatic clearance using literature-reported blood and liver concentration-time data in rats for compounds mainly cleared by the liver. Further handled were various complexities, including nonlinear hepatic binding and metabolism, differing absorption kinetics, and sites of administration. The SCM containing one to five (n) liver subcompartments yields similar fittings and provides comparable estimates for hepatic extraction ratio (ER), prehepatic availability (Fg ), and first-order absorption rate constants (ka ). However, they produce decreased intrinsic clearances (CLint ) and liver-to-plasma partition coefficients (Kph ) with increasing n as expected. Model simulations demonstrated changes in intravenous and oral PK profiles with alterations in Kph and ka and with hepatic metabolic zonation. The permeability (PAMPA P) of the various compounds well explained the fitted fractional distribution (fd ) parameters. The SCM and mPBPK models offer advantages in distinguishing systemic, extrahepatic, and hepatic clearances. The SCM allows for incorporation of liver zonation and is useful in assessing changes in internal concentration gradients potentially masked by similar blood PK profiles. Improved assessment of intraorgan drug concentrations may offer insights into active moieties driving metabolism, biliary excretion, pharmacodynamics, and hepatic toxicity. SIGNIFICANCE STATEMENT: The minimal physiologically based pharmacokinetic model and the series-compartment model are useful in assessing oral absorption and hepatic clearance. They add flexibility in accounting for various drug- or system-specific complexities, including fractional distribution, nonlinear binding and saturable hepatic metabolism, and hepatic zonation. These models can offer improved insights into the intraorgan concentrations that reflect physiologically active moieties often driving disposition, pharmacodynamics, and toxicity.

Successful and Unsuccessful Prediction of Human Hepatic Clearance for Lead Optimization

Development of new chemical entities is costly, time-consuming, and has a low success rate. Accurate prediction of pharmacokinetic properties is critical to progress compounds with favorable drug-like characteristics in lead optimization. Of particular importance is the prediction of hepatic clearance, which determines drug exposure and contributes to projection of dose, half-life, and bioavailability. The most commonly employed methodology to predict hepatic clearance is termed in vitro to in vivo extrapolation (IVIVE) that involves measuring drug metabolism in vitro, scaling-up this in vitro intrinsic clearance to a prediction of in vivo intrinsic clearance by reconciling the enzymatic content between the incubation and an average human liver, and applying a model of hepatic disposition to account for limitations of protein binding and blood flow to predict in vivo clearance. This manuscript reviews common in vitro techniques used to predict hepatic clearance as well as current challenges and recent theoretical advancements in IVIVE.

Intranasal Coadministration of a Diazepam Prodrug with a Converting Enzyme Results in Rapid Absorption of Diazepam in Rats

Intranasal administration is an attractive route for systemic delivery of small, lipophilic drugs because they are rapidly absorbed through the nasal mucosa into systemic circulation. However, the low solubility of lipophilic drugs often precludes aqueous nasal spray formulations. A unique approach to circumvent solubility issues involves coadministration of a hydrophilic prodrug with an exogenous converting enzyme. This strategy not only addresses poor solubility but also leads to an increase in the chemical activity gradient driving drug absorption. Herein, we report plasma and brain concentrations in rats following coadministration of a hydrophilic diazepam prodrug, avizafone, with the converting enzyme human aminopeptidase B Single doses of avizafone equivalent to diazepam at 0.500, 1.00, and 1.50 mg/kg were administered intranasally, resulting in 77.8% ± 6.0%, 112% ± 10%, and 114% ± 7% bioavailability; maximum plasma concentrations 71.5 ± 9.3, 388 ± 31, and 355 ± 187 ng/ml; and times to peak plasma concentration 5, 8, and 5 minutes for each dose level, respectively. Both diazepam and a transient intermediate were absorbed. Enzyme kinetics incorporated into a physiologically based pharmacokinetic model enabled estimation of the first-order absorption rate constants: 0.0689 ± 0.0080 minutes^-1 for diazepam and 0.122 ± 0.022 minutes^-1 for the intermediate. Our results demonstrate that diazepam, which is practically insoluble, can be delivered intranasally with rapid and complete absorption by coadministering avizafone with aminopeptidase B. Furthermore, even faster rates of absorption might be attained simply by increasing the enzyme concentration, potentially supplanting intravenous diazepam or lorazepam or intramuscular midazolam in the treatment of seizure emergencies.

The Constraints, Construction, and Verification of a Strain-Specific Physiologically Based Pharmacokinetic Rat Model

The use of in vitro-in vivo extrapolation (IVIVE) techniques, mechanistically incorporated within physiologically based pharmacokinetic (PBPK) models, can harness in vitro drug data and enhance understanding of in vivo pharmacokinetics. This study's objective was to develop a user-friendly rat (250 g, male Sprague-Dawley) IVIVE-linked PBPK model. A 13-compartment PBPK model including mechanistic absorption models was developed, with required system data (anatomical, physiological, and relevant IVIVE scaling factors) collated from literature and analyzed. Overall, 178 system parameter values for the model are provided. This study also highlights gaps in available system data required for strain-specific rat PBPK model development. The model's functionality and performance were assessed using previous literature-sourced in vitro properties for diazepam, metoprolol, and midazolam. The results of simulations were compared against observed pharmacokinetic rat data. Predicted and observed concentration profiles in 10 tissues for diazepam after a single intravenous (i.v.) dose making use of either observed i.v. clearance (CL_iv) or in vitro hepatocyte intrinsic clearance (CL_int) for simulations generally led to good predictions in various tissue compartments. Overall, all i.v. plasma concentration profiles were successfully predicted. However, there were challenges in predicting oral plasma concentration profiles for metoprolol and midazolam, and the potential reasons and according solutions are discussed.

Prediction of human pharmacokinetics—evaluation of methods for prediction of hepatic metabolic clearance

Methods for prediction of hepatic clearance (CLH) in man have been evaluated. A physiologically-based in-vitro to in-vivo (PB-IVIV) method with human unbound fraction in blood (fu,bl) and hepatocyte intrinsic clearance (CLint)-data has a good rationale and appears to give the best predictions (maximum ∼2-fold errors; < 25% errors for half of CL-predictions; appropriate ranking). Inclusion of an empirical scaling factor is, however, needed, and reasons include the use of cryopreserved hepatocytes with low activity, and inappropriate CLint- and fu,bl-estimation methods. Thus, an improvement of this methodology is possible and required. Neglect of fu,bl or incorporation of incubation binding does not seem appropriate. When microsome CLint-data are used with this approach, the CLH is underpredicted by 5- to 9-fold on average, and a 106-fold underprediction (attrition potential) has been observed. The poor performance could probably be related to permeation, binding and low metabolic activity. Inclusion of scaling factors and neglect of fu,bl for basic and neutral compounds improve microsome predictions. The performance is, however, still not satisfactory. Allometry incorrectly assumes that the determinants for CLH relate to body weight and overpredicts human liver blood flow rate. Consequently, allometric methods have poor predictability. Simple allometry has an average overprediction potential, > 2-fold errors for ∼1/3 of predictions, and 140-fold underprediction to 5800-fold overprediction (potential safety risk) range. In-silico methodologies are available, but these need further development. Acceptable prediction errors for compounds with low and high CLH should be ∼50 and ∼10%, respectively. In conclusion, it is recommended that PB-IVIV with human hepatocyte CLint and fu,bl is applied and improved, limits for acceptable errors are decreased, and that animal CLH-studies and allometry are avoided.

Prediction of human pharmacokinetics — renal metabolic and excretion clearance

The kidneys have the capability to both excrete and metabolise drugs. An understanding of mechanisms that determine these processes is required for the prediction of pharmacokinetics, exposures, doses and interactions of candidate drugs. This is particularly important for compounds predicted to have low or negligible non-renal clearance (CL). Clinically significant interactions in drug transport occur mostly in the kidneys. The main objective was to evaluate methods for prediction of excretion and metabolic renal CL (CLR) in humans. CLR is difficult to predict because of the involvement of bi-directional passive and active tubular transport, differences in uptake capacity, pH and residence time on luminal and blood sides of tubular cells, and limited knowledge about regional tubular residence time, permeability (Pe) and metabolic capacity. Allometry provides poor predictions of excretion CLR because of species differences in unbound fraction, urine pH and active transport. The correlation between fraction excreted unchanged in urine (fe) in humans and animals is also poor, except for compounds with high passive Pe (extensive/complete tubular reabsorption; zero/negligible fe) and/or high non-renal CL. Physiologically based in-vitro/in-vivo methods could potentially be useful for predicting CLR. Filtration could easily be predicted. Prediction of tubular secretion CL requires an in-vitro transport model and establishment of an in-vitro/in-vivo relationship, and does not appear to have been attempted. The relationship between passive Pe and tubular fraction reabsorbed (freabs) for compounds with and without apparent secretion has recently been established and useful equations and limits for prediction were developed. The suggestion that reabsorption has a lipophilicity cut-off does not seem to hold. Instead, compounds with passive Pe that is less than or equal to that of atenolol are expected to have negligible passive freabs. Compounds with passive Pe that is equal to or higher than that of carbamazepine are expected to have complete freabs. For compounds with intermediate Pe the relationship is irregular and freabs is difficult to predict. Tubular cells are comparably impermeable (for passive diffusion), and show regional differences in enzymatic and transporter activities. This limits the usefulness of microsome data and makes microsome-based predictions of metabolic CLR questionable. Renal concentrations and activities of CYP450s are comparably low, suggesting that CYP450 substrates have negligible metabolic CLR. The metabolic CLR of high-Pe UDP-glucuronyltransferase substrates could contribute to the total CL.

Glucuronidation of DRF-6574, hydroxy metabolite of DRF-4367 (a novel COX-2 inhibitor) by pooled human liver, intestinal microsomes and recombinant human UDP-glucuronosyltransferases (UGT): Role of UGT1A1,1A3 and 1A8

DRF-4367 is a novel COX-2 inhibitor, which showed good efficacy in several animal models of inflammation. In a comparative in vitro metabolism in various liver microsomes, DRF-4367 forms a hydroxy metabolite (DRF-6574) mediated by CYP2D6 and 2C19 isoenzymes. DRF-6574 readily undergoes Phase-II metabolism and forms glucuronide and sulfate conjugates both in vitro and in vivo. The objective of the present study was two folds: to study the glucuronidation of DRF-6574 in human liver and intestinal microsomes and to identify the recombinant human liver and intestinal UDP-glucuronosyltransferase (UGT) enzymes responsible for glucuronidation of DRF-6574. Of twelve recombinant UGTs tested, two hepatic UGTs viz., UGT1A1 and 1A3 and an extra hepatic UGT i.e., UGT1A8 showed the catalytic activity. The enzyme kinetics in pooled human liver, intestinal and recombinant UGT microsomes showed a typical Michaelis-Menten plot. The apparent Km and Vmax value for DRF-6574 was found to be 116 +/- 24 microM and 2.07 +/- 0.12 microg/min/mg protein and 142 +/- 17 microM and 3.83 +/- 0.15 microg/min/mg protein in pooled human liver and intestinal microsomes, respectively. The intrinsic clearance (Vmax/Km) value for DRF-6574 was estimated to be 0.043 and 0.065 ml/min/mg protein, respectively in pooled human liver and intestinal microsomes. Moreover we have determined the Km and Vmax and intrinsic clearance values for specific UGTs viz., UGT 1A1, 1A3 and 1A8. The apparent Km and Vmax values are 23 +/- 7.2 microM, 3.44 +/- 0.17 microg/min/mg protein for UGT1A1, 60 +/- 7.9 microM, 3.67 +/- 0.11 microg/min/mg protein for UGT1A3, 96 +/- 8.0 microM, 2.95 +/- 0.06 microg/min/mg protein for UGT1A8. The intrinsic clearance values (Vmax/Km) estimated were 0.367, 0.148, 0.074 ml/min/mg protein for UGT1A1, 1A3 and 1A8, respectively. The intrinsic clearance value in UGT1A8 was very close to that in human intestinal and liver microsomes. The formation of DRF-6574 glucuronide by human liver, intestinal and UGT1A1, 1A3 and 1A8 microsomes was effectively inhibited by phenylbutazone.

Metabolic fate of a new histamine H2-receptor antagonist, Z-300, in rat

1. The metabolisms of (N-[3-[3-(piperidinomethyl)-phenoxy]-propyl]-2-(2-hydroxyethyl-1-t hio) acetamide.2-(4-hydroxy benzoyl) benzoate) (Z-300) in rat have been studied. 2. Five metabolites were identified, namely as sulphoxide (SO), sulphone (SO2), phenol (PH), propionic acid (PA) and thioglycolic acid (TH) derivatives. Furthermore, two metabolites, M1 and M2, were proposed as PH conjugated to sulphate and glucuronic acid respectively. 3. There were several metabolites in the plasma samples after oral administration of 14C-Z-300. The main metabolites were PH (mainly existing as the conjugates), PA and SO, and Z-300 was present as <10% of plasma radioactivity. SO2 and TH were barely detected in plasma. In the urine, the composition of metabolites was similar to that in the plasma. In the bile, the main metabolites were PH (mainly existing as the conjugates) and SO, and very little Z-300 was detected. In the faeces, the main compounds were Z-300 and TH, and very little SO. 4. During in vitro metabolism experiments, Z-300 was metabolized mainly in the liver, and a little in the kidney and lung. PH or its conjugates were produced little on the in vitro metabolism; this was different from that found in the in vivo experiment. The metabolism of Z-300 to SO appeared to be catalysed by cytochrome P450 rather than the flavin-containing monooxygenase. Both hepatic clearance calculated from the in vitro metabolism experiment (28.9 ml/min/kg) and from the liver perfusion metabolism experiment (33.2 ml/min/kg) were lower than the total clearance from the in vivo experiment (119 ml/min/kg). 5. After administration of 14C-SO into the caecum of rat, SO was converted to Z-300 and excreted in the faeces. 14C-SO was reduced by incubation for 4 h within the caecal content of rat under an anaerobic condition in vitro. The results suggest that SO was reduced to Z-300 by gut flora in the lower intestine.

Kinetics of diazepam metabolism in rat hepatic microsomes and hepatocytes and their use in predicting in vivo hepatic clearance

1. The rates of diazepam (DZ) metabolism to the primary metabolites 3-hydroxydiazepam, 4'-hydroxydiazepam and nordiazepam were studied in vitro using rat hepatic microsomes and hepatocytes. 4'-hydroxydiazepam had the largest intrinsic clearance (Vmax/Km ratio, CL(int)) in both microsomes and hepatocytes representing 49 and 70% of total metabolism respectively. Whereas the contribution of 3-hydroxydiazepam was similar in both systems (21-24%), the N-demethylation pathway was greater in microsomes (27%) than hepatocytes (9%). 2. The pharmacokinetics of DZ were determined in vivo using the intraportal route to avoid blood flow limitations due to the high clearance of DZ. No dose dependency was observed in either clearance or steady state volume of distribution, which were estimated to be 38 ml/min/SRW (where SRW is a standard rat weight of 250 g) and 1.3 L/SRW respectively. Blood binding of DZ was concentration independent, the unbound fraction being 0.22. 3. Scaling factors were used to relate the in vitro CL(int) to the in vivo unbound clearance. Hepatocytes (123 ml/min/SRW) produced a more realistic prediction for the in vivo value (174 ml/min/SRW) than microsomes (41 ml/min/SRW). This situation is believed to arise from the quantitative differences in the three metabolic pathways in the two in vitro systems. It is speculated that end product inhibition is responsible for reduced total metabolism in microsomes whereas hepatocytes operate kinetically in a manner close to in vivo.

Concentration-dependent metabolism of diazepam in mouse liver

Previous mouse liver studies with diazepam (DZ), N-desmethyldiazepam (NZ), and temazepam (TZ) confirmed that under first-order conditions, DZ formed NZ and TZ in parallel. Oxazepam (OZ) was generated via NZ and not TZ despite that preformed NZ and TZ were both capable of forming OZ. In the present studies, the concentration-dependent sequential metabolism of DZ was studied in perfused mouse livers and microsomes, with the aim of distinguishing the relative importance of NZ and TZ as precursors of OZ. In microsomal studies, the Kms and Vmaxs, corrected for binding to microsomal proteins, were 34 microM and 3.6 nmole/min per mg and 239 microM and 18 nmole/min per mg, respectively, for N-demethylation and C3-hydroxylation of DZ. The Kms and Vmaxs for N-demethylation and C3-hydroxylation of TZ and NZ, respectively, to form OZ, were 58 microM and 2.5 nmole/min per mg and 311 microM and 2 nmole/min per mg, respectively. The constants suggest that at low DZ concentrations, NZ formation predominates and is a major source of OZ, whereas at higher DZ concentrations, TZ is the important source of OZ. In livers perfused with DZ at input concentrations of 13 to 35 microM, the extraction ratio of DZ (E[DZ]) decreased from 0.83 to 0.60. NZ was the major metabolite formed although its appearance was less than proportionate with increasing DZ input concentration. By contrast, the formation of TZ increased disproportionately with increasing DZ concentration, whereas that for OZ decreased and paralleled the behavior of NZ. Computer simulations based on a tubular flow model and the in vitro enzymatic parameters provided a poor in vitro-organ correlation. The E[DZ], appearance rates of the metabolites, and the extraction ratio of formed NZ (E[NZ, DZ]) were poorly predicted; TZ was incorrectly identified as the major precursor of OZ. Simulations with optimized parameters improved the correlations and identified NZ as the major contributor of OZ. Saturation of DZ N-demethylation at higher DZ concentrations increased the role of TZ in the formation of OZ. The poor aqueous solubility (limiting the concentration range of substrates used in vitro), avid tissue binding and the coupling of enzymatic reactions in liver, favoring sequential metabolism, are possible explanations for the poor in vitro-organ correlation. This work emphasizes the complexity of the hepatic intracellular milieu for drug metabolism and the need for additional modeling efforts to adequately describe metabolite kinetics.

Prediction of ACNU plasma concentration-time profiles in humans by animal scale-up

Plasma concentration-time profiles of nimustine hydrochloride, 1-[(4-amino-2-methyl-5-pyrimidinyl)methyl]-3-(2-chloroethyl)-3-nitrosour ea hydrochloride (ACNU), in the mouse, rat, rabbit, and dog were determined by high-performance liquid chromatographic analysis. The pharmacokinetic parameters for these four animal species and previously reported clinical data were analyzed for investigation of interspecies correlation. Log-log plots of body weight (W; kg) vs total plasma clearance (CLtot,p; ml/min) and steady-state distribution volume (Vd,ss; l) for the four animal species were linear, with high correlation coefficients (r 0.996 for both parameters), despite the fact that the nonrenal clearance was greater than 97% in these species. Linear regression on the plots excluding human data yielded allometric equations (CLtot,p = 50.6 W0.957; Vd, ss = 1.29 W1.03) that were extrapolated to predict ACNU pharmacokinetic parameters in humans. For both parameters, however, there were 3-fold differences between the predicted and observed parametric values. To investigate these discrepancies, we measured serum protein binding of ACNU in these animal species and in humans. The values of CLtot,p and Vd,ss were converted into those of CLutot,p and Vd,uss, which correspond to the parameters for unbound ACNU. In this case, correlation coefficients of the log-log plots excluding human data (CLutot,p = 71.7 W0.891; Vd,uss = 1.82 W0.966) were also high (r greater than or equal to 0.991). The extrapolated values vs those observed in a 70-kg human were the following: CLutot,p, 3,160 vs 2,290 ml/min; Vd,uss, 110 vs 106 l. Thus, the animal data were successfully extrapolated to yield better predictions of human pharmacokinetic parameters if the analysis was based on the unbound plasma concentration of ACNU. In addition, the predicted plasma concentration-time profile for humans also showed good agreement with the observed ones. These results suggest the importance of measuring unbound fractions of drugs for more accurate prediction of human pharmacokinetic parameters by extrapolation of animal data to the human situation.

Mephenytoin stereoselective elimination in the rat: II. Comparison of mephenytoin stereoselective clearance during chronic intravenous and hepatic portal vein administration

The stereoselective clearances of R- and S-mephenytoin were determined in rats receiving either an intravenous or hepatic portal vein infusion of racemic mephenytoin. The mean +/- SD intravenous clearances of R- and S-mephenytoin were 1630 +/- 250 ml/hr and 630 +/- 250 ml/hr, respectively. The corresponding portal vein clearances for these enantiomers were 2560 +/- 1230 ml/hr (R-mephenytoin) and 540 +/- 230 ml/hr (S-mephenytoin). In spite of the slightly higher clearance for R-mephenytoin following portal vein administration, the difference between the intravenous and portal vein clearances for R- or S-mephenytoin were not found to be significant. Subsequent computer simulations of the data indicated there was less than a 5% probability that this result could be attributed solely to interanimal variability in drug clearance. The estimated extraction ratio of R-mephenytoin by the liver was modest and suggested mephenytoin may undergo a substantial degree of extrahepatic elimination in the rat.

Kinetic analysis of the elimination process of human epidermal growth factor (hEGF) in rats

Pharmacokinetic study of human epidermal growth factor (hEGF) in rats was performed in vivo. The hepatic extraction ratio (EH) of [125I]hEGF, determined from the difference between the artery and the hepatic vein plasma concentrations at steady state, was 0.19. The hepatic clearance (CLH:7.56 ml/min/kg body wt), calculated by multiplying EH by the hepatic plasma flow rate (QP,H), was approximately 70% of the total body clearance (CLtot: 10.8 ml/min/kg body wt), which was determined from the steady-state arterial plasma concentration and the infusion rate. These results indicated that the liver is the main organ responsible for the removal of [125I]hEGF from the systemic circulation in rats. The renal extraction ratio (ER) of [125I]hEGF was half of that of [14C]inulin; this may have resulted from the plasma protein binding of [125I]hEGF, which was approximately 50% as determined by the charcoal adsorption method and the equilibrium gel-filtration method. The renal clearance (CLR:2.65 ml/min/kg body wt), calculated by multiplying ER by the renal plasma flow rate (QPR), was approximately 17% of the CLtot (15.6 ml/min/kg body wt), indicating a minor contribution of CLR to CLtot compared with that of CLH to CLtot. The CLR of [125I]hEGF calculated from the urinary excretion data was one-tenth of that calculated from the plasma concentration difference between the femoral artery and the renal vein at steady state. These results suggest that the bulk of [125I]hEGF cleared from the plasma by the kidney may have been metabolized further in the renal tubules before appearing in the urine.

Man versus beast: pharmacokinetic scaling in mammals

Land mammals range in size from the 3-g shrew to the 3000-kg elephant. Despite this 10(6) range in weight, most land mammals have similar anatomy, physiology, biochemistry, and cellular structure. This similarity has allowed interspecies scaling of physiologic properties such as heart rate, blood flow, blood volume, organ size, and longevity. The equation that is the basis for scaling physiologic properties among mammals is the power equation Y = aWb, where Y is the physiologic variable of interest, W is body weight, and log a is the y-intercept and b is the slope obtained from the plot of log Y versus log W. Animals commonly used in preclinical drug studies (i.e., mice, rats, rabbits, monkeys, and dogs) do not eliminate drugs at the same rate that humans eliminate drugs; small mammals usually eliminate drugs faster than large mammals. Since drug elimination is intimately associated with physiologic properties that are well described among species, it seems reasonable to surmise that drug elimination can be scaled among mammals. Analysis of drug pharmacokinetics in numerous species demonstrates that drug elimination among species is predictable and, in general, obeys the power equation Y = aWb. Early papers on interspecies pharmacokinetic scaling normalized the x- and y-axes to illustrate the superimpossibility of pharmacokinetic curves from different species. More recently, the x- and y-axes have been left in the common units of concentration and time, and individual pharmacokinetic variables have been adjusted to predict pharmacokinetic profiles in an untested species, usually humans.

A simulation study on the effect of a uniform diffusional barrier across hepatocytes on drug metabolism by evenly or unevenly distributed uni-enzyme in the liver

The effect of a uniform diffusional barrier on hepatic extraction of the parent drug by evenly or unevenly distributed uni-enzyme was quantitatively determined by the present simulation study. Five models of enzymic distribution were defined with regard to the hepatic blood flow path, and the extraction ratios were calculated or simulated under the various conditions of average intrinsic clearances and diffusion clearances across hepatocytes. Differences in the extraction ratios among the five models were evaluated by the "relative extraction ratios," which are the extraction ratios in each model divided by that in the model where the enzymatic activity is evenly distributed. It was found that when a diffusion clearance was high compared to the intrinsic clearance, enzymic distribution was not an important determinant of the extent of hepatic extraction. By contrast, when a diffusional barrier across hepatocytes exists, i.e., the diffusion clearance is low or intermediate compared to the intrinsic clearance, extraction ratios differed widely among the models of enzymic distribution, especially at intermediate average intrinsic clearances. In the presence of a diffusional barrier, the more skewed the distribution of the enzymatic activity is, the lesser the amount of drug eliminated at steady state. The most efficient metabolism occurred when the enzymatic activity was evenly distributed.

Prediction of the disposition of nine weakly acidic and six weakly basic drugs in humans from pharmacokinetic parameters in rats

Various pharmacokinetic parameters--disposition half-life, t1/2,z, metabolic clearance CLm, volume of distribution V, intrinsic clearance of unbound drug CLuint, and unbound volume of distribution of tissues (distributive tissue volume/fraction of drug in tissue unbound, VT/fuT--are compared in rat and human for nine weakly acidic drugs, phenytoin, hexobarbital, pentobarbital, phenylbutazone, warfarin, tolbutamide, valproate, phenobarbital, and amobarbital, and six weakly basic drugs, quinidine, chlorpromazine, propranolol, pentazocin, antipyrine, and diazepam. With regard to all parameters, statistically significant correlations are obtained when parameters are plotted on a log-log plot. Correlation coefficients between the intrinsic parameters (CLuint or VT/fuT) were higher than those between the hybrid parameters (t1/2,z, CLm or V). In general, these drugs were metabolized ten times more rapidly in rat than in human. With regard to the tissue distribution of these drugs, there was little difference between rat and human. Predictions of CLm, V, and t1/2, in humans using rat data were successful for most drugs, with a few marked exceptions.

Prediction of diazepam disposition in the rat and man by a physiologically based pharmacokinetic model

A physiologically based pharmacokinetic model for diazepam disposition was developed in the rat, incorporating anatomical, physiological, and biochemical parameters, i.e., tissue volume, blood flow rate, serum free fraction, distribution of diazepam into red blood cells, drug metabolism and tissue-to-blood distribution ratio. The serum free fraction of diazepam was determined by equilibrium dialysis at 37 degrees C and was constant over a wide concentration range. Partition of diazepam between plasma and erythrocytes was determined in vitro at 37 degrees C, and the resultant blood-to-plasma concentration ratio was constant over a wide concentration range. The enzymatic parameters (Km, Vmax) of the eliminating organs, i.e., liver, kidney, and lung, previously determined using microsomes, were used for the prediction. The tissue-to-blood distribution ratios inferred by inspection of the data when pseudoequilibrium is reached after i.v. bolus injection of 1.2 mg/kg diazepam were corrected according to the method of Chen and Gross. Predicted diazepam concentration time-course profiles in plasma and various organs or tissues, using an 11-compartmental model, were compared with those observed. Prediction was successful in all compartments including brain, the target organ of diazepam. Scale-up of the disposition kinetics of diazepam from rat to man was also successful.