QUICK SEARCH:   [advanced] Author: Keyword(s): Year:  Vol:  Page:

This Article Abstract Full Text (PDF) All Versions of this Article: dmd.107.017301v1 36/2/452    most recent Submit a response Alert me when this article is cited Alert me when eLetters are posted Alert me if a correction is posted Services Similar articles in this journal Similar articles in PubMed Alert me to new issues of the journal Download to citation manager Citing Articles Citing Articles via HighWire Citing Articles via Google Scholar Google Scholar Articles by Acharya, P. Articles by Bentz, J. Search for Related Content PubMed PubMed Citation Articles by Acharya, P. Articles by Bentz, J.

Kinetic Identification of Membrane Transporters That Assist P-glycoprotein-Mediated Transport of Digoxin and Loperamide through a Confluent Monolayer of MDCKII-hMDR1 Cells

Department of Bioscience & Biotechnology, Drexel University, Philadelphia, Pennsylvania (P.A., M.P.O., J.B.); Preclinical Drug Metabolism and Pharmacokinetics, GlaxoSmithKline, King of Prussia, Pennsylvania (P.A., H.E.), Preclinical Drug Metabolism and Pharmacokinetics, GlaxoSmithKline, Research Triangle Park, North Carolina (J.W.P.); and Preclinical Drug Metabolism and Pharmacokinetics, GlaxoSmithKline, Welwyn, England (A.A.)

(Received June 29, 2007; Accepted October 23, 2007)

	Abstract

Top Abstract Materials and Methods Results Discussion References

 
A robust screen for compound interaction with P-glycoprotein (P-gp) has some obvious requirements, such as a cell line expressing P-gp and a probe substrate that is transported solely by P-gp and passive permeability. It is actually difficult to prove that a particular probe substrate interacts only with P-gp in the chosen cell line. Using a confluent monolayer of MDCKII-hMDR1 cells, we have determined the elementary rate constants for the P-gp efflux of amprenavir, digoxin, loperamide, and quinidine. For amprenavir and quinidine, transport was fitted with just P-gp and passive permeability. For digoxin and loperamide, fitting required a basolateral transporter (p < 0.01), which was inhibited by the P-gp inhibitor 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 (GF120918). This means that when digoxin is used as a probe substrate and a compound is shown to inhibit digoxin flux, it could be that the inhibition occurs at the basolateral transporter rather than at P-gp. Digoxin basolateral>apical efflux also required an apical importer (p < 0.05). We propose that amprenavir and quinidine are robust probe substrates for assessing P-gp interactions using the MDCKII-hMDR1 confluent cell monolayer. Usage of another cell line, e.g., LLC-hMDR1 or Caco-2, would require the same kinetic validation to ensure that the probe substrate interacts only with P-gp. Attempts to identify the additional digoxin and loperamide transporters using a wide range of substrates/inhibitors of known epithelial transporters (organic cation transporters, organic anion transporters, organic ion-transporting polypeptide, uric acid transporter, or multidrug resistance-associated protein) failed to inhibit the digoxin or loperamide transport through their basolateral transporter.

The human multidrug resistance transporter P-glycoprotein (P-gp) (Juliano and Ling, 1976) is the product of the hMDR1 (ABCB1) gene and is widely expressed in human epithelial tissue as a protection against xenobiotics (Dean et al., 2001). Polarized confluent cell monolayers overexpressing P-gp have been used extensively as a model system to study P-gp transport mechanisms and to assess the risk of P-gp-mediated drug-drug interactions (Tang et al., 2002; Rautio et al., 2006; Bartholomé et al., 2007; Korjamo et al., 2007).

Using the polarized Madin-Darby canine kidney II cell line that overexpresses human MDR1 (MDCKII-hMDR1) confluent cell monolayer, which overexpresses human P-gp in the apical plasma membrane and a detailed mass action analysis of P-gp efflux kinetics, we have been able to answer several fundamental questions about P-gp function (Bentz et al., 2005; Tran et al., 2005; Acharya et al., 2006). With P-gp as the only transporter in the kinetic model, we have found that amprenavir and quinidine were well fitted, whereas loperamide showed more efflux than could be explained by just P-gp (Tran et al., 2005; Acharya et al., 2006). Previous work had shown other transporters in model cell lines. Secretory transport of rhodamine 123 across Caco-2 cell monolayers involves a basolateral uptake transporter in addition to P-gp (Troutman and Thakker, 2003). MDCKII wild-type cells and MDCKII cells overexpressing bcrp1 contain an apical influx transporter for the BCRP substrate mitoxantrone (Pan and Elmquist, 2007).

These observations led 1) to the speculation that there was another transporter involved for loperamide in the MDCKII-hMDR1 cells and 2) to the question of the extent to which other transporters in these cells might affect inhibition of a probe substrate in screening for drug-drug interactions with P-gp. A "robust" probe substrate for P-gp inhibition screening should be transported across the cell monolayer solely by P-gp and the passive permeability, which is measured in the presence of a potent P-gp inhibitor, such as GF120918. In this case, the P-gp-mediated efflux can be roughly estimated by the difference between the two curves, i.e., the total efflux (–GF120918) – the passive efflux (+GF120918). We say "roughly" because backflow, from the receiver side to the donor side, must be ignored (Tran et al., 2005).

If a probe substrate also uses another transporter, then there are two possibilities. The first case is easy. If this new transporter is not inhibited by GF120918, then the +GF120918 efflux would contain both the contributions from the other transporter and from passive permeability. In this case, the difference between the two curves would still be due to P-gp alone. One just needs to be clear that "passive permeability," +GF120918, measurements can include transporters. The second case is the problem. If the novel transporter is inhibited by GF120918, then the passive permeability, +GF120918, will underestimate the amount of probe substrate actually reaching P-gp in the –GF120918 experiment. Then the difference between the curves yields an overestimate of the P-gp-mediated efflux. Now, if a novel compound binds to and inhibits the other transporter, –GF120918, then less probe substrate will reach P-gp, thus reducing its efflux kinetics. Such a reduction would probably be misinterpreted as the novel compound being a P-gp inhibitor.

Here we have used our kinetic analysis to detect the presence of novel transporters of loperamide and digoxin in the basolateral membrane of this cell line, as well as a digoxin transporter in the apical membrane. Digoxin has been used as a probe substrate with MDCKII-hMDR1 confluent cell monolayers to screen novel compounds for P-gp inhibition (Rautio et al., 2006). A rigorous statistical analysis confirms the need for these additional transporters and thus supports our kinetic findings. For amprenavir and quinidine, only P-gp and passive efflux (–GF120918) are required to fit the data, suggesting that they are robust probe substrates for this cell line.

	Materials and Methods

