The use of a physiologically based pharmacokinetic model to evaluate deconvolution measurements of systemic absorption

PKQuest: PBPK modeling of highly lipid soluble and extracellular solutes

One of the primary objectives of physiologically based pharmacokinetics (PBPK) is the prediction of a drug’s pharmacokinetics just from knowledge of its physicochemical structure. Unfortunately, at present, the accuracy of this prediction is limited for most drugs because of uncertainty about the drug’s organ/blood partition coefficient (K). However, there are two classes of solutes which are exceptions to this: 1) the highly lipid soluble (HLS) solutes, and 2) the extracellular (ECS) solutes. Since the HLS drugs (eg, volatile anesthetics, propofol, cannabinol) have lipid/water partition coefficients (P_L/W) of 100 or greater, their K is dominated by the tissue fat fraction and one can accurately predict K just from in vitro measurements of P_L/W along with prior anatomic measurements of the fat fraction of the organs in the PBPK model. Since the ECS drugs, such as most antibiotics, cannot penetrate cells, they are not subject to the intracellular binding that complicates the prediction of K for the weak bases and acids. The ECS K is determined primarily by plasma and interstitial albumin binding and can be predicted from in vitro measurements of plasma albumin binding along with prior measurements of interstitial tissue volume and albumin concentrations. This review provides an in depth discussion of the PBPK modeling of these two drug classes along with many specific clinical examples illustrating the good PBPK predictions possible with just zero (volatile anesthetics) or 1 (the clearance) adjustable parameter. The PBPK analysis uses PKQuest, a freely distributed, general purpose pharmacokinetic program. PKQuest is designed so that application to the HLS and ECS solute classes is especially easy. The user only needs to enter the specific parameters that are required to characterize the drug (eg, P_L/W for HLS or plasma albumin binding for ECS) with all the other PBPK parameters (organ blood flow, fat fraction, extracellular volumes, etc.) are set by default.

Using deconvolution to understand the mechanism for variable plasma concentration-time profiles after intramuscular injection

To introduce better antibiotics for the treatment of some infectious diseases in sheep and to expand the range of antibiotics available for veterinary medicine, pharmacokinetics of two antibiotics marbofloxacin (MBX) and trovafloxacin (TVX) were investigated in sheep after intramuscular injection. Variable and irregular plasma concentration-time profiles were observed for TVX but not for MBX. To understand the mechanism of this phenomenon, intravenous studies were performed for both drugs and data were analyzed using a population approach. Deconvolution was then performed using various approaches to obtain absorption profiles of both drugs in sheep after intramuscular injection. The Loo-Riegelman and staircase deconvolution function methods were found to provide more reliable estimates of absorption rate than the Spath-spline and B-spline constraining break points deconvolution methods. The absorption profiles resulting from deconvolution indicated a zero-order absorption process for TVX and a first-order process for MBX. Precipitation of TVX at the injection site was suspected to cause the pseudo zero-order absorption. This hypothesis was supported by the observation of crystalline deposits of TVX in sheep meat after direct injection, using reflectance confocal microscopy.

Quantitation of small intestinal permeability during normal human drug absorption

BACKGROUND

Understanding the quantitative relationship between a drug's physical chemical properties and its rate of intestinal absorption (QSAR) is critical for selecting candidate drugs. Because of limited experimental human small intestinal permeability data, approximate surrogates such as the fraction absorbed or Caco-2 permeability are used, both of which have limitations.

METHODS

