Data-Driven Ensemble Docking to Map Molecular Interactions of Steroid Analogs with Hepatic Organic Anion Transporting Polypeptides

Hepatic organic anion transporting polypeptides—OATP1B1, OATP1B3, and OATP2B1—are expressed at the basolateral membrane of hepatocytes, being responsible for the uptake of a wide range of natural substrates and structurally unrelated pharmaceuticals. Impaired function of hepatic OATPs has been linked to clinically relevant drug–drug interactions leading to altered pharmacokinetics of administered drugs. Therefore, understanding the commonalities and differences across the three transporters represents useful knowledge to guide the drug discovery process at an early stage. Unfortunately, such efforts remain challenging because of the lack of experimentally resolved protein structures for any member of the OATP family. In this study, we established a rigorous computational protocol to generate and validate structural models for hepatic OATPs. The multistep procedure is based on the systematic exploration of available protein structures with shared protein folding using normal-mode analysis, the calculation of multiple template backbones from elastic network models, the utilization of multiple template conformations to generate OATP structural models with various degrees of conformational flexibility, and the prioritization of models on the basis of enrichment docking. We employed the resulting OATP models of OATP1B1, OATP1B3, and OATP2B1 to elucidate binding modes of steroid analogs in the three transporters. Steroid conjugates have been recognized as endogenous substrates of these transporters. Thus, investigating this data set delivers insights into mechanisms of substrate recognition. In silico predictions were complemented with in vitro studies measuring the bioactivity of a compound set on OATP expressing cell lines. Important structural determinants conferring shared and distinct binding patterns of steroid analogs in the three transporters have been identified. Overall, this comparative study provides novel insights into hepatic OATP-ligand interactions and selectivity. Furthermore, the integrative computational workflow for structure-based modeling can be leveraged for other pharmaceutical targets of interest.

Current Advances in Studying Clinically Relevant Transporters of the Solute Carrier (SLC) Family by Connecting Computational Modeling and Data Science

Organic anion and cation transporting proteins (OATs, OATPs, and OCTs), as well as the Multidrug and Toxin Extrusion (MATE) transporters of the Solute Carrier (SLC) family are playing a pivotal role in the discovery and development of new drugs due to their involvement in drug disposition, drug-drug interactions, adverse drug effects and related toxicity. Computational methods to understand and predict clinically relevant transporter interactions can provide useful guidance at early stages in drug discovery and design, especially if they include contemporary data science approaches. In this review, we summarize the current state-of-the-art of computational approaches for exploring ligand interactions and selectivity for these drug (uptake) transporters. The computational methods discussed here by highlighting interesting examples from the current literature are ranging from semiautomatic data mining and integration, to ligand-based methods (such as quantitative structure-activity relationships, and combinatorial pharmacophore modeling), and finally structure-based methods (such as comparative modeling, molecular docking, and molecular dynamics simulations). We are focusing on promising computational techniques such as fold-recognition methods, proteochemometric modeling or techniques for enhanced sampling of protein conformations used in the context of these ADMET-relevant SLC transporters with a special focus on methods useful for studying ligand selectivity.

Prediction and Quantification of Hepatic Transporter-Mediated Uptake of Pitavastatin Utilizing a Combination of the Relative Activity Factor Approach and Mechanistic Modeling