Top Abstract Materials and Methods Results Discussion References

 
Materials. Amprenavir and GF120918 were from GlaxoSmithKline Pharmaceuticals (Research Triangle Park, NC), and loperamide and quinidine were from Sigma-Aldrich (St. Louis, MO). [³H]Loperamide (10 Ci/mmol) and [³H]amprenavir (21 Ci/mmol) were custom synthesized by GE Healthcare (Little Chalfont, Buckinghamshire, UK). [³H]Quinidine (20 Ci/mmol) was from MP Biomedicals, Inc. (Solon, OH). Dulbecco's modified Eagle's medium (DMEM) was from MediaTech, VWR International, Inc. (West Chester, PA). DMEM with 25 mM HEPES buffer, high glucose (4.5 g/l), L-glutamine, fetal bovine serum, pyridoxine hydrochloride, without sodium pyruvate, and with phenol red was from Invitrogen (Carlsbad, CA). The same medium without phenol red, denoted transport media, was used for transport experiments. Transwell 12-well plates with polycarbonate inserts (0.4-µm pore size and 12-mm diameter) were obtained from Corning Life Sciences (Acton, MA). All compounds tested in Table 2 were from Sigma Chemical.

Cell Line and Culture Conditions. The MDCKII-hMDR1 was purchased from the Netherlands Cancer Institute (Amsterdam, The Netherlands). MDCKII-hMDR1 cells were grown in 175-cm² culture flasks using DMEM with 10% fetal bovine serum, 1% L-glutamine, and 50 U/ml penicillin and 50 mg/ml streptomycin at 37°C in 5% CO₂ atmosphere. Cells were split twice a week at 70 to 80% confluence in a ratio of 1:40, after at least two washes in phosphate-buffered saline and trypsinization with 0.25% trypsin-EDTA. All transport assays were done with cells from passages 30 to 55. Cells were kept at 37°C in 5% CO₂.

Single Substrate Efflux Assay. Cells were seeded in 12-well Costar transwell plates with polycarbonate membrane inserts at a density of 175,000 cells per insert and grown for 4 days in culture medium. Cells were given fresh media 1 day after seeding. On the day of the experiment, culture medium was removed, and cells were preincubated for 30 min with either transport media alone or transport media with 2 µM GF120918, an inhibitor of P-gp. Efflux of amprenavir, digoxin, loperamide, and quinidine across the confluent monolayer of cells was measured in both directions, i.e., apical to basolateral (A>B) and basolateral to apical (B>A) in the presence and absence of GF120918. [³H]Amprenavir, [³H]digoxin, [³H]loperamide, or [³H]quinidine (0.5 mCi/ml) was added to each respective substrate concentration to allow quantitation of efflux from donor to receiver chambers. Lucifer yellow (100 µM) was added to the donor chamber to monitor integrity of the confluent cell monolayer.

Samples (25 µl) were taken over a period of 4 or 6 h, as shown, from both donor and receiver chambers into 96-well Lumaplates and dried overnight, after which the radioactivity was counted by using a TopCount model 9912 counter (PerkinElmer Life and Analytical Sciences, Boston, MA). The initial concentration measurement was taken at 6 min after the addition of substrate in the first well (Tran et al., 2005), and subsequent measurements were taken at multiple time points up to 6 h for all experiments (Tran et al., 2005). After each aliquot was taken, the transwell plates were placed in a shaker at a speed of 30 rpm at 37°C in 5% CO₂. The fluorescence of lucifer yellow (excitation maximum, 438 nm; emission maximum, 530 nm) was measured at time 0 from aliquots taken directly from the vials and compared with that of samples aliquoted at the end of the experiment from both the basolateral and apical chambers; the fluorescence was analyzed using a SpectraMax microplate reader. Passive permeability of lucifer yellow was always <10 nm/s over the entire experiment.

Cell Stability and Substrate Metabolism. The cells do not lose any capacity for drug transport over 6-h incubations (Tran et al., 2005; Acharya et al., 2006), implying that cellular ATP levels remain adequate for full P-gp function throughout the experiment. We showed previously that the stability of the cell monolayer and plasma membrane with respect to passive and active transport was not affected by the prolonged exposure times to amprenavir for at least 6 h (data not shown) (Tran et al., 2005). It was also shown that metabolism or decomposition was insignificant for amprenavir, quinidine, and loperamide on this time scale using radio-high-performance liquid chromatography (data not shown) (Tran et al., 2005).

Inhibition Studies. Cells were seeded and fed as described above. On the day of the experiment culture medium was removed by aspiration. The cells were then preincubated with the inhibitor/substrate drugs in both chambers for 30 min. During the preincubation, half of the wells received inhibitor solution without GF120918 (to study active transport) and the other half received inhibitor solution with GF120918 (to study passive transport). After 30 min, the preincubation solutions were removed by aspiration and the labeled substrate was added to the donor chamber, and both chambers were refilled with fresh inhibitor drug solution, again one half without GF120918 and the other half with it.

Digoxin Cell Concentration Exclusion Assay. Cells were seeded and fed as described above. On the day of the experiment, culture medium was removed by aspiration. The cells were then preincubated for 30 min with transport medium (TM), TM plus quinidine (20 µM), TM plus GF120918 (2 µM), or TM plus quinidine (20 µM) plus GF120918 (2 µM) in both chambers. The same protocol was repeated with an additional 1 µM benzbromarone. All of these preincubation conditions were repeated in triplicate. Preincubation solutions were then removed by aspiration. Donor and receiver solutions contained the radiolabeled substrate digoxin (5 µM), in addition to the respective preincubation condition. Basolateral and apical chambers were filled with donor and receiver solutions, respectively. Transport from B>A was measured over 4 h. Because the same concentration of digoxin was used on both sides, a flat curve showing steady state was obtained. At the 4th h, all solutions were removed. Both chambers were washed three times with cold transport media containing 2 µM GF120918. The polycarbonate membrane inserts with the cells were carefully cut out, and radioactivity associated to the cells on the membrane was counted by liquid scintillation.

