Kinetic considerations for the quantitative assessment of efflux activity and inhibition: implications for understanding and predicting the effects of efflux inhibition

Applicability of MDR1 Overexpressing Abcb1KO-MDCKII Cell Lines for Investigating In Vitro Species Differences and Brain Penetration Prediction

Implementing the 3R initiative to reduce animal experiments in brain penetration prediction for CNS-targeting drugs requires more predictive in vitro and in silico models. However, animal studies are still indispensable to obtaining brain concentration and determining the prediction performance of in vitro models. To reveal species differences and provide reliable data for IVIVE, in vitro models are required. Systems overexpressing MDR1 and BCRP are widely used to predict BBB penetration, highlighting the impact of the in vitro system on predictive performance. In this study, endogenous Abcb1 knock-out MDCKII cells overexpressing MDR1 of human, mouse, rat or cynomolgus monkey origin were used. Good correlations between ERs of 83 drugs determined in each cell line suggest limited species specificities. All cell lines differentiated CNS-penetrating compounds based on ERs with high efficiency and sensitivity. The correlation between in vivo and predicted Kp,uu,brain was the highest using total ER of human MDR1 and BCRP and optimized scaling factors. MDR1 interactors were tested on all MDR1 orthologs using digoxin and quinidine as substrates. We found several examples of inhibition dependent on either substrate or transporter abundance. In summary, this assay system has the potential for early-stage brain penetration screening. IC50 comparison between orthologs is complex; correlation with transporter abundance data is not necessarily proportional and requires the understanding of modes of transporter inhibition.

Translatability of in vitro Inhibition Potency to in vivo P-Glycoprotein Mediated Drug Interaction Risk

P-glycoprotein (P-gp) may limit oral drug absorption of substrate drugs due to intestinal efflux. Therefore, regulatory agencies require investigation of new chemical entities as possible inhibitors of P-gp in vitro. Unfortunately, inter-laboratory and inter-assay variability have hindered the translatability of in vitro P-gp inhibition data to predict clinical drug interaction risk. The current study was designed to evaluate the impact of potential IC₅₀ discrepancies between two commonly utilized assays, i.e., bi-directional Madin-Darby Canine Kidney-MDR1 cell-based and MDR1 membrane vesicle-based assays. When comparing vesicle- to cell-based IC₅₀ values (n = 28 inhibitors), non-P-gp substrates presented good correlation between assay formats, whereas IC₅₀s of P-gp substrates were similar or lower in the vesicle assays. The IC₅₀s obtained with a cell line expressing relatively low P-gp aligned more closely to those obtained from the vesicle assay, but passive permeability of the inhibitors did not appear to influence the correlation of IC₅₀s, suggesting that efflux activity reduces intracellular inhibitor concentrations. IC₅₀s obtained between two independent laboratories using the same assay type showed good correlation. Using the G-value (i.e., ratio of estimated gut concentration-to-inhibition potency) >10 cutoff recommended by regulatory agencies resulted in minimal differences in predictive performance, suggesting this cutoff is appropriate for either assay format.

Examining the efficacy of localised gemcitabine therapy for the treatment of pancreatic cancer using a hybrid agent-based model

The prognosis for pancreatic ductal adenocarcinoma (PDAC) patients has not significantly improved in the past 3 decades, highlighting the need for more effective treatment approaches. Poor patient outcomes and lack of response to therapy can be attributed, in part, to a lack of uptake of perfusion of systemically administered chemotherapeutic drugs into the tumour. Wet-spun alginate fibres loaded with the chemotherapeutic agent gemcitabine have been developed as a potential tool for overcoming the barriers in delivery of systemically administrated drugs to the PDAC tumour microenvironment by delivering high concentrations of drug to the tumour directly over an extended period. While exciting, the practicality, safety, and effectiveness of these devices in a clinical setting requires further investigation. Furthermore, an in-depth assessment of the drug-release rate from these devices needs to be undertaken to determine whether an optimal release profile exists. Using a hybrid computational model (agent-based model and partial differential equation system), we developed a simulation of pancreatic tumour growth and response to treatment with gemcitabine loaded alginate fibres. The model was calibrated using in vitro and in vivo data and simulated using a finite volume method discretisation. We then used the model to compare different intratumoural implantation protocols and gemcitabine-release rates. In our model, the primary driver of pancreatic tumour growth was the rate of tumour cell division. We were able to demonstrate that intratumoural placement of gemcitabine loaded fibres was more effective than peritumoural placement. Additionally, we quantified the efficacy of different release profiles from the implanted fibres that have not yet been tested experimentally. Altogether, the model developed here is a tool that can be used to investigate other drug delivery devices to improve the arsenal of treatments available for PDAC and other difficult-to-treat cancers in the future.

Mechanistic Modeling of Central Nervous System Pharmacokinetics and Target Engagement of HER2 Tyrosine Kinase Inhibitors to Inform Treatment of Breast Cancer Brain Metastases

PURPOSE

This study evaluated the central nervous system (CNS) pharmacokinetics and target engagement of lapatinib, neratinib, and tucatinib in patients with cancer, using a physiologically based pharmacokinetic (PBPK) modeling approach.

EXPERIMENTAL DESIGN

Drug-specific parameters for in vitro metabolism, binding to plasma proteins and brain tissues, transcellular passive permeability, and interactions with efflux transporters were determined. Whole-body PBPK models integrated with a 4-compartment permeability-limited brain model was developed and verified for predicting plasma and CNS pharmacokinetics. Target engagement ratio (TER), defined as the ratio of the average steady-state unbound drug brain concentration (Css,ave,br) to in vitro IC50 for HER2 inhibition, was used as a predictor of intracranial efficacy.

RESULTS

PBPK models predicted that following 1 cycle of standard dosing, tucatinib and lapatinib achieved similar Css,ave,br (14.5 vs. 16.8 nmol/L), while neratinib Css,ave,br (0.68 nmol/L) was 20-fold lower. Tucatinib and neratinib were equally potent for HER2 inhibition (IC50, 6.9 vs. 5.6 nmol/L), while lapatinib was less potent (IC50, 109 nmol/L). The model-predicted population mean TER in the human normal brain was 2.1 for tucatinib, but < 0.20 for lapatinib and neratinib.

CONCLUSIONS

The PBPK modeling suggests that tucatinib induces sufficient HER2 inhibition (TER > 2.0) in not only brain metastases with a disrupted blood-brain barrier (BBB), but also micrometastases where the BBB largely remains intact. These findings, in line with available clinical pharmacokinetics and efficacy data, support the therapeutic value of tucatinib for treatment of brain metastases and warrant further clinical investigation for the prevention of brain metastases in patients with HER2-positive breast cancer.

In Vitro–In Silico Modeling of Caffeine and Diclofenac Permeation in Static and Fluidic Systems with a 16HBE Lung Cell Barrier