Quantification of the fraction transported (f_t) by a particular transporter will facilitate more robust estimations of transporter interactions. Using pitavastatin as a model uptake transporter substrate, we investigated the utility of the relative activity factor (RAF) approach and mechanistic modeling to estimate f_t in hepatocytes. The transporters evaluated were organic anion-transporting polypeptides OATP1B1 and OATP1B3 and sodium-taurocholate cotransporting polypeptide. Transporter-expressing human embryonic kidney 293 cells and human hepatocytes were used for determining RAF values, which were then incorporated into the mechanistic model to simulate hepatocyte uptake of pitavastatin over time. There was excellent agreement between simulated and observed hepatocyte uptake of pitavastatin, indicating the suitability of this approach for translation of uptake from individual transporter-expressing cells to more holistic in vitro models. Subsequently, f_t values were determined. The largest contributor to hepatocyte uptake of pitavastatin was OATP1B1, which correlates with what is known about the in vivo disposition of pitavastatin. The f_t values were then used for evaluating in vitro-in vivo correlations of hepatic uptake inhibition with OATP inhibitors rifampicin and cyclosporine. Predictions were compared with previously reported plasma exposure changes of pitavastatin with these inhibitors. Although hepatic uptake inhibition of pitavastatin was 2-3-fold underpredicted, incorporation of scaling factors (SFs) into RAF values significantly improved the predictive ability. We propose that calibration of hepatocytes with standard transporter substrates and inhibitors would allow for determination of system-specific SFs, which could subsequently be used for refining predictions of clinical DDI potential for new chemical entities that undergo active hepatic uptake.

Cell cultures in drug discovery and development: The need of reliable in vitro-in vivo extrapolation for pharmacodynamics and pharmacokinetics assessment

For ethical and cost-related reasons, use of animals for the assessment of mode of action, metabolism and/or toxicity of new drug candidates has been increasingly scrutinized in research and industrial applications. Implementation of the 3 "Rs"¹; rule (Reduction, Replacement, Refinement) through development of in silico or in vitro assays has become an essential element of risk assessment. Physiologically based pharmacokinetic (PBPK²) modeling is the most potent in silico tool used for extrapolation of pharmacokinetic parameters to animal or human models from results obtained in vitro. Although, many types of in vitro assays are conducted during drug development, use of cell cultures is the most reliable one. Two-dimensional (2D) cell cultures have been a part of drug development for many years. Nowadays, their role is decreasing in favor of three-dimensional (3D) cell cultures and co-cultures. 3D cultures exhibit protein expression patterns and intercellular junctions that are closer to in vivo states in comparison to classical monolayer cultures. Co-cultures allow for examinations of the mutual influence of different cell lines. However, the complexity and high costs of co-cultures and 3D equipment exclude such methods from high-throughput screening (HTS).³In vitro absorption, distribution, metabolism, and excretion assessment, as well as drug-drug interaction (DDI), are usually performed with the use of various cell culture based assays. Progress in in silico and in vitro methods can lead to better in vitro-in vivo extrapolation (IVIVE⁴) outcomes and have a potential to contribute towards a significant reduction in the number of laboratory animals needed for drug research. As such, concentrated efforts need to be spent towards the development of an HTS in vitro platform with satisfactory IVIVE features.

Evaluation of transporters in drug development: Current status and contemporary issues

Transporters govern the access of molecules to cells or their exit from cells, thereby controlling the overall distribution of drugs to their intracellular site of action. Clinically relevant drug-drug interactions mediated by transporters are of increasing interest in drug development. Drug transporters, acting alone or in concert with drug metabolizing enzymes, can play an important role in modulating drug absorption, distribution, metabolism and excretion, thus affecting the pharmacokinetics and/or pharmacodynamics of a drug. The drug interaction guidance documents from regulatory agencies include various decision criteria that may be used to predict the need for in vivo assessment of transporter-mediated drug-drug interactions. Regulatory science research continues to assess the prediction performances of various criteria as well as to examine the strength and limitations of each prediction criterion to foster discussions related to harmonized decision criteria that may be used to facilitate global drug development. This review discusses the role of transporters in drug development with a focus on methodologies in assessing transporter-mediated drug-drug interactions, challenges in both in vitro and in vivo assessments of transporters, and emerging transporter research areas including biomarkers, assessment of tissue concentrations, and effect of diseases on transporters.

An S-warfarin and AZD1981 interaction: in vitro and clinical pilot data suggest the N-deacetylated amino acid metabolite as the primary perpetrator

AIM

AZD1981 is an orally bioavailable chemoattractant receptor-homologous molecule expressed on Th2 cells (CRTh2) receptor antagonist progressed to phase II trials for the treatment of allergic asthma. Previously performed in vitro human hepatocyte incubations identified N-deacetylated AZD1981 as a primary metabolite. We report on metabolite exposure from a clinical excretion balance, on in vitro studies performed to determine the likelihood of a metabolite-dependent drug-drug interaction (DDI) and on a clinical warfarin DDI study. The aim was to demonstrate that N-deacetylated AZD1981 is responsible for the observed interaction.

