Skip to main content


The pharmacokinetics of the interstitial space in humans

Article metrics



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.


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 β-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 β-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.


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.


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


The distribution of most hydrophilic solutes (e.g. β-lactam antibiotics) is limited to the extracellular space because of their cell membrane impermeability. Thus, the pharmacokinetics of this large class of compounds is determined primarily by the volume of distribution and the kinetics of exchange with the interstitial space. Given the general importance of this class of drugs, there has been surprisingly little analysis of the pharmacokinetics of the interstitial space in humans. The 1959 review article by Edelman et. al. [1] remains the standard reference on the extracellular space. In addition, there are two extensive reviews of the structure and function of the interstitial space [2, 3].

An accurate pharmacokinetic description of the interstitial space is essential for the development of a physiologically based pharmacokinetic model (PBPK) for extracellular solutes. Although PBPK models have been used extensively to describe human pharmacokinetics, nearly all of these studies have involved solutes that have intracellular distributions, and, thus, do not require detailed modeling of the interstitial space. One exception is the recent application of PKQuest to the extracellular solutes inulin and the beta-lactam antibiotics [4]. PKQuest is a new PBPK program that has been applied to the human pharmacokinetics of a large number of solutes [48]. Although the agreement between the PKQuest PBPK model predictions and the experimental data for these extracellular solutes was satisfactory, subsequent application of PKQuest to other extracellular solutes demonstrated that there was a small, but systematic, error in the PKQuest predictions. This error is the stimulus for this more in depth analysis of interstitial pharmacokinetics.

Probably the most detailed application of a PBPK model to an extracellular solute is the study of Tsuji et. al. [9] of the pharmacokinetics of cefazolin in the rat. In this investigation, nearly every PBPK parameter required by the model was directly measured, including organ blood flows and weights, renal and hepatic clearance and the time dependence of cefazolin concentration in the different organs. Tsuji et. al. [9] estimated interstitial volumes from measurements of the steady state inulin tissue concentrations. The interstitial volumes in the original PKQuest PBPK model were based on these inulin volumes of Tsuji et. al. [9] scaled to give the correct total volume in humans. As will be shown below, inulin is not a good solute to use to measure interstitial volume and this probably contributed to the error in the earlier PKQuest PBPK predictions.


All the experimental data were obtained from earlier publications. In most cases, the published data represented the average of the experimental measurements in a number of subjects and it was assumed that it represented one "average" subject. In a few cases (indicated by a N = 1 in Table), data for a single subject were published and used in the calculations. The program UN-SCAN-IT (Silk Scientific Corporation) was used to read the data from the published figures. Most of the figures shown in this paper are direct copies (in jpeg format) of standard PKQuest output.

Measurement of steady state volume of distribution (Vss)

To ensure that a consistent calculation of Vss was used, all the values of Vss were recalculated by applying PKQuest [8] to the published venous concentration data. In this calculation, deconvolution of the intravenous input is used to find the 2 or 3 exponential unit dose bolus response function r(t):

The expression for Vss is then found from the area under the curve (AUC) and mean residence time (MRT) for this r(t) (for unit dose):

Measurement of equilibrium volume of distribution (Veq)

The direct approach to the determination of Veq is to give a constant IV infusion for a period long enough to establish a steady state plasma concentration (Ceq). For solutes in which the systemic clearance results only from renal clearance, the peripheral (tissue) and central concentrations must be in equilibrium at this steady state and Veq is defined by:

(3) V eq = A tot / C eq

where Atot is the total amount of solute present. For the special case where there is no metabolism and renal clearance in the only excretion route, Atot can be determined by quantitative urine collection after stopping the infusion. If the clearance is from just the central compartment (e.g. renal excretion) Veq (also referred to as Vdrug) is equal to Vss [10].

The value of Veq can also be estimated for subjects in renal failure that have a very low rate of renal clearance. As shown in Appendix I, for the case where the clearance is very small compared to the rate of exchange between compartments, Veq is approximately equal to Vdext (using the notation of Wagner [11]):

(4) V eq V dext = D / B β

where Bβ is the coefficient of the slow, terminal exponential term of the bolus response function (eq. (1)) and D is the bolus dose. The value of Bβ was determined using the deconvolution feature in PKQuest to find the 2-exponential response function.

PBPK model for extracellular solutes

A new PBPK model was developed that was consistent with the experimental data in Tables 1,2,3,4,5. Table 6 summarizes the new set of PBPK parameters for the "standardhuman" with a total weight of 70 Kg and a 20% fat content. The parameters are scaled for different body weight and fat content [7]. The major difference between this parameter set and that used previously [7] is the distribution of the extracellular water (ecf_water) among the different organs. The column labeled ecf_fract is the fraction of the extravascular organ water that is extracellular- i.e. interstitial. The "adipose" ecf_fract is 1 because it is assumed that all adipose tissue water is extracellular and the cells are pure lipid. The brain ecf_fract is zero because the blood brain barrier will prevent the distribution into the interstitial space of any of the test solutes used in Table 1. The term "solid" refers to the non-fat, organ solids It is assumed that the entire animal's fat is in the single "lipid" organ and that this organ is 80% fat. This is a major oversimplification since adipose blood flow may show large variations in different adipose tissues [12] and two different "fat" compartments have been used in some PBPK models [13]. An attempt has been made to minimize the number of adjustable parameters and a single fat compartment has provided satisfactory fits to the solutes previously investigated with PKQuest, including the highly fat soluble volatile anesthetics [5]. "Portal" refers to the organs drained by the portal vein (stomach, small and large intestine, pancreas and spleen). "Bone" refers to the inert solid component of bone that has no volume of distribution or blood flow. The blood flow and water component of the skeleton is distributed among the rest of the organs.

Table 1 Steady-state Volume of Distribution relative to body weight (VSS) and total body water (VSSW) for selected solutes.
Table 2 Equilibrium Volume of Distribution in Humans relative to body weight (VEQ) for selected solutes.
Table 3 Steady state volume of distribution (VSS) and interstitial volume (VI) relative to body weight for β-lactam antibiotics as a function of degree of protein binding.
Table 4 Interstitial volume (VI) – fraction of total extravascular water
Table 5 Interstitial albumin volume as fraction of interstitial EDTA volume (VI A/VI E); the interstitial albumin concentration relative to plasma albumin (CI/CP); and the product KA = (VI A/VI E)(CI/CP).
Table 6 Spreadsheat description of PBPK parameters for the 70 Kg, 20% fat "standardhuman".

The organ weights in Table 6 are close to the values in the Report on the Task Group on Reference Man [12]. The "skin" is the sum of the weight of the dermis and epidermis. The organ "other" represents primarily the loose connective tissue and the organ "tendon" represents the denser connective tissue components. The total body water is 41.84 liters, in agreement with the results of Chumlea et. al. [14]. It is assumed in Table 6 that the "blood" volume includes the intra-organ blood volume, so that the organ parameters refer to the extravascular composition.