Static in vitro permeation experiments are commonly used to gain insights into the permeation properties of drug substances but exhibit limitations due to missing physiologic cell stimuli. Thus, fluidic systems integrating stimuli, such as physicochemical fluxes, have been developed. However, as fluidic in vitro studies display higher complexity compared to static systems, analysis of experimental readouts is challenging. Here, the integration of in silico tools holds the potential to evaluate fluidic experiments and to investigate specific simulation scenarios. This study aimed to develop in silico models that describe and predict the permeation and disposition of two model substances in a static and fluidic in vitro system. For this, in vitro permeation studies with a 16HBE cellular barrier under both static and fluidic conditions were performed over 72 h. In silico models were implemented and employed to describe and predict concentration–time profiles of caffeine and diclofenac in various experimental setups. For both substances, in silico modeling identified reduced apparent permeabilities in the fluidic compared to the static cellular setting. The developed in vitro–in silico modeling framework can be expanded further, integrating additional cell tissues in the fluidic system, and can be employed in future studies to model pharmacokinetic and pharmacodynamic drug behavior.

A Bidirectional Permeability Assay for beyond Rule of 5 Compounds

Bidirectional permeability measurement with cellular models grown on Transwell inserts is widely used in pharmaceutical research since it not only provides information about the passive permeability of a drug, but also about transport proteins involved in the active transport of drug substances across physiological barriers. With the increasing number of investigative drugs coming from chemical space beyond Lipinski’s Rule of 5, it becomes more and more challenging to provide meaningful data with the standard permeability assay. This is exemplified here by the difficulties we encountered with the cyclic depsipeptides emodepside and its close analogs with molecular weight beyond 1000 daltons and cLogP beyond 5. The aim of this study is to identify potential reasons for these challenges and modify the permeability assays accordingly. With the modified assay, intrinsic permeability and in vitro efflux of depsipeptides could be measured reliably. The improved correlation to in vivo bioavailability and tissue distribution data indicated the usefulness of the modified permeability assay for the in vitro screening of compounds beyond the Rule of 5.

Metformin Transport Rates Between Plasma and Red Blood Cells in Humans

Background

Metformin has been used for the treatment of type 2 diabetes for over 60 years; however, its mechanism of pharmacological action is not fully clear. Different hypotheses exist regarding metformin distribution and redistribution mechanisms between plasma and erythrocytes/red blood cells (RBCs).

Objective

We aimed to test the hypothesis that the metformin distribution between plasma and RBC occurs via concentration difference-driven passive transport and estimated transport rate coefficient values based on metformin concentration time series in plasma and RBCs from in vivo studies.

Methods

An ordinary differential equation (ODE) system with two compartments was used to describe diffusion-based passive transport between plasma and RBCs. Metformin concentration time series in plasma and RBCs of 35 individuals were used for metformin transport parametrization. Plasma concentration has been approximated by biexponential decline.

Results

A single passive transport coefficient, k = 0.044 ± 0.014 (h^–1), can be applied, describing the uptake and release transport rate versus the linear equation v = k × (M_pl − M_RBC), where M_pl is the metformin concentration in plasma and M_RBC is the metformin concentration in RBCs.

Conclusions

Our research suggests that passive transport can explain metformin distribution dynamics between plasma and RBCs because transport speed is proportional to the metformin concentration difference and independent of the transport direction. Concentration difference-driven passive transport can explain the mechanism of faster metformin distribution to RBCs the first few hours after administration, and faster release and domination of the redistribution transport rate after metformin concentration in plasma becomes smaller than in RBCs.

Supplementary Information

The online version contains supplementary material available at 10.1007/s40262-021-01058-2.

Mechanistic Model for the Prediction of Small-Molecule Vitreal Clearance Combining Diffusion-Limited and Permeability-Limited Clearance

The discovery of new small-molecule drugs for intravitreal administration would benefit from simple models to predict vitreal clearance (CL). The current models available have limitations in their applicability to small-molecule drugs and translatability to humans. We developed a mechanistic model combining the diffusion rate of the molecule in the vitreous and permeability across posterior segment tissues and applied it to 30 small molecules with observed CL available mostly from literature. We used Caco-2 permeability as a surrogate for ocular tissue permeability. The model predicted rabbit vitreal CL well, with 80% of the predictions being within a 2-fold range of the observed CL. For an accurate prediction, it was crucial to consider the anterior diffusion CL from the vitreous to the aqueous and a limiting diffusion CL for the whole eye. We observed no major differences in model accuracy when using literature permeability values from retinal pigment epithelial cell models. Importantly, by adopting the specific dimensions of the human eye, the model was able to accurately predict vitreal CL of four compounds for which human vitreal CL data are available. In summary, this mechanistic model enables a simple, accurate, and translatable estimation of small-molecule vitreal CL.

The role of drug efflux and uptake transporters ABCB1 (P-gp), ABCG2 (BCRP) and OATP1A/1B and of CYP3A4 in the pharmacokinetics of the CDK inhibitor milciclib

The promising anticancer drug milciclib potently inhibits cyclin-dependent kinase (CDK) 2 and tropomyosin receptor kinase (TRK) A, and is currently in phase II clinical studies. To characterize factors affecting milciclib pharmacokinetics, we investigated whether milciclib is a substrate of the multidrug efflux and uptake transporters ABCB1 (P-gp), ABCG2 (BCRP), and OATP1A/1B, and the drug-metabolizing enzyme CYP3A, using genetically-modified mouse models and Madin-Darby Canine Kidney (MDCK-II) cells. In vitro, milciclib was transported by mAbcg2, and this was inhibited by the ABCG2 inhibitor Ko143. Upon oral administration of milciclib, its plasma exposure in Abcb1a/1b^-/-, Abcg2^-/-, and Abcb1a/1b;Abcg2^-/- mice was similar to that found in wild-type mice. Milciclib showed good brain penetration even in wild-type mice (brain-to-plasma ratio of 1.2), but this was further increased by 5.2-fold when both Abcb1 and Abcg2 were ablated, and to a lesser extent in single Abcb1- or Abcg2-deficient mice. Oatp1a/1b deficiency had only a minor impact on the milciclib plasma AUC_0-24h and C_max. The milciclib AUC_0-8h increased 1.9-fold in Cyp3a^-/- mice but decreased only 1.3-fold upon overexpression of human CYP3A4. Thus, ABCB1 and ABCG2 cooperatively limit milciclib brain penetration. The low impact of OATP1 and CYP3A could be clinically favorable for milciclib, reducing the risks of unintended drug-drug interactions or interindividual variation in CYP3A4 activity.

Improved In Vitro-In Vivo Correlation by Using the Unbound-Fraction-Adjusted IC50 for Breast Cancer Resistance Protein Inhibition

PURPOSE

One function of the blood-brain barrier (BBB) is the efflux of xenobiotics by breast cancer resistance protein (BCRP), and inhibition of BCRP can cause unexpected central nervous system toxicity. Despite the importance of BCRP inhibition and the associated risk of BBB penetration in vivo, there has been little investigation of it to date. In this study, inhibition of BCRP-mediated transport was assessed by in vitro assay in the presence of bovine serum albumin (BSA) to change the unbound inhibitor concentrations, and the in vitro-in vivo correlation (IVIVC) at the BBB was evaluated.

METHODS AND RESULTS