METHODS

The excretion and biotransformation of [¹⁴ C]-AZD1981 were studied in healthy male volunteers, and subsequently in vitro cytochrome P450 (CYP) inhibition and hepatocyte uptake investigations were carried out with metabolites and the parent drug. A clinical DDI study using coadministered twice-daily 100 mg and 400 mg AZD1981 with 25 mg warfarin was performed.

RESULTS

The excretion balance study showed N-deacetylated AZD1981 to be the most abundant metabolite present in plasma. In vitro data revealed the metabolite to be a weak CYP2C9 time-dependent inhibitor, subject to more active hepatic uptake than the parent molecule. Clinically, the S-warfarin area under the plasma concentration-time curve increased, on average, 1.4-fold [95% confidence interval (CI) 1.22, 1.50] and 2.4-fold (95% CI 2.11, 2.64) after 100 mg (n = 13) and 400 mg (n = 11) AZD1981 administration, respectively. In vitro CYP inhibition and hepatocyte uptake data were used to explain the interaction.

CONCLUSIONS

N-deacetylated AZD1981 can be added to the small list of drug metabolites reported as sole contributors to clinical drug-drug interactions, with weak time-dependent inhibition exacerbated by efficient hepatic uptake being the cause.

Development of a high-throughput mass spectrometry based analytical method to support an in vitro OATP1B1 inhibition screening assay

RATIONALE

It is well known that the organic anion transporting polypeptide 1B1 (OATP1B1) plays a major role in the hepatic uptake of a range of drugs. To this end, it is pivotal that the potential for new molecular entities (NMEs) to inhibit OATP1B1 activity be assessed during early drug discovery. The work reported herein describes the development of a high-throughput analytical method to measure the clinically relevant probe substrate, pitavastatin, for the in vitro assessment of OATP1B1 inhibition.

METHODS

Development of an analytical method capable of very fast throughput was crucial for the success of this assay and was accomplished using a system which combines direct, on-line solid-phase extraction (SPE) with highly sensitive, label-free tandem mass spectrometry (MS/MS)-based detection. Mass spectrometry analysis of pitavastatin, along with the stable isotopically labeled internal standard d5-pitavastatin, was conducted using positive electrospray ionization (ESI) in selected reaction monitoring (SRM) mode.

RESULTS

The on-line SPE-MS/MS platform demonstrated similar sensitivity, selectivity, reproducibility, linearity and robustness to existing methodologies while achieving analytical cycle times of 10.4 seconds per well. Sensitivity exceeded what was necessary for our assay conditions, with a determined lower limit of quantification (LLOQ) for pitavastatin of 10 pM (picomolar) in assay matrix. Furthermore, the potency of multiple reference compounds was shown to be within 2-fold of IC50 values generated from liquid chromatography (LC)/MS/MS-based literature values.

CONCLUSIONS

A very fast and robust analytical method was successfully developed for the measurement of the clinically relevant OATP1B1 substrate, pitavastatin. The successful development and implementation of this very important early liability screen has helped to facilitate judicious lead candidate progression and will ultimately help build a greater understanding of OATP1B1-NME interactions, in general. Copyright © 2016 John Wiley & Sons, Ltd.

Effect of OATP-binding on the prediction of biliary excretion

1. Biliary excretion of compounds is dependant on several transporter proteins for the active uptake of compounds from the blood into the hepatocytes. Organic anion-transporting polypeptides (OATPs) are some of the most abundant transporter proteins in the sinusoidal membrane and have been shown to have substrate specificity similar to the structural characteristics of cholephilic compounds. 2. In this study, we sought to use measures of OATP binding as predictors of biliary excretion in conjunction with molecular descriptors in a quantitative structure-activity relationship (QSAR) study. Percentage inhibitions of three subtypes of OATPs were used as surrogate indicators of OATP substrates. Several statistical modelling techniques were incorporated including classification and regression trees, boosted trees, random forest and multivariate adaptive regression splines (MARS) in order to first develop QSARs for the prediction of OATP inhibition of compounds. The predicted OATP percentage inhibition using selected models were then used as features of the QSAR models for the prediction of biliary excretion of compounds in rat. 3. The results indicated that incorporation of predicted OATP inhibition improves accuracy of biliary excretion models. The best result was obtained from a simple regression tree that used predicted OATP1B1 percentage inhibition at the root node of the tree.