Given the blood concentration following an oral and intravenous dose, the time course of intestinal absorption in humans was determined by deconvolution and related to the intestinal permeability by the use of a new 3 parameter model function ("Averaged Model" (AM)). The theoretical validity of this AM model was evaluated by comparing it to the standard diffusion-convection model (DC). This analysis was applied to 90 drugs using previously published data. Only drugs that were administered in oral solution form to fasting subjects were considered so that the rate of gastric emptying was approximately known. All the calculations are carried out using the freely available routine PKQuest Java (http://www.pkquest.com) which has an easy to use, simple interface.

RESULTS

Theoretically, the AM permeability provides an accurate estimate of the intestinal DC permeability for solutes whose absorption ranges from 1% to 99%. The experimental human AM permeabilities determined by deconvolution are similar to those determined by direct human jejunal perfusion. The small intestinal pH varies with position and the results are interpreted in terms of the pH dependent octanol partition. The permeability versus partition relations are presented separately for the uncharged, basic, acidic and charged solutes. The small uncharged solutes caffeine, acetaminophen and antipyrine have very high permeabilities (about 20 x 10-4 cm/sec) corresponding to an unstirred layer of only 45 μm. The weak acid aspirin also has a large AM permeability despite its low octanol partition at pH 7.4, suggesting that it is nearly completely absorbed in the first part of the intestine where the pH is about 5.4.

CONCLUSIONS

The AM deconvolution method provides an accurate estimate of the human intestinal permeability. The results for these 90 drugs should provide a useful benchmark for evaluating QSAR models.

On the search for in vitro in vivo correlations in the field of intra-articular drug delivery: administration of sodium diatrizoate to the horse

Development of suitable in vitro release models for formulation development as well as quality control purposes has to be initiated in the early design phase of injectable depots. Optimally, construction of an in vitro release model may lead to the establishment of in vitro in vivo correlations. By using a model compound (sodium diatrizoate, DTZ), the purpose of this study was to investigate the possibility of establishing in vitro in vivo relations between the DTZ disappearance profile obtained from the donor compartment of the rotating dialysis cell model and the joint disappearance profile following intra-articular administration. In vitro experiments were conducted by applying solutions of DTZ to the donor compartment. In the in vivo experiments, five horses were subjected to both intravenous and intra-articular administration of an aqueous solution of 3.9 mg DTZ/kg. A strong relation (R(2)=0.99) was obtained between the disappearance data from the donor compartment of the in vitro model and the disappearance data from the synovial fluid after intra-articular administration of DTZ. Furthermore, a relation (R(2)=0.91) between the appearance data obtained from the acceptor compartment and the deconvolved appearance serum data upon intra-articular administration of DTZ was obtained. The correlations obtained in this study hold promise that the rotating dialysis cell model has a role in the prediction of the intra-articular fate of drugs injected as solutions.

Comparison of spectroscopically measured tissue alcohol concentration to blood and breath alcohol measurements

Alcohol testing is an expanding area of interest due to the impacts of alcohol abuse that extend well beyond drunk driving. However, existing approaches such as blood and urine assays are hampered in some testing environments by biohazard risks. A noninvasive, in vivo spectroscopic technique offers a promising alternative, as no body fluids are required. The purpose of this work is to report the results of a 36-subject clinical study designed to characterize tissue alcohol measured using near-infrared spectroscopy relative to venous blood, capillary blood, and breath alcohol. Comparison of blood and breath alcohol concentrations demonstrated significant differences in alcohol concentration [root mean square of 9.0 to 13.5 mg/dL] that were attributable to both assay accuracy and precision as well as alcohol pharmacokinetics. A first-order kinetic model was used to estimate the contribution of alcohol pharmacokinetics to the differences in concentration observed between the blood, breath, and tissue assays. All pair-wise combinations of alcohol assays were investigated, and the fraction of the alcohol concentration variance explained by pharmacokinetics ranged from 41.0% to 83.5%. Accounting for pharmacokinetic concentration differences, the accuracy and precision of the spectroscopic tissue assay were found to be comparable to those of the blood and breath assays.

PKQuest_Java: free, interactive physiologically based pharmacokinetic software package and tutorial

Background

Physiologically based pharmacokinetics (PBPK) uses a realistic organ model to describe drug kinetics. The blood-tissue exchange of each organ is characterized by its volume, perfusion, metabolism, capillary permeability and blood/tissue partition coefficient. PBPK applications require both sophisticated mathematical modeling software and a reliable complete set of physiological parameters. Currently there are no software packages available that combine ease of use with the versatility that is required of a general PBPK program.

Findings

The program is written in Java and is available for free download at . Included in the download is a detailed tutorial that discusses the pharmacokinetics of 6 solutes (D₂O, amoxicillin, desflurane, propofol, ethanol and thiopental) illustrated using experimental human pharmacokinetic data. The complete PBPK description for each solute is stored in Excel spreadsheets that are included in the download. The main features of the program are: 1) Intuitive and versatile interactive interface; 2) Absolute and semi-logarithmic graphical output; 3) Pre-programmed optimized human parameter data set (but, arbitrary values can be input); 4) Time dependent changes in the PBPK parameters; 5) Non-linear parameter optimization; 6) Unique approach to determine the oral "first pass metabolism" of non-linear solutes (e.g. ethanol); 7) Pulmonary perfusion/ventilation heterogeneity for volatile solutes; 8) Input and output of Excel spreadsheet data; 9) Antecubital vein sampling.