Numerical Integrations. We used the stiffest integrator in MATLAB, ode23s, with absolute and relative tolerances set to 10^–8. Other MATLAB integrators, although faster, were not accurate enough at the later times of simulations. In data fitting, all concentration curves are simultaneously fitted, so that despite the fact that the A:B>A curve, i.e., substrate concentration in the apical chamber when the basolateral chamber is the donor, is the most visually striking, all curves contribute to minimizing the difference between data and simulated curves. MATLAB fminsearch minimizes the coefficient of variation between the data and the simulated curves. Further details can be found in Tran et al. (2005) and Bentz et al. (2005).

When the fitting included the novel basolateral and apical transporters, both V_max and K_m were fitted independently using eq. 2 in the mass action kinetic model for transport (Tran et al., 2005; Acharya et al., 2006). For the best fits, the V_max and K_m values did not converge, but their ratio was constant for all concentrations of loperamide (V_max/K_m = 100 s^–1) and digoxin (V_max/K_m = 30 s^–1) (Table 1). For loperamide, this value was the average of the best fits for a concentration range of 0.012 to 10 µM, i.e., 3 orders of magnitude, suggesting that a single transporter is responsible for the increased flux.

Statistical Analysis. For the digoxin and loperamide data, we sought statistical validation for the proposed basolateral (BT) and apical (AT) transporters. To compare the predictions by the kinetic model and actual measured amounts of digoxin in apical and basal compartments over time, we sought to evaluate four hypotheses: that predictions from adequate models would 1) be highly correlated with empirical measurements, 2) have small root-mean-square (RMS) prediction errors, 3) be linearly related to empirical measurements, and 4) have a slope of 1 for the regression of empirical data on the predictions. Regressions were performed using singular value decomposition (Press et al., 1992). Refinement of predictions by allowing curvilinear regressions was measured by the contribution of second-order terms to minimize prediction errors (Rao, 1998; Quinn and Keough, 2002). Regression slopes were compared with unity using t tests comparing the estimated slopes to the hypothetical values (Sokal and Rohlf, 1995; Rao, 1998).

Because the amounts of digoxin remaining in the basal compartment or transported to the apical compartment were measured, sequential measurements were potentially autocorrelated and hence, not independent, violating an assumption of regression analyses. To circumvent this violation, we performed autocorrelation analyses and showed that with lag times of 2 to 3 h, autocorrelation (r²) fell to <0.01 with p > 0.70. Hence, values were independent after this time lag. All analyses reported here used sequential time points every 3 h.

Loperamide trials were shorter than the stitched digoxin trials (6 versus 30 h), but transport experiments were run for several initial concentrations of loperamide. Assuming that parameters remained constant at all loperamide concentrations allowed us to construct an analysis analogous to a multiple linear regression analysis, treating each transporter (P-gp-based transport, basolateral transport by a BT, and apical transport by an AT) as a new predictor using 1 df. The need to add further transporters to a minimal model could be tested analogously to the need to add new predictors in a multiple linear regression model, by testing the improvement of the prediction via the decrease in the residual error of the prediction versus measurement regression (Rao, 1998; Quinn and Keough, 2002). Because adding a new transporter yielded an entirely new simulation and an entirely new prediction, this approach technically does violate one of the assumptions of the linear analysis, i.e., independence of the transporters. Therefore, the probability estimates must be regarded as approximate.

Mass Action Kinetic Model. Fig. 1 is a cartoon model of a confluent cell monolayer, featuring the polarized MDCKII-hMDR1 cells, where the basolateral membrane is attached to the polycarbonate filters and P-gp (upward arrows) expressed on the apical surface. The apical and basolateral chambers are kept separate by the tight junctions. Active transport by P-gp occurs vectorially, with substrate binding to a site on P-gp within the apical membrane inner monolayer and with efflux into the apical chamber (Loo and Clarke, 2005; Lugo and Sharom, 2005).

View larger version (37K): [in this window] [in a new window]   FIG. 1. Cartoon model of an MDCKII-hMDR1 confluent cell monolayer. The apical membrane with microvilli is on top and the basolateral membrane is on the bottom, where the cells attach to the polycarbonate insert. Passive permeability occurs in both directions and is measured when P-gp (open arrows) is completely inhibited by GF120918. P-gp expressed on the apical membrane transports substrate from the inner apical monolayer into the apical chamber. The concentrations of substrate in the apical and basolateral chambers, CA and CB, are measured. The concentrations of substrate in the inner plasma membrane, CPC, and the cytosol, CC, are predicted as part of the data-fitting analysis. Transporters other than P-gp are known to exist in the basolateral and apical membranes, shown by the closed double arrows, which is not meant to suggest that they are all bidirectional.

The simplest Michaelis-Menten mass action reaction to model P-gp is

(1)

where T₀ is the empty transporter, C_PC is the substrate in the apical membrane inner monolayer, T_C is the transporter bound by substrate, and C_A is the substrate after efflux into the apical chamber. Although P-gp has more than one efflux pathway with cooperative interactions at low probe substrate and inhibitor substrate concentrations, when the inhibitor concentrations are near the IC₅₀ of the inhibitor, P-gp behaves like a simple single competitive site transporter (Acharya et al., 2006). The additional basolateral and apical transporters are modeled by Michaelis-Menten steady-state equations

(2)

This is a compromise between rigor and flexibility, because this model can assess roughly whether there is a kinetically significant binding step to the transporter before influx into the cell, without having to fit real binding constants. For each transporter, both V_max and K_m are used to fit the data, but we have found, thus far, that only their ratio, V_max/K_m (seconds^–1), was constant. From eq. 2, this would imply that the K_m of the transporter is much larger than the substrate concentration in the cytosol. For this modeling, we have assumed that the unknown transporters are facilitated or bidirectional transporters, i.e., they support transport in both directions. This was not a strong assumption, because substrate flux was mostly unidirectional import during these experiments. Once the transporters are identified, more rigorous modeling can be applied.

For many substrates, including those studied here, passive permeability is a significant fraction of total transport and is quantitatively analyzed separately using the P-gp inhibitor, GF120918 (Hyafil et al., 1993; Tran et al., 2005; Acharya et al., 2006). Our operational definition of passive permeability is +GF120918, which completely inhibits P-gp and some other transporters, both known and novel. For transporters that are not inhibited by GF120918, their transport would be accounted for as part of the passive permeability component, regardless of their actual mechanism of transport. The other part of the passive permeability is the passive diffusion through the lipid bilayers of both membranes and through the tight junctions. Whereas the relative contribution for each "passive" pathway is not known, that information is not essential for our analysis of P-gp function. What is essential is that the passive permeability is accounted for as accurately as possible, so that we can accurately measure how much of the total transport, –GF120918, is due to P-gp and any GF120918 sensitive transporters (Tran et al., 2005).