Hepatic drug transporters: the journey so far

INTRODUCTION

The key role of transporter biology in both the manifestation and treatment of disease is now firmly established. Experiences of sub-optimal drug exposure due to drug-transporter interplay have supported incorporation of studies aimed at understanding the interactions between compounds and drug transporters much earlier in drug discovery. While drug transporters can impact the most pivotal pharmacokinetic parameter with respect to human dose and exposure projections, clearance, at a renal or hepatobiliary level, the latter will form the focus of this perspective.

AREAS COVERED

A synopsis of guidelines on which transporters to study together with an overview of the currently available toolkit is presented. A perspective on when to conduct studies with various hepatic transporters is also provided together with structural "alerts" which should prompt early investigation.

EXPERT OPINION

Great progress has been made in individual laboratories and via consortia to understand the role of drug transporters in disease, drug disposition, drug-drug interactions and toxicity. A systematic analysis of the value posed by the available approaches and an inter-lab comparison now seems warranted. The emerging ability to use physico-chemical properties to guide future screening cascades promises to revolutionise the efficiency of early drug discovery.

Montelukast Disposition: No Indication of Transporter-Mediated Uptake in OATP2B1 and OATP1B1 Expressing HEK293 Cells

Clinical studies with montelukast show variability in effect and polymorphic OATP2B1-dependent absorption has previously been implicated as a possible cause. This claim has been challenged with conflicting data and here we used OATP2B1-transfected HEK293 cells to clarify the mechanisms involved. For montelukast, no significant difference in cell uptake between HEK-OATP2B1 and empty vector cell lines was observed at pH 6.5 or pH 7.4, and no concentration-dependent uptake was detected. Montelukast is a carboxylic acid, a relatively potent inhibitor of OATP1B1, OATP1B3, and OATP2B1, and has previously been postulated to be actively transported into human hepatocytes. Using OATP1B1-transfected HEK293 cells and primary human hepatocytes in the presence of OATP inhibitors we demonstrate for the first time that active OATP-dependent transport is unlikely to play a significant role in the human disposition of montelukast.

Investigation of the impact of substrate selection on in vitro organic anion transporting polypeptide 1B1 inhibition profiles for the prediction of drug-drug interactions

The risk assessment of organic anion transporting polypeptide (OATP) 1B1-mediated drug-drug interactions (DDIs) is an indispensable part of drug development. We previously reported that in vitro inhibitory potencies of several inhibitors on OATP1B1 depend on the substrates when prototypical substrates, estradiol-17β-glucuronide (E₂G), estrone-3-sulfate, and sulfobromophthalein were used as test substrates. The purpose of this study was to comprehensively investigate this substrate-dependent inhibition of OATP1B1 using clinically relevant OATP1B1 inhibitors and substrate drugs. Effects of cyclosporine A (CsA), rifampin, and gemfibrozil on OATP1B1-mediated uptake of 12 substrate drugs were examined in OATP1B1-expressing human embryonic kidney 293 cells. The Ki values (μM) for CsA varied from 0.0771 to 0.486 (6.3-fold), for rifampin from 0.358 to 1.23 (3.4-fold), and for gemfibrozil from 9.65 to 252 (26-fold). Except for the inhibition of torasemide uptake by CsA and that of nateglinide uptake by gemfibrozil, the Ki values were within 2.8-fold of those obtained using E₂G as a substrate. Preincubation potentiated the inhibitory effect of CsA on OATP1B1 with similar magnitude regardless of the substrates. R values calculated based on a static model showed some variation depending on the Ki values determined with various substrates, and such variability could have an impact on the DDI predictions particularly for a weak-to-moderate inhibitor (gemfibrozil). OATP1B1 substrate drugs except for torasemide and nateglinide, or E₂G as a surrogate, is recommended as an in vitro probe in the inhibition experiments, which will help mitigate the risk of false-negative DDI predictions potentially caused by substrate-dependent Ki variation.

