Email updates

Keep up to date with the latest news and content from BMC Pharmacology and Toxicology and BioMed Central.

Open Access Highly Accessed Research article

Quantitative relationship between the octanol/water partition coefficient and the diffusion limitation of the exchange between adipose and blood

David G Levitt

Author Affiliations

Department of Integrative Biology and Physiology, University of Minnesota, 6-125 Jackson Hall, 321 Church St. S. E., Minneapolis, MN 55455, USA

BMC Clinical Pharmacology 2010, 10:1  doi:10.1186/1472-6904-10-1

The electronic version of this article is the complete one and can be found online at: http://www.biomedcentral.com/1472-6904/10/1


Received:15 September 2009
Accepted:7 January 2010
Published:7 January 2010

© 2010 Levitt; licensee BioMed Central Ltd.

This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Abstract

Background

The goal of physiologically based pharmacokinetics (PBPK) is to predict drug kinetics from an understanding of the organ/blood exchange. The standard approach is to assume that the organ is "flow limited" which means that the venous blood leaving the organ equilibrates with the well-stirred tissue compartment. Although this assumption is valid for most solutes, it has been shown to be incorrect for several very highly fat soluble compounds which appear to be "diffusion limited". This paper describes the physical basis of this adipose diffusion limitation and its quantitative dependence on the blood/water (Kbld-wat) and octanol/water (Kow) partition coefficient.

Methods

Experimental measurements of the time dependent rat blood and adipose concentration following either intravenous or oral input were used to estimate the "apparent" adipose perfusion rate (FA) assuming that the tissue is flow limited. It is shown that the ratio of FA to the anatomic perfusion rate (F) provides a measure of the diffusion limitation. A quantitative relationship between this diffusion limitation and Kbld-wat and Kow is derived. This analysis was applied to previously published data, including the Oberg et. al. measurements of the rat plasma and adipose tissue concentration following an oral dose of a mixture of 13 different polychlorinated biphenyls.

Results

Solutes become diffusion limited at values of log Kow greater than about 5.6, with the adipose-blood exchange rate reduced by a factor of about 30 for a solute with a log Kow of 7.36. Quantitatively, a plot of FA/F versus Kow is well described assuming an adipose permeability-surface area product (PS) of 750/min. This PS corresponds to a 0.14 micron aqueous layer separating the well-stirred blood from the adipose lipid. This is approximately equal to the thickness of the rat adipose capillary endothelium.

Conclusions

These results can be used to quantitate the adipose-blood diffusion limitation as a function of Kow. This is especially important for the highly fat soluble persistent organic chemicals (e.g. polychlorinated biphenyls, dioxins) whose pharmacokinetics are primarily determined by the adipose-blood exchange kinetics.

Background

Physiologically Based Pharmacokinetics (PBPK) refers to the approach of modeling drug kinetics using a realistic physiological description of the animal [1,2]. A central feature of this approach is the quantitative description of the tissue-blood exchange. The most basic approach (and the standard one) is to assume that this exchange is "flow limited" - that is, the venous blood leaving the capillary has equilibrated with the well mixed tissue space. For the flow limited model this tissue-blood exchange depends on only two parameters: a) the tissue perfusion rate (kg/min/kg); and b) the blood/tissues partition coefficient. This model has the major advantage that if one has determined the perfusion rate using one solute, the tissue-blood kinetic exchange can be predicted for any other solute for which the blood/tissue partition is known.

There are a number of solutes for which this flow limited model is clearly not valid and for which a "diffusion limited" model must be used. One such class are the large (e.g. inulin) or highly protein bound (e.g. dicloxacillin and ceftriaxone) extracellular solutes which have a significant capillary permeability limitation [3]. Similarly, a number of water soluble molecules such as actinomycin-D and methotrexate have been shown to have cell membrane limited uptake [4]. It is not surprising that these highly water soluble solutes with their low cell membrane permeability might be diffusion limited. These solutes represent relatively rare exceptions and the great majority of solutes that are used in PBPK modeling are highly lipid soluble with high (nearly infinite) cell membrane permeability and the predictions of the flow limited PBPK model generally provides good agreement with the experimental pharmacokinetics.

However, there are some notable exceptions to this general rule where the adipose-blood exchange of several highly lipid soluble molecules seems to be diffusion limited [5-7]. The purpose of this paper is to provide a detailed mechanistic analysis of the origin of the diffusion limitation of these highly lipid soluble molecules. It will be shown that the magnitude of the diffusion limitation can be directly related to and predicted by the lipid/water partition coefficient. This analysis provides general criteria for predicting the degree of adipose tissue diffusion limitation just from knowledge of the lipid/water (or octanol/water, see below) partition coefficient. This is the first detailed discussion of the quantitative relationship between lipid/water partition and diffusion limitation that I am aware of.

This analysis is especially important for the persistent organic pollutants (e.g. dioxins and polychlorinated biphenyls) whose pharmacokinetics are dominated by the kinetics of adipose-blood exchange [8]. The use of PBPK model predictions is essential for this solute class because it is not possible to accurately measure their pharmacokinetics in humans. As the following analysis shows, the adipose-blood exchange rate for the most highly lipid soluble solutes can be as much as 30 times slower than is predicted assuming flow limited kinetics.