We have fitted the drug-independent elementary kinetic rate constants of P-gp, i.e., the density of efflux of active P-gp in the apical membrane and the drug-P-gp association rate constant, using drugs that have no other transporters in this cell line, i.e., amprenavir and quinidine (Tran et al., 2005; Acharya et al., 2006). Thus, it is possible now to fit the kinetic parameters of these novel transporters for digoxin and loperamide independently, as shown below.

	Results

Top Abstract Materials and Methods Results Discussion References

 
The kinetic parameters for P-gp transport published in Acharya et al. (2006) were determined under the assumption that P-gp was the sole transporter. This model worked very well for amprenavir and quinidine. However, the "P-gp alone" model gave adequate fits for loperamide only with concentrations >3 µM (Tran et al., 2005; Acharya et al., 2006). At lower concentrations, more loperamide efflux was observed than could be fitted assuming that the passive permeability coefficients fitted +GF120918 and that the P-gp kinetic parameters fitted at high (3 µM) loperamide concentrations.

The loperamide data have been fitted again, now allowing a basolateral transporter, modeled by the Michaelis-Menten equation, eq. 2. Both B>A and A>B transport data were used to find the consensus value of V_max/K_m. Here, for economy, we will only show the data and fits for the B>A transport of loperamide across MDCKII-hMDR1 monolayers. The A>B transport data fitted just as well.

Figure 2 shows the concentration of loperamide over time in each chamber and the fits for three of the eight loperamide concentrations we have examined, from 12 nM to 30 µM. The B:B>A and A:B>A curves denote the concentrations of drug in the basolateral and apical chambers when the basolateral chamber was the donor, i.e., transport runs B>A. Data are shown by when the receiver apical chamber is sampled and by when the donor basolateral chamber is sampled. All data points represent triplicate measurements at each time point with the corresponding S.D. With the best fits, the V_max and K_m values did not converge, but their ratio was constant for all concentrations of loperamide <µM, where P-gp saturated and contributes little to total transport. V_max/K_m is 100 ± 20 s^–1 (Table 1) was the average of the best fits for loperamide concentrations from 0.012 to 10 µM, i.e., 3 orders of magnitude, suggesting that a single transporter is responsible for the increased basolateral flux. The fits are shown by a dotted line for P-gp alone (Acharya et al., 2006) and by a solid line when the basolateral facilitated transporter, denoted BT, is added using eq. 2. In addition, all curves contain the predicted passive efflux, +GF120918.

View larger version (18K): [in this window] [in a new window]   FIG. 2. Loperamide B>A transport data in MDCKII-MDR1 cells fitted with and without a basolateral transporter. The concentration of loperamide over time across the MDCKII-MDR1 cell monolayer in each chamber is shown as a function of different initial loperamide concentrations in the basolateral chamber. The B:B>A and A:B>A curves denote the concentration of drug in the basolateral and apical chambers, respectively, when the basolateral chamber was the donor, i.e., transport runs B>A. , concentration in the apical chamber; , concentration in the basolateral chamber. All data points represent triplicate measurements at each time point with the S.D. shown by the error bars. The fits for P-gp alone are shown by a dotted line (Acharya et al., 2006), whereas the solid line shows the effect of adding a basolateral facilitated transporter, with Vmax/Km = 100 s–1. All curves also contain the passive permeability (+GP120918) contribution. A, starting with 0.03 µM loperamide shows a clear underestimate of the data for P-gp alone. B, starting with 1 µM loperamide shows that the two models are converging, but the curve for P-gp and the basolateral transporter is clearly superior. C, starting with 10 µM loperamide shows that the two models have largely converged, as the P-gp is nearly saturated and most transport is through passive permeability. Yet, the model with the basolateral transporter remains superior: A:B>A after 3 h.

Although the fits looked good, we wanted to quantitate their statistical significance. We compared predictions separately for apical and basolateral compartments. In each case, predictions were significantly improved by adding the BT to the model with only the P-gp transporter (p < 5 x 10^–7). For neither the apical nor the basolateral compartment did adding an AT to the P-gp plus BT model further improve the predicted concentrations of loperamide (p > 0.7).

We also fitted the transport data of loperamide in the A>B direction. The fits were just as good (data not shown). This finding suggests that the basolateral transporter is bidirectional, as it was modeled.

Figure 3 shows the B>A transport of 1 µM digoxin across the MDCKII-MDR1 confluent cell monolayers. Figure 3A shows that for the first 6 h the concentrations change roughly linearly with time, yielding no "fittable" data. The approach to the true steady-state is required to fit the kinetic parameters. When the cells were incubated for an extended period of time (>12 h), the transport curves showed toxic effects, e.g., lucifer yellow leakage or an A:B>A curve first increasing and then decreasing (data not shown).

View larger version (20K): [in this window] [in a new window]   FIG. 3. Digoxin B>A transport curve is stitched together by measuring transport across the MDCKII-MDR1 cell monolayers over consecutive 6-h experiments. A, concentration of digoxin over time in each chamber starting with 1 µM initially in the basolateral donor chamber. The increase in digoxin concentration in the apical chamber is linear, so fitting was not possible. B, to measure transport for a longer time course, we constructed a "data stitching" method. Once we know the concentrations of digoxin in the apical and basolateral chambers at the 6th h of A, we start the next experiment such that the initial concentration of digoxin in both chambers match as closely as possible. The concentration of digoxin in the apical chamber increases over time and nearly reaches a steady state. C, final stitched data. The fits for P-gp alone are shown by a dashed line (Acharya et al., 2006), whereas the solid line shows the effect of adding a basolateral facilitated transporter, with Vmax/Km = 30 s–1. All curves also contain the passive permeability (+GP120918) contribution. P-gp alone clearly underestimates transport by at least 50%. The solid line shows the fit from adding a basolateral transporter with Vmax/Km = 30 s–1, in addition to P-gp. The solid line fits up to 12 h and beyond that, it predicts more transport than we actually measure. An apical facilitated transporter or active importer for digoxin, in addition to P-gp, would cause the final steady-state concentration for digoxin in the apical chamber to be smaller. This is shown by the dashed line, with Vmax/Km = 2 s–1 for the additional apical facilitator transporter. Statistical analysis of these data shows that the basolateral transporter is significant with p < 0.01 and the apical transporter is significant with p < 0.05, as explained in the text.

