DMD

Home Help [Feedback] [For Subscribers] [Archive] [Search] [Contents]
 QUICK SEARCH:   [advanced]


     


Drug Metabolism and Disposition Fast Forward
First published on June 4, 2007; DOI: 10.1124/dmd.107.016014


0090-9556/07/3509-1574-1579$20.00
DMD 35:1574-1579, 2007

This Article
Right arrow Abstract Freely available
Right arrow Full Text (PDF)
Right arrow All Versions of this Article:
dmd.107.016014v1
35/9/1574    most recent
Right arrow Submit a response
Right arrow Alert me when this article is cited
Right arrow Alert me when eLetters are posted
Right arrow Alert me if a correction is posted
Services
Right arrow Similar articles in this journal
Right arrow Similar articles in PubMed
Right arrow Alert me to new issues of the journal
Right arrow Download to citation manager
Citing Articles
Right arrow Citing Articles via HighWire
Right arrow Citing Articles via Google Scholar
Google Scholar
Right arrow Articles by Korjamo, T.
Right arrow Articles by Mönkkönen, J.
Right arrow Search for Related Content
PubMed
Right arrow PubMed Citation
Right arrow Articles by Korjamo, T.
Right arrow Articles by Mönkkönen, J.

Decrease in Intracellular Concentration Causes the Shift in Km Value of Efflux Pump Substrates

Timo Korjamo, Heidi Kemiläinen, Aki T. Heikkinen, and Jukka Mönkkönen

Department of Pharmaceutics, University of Kuopio, Finland

(Received April 3, 2007; accepted May 31, 2007)


    Abstract
 Top
 Abstract
 Materials and Methods
 Results
 Discussion
 References
 
Passive permeability and active efflux are parallel processes in transcellular flux. Therefore, the observed kinetics of a transporter substrate depends on both of these factors. The transporter expression has been shown to affect both the apparent Km and Vmax values. Kinetic parameters can be obtained from various experimental settings, but these do not necessarily reflect the situation in transcellular flux. Kinetic absorption models need reliable estimates of saturable kinetics when accurate in silico predictions are to be made. The effect of increasing P-glycoprotein expression on apparent transport kinetics was studied using quinidine and digoxin as model compounds. The intracellular concentrations of drugs during the transport process were also measured. A dynamic simulation model was constructed to study the observed data. The apparent Km and Vmax values increased as the P-glycoprotein expression increased. Simulations reproduced the shift in both kinetic parameters as a function of efflux pump expression. In addition, the apparent Km value showed a strong inverse relationship to the passive permeability. In contrast, the apparent Vmax value reached a maximum at intermediate passive permeability and declined above and below this passive permeability. The true Vmax and Km values were never reached. The shift in Km was assigned to a decrease in intracellular concentration at the P-glycoprotein interaction site with both experimental and simulation data. In conclusion, the apparent kinetic parameters in transcellular permeability assays depend on passive permeability and efflux pump activity. Therefore, parameters that are obtained from in vitro assays should be cautiously applied to in vivo predictions.


The determinants of oral drug bioavailability are of major interest in drug research. Absorption, distribution, metabolism, and excretion have been intensively studied, and in vitro-in vivo correlation has become an important issue in these studies. There are several applications for prediction of intestinal drug absorption from in vitro data. The compartmental absorption and transit model (Yu and Amidon, 1999Go) and other physiologically based computer models (Willmann et al., 2004Go; Usansky and Sinko, 2005Go) simulate drug absorption using estimates of permeability and solubility properties of drug molecules and physiological parameters of the gastrointestinal tract. The present models usually work well for passively absorbed drugs, but laborious in vivo experiments are required for accurate simulations on transporter substrates (Kwon et al., 2004Go).

There are several methods available for creating input parameters for these simulation models. Passive permeability can be measured relatively easily with parallel artificial membrane permeability assay or cell cultures (Hidalgo, 2001Go) or predicted from the molecular properties of the chemical (Bergström et al., 2003Go). Kinetic parameters of a transporter, e.g., P-glycoprotein (P-gp), may be evaluated based on binding to the protein (Döppenschmitt et al., 1999Go), ATP hydrolysis in microsomal preparations (Boulton et al., 2002Go), inhibition studies (Gao et al., 2001Go), or permeability experiments with suitable cell monolayers (Horie et al., 2003Go; Troutman and Thakker, 2003Go; Balakrishnan et al., 2007Go). Because the methods usually give different results, the question of which these parameters should be used in a simulation model is raised. To use the kinetic parameters from studies with isolated proteins in a simulation model, one would need a realistic estimate of the concentration at the intracellular interaction site. This is not an easy task because many processes are involved in the intracellular kinetics of drug molecules, as seen, e.g., for paclitaxel (Kuh et al., 2000Go). Monolayer studies combine the expression level of the transporter of interest and the passive transport, meaning that the kinetic values obtained vary, depending on the cell line. Furthermore, the choice of experimental parameters such as pH conditions and stirring rate influence the passive permeability of drugs and, thus, different kinetic parameters can be obtained despite the use of identical cell lines.