Although the blood flows to the different organs in Table 6 are in general agreement with the "proposed reference" values of Williams and Leggett [15], they have been fine-tuned in order to optimize the PBPK model predictions for selected solutes. The resting muscle blood flow was adjusted to a value of 2.25 ml/min/100 g (0.0225 L/Kg/min) to fit the D2O data of Schloerb et. al. [16] as described previously [6]. This value for muscle blood flow is similar to the recent measurements in man of resting muscle blood flow of 2.3 [17] and 3.12 ml/min/100 gm [Ruotsalainen, #239] using [15O]H2O PET scans and of 2.5 ml/min/100 gm [18] using plethysmographic measurements of calf blood flow. This value is somewhat higher then 133Xe measurements, which vary from 1.4 to 1.9 ml/min/100 gm [19]. Similarly, the adipose blood flow in Table 6 was adjusted [5] to a value of 4.4 ml/min/100 gm to fit the enflurane data of Munson et. al. [20]. This value is close to the recent human [15O] H2O measurements of subcutaneous abdominal (4.8 ml/min/100 gm), visceral (5.9) and perirenal (4.9) fat blood flow [21]. The "tendon" (e.g. connective tissue) blood flow is poorly characterized and has a large variation. The value of 1 ml/100 gm/min is representative of direct measurements of tendon blood flow using microspheres [22, 23]. The total cardiac output of 5.62 L/min is close to the value of 5.82 L/min in normal young men reported by Grimby et. al. [24].

The volume of distribution in the different organs for solutes that bind to albumin is characterized by the parameter KA i, the ratio of the albumin concentration in the EDTA interstitial space to the plasma albumin concentration for organ i (see Appendix, II and III). The PBPK values of KA i for the different organs are listed in Table 6. In the analysis of the pharmacokinetics of the protein bound β-lactam antibiotics, an average value of KA (see eq. (10)) for the entire human is estimated. This corresponds to the organ-weighted sum:

where ecfi is the interstitial volume of organ i. The PBPK organ weighted values are listed in the last column of Table 6, and correspond to an average value of KA of 4.94/17.45 = 0.28. This is identical to the value that was estimated using the β-lactam antibiotics (see fig. 15).

As described previously [4], capillary permeability limited solutes are characterized by the organ parameter fclear [i], the fraction of solute that equilibrates with the tissue in one pass through the capillary. The fclear of organ i is related to the capillary permeability-surface area product (PSi) by the relation:

(6) fclear [i] = 1 - exp (-f P PS i / F i )

where fP is the fraction free in the plasma and Fi is the organ plasma volume flow (liters/min/Kg) [7]. The parameter fclear ranges from 0 (impermeable capillary) to 1 (flow limited). In PKQuest, the default procedure is to input just one parameter, the fclear [muscle], and the permeability of all the other organs are then automatically determined using the default values of PSi/PSmuscle programmed into PKQuest. This accounts for the other physiological factors that determine fclear (organ plasma flow and plasma protein binding). The following set of default PS values are assumed: liver, kidney and portal organs are flow limited, heart and lung have a PS 50 times that of skeletal muscle, brain has a PS = 0 for solutes that have fclear [muscle]<1 due to the blood brain barrier, and the rest of the organs have a PS equal to that of muscle [25]. As can be seen from eq. (6), solutes that have a high intrinsic capillary permeability (PS) can have a functional capillary permeability limitation (fclear < 1) if there is a high degree of plasma protein binding (i.e. fP is small). This dependence of fclear on protein binding is applied below in the application of PKQuest to a series of β-lactam antibiotics with varying protein binding.

All the PBPK calculations use PKQuest [7]. A small modification of the earlier version of the software is required to account for the new situation where the interstitial volume of distribution of protein is less than that of EDTA (see Appendix, III). Four parameters characterize the interstitial volume of distribution for each organ i: 1) ecf [i] – the interstitial volume of EDTA as the fraction of the total organ water; 2) mecf – a multiplicative constant that scales all the values of ecf [i]; 3) frecf [i] – interstitial volume of the solute relative to that of EDTA; and 4) (for solutes that are protein bound) cProt [i], which is the PKQuest parameter corresponding to KA i. Standard values of ecf (Table 6), frecf (=1) and cProt (KA i, Table 6) are preprogrammed in PKQuest. The only parameter that needs to be input by the user is mecf. Setting mecf =1 indicates that the solute distributes in the standard ecf space. (The default value of mecf is -1, indicating that the solutes distributes in all the tissue (interstitial and intracellular) space.) The complete PBPK pharmacokinetics of each solute is characterized by a short Maple worksheet that lists the parameters that differ from the default "standardhuman" parameters. The complete worksheets for the 11 solutes are included in the additional file "PKQuest_worksheets.doc". For all of the solutes investigated in this paper, the kidney is the major excretion pathway. The renal clearance is described two different ways in these worksheets. If the clearance is close to the free plasma glomerular clearance, than the parameter "rclr" (glomerular filtration rate) is used. For solutes that are cleared both by glomerular filtration and tubular secretion, the single parameter Tclr [kidney] is used (the rate of renal clearance of the free (unbound) plasma solute).

PKQuest is used to find the values of the adjustable parameters that provide the best fit to the experimental data by minimizing the error function:

The parameter "noise" adds a weighting factor that reduces the contribution of the low concentration data points. The default value for noise is 10% of the average concentration. (Optionally, PKQuest allows the choice of a mean square error term, and an arbitrary "noise" can be input).

Correction for antecubital vein sampling

Ideally, one would like to be able to sample arterial blood for the calculations of steady state volume of distribution and for PBPK model fitting. However, all of the data used in this paper was obtained by measuring the plasma concentrations in antecubital venous blood which differs from the arterial concentration due to the exchange with the tissues drained by the antecubital vein. A new feature has been added to PKQuest that allows one to directly correct for this effect in the PBPK model fitting. The tissue/flow distribution supplying the antecubital vein was established by applying PKQuest to the experimental data for a number of solutes for which simultaneous arterial and antecubital vein concentrations are available (Levitt, in preparation). This analysis indicated that antecubital vein blood represents approximately 10% muscle, 20% skin, 5% "other", 5% adipose and the rest (60%) A-V anastomoses. In PKQuest, specifying "arm" as the sample site outputs the antecubital vein concentration and uses this concentration to optimize the PBPK parameters using the PKQuest minimization routines [7]. This allows one to use antecubital vein blood samples when adjusting PBPK parameters for an arbitrary uncharacterized solute.

Figure 1 shows the results of this analysis for two of the extracellular solutes (DTPA and inulin) that were investigated to determine the organ composition of antecubital blood. The squares show the experimental measurements of the 99Tcm-DTPA (left column) and inulin (right column) sampled simultaneously from the artery (black) and antecubital vein (red) after a bolus venous input. The lines show the PBPK model fits for the arterial (black) and antecubital vein (red) concentrations. The difference between the arterial and antecubital vein concentrations become negligible after the first 5 to 10 minutes because the antecubital vein blood is supplied primarily by high blood flow skin and A-V anastomoses in the hand [26]. For the solutes investigated in this paper, the differences between arterial and antecubital vein concentration are very small (e.g. see fig. 2) and this correction is of minor importance. It is fortunate that this is the case. Otherwise, one could not use antecubital vein data for PBPK modeling. Unless otherwise stated, all the PBPK model data shown in the figures in this paper correspond to antecubital vein blood.

Figure 1

Arterial (black) and antecubital vein (red) concentrations of 99Tcm-DTPA (right column) and inulin (left column) after a bolus injection using the experimental data (squares) of Cousins et. al. [56]. The solid lines indicate the PBPK model predictions for the arterial (black) and antecubital vein (red) concentrations. The top row shows a plot of the absolute concentration and the bottom row is a semi-log plot of the same data. The PBPK parameters for DTPA were identical to those used for EDTA (see fig. 3).

Figure 2

Comparison of PBPK model predictions and experimental data for mannitol plotted using linear (left) or semi-log (right). The predictions for the experimental data of Laker et. al. [57] (top) and Elia et. al. [58] (bottom) are shown using the parameters in Table 6. Both the model arterial (black) and antecubital vein (red) concentration curves are shown.

Chiou [27] has emphasized that the measurement of the steady state volume of distribution (Vss) using eq. 2 is dependent on the sampling site and has shown that the value of Vss determined using antecubital vein blood will be greater than the true value of Vss using arterial blood. This effect can be quantitated for the extracellular solutes investigated in this paper using the experimental or PBPK model data shown in figure 1. Using this data, the value of Vss for the antecubital sampling site is from 4% (PBPK model) to 10% (experimental data) greater than the true arterial Vss. The values of Vss reported below are for the antecubital vein data. A small correction will be applied when these results are used to estimate the true extracellular volumes.


Experimental measurements of whole animal and individual organ volumes of distribution for extracellular solutes

Table 1 summarizes the experimental measurements in humans of the steady state volume of distribution (Vss) for mannitol, EDTA, amoxicillin, morphine-3-glucuronide, morphine-6-glucuronide and inulin. The volume of distribution is expressed in terms of total body weight (Vss, liters/Kg) and total body water (Vssw, liters/liters). The latter calculation was based on an assumption of a value of 41 liters of water for the standard 70 Kg man with 20% body fat [14]. The total body water was adjusted for body fat estimates if information was available. Table 2 lists the values of the equilibrium volume of distribution (Veq) for amoxicillin, morphine-6-glucuronide, sucrose and inulin (see Methods). Table 3 lists the Vss and the interstitial volume (VI) for a series of β-lactam antibiotics with varying amounts of plasma protein binding. VI was estimated by subtracting a plasma water volume of 2.68 liters/70 Kg from Vss. The human serum protein binding values used in Tables 1,2,3 were from the following references: amoxicillin [28, 29]; piperacillin [30]; cefatrizine [31]; ceforanide [32]; dicloxacillin [33]; flucloxacillin [29, 33]; morphine-3-glucuronide and morphine-6-glucoronide [34]. Table 4 lists the interstitial volume as a fraction of total extravascular organ water for selected solutes in rats, rabbits, cats, dogs and man.

The solutes listed in Table 4 have a rapid, possibly flow limited, exchange between the blood and interstitial space so that the interstitial space equilibrates with the blood. In contrast, albumin has a very slow rate of tissue-plasma exchange with a time constant of about 100 hours in the rat [35]. Because of this low permeability, the steady state interstitial albumin concentration in a given tissue (CI) is less than the plasma concentration and is determined by the balance between the rate of trans-capillary exchange and the rate of removal by lymph flow [36]. The standard procedure used to measure this concentration CI is to sample the lymph draining that tissue, or to use the wick method to sample the free tissue albumin. The interstitial albumin volume of distribution VI A for that tissue is then defined by:

The fourth column in Table 5 lists the steady-state albumin tissue/plasma concentration ratios (CI/CP) and the fifth column lists the corresponding ratio VI A/VI E (where VI E is the interstitial volume of EDTA) for skeletal muscle, skin and tendon. VI A is less than VI E because the interstitial matrix behaves like a size exclusion gel, restricting the albumin volume of distribution [3]. This excluded volume is larger in skin and tendon then in skeletal muscle, resulting, presumably, from varying concentrations of collagen and hyaluronan in the different tissues. The sixth column is the product of these two ratios:

This product for organ i, which will be referred to as KA i, is the relevant parameter determining the pharmacokinetics of albumin bound solutes, such as the β-lactam antibiotics (see Appendix, II).

PBPK model for extracellular solutes

The PKQuest PBPK model [7] was applied to the extracellular solutes using the new set of parameters (Table 6). A comparison of the model predictions with the experimental data is shown in figs. 2,3,4,5 for mannitol, EDTA, morphine-6-glucuronide and morphine-3-glucuronide – solutes with low or negligible protein binding. Since one of the purposes of this analysis is to find a single, unique PBPK model applicable to all extracellular solutes, the identical set of PBPK parameters has been used in all these figures. The only parameter that was varied for each solute is the renal clearance, determined using the optimization feature in PKQuest (see Methods). Figure 2 shows the PBPK model predictions for both the arterial and antecubital vein concentration (see Methods). It can be seen that the difference is small. In all other figures, only the antecubital vein concentration is shown.

Figure 3

Comparison of antecubital PBPK model predictions (line) and experimental data (squares) for EDTA [58] plotted using linear (left) or semi-log (right).

Figure 4

Comparison of PBPK linear (left) or semi-log (right) plots of model predictions and experimental data for morphine-6-glucuronide for experimental data of Lotsch et. al. [59] (top row) and Penson et. al. [42] (bottom row).

Figure 5

Comparison of PBPK linear (left) or semi-log (right) plots of model predictions and experimental data for morphine-3-glucuronide for experimental data of Penson et. al. [60].

All of the solutes shown in figs. 2,3,4,5 have been assumed to be flow limited. This is clearly not the case for inulin, which is the prototypical capillary permeability limited solute. In addition, the interstitial volume of distribution of inulin (VI) is clearly significantly smaller than that of EDTA (Table 1). In PKQuest, the parameter mecf scales the interstitial volume of all non-blood organs by the same amount. Thus, for inulin, there are 3 adjustable parameters: 1) the renal clearance ("rclr"); 2) the interstitial volume relative to that for EDTA ("mecf"); and 3) the capillary permeability for muscle ("fclear [muscle]") which sets the permeability for all the other organs (see Methods). The parameter fclear indicates the fraction of the solute that equilibrates with the tissue in one pass through the capillary. The flow-limited case corresponds to an fclear equal to 1. Figure 8 shows the PBPK model results for four different sets of experimental inulin data. The values of mecf (= 0.42) and fclear [muscle] (= 0.409, corresponding to PSmuscle = 0.61 ml/min/100 g, see eq. (6)) were first determined by optimizing the fit to the data of Odeh et. al. [37] (top row, fig. 6). A good fit (average error less than 2%) was obtained with the 3 adjustable parameters. To test the validity of this result, 3 other sets of experimental data were fitted, using these identical values of mecf and fclear [muscle] and only adjusting the renal clearance. The result of this procedure are shown in figure 6 for the experimental data of Orlando et. al. [38] (second row), Ladegaard-Pedersen et. al. [39] (third row) and Prescott et. al. [40] (bottom row). The agreement between the PBPK model and the experimental data is quite good, except for the long time data of Prescott et. al. [40]. This discrepancy is expected, since the inulin clearance in the experiments of Prescott et. al. was non-linear, decreasing by about 50% at low concentrations (long times). It is not clear why the data in the other three experiments in fig. 6 do not show this marked non-linearity.

Figure 6

Comparison of PBPK linear (left) or semi-log (right) plots of model predictions and experimental data for inulin. The parameters rclr, mecf, and fclear [muscle] were adjusted to give the optimal fit to the data of Odeh et. al. [37] (top row). These values of mecf and fclear [muscle] were then used to fit the experimental data of Orlando et. al. [38] (second row), Ladegaard-Pedersen et. al. [39] (third row) and Prescott et. al. [40] (bottom row), varying only the renal clearance to obtain the best fit.

Figures 7,8,9,10,11,12 compare the PBPK model predictions with experimental data for the six β-lactam antibiotics listed in Table 3 (amoxicillin, piperacillin, cefatrizine, ceforanide, flucloxacillin and dicloxacillin). The influence of the protein (albumin) binding on the pharmacokinetics of these solutes is characterized by the default PKQuest values for KA i in the different organs (Table 6, see Appendix, II and III) and the experimental values for the fraction bound in plasma (fp, Table 1). The effect of including a capillary permeability limitation was investigated for these solutes because, as discussed above (see eq. (6)), one would predict that solutes with a high amount of protein binding should be capillary permeability limited. (All the other parameters were identical to those used in figs. 2,3,4,5). The addition of a capillary permeability limitation provides a significant improvement in the fit between the PBPK model and the experimental data for the two antibiotics with the highest degree of binding, dicloxacillin (97% bound, fig. 12) and flucloxacillin (93% bound, fig. 11). For ceforanide (fig. 10), which is 82% bound, there is a slight improvement in the fit with the addition of a small permeability limitation. The value of the PS product for these three permeability limited solutes is 11.3 ml/min/100 g (fclear= 0.25) for dicloxacillin, 7.0 ml/min/100 g (fclear = 0.344) for flucloxacillin and 5 ml/min/100 g (fclear = 0.52) for ceforanide. For the other antibiotics (amoxicillin, piperacillin and cefatrizine) with protein binding of 62% or less, the addition of a capillary permeability limitation does not significantly improve the fit. For these flow limited antibiotics (amoxicillin, piperacillin, cefatrizine), there is only 1 adjustable parameter – the renal clearance. For the capillary permeability limited solutes (ceforanide, flucloxacillin, dicloxacillin) there are two adjustable parameters, the renal clearance and fclear [muscle].

Figure 7

Comparison of PBPK linear (left) or semi-log (right) plots of model predictions and experimental data for amoxicilllin for experimental data of Sjovall et. al. [61] (top row) and Arancibia et. al. [62] (bottom row). The flow limited PBPK model is used (no capillary permeability limitation).

Figure 8

Comparison of PBPK linear (left) or semi-log (right) plots of model predictions and experimental data for piperacillin from Lode et. al. [30]. The flow limited PBPK model is used (no capillary permeability limitation).

Figure 9

Comparison of PBPK linear (left) or semi-log (right) plots of model predictions and experimental data for cefatrizine from Pfeffer et. al. [31]. The top row is for the flow limited model (fclear = 1) with an average error between PBPK model and experimental data of 5.4%. The addition of a small capillary permeability limitation does not significantly improve the fit (not shown).

Figure 10

Comparison of PBPK linear (left) or semi-log (right) plots of model predictions and experimental data for ceforanide from Pfeffer et. al. [32]. The top row is the optimal fit for the case of a capillary permeability limitation (fclear [muscle] = 0.518, PSmuscle = 5 ml/min/100 gm, average error = 4.3%). The bottom row is the best fit for the flow limited case (fclear = 1, average error = 8.6%).

Figure 11

Comparison of PBPK linear (left) or semi-log (right) plots of model predictions and experimental data for flucloxacillin from Wise et. al.[29]. The top row is the optimal fit for the case of a capillary permeability limitation (fclear [muscle] = 0.344, PSmuscle = 7.0 ml/min/100 gm, average error = 6.4%). The bottom row is the best fit for the flow limited case (fclear = 1, average error = 11%).

Figure 12