To extend the time of the experiment we constructed a data stitching experiment. Figure 3B illustrates this construction. Once the concentrations of digoxin in the apical and basolateral chambers at the 6th h are known, the next experiment begins with these initial concentrations of digoxin in the appropriate chambers. Using this approach, data are collected for consecutive stretches of 6 h and stitched together to create a time course up to 30 h. The only artifact in stitching is that for the next experiment the cytosol is initially empty of substrate, which it would not be at the end of the prior experiment. However, because the volume of the confluent cell monolayer cytosol is very small, the cells will fill within the first few minutes of substrate addition, which is insignificant on the time scale of these experiments (Acharya et al., 2006).

The kinetic fits for the digoxin concentration time curve are shown in Figure 3C. The fits are shown by a dotted line for P-gp alone (Acharya et al., 2006), a solid line when the basolateral facilitated transporter, denoted BT, is added and a dashed line when the basolateral facilitated transporter and the apical facilitated transporter, denoted AT, are added using eq. 2. In addition, all curves contain the predicted passive efflux, +GF120918. The P-gp alone model underestimates transport by at least 50% at all times. Adding a basolateral transporter yields a good fit up to 12 h when V_max/K_m is 30 s^–1, and the parameters are given in Table 1. As mentioned above V_max and K_m were fitted independently, but only their ratio was constant. However, beyond 12 h it predicts too much transport, well away from the S.D. of the data. Addition of an apical facilitated transporter with a V_max/K_m = 2s^–1 for digoxin, in addition to P-gp and the basolateral transporter, yields a lower true steady-state concentration for digoxin in the apical chamber and a good fit to the data.

To estimate the statistical significance of these fits, we compared model predictions to measurements separately for apical and basolateral compartments. In all cases, model predictions were highly correlated with measured amounts of digoxin in each compartment (r² > 0.90 for linear regressions). Correlations increased and RMS prediction errors decreased, in general, when a BT and then an AT were added to the models (r² > 0.995, RMS error < 0.03 for full models), Adding an apical transporter to the P-gp plus BT model neither increased the correlation nor decreased the prediction error for the loss of digoxin from the basolateral chamber.

Regressing measured basolateral compartment digoxin data versus model predictions yielded slopes significantly >1.0 for the P-gp alone model (slope = 1.86, p < 0.01). Both the P-gp plus BT model and the P-gp plus BT plus AT model yielded virtually indistinguishable prediction with slopes (0.974 and 0.984, respectively) not different from 1.0 (p > 0.05). Regressions of apical compartment digoxin on predictions gave slopes of 1.58, 0.822, and 0.966 for the P-gp alone model, P-gp plus BT model, and P-gp plus BT plus AT model, respectively. Only the P-gp plus BT plus AT model slope was indistinguishable from 1.0 at p = 0.05.

Analyses of the curvilinearity of relationships gave similar results. For basal compartment digoxin, the P-gp alone model regression was strongly curved (p < 0.01), but the other two models had no significant curvature (p > 0.05) to the regressions. For apical compartment digoxin, both the P-gp and P-gp plus BT model regressions were strongly curved (p < 0.0001). The P-gp plus BT plus AT model regression was significantly curved (p = 0.01) but with a curvature so small that it did not affect the biological predictions.

An attempt was made to identify the canine kidney transporters involved with digoxin and loperamide transport. By preincubating cells with a "high" concentration of a substrate/inhibitor for the basolateral transporter, the transport of digoxin should decrease to the levels predicted by the P-gp efflux alone curves (e.g., a 50% reduction in digoxin transport). This makes the kinetic analysis a very sensitive instrument for determining substrates or inhibitors of the basolateral transporter. Discovering a substrate/inhibitor of the apical digoxin transporter would be more difficult, given its smaller effect on total transport.

Inhibition studies were completed with a variety of prototypical substrates or inhibitors for OAT, OATP, OCT, and MRP transporters (Table 2). The compounds in Table 2 were selected because they cover a broad range of possibilities for known transporters on the basolateral or apical membranes of human and mouse kidney cells. To evaluate whether these compounds have any effect on loperamide or digoxin transport, we used concentrations at or higher than the maximum concentration used for the experiments in the respective references. Clearly the human isoform of a particular transporter could have a substrate or inhibitor not shared by its canine isoform. So this strategy need not succeed, but it is certainly the obvious first step.

Except for loperamide, incubations with the substrates/inhibitors failed to inhibit digoxin transport, B>AorA>B, within 3 h; the canine transporters do not interact with these substrates/inhibitors. The inhibition of digoxin transport by loperamide and vice versa certainly involves P-gp, but to determine whether they share a single basolateral transporter we will need a protocol for assaying digoxin transport kinetic parameters more practical than the cumbersome "stitching" approach used here.

We did find that 100 and 200 µM benzbromarone inhibited 1 µM digoxin transport in the B>A direction by 20 and 50%, respectively (data not shown), but these concentrations also showed toxic effects on the cells, e.g., the transport curves stopped abruptly after a few hours and the cell monolayer became leaky. These results led us to perform the cell concentration exclusion assay using 1 µM benzbromarone with the aim to determine whether it actually caused less accumulation of digoxin in the cells. The amount of radioactivity from digoxin associated with the cells was measured in the presence and absence of benzbromarone. Presumably if benzbromarone inhibits the basolateral and/or apical transporter(s), then less digoxin would be accumulated within the cells.

Figure 4 shows cellular digoxin radioactivity by itself and in the presence of quinidine (a P-gp substrate inhibitor), GF120918 (inhibits P-gp and the additional basolateral and apical transporters), and benzbromarone. In the presence of 20 µM quinidine, P-gp is relatively saturated (Acharya et al., 2006), and there is more digoxin present within the cells compared with digoxin alone. In the presence of 2 µM GF120918, P-gp and the additional transporters are completely inactive, and there is even more digoxin present within the cells as a result of the lack of P-gp efflux of digoxin. In contrast, in presence of 1 µM benzbromarone there is significantly less digoxin accumulation, suggesting inhibition of a basolateral and/or apical uptake transporter(s), whereas P-gp could actively efflux digoxin.

View larger version (13K): [in this window] [in a new window]   FIG. 4. Digoxin cell concentration exclusion by benzbromarone. This figure shows radioactivity in disintegrations per minute from radiolabeled digoxin (DGX) in MDCKII-hMDR1 cells by itself, in the presence of 20 µM quinidine (QND), 2 µM GF120918 (918), and 1 µM benzbromarone (BZD).

	Discussion