Some previous works described the interesting phenomenon that apparent Km and/or Vmax values in a transport experiment increase when the transporter expression increases. This result has been observed for both influx (Irie et al., 2006Go; Balakrishnan et al., 2007Go) and efflux (Zhang and Benet, 1998Go; Soldner et al., 2000Go; Horie et al., 2003Go). Recent computational analyses suggest that transport assays cannot give accurate estimates of the molecular Michaelis-Menten parameters (Bentz et al., 2005Go; Balakrishnan et al., 2007Go). A very recent theoretical work using a basic three-compartment transport model suggested that the apparent kinetics of efflux pump substrates are direction-dependent and are also affected by both passive permeability and efflux pump expression (Kalvass and Pollack, 2007Go). Therefore, there is a need to find useful relationships for kinetic values obtained from in vitro assays and their significance in observed permeability.


Figure 1
View larger version (10K):
[in this window]
[in a new window]

 
FIG. 1. An overview of the three compartmental simulation model. Virtual compound was applied to either the apical or the basolateral chamber. The mass was transferred according to equations next to the arrows. Distribution into intracellular structures was modeled with apparent cellular distribution coefficient K = Mbound/Mfree. Passive transport (Pa = 20-500 x 10-6 cm/s and Pb = 20-500 x 10-6 cm/s) proceeded bidirectionally and was driven by the concentration gradient of free molecules across the cell membrane. Active efflux (Km = 5 µM, Vmax = 0-3600 fmol/s) took its substrate from the free intracellular pool. The simulation was run up to 120 min.

 
P-gp is of considerable interest because it affects the absorption of several drugs. However, it is not known how the expression level of P-gp affects the apparent Michaelis-Menten constants in a permeability assay. We have previously described P-gp expression levels in three cell lines, Caco-2/WT < Caco-2/hPXR < MDCK-MDR1 (Korjamo et al., 2006Go). In this report, the changes in the apparent kinetic parameters of the known P-gp substrates quinidine and digoxin were studied with these cell lines expressing various amounts of P-gp. Direct experimental evidence showing that the changes in the intracellular drug concentrations lead to a shift in the apparent kinetic parameters was obtained. The effect of passive permeability on the apparent transport kinetics was investigated with a novel simulation model.


    Materials and Methods
 Top
 Abstract
 Materials and Methods
 Results
 Discussion
 References
 
Materials. Digoxin and quinidine were from Sigma-Aldrich (St. Louis, MO). [3H]Digoxin and [3H]mannitol were from PerkinElmer Life and Analytical Sciences (Boston, MA). [3H]Quinidine was from American Radiolabeled Chemicals (St. Louis, MO). The P-gp inhibitor GF120918 (Hyafil et al., 1993Go) was a kind gift from GlaxoSmithKline (Research Triangle Park, NC).

Cell Cultures and Filters. Caco-2/WT (HTB-37; American Type Culture Collection, Manassas, VA) and Caco-2/hPXR (Korjamo et al., 2005Go) cell lines were cultured as described previously (Korjamo et al., 2005Go).

MDCKII-MDR1 cells were purchased from The Netherlands Cancer Institute (Amsterdam, The Netherlands). Growth medium was Dulbecco's modified Eagle's medium (Gibco 61965-026; Invitrogen, Carlsbad, CA) with 10% heat-inactivated fetal bovine serum and 100 U/ml penicillin and 100 µg/ml streptomycin. The cells were grown in 75-cm3 flasks at +37°C in a 95% air/5% CO2 atmosphere and subcultured twice a week.

For permeability experiments, 82,000 cells/cm2 (Caco-2/WT and Caco-2/hPXR) or 390,000 cells/cm2 (MDCKII-MDR1) were seeded onto uncoated polycarbonate filters (pore size 0.4 µm, area 1.1 cm2; Corning Life Sciences, Corning, NY). Cell passage numbers were 37 to 46, 25 to 31, and 29 to 46, respectively. For both Caco-2 cell lines, the medium was changed three times per week and always on the day before the experiment. The monolayer was used for the experiments on the 21st day. For MDCKII-MDR1 cells, the growth medium was changed daily and the monolayer experiments were performed on the 4th day. Standard growth medium was used for all the cell lines.

Permeability Experiments. Apparent absorptive and secretory permeability coefficients (Papp) were determined for P-glycoprotein substrates digoxin and quinidine at various donor concentrations using radiolabeled compound (1 µl/ml) in the donor chamber as a marker. The permeability experiments were performed in both apical-to-basolateral (AB) and basolateral-to-apical (BA) directions as described previously (Korjamo et al., 2005Go). Passive permeability (Ppass) was studied at each drug concentration with assay buffer containing the P-glycoprotein inhibitor, GF120918 (2 µM). The stirring rate was 450 rpm. The experiments were conducted for up to 45 min (quinidine) or 90 min (digoxin). All assays were done in triplicate. Transepithelial electrical resistance values were measured at the beginning and at the end of the assay, and the filters with too low values [<300 {Omega}·cm2 (Caco-2/WT), <200 {Omega}·cm2 (Caco-2/hPXR), and <150{Omega}·cm2 (MDCKII-MDR1)] were discarded. Additionally, [3H]mannitol was used as a control substance for the monolayer integrity in all the assays. Mass balance was calculated as (amount in donor at last time point + amount in receiver at last time point)/amount initially added x 100%.

