Case Study 8: Status of the Structural Mass Action Kinetic Model of P-gp-Mediated Transport Through Confluent Cell Monolayers.

In the first edition of this book, we presented the basics of explicitly incorporating the lipid biochemistry into a confluent cell monolayer transport model and the novel findings of this model up to 2013, including the use of global optimization to fit the elementary rate constants and the efflux active P-glycoprotein (P-gp) membrane concentrations for the transport of four P-gp substrates across MDCKII-hMDR1-NKI confluent cell monolayers. This chapter is an update on that model, which has been focused primarily on discovering how microvilli morphology regulates the efflux active P-gp and the existence of, as yet, unidentified uptake transporters of P-gp substrates in all of the commonly used P-gp expressing cell lines used in the pharmaceutical industry, thereby adding new players to DDI predictions and IVIVE. The structural mass action kinetic model uses the general mass action reactions for P-gp binding and efflux, with the membrane structural parameters for the confluent cell monolayer to predict drug transport over time.

Mechanistic kinetic modeling generates system-independent P-glycoprotein mediated transport elementary rate constants for inhibition and, in combination with 3D SIM microscopy, elucidates the importance of microvilli morphology on P-glycoprotein mediated efflux activity.

In vitro transporter kinetics are typically analyzed by steady-state Michaelis-Menten approximations. However, no clear evidence exists that these approximations, applied to multiple transporters in biological membranes, yield system-independent mechanistic parameters needed for reliable in vivo hypothesis generation and testing. Areas covered: The classical mass action model has been developed for P-glycoprotein (P-gp) mediated transport across confluent polarized cell monolayers.

Hepatobiliary Disposition of Atovaquone: A Case of Mechanistically Unusual Biliary Clearance.

Atovaquone, an antiprotozoal and antipneumocystic agent, is predominantly cleared by biliary excretion of unchanged parent drug. Atovaquone is ≥10,000-fold concentrated in human bile relative to unbound plasma. Even after correcting for apparent nonspecific binding and incomplete solubility in bile, atovaquone is still concentrated ≥100-fold in bile, consistent with active biliary excretion.

Derivation of a System-Independent for P-glycoprotein Mediated Digoxin Transport from System-Dependent IC Data.

It has been previously demonstrated that IC values for inhibition of digoxin transport across confluent polarized cell monolayers are system-dependent. Digoxin IC data from five laboratories participating in the P-glycoprotein (P-gp) IC Initiative, using Caco-2, MDCKII-hMDR1 or LLC-PK1-hMDR1 cells, were fitted by the structural mass action kinetic model for P-gp-mediated transport across confluent cell monolayers. We determined their efflux-active P-gp concentration [T(0)], inhibitor elementary dissociation rate constant from P-gp (), digoxin basolateral uptake clearance (), and inhibitor binding affinity to the digoxin basolateral uptake transporter ().

Prediction of the Transporter-Mediated Drug-Drug Interaction Potential of Dabrafenib and Its Major Circulating Metabolites.

The BRAF inhibitor dabrafenib was recently approved for the treatment of certain BRAF V600 mutation-positive tumors, either alone or in combination therapy with the mitogen-activated extracellular signal regulated kinase 1 (MEK1) and MEK2 inhibitor, trametinib. This article presents the dabrafenib transporter-mediated drug-drug interaction (DDI) risk assessment, which is currently an important part of drug development, regulatory submission, and drug registration. Dabrafenib and its major circulating metabolites (hydroxy-, carboxy-, and desmethyl-dabrafenib) were investigated as inhibitors of the clinically relevant transporters P-gp, BCRP, OATP1B1, OATP1B3, OCT2, OAT1, and OAT3.

Extrapolation of Elementary Rate Constants of P-glycoprotein-Mediated Transport from MDCKII-hMDR1-NKI to Caco-2 Cells.

The best parameters for incorporation into mechanistic physiologically based pharmacokinetic models for transporters are system-independent kinetic parameters and active (not total) transporter levels. Previously, we determined the elementary rate constants for P-glycoprotein (P-gp)-mediated transport (on- and off-rate constants from membrane to P-gp binding pocket and efflux rate constant into the apical chamber) using the structural mass action kinetic model in confluent MDCKII-hMDR1-NKI cell monolayers. In the present work, we extended the kinetic analysis to Caco-2 cells for the first time and showed that the elementary rate constants are very similar compared with MDCKII-hMDR1-NKI cells, suggesting they primarily depend on the interaction of the compound with P-gp and are therefore mostly independent of the in vitro system used.

Microvilli Morphology Can Affect Efflux Active P-Glycoprotein in Confluent MDCKII -hMDR1-NKI and Caco-2 Cell Monolayers.