Mechanistic interpretation of conventional Michaelis-Menten parameters in a transporter system

The aim was to elucidate how steps in drug translocation by a solute carrier transporter impact Michaelis-Menten parameters Km, Ki, and Vmax. The first objective was to derive a model for carrier-mediated substrate translocation and perform sensitivity analysis with regard to the impact of individual microrate constants on Km, Ki, and Vmax. The second objective was to compare underpinning microrate constants between compounds translocated by the same transporter. Equations for Km, Ki, and Vmax were derived from a six-state model involving unidirectional transporter flipping and reconfiguration. This unidirectional model is applicable to co-transporter type solute carriers, like the apical sodium-dependent bile acid transporter (ASBT) and the proton-coupled peptide cotransporter (PEPT1). Sensitivity analysis identified the microrate constants that impacted Km, Ki, and Vmax. Compound comparison using the six-state model employed regression to identify microrate constant values that can explain observed Km and Vmax values. Results yielded some expected findings, as well as some unanticipated effects of microrate constants on Km, Ki, and Vmax. Km and Ki were found to be equal for inhibitors that are also substrates. Additionally, microrate constant values for certain steps in transporter functioning influenced Km and Vmax to be low or high.

Interaction of digitalis-like compounds with liver uptake transporters NTCP, OATP1B1, and OATP1B3

Digitalis-like compounds (DLCs) such as digoxin, digitoxin, and ouabain, also known as cardiac glycosides, are among the oldest pharmacological treatments for heart failure. The compounds have a narrow therapeutic window, while at the same time, DLC pharmacokinetics is prone to drug-drug interactions at the transport level. Hepatic transporters organic anion transporting polypeptide (OATP) 1B1, OATP1B3, and Na(+)-dependent taurocholate co-transporting polypeptide (NTCP) influence the disposition of a variety of drugs by mediating their uptake from blood into hepatocytes. The interaction of digoxin, digitoxin, and ouabain with hepatic uptake transporters has been studied before. However, here, we systematically investigated a much wider range of structurally related DLCs for their capability to inhibit or to be transported by these transporters in order to better understand the relation between the activity and chemical structure of this compound type. We studied the uptake and inhibitory potency of a series of 14 structurally related DLCs in Chinese hamster ovary cells expressing NTCP (CHO-NTCP) and human embryonic kidney cells expressing OATP1B1 and OATP1B3 (HEK-OATP1B1 and HEK-OATP1B3). The inhibitory effect of the DLCs was measured against taurocholic acid (TCA) uptake in CHO-NTCP cells and against uptake of β-estradiol 17-β-d-glucuronide (E217βG) in HEK-OATP1B1 and HEK-OATP1B3 cells. Proscillaridin A was the most effective inhibitor of NTCP-mediated TCA transport (IC50 = 22 μM), whereas digitoxin and digitoxigenin were the most potent inhibitors of OATP1B1 and OAPTP1B3, with IC50 values of 14.2 and 36 μM, respectively. Additionally, we found that the sugar moiety and hydroxyl groups of the DLCs play different roles in their interaction with NTCP, OATP1B1, and OATP1B3. The sugar moiety decreases the inhibition of NTCP and OATP1B3 transport activity, whereas it enhances the inhibitory potency against OATP1B1. Moreover, the hydroxyl group at position 12 reinforces the inhibition of NTCP but decreases the inhibition of OATP1B1 and OATP1B3. To investigate whether DLCs can be translocated, we quantified their uptake in transporter-expressing cells by LC-MS. We demonstrated that convallatoxin, ouabain, dihydroouabain, and ouabagenin are substrates of OATP1B3. No transport was observed for the other compounds in any of the studied transporters. In summary, this work provides a step toward an improved understanding of the interaction of DLCs with three major hepatic uptake transporters. Ultimately, this can be of use in the development of DLCs that are less prone to transporter-mediated drug-drug interactions.