Unstirred Water Layer. We studied the effect of the unstirred water layer (UWL) on quinidine permeability with Caco-2/WT cells at a 0.3 µM concentration using various stirring rates (250-450 rpm). The assay was performed in both directions with and without 5 µM GF120918 at each stirring rate.

Cellular Retention. After the kinetic parameters were determined, the cellular accumulation of drugs during a permeability experiment was determined at two concentrations (below and above the observed Km value). The concentrations were quinidine Caco/WT 1 and 10 µM, Caco/hPXR and MDCK-MDR1 10 and 100 µM, and digoxin 50 µM (only because of solubility problems). Altogether 12 inserts were prepared, and the permeability experiment was started as described above. At time points 5, 15, 30, and 60 min, three inserts were sampled on both chambers. These same inserts were washed from both sides with ice-cold Hanks' balanced salt solution buffered with 25 mM Hepes. Thereafter, 300 µl of Hanks' balanced salt solution-Hepes containing 1% Triton X-100 was added, and the monolayers were lysed on an orbital shaker for 20 to 30 min. A 100-µl sample was withdrawn to liquid scintillation counting.

Sample Detection and Result Analysis. Radioactivity of the assay samples was measured with a MicroBeta liquid scintillation counter (PerkinElmer Wallac, Turku, Finland) using OptiPhase HiSafe scintillation cocktail (PerkinElmer Wallac). Apparent permeability coefficients under sink conditions were calculated for each drug concentration according to the equation

Formula(1)

where J is flux across the cell monolayer (nanomoles per second), A is area of the filter (square centimeters), and C0 is the initial donor concentration (micromolar). Each drug concentration was plotted against the flux (J), and Michaelis-Menten kinetic parameter values Vmax and Km were calculated from these curves according to the equation

Formula(2)

where Ppass is the passive permeability in the presence of P-gp inhibitor GF120918 (centimeters per second), A is the area of the filter (square centimeters), and C0 is the initial donor concentration (micromolar). GraphPad Prism software (GraphPad Software Inc., San Diego, CA) was applied for curve fitting. Weighting of 1/concentration was used in the fitting procedure.

Simulations. A novel three-compartment model (Fig. 1) was constructed with Stella 8.1.4 software (isee systems, Lebanon, NH) to simulate mass flow in a permeability experiment. Three-compartment models have been used to simulate in vitro permeability experiments (Ito et al., 1999Go; González-Alvarez et al., 2005Go). We introduce an apparent cellular distribution coefficient (K = Mbound/Mfree) to the model. Mfree represents the drug amount in cytosol that is ready to diffuse across cell membranes, and Mbound is the amount of drug that is bound to intracellular structures, which leads to high cellular drug accumulation. This approach can accurately simulate the time course of intracellular accumulation of drug during a permeability experiment (unpublished results). The flux across the basolateral membrane was assumed to be purely passive. In addition to the passive component, an efflux transporter obeying regular Michaelis-Menten kinetics affected the flux across the apical membrane. A virtual efflux pump substrate with high cellular retention (K = 200, corresponding to a fitted value for quinidine, unpublished results) was created. The Km value was fixed to 5 µM. The volumes of the compartments were 0.5 ml (apical), 1.5 ml (basolateral), and 0.0022 ml (cellular assuming cell height of 20 µm). The virtual compounds accumulated into cells according to apparent cellular distribution coefficient. The efflux pump used the "free" cellular concentration as its substrate. The true Vmax and membrane permeabilities were varied 26-fold (80-2100 fmol/s) and 5-fold (50-250 x 10-6 cm/s), respectively. Each simulation was run at eight initial concentrations (0.1-300 µM). For each parameter combination, the simulation with Vmax = 0 was considered as an efflux pump-inhibited experiment. For each simulation, the highest flux (1-min intervals) to the receiver chamber was searched. The initial applied donor concentration was used for calculations. The passive flux (true Vmax = 0) was subtracted from the total flux in the presence of efflux activity. This "active flux" data were then fitted to regular Michaelis-Menten equation

Formula(3)

where A is the area of the filter (square centimeters) and C0 is the initial donor concentration (micromolar). GraphPad Prism software was applied for curve fitting. Weighting could be omitted because simulation data lack experimental error.


Figure 2
View larger version (12K):
[in this window]
[in a new window]

 
FIG. 2. The apparent permeability of quinidine across Caco/WT cells at different stirring rates in the AB ({triangleup}, {blacktriangleup}) and BA ({square}, {blacksquare}) direction. Experiments were conducted both with (open symbols) and without (filled symbols) a 5 µM concentration of the P-glycoprotein inhibitor GF120918. Symbols represent average values of three filters and bars represent S.D.

 
Statistical Testing. The difference in the intracellular accumulation between P-gp inhibited and uninhibited wells at every time point was tested with Student's t test. The analysis was performed with GraphPad Prism software.


    Results
 Top
 Abstract
 Materials and Methods
 Results
 Discussion
 References
 
Permeability. The passive permeability of quinidine was unpolarized and of similar magnitude in all cell lines (Table 1). The mass balance of quinidine was 75 to 80% in the AB direction and 87 to 90% in the BA direction. The Ppass of digoxin was moderately lower in the MDCK-MDR1 cells than in the Caco-2 cell lines, and the mass balance was practically 100%. P-gp-inhibited digoxin permeability decreased at concentrations 200 µM in all cell lines in both the AB and BA directions (data not shown), and permeability values obtained at higher concentrations were not used for curve-fitting.