Conclusion

PKQuest_Java is a free, easy to use, interactive PBPK software routine. The user can either directly use the pre-programmed optimized human or rat data set, or enter an arbitrary data set. It is designed so that drugs that are classified as "extracellular" or "highly fat soluble" do not require information about tissue/blood partition coefficients and can be modeled by a minimum of user input parameters. PKQuest_Java, along with the included tutorial, could be used as the basis of an interactive, on-line, pharmacokinetic course.

An efficient deconvolution algorithm for estimating oxygen consumption during muscle activities

The reconstruction of an unknown input function from noisy measurements in a biological system is an ill-posed inverse problem. Any computational algorithm for its solution must use some kind of regularization technique to neutralize the disastrous effects of amplified noise components on the computed solution. In this paper, following a hierarchical Bayesian statistical inversion approach, we seek estimates for the input function and regularization parameter (hyperparameter) that maximize the posterior probability density function. We solve the maximization problem simultaneously for all unknowns, hyperparameter included, by a suitably chosen quasi-Newton method. The optimization approach is compared to the sampling-based Bayesian approach. We demonstrate the efficiency and robustness of the deconvolution algorithm by applying it to reconstructing the time courses of mitochondrial oxygen consumption during muscle state transitions (e.g., from resting state to contraction and recovery), from the simulated noisy output of oxygen concentration dynamics on the muscle surface. The model of oxygen transport and metabolism in skeletal muscle assumes an in vitro cylindrical structure of the muscle in which the oxygen from the surrounding oxygenated solution diffuses into the muscle and is then consumed by the muscle mitochondria. The algorithm can be applied to other deconvolution problems by suitably replacing the forward model of the system.

125I used for labelling of proteins in an absorption model changes the absorption rate of insulin aspart

The aim of this study is to validate the ability of the disappearance model to predict absorption rates of insulin aspart in pigs. The disappearance model is used as a screening tool to estimate absorption rates after subcutaneous injections in humans or pigs especially of insulin and insulin analogues. The disappearance model measures remaining radioactivity at the injection site and therefore radioactive labelling of the insulin analogue is necessary. The labelling is done with 125I. One of the assumptions for the disappearance model to be reliable is that absorption rates of the labelled and non-labelled molecules are comparable. In this study, we compared disappearance data with absorption calculated from plasma samples of insulin aspart. The calculated absorption is based on non-labelled insulin aspart. The absorption rate from the disappearance data was statistical significant (p = 0.0028) different from the absorption rate based on plasma samples. A control study was carried out where 125I labelled insulin aspart was compared to 127I (the natural non-radioactive isotope) insulin aspart. In this study, absorption rate from the disappearance data and absorption rate based on plasma samples were similar (p = 0.63).

CONCLUSION

Iodination of insulin aspart changes the subcutaneous absorption rate.

Heterogeneity of human adipose blood flow

BACKGROUND

The long time pharmacokinetics of highly lipid soluble compounds is dominated by blood-adipose tissue exchange and depends on the magnitude and heterogeneity of adipose blood flow. Because the adipose tissue is an infinite sink at short times (hours), the kinetics must be followed for days in order to determine if the adipose perfusion is heterogeneous. The purpose of this paper is to quantitate human adipose blood flow heterogeneity and determine its importance for human pharmacokinetics.

METHODS

The heterogeneity was determined using a physiologically based pharmacokinetic model (PBPK) to describe the 6 day volatile anesthetic data previously published by Yasuda et. al. The analysis uses the freely available software PKQuest and incorporates perfusion-ventilation mismatch and time dependent parameters that varied from the anesthetized to the ambulatory period. This heterogeneous adipose perfusion PBPK model was then tested by applying it to the previously published cannabidiol data of Ohlsson et. al. and the cannabinol data of Johansson et. al.

RESULTS

The volatile anesthetic kinetics at early times have only a weak dependence on adipose blood flow while at long times the pharmacokinetics are dominated by the adipose flow and are independent of muscle blood flow. At least 2 adipose compartments with different perfusion rates (0.074 and 0.014 l/kg/min) were needed to describe the anesthetic data. This heterogeneous adipose PBPK model also provided a good fit to the cannabinol data.

CONCLUSION