The approach described here is based on an analysis of experimental measurements of the "apparent" adipose perfusion rate (FA). The "apparent" rate is the perfusion rate that would be predicted assuming that the exchange is blood flow limited. It will be shown that the ratio of this "apparent" rate and the true anatomic adipose perfusion rate (FA/F) can be used to quantitate the degree of diffusion limitation.

Methods

Dependence of the adipose-blood exchange rate (Clr) on the intrinsic capillary permeability-surface area product (PS) and the blood-water partition coefficient (Kbld-wat)

Figure 1 shows a schematic diagram of the factors involved in the solute exchange between the tissue and blood. The upper case letters in Figure 1 indicate the absolute blood (CB, CA, CV) or tissue (CT) concentration while the lower case indicates the free aqueous blood (cB) or tissue (cT) concentration. It is this free aqueous concentration that determines the diffusional exchange rate between the blood and tissue. The capillary concentration varies as a function of the distance (x) from the arterial end of the capillary. It is assumed that the tissue compartment is well mixed and can be represented by an average value that does not depend on x.

thumbnailFigure 1. Schematic diagram of the concentrations in capillary blood and adipose tissue and the factors involved in the tissue-blood solute exchange. The upper case letters indicate the absolute blood (CB, CA, CV) or tissue (CT) concentration while the lower case (cB, cT) indicates the free aqueous blood (cB) or tissue (cT) concentration.

It is assumed that the tissue consists of N capillaries/cm3, all with exactly the same geometry, blood flow, permeability, etc. It is also assumed that the relation between the capillary and tissue concentration is in a pseudo steady state. The steady state differential equation for the concentration in the capillary as a function of position is

(1)

where F is the tissue perfusion rate (kg/min/kg), P is the intrinsic capillary permeability defined in terms of the free water concentration (cm/min), N is the number of capillaries per cm3 tissue, and a is the capillary radius (cm). The total blood concentration CB can be related to the free water concentration (cB) using the blood/water partition coefficient (CB = Kbld-wat cB). Integrating eq [1] over the length of the capillary (L) and solving for the venous concentration leaving the capillary:

(2)

where PS (min-1) is the permeability-surface area product per tissue weight (S = 2πaNL). Equation [2] can be related to the fractional clearance or equilibration that occurs in one pass through the capillary:

(3)

For κ>>1, the venous concentration leaving the tissue (cv) is nearly equal to the tissue concentration (cT), the clearance (Clr) approaches 1 and the solute is flow limited. The clearance depends on both the intrinsic permeability (PS) and the blood-water partition coefficient. Solutes that have a high intrinsic permeability may be permeability (i.e. diffusion) limited if a large fraction in the blood is solute bound (i.e. large Kbld-wat).

Relationship between the anatomic adipose perfusion rate (F) and the "apparent" perfusion rate (FA)

The general, diffusion limited, differential equation describing the change in the well mixed tissue concentration (CT(t)) produced by a time dependent arterial input (CA(t)) is:

(4)

where CV(t) is the time dependent venous concentration leaving the tissue and cA and cV are the unbound free blood water concentration. By definition, if the tissue exchange is flow limited, then the venous concentration equilibrates with the tissue concentration (cv(t) = cT(t)):

(5)

Using the definition of clearance (eq. [6]), the general diffusion limited equation (eq. [5]) can be written in the form:

(6)

where FA is defined as the "apparent" perfusion rate. Equation [6] is identical to the flow limited expression (eq. [5]), the only difference is that FA is used in place of the true anatomic perfusion rate (F). Using eq [6], the general differential equation that describes the time dependent adipose tissue concentration (CT(t)) as a function of the arterial concentration (CA(t)) is then:

(7)

The adipose tissue concentration (CT) depends on just two parameters - the apparent perfusion rate (FA) and the adipose/blood partition coefficient (Kad-bld).

Relationship between Kbld-wat and octanol/water partition coefficient (Kow)

For the solutes used here the experimental value of Kbld-wat is not available and the following procedure was used to estimate it. It is well recognized that the tissue/blood partition for the highly lipid soluble molecules is roughly equal to the tissue/blood lipid concentration ratio [9-11], indicating that, as a first approximation, solutes in the blood are bound as if they were binding to an equivalent blood lipid (i.e. oil) fraction. Thus, the blood/water concentration can be approximated by:

(8)

where Koil-wat is the oil (e.g. olive oil)/water partition coefficient and fL is the equivalent blood lipid concentration. The last equality in eq. [8] results from the observation that, for non-polar solutes, Koil-wat is approximately equal to the octanol/water partition coefficient (Kow). A detailed analysis of the dependence of the relationship between Koil-wat and Kow on chemical structure is described in section 2 of the supplemental file (Additional File 1). Equation [8] is only an approximate estimate of Kbld-wat since there is clearly some specific binding to the cellular and protein blood components [12,13]. The analysis described below relies on using compounds with a very large range of Kow (103) so that small errors in the absolute value of Kbld-wat are relatively insignificant.

Additional file 1. I. Mathematical representation of the time dependent plasma concentration.

II. Prediction of oil/water partition coefficient using octanol/water coefficient.

Format: DOC Size: 230KB Download file

This file can be viewed with: Microsoft Word ViewerOpen Data