The IC₅₀ values of BCRP inhibitors were determined in vitro with and without BSA and the inhibitors were categorized into two groups. One group of compounds had little risk of inhibiting BCRP because of their low unbound concentrations. In contrast, the other group has the potential to facilitate BBB penetration by inhibiting BCRP. In the IVIVC approach, brain concentrations and the brain-to-plasma ratio were better correlated with the ratio of the unbound plasma concentration at steady-state to the unbound-fraction-adjusted IC₅₀.

CONCLUSION

We have found a way to obtain a better in vitro-in vivo correlation for BCRP-mediated transport.

Novel insights into the organic solute transporter alpha/beta, OSTα/β: From the bench to the bedside

Organic solute transporter alpha/beta (OSTα/β) is a heteromeric solute carrier protein that transports bile acids, steroid metabolites and drugs into and out of cells. OSTα/β protein is expressed in various tissues, but its expression is highest in the gastrointestinal tract where it facilitates the recirculation of bile acids from the gut to the liver. Previous studies established that OSTα/β is upregulated in liver tissue of patients with extrahepatic cholestasis, obstructive cholestasis, and primary biliary cholangitis (PBC), conditions that are characterized by elevated bile acid concentrations in the liver and/or systemic circulation. The discovery that OSTα/β is highly upregulated in the liver of patients with nonalcoholic steatohepatitis (NASH) further highlights the clinical relevance of this transporter because the incidence of NASH is increasing at an alarming rate with the obesity epidemic. Since OSTα/β is closely linked to the homeostasis of bile acids, and tightly regulated by the nuclear receptor farnesoid X receptor, OSTα/β is a potential drug target for treatment of cholestatic liver disease, and other bile acid-related metabolic disorders such as obesity and diabetes. Obeticholic acid, a semi-synthetic bile acid used to treat PBC, under review for the treatment of NASH, and in development for the treatment of other metabolic disorders, induces OSTα/β. Some drugs associated with hepatotoxicity inhibit OSTα/β, suggesting a possible role for OSTα/β in drug-induced liver injury (DILI). Furthermore, clinical cases of homozygous genetic defects in both OSTα/β subunits resulting in diarrhea and features of cholestasis have been reported. This review article has been compiled to comprehensively summarize the recent data emerging on OSTα/β, recapitulating the available literature on the structure-function and expression-function relationships of OSTα/β, the regulation of this important transporter, the interaction of drugs and other compounds with OSTα/β, and the comparison of OSTα/β with other solute carrier transporters as well as adenosine triphosphate-binding cassette transporters. Findings from basic to more clinically focused research efforts are described and discussed.

Facts, Patterns, and Principles in Drug Discovery: Appraising the Rule of 5 with Measured Physicochemical Data

The rule of 5 was designed to estimate the likelihood of poor absorption or permeation, noting the impact of poor solubility. This Perspective explores the impact of various physicochemical descriptors and contemporary lipophilicity measurements on permeability and solubility, showing that the distribution coefficient log D_7.4 (rather than log P) is the most impactful parameter. Molecular weight, almost invariably the defining characteristic of "beyond the rule of 5" compounds, has little impact on solubility when log D_7.4 measurements and aromaticity are considered. Predicting permeation is more complex, given passive and carrier transport mechanisms; however, notable patterns of behavior are apparent, giving insight even "beyond the rule of 5". Recommended best practices should involve using the facts (measurements) and the patterns they reveal to establish informative principles rather than fastidious rules.

Simplified Method to Determine the Efflux Ratio on P-Glycoprotein Substrates Using Three-Compartment Model Analysis for Caco-2 Cell Assay Data

PURPOSE

Multiple time-point sampling is required in transcellular transport studies to accurately calculate the appropriate efflux ratio (ER). Our study sought to develop a simplified method to determine the ER in Caco-2 cells.

METHODS

The equation for the ER was derived from a three-compartment model of apical to basal and basal to apical transport. Transcellular transport studies were conducted with 10 non-P-glycoprotein (P-gp) and 6 P-gp substrates in Caco-2 cells, and the ER was calculated using this equation.

RESULTS

The equation for the ER used the concentration ratio in the receiver compartment at the same time-point; therefore, the ER can theoretically be calculated using only a single point. The ER of all non-P-gp substrates tested was close to 1 at all sampling times. The ERs of cyclosporine A calculated from the concentration ratio at 30, 60, 90, and 120 min incubation were 2.93, 6.43, 7.12, and 9.57, respectively, and the ER at 120 min was almost identical to the theoretical value (9.62) calculated using three-compartment model analysis. The other 5 P-gp substrates showed a similar tendency. Single-point sampling can be used to accurately calculate ER at 120 min.

CONCLUSIONS

Single-point sampling is a promising approach for calculating appropriate ERs in the drug discovery stage.

Examining P-gp efflux kinetics guided by the BDDCS - Rational selection of in vitro assay designs and mathematical models

The generation of reliable kinetic parameters to describe P-glycoprotein (P-gp) activity is essential for predicting the impact of efflux transport on gastrointestinal drug absorption. The compound-specific selection of in vitro assay designs and ensuing data analysis methods is explored in this manuscript. We measured transcellular permeability and cellular uptake of five P-gp substrates in Caco-2 and LLC-PK1 MDR1 cells. Kinetic parameters of P-gp-mediated efflux transport (K_m, V_max) were derived from conventional and mechanistic compartmental models. The estimated apparent K_m values based on medium concentrations in the conventional permeability model indicated significant differences between the cell lines. The respective intrinsic K_m values based on unbound intracellular concentrations in the mechanistic compartmental models were significantly lower and comparable between cell lines and assay formats. Non-specific binding or lysosomal trapping were shown to cause discrepancies in the kinetic parameters obtained from different assay formats. A guidance for the selection of in vitro assays and kinetic assessment methods is proposed in line with the Biopharmaceutics Drug Disposition Classification System (BDDCS). The recommendations are expected to aid the acquisition of robust and reproducible kinetic parameters of P-gp-mediated efflux transport.

Novel in Vitro Method Reveals Drugs That Inhibit Organic Solute Transporter Alpha/Beta (OSTα/β)

Drug interactions with the organic solute transporter alpha/beta (OSTα/β) are understudied even though OSTα/β is an important transporter that is expressed in multiple human tissues including the intestine, kidneys, and liver. In this study, an in vitro method to identify novel OSTα/β inhibitors was first developed using OSTα/β-overexpressing Flp-In 293 cells. Incubation conditions were optimized using previously reported OSTα/β inhibitors. A method including a 10 min preincubation step with the test compound was used to screen for OSTα/β inhibition by 77 structurally diverse compounds and fixed-dose combinations. Seven compounds and one fixed-dose combination (100 μM final concentration) inhibited OSTα/β-mediated dehydroepiandrosterone sulfate (DHEAS) uptake by >25%. Concentration-dependent OSTα/β inhibition was evaluated for all putative inhibitors (atorvastatin, ethinylestradiol, fidaxomicin, glycochenodeoxycholate, norgestimate, troglitazone, and troglitazone sulfate). Ethinylestradiol, fidaxomicin, and troglitazone sulfate yielded a clear concentration-inhibition response with IC50 values <200 μM. Among all tested compounds, there was no clear association between physicochemical properties, the severity of hepatotoxicity, and the degree of OSTα/β inhibition. This study utilized a novel in vitro method to identify OSTα/β inhibitors and, for the first time, provided IC50 values for OSTα/β inhibition. These data provide evidence that several drugs, some of which are associated with cholestatic drug-induced liver injury, may impair the function of the OSTα/β transporter.