Human adipose blood flow is markedly heterogeneous, varying by at least 5 fold. This heterogeneity significantly influences the long time pharmacokinetics of the volatile anesthetics and tetrahydrocannabinol. In contrast, using this same PBPK model it can be shown that the long time pharmacokinetics of the persistent lipophilic compounds (dioxins, PCBs) do not depend on adipose blood flow. The ability of the same PBPK model to describe both the anesthetic and cannabinol kinetics provides direct qualitative evidence that their kinetics are flow limited and that there is no significant adipose tissue diffusion limitation.

Human physiologically based pharmacokinetic model for ACE inhibitors: ramipril and ramiprilat

Background

The angiotensin-converting enzyme (ACE) inhibitors have complicated and poorly characterized pharmacokinetics. There are two binding sites per ACE (high affinity "C", lower affinity "N") that have sub-nanomolar affinities and dissociation rates of hours. Most inhibitors are given orally in a prodrug form that is systemically converted to the active form. This paper describes the first human physiologically based pharmacokinetic (PBPK) model of this drug class.

Methods

The model was applied to the experimental data of van Griensven et. al for the pharmacokinetics of ramiprilat and its prodrug ramipril. It describes the time course of the inhibition of the N and C ACE sites in plasma and the different tissues. The model includes: 1) two independent ACE binding sites; 2) non-equilibrium time dependent binding; 3) liver and kidney ramipril intracellular uptake, conversion to ramiprilat and extrusion from the cell; 4) intestinal ramipril absorption. The experimental in vitro ramiprilat/ACE binding kinetics at 4°C and 300 mM NaCl were assumed for most of the PBPK calculations. The model was incorporated into the freely distributed PBPK program PKQuest.

Results

The PBPK model provides an accurate description of the individual variation of the plasma ramipril and ramiprilat and the ramiprilat renal clearance following IV ramiprilat and IV and oral ramipril. Summary of model features: Less than 2% of total body ACE is in plasma; 35% of the oral dose is absorbed; 75% of the ramipril metabolism is hepatic and 25% of this is converted to systemic ramiprilat; 100% of renal ramipril metabolism is converted to systemic ramiprilat. The inhibition was long lasting, with 80% of the C site and 33% of the N site inhibited 24 hours following a 2.5 mg oral ramipril dose. The plasma ACE inhibition determined by the standard assay is significantly less than the true in vivo inhibition because of assay dilution.

Conclusion

If the in vitro plasma binding kinetics of the ACE inhibitor for the two binding sites are known, a unique PBPK model description of the Griensven et. al. experimental data can be obtained.

Human physiologically based pharmacokinetic model for propofol

Background

Propofol is widely used for both short-term anesthesia and long-term sedation. It has unusual pharmacokinetics because of its high lipid solubility. The standard approach to describing the pharmacokinetics is by a multi-compartmental model. This paper presents the first detailed human physiologically based pharmacokinetic (PBPK) model for propofol.

Methods

PKQuest, a freely distributed software routine , was used for all the calculations. The "standard human" PBPK parameters developed in previous applications is used. It is assumed that the blood and tissue binding is determined by simple partition into the tissue lipid, which is characterized by two previously determined set of parameters: 1) the value of the propofol oil/water partition coefficient; 2) the lipid fraction in the blood and tissues. The model was fit to the individual experimental data of Schnider et. al., Anesthesiology, 1998; 88:1170 in which an initial bolus dose was followed 60 minutes later by a one hour constant infusion.

Results

The PBPK model provides a good description of the experimental data over a large range of input dosage, subject age and fat fraction. Only one adjustable parameter (the liver clearance) is required to describe the constant infusion phase for each individual subject. In order to fit the bolus injection phase, for 10 or the 24 subjects it was necessary to assume that a fraction of the bolus dose was sequestered and then slowly released from the lungs (characterized by two additional parameters). The average weighted residual error (WRE) of the PBPK model fit to the both the bolus and infusion phases was 15%; similar to the WRE for just the constant infusion phase obtained by Schnider et. al. using a 6-parameter NONMEM compartmental model.

Conclusion

A PBPK model using standard human parameters and a simple description of tissue binding provides a good description of human propofol kinetics. The major advantage of a PBPK model is that it can be used to predict the changes in kinetics produced by variations in physiological parameters. As one example, the model simulation of the changes in pharmacokinetics for morbidly obese subjects is discussed.

Monitoring the initial pulmonary absorption of two different beclomethasone dipropionate aerosols employing a human lung reperfusion model

BACKGROUND