View this table:
[in this window]
[in a new window]

 
TABLE 1 Passive permeability and kinetic parameters of quinidine and digoxin

The apparent permeabilities of quinidine and digoxin across cell monolayers were determined at a stirring rate of 450 rpm. Passive permeability was determined in the presence of a 5 µM concentration of the P-glycoprotein inhibitor GF120918. Total fluxes at different concentrations were calculated, and the results were fitted to eq. 2. Permeabilities are expressed as means ± S.D. and kinetic parameters as means ± S.E.

 

The unstirred water layer forms a considerable permeability barrier for compounds that traverse the cell monolayer very rapidly (Karlsson and Artursson, 1991Go). Therefore, we determined the Papp values at a low quinidine concentration (0.3 µM) in both directions at various stirring rates. The passive permeability in both directions increased as the stirring became more vigorous (Fig. 2). The uninhibited transport, however, was increased only in the BA direction. Thus, the apparent P-gp-mediated transport (Papp - Ppass) was relatively constant in the BA direction but increased considerably in the AB direction at higher stirring rates. Therefore, the kinetic studies were conducted at the highest stirring rate (450 rpm). The apparent permeability of digoxin did not change at different stirring speeds (data not shown).

Cellular Retention. P-glycoprotein expression affected intracellular drug concentrations. The P-gp inhibitor GF120918 increased the cellular retention of quinidine in Caco/WT at early time points at 1 µM concentration (Fig. 3a). In contrast, no difference in the cellular quinidine amount between inhibited and uninhibited cells was observed when the concentration was increased well above the apparent Km values (10 µM). Caco/hPXR and MDCK-MDR1 cells have sufficient P-gp expression to maintain decreased cellular quinidine levels at 10 µM (Fig. 3, b and c). However, a considerable increase (to 100 µM) in donor concentration removes the difference in cellular retention between inhibited and uninhibited wells. The results were similar in the BA direction even though the time courses of accumulation were somewhat different (data not shown). The overall mass balance (apical + basolateral + cellular) remained constant after initial drop of 3 to 9%, which is probably due to binding to the apparatus (Palmgrén et al., 2006Go) or small pipetting errors. Similar observations on cellular accumulation were made in the simulation experiments when the efflux pump expression level was altered (Fig. 4).


Figure 3
View larger version (23K):
[in this window]
[in a new window]

 
FIG. 3. The accumulation of quinidine into Caco/WT (a), Caco/hPXR (b), and MDCK-MDR1 (c) cell monolayers during an apical-to-basolateral permeability experiment. At each indicated time point, the cell monolayers were washed with ice-cold buffer and lysed with 1% Triton X-100. Drug amount was determined from the lysate and normalized to the initial drug amount in the donor. Initial substrate concentrations were 1 ({circ}, bullet), 10 ({triangleup}, {blacktriangleup}), and 100 ({square}, {blacksquare}) µM. The experiments were performed in the absence (filled symbols, solid lines) or presence (open symbols, dashed lines) of a 5 µM concentration of the P-glycoprotein inhibitor GF120918. Symbols represent average values of 3 to 12 filters. Error bars were omitted for clarity but the coefficient of variation was generally <10%. The differences in intracellular quinidine amounts between inhibited and uninhibited wells were statistically tested with Student's t test at each time point. Statistical significance is indicated as n.s. = p > 0.05, *, p < 0.05, **, p < 0.01, ***, p < 0.001. The symbols at the top of each graph refer to lower and symbols at the bottom to higher initial quinidine concentration.

 

Figure 4
View larger version (16K):
[in this window]
[in a new window]

 
FIG. 4. The cellular accumulation of a quinidine-like virtual compound (Pa = Pb = 250 x 10-6 cm/s, Kbound/free = 200) during an apical-to-basolateral simulation at 10 µM initial donor concentration. The total cellular amount of the virtual compound after each time interval was normalized to the initial donor amount. The true Vmax of the efflux pump was set to 0 ({square}), 400 ({blacksquare}), 1800 ({blacktriangleup}), or 3600 (bullet) fmol/s, which corresponded to the same fold differences of P-gp expression that were observed in our cell lines (Table 2).

 
Cellular retention of digoxin was clearly lower than that of quinidine (Fig. 5), and the mass balance was nearly perfect. P-gp could reduce the cellular accumulation to a similar extent in the two Caco-2 cell lines. However, MDCK-MDR1 cells retained much less digoxin than the Caco-2 cells, and P-gp could only modestly decrease its accumulation.