P-glycoprotein and breast cancer resistance protein restrict brigatinib brain accumulation and toxicity, and, alongside CYP3A, limit its oral availability

Brigatinib is an FDA-approved oral anaplastic lymphoma kinase (ALK) inhibitor for treatment of metastatic non-small cell lung cancer (NSCLC). Using genetically modified mouse models, we investigated the roles of the multidrug efflux transporters ABCB1 and ABCG2, and the multispecific drug-metabolizing enzyme CYP3 A in plasma pharmacokinetics and tissue distribution of brigatinib. In vitro, brigatinib was exceptionally well transported by human ABCB1 and mouse Abcg2, and efficiently by human ABCG2. Following oral brigatinib administration (10 mg/kg), brain accumulation was dramatically increased in Abcb1a/1b^-/- (19.3-fold) and Abcb1a/1b;Abcg2^-/-(41.8-fold), but not in single Abcg2^-/- mice compared to wild-type mice. Brigatinib testis accumulation showed qualitatively similar behavior. mAbcb1a/1b and mAbcg2 together restricted systemic exposure of brigatinib: with both systems absent oral availability increased 1.9-fold. Coadministration of elacridar, an ABCB1/ABCG2 inhibitor, caused a pronounced increase (36-fold) in brain-to-plasma ratios of brigatinib, approaching the levels seen in Abcb1a/1b;Abcg2^-/- mice. Unexpectedly, lethal toxicity of oral brigatinib was observed in mice with genetic knockout or pharmacological inhibition of mAbcb1a/1b and mAbcg2, indicating a pronounced protective role for these transporters. In Cyp3a^-/- mice, brigatinib plasma exposure increased 1.3-fold, and was subsequently 1.8-fold reduced by transgenic overexpression of human CYP3 A4 in liver and intestine. The relative tissue distribution of brigatinib, however, remained unaltered. ABCB1 and ABCG2 thus limit brain accumulation, toxicity, and systemic exposure of brigatinib, whereas CYP3 A also markedly restricts its oral availability. Unexpected toxicities should therefore be carefully monitored when brigatinib is coadministered with ABCB1/ABCG2 inhibitors in patients. Collectively, these insights may support the clinical application of brigatinib.

Prediction of brain:blood unbound concentration ratios in CNS drug discovery employing in silico and in vitro model systems

Recent years have seen a paradigm shift away from optimizing the brain:blood concentration ratio toward the more relevant brain:blood unbound concentration ratio (K_p,uu,br) in CNS drug discovery. Here, we review the recent developments in the in silico and in vitro model systems to predict the K_p,uu,br of discovery compounds with special emphasis on the in-vitro-in-vivo correlation. We also discuss clinical 'translation' of rodent K_p,uu,br and highlight the future directions for improvement in brain penetration prediction. Important in this regard are in silico K_p,uu,br models built on larger datasets of high quality, calibration and deeper understanding of experimental in vitro transporter systems, and better understanding of blood-brain barrier transporters and their in vivo relevance aside from P-gp and BCRP.

P-glycoprotein (MDR1/ABCB1) and Breast Cancer Resistance Protein (BCRP/ABCG2) affect brain accumulation and intestinal disposition of encorafenib in mice

Encorafenib (LGX818) is a promising BRAF^V600E inhibitor that has efficacy against metastatic melanoma. To better understand its pharmacokinetics, we studied its interactions with the multidrug efflux transporters ABCB1 and ABCG2 and the multidrug metabolizing enzyme CYP3A. In polarized MDCK-II cells, encorafenib was efficiently transported by canine and human ABCB1 and ABCG2 and by mouse Abcg2. Upon oral administration to wild-type, Abcb1a/1b^-/-, Abcg2^-/-, and Abcb1a/1b;Abcg2^-/- mice, encorafenib was absorbed very quickly and to very high plasma levels, but without clear changes in oral availability between the strains. Upon oral or intravenous administration, encorafenib brain accumulation was markedly increased in Abcb1a/1b;Abcg2^-/- mice and to a lesser extent in Abcb1a/1b^-/- mice. However, absolute brain concentrations and brain-to-plasma ratios remained very low in all strains, indicating intrinsically poor brain penetration of encorafenib. Upon intravenous administration, Abcb1a/1b;Abcg2^-/- mice showed somewhat reduced plasma elimination of encorafenib compared to wild-type mice, and lower accumulation of the drug in the intestinal tract, suggesting a limited role for these transporters in intestinal elimination of the drug. In Cyp3a^-/- mice plasma levels of encorafenib were not markedly increased, suggesting a limited impact of Cyp3a on encorafenib oral availability. The low brain penetration of encorafenib might limit its efficacy against malignancies positioned behind a functional blood-brain barrier, but its oral bioavailability and distribution to other tested organs (liver, kidney, spleen, testis) was high.

Quantitative and Mechanistic Understanding of AZD1775 Penetration across Human Blood-Brain Barrier in Glioblastoma Patients Using an IVIVE-PBPK Modeling Approach

Purpose: AZD1775, a first-in-class, small-molecule inhibitor of the Wee1 tyrosine kinase, is under evaluation as a potential chemo- and radiosensitizer for treating glioblastoma. This study was to prospectively, quantitatively, and mechanistically investigate the penetration of AZD1775 across the human blood-brain barrier (BBB).Experimental Design: AZD1775 plasma and tumor pharmacokinetics were evaluated in 20 patients with glioblastoma. The drug metabolism, transcellular passive permeability, and interactions with efflux and uptake transporters were determined using human derived in vitro systems. A whole-body physiologically based pharmacokinetic (PBPK) model integrated with a four-compartment permeability-limited brain model was developed for predicting the kinetics of AZD1775 BBB penetration and assessing the factors modulating this process.Results: AZD1775 exhibited good tumor penetration in patients with glioblastoma, with the unbound tumor-to-plasma concentration ratio ranging from 1.3 to 24.4 (median, 3.2). It was a substrate for ABCB1, ABCG2, and OATP1A2, but not for OATP2B1 or OAT3. AZD1775 transcellular passive permeability and active efflux clearance across MDCKII-ABCB1 or MDCKII-ABCG2 cell monolayers were dependent on the basolateral pH. The PBPK model well predicted observed drug plasma and tumor concentrations in patients. The extent and rate of drug BBB penetration were influenced by BBB integrity, efflux and uptake active transporter activity, and drug binding to brain tissue.Conclusions: In the relatively acidic tumor microenvironment where ABCB1/ABCG2 transporter-mediated efflux clearance is reduced, OATP1A2-mediated active uptake becomes dominant, driving AZD1775 penetration into brain tumor. Variations in the brain tumor regional pH, transporter expression/activity, and BBB integrity collectively contribute to the heterogeneity of AZD1775 penetration into brain tumors. Clin Cancer Res; 23(24); 7454-66. ©2017 AACRSee related commentary by Peer et al., p. 7437.