Top Abstract Materials and Methods Results Discussion References

 
The rigorous kinetic analysis of P-gp-mediated transport through the confluent MDCKII-hMDR1 cell monolayer has provided a quantitative understanding of P-gp function, because the elementary rate constants so clearly define the function of P-gp and the selective pressure on it (Tran et al., 2005; Acharya et al., 2006). The transports of amprenavir and quinidine across the MDCKII-hMDR1 confluent cell monolayer are quantitatively fitted by P-gp and passive permeability alone. However, for loperamide and digoxin the kinetic modeling shows other novel transporter(s) in the MDCKII cells that facilitate passage through the basolateral membrane and also through the apical membrane for digoxin.

In Tran et al. (2005) and Acharya et al. (2006) it was noted that loperamide transport, especially at low concentrations, was being affected by another process. Adding a facilitated transporter on the basolateral membrane allows loperamide to enter the cytosol from the basolateral chamber more rapidly than predicted by the +GF120918 passive permeability coefficients, i.e., the loperamide basolateral transporter is inhibited by GF120918. We obtained good fits from 12 nM to 30 µM loperamide with the transporter over a 3 order of magnitude range in concentration (Table 1). Statistical analysis supported the existence of the basolateral transporter (p < 0.000001). There is no statistical evidence for an apical loperamide transporter. Because we fitted both B>A and A>B transport data for loperamide using the basolateral transporter, this result suggests that transport is bidirectional.

This addition to the kinetic model for loperamide has changed the estimate for the binding constant of loperamide to P-gp from the apical membrane, K_C, from 4000 M^–1 (Acharya et al., 2006) to 20,000 M^–1 (Table 1). The binding constant is fit primarily from A>B transport (Tran et al., 2005), and the basolateral transporter for loperamide speeds up permeation to the basolateral chamber. Thus, a stronger binding constant to P-gp is required to fit the same A>B data.

For digoxin, the passive permeability was so slow that transport was linear with time for the first 6 h. Without some degree of curvature, i.e., approaching steady state, fitting cannot be done. We applied a stitching method to approximate a 30-h transport experiment. The end of one 6-h experiment became the starting point for the next 6-h experiment. Both a basolateral transporter and an apical transporter were needed to fit this data. We did not attempt the A>B experiment, because the transport in this direction would be much smaller, thereby making the stitching approach more problematic. We are currently working on a new protocol under which steady state should be reached faster. This will be required to determine whether digoxin and loperamide use the same basolateral transporter or not. They partially inhibit each other's transport, but we cannot yet prove whether or not this inhibition is beyond what would be predicted by inhibition of P-gp alone.

For digoxin, there is a basolateral importer and an apical importer that allow its reentry into the cells after efflux into the apical chamber. Both were modeled as bidirectional, but only B>A data were fitted, so only import was kinetically important. Statistical analyses have shown that the basolateral transporter was required to precisely predict the amount of digoxin remaining in the basolateral compartment (p < 0.01). Both the basolateral and apical transporters were needed to accurately predict the amount of digoxin transported to the apical compartment (p < 0.05).

The magnitude of the increased flux of digoxin and loperamide into the cells (–GF120918) proves that the basolateral transporter(s) are inhibited by GF120918. For the transport of loperamide from the basolateral chamber into the cytosol, 60% is due to passive permeability (+GF120918) and 40% is due to the transporter (–GF120918) (Table 1). For the transport of digoxin from the basolateral chamber into the cytosol, 30% is due to passive permeability (+GF120918) and 70% is due to the basolateral transporter (–GF120918) (Table 1).

On the other hand, for the digoxin transporter from the apical chamber into the cytosol, 95% of the flux is due to passive permeability (+GF120918) and 5% is due to the apical transporter (–GF120918). Five percent digoxin transport by the apical transporter is too small to determine whether GF120918 inhibits this transporter, which will require performing the A>B transport experiment with an improved experimental protocol.

As a compromise between computational cost and rigor, we modeled the transporters by the Michaelis-Menten equation, with both V_max and K_m independently fitted. For all loperamide concentrations, the best fits for the basolateral transporter yielded 1 to 2 order of magnitude ranges for K_m and V_max, but their ratio was constant for each fit. This result means that the fitting predicted a first-order rate constant for this transporter, just like a passive permeability. So there was no evidence for saturation of the basolateral transporter up to 10 µM loperamide. The convergence between the P-gp alone curve and the P-gp plus basolateral transporter curve is due to the saturation of P-gp at a high substrate concentration.

We have tried to identify these transporters by adding prototypical substrates/inhibitors of known transporters. If we found one, then the loperamide or digoxin transport curves would drop to the levels predicted by the "P-gp only" curves (Figs. 2 and 3). The change in the predicted curves would be large, making the kinetic analysis an extremely sensitive assay for discovering functional transporters in the confluent cell monolayer. Although this approach is reasonable, it is always possible that the human or mouse transporter would not share a substrate/inhibitor with its canine counterpart.

None of the compounds tested in Table 2 showed simple inhibition of transport of digoxin or loperamide during the first 3 h, except loperamide and digoxin, respectively. There were some changes noted after 3 h, but these appeared to be related to toxicity. Because the other transporter(s) face the basolateral or apical chambers, there is no reason that simple inhibition would require more than 3 h (Acharya et al., 2006).

There may be several reasons for the lack of an inhibitory effect by the P-gp substrate fexofenadine. Fexofenadine has a very low passive permeability through Caco-2 cells (Petri et al., 2004). We have shown that at "low" concentrations of probe substrate and of inhibitor substrate, there is no inhibition of probe transport because P-gp has more than one efflux pathway (Acharya et al., 2006). Also, it may be that the binding constant of fexofenadine is much smaller than that of digoxin, so a lower digoxin concentration may be needed.

A cell concentration exclusion assay was also tried, wherein the concentration of digoxin within the confluent cell monolayer was assayed as a function of benzbromarone concentration. Addition of benzbromarone resulted in a very significant reduction of digoxin accumulation within the cell monolayer, as expected for a substrate/inhibitor of the basolateral digoxin transporter. However, benzbromarone had no significant effect on digoxin B>AorA>B transport.

The reduction in cell-associated digoxin by benzbromarone is due to another mechanism besides inhibition of the basolateral or apical transporter. The simplest hypothesis is that benzbromarone reduces the cytosolic volume of the cell, so less digoxin can accumulate. This reduction would cause no predicted change in the transport of digoxin, because the kinetic model predicts that nearly all of it diffuses along the inner plasma membrane (Tran et al., 2005). Thus, cell concentration exclusion assays must be carefully controlled when they are used with amphipathic compounds effluxed from a transporter whose binding sites are within the plasma membrane, such as P-gp (Lugo and Sharom, 2005).

