In vitro and in silico strategies to identify OATP1B1 inhibitors and predict clinical drug-drug interactions.

PURPOSE: To establish in vitro and in silico models that predict clinical drug-drug interactions (DDIs) with the OATP1B1 (SLCO1B1) transporter. METHODS: The inhibitory effect of 146 drugs and drug-like compounds on OATP1B1-mediated transport was studied in HEK293 cells. A computational model was developed to predict OATP1B1 inhibition. Concentration-dependent effects were investigated for six compounds; clinical DDIs were predicted by calculating change in exposure (i.e. R-values) in eight different ways. RESULTS: Sixty-five compounds were identified as OATP1B1 inhibitors at 20 µM. The computational model predicted the test set with 80% accuracy for inhibitors and 91% for non-inhibitors. In vitro-in vivo comparisons underscored the importance of using drugs with known clinical effects as references. Thus, reference drugs, cyclosporin A, gemfibrozil, and fenofibrate, provided an inhibition interval to which three antiviral drugs, atazanavir, lopinavir, and amprenavir, could be compared and their clinical DDIs with OATP1B1 classified. CONCLUSIONS: Twenty-two new OATP1B1 inhibitors were identified, a predictive OATP1B1 inhibition in silico model was developed, and successful predictions of clinical DDIs were obtained with OATP1B1.

Measurement of unbound drug exposure in brain: modeling of pH partitioning explains diverging results between the brain slice and brain homogenate methods.

Currently used methodology for determining unbound drug exposure in brain combines measurement of the total drug concentration in the whole brain in vivo with estimation of brain tissue binding from one of two available in vitro methods: equilibrium dialysis of brain homogenate and the brain slice uptake method. This study of 56 compounds compares the fraction of unbound drug in brain (f(u,brain)), determined using the brain homogenate method, with the unbound volume of distribution in brain (V(u,brain)), determined using the brain slice method. Discrepancies were frequent and were primarily related to drug pH partitioning, attributable to the preservation of cellular structures in the slice that are absent in the homogenate. A mathematical model for pH partitioning into acidic intracellular compartments was derived to predict the slice V(u,brain) from measurements of f(u,brain) and drug pK(a). This model allowed prediction of V(u,brain) from f(u,brain) within a 2.2-fold error range for 95% of the drugs compared with a 4.5-fold error range using the brain homogenate f(u,brain) method alone. The greatest discrepancies between the methods occurred with compounds that are actively transported into brain cells, including gabapentin, metformin, and prototypic organic cation transporter substrates. It was concluded that intrabrain drug distribution is governed by several diverse mechanisms in addition to nonspecific binding and that the slice method is therefore more reliable than the homogenate method. As an alternative, predictions of V(u,brain) can be made from homogenate f(u,brain) using the pH partition model presented, although this model does not take into consideration possible active brain cell uptake.

Endogenous gene and protein expression of drug-transporting proteins in cell lines routinely used in drug discovery programs.

The aim of this study was to investigate the gene and protein expression profiles of important drug-transporting proteins in human cell lines commonly used for studies of drug transport mechanisms. Human cell lines used to transiently or stably express single transporters [HeLa, human embryonic kidney (HEK) 293] and leukemia cell lines used to study drug resistance by ATP-binding cassette transporters (HL-60, K562) were investigated and compared with organotypic cell lines (HepG2, Saos-2, Caco-2, and Caco-2 TC7). For gene expression studies, real-time polymerase chain reaction was used, whereas monospecific polyclonal antibodies were generated and used to investigate protein expression by immunohistochemistry. Thirty-six transporters were studied for gene expression, and nine were studied for protein expression. The antibodies were validated using expression patterns in human tissues. Finally, the function of one ubiquitously expressed transporter, MCT1/SLC16A1, was investigated using [(14)C]lactic acid as a substrate. In general, the adherent cell lines (HeLa, HEK293) displayed low transporter expression, and the expression patterns were barely affected by transfection. The leukemia cell lines (K562, HL-60) and Saos-2 also had low endogenous transporter expression, whereas the organotypic cell lines (HepG2 and Caco-2) showed higher expression of some transporters. Comparison of gene and protein expression profiles gave poor correlations, but better agreement was obtained for antibodies with a good validation score, indicating that antibody quality was a significant variable. It is noteworthy that the monocarboxylic acid-transporting protein MCT1 was significantly expressed in all and was functional in most of the cell lines, indicating that MCT1 may be a confounding factor when the transport of small anionic drugs is investigated.