From fits of drug transport kinetics across confluent MDCKII-hMDR1-NKI and Caco-2 cell monolayers we estimated the levels of efflux active P-glycoprotein (P-gp) in these two cell lines (companion paper). In the present work, we compared the efflux active P-gp number to the total P-gp level, using liquid chromatography-tandem mass spectrometry, and showed that in Caco-2 cells total P-gp is about 10-fold greater than efflux active P-gp, whereas in MDCKII-hMDR1-NKI cells these values are within twofold. We further visualized the microvilli in MDCKII-hMDR1-NKI and Caco-2 cells using three-dimensional structured illumination super-resolution microscopy and found that the microvilli in Caco-2 cells are taller and more densely packed than those in MDCK-hMDR1-NKI cells.

Inhibition of Intestinal OATP2B1 by the Calcium Receptor Antagonist Ronacaleret Results in a Significant Drug-Drug Interaction by Causing a 2-Fold Decrease in Exposure of Rosuvastatin.

Rosuvastatin is a widely prescribed antihyperlipidemic which undergoes limited metabolism, but is an in vitro substrate of multiple transporters [organic anion transporting polypeptide 1B1 (OATP1B1), OATP1B3, OATP1A2, OATP2B1, sodium-taurocholate cotransporting polypeptide, breast cancer resistance protein (BCRP), multidrug resistance protein 2 (MRP2), MRP4, organic anion transporter 3]. It is therefore frequently used as a probe substrate in clinical drug-drug interaction (DDI) studies to investigate transporter inhibition. Although each of these transporters is believed to play a role in rosuvastatin disposition, multiple pharmacogenetic studies confirm that OATP1B1 and BCRP play an important role in vivo.

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 IC50 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 IC50 values, which can be easily adopted across the industry. The current work provides a quantitative examination of critical factors affecting the quality of IC50 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 IC50 initiative.

Breast cancer resistance protein (ABCG2) in clinical pharmacokinetics and drug interactions: practical recommendations for clinical victim and perpetrator drug-drug interaction study design.

Breast cancer resistance protein (BCRP; ABCG2) limits intestinal absorption of low-permeability substrate drugs and mediates biliary excretion of drugs and metabolites. Based on clinical evidence of BCRP-mediated drug-drug interactions (DDIs) and the c.421C>A functional polymorphism affecting drug efficacy and safety, both the US Food and Drug Administration and European Medicines Agency recommend preclinical evaluation and, when appropriate, clinical assessment of BCRP-mediated DDIs.

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).

Disposition and metabolism of darapladib, a lipoprotein-associated phospholipase A2 inhibitor, in humans.

The absorption, metabolism, and excretion of darapladib, a novel inhibitor of lipoprotein-associated phospholipase A2, was investigated in healthy male subjects using [(14)C]-radiolabeled material in a bespoke study design. Disposition of darapladib was compared following single i.v.

The importance of villous physiology and morphology in mechanistic physiologically-based pharmacokinetic models.

Purpose: Existing PBPK models incorporating intestinal first-pass metabolism account for effect of drug permeability on accessible absorption surface area by use of "effective" permeability, P eff , without adjusting number of enterocytes involved in absorption or proportion of intestinal CYP3A involved in metabolism. The current model expands on existing models by accounting for these factors.

Methods: The PBPK model was developed using SAAM II.

Transport inhibition of digoxin using several common P-gp expressing cell lines is not necessarily reporting only on inhibitor binding to P-gp.

We have reported that the P-gp substrate digoxin required basolateral and apical uptake transport in excess of that allowed by digoxin passive permeability (as measured in the presence of GF120918) to achieve the observed efflux kinetics across MDCK-MDR1-NKI (The Netherlands Cancer Institute) confluent cell monolayers. That is, GF120918 inhibitable uptake transport was kinetically required. Therefore, IC50 measurements using digoxin as a probe substrate in this cell line could be due to inhibition of P-gp, of digoxin uptake transport, or both.

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.

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).

ITC recommendations for transporter kinetic parameter estimation and translational modeling of transport-mediated PK and DDIs in humans.

This white paper provides a critical analysis of methods for estimating transporter kinetics and recommendations on proper parameter calculation in various experimental systems. Rational interpretation of transporter-knockout animal findings and application of static and dynamic physiologically based modeling approaches for prediction of human transporter-mediated pharmacokinetics and drug-drug interactions (DDIs) are presented. The objective is to provide appropriate guidance for the use of in vitro, in vivo, and modeling tools in translational transporter science.

Optimizing the in vitro and clinical assessment of drug interaction risk by understanding co-medications in patient populations.