Figure 5
View larger version (17K):
[in this window]
[in a new window]

 
FIG. 5. The accumulation of digoxin (50 µM) into Caco/WT ({diamond}, {diamondsuit}), Caco/hPXR ({triangleup}. {blacktriangleup}), and MDCK-MDR1 ({square}, {blacksquare}) cells during a basolateral-to-apical permeability experiment. At each indicated time point, the cell monolayers were washed with ice-cold buffer and lysed with 1% Triton X-100. The drug amount was determined from the lysate and normalized to the initial drug amount in the donor. The experiments were performed in the absence (filled symbols, solid lines) or presence (open symbols, dashed lines) of 5 µM concentration of the P-glycoprotein inhibitor GF120918. Symbols represent average values of 3 to 12 filters. Error bars were omitted for clarity, but the coefficient of variation was generally <10%. The difference in intracellular digoxin amounts between inhibited and uninhibited wells was statistically tested with Student's t test. Intracellular accumulation was statistically lower in the absence than in the presence of GF120918 at all time points (p < 0.05).

 
Experimental Kinetic Parameters. The apparent kinetic transport parameters are presented in Table 1. The Km and Vmax values increased in the order Caco/WT < Caco/hPXR < MDCK-MDR1 with both substrates. However, parameters obtained for digoxin from MDCK-MDR1 cells should be treated with caution because the Km value is 5-fold higher than the highest experimental concentration. A practically linear relationship was obtained when Km and Vmax values of quinidine in the AB direction are plotted against each other (Fig. 6). A similar but not equally linear relationship was seen for both drug molecules in the BA direction. Some literature data are also presented for comparison in these figures.


Figure 6
View larger version (17K):
[in this window]
[in a new window]

 
FIG. 6. Correlation of the apparent kinetic parameters (Km and Vmax) of quinidine (a) and digoxin (b) in cells expressing different amounts of P-glycoprotein. Figures include our experimental data in the AB ({blacksquare}) and/or BA ({blacktriangleup}) direction and literature data ({square}, {triangleup}) from references (Stephens et al., 2001Go; Lowes et al., 2003Go; Troutman and Thakker, 2003Go).

 

The total P-gp levels in the cell lines were obtained from the background corrected optical density values of the Western blot figure published recently (Korjamo et al., 2006Go). Table 2 presents the protein levels and kinetic parameters in relation to Caco/WT cells. The kinetic parameters of quinidine in the AB direction correlate almost perfectly to the protein data. Other apparent kinetic parameters of Caco/hPXR cells change less than expected from the protein data. The Km values from MDCK-MDR1 cells show consistency with protein data, but there is more variation in Vmax.


View this table:
[in this window]
[in a new window]

 
TABLE 2 Comparison of P-gp protein level and kinetic parameters with those in Caco/WT cells

 

Simulations. The apparent Km and Vmax values both increased as the efflux pump expression increased. This shift in the apparent Km value increased as the passive permeability decreased (Fig. 7a). At low passive permeability, the apparent Vmax was low (Fig. 7b). When the permeability increased to medium levels, the apparent Vmax reached its peak value. Finally, apparent Vmax decreased again at high passive permeability. Similar to the experimental observations, an increase in efflux pump expression decreased intracellular drug concentration (Fig. 4). The apparent Km was always considerably higher in the AB than in the BA direction (Table 3).


Figure 7
View larger version (22K):
[in this window]
[in a new window]

 
FIG. 7. The effect of efflux pump expression and passive permeability on the apparent kinetic parameters of a virtual compound (Kbound/free = 200) in apical-to-basolateral simulations at different expression levels of efflux pump. a, the apparent Km is plotted against the true Vmax (80-2100 fmol/s). Membrane permeabilities were set to equal values of 20 (x), 50 ({square}), 125 ({triangleup}), 250 ({diamond}) or 500 ({circ}) x 10-6 cm/s. b, the apparent Vmax is plotted against passive membrane permeability (20-500 x 10-6 cm/s). The true Vmax values were 80 ({diamond}), 400 ({circ}), 1250 ({triangleup}), or 2100 ({square}) fmol/s.

 

View this table:
[in this window]
[in a new window]

 
TABLE 3 Apparent kinetic parameters of a quinidine-like (Pa = Pb = 250 x 10–6 cm/s, K = 200) virtual compound in AB and BA directions at various efflux pump expression levels

 


    Discussion
 Top
 Abstract
 Materials and Methods
 Results
 Discussion
 References
 
The apparent transport kinetics of both influx and efflux transporter substrates depend on the access of the substrate to the interaction site. The increase in transporter expression leads to biased apparent kinetic parameters in monolayer transport assays. These observed shifts have been hypothesized to result either from a decreased concentration in the UWL next to the plasma membrane (Winne, 1977Go; Balakrishnan et al., 2007Go) or from a decreased intracellular concentration (Bentz et al., 2005Go). The previous works have based their hypotheses mainly on computational data. This study is the first that shows direct experimental evidence of how the increase in efflux pump expression decreases the intracellular drug amount, which then leads to larger shifts in the apparent kinetic parameters. Our novel simulation model can reproduce the time course of intracellular accumulation. Simulations show that the passive membrane permeability altered the bias in kinetic parameters. The shift in the apparent Km value is inversely proportional to the passive permeability. In contrast, the apparent Vmax first increases and then decreases as the passive permeability increases.