Comparison of PBPK linear (left) or semi-log (right) plots of model predictions and experimental data for dicloxacillin from Lofgren et. al. [63]. The top row is the optimal fit for the case of a capillary permeability limitation (fclear [muscle] = 0.25, PSmuscle = 11 ml/min/100 gm, average error = 12%. The bottom row is the best fit for the flow limited case (fclear = 1, average error = 18%).

This new PBPK parameter set (Table 6) also provided good fits (not shown) to the other solutes that have been previously investigated with PKQuest (propranolol [7], D2O and ethanol [6], anesthetic gases and toluene [5]). The slightly modified Maple worksheets describing these solutes are available on the PKQuest web site


Volume of distribution of selected extracellular solutes in humans

Table 1 summarizes the experimental measurements in humans of the steady state volume of distribution (Vss) for mannitol, EDTA, amoxicillin, morphine-3-glucuronide, morphine-6-glucuronide and inulin. In order to be included in this table, solutes had to meet the following criteria: 1) extracellular distribution; 2) low level of plasma protein binding, 3) no evidence of any unusual tissue binding or accumulation; 4) major excretory pathway is renal; 5) linear pharmacokinetics; and 6) published data describes the plasma time course in humans in enough detail that the value of Vss could be recalculated using PKQuest (see Methods). These 6 conditions are quite restrictive and most other drug classes do not meet all of these conditions. For example, although the aminoglycosides (e.g. gentamicin) satisfy conditions 1, 2, 4 and 5 and 6, they are known to be actively concentrated in the renal cortex and in several other organs [41]. There is even a suggestion that inulin, the classic extracellular solute, may not satisfy condition 5 because of non-linear renal clearance [40].

There are two ambiguities in this calculation of Vss. First, as discussed above in the Methods section, the use of antecubital vein samples should lead to a value of Vss that is about 5 to 10% greater than the true value. Second, the calculation of Vss using eq. (2) depends on an accurate extrapolation of the venous concentration to long times to determine the integrated mean residence time. Errors in the long time, low concentration venous concentrations can produce large errors in the extrapolated MRT. The procedure used in Table 1 involves using deconvolution to find a multi-exponential unit bolus response function, and then integrating this function from 0 to infinity (see eq. 2). The values of Vss were calculated using either a 2 or 3 (if there were enough experimental data points) exponential function. Some indication of the accuracy of this extrapolation can be obtained by comparing the results for the 2 and 3 exponential functions, listed in Table 1. The calculation with the largest difference between the Vss calculated for the 2 and 3-exponential function is for the morphine-6-glucuronide data of Penson et. al. [42]. Figure 13 compares the semi-log fits for the 2 (black line) and 3-exponential (red line) functions to the experimental venous data. In the figure on the left, which displays the two fits over the time of the experimental measurements (600 minutes), the fits look almost identical. However, as can be seen in the right figure, which extends the fits out to 2000 minutes, the 3-exponential fit has a markedly different extrapolation, producing the 33% greater value of Vss for the 3-exponential fit (Table 1). As this figure illustrate, the value of Vss will be underestimated if there is a compartment with a very long time constant (e.g. 690 minutes for the 3-exponential fit, fig. 1). For each solute in Table 1, the 2-exponential fit was superior using the standard statistical tests [43], and this is the value that was used to determine Vssw (Table 1) and was used in this paper.

Figure 13

Comparison of morphine-6-glucuronide experimental data [42] with venous concentration predicted using either a 2 (black) or 3 (red) exponential bolus response function.

Both of these ambiguities can be avoided if one can measure the true equilibrium volume of distribution (Veq). One approach that has been used in humans is to first establish a steady state concentration using a long, constant IV infusion. The equilibrium volume of distribution (Veq) is then determined using eq. 3, where Ceq is the steady state plasma concentration, and Atot is the total amount of solute in the animal. If the solute is not metabolized and the kidney is the only route of excretion, Atot can be determined by quantitative urine collection after stopping the infusion. The measurements of Schwartz et. al. [44] of the Veq of inulin represent the most detailed investigation of this type in humans (Table 2). Schwartz et. al. [44] infused inulin for periods up to 30 hours and found that equilibration was at least 95% complete by 5 hours. This indicates that there are no compartments with significant volumes of distribution with time constants longer than 300 minutes. There is some evidence that the renal clearance of inulin is non-linear [40] which, if correct, could produce an error in the Vss calculation for inulin in Table 1. However, this non-linearity would have no effect on the Veq calculations in Table 2. Deane et. al. [45] used this constant infusion approach to compare the Veq of inulin and sucrose and reported that the sucrose Veq is slightly smaller than that of inulin (Table 2). This result is not consistent with the other measurement of Veq or Vss (Tables 1 and 2), which indicate that solutes of the size of sucrose should have significantly larger volumes of distribution than inulin. One explanation for this inconsistency is that there is significant extrarenal clearance of sucrose so that the value of Atot (determine by 24 hour urine collection) is underestimated.

For animal studies, the most direct measurement of Veq for solutes with no extrarenal clearance is to nephrectomize the animal and determine the equilibrium concentration after a bolus IV injection. Although this is not an option for humans, a good approximation to this situation can be obtained in subjects with renal failure and very low rates of clearance. For this case, Veq can be estimated from the coefficient of the slow exponential term in the bolus response function (see Methods). Figure 14 shows the good agreement between the predicted venous concentration using a 2-exponential bolus response function and the experimental data of Arancibia et. al. [46] for amoxicillin and Hanna et. al. [47] for morphine-6-glucuronide. Because of the low clearance rate in these subjects, measurements can be extended for longer times (720 minutes for morphine-6-glucuronide and 1440 minutes for amoxicillin) then those used in Table 1. The value of Veq and the time constants for the slow and fast compartments are listed in Table 2. For both solutes, the time constant for the slow compartment is about 100 times longer that that for the fast compartment. The time constant of the slow compartment (909 and 2083 minutes) reflects the slow rate of systemic clearance in these subjects. If there were a slowly exchanging compartment that was missed in these subjects, it must have a time constant longer than 900 minutes, which seems unlikely. The values of the "equilibrium" volume of distribution listed in Table 2 are consistent with the values of Vss listed in Table 1 which, as discussed above, should be 5 to 10% too large because of antecubital vein sampling.

Figure 14

Comparison of experimental data with venous concentration predicted using a 2-exponential bolus response function for amoxicillin (left) and morphine-6-glucuronide (right) in subjects with renal failure.

Figure 15

Plot of VI i/VI E versus fp i for the β-lactam antibiotics listed in Table 3. The red line, which is the least squares linear fit to the data, has a slope of 0.83 and intercept of 0.28.

The results in Tables 1 and 2 indicate that the extracellular volume of distribution is surprisingly large, about 42% of the total body water. The distribution of this space among the different organs will be discussed in the next section. Tables 1 and 2 also show that the volume of distribution of inulin is about one half the value of the other solutes. This difference is supported by the results of Ladegaard-Pedersen and Engell [48] who found that the inulin volume was, on average, 48% smaller than the EDTA volume in a series of 22 patients. It is ironic that inulin, which is usually regarded as the classical, standard, extracellular marker, seems to be a poor representative of the extracellular space for the typical pharmacokinetic agent. The reduced volume for inulin is presumably produced by the "excluded volume" associated with the interstitial matrix [3], similar to the excluded volume of albumin.

Tissue distribution of interstitial volume (VI)

Table 4 summarizes the measurements in the literature of the interstitial volumes (VI) of different organs for a series of extracellular solutes. The value of VI is expressed as the fraction of the total extravascular organ water. The most striking aspect of the results in this table is the marked variation in VI for different tissues – varying from about 10 to 15% of total tissue water for skeletal muscle, to 50 to 70% for skin to 80% for tendon. (The high value of 23% for VI obtained by Bell et. al. [49] in the dog is probably associated with the manipulations required to collect the lymph [50]). The measurements in the dog by Nichols et. al. [51] in 1953 (Table 4) represent the first detailed investigation of the tissue distribution of interstitial volume and form the basis of the current understanding of the extracellular space. These early measurements of Nichols et. al. [51] emphasized the important contribution of connective tissue (represented by tendon) to the total interstitial kinetics [51, 52]. In general, the more recent results in Table 4 are in agreement with these earlier measurements.

Volume of distribution of protein bound β-lactam antibiotics

The above discussion has assumed that the test solutes distribute freely in the extracellular blood and tissue water without any binding so that, at equilibrium, the interstitial and plasma concentrations are equal. This is not true if there is blood and/or tissue binding. In this case, the volume of distribution (see eq. (3)) will depend on the degree of tissue binding relative to plasma binding. For the general case, the binding in the different tissues will be unique for each solute, requiring direct tissue binding measurements. However, for the case where the binding is linear and is only due to the plasma and tissue albumin, one can derive (see Appendix, II) the following general relationship for the ratio of the interstitial volume of distribution of solute i (VI i) relative to the interstitial EDTA volume of distribution (VI E) as a function of fP i, the fraction of solute i that is free in plasma:

The parameter KA is the ratio of the albumin concentration in the EDTA interstitial space to the plasma albumin concentration (eq. (9)). The KA in this equation is for the whole animal and is the weighted average of the KA i (see eq. (5)) for the individual organs (see Tables 5 and 6). The parameter αi is the ratio of the interstitial space for solute i relative to EDTA in the absence of binding. Since EDTA and the β-lactam antibiotics are of similar size, αi should be close to 1. In the limit of no protein binding (fP = 1), VI i/VI E = αi. In the limit of very strong binding (fP = 0), VI i/VI E = KA (the ratio of tissue to plasma albumin).

