Variability in P-glycoprotein inhibitory potency (IC₅₀) using various in vitro experimental systems: implications for universal digoxin drug-drug interaction risk assessment decision criteria

Prediction of drug-ABC-transporter interaction--Recent advances and future challenges

With the discovery of P-glycoprotein (P-gp), it became evident that ABC-transporters play a vital role in bioavailability and toxicity of drugs. They prevent intracellular accumulation of toxic compounds, which renders them a major defense mechanism against xenotoxic compounds. Their expression in cells of all major barriers (intestine, blood–brain barrier, blood–placenta barrier) as well as in metabolic organs (liver, kidney) also explains their influence on the ADMET properties of drugs and drug candidates. Thus, in silico models for the prediction of the probability of a compound to interact with P-gp or analogous transporters are of high value in the early phase of the drug discovery process. Within this review, we highlight recent developments in the area, with a special focus on the molecular basis of drug–transporter interaction. In addition, with the recent availability of X-ray structures of several ABC-transporters, also structure-based design methods have been applied and will be addressed.

Nature and uses of fluorescent dyes for drug transporter studies

INTRODUCTION

Drug transporters are now recognized as major players involved in pharmacokinetics and toxicology. Methods for assessing their activity are important to consider, particularly owing to regulatory requirements with respect to inhibition of drug transporter activity and prediction of drug-drug interactions. In this context, the use of fluorescent-dye-based transport assays is likely to deserve attention.

AREAS COVERED

This review provides an overview of the nature of fluorescent dye substrates for ATP-binding cassette and solute carrier drug transporters. Their use for investigating drug transporter activity in cultured cells and clinical hematological samples, drug transporter inhibition, drug transporter imaging and drug transport at the organ level are summarized.

EXPERT OPINION

A wide range of fluorescent dyes is now available for use in various aspects of drug transporter studies. The use of these dyes for transporter analyses may, however, be hampered by classic pitfalls of fluorescence technology, such as quenching. Transporter-independent processes such as passive diffusion of dyes through plasma membrane or dye sequestration into subcellular compartments must also be considered, as well as the redundant handling by various distinct transporters of some fluorescent probes. Finally, standardization of dye-based transport assays remains an important on-going issue.

The Use of Transporter Probe Drug Cocktails for the Assessment of Transporter-Based Drug-Drug Interactions in a Clinical Setting-Proposal of a Four Component Transporter Cocktail

Probe drug cocktails are used clinically to assess the potential for drug-drug interactions (DDIs), and in particular, DDIs resulting from coadministration of substrates and inhibitors of cytochrome P450 enzymes. However, a probe drug cocktail has not been identified to assess DDIs involving inhibition of drug transporters. We propose a cocktail consisting of the following substrates to explore the potential for DDIs caused by inhibition of key transporters: digoxin (P-glycoprotein, P-gp), rosuvastatin (breast cancer resistance protein, BCRP; organic anion transporting polypeptides, OATP), metformin (organic cation transporter, OCT; multidrug and toxin extrusion transporters, MATE), and furosemide (organic anion transporter, OAT). Furosemide was evaluated in vitro, and is a substrate of OAT1 and OAT3, with Km values of 38.9 and 21.5 μM, respectively. Furosemide was also identified as a substrate of BCRP, OATP1B1, and OATP1B3. Furosemide inhibited BCRP (50% inhibition of drug transport: 170 μM), but did not inhibit OATP1B1, OATP1B3, OCT2, MATE1, and MATE2-K at concentrations below 300 μM, and P-gp at concentrations below 2000 μM. Conservative approaches for the estimation of the likelihood of in vivo DDIs indicate a remote chance of in vivo transporter inhibition by these probe drugs when administered at low single oral doses. This four component probe drug cocktail is therefore proposed for clinical evaluation.

Unstirred Water Layers and the Kinetics of Organic Cation Transport

PURPOSE

Unstirred water layers (UWLs) present an unavoidable complication to the measurement of transport kinetics in cultured cells, and the high rates of transport achieved by overexpressing heterologous transporters exacerbate the UWL effect. This study examined the correlation between measured Jmax and Kt values and the effect of manipulating UWL thickness or transport Jmax on the accuracy of experimentally determined kinetics of the multidrug transporters, OCT2 and MATE1.

METHODS

Transport of TEA and MPP was measured in CHO cells that stably expressed human OCT2 or MATE1. UWL thickness was manipulated by vigorous reciprocal shaking. Several methods were used to manipulate maximal transport rates.

RESULTS