Digoxin permeabilities up to 500 µM has been reported (Troutman and Thakker, 2003Go). However, we detected that the Ppass decreased above 200 µM (data not shown), probably due to precipitation of the drug or saturation of some uptake transporter. Indeed, organic anion transporting polypeptide E has been shown to transport digoxin (Mikkaichi et al., 2004Go). However, its expression is low both in Caco-2 (Hilgendorf et al., 2007Go) and MDCK cells (Goh et al., 2002Go). UWL affects the kinetic parameters of absorptive transporters (Winne, 1977Go; Balakrishnan et al., 2007Go). It has been hypothesized but not experimentally proven that the apparent Km value increases because the concentration next to plasma membrane decreases as the UWL thickness increases. In the present quinidine experiments, the UWL affected the Papp in the BA direction and the Ppass in both directions (Fig. 2), but the Papp in the AB direction remained fairly constant. Therefore, the UWL influences the kinetic parameters also for efflux transporter substrates. The kinetic parameters in the AB direction may be more sensitive to changes in the stirring rate than those in the BA direction because the P-gp effect (Papp - Ppass) changes considerably in the AB but not in the BA direction. Digoxin permeability and, consequently, the kinetic parameters are fairly insensitive to the stirring rate because its transport is not rate-limited by the UWL.

The apparent kinetic values (Table 1) show that both Km and Vmax increase as the P-gp activity increases and that these values are strongly correlated (Fig. 6). Moreover, there is often a good correlation between the kinetic parameters and the expression of P-gp (Table 2). A plausible explanation for this is that P-gp reduces the drug concentration at the interaction site below saturation concentrations. Therefore, the donor concentration does not reflect the concentration at the P-gp interaction site. Indeed, the total cellular amounts of quinidine (Fig. 3) and digoxin (Fig. 5) are smaller in the absence than in the presence of the P-gp inhibitor when the donor concentration does not saturate P-gp. In addition, the simulation model reproduced the experimental data for quinidine (Fig. 4). Together these data confirm the hypothesis that the reduced intracellular concentration causes the shift in observed kinetic parameters in monolayer permeability assays. The actual concentration at the interaction site may be even more reduced if the efflux pump is able to create a concentration gradient inside the cell.

Recent results from a simulation model suggest that apparent Vmax of P-gp from permeability experiments follows the molecular Vmax and, thus, the protein expression level reasonably well (Bentz et al., 2005Go). However, the apparent Km value is always expected to be higher than the molecular Km, and it was shown to depend on the passive permeability and the maximal pumping capacity of P-gp. Supporting this result, our simulations show that the decrease in passive permeability increases the observed Km shift when efflux pump expression increases (Fig. 7a). In contrast, the apparent Vmax value at a particular transporter expression level increases up to a certain passive permeability and starts to decline thereafter (Fig. 7b). The apparent Vmax never reached the true Vmax but only attained <50% of this value. The decrease in apparent Vmax at low passive permeability may be caused by the combination of slow entrance and high accumulation into cellular structures. On the other hand, the high passive permeability masks the active transport and the apparent Vmax decreases. Thus, the parallel active and passive transport processes prevent obtaining of true kinetic values of transporters, at least with simple approaches (eqs. 2 and 3) that are used to treat such permeability data. Recently a new kinetic approach for calculating the parameters in transport experiments has been suggested (Kalvass and Pollack, 2007Go). Nevertheless, the present results give an idea of how different parameters affect the apparent transport kinetics viewed from the donor solution. Actually, the UWL and the plasma membrane both form a diffusion barrier to the interaction site, and the Km shift is expected to be more pronounced for efflux than for influx transporters with an extracellular interaction site (Balakrishnan et al., 2007Go).

The relationship between P-gp expression and kinetic parameters is not simple, and it can vary, depending on the substrate and transport direction (Table 2). A previous study reported that the apparent Km value is always higher in the AB than in the BA direction (Troutman and Thakker, 2003Go). This observation was explained by rigidity of the apical membrane, which results in lower permeability than at the basolateral membrane. In this study, digoxin behaved similarly, because the Km value in the AB direction is much higher than 200 µM (Table 1). The fact that our simulations with equal membrane permeabilities gave higher Km values in the AB than in the BA direction (Table 3) suggests that the observed difference may partly be an artifact caused by the Transwell apparatus. Because the donor and receiver volumes are unequal and small, the donor concentration decreases more rapidly in the AB than in the BA direction. Thus, the cellular concentrations remain lower in the AB than in the BA direction at any given donor concentrations, and the Km shift increases. Moreover, the donor concentration is considerably below the initial concentration at steady state in which the maximal flux is taken. The directional difference in the apparent Km value was practically abolished when both apical and basolateral chambers were increased to 5000 ml in simulations (data not shown). Now the sink conditions could be maintained, and the donor concentration remained at initial level. The present experimental data suggest that the apparent Km value may also be higher in the secretory than in the absorptive direction (quinidine) (Table 1). This result requires that the resistance at the apical membrane be considerably lower than that at the basolateral membrane. Our hypothesis is that the applied vigorous stirring decreases the total passive resistance (UWL + cell membrane) at the apical membrane to a lower value than at the basolateral membrane. This hypothesis needs further experimental evidence.

In conclusion, The apparent kinetic parameters shift when a diffusion barrier is placed between the donor and the interaction site. The access to the interaction site may be very different in simple laboratory experiments and in whole organisms. Therefore, scaling in vitro parameters to permeability at different locations in vivo is not straightforward. Further studies are needed to understand the relationships between permeability barriers and observed kinetics. Ultimately, a predictive permeability model may be constructed using simple in vitro parameters reflecting the passive permeability, intracellular distribution, and transporter affinity.


    Acknowledgments
 
Technical assistance from Paula Nyyssönen and Markku Taskinen is gratefully acknowledged.


    Footnotes
 