The pulmonary residence time of inhaled glucocorticoids as well as their rate and extend of absorption into systemic circulation are important facets of their efficacy-safety profile. We evaluated a novel approach to elucidate the pulmonary absorption of an inhaled glucocorticoid. Our objective was to monitor and compare the combined process of drug particle dissolution, pro-drug activation and time course of initial distribution from human lung tissue into plasma for two different glucocorticoid formulations.

METHODS

We chose beclomethasone dipropionate (BDP) delivered by two different commercially available HFA-propelled metered dose inhalers (Sanasthmax/Becloforte and Ventolair/Qvar). Initially we developed a simple dialysis model to assess the transfer of BDP and its active metabolite from human lung homogenate into human plasma. In a novel experimental setting we then administered the aerosols into the bronchus of an extracorporally ventilated and reperfused human lung lobe and monitored the concentrations of BDP and its metabolites in the reperfusion fluid.

RESULTS

Unexpectedly, we observed differences between the two aerosol formulations Sanasthmax/Becloforte and Ventolair/Qvar in both the dialysis as well as in the human reperfusion model. The HFA-BDP formulated as Ventolair/Qvar displayed a more rapid release from lung tissue compared to Sanasthmax/Becloforte. We succeeded to explain and illustrate the observed differences between the two aerosols with their unique particle topology and divergent dissolution behaviour in human bronchial fluid.

CONCLUSION

We conclude that though the ultrafine particles of Ventolair/Qvar are beneficial for high lung deposition, they also yield a less desired more rapid systemic drug delivery. While the differences between Sanasthmax/Becloforte and Ventolair/Qvar were obvious in both the dialysis and lung perfusion experiments, the latter allowed to record time courses of pro-drug activation and distribution that were more consistent with results of comparable clinical trials. Thus, the extracorporally reperfused and ventilated human lung is a highly valuable physiological model to explore the lung pharmacokinetics of inhaled drugs.

Percutaneous absorption of m-xylene vapour in volunteers during pre-steady and steady state

Percutaneous absorption of m-xylene (XYL) was determined in volunteers exposed to 29.4 microg cm(-3) XYL vapour on the forearm and hand for 20, 45, 120 and 180 min. The internal exposure was assessed by measuring the concentration of XYL in exhaled air. The systemic kinetics were determined using a reference exposure by inhalation. The dermal permeation rate and the cumulative absorption of XYL as a function of time were calculated using mathematical deconvolution. From these relationships, the average flux into the skin throughout the exposure (J(skin, average)) and the maximal flux into the blood (J(blood, max)) were derived. Both fluxes were dependent on the duration of exposure, approaching each other at longer exposure durations. The values of J(skin, average), adjusted to a concentration of 1 microg cm(-3), were 0.091 microg cm(-2) h(-1) during 20-min exposure falling to 0.072, 0.066 and 0.061 microg cm(-2) h(-1) for 45, 120 and 180 min, respectively. The values of J(blood, max) showed an opposite trend, gradually increasing from 0.034 microg cm(-2) h(-1) at an exposure duration of 20 min to 0.042, 0.059 and 0.063 microg cm(-2) h(-1) for 45, 120 and 180 min of exposure durations, respectively.

Physiologically based pharmacokinetic modeling of arterial - antecubital vein concentration difference

BACKGROUND

Modeling of pharmacokinetic parameters and pharmacodynamic actions requires knowledge of the arterial blood concentration. In most cases, experimental measurements are only available for a peripheral vein (usually antecubital) whose concentration may differ significantly from both arterial and central vein concentration.

METHODS

A physiologically based pharmacokinetic (PBPK) model for the tissues drained by the antecubital vein (referred to as "arm") is developed. It is assumed that the "arm" is composed of tissues with identical properties (partition coefficient, blood flow/gm) as the whole body tissues plus a new "tissue" representing skin arteriovenous shunts. The antecubital vein concentration depends on the following parameters: the fraction of "arm" blood flow contributed by muscle, skin, adipose, connective tissue and arteriovenous shunts, and the flow per gram of the arteriovenous shunt. The value of these parameters was investigated using simultaneous experimental measurements of arterial and antecubital concentrations for eight solutes: ethanol, thiopental, 99Tcm-diethylene triamine pentaacetate (DTPA), ketamine, D2O, acetone, methylene chloride and toluene. A new procedure is described that can be used to determine the arterial concentration for an arbitrary solute by deconvolution of the antecubital concentration. These procedures are implemented in PKQuest, a general PBPK program that is freely distributed http://www.pkquest.com.