Vigorous stirring stimulated uptake of OCT2-mediated transport by decreasing apparent Kt (Ktapp) values. Systematic reduction in transport rates was correlated with reduction in Ktapp values. The slope of these relationships indicated a 1500 μm UWL in multiwell plates. Reducing the influence of UWLs (by decreasing either their thickness or the Jmax of substrate transport) reduced Ktapp by 2-fold to >10-fold.

CONCLUSIONS

Failure to take into account the presence of UWLs in experiments using cultured cells to measure transport kinetics can result in significant underestimates of the apparent affinity of multidrug transporters for substrates.

A novel application of t-statistics to objectively assess the quality of IC50 fits for P-glycoprotein and other transporters

Current USFDA and EMA guidance for drug transporter interactions is dependent on IC₅₀ measurements as these are utilized in determining whether a clinical interaction study is warranted. It is therefore important not only to standardize transport inhibition assay systems but also to develop uniform statistical criteria with associated probability statements for generation of robust IC₅₀ values, which can be easily adopted across the industry. The current work provides a quantitative examination of critical factors affecting the quality of IC₅₀ fits for P-gp inhibition through simulations of perfect data with randomly added error as commonly observed in the large data set collected by the P-gp IC₅₀ initiative. The types of errors simulated were (1) variability in replicate measures of transport activity; (2) transformations of error-contaminated transport activity data prior to IC₅₀ fitting (such as performed when determining an IC₅₀ for inhibition of P-gp based on efflux ratio); and (3) the lack of well defined “no inhibition” and “complete inhibition” plateaus. The effect of the algorithm used in fitting the inhibition curve (e.g., two or three parameter fits) was also investigated. These simulations provide strong quantitative support for the recommendations provided in Bentz et al. (2013) for the determination of IC₅₀ values for P-gp and demonstrate the adverse effect of data transformation prior to fitting. Furthermore, the simulations validate uniform statistical criteria for robust IC₅₀ fits in general, which can be easily implemented across the industry. A calibration of the t-statistic is provided through calculation of confidence intervals associated with the t-statistic.

ITC commentary on the prediction of digoxin clinical drug-drug interactions from in vitro transporter assays

The "P-glycoprotein" IC50 working group reported an 18- to 796-fold interlaboratory range in digoxin transport IC50 (inhibitor concentration achieving 50% of maximal inhibition), raising concerns about the predictability of clinical transporter-based drug-drug interactions (DDIs) from in vitro data. This Commentary describes complexities of digoxin transport, which involve both uptake and efflux processes. We caution against attributing digoxin transport IC50 specifically to P-glycoprotein (P-gp) or extending this composite uptake/efflux IC50 variability to individual transporters. Clinical digoxin interaction studies should be interpreted as evaluation of digoxin safety, not P-gp DDIs.

Physiologically based pharmacokinetic modelling and in vivo [I]/K(i) accurately predict P-glycoprotein-mediated drug-drug interactions with dabigatran etexilate

BACKGROUND AND PURPOSE

In vitro inhibitory potency (Ki )-based predictions of P-glycoprotein (P-gp)-mediated drug-drug interactions (DDIs) are hampered by the substantial variability in inhibitory potency. In this study, in vivo-based [I]/Ki values were used to predict the DDI risks of a P-gp substrate dabigatran etexilate (DABE) using physiologically based pharmacokinetic (PBPK) modelling.

EXPERIMENTAL APPROACH

A baseline PBPK model was established with digoxin, a known P-gp substrate. The Km (P-gp transport) of digoxin in the baseline PBPK model was adjusted to Km (i) to fit the change of digoxin pharmacokinetics in the presence of a P-gp inhibitor. Then 'in vivo' [I]/Ki of this P-gp inhibitor was calculated using Km (i) /Km . Baseline PBPK model was developed for DABE, and the 'in vivo' [I]/Ki was incorporated into this model to simulate the static effect of P-gp inhibitor on DABE pharmacokinetics. This approach was verified by comparing the observed and the simulated DABE pharmacokinetics in the presence of five different P-gp inhibitors.

KEY RESULTS

This approach accurately predicted the effects of five P-gp inhibitors on DABE pharmacokinetics (98-133% and 89-104% for the ratios of AUC and Cmax respectively). The effects of 16 other P-gp inhibitors on the pharmacokinetics of DABE were also confidently simulated.

CONCLUSIONS AND IMPLICATIONS

'In vivo' [I]/Ki and PBPK modelling, used in combination, can accurately predict P-gp-mediated DDIs. The described framework provides a mechanistic basis for the proper design of clinical DDI studies, as well as avoiding unnecessary clinical DDI studies.

The need for human breast cancer resistance protein substrate and inhibition evaluation in drug discovery and development: why, when, and how?