Blood-Brain Barrier Driven Pharmacoresistance in Amyotrophic Lateral Sclerosis and Challenges for Effective Drug Therapies

The blood-brain barrier (BBB) is essential for proper neuronal function, homeostasis, and protection of the central nervous system (CNS) microenvironment from blood-borne pathogens and neurotoxins. The BBB is also an impediment for CNS penetration of drugs. In some neurologic conditions, such as epilepsy and brain tumors, overexpression of P-glycoprotein, an efflux transporter whose physiological function is to expel catabolites and xenobiotics from the CNS into the blood stream, has been reported. Recent studies reported that overexpression of P-glycoprotein and increase in its activity at the BBB drives a progressive resistance to CNS penetration and persistence of riluzole, the only drug approved thus far for treatment of amyotrophic lateral sclerosis (ALS), rapidly progressive and mostly fatal neurologic disease. This review will discuss the impact of transporter-mediated pharmacoresistance for ALS drug therapy and the potential therapeutic strategies to improve the outcome of ALS clinical trials and efficacy of current and future drug treatments.

A CRISPR-Cas9 Generated MDCK Cell Line Expressing Human MDR1 Without Endogenous Canine MDR1 (cABCB1): An Improved Tool for Drug Efflux Studies

Madin-Darby canine kidney (MDCK) II cells stably transfected with transport proteins are commonly used models for drug transport studies. However, endogenous expression of especially canine MDR1 (cMDR1) confounds the interpretation of such studies. Here we have established an MDCK cell line stably overexpressing the human MDR1 transporter (hMDR1; P-glycoprotein), and used CRISPR-Cas9 gene editing to knockout the endogenous cMDR1. Genomic screening revealed the generation of a clonal cell line homozygous for a 4-nucleotide deletion in the canine ABCB1 gene leading to a frameshift and a premature stop codon. Knockout of cMDR1 expression was verified by quantitative protein analysis and functional studies showing retained activity of the human MDR1 transporter. Application of this cell line allowed unbiased reclassification of drugs previously defined as both substrates and non-substrates in different studies using commonly used MDCK-MDR1 clones. Our new MDCK-hMDR1 cell line, together with a previously developed control cell line, both with identical deletions in the canine ABCB1 gene and lack of cMDR1 expression represent excellent in vitro tools for use in drug discovery.

Breast cancer resistance protein (BCRP/ABCG2) and P-glycoprotein (P-gp/ABCB1) transport afatinib and restrict its oral availability and brain accumulation

Afatinib is a highly selective, irreversible inhibitor of EGFR and HER-2. It is orally administered for the treatment of patients with EGFR mutation-positive types of metastatic NSCLC. We investigated whether afatinib is a substrate for the multidrug efflux transporters ABCB1 and ABCG2 and whether these transporters influence oral availability and brain and other tissue accumulation of afatinib. We used in vitro transport assays to assess human (h)ABCB1-, hABCG2- or murine (m)Abcg2-mediated transport of afatinib. To study the single and combined roles of Abcg2 and Abcb1a/1b in oral afatinib disposition, we used appropriate knockout mouse strains. Afatinib was transported well by hABCB1, hABCG2 and mAbcg2 in vitro. Upon oral administration of afatinib, Abcg2^-/-, Abcb1a/1b^-/- and Abcb1a/1b^-/-;Abcg2^-/- mice displayed a 4.2-, 2.4- and 7-fold increased afatinib plasma AUC_0-24 compared with wild-type mice. Abcg2-deficient strains also displayed decreased afatinib plasma clearance. At 2h, relative brain accumulation of afatinib was not significantly altered in the single knockout strains, but 23.8-fold increased in Abcb1a/1b^-/-;Abcg2^-/- mice compared to wild-type mice. Abcg2 and Abcb1a/1b restrict oral availability and brain accumulation of afatinib. Inhibition of these transporters may therefore be of clinical importance for patients with brain (micro)metastases positioned behind an intact blood-brain barrier.

Key to Opening Kidney for In Vitro-In Vivo Extrapolation Entrance in Health and Disease: Part I: In Vitro Systems and Physiological Data

The programme for the 2015 AAPS Annual Meeting and Exhibition (Orlando, FL; 25-29 October 2015) included a sunrise session presenting an overview of the state-of-the-art tools for in vitro-in vivo extrapolation (IVIVE) and mechanistic prediction of renal drug disposition. These concepts are based on approaches developed for prediction of hepatic clearance, with consideration of scaling factors physiologically relevant to kidney and the unique and complex structural organisation of this organ. Physiologically relevant kidney models require a number of parameters for mechanistic description of processes, supported by quantitative information on renal physiology (system parameters) and in vitro/in silico drug-related data. This review expands upon the themes raised during the session and highlights the importance of high quality in vitro drug data generated in appropriate experimental setup and robust system-related information for successful IVIVE of renal drug disposition. The different in vitro systems available for studying renal drug metabolism and transport are summarised and recent developments involving state-of-the-art technologies highlighted. Current gaps and uncertainties associated with system parameters related to human kidney for the development of physiologically based pharmacokinetic (PBPK) model and quantitative prediction of renal drug disposition, excretion, and/or metabolism are identified.

Cell permeability beyond the rule of 5

Drug discovery for difficult targets that have large and flat binding sites is often better suited to compounds beyond the "rule of 5" (bRo5). However, such compounds carry higher pharmacokinetic risks, such as low solubility and permeability, and increased efflux and metabolism. Interestingly, recent drug approvals and studies suggest that cell permeable and orally bioavailable drugs can be discovered far into bRo5 space. Tactics such as reduction or shielding of polarity by N-methylation, bulky side chains and intramolecular hydrogen bonds may be used to increase cell permeability in this space, but often results in decreased solubility. Conformationally flexible compounds can, however, combine high permeability and solubility, properties that are keys for cell permeability and intestinal absorption. Recent developments in computational conformational analysis will aid design of such compounds and hence prediction of cell permeability. Transporter mediated efflux occurs for most investigated drugs in bRo5 space, however it is commonly overcome by high local intestinal concentrations on oral administration. In contrast, there is little data to support significant impact of transporter-mediated intestinal absorption in bRo5 space. Current knowledge of compound properties that govern transporter effects of bRo5 drugs is limited and requires further fundamental and comprehensive studies.

Blood-brain and retinal barriers show dissimilar ABC transporter impacts and concealed effect of P-glycoprotein on a novel verapamil influx carrier

BACKGROUND AND PURPOSE

The respective impact and interplay between ABC (P-glycoprotein/P-gp/Abcb1a, BCRP/ABCG2, MRP/ABCC) and SLC transporter functions at the blood-brain barrier (BBB) and blood-retinal barriers (BRB) are incompletely understood.

EXPERIMENTAL APPROACH

We measured the initial cerebral and retinal distribution of selected ABC substrates by in situ carotid perfusion using P-gp/Bcrp knockout mice and chemical ABC/SLC modulation strategies. P-gp, Bcrp, Mrp1 and Mrp4 were studied by confocal retina imaging.