Introduction: Pharmacokinetic drug interactions resulting from the inhibition of drug elimination mechanisms are of concern in drug development due to the clinical risk associated with elevated drug concentrations. Advances in understanding the mechanistic basis of these drug interactions has resulted in the widespread application of mechanistic in vitro assays and the conduct of clinical drug interaction studies to predict and quantify the risks in drug development.

Areas Covered: The authors investigate co-medication prescriptions in target patient populations and characterize the mechanistic basis and clinical impact of known co-medication drug interactions.

Preclinical strategy to reduce clinical hepatotoxicity using in vitro bioactivation data for >200 compounds.

Drug-induced liver injury is the most common cause of market withdrawal of pharmaceuticals, and thus, there is considerable need for better prediction models for DILI early in drug discovery. We present a study involving 223 marketed drugs (51% associated with clinical hepatotoxicity; 49% non-hepatotoxic) to assess the concordance of in vitro bioactivation data with clinical hepatotoxicity and have used these data to develop a decision tree to help reduce late-stage candidate attrition. Data to assess P450 metabolism-dependent inhibition (MDI) for all common drug-metabolizing P450 enzymes were generated for 179 of these compounds, GSH adduct data generated for 190 compounds, covalent binding data obtained for 53 compounds, and clinical dose data obtained for all compounds.

Fitting the elementary rate constants of the P-gp transporter network in the hMDR1-MDCK confluent cell monolayer using a particle swarm algorithm.

P-glycoprotein, a human multidrug resistance transporter, has been extensively studied due to its importance to human health and disease. In order to understand transport kinetics via P-gp, confluent cell monolayers overexpressing P-gp are widely used. The purpose of this study is to obtain the mass action elementary rate constants for P-gp's transport and to functionally characterize members of P-gp's network, i.

Digoxin is not a substrate for organic anion-transporting polypeptide transporters OATP1A2, OATP1B1, OATP1B3, and OATP2B1 but is a substrate for a sodium-dependent transporter expressed in HEK293 cells.

Digoxin, an orally administered cardiac glycoside cardiovascular drug, has a narrow therapeutic window. Circulating digoxin levels (maximal concentration of ∼1.5 ng/ml) require careful monitoring, and the potential for drug-drug interactions (DDI) is a concern.

A mechanism-based mathematical model of aryl hydrocarbon receptor-mediated CYP1A induction in rats using beta-naphthoflavone as a tool compound.

β-Naphthoflavone (BNF) is a synthetic flavone that selectively and potently induces CYP1A enzymes via aryl hydrocarbon receptor activation. Mechanism-based mathematical models of CYP1A enzyme induction were developed to predict the time course of enzyme induction and quantitatively evaluate the interrelationship between BNF plasma concentrations, hepatic CYP1A1 and CYP1A2 mRNA levels, and CYP1A enzyme activity in rats in vivo. Male Sprague-Dawley rats received a continuous intravenous infusion of vehicle or 1.

Pazopanib-induced hyperbilirubinemia is associated with Gilbert's syndrome UGT1A1 polymorphism.

Background: Pazopanib has shown clinical activity against multiple tumour types and is generally well tolerated. However, isolated elevations in transaminases and bilirubin have been observed. This study examined polymorphisms in molecules involved in pharmacokinetic and pharmacodynamic pathways of pazopanib and their association with hepatic dysfunction.

If the KI is defined by the free energy of binding to P-glycoprotein, which kinetic parameters define the IC50 for the Madin-Darby canine kidney II cell line overexpressing human multidrug resistance 1 confluent cell monolayer?

From previous fits of drug transport kinetics across confluent Madin-Darby canine kidney II cell line overexpressing human multidrug resistance 1 cell monolayers, we found that a drug's binding constant to P-glycoprotein (P-gp) was significantly smaller than its IC(50) when that drug was used as an inhibitor against another P-gp substrate. We tested several IC(50) candidate functions, including the standard function, the Kalvass-Pollack function, and the efflux ratio, to determine whether any of them yielded an IC(50) = K(I), as would be expected for water-soluble enzymes. For the confluent cell monolayer, the IC(50)/K(I) ratio is greater than 1 for all candidate functions tested.

Biopharmaceutics classification system: validation and learnings of an in vitro permeability assay.

The Biopharmaceutics Classification System (BCS) is the scientific basis for classifying drugs based on their aqueous solubility and intestinal permeability that supports in vivo bioavailability and bioequivalence waivers for immediate-release solid dosage form drugs. One requirement of the BCS is that the permeability method must be validated. In order to accommodate the variety of in vitro/in situ permeability models, the BCS Guidance gives a general framework for the validation requirements, necessitating implemented experimental details to be selected by the applicant laboratory.