Although the multiplicity in transport proteins assessed during drug development is continuously increasing, the clinical relevance of the breast cancer resistance protein (BCRP) is still under debate. Here, our aim is to rationalize the need to consider BCRP substrate and inhibitor interactions and to define optimum selection and acceptance criteria between cell-based and vesicle-based assays in vitro. Information on the preclinical and clinical pharmacokinetics (PK), drug-drug interactions, and pharmacogenomics data was collated for 13 marketed drugs whose PK is reportedly associated with BCRP interaction. Clinical examples where BCRP impacts drug PK and efficacy appear to be rare and confounded by interactions with other transporters. Thirty-seven compounds were selected to be tested as BCRP substrates in a cell-based assay using MDCKII cells (Madin-Darby canine kidney cells) and 18 in membrane vesicles. Depending on the physicochemical compound properties, we observed both in vitro systems to give false-negative readouts. In addition, the inhibition potential of 19 compounds against BCRP was assessed in vesicles and in MDCKII cells, where we observed significant system and substrate-dependent IC50 values. Therefore, neither of the two test systems is superior to the other. Instead, one system may offer advantages under certain situations (e.g., low permeability) and thus should be selected based on the physicochemical compound properties. Finally, given the clinical relevance of BCRP, we propose that its evaluation should remain issue-driven: for low permeable, low bioavailable drugs, in particular when other more common processes do not allow a mechanistic understanding of any unexpected absorption or brain disposition, and for drugs with a low therapeutic window.

Characterization of gastrointestinal absorption of digoxin involving influx and efflux transporter in rats: application of mdr1a knockout (-/-) rats into absorption study of multiple transporter substrate

1.  This study was aimed to characterize gastrointestinal absorption of digoxin using wild-type (WT) and multidrug resistance protein 1a [mdr1a; P-glycoprotein (P-gp)] knockout (-/-) rats. 2.  In WT rats, the area under the plasma concentration-time curve (AUC) of oral digoxin increased after oral pretreatment with quinidine at 30 mg/kg compared with non-treatment, but the increasing ratio tended to decrease at a high dose of 100 mg/kg. In mdr1a (-/-) rats, however, quinidine pretreatment caused a dose-dependent decrease in the AUC. 3.  Quinidine pretreatment did not alter the hepatic availability of digoxin, indicating that the changes in the digoxin AUC were attributable to inhibition of the absorption process by quinidine; i.e. inhibition of influx by quinidine in mdr1a (-/-) rats and inhibition of efflux and influx by quinidine in WT rats. 4.  An in situ rat intestinal closed loop study using naringin implied that organic anion transporting peptide (Oatp) 1a5 may be a responsible transporter in the absorption of digoxin. 5.  These findings imply that the rat absorption behavior of digoxin is possibly governed by Oatp1a5-mediated influx and P-gp-mediated efflux. The mdr1a (-/-) rat is therefore a useful in vivo tool to investigate drug absorption associated with multiple transporters including P-gp.

Quantitative targeted proteomics for membrane transporter proteins: method and application

Although global proteomics has shown promise for discovery of many new proteins, biomarkers, protein modifications, and polymorphisms, targeted proteomics is emerging in the proteomics research field as a complement to untargeted shotgun proteomics, particularly when a determined set of low-abundance functional proteins need to be measured. The function and expression of proteins related to drug absorption, distribution, metabolism, and excretion (ADME) such as cytochrome P450 enzymes and membrane transporters are of great interest in biopharmaceutical research. Since the variation in ADME-related protein expression is known to be a major complicating factor encountered during in vitro-in vivo and in vivo-in vivo extrapolations (IVIVE), the accurate quantification of the ADME proteins in complex biological systems becomes a fundamental element in establishing IVIVE for pharmacokinetic predictions. In this review, we provide an overview of relevant methodologies followed by a summary of recent applications encompassing mass spectrometry-based targeted quantifications of membrane transporters.

Cytochrome p450 inhibitory properties of common efflux transporter inhibitors