This study was financially supported by the Finnish Funding Agency for Technology and Innovation.

doi:10.1124/dmd.107.016014.

ABBREVIATIONS: P-gp, P-glycoprotein; WT, wild type; hPXR, human pregnane X receptor; 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; AB, apical-to-basolateral; BA, basolateral-to-apical; UWL, unstirred water layer.

Address correspondence to: Dr. Timo Korjamo, Department of Pharmaceutics, University of Kuopio, Yliopistonranta 1C, FI-70210, Kuopio, Finland. E-mail: Timo.Korjamo{at}uku.fi


    References
 Top
 Abstract
 Materials and Methods
 Results
 Discussion
 References
 


Balakrishnan A, Hussainzada N, Gonzalez P, Bermejo M, Swaan PW, and Polli JE (2007) Bias in estimation of transporter kinetic parameters from over-expression systems: interplay of transporter expression level and substrate affinity. J Pharmacol Exp Ther 320: 133-144.[Abstract/Free Full Text]

Bentz J, Tran TT, Polli JW, Ayrton A, and Ellens H (2005) The steady-state Michaelis-Menten analysis of P-glycoprotein mediated transport through a confluent cell monolayer cannot predict the correct Michaelis constant Km. Pharm Res 22: 1667-1677.[CrossRef][Medline]

Bergström CA, Strafford M, Lazorova L, Avdeef A, Luthman K, and Artursson P (2003) Absorption classification of oral drugs based on molecular surface properties. J Med Chem 46: 558-570.[CrossRef][Medline]

Boulton DW, DeVane CL, Liston HL, and Markowitz JS (2002) In vitro P-glycoprotein affinity for atypical and conventional antipsychotics. Life Sci 71: 163-169.[CrossRef][Medline]

Döppenschmitt S, Langguth P, Regardh CG, Andersson TB, Hilgendorf C, and Spahn-Langguth H (1999) Characterization of binding properties to human P-glycoprotein: development of a [3H]verapamil radioligand-binding assay. J Pharmacol Exp Ther 288: 348-357.[Abstract/Free Full Text]

Gao J, Murase O, Schowen RL, Aube J, and Borchardt RT (2001) A functional assay for quantitation of the apparent affinities of ligands of P-glycoprotein in Caco-2 cells. Pharm Res 18: 171-176.[CrossRef][Medline]

Goh LB, Spears KJ, Yao D, Ayrton A, Morgan P, Roland Wolf C, and Friedberg T (2002) Endogenous drug transporters in in vitro and in vivo models for the prediction of drug disposition in man. Biochem Pharmacol 64: 1569-1578.[CrossRef][Medline]

González-Alvarez I, Fernandez-Teruel C, Garrigues TM, Casabo VG, Ruiz-Garcia A, and Bermejo M (2005) Kinetic modelling of passive transport and active efflux of a fluoroquinolone across Caco-2 cells using a compartmental approach in NONMEM. Xenobiotica 35: 1067-1088.[CrossRef][Medline]

Hidalgo IJ (2001) Assessing the absorption of new pharmaceuticals. Curr Top Med Chem 1: 385-401.[CrossRef][Medline]

Hilgendorf C, Ahlin G, Seithel A, Artursson P, Ungell AL, and Karlsson JE (2007) Expression of thirty-six drug transporter genes in human intestine, liver, kidney, and organotypic cell lines. Drug Metab Dispos 35: 1334-1341.

Horie K, Tang F, and Borchardt RT (2003) Isolation and characterization of Caco-2 subclones expressing high levels of multidrug resistance protein efflux transporter. Pharm Res 20: 161-168.[CrossRef][Medline]

Hyafil F, Vergely C, Du Vignaud P, and Grand-Perret T (1993) In vitro and in vivo reversal of multidrug resistance by GF120918, an acridonecarboxamide derivative. Cancer Res 53: 4595-4602.[Abstract/Free Full Text]

Irie M, Terada T, Tsuda M, Katsura T, and Inui K (2006) Prediction of glycylsarcosine transport in Caco-2 cell lines expressing PEPT1 at different levels. Pflugers Arch 452: 64-70.[CrossRef][Medline]

Ito S, Woodland C, Sarkadi B, Hockmann G, Walker SE, and Koren G (1999) Modeling of P-glycoprotein-involved epithelial drug transport in MDCK cells. Am J Physiol 277: F84-F96.[Medline]

Kalvass JC and Pollack GM (2007) Kinetic considerations for the quantitative assessment of efflux activity and inhibition: implications for understanding and predicting the effects of efflux inhibition. Pharm Res 24: 265-276.[CrossRef][Medline]

Karlsson J and Artursson P (1991) A method for the determination of cellular permeability coefficients and aqueous boundary layer thickness in monolayers of intestinal epithelial (Caco-2) cells grown in permeable filter chambers. Int J Pharm 71: 55-64.[CrossRef]

Korjamo T, Honkakoski P, Toppinen MR, Niva S, Reinisalo M, Palmgrén JJ, and Mönkkönen J (2005) Absorption properties and P-glycoprotein activity of modified Caco-2 cell lines. Eur J Pharm Sci 26: 266-279.[CrossRef][Medline]