RESULTS

One set of "standard arm" parameters provides an adequate description of the arterial/antecubital vein concentration for ethanol, DTPA, thiopental and ketamine. A significantly different set of "arm" parameters was required to describe the data for D2O, acetone, methylene chloride and toluene - probably because the "arm" is in a different physiological state.

CONCLUSIONS

Using the set of "standard arm" parameters, the antecubital vein concentration can be used to determine the whole body PBPK model parameters for an arbitrary solute without any additional adjustable parameters. Also, the antecubital vein concentration can be used to estimate the arterial concentration for an arbitrary input for solutes for which no arterial concentration data is available.

Physiologically based pharmacokinetic modeling of arterial - antecubital vein concentration difference

BACKGROUND

Modeling of pharmacokinetic parameters and pharmacodynamic actions requires knowledge of the arterial blood concentration. In most cases, experimental measurements are only available for a peripheral vein (usually antecubital) whose concentration may differ significantly from both arterial and central vein concentration.

METHODS

A physiologically based pharmacokinetic (PBPK) model for the tissues drained by the antecubital vein (referred to as "arm") is developed. It is assumed that the "arm" is composed of tissues with identical properties (partition coefficient, blood flow/gm) as the whole body tissues plus a new "tissue" representing skin arteriovenous shunts. The antecubital vein concentration depends on the following parameters: the fraction of "arm" blood flow contributed by muscle, skin, adipose, connective tissue and arteriovenous shunts, and the flow per gram of the arteriovenous shunt. The value of these parameters was investigated using simultaneous experimental measurements of arterial and antecubital concentrations for eight solutes: ethanol, thiopental, 99Tcm-diethylene triamine pentaacetate (DTPA), ketamine, D2O, acetone, methylene chloride and toluene. A new procedure is described that can be used to determine the arterial concentration for an arbitrary solute by deconvolution of the antecubital concentration. These procedures are implemented in PKQuest, a general PBPK program that is freely distributed http://www.pkquest.com.

RESULTS

One set of "standard arm" parameters provides an adequate description of the arterial/antecubital vein concentration for ethanol, DTPA, thiopental and ketamine. A significantly different set of "arm" parameters was required to describe the data for D2O, acetone, methylene chloride and toluene - probably because the "arm" is in a different physiological state.

CONCLUSIONS

Using the set of "standard arm" parameters, the antecubital vein concentration can be used to determine the whole body PBPK model parameters for an arbitrary solute without any additional adjustable parameters. Also, the antecubital vein concentration can be used to estimate the arterial concentration for an arbitrary input for solutes for which no arterial concentration data is available.

The pharmacokinetics of the interstitial space in humans

BACKGROUND

The pharmacokinetics of extracellular solutes is determined by the blood-tissue exchange kinetics and the volume of distribution in the interstitial space in the different organs. This information can be used to develop a general physiologically based pharmacokinetic (PBPK) model applicable to most extracellular solutes.

METHODS

The human pharmacokinetic literature was surveyed to tabulate the steady state and equilibrium volume of distribution of the solutes mannitol, EDTA, morphine-6-glucuronide, morphine-3-glucuronide, inulin and beta-lactam antibiotics with a range of protein binding (amoxicillin, piperacillin, cefatrizine, ceforanide, flucloxacillin, dicloxacillin). A PBPK data set was developed for extracellular solutes based on the literature for interstitial organ volumes. The program PKQuest was used to generate the PBPK model predictions. The pharmacokinetics of the protein (albumin) bound beta-lactam antibiotics were characterized by two parameters: 1) the free fraction of the solute in plasma; 2) the interstitial albumin concentration. A new approach to estimating the capillary permeability is described, based on the pharmacokinetics of the highly protein bound antibiotics.

RESULTS

About 42% of the total body water is extracellular. There is a large variation in the organ distribution of this water - varying from about 13% of total tissue water for skeletal muscle, up to 70% for skin and connective tissue. The weakly bound antibiotics have flow limited capillary-tissue exchange kinetics. The highly protein bound antibiotics have a significant capillary permeability limitation. The experimental pharmacokinetics of the 11 solutes is well described using the new PBPK data set and PKQuest.

CONCLUSIONS

Only one adjustable parameter (systemic clearance) is required to completely characterize the PBPK for these extracellular solutes. Knowledge of just this systemic clearance allows one to predict the complete time course of the absolute drug concentrations in the major organs. PKQuest is freely available http://www.pkquest.com.