Drug transporter inhibitors are important tools to elucidate the contribution of transporters to drug disposition both in vitro and in vivo. These inhibitors are often unselective and affect several transporters as well as drug metabolizing enzymes, which can make experimental results difficult to interpret with confidence. We therefore tested 14 commonly used P-glycoprotein (P-gp), breast cancer resistance protein (BCRP), and multidrug-resistance associated protein (MRP) inhibitors as inhibitors of cytochrome P450 (P450) enzyme activities using recombinant enzymes. A subset of P-gp and/or CYP3A inhibitors were selected (cyclosporin A, elacridar, ketoconazole, quinidine, reserpine, and tacrolimus) for a comparison of P450 inhibition in human microsomes and hepatocytes. Most P-gp inhibitors showed CYP3A4 inhibition, with potencies often in a similar range as their P-gp inhibition, as well as less potent CYP2C19 inhibition. Other P450 enzymes were not strongly inhibited except a few cases of CYP2D6 inhibition. MRP and BCRP inhibitors showed limited P450 inhibition. Some inhibitors showed less P450 inhibition in human hepatocytes than human liver microsomes, for example, elacridar, probably due to differences in binding, permeability limitations, or active, P-gp mediated efflux of the inhibitor from the hepatocytes. Quinidine was a potent P450 inhibitor in hepatocytes but only showed weak inhibition in microsomes. Quinidine shows an extensive cellular uptake, which may potentiate intracellular P450 inhibition. Elacridar, described as a potent and selective P-gp inhibitor, displayed modest P450 inhibition in this study and is thus a useful model inhibitor to define the role of P-gp in drug disposition without interference with other processes.

A structural model for the mass action kinetic analysis of P-gp mediated transport through confluent cell monolayers

The structural model for P-gp mediated transport across confluent cell monolayers uses the generally accepted mass action reactions for P-gp binding and efflux, together with the known structural parameters for P-gp (large substrate binding site accessible from the membrane) and the apical plasma membrane in which it resides (lipid bilayer partition coefficient of substrate and volume of apical plasma membrane allow estimation of substrate concentration at binding site). The model considers binding of substrate to P-gp from within the inner leaflet of the apical membrane, with an on rate constant, k 1 (M(-1)s(-1)), and off rate constant k r (s(-1)), as well as an efflux rate constant from P-gp into the apical chamber, k 2 (s(-1)). The model also explicitly estimates the active P-gp protein level, known as P-gp efflux active surface density T(0). For each new drug, fitting these parameters requires use of multiple initial drug concentrations and multiple time points at each concentration, until steady state is reached between P-gp-mediated efflux into the apical chamber and passive permeability from apical chamber back into the cytosol. Although this model optimally requires a larger than usual dataset for analysis, it does provide important mechanistic information through estimates of these on, off and efflux rate constants, as well as efflux active P-gp surface density. This more detailed description of efflux from polarized confluent cell monolayers has (1) provided insight into the unexpected relationship between P-gp IC50 and K i in this system, (2) highlighted the kinetic need for GF120918 inhibitable apical and basolateral uptake transporters for digoxin, and (3) provided possible explanations for the extreme lab-to-lab variability in P-gp IC50 values observed for inhibition of digoxin transport. This model can also be used to distinguish between efflux active P-gp and total apical plasma membrane P-gp, which may be important when P-gp is expressed in a microvillous membrane.

Use of different parameters and equations for calculation of IC₅₀ values in efflux assays: potential sources of variability in IC₅₀ determination

Drug interactions due to efflux transporters may result in one drug increasing or decreasing the systemic exposure of a second drug. The potential for in vivo drug interactions is estimated through in vitro cell assays. Variability in in vitro parameter determination (e.g., IC₅₀ values) among laboratories may lead to different conclusions in in vivo interaction predictions. The objective of this study was to investigate variability in in vitro inhibition potency determination that may be due to calculation methods. In a Caco-2 cell assay, the absorptive and secretive permeability of digoxin was measured in the presence of spironolactone, itraconazole and vardenafil. From the permeability data, the efflux ratio and net secretory flux where calculated for each inhibitor. IC₅₀ values were then calculated using a variety of equations and software programs. All three drugs decreased the secretory transport of digoxin in a concentration-dependent manner while increasing digoxin's absorption to a lesser extent. The resulting IC₅₀ values varied according to the parameter evaluated, whether percent inhibition or percent control was applied, and the computational IC₅₀ equation. This study has shown that multiple methods used to quantitate the inhibition of drug efflux in a cell assay can result in different IC₅₀ values. The variability in the results in this study points to a need to standardize any transporter assay and calculation methods within a laboratory and to validate the assay with a set of known inhibitors and non-inhibitors against a clinically relevant substrate.

In vitro predictability of drug-drug interaction likelihood of P-glycoprotein-mediated efflux of dabigatran etexilate based on [I]2/IC50 threshold