KEY RESULTS

Chemical or physical disruption of P-gp increased [(3) H]-verapamil transport by ~10-fold at the BBB and ~1.5-fold at the BRB. [(3) H]-Verapamil transport involved influx-mediated by an organic cation clonidine-sensitive/diphenhydramine-sensitive proton antiporter at both barriers; this effect was unmasked when P-gp was partially or fully inhibited/disrupted at the BBB. Studies of [(3) H]-mitoxantrone and [(3) H]-zidovudine transport suggested, respectively, that Bcrp efflux was less involved at the BRB than BBB, whereas Mrps were significantly and similarly involved at both barriers. Confocal imaging showed that P-gp and Bcrp were expressed in intra-retinal vessels (inner BRB/iBRB) but absent from the blood/basal membrane of cells of the retinal pigment epithelium (outer BRB/oBRB/RPE) where, in contrast, Mrp1 and Mrp4 were localized.

CONCLUSIONS AND IMPLICATIONS

P-gp, Bcrp, Mrp1 and Mrp4 are differentially expressed at the outer and inner BRB, resulting in an altered ability to limit substrate distribution at the retina as compared with the BBB. [(3) H]-Verapamil distribution is not P-gp-specific and involves a proton antiporter at both the BBB and BRB. However, this transport is concealed by P-gp at the BBB, but not at the BRB, where P-gp activity is reduced.

Brain accumulation of the EML4-ALK inhibitor ceritinib is restricted by P-glycoprotein (P-GP/ABCB1) and breast cancer resistance protein (BCRP/ABCG2)

We aimed to clarify the roles of the multidrug transporters ABCB1 and ABCG2 in oral availability and brain accumulation of ceritinib, an oral anaplastic lymphoma kinase (ALK) inhibitor used to treat metastatic non-small cell lung cancer (NSCLC) after progression on crizotinib. Importantly, NSCLC is prone to form brain metastases. Transport of ceritinib by human (h) ABCB1 or hABCG2 or mouse (m) Abcg2 was assessed in vitro. To study the single and combined roles of Abcb1a/1b and Abcg2 in ceritinib disposition in vivo, we used appropriate knockout mouse strains. Ceritinib was very efficiently transported by hABCB1, and efficiently by hABCG2 and mAbcg2 in vitro, and transport was specifically inhibited by the ABCB1 inhibitor zosuquidar and ABCG2 inhibitor Ko143, respectively. Absorption and 24-h oral availability were not significantly affected by the absence of Abcb1 and/or Abcg2, but the brain concentrations were greatly increased (>38-fold) in Abcb1a/1b(-/-) mice at 3 and 24h after oral administration of 20mg/kg ceritinib. The brain concentrations increased another ∼ 3-fold (to >90-fold) in Abcb1a/1b;Abcg2(-/-) mice, indicating that there was a significant additional effect of Abcg2-mediated transport of ceritinib as well in vivo. Overall, brain accumulation, but not the 24-h oral availability of ceritinib were profoundly restricted by Abcb1a/1b and Abcg2, with Abcb1a/1b being the dominant efflux protein. Our data suggest that coadministration of ceritinib with a dual ABCB1 and ABCG2 inhibitor may improve treatment of brain (micro) metastases positioned behind a functionally intact blood-brain barrier, and possibly also of tumors resistant to ceritinib due to ABCB1 or ABCG2 overexpression.

Metformin's Intrinsic Blood-to-Plasma Partition Ratio (B/P): Reconciling the Perceived High In Vivo B/P > 10 with the In Vitro Equilibrium Value of Unity

Blood cells are considered an important distributional compartment for metformin based on the high blood-to-plasma partition ratio (B/P) in humans (>10 at Cmin). However, literature reports of metformin's intrinsic in vitro B/P values are lacking. At present, the extent and rate of metformin cellular partitioning was determined in incubations of fresh human and rat blood with [(14)C]metformin for up to 1 week at concentrations spanning steady-state plasma Cmin, Cmax, and a concentration associated with lactic acidosis. The results showed that metformin's intrinsic equilibrium B/P was ∼0.8-1.4 in blood, which is <10% of the reported clinical value. Kinetics of metformin partitioning into human blood cells and repartitioning back into plasma were slow (repartitioning half-life ∼32-39 hours). These data, along with in vivo rapid and efficient renal clearance of plasma metformin (plasma renal extraction ratio ∼90%-100%), explain why the clinical terminal half-life of metformin in plasma (6 hours) is 3- to 4-fold shorter than the half-life in whole blood (18 hours) and erythrocytes (23 hours). The rate constant for metformin repartitioning from blood cells to plasma (∼0.02 h(-1)) is far slower than the clinical renal elimination rate constant (0.3 h(-1)). Blood distributional rate constants were incorporated into a metformin physiologically-based pharmacokinetic model, which predicted the differential elimination half-life in plasma and blood. The present study demonstrates that the extent of cellular drug partitioning in blood observed in a dynamic in vivo system may be very different from the static in vitro values when repartitioning from blood cells is far slower than clearance of drug in plasma.

Apical ABC transporters and cancer chemotherapeutic drug disposition

ATP-binding cassette (ABC) transporters are transmembrane efflux transporters that mediate cellular extrusion of a broad range of substrates ranging from amino acids, lipids, and ions to xenobiotics including many anticancer drugs. ABCB1 (P-GP) and ABCG2 (BCRP) are the most extensively studied apical ABC drug efflux transporters. They are highly expressed in apical membranes of many pharmacokinetically relevant tissues such as epithelial cells of the small intestine and endothelial cells of the blood capillaries in brain and testis, and in the placental maternal-fetal barrier. In these tissues, they have a protective function as they efflux their substrates back to the intestinal lumen or blood and thus restrict the intestinal uptake and tissue disposition of many compounds. This presents a major challenge for the use of many (anticancer) drugs, as most currently used anticancer drugs are substrates of these transporters. Herein, we review the latest findings on the role of apical ABC transporters in the disposition of anticancer drugs. We discuss that many new, rationally designed anticancer drugs are substrates of these transporters and that their oral availability and/or brain disposition are affected by this interaction. We also summarize studies that investigate the improvement of oral availability and brain disposition of many cytotoxic (e.g., taxanes) and rationally designed (e.g., tyrosine kinase inhibitor) anticancer drugs, using chemical inhibitors of these transporters. These findings provide a better understanding of the importance of apical ABC transporters in chemotherapy and may therefore advance translation of promising preclinical insights and approaches to clinical studies.

Increased permeability-glycoprotein inhibition at the human blood-brain barrier can be safely achieved by performing PET during peak plasma concentrations of tariquidar

The permeability-glycoprotein (P-gp) efflux transporter is densely expressed at the blood-brain barrier, and its resultant spare capacity requires substantial blockade to increase the uptake of avid substrates, blunting the ability of investigators to measure clinically meaningful alterations in P-gp function. This study, conducted in humans, examined 2 P-gp inhibitors (tariquidar, a known inhibitor, and disulfiram, a putative inhibitor) and 2 routes of administration (intravenous and oral) to maximally increase brain uptake of the avid and selective P-gp substrate (11)C-N-desmethyl-loperamide (dLop) while avoiding side effects associated with high doses of tariquidar.