Structural requirements for drug inhibition of the liver specific human organic cation transport protein 1.

The liver-specific organic cation transport protein (OCT1; SLC22A1) transports several cationic drugs including the antidiabetic drug metformin and the anticancer agents oxaliplatin and imatinib. In this study, we explored the chemical space of registered oral drugs with the aim of studying the inhibition pattern of OCT1 and of developing predictive computational models of OCT1 inhibition. In total, 191 structurally diverse compounds were examined in HEK293-OCT1 cells. The assay identified 47 novel inhibitors and confirmed 15 previously known inhibitors. The enrichment of OCT1 inhibitors was seen in several drug classes including antidepressants. High lipophilicity and a positive net charge were found to be the key physicochemical properties for OCT1 inhibition, whereas a high molecular dipole moment and many hydrogen bonds were negatively correlated to OCT1 inhibition. The data were used to generate OPLS-DA models for OCT1 inhibitors; the final model correctly predicted 82% of the inhibitors and 88% of the noninhibitors of the test set.

Expression of thirty-six drug transporter genes in human intestine, liver, kidney, and organotypic cell lines.

This study was designed to quantitatively assess the mRNA expression of 36 important drug transporters in human jejunum, colon, liver, and kidney. Expression of these transporters in human organs was compared with expression in commonly used cell lines (Caco-2, HepG2, and Caki-1) originating from these organs to assess their value as in vitro transporter system models, and was also compared with data obtained from the literature on expression in rat tissues to assess species differences. Transporters that were highly expressed in the intestine included HPT1, PEPT1, BCRP, MRP2, and MDR1, whereas, in the liver, OCT1, MRP2, OATP-C, NTCP and BSEP were the main transporters. In the kidney, OAT1 was expressed at the highest levels, followed by OAT3, OAT4, MCT5, MDR1, MRP2, OCT2, and OCTN2. The best agreement between human tissue and the representative cell line was observed for human jejunum and Caco-2 cells. Expression in liver and kidney ortholog cell lines was not correlated with that in the associated tissue. Comparisons with rat transporter gene expression revealed significant species differences. Our results allowed a comprehensive quantitative comparison of drug transporter expression in human intestine, liver, and kidney. We suggest that it would be beneficial for predictive pharmacokinetic research to focus on the most highly expressed transporters. We hope that our comparison of rat and human tissue will help to explain the observed species differences in in vivo models, increase understanding of the impact of active transport processes on pharmacokinetics and distribution, and improve the quality of predictions from animal studies to humans.

A global drug inhibition pattern for the human ATP-binding cassette transporter breast cancer resistance protein (ABCG2).

In this article, we explore the entire structural space of registered drugs to obtain a global model for the inhibition of the drug efflux transporter breast cancer resistance protein (BCRP; ABCG2). For this purpose, the inhibitory effect of 123 structurally diverse drugs and drug-like compounds on mitoxantrone efflux was studied in Saos-2 cells transfected with human wild-type (Arg482) BCRP. The search for BCRP inhibitors throughout the drug-like chemical space resulted in the identification of 29 previously unknown inhibitors. The frequency of BCRP inhibition was 3 times higher for compounds reported to interact with other ATP-binding cassette (ABC) transporters than for compounds without reported ABC transporter affinity. An easily interpreted computational model capable of discriminating inhibitors from noninhibitors using only two molecular descriptors, octanol-water partition coefficient at pH 7.4 and molecular polarizability, was constructed. The discriminating power of this two-descriptor model was 93% for the training set and 79% for the test set, respectively. The results were supported by a global pharmacophore model and are in agreement with a two-step mechanism for the inhibition of BCRP, where both the drug's capacity to insert into the cell membrane and to interact with the inhibitory binding site of the transporter are important.