The value of fL (= Kbld-wat/Koil-wat = Kbld-oil ≈ 1/Kad-bld) can be estimated from the equilibrium partition between blood and adipose tissue (i.e.,oil). In humans, this ratio is about 0.005 for a large series of polychlorinated biphenyls [14]. In addition, in a recent large compilation of Ktis-bld by deBruyn and Gobas [15] the value of 1/Kad-bld is in the range of 0.005 for most of the very highly fat soluble solutes. This value of 0.005 is close to the directly measured values of blood lipid fraction in humans [9] and rats [10] and a value of fL = 0.005 will be assumed here. With the exception of PCB 2,2',5,5' (Table 1), all the Oberg et. al. [16] solutes had values of fL (= 1/Kad-bld) of about 0.005.

Table 1. Summary of the experimental data and the model analysis.

Experimental measurement of "apparent" perfusion rate (FA)

FA was determined from experimental measurements in the rat of the blood and adipose tissue concentration as a function of time. Given the time dependence of the blood concentration entering the tissue and assuming that the tissue is "apparently" flow limited, the adipose tissue concentration depends on only two parameters: FA and Kad-bld (eq. [7]). A continuous smooth plasma concentration curve was generated from the experimental plasma data points (see additional file 1: section I for details) and used as the arterial input (CA(t)). Using this time dependent arterial input, the standard flow limited organ equations (eq. [7]) were solved for the adipose concentration as a function of time. The two adipose parameters FA and Kad-bld were then adjusted to give the optimal fit to the experimental adipose concentration measurements. In most cases the experiments were carried out for long enough times (up to 132 days) that the adipose tissue had come close to equilibrating with the plasma and Kad-bld could be determined directly from this equilibrium value.

The numerical calculations were obtained using PKQuest_Java [17], a freely distributed PBPK software routine that can be downloaded from http://www.pkquest.com webcite. In addition, the free file "Persistent organics" contains the complete experimental data sets for most of the solutes discussed here and a detailed tutorial for reproducing these calculations.

Value of the anatomic adipose perfusion rate (F)

The determination of PS from eq. [9] requires an estimate of the anatomic adipose perfusion rate (F). The reported values for F in rat adipose tissue vary over a large range depending on methodologies, age, strain and condition (i.e. conscious or anesthetized). In addition, there are large differences in blood flow at different anatomic locations, varying from about 0.15 for epididymal to 0.55 kg/min/kg for mesenteric fat [18,19]. Most recent studies using labeled microspheres in conscious rats report values in the range of 0.18 to 0.25 kg/min/kg [20-22]. A value for F of 0.2 kg/min/kg was assumed in the following calculations.

Experimental data

The experimental data was obtained from the literature. The central results are those of Oberg et. al. [16] who simultaneously measured the rat plasma and adipose tissue concentration following an oral dose of a mixture of 13 different polychlorinated biphenyls, varying from 3 to 7 substituted chlorines with log Kow varying from 5.67 to 7.36 (table 1). These simultaneous chemical measurements represent the ideal data set for this analysis and this is the only measurement of this type that I am aware of.

In addition to the Oberg et. al. data set, the literature was screened and a number of other solutes were also modeled. One important experimental limitation for many solutes is that only the total C14 or H3 labeled equivalent was measured and the parent and metabolite compounds were not distinguished. Some of these tracer measurements were used in this analysis (see table 1) if there was supporting information that the labeled metabolite concentration in the plasma and tissue is relatively low. All of the experimental data that was used, along with comments about the experimental limitations, are summarized in table 1.

For published data that was available only in graphical form, the values were read off the graph using the program UN-SCAN-IT (version 6.0, Silk Scientific Corporation).

Summary of procedure used to estimate the adipose-blood diffusion limitation (Clr) and the intrinsic adipose capillary permeability-surface area product (PS)