Korjamo T, Mönkkönen J, Uusitalo J, Turpeinen M, Pelkonen O, and Honkakoski P (2006) Metabolic and efflux properties of Caco-2 cells stably transfected with nuclear receptors. Pharm Res 23: 1991-2001.[CrossRef][Medline]

Kuh HJ, Jang SH, Wientjes MG, and Au JL (2000) Computational model of intracellular pharmacokinetics of paclitaxel. J Pharmacol Exp Ther 293: 761-770.[Abstract/Free Full Text]

Kwon H, Lionberger RA, and Yu LX (2004) Impact of P-glycoprotein-mediated intestinal efflux kinetics on oral bioavailability of P-glycoprotein substrates. Mol Pharm 1: 455-465.[CrossRef][Medline]

Lowes S, Cavet ME, and Simmons NL (2003) Evidence for a non-MDR1 component in digoxin secretion by human intestinal Caco-2 epithelial layers. Eur J Pharmacol 458: 49-56.[CrossRef][Medline]

Mikkaichi T, Suzuki T, Onogawa T, Tanemoto M, Mizutamari H, Okada M, Chaki T, Masuda S, Tokui T, Eto N, et al. (2004) Isolation and characterization of a digoxin transporter and its rat homologue expressed in the kidney. Proc Natl Acad SciUSA 101: 3569-3574.[Abstract/Free Full Text]

Palmgrén JJ, Monkkonen J, Korjamo T, Hassinen A, and Auriola S (2006) Drug adsorption to plastic containers and retention of drugs in cultured cells under in vitro conditions. Eur J Pharm Biopharm 64: 369-378.[CrossRef][Medline]

Soldner A, Benet LZ, Mutschler E, and Christians U (2000) Active transport of the angiotensin-II antagonist losartan and its main metabolite EXP 3174 across MDCK-MDR1 and Caco-2 cell monolayers. Br J Pharmacol 129: 1235-1243.[CrossRef]

Stephens RH, O'Neill CA, Warhurst A, Carlson GL, Rowland M, and Warhurst G (2001) Kinetic profiling of P-glycoprotein-mediated drug efflux in rat and human intestinal epithelia. J Pharmacol Exp Ther 296: 584-591.[Abstract/Free Full Text]

Troutman MD and Thakker DR (2003) Efflux ratio cannot assess P-glycoprotein-mediated attenuation of absorptive transport: asymmetric effect of P-glycoprotein on absorptive and secretory transport across Caco-2 cell monolayers. Pharm Res 20: 1200-1209.[CrossRef][Medline]

Usansky HH and Sinko PJ (2005) Estimating human drug oral absorption kinetics from Caco-2 permeability using an absorption-disposition model: model development and evaluation and derivation of analytical solutions for ka and Fa. J Pharmacol Exp Ther 314: 391-399.[Abstract/Free Full Text]

Willmann S, Schmitt W, Keldenich J, Lippert J, and Dressman JB (2004) A physiological model for the estimation of the fraction dose absorbed in humans. J Med Chem 47: 4022-4031.[CrossRef][Medline]

Winne D (1977) Correction of the apparent Michaelis constant, biased by an unstirred layer, if a passive transport component is present. Biochim Biophys Acta 464: 118-126.[Medline]

Yu LX and Amidon GL (1999) A compartmental absorption and transit model for estimating oral drug absorption. Int J Pharm 186: 119-125.[CrossRef][Medline]

Zhang Y and Benet LZ (1998) Characterization of P-glycoprotein mediated transport of K02, a novel vinylsulfone peptidomimetic cysteine protease inhibitor, across MDR1-MDCK and Caco-2 cell monolayers. Pharm Res 15: 1520-1524.[CrossRef][Medline]


This article has been cited by other articles:


Home page
Drug Metab. Dispos.Home page
P. Acharya, M. P. O'Connor, J. W. Polli, A. Ayrton, H. Ellens, and J. Bentz
Kinetic Identification of Membrane Transporters That Assist P-glycoprotein-Mediated Transport of Digoxin and Loperamide through a Confluent Monolayer of MDCKII-hMDR1 Cells
Drug Metab. Dispos., February 1, 2008; 36(2): 452 - 460.
[Abstract] [Full Text] [PDF]


This Article
Right arrow Abstract Freely available
Right arrow Full Text (PDF)
Right arrow All Versions of this Article:
dmd.107.016014v1
35/9/1574    most recent
Right arrow Submit a response
Right arrow Alert me when this article is cited
Right arrow Alert me when eLetters are posted
Right arrow Alert me if a correction is posted
Services
Right arrow Similar articles in this journal
Right arrow Similar articles in PubMed
Right arrow Alert me to new issues of the journal
Right arrow Download to citation manager
Citing Articles
Right arrow Citing Articles via HighWire
Right arrow Citing Articles via Google Scholar
Google Scholar
Right arrow Articles by Korjamo, T.
Right arrow Articles by Mönkkönen, J.
Right arrow Search for Related Content
PubMed
Right arrow PubMed Citation
Right arrow Articles by Korjamo, T.
Right arrow Articles by Mönkkönen, J.


Home Help [Feedback] [For Subscribers] [Archive] [Search] [Contents]
All ASPET Journals Molecular Pharmacology Pharmacological Reviews
 Molecular Interventions Drug Metabolism and Disposition