Table 3 lists the Vss and VI of a series of β-lactam antibiotics with varying amounts of plasma protein binding. These antibiotics should satisfy the assumptions required for eq. (10) since they all have linear kinetics and the protein binding results primarily from binding to albumin. Figure 15 shows a plot of VI i/VI E versus fp i for these solutes. (The Vss for flucloxacillin in Table 3 was scaled for an amoxicillin value of 0.26). According to eq. (10), the data points in this plot should fall on a straight line, with a slope of (αi-KA) and intercept of KA. The line in fig. 15, which is the least squares fit to the data points, corresponds to a KA of 0.28 and an αi of 1.1. Comparing this value of KA, which represents a weighted average for the whole animal (eq, (5)), to the value of KA i for the individual organs (Table 5) suggest that organs such as skin, tendon and connective tissue dominate the contribution to the total KA.

PBPK model parameters for extracellular solutes

Based on the physiological data in Tables 1,2,3,4,5, a new set of PBPK parameters was developed that should provide a more accurate description of extracellular solutes, while remaining valid for the other solutes modeled previously by PKQuest. Table 6 (see Methods) summarizes the new PBPK data set. The comparison of the PBPK model predictions using this new parameter set with the experimental data is shown in figures 2,3,4,5,6,7,8,9,10,11,12 for the 11 solutes listed in Tables 1 and 3.

Distributing the extracellular water among the different organs presents something of a conundrum. If one uses the data from Table 4 for the organs that are usually regarded as making up most of the non-fat, non-bone body weight (blood, muscle, skin, gastrointestinal (portal) and liver), one can only account for about 7.3 liters of extracellular water, less than half the total value of about 17.3 liters based on the VSS data in Table 1 (extracellular water equal to 42% of total body water). Since the only other organs of any significant weight are the adipose and connective tissue, these organs must represent a large fraction of the total extracellular water. Connective tissue is very heterogeneous, – consisting of tendons, cartilage, subcutaneous, etc. It is poorly characterized and values for its weight, composition, and blood flow represent, at best, crude estimates. In the classic reference on the extracellular space, Edelman and Leibman [1] estimated that the connective tissue had a mass of 6 Kg and was 70% water (4.2 liters). The fraction of adipose tissue that is extracellular water varies for different locations. One component of the adipose tissue is the subcutaneous component of skin ("subcutis"). In the measurements of Reed et. al. [50] in the rat, subcutis is about 70% water, of which about 51% is extracellular (see Table 4). Adipose tissue in other locations can have much smaller extracellular components. If one assumes that, on average, adipose tissues is 70% lipid and 30% water, then, for the standard human (70 Kg, 20% fat) adipose tissue will contain 6 liters of water.

Based on the above estimates, adipose and connective tissue contain about 10.2 liters of water. For the PBPK model in Table 6, it is assumed that most of this water is extracellular and it is distributed between "adipose" tissue and two connective tissue organs: "tendons" with a relatively low water fraction and low blood flow, and "other" with a higher water fraction and higher blood flow. The total extracellular volume for this PBPK data set is 17.45 liters, or 42% of the total water (Table 6).

As described above (see eq. (10), the interstitial volume of distribution of the β-lactam antibiotics depends on the parameter KA, the ratio of the albumin concentration in the EDTA interstitial space to the plasma albumin concentration. This value represents the sum over all the different organs. Each individual organ i is characterized by a value of KA i. The assumed PBPK values of KA i for the different organs are also listed in Table 6. It is assumed that the organs "skin" and "other" have a smaller value of KA i then the other organs, consistent with the experimental data in Table 5. The organ-weighted total value of KA is 0.28 (see Table 6 and eq. (5)), identical to the value obtained by fitting the data in Table 3 (fig. 15).

Capillary Permeability Limitation

Although nearly all PBPK models assume that the blood-tissue exchange for solutes such as EDTA or the β-lactam antibiotics is "flow-limited" (i.e., infinite capillary permeability), there is only indirect evidence to support this assumption. The standard approach that is used to experimentally determine the capillary permeability-surface area product (PS) is the organ perfused early extraction (E0) method [25]. For an ideal, homogeneous organ, E0 is equal to the PKQuest parameter "fclear" (see, eq. (6)), the fraction of solute that equilibrates with the tissue in one pass through the capillary [4]. A flow-limited solute has an fclear of 1. As shown in eq. (6), fclear is a function of 3 different physiological parameters: 1) the capillary permeability-surface area product (PS); 2) the organ plasma flow (F); and 3) the fraction of the solute that is free in the plasma (fP). The standard reference value for the PS of EDTA in skeletal muscle, skin and subcutaneous tissue is about 5 ml/min/100 g [25]. Substituting this value for PS in eq. (6) and using the resting muscle blood flow of 2.25 ml/min/100 gm (Table 6) yields an fclear of 0.89. That is, there is 89% equilibration of EDTA in one pass through the capillary. Since it is difficult to experimentally distinguish an fclear of 0.89 from an fclear of 1 – this result is consistent with the flow limited assumption for EDTA.

There are several uncertainties in this estimate of fclear. To accurately measure PS with the E0 method it is necessary to use relatively high rates of perfused organ flows so that E0 (and fclear) are much less than 1. In addition, the E0 method tends to underestimate the true value of fclear because of factors such as back diffusion and capillary heterogeneity [53, 54]. This underestimate is greater for solutes with relatively high permeability, such as EDTA. The two most recent measurements of the skeletal muscle PS of EDTA are 5.0 ml/min/100 g in the cat [55] and 12.9 in rat [54]. Both of these studies have emphasized the importance of using very high blood flows to eliminate the effects of heterogeneity and back diffusion. However, the PS at these high flows (at least 10 times resting flows) may not be relevant to the value of PS at the resting muscle blood flows that are required for the human pharmacokinetic studies. In particular, one would expect that the PS product would be greater at the high flows because the value of S (capillary surface area) would be increased as a result of recruitment of additional capillaries. In addition, the capillaries that are recruited at high flows may not be representative of the low flow capillaries.

The use of the β-lactam antibiotics provides a new approach to the measurement of capillary permeability at resting blood flows [4]. Since the factor fP multiplies the exponential term in eq. (6), a high protein binding (small fP) is kinetically equivalent to increasing the blood flow. For example, flucloxacillin, which is 93% bound, has an fP of 0.07 and a corresponding fclear of 0.15 for a PS of 5 ml/min/100 g and resting muscle blood flow. That is, at resting blood flows only 15% of the plasma flucloxacillin will equilibrate with the tissue in one pass – a large capillary permeability limitation. The use of these protein bound antibiotics allows one to directly measure the capillary permeability at low flows without the uncertainty introduced by heterogeneity and back diffusion.

The approach used in this paper to establish if the antibiotics are flow limited is to determine whether the agreement between the PBPK model and the experimental data was improved by inputting an fclear [muscle] that is less than 1, i.e. a capillary permeability limitation. This is an indirect approach because one would expect that the fit should improve just because another adjustable parameter is introduced. Despite this proviso, there is little doubt that the two antibiotics with the highest degree of protein binding, dicloxacillin (fig. 14, fP = 0.03; fclear = 0.25) and flucloxacillin (fig. 13, fP = 0.07; fclear = 0.344) have a significant capillary permeability limitation. Comparing the PBPK model output for the flow limited and permeability limited conditions in figs. 13 and 14, it can be seen that the primary effect of the permeability limitation is at early times, increasing the initial peak in the plasma concentration. For ceforanide, the third highest binding antibiotic in Table 3 (fP = 0.18, fig. 12), there is a slight improvement in the fit using a fclear [muscle] of 0.52. For the 3 other antibiotics in Table 3, no significant improvement in the fit was obtained by the addition of an fclear, suggesting that they are flow limited. This is consistent with eq. (6) since these antibiotics have an fP of 0.38 or greater and should have an fclear of greater than 0.7 (assuming a PS = 7, the value obtained for flucloxacillin).

The PBPK estimates of PS of muscle for the β-lactam antibiotics PS are 11.3 ml/min/100 g for dicloxacillin, 7 for flucloxacillin and 5 for ceforanide. Since these three antibiotics have similar structures, one would expect them to have similar PS values. This suggests that the difference in these three PS values provides an indication of the uncertainty in the approach of using PBPK fitting to determine PS. These values are somewhat larger than the recent estimate of Watson [55] for cat skeletal muscle PS of 5 and 7.6 ml/min/100 g for EDTA and mannitol, respectively. These values of Watson [55] were obtained at very high isoproterenol induced blood flows (more than 10 times resting flows) and one would expect these PS values should be greater than the resting PBPK values because of the recruitment of new capillaries and the accompanying increase in capillary surface area. The PBPK value estimated in this paper for the PS of inulin in skeletal muscle is 0.61 ml/min/100 g (see fig. 6) which is in the range of values reported in the literature (0.59 to 1.62 [25]).