Experimental measurements of the rat arterial blood (CA(t)) and adipose tissue (CT(t)) concentrations as a function of time were obtained from the literature. From a numerical solution of the differential eq. [7], the values of the two parameters FA and Kad-bld that gave the best agreement between the predicted and experimental adipose tissue concentration is then obtained. From the ratio of FA/F, the value of the adipose clearance (Clr) and the parameter κ for a given solute can be obtained from eqs. [3] and [6]. Using the relationship between Kbld-wat and Kow (eq. [8], κ can be expressed in terms of Kow:

(9)

As a first approximation, the value of the intrinsic PS (and β) for a given organ and species can be regarded as a constant, independent of the solute. This is because PS is proportional to the aqueous diffusion coefficient which is roughly proportional to 1/radius (Stokes-Einstein relation). For all the solutes considered here the molecular radius varies be less than a factor of 2, a variation that is negligible compared to the variations in Kow (≈ 103). Thus, assuming that β is constant, measurements of κ for a series of solutes with a wide range of values of Kow can be used to estimate the value of PS using the known values of fL (= 0.005) and F (0.2 kg/min/kg).

Results

The three panels in Figures 2 and 3 show the model fits to the experimental data of Oberg et. al. [16] for the solutes with the lowest (2,4,4' PCB, log Kow = 5.67), intermediate (2,3,3',4,4' PCB, log Kow = 6.65) and highest (2,2',3,4,4',5,5', log Kow = 7.36) value of Kow. Figure 2 shows the absolute plot and the insets shows the early time data. Figure 3 shows the semi-log plot of the same data. The red line is the smoothed fit to the plasma data which is used as the arterial input concentration to the organ. For all the solutes investigated, this smoothed curve provided a nearly perfect fit to the experimental blood concentration data points.

thumbnailFigure 2. Absolute plot of the time dependence of the experimental adipose (green triangle) and blood (yellow diamond) concentration and the model predictions for adipose (blue line) and blood (red line). The three panels show the results of Oberg et. al. [16] for a PCB with low (log Kow = 5.67), moderate (log Kow = 6.65) and high (log Kow = 7.36) octanol/water partition coefficient. The insets show the results at short times.

thumbnailFigure 3. Semi-log plot of the time dependence of the experimental adipose (green triangle) and blood (yellow diamond) concentration and the model predictions for adipose (blue line) and blood (red line). The three panels show the results of Oberg et. al. [16] for a PCB with low (log Kow = 5.67), moderate (log Kow = 6.65) and high (log Kow = 7.36) octanol/water partition coefficient. The insets show the results at short times.

The model fit to the adipose concentration data points is a function of two parameters: FA, and Kad-bld. These two parameters were adjusted to give the optimal fit to the adipose data, and these are the results that are shown in Figures 2 and 3 and are listed in Table 1. The adipose tissue fits are clearly not perfect. However, because of the large qualitative difference in the kinetics for the 3 solutes shown in Figure 2, there is no question that there are significant differences in the corresponding values of FA. This is dramatically illustrated in Figure 4 where the plots of the optimal values of FA are compared with the best fits that can be obtained using the value of FA for the neighboring solute in Figure 2. For example, in the middle panel of Figure 4, the red line shows the optimal fit (FA = 0.04) to the experimental adipose concentration data points for PCB 2,3,3',4,4' (solid circles). The green line is the predicted fit if the FA were equal to that of PCB 2,4,4' (FA = 0.15) and the blue line is the optimal fit if the FA was that of PCB 2,2',3,4,4',5,5' (FA = 0.005). The values of Kad-bld for the green and blue lines have been adjusted so that the curves fit the long time (i.e. equilibrium) data point.

thumbnailFigure 4. Comparison of the model fits to the adipose concentrations as a function of the "apparent" perfusion rate (FA) for the three solutes described in Figures 1 and 2. For each solute, the optimal value of FA for that solute (red line) is compared with the optimal fit that can be obtained using the value of FA for the other solutes (green and blue line). For the green and blue fits, the value of the adipose/blood partition (K) was adjusted to approximately fit the last time point. The value of FA and K for each line is listed at the bottom of each panel.

Table 1 lists the model values for FA and Kad-bld for the 13 PCBs studied by Oberg et. al. [16] along with a number of other solutes. One limitation of the Oberg et. al. data is that the first data point is at 6 hours and the early time adipose kinetics are missed. To address this, results are shown in Figure 5 for the very slowly metabolized PCB 2,2',4,4',5,5' which has been carefully investigated by Muhlebach et. al.[23] and has adipose data points at early times (5, 15, 30, 60 and 240 minutes). This data (including the early times) is well described by an FA of 0.02, which is in good agreement with the Oberg. et. al. data for PCBs with a similar Kow. Three of the solutes in Table 1 (thiopental, perchloroethylene and dichlorodiphenylsulfone) have low values of Kow (and Kbld-wat) and would be expected to be flow limited. Also listed in Table 1 are two pairs of solutes that are in the same chemical class and used similar experimental measurements. The first pair isTCDD and TBDD, for which both the pharmacokinetic and Kow measurements were made by the same laboratory. The second is the hexachlorobenzene and hexabromobenzene pair whose kinetics were described in the same publication. For both of these pairs, there is a significant decrease in the values for FA with an increase in Kow. Figure 6 compares the model FA dependence for TCDD and TBDD, similar to the plots in Figure 4.

thumbnailFigure 5. Absolute plot of the time dependence of the experimental adipose (green triangle) and blood (yellow diamond) concentration and the model predictions for adipose (blue line) and blood (red line) for the PCB 2,2',4,4',5,5' data of Muhlebach et. al. [23]. The insets show the results at short times.

thumbnailFigure 6. Perfusion rate dependence (FA) of the model fits to the experimental adipose concentration measurements (solid circles) for tetrachlorodibenzo-p-dioxin (TCDD, Wang et. al. [5]) and tetrabromodibenzo-p-dioxin (TBDD, Kedderis et. al. [35]). The red line is the predicted fit using the optimal value of FA for that solute and the green line is the fit using the FA of the other solute. For the green lines, the value of the adipose/blood partition (K) was adjusted to approximately fit the last time point.

The crucial test of this model of diffusion limitation is to determine if the fractional clearance (= Clr = FA/F) is described by eq. [9] Figure 7 shows the plot of the Oberg et. al. [16] values of FA/F (solid circles) as a function of Kow. The red line is the plot of 1-exp(-β/Kow), with β (eq. [9]) adjusted to provide the best fit to the Oberg et. al. data. For the value of β (= 0.75 × 106) that provides the best fit to the data in Figure 6, PS = βfLF = 750 min-1 (assuming fL = 0.005 and F = 0.2).

thumbnailFigure 7. Plot of the degree of diffusion limitation (FA/F) as a function of the octanol/water partition coefficient for the 13 PCBs studied by Oberg et. al. [16] (see Table 1 for numerical values). The red line is a plot of eq. [9] with a β (eq. [9]) of 0.75 × 106, adjusted to provide the best fit to the data.

Discussion

As described by eq. [3], the blood-water partition coefficient (Kbld-wat) has a critical role in determining whether tissue-solute exchange is diffusion limited. The capillary permeability limitation produced by this blood lipid solute binding is directly analogous to the carbon monoxide (CO) pulmonary diffusion limitation that results from the high affinity of hemoglobin for CO [24]. This same effect was also invoked to quantitate the capillary diffusion limitation of albumin bound β-lactam antibiotics [3]. As seen in eq. [3], given a finite intrinsic permeability (PS), the solute must become diffusion limited if Kbld-wat is large enough. The only question is at what value of Kbld-wat does the solute become diffusion limited?

To answer this question, the degree of diffusion limitation was determined for a series of solutes with varying values of octanol/water partition coefficient (Kow). The "diffusion limitation" was quantitated in rat adipose tissue from experimental measurements of the ratio of the "apparent" perfusion rate (FA) relative to the anatomic perfusion rate (F). Assuming that Kbld-wat is roughly proportional to Kow (eq. [8]), one should observe an increase in diffusion limitation as Kow increases. One would expect that this relationship would only be roughly satisfied because the exact relationship between Kow and Kbld-wat should depend on the specific blood protein binding and the detailed chemical structure [12,13]. This is the reason that it is especially important to have the Oberg et. al. [16] experimental data set for a large series of solutes with similar chemical structures and with Kow varying over a 50 fold range. Another factor that introduces variations in the relation between diffusion limitation and Kow is that there are large variations in the experimental measurements of Kow, depending on the experimental technique that is used. An advantage of the Oberg et. al. data set is that the values of Kow for all the solutes was determined using the same approach [25]. As shown in Table 1, for the Oberg et. al. data there is a qualitative increase in the degree of diffusion limitation (i.e. reduction in "apparent" perfusion rate FA) as Kow increases. Quantitatively, this diffusion limitation is quite will described by eq. [9] (Figure 7).

Also listed in Table 1 are the values of FA for some other solutes. Three of these solutes (thiopental, perchloroethylene, and dichlorodiphenylsulfone) have very small values of Kow and one would expect them to be flow limited. Consistent with this prediction, their values of FA are close to the assumed anatomic value (F) of 0.2 kg/min/kg.

Results are also shown in table 1 for hexachlorobenzene, hexabromobenzene, TCDD, TBDD and PCDE. The hexachlorobenzene and hexabromobenzene pair are directly comparable since their kinetics were described in the same publication using identical procedures. Similarly, the TCDD and TBDD pair results were carried out by the same lab using similar procedures. For both of these solute pairs with very similar structures, there is a significant increase in diffusion limitation (i.e. decrease in FA) with an increase in Kow (see also Figure 6 for a sensitivity analysis of the TCDD and TBDD results). These five solutes have values of Kow in the same range as the flow limited solutes studied by Oberg et. al. [16] and, qualitatively, have a similar increase in diffusion limitation (decrease in FA) with increasing Kow. Quantitatively, these solutes do not fall on exactly the same FA/F versus Kow curve (Figure 7) that was obtained for the Oberg solutes. As discussed above, this is not surprising because one would predict that the proportionality between Kow and Kbld-wat should have some solute dependence and should also depend on the specific details of, for example, the Kow measurements.

From the plot of FA/F versus Kow for the Oberg et. al. data (Figure 7), one can estimate that the corresponding value of the "intrinsic" PS for rat adipose tissue is about 750/min (see Results). To put this very large value in perspective, the highest capillary permeability that has been directly measured is a PS of about 1/min for Na+ in heart capillaries [26]. Solutes with a PS of 750/min only become diffusion limited because they have a very large value of Kbld-wat (≈ 5000).

One can use this value of PS to estimate the "equivalent" thickness of the rat blood-adipose diffusion limiting barrier. The permeability is equal to:

(10)

where DW is the aqueous diffusion coefficient and W is the thickness. Thus, relating W to PS:

(11)

where is S is the rat adipose capillary surface area. Using values of 5 × 10-6 cm2/sec [27] for PCB DW and 35 cm2/cm3 [28] for rat S, the value of W is about 1.4 × 10-5 cm for a PS of 750/min. This value for W (0.14 μm) is approximately equal to the thickness of the rat adipose capillary endothelium (0.25 μm) [29] and seems a reasonable estimate of the aqueous diffusive barrier between the blood (where the solutes are well-stirred and bound to albumin and the lipoproteins) and the adipose lipid.

An inherent assumption in this analysis is that the rate limiting step is the aqueous diffusion across the capillary wall and that the diffusion in the capillary blood and adipose tissue is not limiting. Diffusion in the blood should not be rate limiting because it is well stirred and the solute is bound at a high concentration relative to the water. The relative rates of diffusion in the aqueous and fat tissue are described by:

(12)

and the ratio of lipid to aqueous diffusion is:

(13)

Since the lipid (i.e. olive oil) viscosity is about 85 times greater than water, the Stokes-Einstein relation predicts that the lipid diffusion coefficient is about 85 times smaller than in water. However, this relation is valid only if the solute is much larger than the solvent, which is not true for diffusion in olive oil. For small solutes, such as O2 or N2, the diffusion coefficient in olive oil is only about 3 times smaller than in water [30]. Rogacheva et. al. [31] measured diffusion coefficients of 2-nonanone and benzaldehyde in oil that were about 10 times less than the water value. As a rough estimate, it is assumed that the diffusion coefficient in olive oil for the solutes studied here is about 20 times less than in water (DL/DW = 0.05). Since Koil-wat for the solutes considered here is 105 or greater, aqueous diffusion (eq. 13) is clearly the rate limiting step.

A novel feature of this analysis is the approach that was used to determine the "apparent" adipose perfusion. In the previous publications from which this experimental data was obtained, the blood and adipose concentrations were simultaneously modeled using the complete PBPK multi-tissue model. In many cases the model fit to the blood concentration is only poorly approximated by this PBPK model, presumably because of errors in the model assumptions (e.g. non-linearity). In contrast, in this analysis the exact experimental blood concentration is fit by a smooth curve so that the model input to the adipose tissue is identical to the experimental arterial concentration. The adipose perfusion rate is then adjusted to give the best fit to the experimental adipose concentration. This provides a significantly more accurate estimate of adipose perfusion then the total PBPK model fit with the incorrect blood concentration.

It is of interest to try to extrapolate these rat results to humans. The degree of diffusion limitation for a given Kow is determined by the parameter β = PS/(fLF) (eq. [9]). P corresponds to the diffusion limiting aqueous layer and should be similar for rat and human. Since rats [32] and humans [33] have similar adipose capillary densities, they should have similar capillary surface areas (S). Also, the blood lipid fraction (fL) is similar for rats and humans. One factor that does differ between rats and humans is the adipose blood flow (F), with the average human value of about 0.05 kg/kg/min [34] about 1/4 the rat value (0.2 kg/kg/min). Thus, the human value for β is about 4 times larger than the rat value and there should be corresponding less diffusion limitation for the same value of Kow.

Conclusions

Although it has been previously recognized that some highly fat soluble persistent organic chemicals are diffusion limited [5-7], this analysis provides the first physical explanation of this diffusion limitation along with its quantitative dependence on Kow. This diffusion limitation follows directly from the basic physiology of the blood tissue exchange. There must be some finite aqueous unstirred layer between the blood and the adipose lipid and this layer will become rate limiting if Kbld-wat is large enough. The results described here show that this diffusion limitation reduces the apparent rat adipose perfusion rate from the anatomic value of about 0.2 kg/min/kg for the flow limited solutes (log Kow < 5) to a value of about 0.005 kg/min/kg for a PCB with a log Kow of 7.36. The thickness of the limiting layer estimated from this analysis is about 0.14 μm, approximately equal to the thickness of the adipose capillary epithelial cell.

List of abbreviations

PBPK: Physiologically based pharmacokinetics; F: anatomic organ perfusion rate (kg/min/kg); FA: "apparent" perfusion rate assuming the organ kinetics are flow limited; a: capillary radius (cm); L: capillary length (cm); N: capillary density (#/cm3); S: 2πaNL: capillary surface area surface area per cm3 tissue (1/cm); P: intrinsic permeability (cm/min); PS: intrinsic organ permeability-surface area product (min-1); Clr: fraction of solute that equilibrates across capillary in one pass; Ci: total concentration (mole/kg) in organ i. i: B(blood), A(artery),V(vein), T(adipose); ci: free concentration in water (moles/cm3); fL: "equivalent" lipid fraction of blood; Kbld-wat: blood/water partition coefficient; Kad-wat: adipose/water partition coefficient; Kad-bld: adipose/blood partition coefficient; Koil-wat: oil/water partition coefficient; Kow: octanol/water partition coefficient; Dw and DL: Diffusion coefficient in water and lipid (cm2/sec); W: thickness (cm) of equivalent aqueous layer corresponding to intrinsic permeability; PCB: polychlorinated biphenyl; PCDE: pentachlorodiphenyl ether; TCDD: tetrachlorodibenzo-p-dioxin; TBDD: tetrabromodibenzo-p-dioxin.

Competing interests

The author declares that they have no competing interests.

Authors' contributions

DGL is entirely responsible for the contents of this paper and has read and approved the final manuscript.

References

  1. Reddy MB, Yang RSH, Clewell HJ, Andersen ME: Physiologically Based Pharmacokinetic Modeling. Hoboken, New Jersey: Wiley-Interscience; 2005.

  2. Edginton AN, Theil FP, Schmitt W, Willmann S: Whole body physiologically-based pharmacokinetic models: their use in clinical drug development.

    Expert Opin Drug Metab Toxicol 2008, 4:1143-1152. PubMed Abstract | Publisher Full Text OpenURL

  3. Levitt DG: PKQuest: capillary permeability limitation and plasma protein binding - application to human inulin, dicloxacillin and ceftriaxone pharmacokinetics.

    BMC Clin Pharmacol 2002, 2:7. PubMed Abstract | BioMed Central Full Text | PubMed Central Full Text OpenURL

  4. Lutz RJ, Dedrick RL, Zaharko DS: Physiological pharmacokinetics: an in vivo approach to membrane transport.

    Pharmacol Ther 1980, 11:559-592. PubMed Abstract | Publisher Full Text OpenURL

  5. Wang X, Santostefano MJ, Evans MV, Richardson VM, Diliberto JJ, Birnbaum LS: Determination of parameters responsible for pharmacokinetic behavior of TCDD in female Sprague-Dawley rats.

    Toxicol Appl Pharmacol 1997, 147:151-168. PubMed Abstract | Publisher Full Text OpenURL

  6. Kedderis LB, Mills JJ, Andersen ME, Birnbaum LS: A physiologically based pharmacokinetic model for 2,3,7,8-tetrabromodibenzo-p-dioxin (TBDD) in the rat: tissue distribution and CYP1A induction.

    Toxicol Appl Pharmacol 1993, 121:87-98. PubMed Abstract | Publisher Full Text OpenURL

  7. Kramer HJ, Drenth H, vandenBerg M, Seinen W, DeJongh J: Physiologically based pharmacokinetic model for tetrachlorobenzyltoluenes in rat: comparison of in vitro and in vivo metabolic rates.

    Toxicol Sci 2001, 63:22-28. PubMed Abstract | Publisher Full Text OpenURL

  8. Persistent organic pollutants Berlin Heidelberg: Springer-Verlag; 2003.

  9. Molen GW, Kooijman SA, Slob W: A generic toxicokinetic model for persistent lipophilic compounds in humans: an application to TCDD.

    Fundam Appl Toxicol 1996, 31:83-94. PubMed Abstract | Publisher Full Text OpenURL

  10. Emond C, Charbonneau M, Krishnan K: Physiologically based modeling of the accumulation in plasma and tissue lipids of a mixture of PCB congeners in female Sprague-Dawley rats.

    J Toxicol Environ Health A 2005, 68:1393-1412. PubMed Abstract | Publisher Full Text OpenURL

  11. Brown JF Jr, Lawton RW: Polychlorinated biphenyl (PCB) partitioning between adipose tissue and serum.

    Bull Environ Contam Toxicol 1984, 33:277-280. PubMed Abstract | Publisher Full Text OpenURL

  12. Parham FM, Kohn MC, Matthews HB, DeRosa C, Portier CJ: Using structural information to create physiologically based pharmacokinetic models for all polychlorinated biphenyls.

    Toxicol Appl Pharmacol 1997, 144:340-347. PubMed Abstract | Publisher Full Text OpenURL

  13. Noren K, Weistrand C, Karpe F: Distribution of PCB congeners, DDE, hexachlorobenzene, and methylsulfonyl metabolites of PCB and DDE among various fractions of human blood plasma.

    Arch Environ Contam Toxicol 1999, 37:408-414. PubMed Abstract | Publisher Full Text OpenURL

  14. Wolff MS, Thornton J, Fischbein A, Lilis R, Selikoff IJ: Disposition of polychlorinated biphenyl congeners in occupationally exposed persons.

    Toxicol Appl Pharmacol 1982, 62:294-306. PubMed Abstract | Publisher Full Text OpenURL

  15. deBruyn AM, Gobas FA: The sorptive capacity of animal protein.

    Environ Toxicol Chem 2007, 26:1803-1808. PubMed Abstract | Publisher Full Text OpenURL

  16. Oberg M, Sjodin A, Casabona H, Nordgren I, Klasson-Wehler E, Hakansson H: Tissue distribution and half-lives of individual polychlorinated biphenyls and serum levels of 4-hydroxy-2,3,3',4',5-pentachlorobiphenyl in the rat.

    Toxicol Sci 2002, 70:171-182. PubMed Abstract | Publisher Full Text OpenURL

  17. Levitt DG: PKQuest_Java: free, interactive physiologically based pharmacokinetic software package and tutorial.

    BMC Res Notes 2009, 2:158. PubMed Abstract | BioMed Central Full Text | PubMed Central Full Text OpenURL

  18. West DB, Prinz WA, Greenwood MR: Regional changes in adipose tissue blood flow and metabolism in rats after a meal.

    Am J Physiol 1989, 257:R711-716. PubMed Abstract | Publisher Full Text OpenURL

  19. Crandall DL, Goldstein BM, Huggins F, Cervoni P: Adipocyte blood flow: influence of age, anatomic location, and dietary manipulation.

    Am J Physiol 1984, 247:R46-51. PubMed Abstract | Publisher Full Text OpenURL

  20. Wada DR, Harashima H, Ebling W, Osaki EW, Stanski DR: Effects of thiopental on regional blood flows in the rat.

    Anesthesiology 1996, 84:596-604. PubMed Abstract | Publisher Full Text OpenURL

  21. Kuwahira I, Gonzalez NC, Heisler N, Piiper J: Regional blood flow in conscious resting rats determined by microsphere distribution.

    J Appl Physiol 1993, 74:203-210. PubMed Abstract | Publisher Full Text OpenURL

  22. Delp MD, Evans MV, Duan C: Effects of aging on cardiac output, regional blood flow, and body composition in Fischer-344 rats.

    J Appl Physiol 1998, 85:1813-1822. PubMed Abstract | Publisher Full Text OpenURL

  23. Muhlebach S, Wyss PA, Bickel MH: Comparative adipose tissue kinetics of thiopental, DDE and 2,4,5,2',4',5'-hexachlorobiphenyl in the rat.

    Xenobiotica 1985, 15:485-491. PubMed Abstract | Publisher Full Text OpenURL

  24. Forster RE: Exchange of gases between alveolar air and pulmonary capillary blood: pulmonary diffusing capacity.

    Physiol Rev 1957, 37:391-452. PubMed Abstract | Publisher Full Text OpenURL

  25. Hawker DW, Connell DW: Octanol-water partition coefficients of polychlorinated biphenyl congeners.

    Environ Sci Technol 1988, 22:382-387. Publisher Full Text OpenURL

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

  27. Werner D, Ghosh U, Luthy RG: Modeling polychlorinated biphenyl mass transfer after amendment of contaminated sediment with activated carbon.

    Environ Sci Technol 2006, 40:4211-4218. PubMed Abstract | Publisher Full Text OpenURL

  28. Paaske WP: Absence of restricted diffusion in adipose tissue capillaries.

    Acta Physiol Scand 1977, 100:430-436. PubMed Abstract | Publisher Full Text OpenURL

  29. Fraselle-Jacobs A, Jeanjean M, Hauser N, Remacle C: Effect of aging on the morphology of epididymal adipose tissue in the rat.

    Exp Gerontol 1987, 22:389-408. PubMed Abstract | Publisher Full Text OpenURL

  30. Davidson D, Eggleton P, Foggie P: The diffusion of atmospheric gases through fats and oils.

    Q J Exp Physiol Cogn Med Sci 1952, 37:91-105. PubMed Abstract | Publisher Full Text OpenURL

  31. Rogacheva S, Espinosa-Diaz MA, Voilley A: Transfer of aroma compounds in water-lipid systems: binding tendency of beta-lactoglobulin.

    J Agric Food Chem 1999, 47:259-263. PubMed Abstract | Publisher Full Text OpenURL

  32. Kampf C, Bodin B, Kallskog O, Carlsson C, Jansson L: Marked increase in white adipose tissue blood perfusion in the type 2 diabetic GK rat.

    Diabetes 2005, 54:2620-2627. PubMed Abstract | Publisher Full Text OpenURL

  33. Pasarica M, Sereda OR, Redman LM, Albarado DC, Hymel DT, Roan LE, Rood JC, Burk DH, Smith SR: Reduced adipose tissue oxygenation in human obesity: evidence for rarefaction, macrophage chemotaxis, and inflammation without an angiogenic response.

    Diabetes 2009, 58:718-725. PubMed Abstract | Publisher Full Text OpenURL

  34. Levitt DG: Heterogeneity of human adipose blood flow.

    BMC Clin Pharmacol 2007, 7:1. PubMed Abstract | BioMed Central Full Text | PubMed Central Full Text OpenURL

  35. Kedderis LB, Diliberto JJ, Birnbaum LS: Disposition and excretion of intravenous 2,3,7,8-tetrabromodibenzo-p-dioxin (TBDD) in rats.

    Toxicol Appl Pharmacol 1991, 108:397-406. PubMed Abstract | Publisher Full Text OpenURL

  36. Ebling WF, Wada DR, Stanski DR: From piecewise to full physiologic pharmacokinetic modeling: applied to thiopental disposition in the rat.

    J Pharmacokinet Biopharm 1994, 22:259-292. PubMed Abstract | Publisher Full Text OpenURL

  37. Freese E, Levin BC, Pearce R, Sreevalsan T, Kaufman JJ, Koski WS, Semo NM: Correlation between the growth inhibitory effects, partition coefficients and teratogenic effects of lipophilic acids.

    Teratology 1979, 20:413-440. PubMed Abstract | Publisher Full Text OpenURL

  38. Dallas CE, Muralidhara S, Chen XM, Ramanathan R, Varkonyi P, Gallo JM, Bruckner JV: Use of a physiologically based model to predict systemic uptake and respiratory elimination of perchloroethylene.

    Toxicol Appl Pharmacol 1994, 128:60-68. PubMed Abstract | Publisher Full Text OpenURL

  39. Leo A, Hansch C, Elkins D: Partition coefficients and their uses.

    Chemical Reviews 1971, 71:525-616. Publisher Full Text OpenURL

  40. Parham FM, Matthews HB, Portier CJ: A physiologically based pharmacokinetic model of p,p'-dichlorodiphenylsulfone.

    Toxicol Appl Pharmacol 2002, 181:153-163. PubMed Abstract | Publisher Full Text OpenURL

  41. Yamaguchi Y, Kaswano M, Tatsukawa R: Tissue distribution and excretion of hexabromobenzene and hexachlorobenzene administered to rats.

    Chemosphere 1986, 15:453-459. Publisher Full Text OpenURL

  42. Bahudur NP, Shiu W, Boocock DGB, Mackay D: Temperature dependence of octanol-water partition coefficient for selected chlorobenzenes.

    J Chem Eng Data 1997, 42:685-688. Publisher Full Text OpenURL

  43. Kuramochi H, Maeda K, Kawamoto K: Measurements of water solubilities and 1-octanol/water partition coefficients and estimations of Henry's lay constants for brominated benzenes.

    J Chem Eng Data 2004, 49:720-724. Publisher Full Text OpenURL

  44. Jackson JA, Diliberto JJ, Birnbaum LS: Estimation of octanol-water partition coefficients and correlation with dermal absorption for several polyhalogenated aromatic hydrocarbons.

    Fundam Appl Toxicol 1993, 21:334-344. PubMed Abstract | Publisher Full Text OpenURL

  45. Komsta E, Chu I, Villeneuve DC, Benoit FM, Murdoch D: Tissue distribution metabolism and excretion of 2,2',4,4',5-pentachlorodiphenyl ether in the rat.

    Arch Toxicol 1988, 62:258-262. PubMed Abstract | Publisher Full Text OpenURL

  46. Kurz J, Ballschmiter K: Vapour pressures, aqueous solubilities, Henry's law constants, partition coefficients between gas/water (Kgw), n-octanol/water (Kow) and gas/n-octanol (Kgo) of 106 polychlorinated diphenyl ethers (PCDE).

    Chemosphere 1999, 38:573-586. PubMed Abstract | Publisher Full Text OpenURL

Pre-publication history

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

http://www.biomedcentral.com/1472-6904/10/1/prepub