Dabigatran etexilate, an oral, reversible, competitive, and direct thrombin inhibitor, is an in vitro and in vivo substrate of P-glycoprotein (P-gp). Dabigatran etexilate was proposed as an in vivo probe substrate for intestinal P-gp inhibition in a recent guidance on drug-drug interactions (DDI) from the European Medicines Agency (EMA) and the Food and Drug Administration (FDA). We conducted transcellular transport studies across Caco-2 cell monolayers with dabigatran etexilate in the presence of various P-gp inhibitors to examine how well in vitro IC50 data, in combination with mathematical equations provided by regulatory guidances, predict DDI likelihood. From a set of potential P-gp inhibitors, clarithromycin, cyclosporin A, itraconazole, ketoconazole, quinidine, and ritonavir inhibited P-gp-mediated transport of dabigatran etexilate over a concentration range that may hypothetically occur in the intestine. IC50 values of P-gp inhibitors for dabigatran etexilate transport were comparable to those of digoxin, a well established in vitro and in vivo P-gp substrate. However, IC50 values varied depending whether they were calculated from efflux ratios or permeability coefficients. Prediction of DDI likelihood of P-gp inhibitors using IC50 values, the hypothetical concentration of P-gp inhibitors, and the cut-off value recommended by both the FDA and EMA were in line with the DDI occurrence in clinical studies with dabigatran etexilate. However, it has to be kept in mind that validity of the cut-off criteria proposed by the FDA and EMA depends on in vitro experimental systems and the IC50-calculation methods that are employed, as IC50 values are substantially influenced by these factors.

Application of permeability-limited physiologically-based pharmacokinetic models: part II - prediction of P-glycoprotein mediated drug-drug interactions with digoxin

Digoxin is the recommended substrate for assessment of P-glycoprotein (P-gp)-mediated drug-drug interactions (DDIs) in vivo. The overall aim of our study was to investigate the inhibitory potential of both verapamil and norverapamil on the P-gp-mediated efflux of digoxin in both gut and liver. Therefore, a physiologically-based pharmacokinetic (PBPK) model for verapamil and its primary metabolite was developed and validated through the recovery of observed clinical plasma concentration data for both moieties and the reported interaction with midazolam, albeit a cytochrome P450 3A4-mediated DDI. The validated inhibitor model was then used in conjunction with the model developed previously for digoxin. The range of values obtained for the 10 trials indicated that increases in area under the plasma concentration-time curve (AUC) profiles and maximum plasma concentration observed (Cmax ) values of digoxin following administration of verapamil were more comparable with in vivo observations, when P-gp inhibition by the metabolite, norverapamil, was considered as well. The predicted decrease in AUC and Cmax values of digoxin following administration of rifampicin because of P-gp induction was 1.57- (range: 1.42-1.77) and 1.62-fold (range: 1.53-1.70), which were reasonably consistent with observed values of 1.4- and 2.2-fold, respectively. This study demonstrates the application of permeability-limited models of absorption and distribution within a PBPK framework together with relevant in vitro data on transporters to assess the clinical impact of modulated P-gp-mediated efflux by drugs in development.

Application of receiver operating characteristic analysis to refine the prediction of potential digoxin drug interactions

In the 2012 Food and Drug Administration (FDA) draft guidance on drug-drug interactions (DDIs), a new molecular entity that inhibits P-glycoprotein (P-gp) may need a clinical DDI study with a P-gp substrate such as digoxin when the maximum concentration of inhibitor at steady state divided by IC₅₀ ([I₁]/IC₅₀) is ≥0.1 or concentration of inhibitor based on highest approved dose dissolved in 250 ml divide by IC₅₀ ([I₂]/IC₅₀) is ≥10. In this article, refined criteria are presented, determined by receiver operating characteristic analysis, using IC₅₀ values generated by 23 laboratories. P-gp probe substrates were digoxin for polarized cell-lines and N-methyl quinidine or vinblastine for P-gp overexpressed vesicles. Inhibition of probe substrate transport was evaluated using 15 known P-gp inhibitors. Importantly, the criteria derived in this article take into account variability in IC₅₀ values. Moreover, they are statistically derived based on the highest degree of accuracy in predicting true positive and true negative digoxin DDI results. The refined criteria of [I₁]/IC₅₀ ≥ 0.03 and [I₂]/IC₅₀ ≥ 45 and FDA criteria were applied to a test set of 101 in vitro-in vivo digoxin DDI pairs collated from the literature. The number of false negatives (none predicted but DDI observed) were similar, 10 and 12%, whereas the number of false positives (DDI predicted but not observed) substantially decreased from 51 to 40%, relative to the FDA criteria. On the basis of estimated overall variability in IC₅₀ values, a theoretical 95% confidence interval calculation was developed for single laboratory IC₅₀ values, translating into a range of [I₁]/IC₅₀ and [I₂]/IC₅₀ values. The extent by which this range falls above the criteria is a measure of risk associated with the decision, attributable to variability in IC₅₀ values.