METHODS

Forty-two (11)C-dLop PET scans were obtained from 37 healthy volunteers. PET was performed with (11)C-dLop under the following 5 conditions: injected under baseline conditions without P-gp inhibition, injected 1 h after intravenous tariquidar infusion, injected during intravenous tariquidar infusion, injected after oral tariquidar, and injected after disulfiram. (11)C-dLop uptake was quantified with kinetic modeling using metabolite-corrected arterial input function or by measuring the area under the time-activity curve in the brain from 10 to 30 min.

RESULTS

Neither oral tariquidar nor oral disulfiram increased brain uptake of (11)C-dLop. Injecting (11)C-dLop during tariquidar infusion, when plasma tariquidar concentrations reach their peak, resulted in a brain uptake of the radioligand approximately 5-fold greater than baseline. Brain uptake was similar with 2 and 4 mg of intravenous tariquidar per kilogram; however, the lower dose was better tolerated. Injecting (11)C-dLop after tariquidar infusion also increased brain uptake, though higher doses (up to 6 mg/kg) were required. Brain uptake of (11)C-dLop increased fairly linearly with increasing plasma tariquidar concentrations, but we are uncertain whether maximal uptake was achieved.

CONCLUSION

We sought to increase the dynamic range of P-gp function measured after blockade. Performing (11)C-dLop PET during peak plasma concentrations of tariquidar, achieved with concurrent administration of intravenous tariquidar, resulted in greater P-gp inhibition at the human blood-brain barrier than delayed administration and allowed the use of a lower, more tolerable dose of tariquidar. On the basis of prior monkey studies, we suspect that plasma concentrations of tariquidar did not fully block P-gp; however, higher doses of tariquidar would likely be associated with unacceptable side effects.

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.

Challenges of using in vitro data for modeling P-glycoprotein efflux in the blood-brain barrier

The efficacy of central nervous system (CNS) drugs may be limited by their poor ability to cross the bloodbrain barrier (BBB). Transporters, such as p-glycoprotein, may affect the distribution of many drugs into the CNS in conjunction with the restricted paracellular pathway of the BBB. It is therefore important to gain information on unbound drug concentrations in the brain in drug development to ensure sufficient drug exposure from plasma at the target site in the CNS. In vitro methods are routinely used in drug development to study passive permeability and p-glycoprotein efflux of new drugs. This review discusses the challenges in the use of in vitro data as input parameters in physiologically based pharmacokinetic (PBPK) models of CNS drug disposition of p-glycoprotein substrates. Experience with quinidine demonstrates the variability in in vitro parameters of passive permeability and active pglycoprotein efflux. Further work is needed to generate parameter values that are independent of the model and assay. This is a prerequisite for reliable predictions of drug concentrations in the brain in vivo.

Compartmental models for apical efflux by P-glycoprotein: part 2--a theoretical study on transporter kinetic parameters

PURPOSE

The impact of efflux transporters in intracellular concentrations of a drug can be predicted with modeling techniques. In Part 1, several compartmental models were developed and evaluated. The goal of Part 2 was to apply these models to the characterization and interpretation of saturation kinetic data.

METHODS

The compartmental models from Part 1 were used to evaluate a previously published dataset from cell lines expressing varying levels of P-glycoprotein. Kinetic parameters for the transporter were estimated and compared across models.

RESULTS

Fits and errors for all compartmental models were identical. All compartmental models predicted more consistent parameters than the Michaelis-Menten model. The 5-compartment model with efflux out of the membrane predicted differential impact of P-gp upon apical versus basolateral drug exposure. Finally, the saturable kinetics of active efflux along with a permeability barrier was modeled to delineate a relationship between intracellular concentration with or without active efflux versus donor concentration. This relationship was not a rectangular hyperbola, but instead was shown to be a quadratic function.

CONCLUSIONS

One approach to estimate an in vivo transporter effect is to first model an intracellular Km value from in vitro data, and use this value along with the appropriate tissue transporter expression levels and relative surface area to calculate the relevant apparent Km (or Ki) values. Together with the results from Part 1, these studies suggest that compartmental models can provide a path forward to better utilize in vitro transporter data for in vivo predictions such as physiologically based pharmacokinetic modeling.

Compartmental models for apical efflux by P-glycoprotein--part 1: evaluation of model complexity

PURPOSE

With the goal of quantifying P-gp transport kinetics, Part 1 of these manuscripts evaluates different compartmental models and Part 2 applies these models to kinetic data.

METHODS

Models were developed to simulate the effect of apical efflux transporters on intracellular concentrations of six drugs. The effect of experimental variability on model predictions was evaluated. Several models were evaluated, and characteristics including membrane configuration, lipid content, and apical surface area (asa) were varied.

RESULTS

Passive permeabilities from MDCK-MDR1 cells in the presence of cyclosporine gave lower model errors than from MDCK control cells. Consistent with the results in Part 2, model configuration had little impact on calculated model errors. The 5-compartment model was the simplest model that reproduced experimental lag times. Lipid content and asa had minimal effect on model errors, predicted lag times, and intracellular concentrations. Including endogenous basolateral uptake activity can decrease model errors. Models with and without explicit membrane barriers differed markedly in their predicted intracellular concentrations for basolateral drug exposure. Single point data resulted in clearances similar to time course data.

CONCLUSIONS

Compartmental models are useful to evaluate the impact of efflux transporters on intracellular concentrations. Whereas a 3-compartment model may be sufficient to predict the impact of transporters that efflux drugs from the cell, a 5-compartment model with explicit membranes may be required to predict intracellular concentrations when efflux occurs from the membrane. More complex models including additional compartments may be unnecessary.

Computational approaches to analyse and predict small molecule transport and distribution at cellular and subcellular levels

Quantitative structure-activity relationship (QSAR) studies and mechanistic mathematical modeling approaches have been independently employed for analysing and predicting the transport and distribution of small molecule chemical agents in living organisms. Both of these computational approaches have been useful for interpreting experiments measuring the transport properties of small molecule chemical agents, in vitro and in vivo. Nevertheless, mechanistic cell-based pharmacokinetic models have been especially useful to guide the design of experiments probing the molecular pathways underlying small molecule transport phenomena. Unlike QSAR models, mechanistic models can be integrated from microscopic to macroscopic levels, to analyse the spatiotemporal dynamics of small molecule chemical agents from intracellular organelles to whole organs, well beyond the experiments and training data sets upon which the models are based. Based on differential equations, mechanistic models can also be integrated with other differential equations-based systems biology models of biochemical networks or signaling pathways. Although the origin and evolution of mathematical modeling approaches aimed at predicting drug transport and distribution has occurred independently from systems biology, we propose that the incorporation of mechanistic cell-based computational models of drug transport and distribution into a systems biology modeling framework is a logical next step for the advancement of systems pharmacology research.

PBPK models for the prediction of in vivo performance of oral dosage forms

