![]() |
|
|
| |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
Research Service, V.A. Medical Center, Milwaukee, Wisconsin (D.L.R., S.B.A.); Department of Anesthesiology (D.L.R.) and Department of Pulmonary and Critical Care Medicine (S.H.A.), Medical College of Wisconsin, Milwaukee, Wisconsin; and Department of Biomedical Engineering, Marquette University, Milwaukee, Wisconsin (S.H.A.)
(Received March 30, 2004; accepted June 9, 2004)
| Abstract |
|---|
|
|
|---|
P-Glycoprotein has been referred to as an atypical ATP-dependent transporter that is highly promiscuous based on the hundreds of compounds that have been identified as substrates (Sharom, 1997
). The physical chemical properties of P-gp substrates are diverse and most commonly include some degree of lipophilicity and a positive charge, and/or an amine nitrogen (Frezard et al., 2001
). Seelig (1998
) reviewed 100 compounds that interact with P-gp, and more than 60 were basic lipophilic amines, a property shared with compounds that readily accumulate in lung tissue. The lung is well known for its extensive accumulation of lipophilic amines, and this similarity to the substrate specificity of P-gp plus the very nature of this drug efflux transporter strongly suggest that P-gp-mediated drug efflux could play a moderating role in the pulmonary accumulation of amine drugs (Junod, 1976
; Wilson et al., 1979
; Roerig et al., 1983
; 1984
; Bend et al., 1985
). However, being a substrate for P-gp does not ensure that this efflux transporter will significantly moderate drug disposition in lung tissue. For an intact organ, tissue accumulation is dependent on multiple, and often competing, vascular and intracellular dispositional processes. These include rate of delivery to the organ (blood flow), diffusion into or out of specific cell types, interactions such as binding to intracellular sites, or accumulation in subcellular organelles. An additional factor is saturation of one or more of these kinetic processes, which could reduce the impact of the saturable process on overall tissue accumulation of the drug.
This work is aimed at developing a kinetic model of pulmonary drug disposition to identify the contribution of the different kinetic processes, including P-gp-mediated drug efflux, in the lung tissue accumulation of lipophilic amines. Because of its ease of ex vivo perfusion, the ability to manipulate perfusion variables, and expression of P-gp in normal lung tissue, the lung may also serve as an organ model in general. For this study, we selected the lipophilic amine dye, rhodamine 6G (R6G) as the test substrate based on its reported uptake in cells (Loetchutinat et al., 2003
).
| Materials and Methods |
|---|
|
|
|---|
Animals. All studies were carried out with New Zealand White rabbits of either gender with a mean ± S.E. weight of 2.65 ± 0.2 kg. The rabbits were housed one per cage, with free access to food and water, and maintained on a 12-h light/dark cycle. All rabbits were purchased from New Franken Research Rabbits (New Franken, WI). All animal care and treatment procedures were approved by the V.A. Medical Center Research Institutional Animal Care and Use Committee before initiating these studies.
Isolated Lung Perfusion System. The isolated perfused rabbit lung (IPL) system was similar to that previously described (Roerig et al., 1992
; Audi et al., 1995
, 1996
). Briefly, rabbits were anesthetized with sodium pentobarbital (25-40 mg/kg. i.v., ear vein), and the carotid artery was cannulated. The rabbits were heparinized (1000 units/kg) and exsanguinated via the carotid cannula. The chest was opened and the pulmonary artery, pulmonary vein, and trachea were cannulated. The lung was removed from the chest and suspended in the perfusion apparatus by the cannulas, and the perfusion system was primed with a perfusate consisting of a physiological salt solution containing 4.5% bovine serum albumin. Perfusate was initially pumped through the IPL at 200 ml/min (Master Flex roller pump; Cole-Parmer Instrument Co., Vernon Hills, IL), and the uniformity of blanching after the residual blood was washed out of the lung was used to estimate the extent of perfusion of the IPL. Only lungs greater than 90% perfused were used for these studies. The venous outflow was then directed to the reservoir supplying the pump. Additional perfusate was added to the reservoir to give a recirculating perfusate volume of 134 ± 1 ml. The pulmonary artery pressure, referenced to the level of the left atrium, was 8.1 ± 0.3 cm of H2O at end expiration, and IPLs that exhibited an arterial perfusion pressure greater than 10 cm of H2O were not studied. The perfused lung was maintained at 37°C. The lungs were continuously ventilated (10 breaths/min) with a gas mixture containing 15% O2, 5% CO2, balance N2 with end-inspiratory and -expiratory airway pressures of 7.7 ± 0.1 and 2.1 ± 0.1 cm of H2O, respectively. At the end of each experiment the lungs were weighed and lyophilized to a constant weight. The wet/dry weight ratio of all lungs was 6.11 ± 0.06.
Experimental Protocols. The uptake of R6G with time was determined by adding known amounts of R6G to the recirculating reservoir and then removing 1.7 ml of perfusate samples from the reservoir at fixed time intervals out to 120 min, depending on the particular experiment. The sampling interval after R6G addition to the reservoir was every 5 min for the first 30 min, every 10 min between 30 and 60 min, and every 15 min thereafter. Depending on the experiment, the P-gp inhibitor, GF120918, was added to the perfusate reservoir 5 min before the R6G or at various times after R6G addition. A standard concentration curve for R6G was constructed using perfusate that had recirculated through the lung before R6G addition and a portion of the stock R6G solution that was added to the recirculating reservoir. All samples and standards were centrifuged to remove trace red blood cells, and the absorbance at 540 nm was used to determine the R6G perfusate concentration. The amount of R6G taken up by the lung was determined by its disappearance from the recirculating reservoir corrected for the amount of R6G removed from the perfusion system during sampling. Each lung was used only once, and the number (n) of individual rabbit lungs used in the different experiments is given in the figure legends.
| Results |
|---|
|
|
|---|
|
Figure 2 depicts a family of curves showing the increase in the total nanomoles of R6G in the rat lungs versus time with initial perfusate R6G concentrations ranging from 1.25 µM to 12.6 µM. This experiment provides one of the discriminating data sets necessary to estimate values of the kinetic model parameters (see "Kinetic Model and Data Analysis"). This experiment also provides insight into possible saturation of dispositional processes involved in R6G lung accumulation. For example, in the absence of any saturation, the fraction of the initial amount of R6G in the perfusate in the lung at the end of the perfusion should be constant over the R6G perfusate concentration range studied. However, as shown in Fig. 3, the fraction of total R6G in the lung at 90 min increased as a function of total R6G initially present in the reservoir.
|
|
The effect of GF120918 concentration on R6G accumulation in the lung was also determined. Lungs were perfused for 5 min with perfusate concentrations of GF120918 of 0, 2, 10, 20, and 70 µM before recirculation of 2.5 µM R6G. Accumulation of R6G in the lung was then determined at each inhibitor concentration and is shown in Fig. 4. The accumulation of R6G in the lung increased with increasing GF120918 perfusate concentration. This experiment also provides a data set necessary to estimate values of the kinetic model parameters as explained under "Kinetic Model and Data Analysis," described below. It is apparent in Fig. 4 that inhibition of P-gp is approaching a maximum, and the 20 µM GF120918 concentration (Fig. 4, open triangles) used to inhibit P-gp-mediated R6G efflux out of the lung in Figs. 1 and 5 results in near maximal P-gp inhibition.
|
|
In an additional study, 20 µM GF120918 was added to the recirculating perfusate at various times after the start of R6G recirculation through the lung. Figure 5 shows the increase in the total R6G in the lung versus time with GF120918 added at the times indicated by the arrows. Addition of the P-gp inhibitor caused a sudden increase in R6G accumulation that was greater the earlier the GF120918 was added to the recirculating perfusate.
Kinetic Model and Data Analysis. To evaluate the effect of P-gp and other cellular processes on R6G accumulation in the lung (Figs. 2, 3, 4, 5), we developed a parsimonious kinetic model for the pulmonary disposition of R6G. The model (Fig. 6) consists of three compartments, a vascular compartment and two extravascular compartments. The vascular compartment (compartment 1) represents the lung vascular region, reservoir, and tubing. The two extravascular compartments (compartments 2 and 3) represent distribution of R6G within the lung tissue. The model allows for R6G diffusion between compartments 1 and 2 with the rate of diffusion represented by the permeability-surface area product, PS (ml/min); P-gp-mediated efflux of R6G from compartment 2 into compartment 1, described by a Michaelis-Menten-type process with an apparent maximum efflux rate Vmax21 (nmol/min) and a Michaelis constant Km21 (nmol); accumulation of R6G in compartment 3, also described by a Michaelis-Menten-type process with an apparent maximum uptake rate Vmax23 (nmol/min) and a Michaelis constant Km23 (µM); and a passive intercompartmental transfer from compartment 3 to compartment 2, described by the first order rate constant k32 (min-1). Compartment 3 could represent intracellular binding sites or sequestration in organelles. The extensive pulmonary uptake of R6G suggests that this is a high-capacity, high-affinity intracellular pool for R6G. Thus, in eqs. 1 to 3, below, Vmax23 and Km23 are replaced by the first order rate constant k23 = Vmax23/(Km23V2) (min-1). The inhibitor GF120918 is assumed not only to inhibit P-gp-mediated R6G efflux with an inhibition constant Ki1 (µM), but also to inhibit R6G accumulation in compartment 3 with an inhibition constant Ki2 (µM). The kinetic model makes no assumption of the anatomical location of the P-gp protein but does assume that P-gp transports R6G from the lung tissue into the vascular space.
|
The temporal variations in the concentrations of the various species in the three compartments are described by eqs. 1 to 3.
![]() | (1) |
![]() | (2) |
![]() | (3) |
where [C1](t) is the vascular (perfusate) concentration (µM) of R6G; C2(t) and C3(t) are the amounts (nmol) of R6G in compartments 2 and 3, respectively; [i] is the concentration (µM) of the inhibitor GF120918 in the extravascular region;
is the inhibitor's extravascular-to-vascular partition coefficient; k23 (min-1) is an approximation to the Vmax23/(Km23V2) ratio, and V1 and V2 (ml) are the volumes of compartments 1 and 2, respectively.
The identifiable model parameters are the following parameter groups: PS (ml/min); kout = PS/V2 (min-1), a measure of the rate constant for diffusion-mediated lung efflux of R6G; (Km21V2) (nmol), measure of the affinity of P-gp for R6G; Vmax21 (nmol/min), measure of the capacity of P-gp for R6G; k23 and k32 (min-1), first order rate constants for the transfer of R6G in and out of the subcellular pool (compartment 3), respectively; and Ki1/
and Ki2/
, respective measures of the inhibition constants for P-gp-mediated R6G efflux and R6G accumulation in compartment 3. The volume of compartment 1, V1, was set at 134 ml.
The above model parameters are not estimable from a single set of data or experimental condition. In other words, the information content of a single set of data or experimental condition does not provide sufficiently discriminating information about all dominant processes to separately estimate all the model parameters. Thus, the approach that we have taken to estimate values for the kinetic model parameters was to carry out simultaneous nonlinear regression fits of the solution of eqs. 1 to 3 to all the mean data in Figs. 2 and 4 obtained by varying the perfusion conditions as indicated above. The model solution was corrected for the change in reservoir volume due to sampling. The estimated values of the model parameters are given in Table 1. The solid lines imposed on the data in Figs. 2 and 4 are the model fit to the data. The solid lines superimposed on the data in Figs. 3 and 5 are model predictions obtained by solving eqs. 1 to 3 with the model parameters set at the values estimated from the data in Figs. 2 and 4.
|
To examine the contribution of P-gp-mediated efflux to the pulmonary uptake of R6G, we compared the estimated values of the kinetic parameters kpgp = Vmax21/(Km21V2), kout = PS/V2, kin = PS/V1, and k23 = Vmax23/(Km23V2), which represent rate constants (min-1) for the competing dispositional processes, namely P-gp-mediated efflux, non-P-gp-mediated efflux, diffusion into the lung, and intracellular transfer into the subcellular pool (see Table 1).
A Fortran-based computer program (Fortran PowerStation, Version 4.0) was developed for evaluating eqs. 1 to 3 (kinetic model) and for estimating the model parameters descriptive of the various processes (Audi et al., 1998
). As described above, for the appropriate initial conditions, the model equations were solved numerically (Gear's method) using International Mathematics and Statistics Library (IMSL) subroutine DIPAG (IMSL Math/Library). The optimization for fitting the model to the data and parameter estimation was carried out using IMSL subroutine DBCLSF, which finds the best nonlinear least-squares solution using a modified Levenberg-Marquardt algorithm with finite difference Jacobian. The 95% asymptotic parameter confidence intervals (Table 1) were calculated from the sensitivity matrix as previously described (Audi et al., 1998
).
| Discussion |
|---|
|
|
|---|
In normal lung tissue, P-gp is expressed in both airway epithelial cells and lung endothelial cells (Lechapt-Zalcman et al., 1997
; Demeule et al., 2001
; Campbell et al., 2003
). Campbell et al. (2003
) suggested a protective role of epithelial P-gp as a barrier to xenobiotic transfer from alveoli to pulmonary interstitium and blood. We have focused on the role of P-gp in pulmonary drug accumulation from the circulation, where a protective role of P-gp from blood-borne xenobiotics is of great potential importance, as well as an even larger potential role in limiting pulmonary accumulation of a large number of systemic therapeutic agents where the lung is the target organ. The accessibility to this transporter from the vascular space and the reported presence of P-gp in pulmonary endothelial cells suggests that endothelial P-gp moderates pulmonary R6G accumulation (Demeule et al., 2001
). However, pulmonary epithelial P-gp must also be considered. Epithelial P-gp would transport and retain R6G in the airspace rather than transport it out of the lung back into the circulation. Inhibition of epithelial P-gp would therefore decrease net lung R6G accumulation. This is the exact opposite of our results, suggesting that epithelial P-gp contributes little to the observed kinetics of pulmonary R6G accumulation from the vascular space.
The promiscuous nature of P-gp substrate specificity parallels the affinity of lipophilic amine drugs for lung tissue. Depending on the drug, accumulation of lipophilic amines in the intact lung results from both rapidly and slowly equilibrating interactions with lung tissue (Audi et al., 1995
, 1996
, 1998
; Roerig et al., 1999
). To study this differential amine interaction with macromolecules and/or partitioning among subcellular organelles, we developed kinetic models of pulmonary drug disposition to identify and characterize kinetic parameters descriptive of different drug-tissue interactions in the vascular and extravascular space (Audi et al., 1998
, 2002
; Roerig et al., 1999
). We used a similar approach to develop a kinetic model for pulmonary R6G accumulation in the intact lung that includes kinetic parameters descriptive of the role of P-gp-mediated efflux and other drug-dispositional mechanisms inferred by the data in Figs. 1, 2, 3, 4. Based on Fig. 3, P-gp-mediated efflux of R6G is saturable and is assumed to follow Michaelis-Menten kinetics as proposed by others (Weiss and Kang, 2002
). The extensive net accumulation of R6G over the 90 min of recirculation (Figs. 1, 2, 4, and 5) indicates that at steady state, distribution is far in favor of the lung, suggesting a large apparent lung tissue volume of distribution. Since R6G is known and marketed as a mitochondrion-selective dye, compartment 3 (Fig. 6) represents, at least in part, sequestration in a subcellular organelle (Haugland, 2002
). Other slowly equilibrating R6G-lung tissue interactions may also be involved but cannot be separately identified with the present data. Weiss and Kang (2002
) developed a similar compartmental model to describe the role of P-gp in idarubicin disposition in the isolated perfused heart and proposed a carrier-mediated transport process from compartment 2 to 3. Again, this cannot be separately identified from the data presented here. The model (Fig. 6) developed from these considerations was simultaneously fit to the data in Figs. 2 and 4 and, as can be seen, resulted in a reasonably good fit (solid line) to the experimental data points. To explain (fit) all the data in Figs. 2 and 4, it was necessary to include an additional dispositional kinetic parameter, Ki2, that describes inhibition of R6G accumulation in compartment 3 by the P-gp inhibitor, GF120918. The model was further tested by using the model parameter estimates obtained from fitting the data in Figs. 2 and 4 to predict the shape of the outflow curves in Fig. 5 where the P-gp inhibitor GF120918 was added at various times after the start of R6G recirculation. As can be seen from the solid lines in Fig. 5, the model faithfully predicts the increase in R6G accumulation that occurs after GF120918 addition at different times during the recirculation period.
The kinetic model parameters obtained from numerically solving eqs. 1 to 3 are shown in Table 1. The effect of P-gp-mediated R6G efflux relative to the other dispositional kinetic processes can be examined by comparing the rate constants (min-1) for P-gp-mediated efflux (kpgp), non-P-gp-mediated efflux (kout), diffusion into the lung (kin), and sequestration into compartment 3 (k23) using the model parameter values in Table 1. If P-gp were fully inhibited or not present, sequestration in compartment 3 (k23 = 0.22 min-1) would dominate over back-diffusion (kout = 0.011 min-1), resulting in near unidirectional influx and sequestration in the lung. This is supported by the fact that after 90 min of R6G recirculation in the presence of an inhibitor (Fig. 1), 87.5 ± 1.1% (S.E.) of the R6G was in the lung. The value of kpgp of 1.44 min-1 is 2 orders of magnitude greater than that of kout, suggesting that P-gp-mediated efflux is by far the dominant efflux mechanism. In addition, kpgp is 15 times greater than kin = 0.093 min-1 and 6 times the rate constant for sequestration (k23 = 0.22 min-1), resulting in P-gp-mediated efflux significantly moderating the amount of R6G available for disposition into compartment 3 and the lung as a whole.
The inhibition of R6G accumulation into compartment 3 in the presence of GF120918 is interesting since it poses a mechanism that would offset the effect of P-gp inhibition. GF120918 is a lipophilic amine, and inhibition of sequestration of one lipophilic amine by another in the lung has been reported (Anderson et al., 1974
; Wilson et al., 1979
; Roerig et al., 1989
). This mechanism may be of some importance since reversal or modulation therapy using inhibitors of P-gp has been proposed and studied as a method of reversing the development of multidrug resistance in cancer chemotherapy (Ford and Hait, 1990
).
Others have proposed models to estimate kinetic parameters to provide insight into the role of P-gp-mediated drug efflux in drug disposition. Loetchutinat et al. (2003
) estimated the rate constants, ka and k, for P-gp-mediated efflux and passive membrane permeability, respectively, of different rhodamine dyes in a human leukemia cell line. For R6G, the ratio of ka/k was about 29, indicating that P-gp also dominated the rate of R6G efflux out of the cell in their studies. Wielinga et al. (2000
) also determined a kinetic parameter for passive permeation and a measure of P-gp-mediated efflux using four anthracycline analogs in an MDR cell line. They concluded that passive permeation plays a substantial role in determining the drug resistance for these anthracyclines. Wielinga et al. (2000
) also identified a rate constant, ks, for passive transport between an inner compartment and a shell compartment of the cell. Based on their data, like our findings with R6G, the rate constant for transport into an intracellular sequestered compartment was greater than that of the passive permeation parameter. The results of both of the above studies demonstrated that the difference for passive membrane permeability between the different compounds studied did not parallel the differences in the derived rate constants for P-gp-mediated efflux, with high lipophilicity having more effect on passive efflux. Weiss and Kang (2002
) developed a whole organ model for evaluating the role of P-gp in the uptake of idarubicin in the isolated perfused rat heart. Their three-compartment model contained kinetic parameters descriptive of different intracellular dispositional processes. They concluded that the rate of idarubicin uptake in the heart represents a balance between two opposing processes: active P-gp-mediated efflux and uptake into mitochondria. These studies, together with the present work, all serve to emphasize the importance of competing intracellular processes when defining the impact of P-gp-mediated drug efflux on overall tissue disposition.
In conclusion, these studies demonstrate a functional P-gp-mediated drug efflux transporter in the lung that significantly moderates the disposition of R6G between the vascular and extravascular space in the lung. Most importantly, the experimental protocol presented, together with the proposed kinetic model and appropriate experimental data sets, can provide estimates of kinetic parameters that are descriptive of the different dispositional mechanisms involved in pulmonary drug accumulation. Reliable estimates of all of these kinetic parameters will be essential in assessing the role of P-gp-mediated drug efflux in the net pulmonary accumulation of a large variety of therapeutic agents.
| Footnotes |
|---|
Article, publication date, and citation information can be found at http://dmd.aspetjournals.org.
ABBREVIATIONS: P-gp, P-glycoprotein; MDR, multidrug resistance protein; R6G, rhodamine 6G; GF120918, N-{4[2-(1,2,3,4-tetrahydro-6,7-dimethoxy-2-isoquinolinyl)-ethyl]-phenyl}-9,10-dihydro-5-methoxy-9-oxo-4-acridine carboxamide; IPL, isolated perfused rabbit lung; IMSL, International Mathematics and Statistics Library; PS, permeability-surface area product.
Address correspondence to: Dr. David L. Roerig, Research Service 151, V.A. Medical Center, 5000 W. National Ave., Milwaukee, WI 53295. E-mail: droerig{at}mcw.edu
| References |
|---|
|
|
|---|
This article has been cited by other articles:
![]() |
B. Hosten, C. Abbara, B. Petit, A. Dauvin, F. Bourasset, R. Farinotti, P. Gonin, and L. Bonhomme-Faivre Effect of Interleukin-2 Pretreatment on Paclitaxel Absorption and Tissue Disposition after Oral and Intravenous Administration in Mice Drug Metab. Dispos., August 1, 2008; 36(8): 1729 - 1735. [Abstract] [Full Text] [PDF] |
||||
| |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||