Hepatic uptake of atorvastatin: influence of variability in transporter expression on uptake clearance and drug-drug interactions

Differences in the expression and function of the organic anion transporting polypeptide (OATP) transporters contribute to interindividual variability in atorvastatin clearance. However, the importance of the bile acid transporter sodium taurocholate cotransporting polypeptide (NTCP, SLC10A1) in atorvastatin uptake clearance (CLupt) is not yet clarified. To elucidate this issue, we investigated the relative contribution of NTCP, OATP1B1, OATP1B3, and OATP2B1 to atorvastatin CLupt in 12 human liver samples. The impact of inhibition on atorvastatin CLupt was also studied, using inhibitors of different isoform specificities. Expression levels of the four transport proteins were quantified by liquid chromatography tandem mass spectrometry. These data, together with atorvastatin in vitro kinetics, were used to predict the maximal transport activity (MTA) and interindividual differences in CLupt of each transporter in vivo. Subsequently, hepatic uptake impairment on coadministration of five clinically interacting drugs was predicted using in vitro inhibitory potencies. NTCP and OATP protein expression varied 3.7- to 32-fold among the 12 sample donors. The rank order in expression was OATP1B1 > OATP1B3 ≈ NTCP ≈ OATP2B1. NTCP was found to be of minor importance in atorvastatin disposition. Instead, OATP1B1 and OATP1B3 were confirmed as the major atorvastatin uptake transporters. The average contribution to atorvastatin uptake was OATP1B1 > OATP1B3 >> OATP2B1 > NTCP, although this rank order varied among individuals. The interindividual differences in transporter expression and CLupt resulted in marked differences in drug-drug interactions due to isoform-specific inhibition. We conclude that this variation should be considered in in vitro to in vivo extrapolations.

A perspective on the prediction of drug pharmacokinetics and disposition in drug research and development

Prediction of human pharmacokinetics of new drugs, as well as other disposition attributes, has become a routine practice in drug research and development. Prior to the 1990s, drug disposition science was used in a mostly descriptive manner in the drug development phase. With the advent of in vitro methods and availability of human-derived reagents for in vitro studies, drug-disposition scientists became engaged in the compound design phase of drug discovery to optimize and predict human disposition properties prior to nomination of candidate compounds into the drug development phase. This has reaped benefits in that the attrition rate of new drug candidates in drug development for reasons of unacceptable pharmacokinetics has greatly decreased. Attributes that are predicted include clearance, volume of distribution, half-life, absorption, and drug-drug interactions. In this article, we offer our experience-based perspectives on the tools and methods of predicting human drug disposition using in vitro and animal data.

Development of Novel, 384-Well High-Throughput Assay Panels for Human Drug Transporters

Transporter proteins are known to play a critical role in affecting the overall absorption, distribution, metabolism, and excretion characteristics of drug candidates. In addition to efflux transporters (P-gp, BCRP, MRP2, etc.) that limit absorption, there has been a renewed interest in influx transporters at the renal (OATs, OCTs) and hepatic (OATPs, BSEP, NTCP, etc.) organ level that can cause significant clinical drug-drug interactions (DDIs). Several of these transporters are also critical for hepatobiliary disposition of bilirubin and bile acid/salts, and their inhibition is directly implicated in hepatic toxicities. Regulatory agencies took action to address transporter-mediated DDI with the goal of ensuring drug safety in the clinic and on the market. To meet regulatory requirements, advanced bioassay technology and automation solutions were implemented for high-throughput transporter screening to provide structure-activity relationship within lead optimization. To enhance capacity, several functional assay formats were miniaturized to 384-well throughput including novel fluorescence-based uptake and efflux inhibition assays using high-content image analysis as well as cell-based radioactive uptake and vesicle-based efflux inhibition assays. This high-throughput capability enabled a paradigm shift from studying transporter-related issues in the development space to identifying and dialing out these concerns early on in discovery for enhanced mechanism-based efficacy while circumventing DDIs and transporter toxicities.