The capillary permeability measurements reported here using the PBPK approach are limited by fact that they are whole animal studies. The pharmacokinetics for these extracellular solutes and, therefore, the capillary permeability estimates, are dominated by the organs "skin", "tendon" and "other" which represent 62% of the extravascular water (Table 6). In the implementation in PKQuest, the parameter fclear [muscle] is input, and the fclear of all the other organs are then determined using eq. (6) and a preprogrammed value of PS of each organ relative to muscle. It has been assumed in PKQuest that "skin", "tendon" and "other" have the same PS as skeletal muscle (see Methods).

Despite the quantitative uncertainty in these PS measurements, the qualitative implications are clear. For solutes such as EDTA, mannitol, morphine-6-glucluronide and the weakly bound β-lactam antibiotics, adding the additional factor of a capillary permeability limitation does not significantly improve the fit between the experimental data and the PBPK model prediction. That is, a flow-limited model adequately describes these solutes. In contrast, as predicted using eq. (6), a capillary permeability limitation significantly improves the fit for the highly protein bound β-lactam antibiotics.


In general, each solute that is modeled by PBPK requires unique information about the binding or partition in the individual organs – information that can only by obtained by measurements in animals and then extrapolated to humans. However, this step is not necessary for the extracellular solutes investigated here. Since the solutes do not enter cells, all the diverse cellular binding and uptake processes are circumvented. In addition, if the interstitial binding is primarily due to albumin, as in the case of the β-lactam antibiotics, then the interstitial binding can be predicted just from information about the interstitial albumin concentration (KA, Table 6). Thus, one can develop a PBPK model or parameter set (Table 6) that is generally applicable to this class of drugs.

This approach is illustrated here for the 11 solutes listed in Tables 1 and 3. A comparison of the model predictions with the experimental data is shown in figs. 2,3,4,5,6,7,8,9,10,11,12. The same PBPK parameter set (Table 6) was used for all of the solutes. The PBPK model predictions shown in these figures were based on the use of only 1 (renal clearance) or 2 (skeletal muscle PS for drugs with high binding) adjustable parameters. In principal, the PS can also be predicted using eq. (6) if an accurate value of fP is available. Thus, if one knows the systemic clearance for these extracellular solutes, it should be possible to predict the complete time course of the absolute drug concentrations in all the major organs.


I. Volume of distribution of 2-compartment model for case where clearance is slow compared to compartment exchange

A two compartment system with a clearance rate constant of k10 from the central compartment and exchanges rates of k12 and k21 between the central and peripheral compartments has the following kinetic solution for the central (venous) concentration c(t) in response to a unit bolus input (V1 is the volume of the central compartment, D is the bolus dose):

For the special case where the clearance rate constant is small compared to the other exchange rates (k10 << k12 + k21):

One can define an "equilibrium volume of distribution" (Veq) for the case where the clearance rate constant (k10) is set to zero and the exchange between the peripheral and central compartment is in equilibrium (A1 and A2 are the amounts in the two compartments):

Substituting this relation for Veq into eq. (A.12) for B2:

(A.14)     V eq D / B 2

This expression for Veq is identical to volume of distribution that Wagner [11] defined as Vdext:

(A.15)     V eq V dext = D / B β

where Bβ is the coefficient of the slow, terminal exponential, determined by extrapolating the terminal component back to time 0, and Vdext is the "extrapolated volume of distribution".

II. Derivation of expression for steady state volume of distribution for the case of blood and tissue albumin binding


VD Steady state volume of distribution

VP Plasma volume

VE Interstitial space for EDTA.

Vi Interstitial space for solute i.

λE, λi Fraction of VE and Vi that albumin distributes in (λ ≤ 1).

αi = Vi/VEi ≥ λE).

CP Total (free + bound) solute concentration in plasma.

cP, c Free concentration in plasma and tissue.

fP, fi, f0 Fraction of total solute that is free in plasma and tissue and at 0 concentration.

k Albumin binding constant for solute i.

BP Concentration of protein binding sites in plasma.

BI Concentration of protein binding sites in interstitial protein space.

β = BI/BP.

KA = βλE

The steady state volume of distribution of solute i (VD) is define by:

(A.16)     Amount = V D C P = V P C P + V 1 C 1 + V 2 C 2

where Amount is the total amount of solute in the organ, CP,C1 and C2 are the total (free plus bound) solute concentration in plasma space VP, in the interstitial space V1 that albumin distributes in and in the remainder of the interstitial space V2 available to the solute but not to albumin. Using the above definitions of λi, V1 and V2 can be related to the interstitial space of solute i (Vi):

(A.17)     V 1 = λ i V i     V 2 = (1 - λ i ) V i

Since the total concentration C is related to the free concentration c by C = c/f, and at steady state, c has the same value in all three compartments (assuming removal only from the plasma compartment):

where fp and fi are the fraction of total solute that is free in the plasma and interstitial space, respectively, and it has been assumed there is no protein binding in V2. Solving (A.18) for VD:

The fraction free (f) is related to the protein binding constant (k) and the concentration of protein binding sites (B) by the Scatchard equation:

For the case of linear protein binding (no saturation), kc << 1 and δ (eq. (A.19)) is equal to:

(A.21)     δ = f p (1 + λkB I )     (Linear)

where BI is the interstitial protein concentration. For the linear case, the binding constant k in eq. (A.21) can be expressed in terms of fp and Bp (plasma protein concentration):

(A.22)     δ = βλ i + f P (1 - βλ i )     (Linear)

where β is the ratio BI/Bp. Finally, since the parameters λi and Vi depend on the specific solute, they need to be related to some standard value. In PKQuest, the standard is EDTA (λE, VE):

Also, it is assumed that Vi ≥ Valb, and, therefore, αi ≥ λE. Substituting these relations into eqs. (A.22) and (A.19) :

(A.24)     V D = V P + V E [βλ E + f p i - βλ E )]     (Linear)

The interstitial volume of distribution of the bound solute VI (= VD-VP) can then be expressed relative to VE, the interstitial volume of the EDTA (assumed to have no protein binding):

(A.25)     V I / V E = (α i - K A ) f p + K A     K A = βλ E

where KA is the albumin concentration in the EDTA interstitial space divided by the plasma albumin concentration.

III. Modification of differential equations describing PBPK model if the interstitial protein volume is less than the interstitial volume of solute i

In the original derivation for PKQuest [7], it was assumed that the solute and the protein (albumin) distributed in the same interstitial volume. If the protein volume is less than the solute volume, then the differential equation describing the change in total tissue solute concentration (eq. 17 in [7]) must be modified. The total interstitial amount (Ai) of solute i is described by the second two terms on the right side of eq. (A.18):

where ci and cb are the free and bound tissue concentrations of solute i and eq. (A.20) has been substituted for fi. Comparing the expression for cb in eq. (A.26) with the original expression (eq. 16 in [7]), it can be seen that the only modification that is required is the substitution of λiBI in place of BI.

In the implementation in PKQuest the interstitial binding in organ j is characterized by the following set of parameters:

ecf [j]: interstitial EDTA volume as fraction of total organ water.

frecf [j,i]: = αi = interstitial space of solute i as fraction of EDTA space

cProt [j]: = KA j = βλE