In conclusion, we have used a rigorous kinetic analysis to detect the presence of novel transporters of loperamide and digoxin in the basolateral membrane of this cell line, as well as a digoxin transporter in the apical membrane. Based on the current findings, one could question whether digoxin and loperamide are robust substrates for screening compounds for P-gp interactions in the MDCKII-hMDR cell line. The worst case scenario in using them would be a false-positive result if the novel compound did indeed inhibit the basolateral transporter but not P-gp. If a compound inhibits the efflux of digoxin, a follow-up experiment with amprenavir or quinidine would identify the source of the inhibition. Quinidine and amprenavir appear to be robust probe substrates, as their transport across the MDCKII-hMDR1 cell monolayer depends only upon passive permeability and P-gp. Use of another cell line, e.g., LLC-hMDRI or Caco-2, with digoxin or loperamide does not avoid the problem, because the same kinetic analysis as shown here is required to determine which compounds are robust P-gp probe substrates for individual cell lines.

	Footnotes

 
Article, publication date, and citation information can be found at http://dmd.aspetjournals.org.

doi:10.1124/dmd.107.017301.

ABBREVIATIONS: P-gp, the P-glycoprotein product of the hMDR1 gene; MDCKII-hMDR1, Madin-Darby canine kidney II cell line that overexpresses human multidrug resistance 1; bcrp/BRCP, breast cancer resistance protein; 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; DMEM, Dulbecco's modified Eagle's medium; A>B (or B> transport across the confluent cell monolayer when the donor chamber is apical (or basolateral) and the receiver chamber is basolateral (or apical); TM, transport medium; BT, basolateral transporter; AT, apical transporter; RMS, root-mean-square; OAT, organic anion transporters; OATP, organic ion-transporting polypeptide; OCT, organic cation transporters; MRP, multidrug resistance-associated protein.

Address correspondence to: Dr. Joe Bentz, Department of Bioscience & Biotechnology, Drexel University, 32nd and Chestnut Sts., Philadelphia, PA 19104. E-mail: bentzj{at}drexel.edu

	References

Top Abstract Materials and Methods Results Discussion References

 


Acharya P, Tran TT, Polli JW, Ayrton A, Ellens H, and Bentz J (2006) P-gp expressed in a confluent monolayer of hMDR1-MDCKII cells has more than one efflux pathway with cooperative binding sites. Biochemistry 45: 15505–15519.[CrossRef][Medline]

Bartholomé K, Rius M, Letschert K, Keller D, Timmer J, and Keppler D (2007) Data-based mathematical modeling of vectorial transport across double-transfected polarized cells. Drug Metab Dispos 35: 1476–1481.[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 K_m. Pharmacol Res 22: 1667–1677.[CrossRef]

Collett A, Tanianis-Hughes J, Carlson GL, Harwood MD, and Warhurst G (2005) Comparison of P-glycoprotein-mediated drug-digoxin interactions in Caco-2 with human and rodent intestine: relevance to in vivo prediction. Eur J Pharm Sci 26: 386–393.[CrossRef][Medline]

Dean M, Rzhetsky A, and Allikmets R (2001) The human ATP-binding cassette (ABC) transporter superfamily. Genome Res 11: 1156–1166.[Abstract/Free Full Text]

Dresser MJ, Leabman MK, and Giacomini KM (2001) Transporters involved in the elimination of drugs in the kidney: organic anion transporters and organic cation transporters. J Pharm Sci 90: 397–421.[CrossRef][Medline]

Enomoto A and Endou H (2005) Roles of organic anion transporters (OATs) and a urate transporter (URAT1) in the pathophysiology of human disease. Clin Exp Nephrol 9: 195–205.[CrossRef][Medline]

Harris MJ, Kagawa T, Dawson PA, and Arias IM (2004) Taurocholate transport by hepatic and intestinal bile acid transporters is independent of FIC1 overexpression in Madin-Darby canine kidney cells. J Gastroenterol Hepatol 19: 819–825.[CrossRef][Medline]

Horikawa M, Kato Y, Tyson CA, and Sugiyama Y (2002) The potential for an interaction between MRP2 (ABCC2) and various therapeutic agents: probenecid as a candidate inhibitor of the biliary excretion of irinotecan metabolites. Drug Metab Pharmacokinet 17: 23–33.[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 acridone carboxamide derivative. Cancer Res 53: 4595–4602.[Abstract/Free Full Text]

Ismair MG, Stanca C, Ha HR, Renner EL, Meier PJ, and Kullak-Ublick GA (2003) Interactions of glycyrrhizin with organic anion transporting polypeptides of rat and human liver. Hepatol Res 26: 343–347.[CrossRef][Medline]

Iwanaga T, Kobayashi D, Hirayama M, Maeda T, and Tamai I (2005) Involvement of uric acid transporter in increased renal clearance of the xanthine oxidase inhibitor oxypurinol induced by a uricosuric agent, benzbromarone. Drug Metab Dispos 33: 1791–1795.[Abstract/Free Full Text]

Juliano RL and Ling V (1976) A surface glycoprotein modulating drug permeability in Chinese hamster ovary cell mutants. Biochim Biophys Acta 455: 152–162.[Medline]

Kikuchi R, Kusuhara H, Abe T, Endou H, and Sugiyama Y (2004) Involvement of multiple transporters in the efflux of 3-hydroxy-3-methylglutaryl-CoA reductase inhibitors across the blood-brain barrier. J Pharmacol Exp Ther 311: 1147–1153.[Abstract/Free Full Text]

Korjamo T, Kemiläinen H, Heikkinen AT, and Mönkkönen J (2007) Decrease in intracellular concentration causes the shift in K_m value of efflux pump substrates. Drug Metab Dispos 35: 1574–1579.[Abstract/Free Full Text]

Lau YY, Okochi H, Huang Y, and Benet LZ (2006) Multiple transporters affect the disposition of atorvastatin and its two active hydroxy metabolites: application of in vitro and ex situ systems. J Pharmacol Exp Ther 316: 762–771.[Abstract/Free Full Text]

Loo TW and Clarke DM (2005) Recent progress in understanding the mechanism of P-glycoprotein-mediated drug efflux. J Membr Biol 206: 173–185.[CrossRef][Medline]

Lugo MR and Sharom FJ (2005) Interaction of LDS-751 with P-glycoprotein and mapping of the location of the R drug binding site. Biochemistry 44: 643–655.[CrossRef][Medline]

Mita S, Suzuki H, Akita H, Hayashi H, Onuki R, Hofmann AF, and Sugiyama Y (2006) Inhibition of bile acid transport across Na⁺/taurocholate cotransporting polypeptide (SLC10A1) and bile salt export pump (ABCB 11)-coexpressing LLC-PK1 cells by cholestasis-inducing drugs. Drug Metab Dispos 34: 1575–1581.[Abstract/Free Full Text]

Mizuno N, Niwa T, Yotsumoto Y, and Sugiyama Y (2003) Impact of drug transporter studies on drug discovery and development. Pharmacol Rev 55: 425–461.[Abstract/Free Full Text]

Morita N, Kusuhara H, Sekine T, Endou H, and Sugiyama Y (2001) Functional characterization of rat organic anion transporter 2 in LLC-PK1 Cells. J Pharmacol Exp Ther 298: 1179–1184.[Abstract/Free Full Text]

Morrow CS, Peklak-Scott C, Bishwokarma B, Kute TE, Smitherman PK, and Townsend AJ (2006) Multidrug resistance protein 1 (MRP1, ABCC1) mediates resistance to mitoxantrone via glutathione-dependent drug efflux. Mol Pharmacol 69: 1499–1505.[Abstract/Free Full Text]

Pan G and Elmquist WF (2007) Mitoxantrone permeability in MDCKII cells is influenced by active influx transport. Mol Pharmacol 4: 475–483.[CrossRef][Medline]

Petri N, Tannergren C, Rungstad D, and Lennernas H (2004) Transport characteristics of fexofenadine in the Caco-2 cell model. Pharmacol Res 21: 1398–1404.[CrossRef]

Press WH, Teukolsky SA, Vettering WT, and Flannery BP (1992) Numerical Recipes in C: The Art of Scientific Computing, 2nd ed, Cambridge University Press, Cambridge, UK.

Quinn GP and Keough MJ (2002) Experimental Design and Data Analysis for Biologists, Cambridge University Press, Cambridge, UK.

Rao PV (1998) Statistical Research Methods in the Life Sciences, Duxbury Press, Pacific Grove, CA.

Rautio J, Humphreys JE, Webster LO, Balakrishnan A, Keogh JP, Kunta JR, Serabjit-Singh CJ, and Polli JW (2006) In vitro P-glycoprotein inhibition assays for assessment of clinical drug interaction potential of new drug candidates: a recommendation for probe substrates. Drug Metab Dispos 34: 786–792.[Abstract/Free Full Text]

Robertson EE and Rankin GO (2006) Human renal organic anion transporters: characteristics and contributions to drug and drug metabolite excretion. Pharmacol Ther 109: 399–412.[CrossRef][Medline]

Sekine T, Miyazaki H, and Endou H (2006) Molecular physiology of renal organic anion transporters. Am J Physiol Renal Physiol 290: F251–F261.[Abstract/Free Full Text]

Shitara Y, Horiea T, and Sugiyama Y (2006) Transporters as a determinant of drug clearance and tissue distribution. Eur J Pharm Sci 27: 425–446.[CrossRef][Medline]

Sokal RR and Rohlf FJ (1995) Biometry, 3rd ed, W. H. Freeman and Company, New York.

Spears KJ, Ross J, Stenhouse A, Ward CJ, Goh LB, Wolf CR, Morgan P, Ayrton A, and Friedberg TH (2005) Directional trans-epithelial transport of organic anions in porcine LLC-PK1 cells that co-express human OATP1B1 (OATP-C) and MRP2. Biochem Pharmacol 69: 415–423.[CrossRef][Medline]

Sugiyama D, Kusuhara H, Shitara Y, Abe T, Meier PJ, and Sekine T (2001) Characterization of the efflux transport of 17β-estradiold-17β-glucuronide from the brain across the blood-brain barrier. J Pharmacol Exp Ther 298: 316–322.[Abstract/Free Full Text]

Tahara H, Kusuhara H, Maeda K, Koepsell H, Fuse E, and Sugiyama Y (2006) Inhibition of Oat3-mediated renal uptake as a mechanism for drug-drug interaction between fexofenadine and probenecid. Drug Metab Dispos 34: 743–747.[Abstract/Free Full Text]

Tang F, Horie K, and Borchardt RT (2002) Are MDCK cells transfected with the human MDR1 gene a good model of the human intestinal mucosa? Pharmacol Res 19: 773–779.[CrossRef]

Tran TT, Mittal A, Aldinger T, Polli JW, Ayrton A, Ellens H, and Bentz J (2005) The elementary mass action rate constants of P-gp transport for a confluent monolayer of MDCKII-hMDR1 cells. Biophys J 88: 715–738.[CrossRef][Medline]

Troutman MD and Thakker DR (2003) novel experimental parameters to quantify the modulation of absorptive and secretory transport of substrates by P-glycoprotein in cell culture models of intestinal epithelium. Pharmacol Res 220: 1210–1224.

Vavricka SR, Montfoort JV, Ha HR, Meier PJ, and Fattinger K (2002) Interactions of rifamycin SV and rifampicin with organic anion uptake systems of human liver. Hepatology 36: 164–172.[CrossRef][Medline]

This article has been cited by other articles:

				E. Kamiyama, D. Sugiyama, D. Nakai, S.-i. Miura, and O. Okazaki Culture Period-Dependent Change of Function and Expression of ATP-Binding Cassette Transporters in Caco-2 Cells Drug Metab. Dispos., September 1, 2009; 37(9): 1956 - 1962. [Abstract] [Full Text] [PDF]

				L. Wang, M. Leggas, M. Goswami, P. E. Empey, and P. J. McNamara 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 (GF120918) As a Chemical ATP-Binding Cassette Transporter Family G Member 2 (Abcg2) Knockout Model to Study Nitrofurantoin Transfer into Milk Drug Metab. Dispos., December 1, 2008; 36(12): 2591 - 2596. [Abstract] [Full Text] [PDF]

This Article Abstract Full Text (PDF) All Versions of this Article: dmd.107.017301v1 36/2/452    most recent Submit a response Alert me when this article is cited Alert me when eLetters are posted Alert me if a correction is posted Services Similar articles in this journal Similar articles in PubMed Alert me to new issues of the journal Download to citation manager Citing Articles Citing Articles via HighWire Citing Articles via Google Scholar Google Scholar Articles by Acharya, P. Articles by Bentz, J. Search for Related Content PubMed PubMed Citation Articles by Acharya, P. Articles by Bentz, J.