Substrate-dependent inhibition of organic anion transporting polypeptide 1B1: comparative analysis with prototypical probe substrates estradiol-17β-glucuronide, estrone-3-sulfate, and sulfobromophthalein

Organic anion transporting polypeptide (OATP) 1B1 plays an important role in the hepatic uptake of many drugs, and the evaluation of OATP1B1-mediated drug-drug interactions (DDIs) is emphasized in the latest DDI (draft) guidance documents from U.S. and E.U. regulatory agencies. It has been suggested that some OATP1B1 inhibitors show a discrepancy in their inhibitory potential, depending on the substrates used in the cell-based assay. In this study, inhibitory effects of 14 compounds on the OATP1B1-mediated uptake of the prototypical substrates [³H]estradiol-17β-glucuronide (E₂G), [³H]estrone-3-sulfate (E₁S), and [³H]sulfobromophthalein (BSP) were studied in OATP1B1-transfected cells. Inhibitory potencies of tested compounds varied depending on the substrates. Ritonavir, gemfibrozil, and erythromycin caused remarkable substrate-dependent inhibition with up to 117-, 14-, and 13-fold difference in their IC₅₀ values, respectively. Also, the clinically relevant OATP inhibitors rifampin and cyclosporin A exhibited up to 12- and 6-fold variation in their IC₅₀ values, respectively. Regardless of the inhibitors tested, the most potent OATP1B1 inhibition was observed when [³H]E₂G was used as a substrate. Mutual inhibition studies of OATP1B1 indicated that E₂G and E₁S competitively inhibited each other, whereas BSP noncompetitively inhibited E₂G uptake. In addition, BSP inhibited E₁S in a competitive manner, but E₁S caused an atypical kinetics on BSP uptake. This study showed substrate-dependent inhibition of OATP1B1 and demonstrated that E₂G was the most sensitive in vitro OATP1B1 probe substrate among three substrates tested. This will give us an insight into the assessment of clinically relevant OATP1B1-mediated DDI in vitro with minimum potential of false-negative prediction.

Cynomolgus monkey as a potential model to assess drug interactions involving hepatic organic anion transporting polypeptides: in vitro, in vivo, and in vitro-to-in vivo extrapolation

Organic anion-transporting polypeptides (OATP) 1B1, 1B3, and 2B1 can serve as the loci of drug-drug interactions (DDIs). In the present work, the cynomolgus monkey was evaluated as a potential model for studying OATP-mediated DDIs. Three cynomolgus monkey OATPs (cOATPs), with a high degree of amino acid sequence identity (91.9, 93.5, and 96.6% for OATP1B1, 1B3, and 2B1, respectively) to their human counterparts, were cloned, expressed, and characterized. The cOATPs were stably transfected in human embryonic kidney cells and were functionally similar to the corresponding human OATPs (hOATPs), as evident from the similar uptake rate of typical substrates (estradiol-17β-d-glucuronide, cholecystokinin octapeptide, and estrone-3-sulfate). Moreover, six known hOATP inhibitors exhibited similar IC(50) values against cOATPs. To further evaluate the appropriateness of the cynomolgus monkey as a model, a known hOATP substrate [rosuvastatin (RSV)]-inhibitor [rifampicin (RIF)] pair was examined in vitro; the monkey-derived parameters (RSV K(m) and RIF IC(50)) were similar (within 3.5-fold) to those obtained with hOATPs and human primary hepatocytes. In vivo, the area under the plasma concentration-time curve of RSV (3 mg/kg, oral) given 1 hour after a single RIF dose (15 mg/kg, oral) was increased 2.9-fold in cynomolgus monkeys, consistent with the value (3.0-fold) reported in humans. A number of in vitro-in vivo extrapolation approaches, considering the fraction of the pathways affected and free versus total inhibitor concentrations, were also explored. It is concluded that the cynomolgus monkey has the potential to serve as a useful model for the assessment of OATP-mediated DDIs in a nonclinical setting.