Drug absorption from the gastrointestinal (GI) tract is a highly complex process dependent upon numerous factors including the physicochemical properties of the drug, characteristics of the formulation and interplay with the underlying physiological properties of the GI tract. The ability to accurately predict oral drug absorption during drug product development is becoming more relevant given the current challenges facing the pharmaceutical industry. Physiologically-based pharmacokinetic (PBPK) modeling provides an approach that enables the plasma concentration-time profiles to be predicted from preclinical in vitro and in vivo data and can thus provide a valuable resource to support decisions at various stages of the drug development process. Whilst there have been quite a few successes with PBPK models identifying key issues in the development of new drugs in vivo, there are still many aspects that need to be addressed in order to maximize the utility of the PBPK models to predict drug absorption, including improving our understanding of conditions in the lower small intestine and colon, taking the influence of disease on GI physiology into account and further exploring the reasons behind population variability. Importantly, there is also a need to create more appropriate in vitro models for testing dosage form performance and to streamline data input from these into the PBPK models. As part of the Oral Biopharmaceutical Tools (OrBiTo) project, this review provides a summary of the current status of PBPK models available. The current challenges in PBPK set-ups for oral drug absorption including the composition of GI luminal contents, transit and hydrodynamics, permeability and intestinal wall metabolism are discussed in detail. Further, the challenges regarding the appropriate integration of results from in vitro models, such as consideration of appropriate integration/estimation of solubility and the complexity of the in vitro release and precipitation data, are also highlighted as important steps to advancing the application of PBPK models in drug development. It is expected that the "innovative" integration of in vitro data from more appropriate in vitro models and the enhancement of the GI physiology component of PBPK models, arising from the OrBiTo project, will lead to a significant enhancement in the ability of PBPK models to successfully predict oral drug absorption and advance their role in preclinical and clinical development, as well as for regulatory applications.

Impact of P-glycoprotein (ABCB1) and breast cancer resistance protein (ABCG2) gene dosage on plasma pharmacokinetics and brain accumulation of dasatinib, sorafenib, and sunitinib

Low brain accumulation of anticancer drugs due to efflux transporters may limit chemotherapeutic efficacy, necessitating a better understanding of the underlying mechanisms. P-glycoprotein (Abcb1a/1b) and breast cancer resistance protein (Abcg2) combination knockout mice often display disproportionately increased brain accumulation of shared drug substrates compared with single transporter knockout mice. Recently developed pharmacokinetic models could explain this phenomenon. To experimentally test these models and their wider relevance for tyrosine kinase inhibitors and other drugs, we selected dasatinib, sorafenib, and sunitinib because of their divergent oral availability and brain accumulation profiles: the brain accumulation of dasatinib is mainly restricted by Abcb1, that of sorafenib mainly by Abcg2, and that of sunitinib equally by Abcb1 and Abcg2. We analyzed the effect of halving the efflux activity of these transporters at the blood-brain barrier by generating heterozygous Abcb1a/1b;Abcg2 knockout mice and testing the plasma and brain levels of the drugs after oral administration at 10 mg/kg. Real-time reverse transcription-polymerase chain reaction analysis confirmed the ∼2-fold decreased expression of both transporters in brain. Interestingly, whereas complete knockout of the transporters caused 24- to 36-fold increases in brain accumulation of the drugs, the heterozygous mice only displayed 1.6- to 1.9-fold increases of brain accumulation relative to wild-type mice. These results are well in line with the predictions of the pharmacokinetic models and provide strong support for their validity for a wider range of drugs. Moreover, retrospective analysis of fetal accumulation of drugs across the placenta in Abcb1a/1b heterozygous knockout pups suggests that these models equally apply to the maternal-fetal barrier.

Application of permeability-limited physiologically-based pharmacokinetic models: part I-digoxin pharmacokinetics incorporating P-glycoprotein-mediated efflux

A prerequisite for the prediction of the magnitude of P-glycoprotein (P-gp)-mediated drug-drug interactions between digoxin and P-gp inhibitors (e.g. verapamil and its metabolite norverapamil) or P-gp inducers (e.g. rifampicin) is a predictive pharmacokinetic model for digoxin itself. Thus, relevant in vitro metabolic, transporter and inhibitory data incorporated into permeability-limited models, such as the "advanced dissolution, absorption and metabolism" (ADAM) module and the permeability-limited liver (PerL) module, integrated with a mechanistic physiologically-based pharmacokinetic (PBPK) model such as that of the Simcyp Simulator (version 12.2) are necessary. Simulated concentration-time profiles of digoxin generated using the developed model were consistent with observed data across 31 independent studies [13 intravenous single dose (SD), 12 per oral SD and six multiple dose studies]. The fact that predicted tmax (time of maximum plasma concentration observed) and Cmax (maximum plasma concentration observed) of oral digoxin were similar to observed values indicated that the relative contributions of permeation and P-gp-mediated efflux in the model were appropriate. There was no indication of departure from dose proportionality over the dose range studied (0.25-1.5 mg). All dose normalised area under the plasma concentration-time curve profiles (AUCs) for the 0.25, 0.5, 0.75 and 1 mg doses resembled each other. Thus, PBPK modelling in conjunction with mechanistic absorption and distribution models and reliable in vitro transporter data can be used to assess the impact of dose on P-gp-mediated efflux (or otherwise).

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

A P-glycoprotein (P-gp) IC₅₀ working group was established with 23 participating pharmaceutical and contract research laboratories and one academic institution to assess interlaboratory variability in P-gp IC₅₀ determinations. Each laboratory followed its in-house protocol to determine in vitro IC₅₀ values for 16 inhibitors using four different test systems: human colon adenocarcinoma cells (Caco-2; eleven laboratories), Madin-Darby canine kidney cells transfected with MDR1 cDNA (MDCKII-MDR1; six laboratories), and Lilly Laboratories Cells--Porcine Kidney Nr. 1 cells transfected with MDR1 cDNA (LLC-PK1-MDR1; four laboratories), and membrane vesicles containing human P-glycoprotein (P-gp; five laboratories). For cell models, various equations to calculate remaining transport activity (e.g., efflux ratio, unidirectional flux, net-secretory-flux) were also evaluated. The difference in IC₅₀ values for each of the inhibitors across all test systems and equations ranged from a minimum of 20- and 24-fold between lowest and highest IC₅₀ values for sertraline and isradipine, to a maximum of 407- and 796-fold for telmisartan and verapamil, respectively. For telmisartan and verapamil, variability was greatly influenced by data from one laboratory in each case. Excluding these two data sets brings the range in IC₅₀ values for telmisartan and verapamil down to 69- and 159-fold. The efflux ratio-based equation generally resulted in severalfold lower IC₅₀ values compared with unidirectional or net-secretory-flux equations. Statistical analysis indicated that variability in IC₅₀ values was mainly due to interlaboratory variability, rather than an implicit systematic difference between test systems. Potential reasons for variability are discussed and the simplest, most robust experimental design for P-gp IC₅₀ determination proposed. The impact of these findings on drug-drug interaction risk assessment is discussed in the companion article (Ellens et al., 2013) and recommendations are provided.