λiBI = cProt*Bp /frecf


  1. 1.

    Edelman IS, Liebman J: Anatomy of body water and electrolytes. American Journal of Medicine. 1959, 27: 256-275.

  2. 2.

    Aukland K, Nicolaysen G: Interstitial fluid volume: local regulatory mechanisms. Physiol Rev. 1981, 61: 556-643.

  3. 3.

    Bert JL, Pearce RH: The interstitium and microvascular exchange. In: Hanbook of Physiology, The Cardiovascular System Microcirculation, Part I. Edited by: Renkin EM, Michel CC. 1984, Bethesda, MD: American Physiological Society, IV: 521-547.

  4. 4.

    Levitt DG: PKQuest: capillary permeability limitation and plasma protein binding – application to human inulin, dicloxacillin and ceftriaxone pharmacokinetics. BMC Clin Pharmacol. 2002, 2: 7-10.1186/1472-6904-2-7.

  5. 5.

    Levitt DG: PKQuest: volatile solutes – application to enflurane, nitrous oxide, halothane, methoxyflurane and toluene pharmacokinetics. BMC Anesthesiol. 2002, 2: 5-10.1186/1471-2253-2-5.

  6. 6.

    Levitt DG: PKQuest: measurement of intestinal absorption and first pass metabolism – application to human ethanol pharmacokinetics. BMC Clin Pharmacol. 2002, 2: 4-10.1186/1472-6904-2-4.

  7. 7.

    Levitt DG: PKQuest: a general physiologically based pharmacokinetic model. Introduction and application to propranolol. BMC Clin Pharmacol. 2002, 2: 5-10.1186/1472-6904-2-5.

  8. 8.

    Levitt DG: The use of a physiologically based pharmacokinetic model to evaluate deconvolution measurements of systemic absorption. BMC Clin Pharmacol. 2003, 3: 1-10.1186/1472-6904-3-1.

  9. 9.

    Tsuji A, Yoshikawa T, Nishide K, Minami H, Kimura M, Nakashima E, Terasaki T, Miyamoto E, Nightingale CH, Yamana T: Physiologically based pharmacokinetic model for beta-lactam antibiotics I: Tissue distribution and elimination in rats. J Pharm Sci. 1983, 72: 1239-52.

  10. 10.

    Rominger KL, Wolf M: A critical view of distribution volumes. In: Pharmacokinetics during drug development: data analysis and evaluation techniques. Edited by: Bozler G, van Rossum JM. 1982, Stuttgart: Gustav Fischer Verlag, 45-64.

  11. 11.

    Wagner JG: Linear pharmacokinetic equations allowing direct calculation of many needed pharmacokinetic parameters from the coefficients and exponents of polyexponential equations which have been fitted to the data. J Pharmacokinet Biopharm. 1976, 4: 443-67.

  12. 12.

    Snyder WS, ed: Report of the task group on refereance man. 1975, Oxford: Pergamon Press

  13. 13.

    Jonsson F, Johanson G: Physiologically Based Modeling of the Inhalation Kinetics of Styrene in Humans Using a Bayesian Population Approach. Toxicol Appl Pharmacol. 2002, 179: 35-49. 10.1006/taap.2001.9331.

  14. 14.

    Chumlea WC, Guo SS, Zeller CM, Reo NV, Baumgartner RN, Garry PJ, Wang J, Pierson RN, Heymsfield SB, Siervogel RM: Total body water reference values and prediction equations for adults. Kidney Int. 2001, 59: 2250-8. 10.1046/j.1523-1755.2001.0590062250.x.

  15. 15.

    Williams LR, Leggett RW: Reference values for resting blood flow to organs of man. Clin Phys Physiol Meas. 1989, 10: 187-217. 10.1088/0143-0815/10/3/001.

  16. 16.

    Schloerb PR, Friis-Hansen BJ, Edelman IS, Solomon AK, Moore FD: The measurement of total body water in the human subject by deuterium oxide dilution. J Clin Invest. 1950, 29: 1296-310.

  17. 17.

    Laine H, Yki-Jarvinen H, Kirvela O, Tolvanen T, Raitakari M, Solin O, Haaparanta M, Knuuti J, Nuutila P: Insulin resistance of glucose uptake in skeletal muscle cannot be ameliorated by enhancing endothelium-dependent blood flow in obesity. J Clin Invest. 1998, 101: 1156-62.

  18. 18.

    Eriksson JW, Smith U, Waagstein F, Wysocki M, Jansson PA: Glucose turnover and adipose tissue lipolysis are insulin-resistant in healthy relatives of type 2 diabetes patients: is cellular insulin resistance a secondary phenomenon?. Diabetes. 1999, 48: 1572-8.

  19. 19.

    Elia M, Kurpad A: What is the blood flow to resting human muscle?. Clin Sci (Lond). 1993, 84: 559-63.

  20. 20.

    Munson ES, Eger EI, Tham MK, Embro WJ: Increase in anesthetic uptake, excretion, and blood solubility in man after eating. Anesth Analg. 1978, 57: 224-31.

  21. 21.

    Virtanen KA, Lonnroth P, Parkkola R, Peltoniemi P, Asola M, Viljanen T, Tolvanen T, Knuuti J, Ronnemaa T, Huupponen R, Nuutila P: Glucose uptake and perfusion in subcutaneous and visceral adipose tissue during insulin stimulation in nonobese and obese humans. J Clin Endocrinol Metab. 2002, 87: 3902-10. 10.1210/jc.87.8.3902.

  22. 22.

    Tothill P, Hooper G: Measurement of tendon blood flow in rabbits by microsphere uptake and 133Xe washout. J Hand Surg [Br]. 1985, 10: 17-20.

  23. 23.

    Riggi K, Wood MB, Ilstrup DM: Dose-dependent variations in blood flow evaluation of canine nerve, nerve graft, tendon, and ligament tissue by the radiolabeled-microsphere technique. J Orthop Res. 1990, 8: 909-16.

  24. 24.

    Grimby G, Nilsson NJ, Sanne H: Serial determinations of cardiac output at rest. Br Heart J. 1966, 28: 118-21.

  25. 25.

    Crone C, Levitt DG: Capillary permeability to small solutes. 1984, Bethesda, Maryland: American Physiological Society

  26. 26.

    Lossius K, Eriksen M, Walloe L: Fluctuations in blood flow to acral skin in humans: connection with heart rate and blood pressure variability. J Physiol. 1993, 460: 641-55.

  27. 27.

    Chiou WL: The phenomenon and rationale of marked dependence of drug concentration on blood sampling site. Implications in pharmacokinetics, pharmacodynamics, toxicology and therapeutics (Part I). Clin Pharmacokinet. 1989, 17: 175-99.

  28. 28.

    Sutherland R, Croydon EA, Rolinson GN: Amoxycillin: a new semi-synthetic penicillin. Br Med J. 1972, 3: 13-6.

  29. 29.

    Wise R, Gillett AP, Cadge B, Durham SR, Baker S: The influence of protein binding upon tissue fluid levels of six beta-lactam antibiotics. J Infect Dis. 1980, 142: 77-82.

  30. 30.

    Lode H, Elvers A, Koeppe P, Borner K: Comparative pharmacokinetics of apalcillin and piperacillin. Antimicrob Agents Chemother. 1984, 25: 105-8.

  31. 31.

    Pfeffer M, Gaver RC, Ximenez J: Human intravenous pharmacokinetics and absolute oral bioavailability of cefatrizine. Antimicrob Agents Chemother. 1983, 24: 915-20.

  32. 32.

    Pfeffer M, Gaver RC, Van Harken DR: Human pharmacokinetics of a new braod-spectrum parenteral cephalosporin antibiotic, ceforanide. J Pharm Sci. 1980, 69: 398-403.

  33. 33.

    Roder BL, Frimodt-Moller N, Espersen F, Rasmussen SN: Dicloxacillin and flucloxacillin: pharmacokinetics, protein binding and serum bactericidal titers in healthy subjects after oral administration. Infection. 1995, 23: 107-12.

  34. 34.

    Van Crugten JT, Sallustio BC, Nation RL, Somogyi AA: Renal tubular transport of morphine, morphine-6-glucuronide, and morphine-3-glucuronide in the isolated perfused rat kidney. Drug Metab Dispos. 1991, 19: 1087-92.

  35. 35.

    Wiig H, Kaysen GA, al-Bander HA, De Carlo M, Sibley L, Renkin EM: Interstitial exclusion of IgG in rat tissues estimated by continuous infusion. Am J Physiol. 1994, 266: H212-9.

  36. 36.

    Johnson JA: Capillary permeability, extracellular space estimation, and lymph flow. Am J Physiol. 1966, 211: 1261-3.

  37. 37.

    Odeh YK, Wang Z, Ruo TI, Wang T, Frederiksen MC, Pospisil PA, Atkinson AJ: Simultaneous analysis of inulin and 15N2-urea kinetics in humans. Clin Pharmacol Ther. 1993, 53: 419-25.

  38. 38.

    Orlando R, Floreani M, Padrini R, Palatini P: Determination of inulin clearance by bolus intravenous injection in healthy subjects and ascitic patients: equivalence of systemic and renal clearances as glomerular filtration markers. Br J Clin Pharmacol. 1998, 46: 605-9. 10.1046/j.1365-2125.1998.00824.x.

  39. 39.

    Ladegaard-Pedersen HJ: Measurement of extracellular volume and renal clearance by a single injection of inulin. Scand J Clin Lab Invest. 1972, 29: 145-53.

  40. 40.

    Prescott LF, McAuslane JA, Freestone S: The concentration-dependent disposition and kinetics of inulin. Eur J Clin Pharmacol. 1991, 40: 619-24.

  41. 41.

    Nassberger L, DePierre JW: Uptake, distribution and elimination of 3H-gentamicin in different organs of the rat as determined by scintillation counting. Acta Pharmacol Toxicol (Copenh). 1986, 59: 356-61.

  42. 42.

    Penson RT, Joel SP, Roberts M, Gloyne A, Beckwith S, Slevin ML: The bioavailability and pharmacokinetics of subcutaneous, nebulized and oral morphine-6-glucuronide. Br J Clin Pharmacol. 2002, 53: 347-54. 10.1046/j.1365-2125.2002.01554.x.

  43. 43.

    Ludden TM, Beal SL, Sheiner LB: Comparison of the Akaike Information Criterion, the Schwarz criterion and the F test as guides to model selection. J Pharmacokinet Biopharm. 1994, 22: 431-45.

  44. 44.

    Schwartz IL, Schacthter D, Freinkel N: The measurement of extracellular fluid in man by means of a constant infusion technique. Journal of Clinical Investigaton. 1945, 28: 1117-1125.

  45. 45.

    Deane N, Schreiner GE, Robertson JS: The velocity of distribution of sucrose between the plasma and interstitial fluid, with reference to the use of sucrose for the extracellular fluid in man. J Clin Investigation. 1951, 30: 1463-1468.

  46. 46.

    Arancibia A, Drouguett MT, Fuentes G, Gonzalez G, Gonzalez C, Thambo S, Palombo G: Pharmacokinetics of amoxicillin in subjects with normal and impaired renal function. Int J Clin Pharmacol Ther Toxicol. 1982, 20: 447-53.

  47. 47.

    Hanna MH, D'Costa F, Peat SJ, Fung C, Venkat N, Zilkha TR, Davies S: Morphine-6-glucuronide disposition in renal impairment. Br J Anaesth. 1993, 70: 511-4.

  48. 48.

    Ladegaard-Pedersen HJ, Engell HC: A comparison of the distribution volumes of inulin and (51 Cr)EDTA in man and nephrectomized dogs. Scand J Clin Lab Invest. 1972, 30: 267-70.

  49. 49.

    Bell DR, Watson PD, Renkin EM: Exclusion of plasma proteins in interstitium of tissues from the dog hind paw. Am J Physiol. 1980, 239: H532-H538.

  50. 50.

    Reed RK, Lepsoe S, Wiig H: Interstitial exclusion of albumin in rat dermis and subcutis in over- and dehydration. Am J Physiol. 1989, 257: H1819-27.

  51. 51.

    Nichols G, Nichols N, Weil WB, Wallace WM: The direct measurement of the extracellular phase of tissues. J Clinical Investigation. 1953, 32: 1299-1308.

  52. 52.

    Kruhoffer P: Inullin as an indicator for the extracellular space. Acta Physiol Scand. 1946, 11: 16-36.

  53. 53.

    Bass L: Flow dependence of first-order uptake of substances by heterogeneous perfused organs. J Theor Biol. 1980, 86: 365-76.

  54. 54.

    Haraldsson B, Rippe B: Restricted diffusion of CrEDTA and cyanocobalamine across the exchange vessels in rat hindquarters. Acta Physiol Scand. 1986, 127: 359-72.

  55. 55.

    Watson PD: Permeability of cat skeletal muscle capillaries to small solutes. Am J Physiol. 1995, 268: H184-93.

  56. 56.

    Cousins C, Mohammadtaghi S, Mubashar M, Strong R, Gunasekera RD, Myers MJ, Peters AM: Clearance kinetics of solutes used to measure glomerular filtration rate. Nucl Med Commun. 1999, 20: 1047-54.

  57. 57.

    Laker MF, Bull HJ, Menzies IS: Evaluation of mannitol for use as a probe marker of gastrointestinal permeability in man. Eur J Clin Invest. 1982, 12: 485-91.

  58. 58.

    Elia M, Behrens R, Northrop C, Wraight P, Neale G: Evaluation of mannitol, lactulose and 51Cr-labelled ethylenediaminetetra-acetate as markers of intestinal permeability in man. Clin Sci (Lond). 1987, 73: 197-204.

  59. 59.

    Lotsch J, Weiss M, Kobal G, Geisslinger G: Pharmacokinetics of morphine-6-glucuronide and its formation from morphine after intravenous administration. Clin Pharmacol Ther. 1998, 63: 629-39.

  60. 60.

    Penson RT, Joel SP, Clark S, Gloyne A, Slevin ML: Limited phase I study of morphine-3-glucuronide. J Pharm Sci. 2001, 90: 1810-6.

  61. 61.

    Sjovall J, Westerlund D, Alvan G: Renal excretion of intravenously infused amoxycillin and ampicillin. Br J Clin Pharmacol. 1985, 19: 191-201.

  62. 62.

    Arancibia A, Guttmann J, Gonzalez G, Gonzalez C: Absorption and disposition kinetics of amoxicillin in normal human subjects. Antimicrob Agents Chemother. 1980, 17: 199-202.

  63. 63.

    Lofgren S, Bucht G, Hermansson B, Holm SE, Winblad B, Norrby SR: Single-dose pharmacokinetics of dicloxacillin in healthy subjects of young and old age. Scand J Infect Dis. 1986, 18: 365-9.

  64. 64.

    Wiig H, DeCarlo M, Sibley L, Renkin EM: Interstitial exclusion of albumin in rat tissues measured by a continuous infusion method. Am J Physiol. 1992, 263: H1222-33.

  65. 65.

    Wiig H, Reed RK, Tenstad O: Interstitial fluid pressure, composition of interstitium, and interstitial exclusion of albumin in hypothyroid rats. Am J Physiol Heart Circ Physiol. 2000, 278: H1627-39.

  66. 66.

    Wiig H, Lund T: Relationship between interstitial fluid volume and pressure (compliance) in hypothyroid rats. Am J Physiol Heart Circ Physiol. 2001, 281: H1085-92.

  67. 67.

    Rippe C, Rippe B, Erlanson-Albertsson C: Capillary diffusion capacity and tissue distribution of pancreatic procolipase in rat. Am J Physiol. 1998, 275: G1179-84.

  68. 68.

    Larsson M, Johnson L, Nylander G, Ohman U: Plasma water and 51Cr EDTA equilibration volumes of different tissues in the rat. Acta Physiol Scand. 1980, 110: 53-7.

  69. 69.

    Katz J, Bonorris G, Golden S, Sellers AL: Extravascular albumin mass and exchange in rat tissues. Clin Sci. 1970, 39: 705-24.

  70. 70.

    Barratt TM, Walser M: Extracellular volume in skeletal muscle of the rat and dog: a comparison of radiosulphate and radiobromide spaces. Clin Sci. 1968, 35: 525-36.

  71. 71.

    Holmang A, Bjorntorp P, Rippe B: Tissue uptake of insulin and inulin in red and white skeletal muscle in vivo. Am J Physiol. 1992, 263: H1170-6.

  72. 72.

    Mullins RJ, Powers MR, Bell DR: Albumin and IgG in skin and skeletal muscle after plasmapheresis with saline loading. Am J Physiol. 1987, 252: H71-9.

  73. 73.

    Poole-Wilson PA, Cameron IR: ECS, intracellular pH, and electrolytes of cardiac and skeletal muscle. Am J Physiol. 1975, 229: 1299-1304.

  74. 74.

    Wittmers LE, Bartlett M, Johnson JA: Estimation of the capillary permeability coefficients of inulin in various tissues of the rabbit. Microvasc Res. 1976, 11: 67-78.

  75. 75.

    Wiig H, Reed RK: Interstitial compliance and transcapillary Starling pressures in cat skin and skeletal muscle. Am J Physiol. 1985, 248: H666-73.

  76. 76.

    Wiig H, Reed RK: Volume-pressure relationship (compliance) of interstitium in dog skin and muscle. Am J Physiol. 1987, 253: H291-8.

  77. 77.

    Sjogaard G, Saltin B: Extra- and intracellular water spaces in muscles of man at rest and with dynamic exercise. Am J Physiol. 1982, 243: R271-80.

  78. 78.

    Reinoso RF, Telfer BA, Rowland M: Tissue water content in rats measured by desiccation. J Pharmacol Toxicol Methods. 1997, 38: 87-92. 10.1016/S1056-8719(97)00053-1.

  79. 79.

    Mullins RJ, Bell DR: Changes in interstitial volume and masses of albumin and IgG in rabbit skin and skeletal muscle after saline volume loading. Circ Res. 1982, 51: 305-13.

Pre-publication history

  1. The pre-publication history for this paper can be accessed here:

Download references


I wish to thank Dr Mike Laker for providing additional information about an experimental subject in his published paper.

Author information

Correspondence to David G Levitt.

Additional information

Competing Interests

None declared.

Electronic supplementary material

Authors’ original submitted files for images

Rights and permissions

Reprints and Permissions

About this article


  • Inulin
  • PBPK Model
  • Total Body Water
  • Capillary Permeability
  • Extracellular Solute