Substantially Increased Plasma Coproporphyrin-I Concentrations Associated With OATP1B1*15 Allele in Japanese General Population

Utility of Coproporphyrin-I Determination in First-in-Human Study for Early Evaluation of OATP1B Inhibitory Potential Based on Investigation of Ensitrelvir, an Oral SARS-CoV-2 3C-Like Protease Inhibitor

Coproporphyrin-I (CP-I) has been investigated as an endogenous biomarker of organic anion transporting polypeptide (OATP) 1B. Here, we determined the CP-I concentrations in a cocktail drug-drug interaction (DDI) study of ensitrelvir to evaluate the OATP1B inhibitory potential because ensitrelvir had increased plasma concentrations of rosuvastatin in this study, raising concerns about breast cancer resistance protein and OATP1B inhibition. Furthermore, CP-I concentrations were compared between active and placebo groups in a first-in-human (FIH) study of ensitrelvir to verify whether the OATP1B inhibitory potential could be estimated at an early drug development stage. In the cocktail DDI study, CP-I did not differ between with/without administration of ensitrelvir, indicating that ensitrelvir has no OATP1B inhibitory effect. Although there were some individual variabilities in CP-I concentrations among the treatment groups in the FIH study, the normalization of CP-I concentrations with pre-dose values minimized these variabilities, suggesting that this normalized method would be helpful for comparing the CP-I from different participants. Finally, we concluded that CP-I concentrations were not affected by ensitrelvir in the FIH study. These results suggested that the CP-I determination in an FIH study and its normalized method can be useful for an early evaluation of the OATP1B-mediated DDI potential in humans.

Relationship of plasma 3-carboxy-4-methyl-5-propyl-2-furanpropanoic acid concentration with OATP1B activity in patients with chronic kidney disease

Organic anion-transporting polypeptides (OATP)1B are drug transporters mainly expressed in the sinusoidal membrane. Many studies have suggested that OATP1B activity is affected by genetic factor, the uremic toxin 3-carboxy-4-methyl-5-propyl-2-furanpropanoic acid (CMPF), and inflammatory cytokines, such as tumor necrosis factor-α (TNF-α) and interleukin-6 (IL-6). Coproporphyrin-I (CP-I) is spotlighted as a highly accurate endogenous substrate of OATP1B. We previously reported a positive correlation between plasma CMPF and CP-I concentrations in patients with chronic kidney disease (CKD). The present study evaluated the impact of genetic polymorphisms, CMPF, IL-6, TNF-α, and estimated glomerular filtration rate (eGFR) on individual differences in OATP1B activity in patients with CKD. Seventy-three patients with CKD who received kidney transplant at least 3 months earlier were analyzed. Plasma CP-I concentration was higher in OATP1B1*15 carriers than in non-carriers. In all patients, CP-I did not correlate significantly with CMPF, IL-6, TNF-α, or eGFR. However, when the dataset was cut off at CMPF concentration of 8 and 7 μg/mL, 4 μg/mL, 3 μg/mL or 2 μg/mL, CMPF correlated positively with CP-I, and correlation coefficient tended to be higher as plasma CMPF concentration was lower. In conclusion, OATP1B1*15 impacted OATP1B activity in patients with CKD, but IL-6 and TNF-α did not. However, the impact of CMPF on OATP1B activity was limited to low CMPF concentrations, and the effect could be saturated at high concentrations. When prescribing an OATP1B substrate drug for patients with CKD, the OATP1B1*15 carrier status and plasma CMPF concentration may need to be considered to decide the dose regimen.

Progress in the Quantitative Assessment of Transporter-Mediated Drug-Drug Interactions Using Endogenous Substrates in Clinical Studies

Variations in drug transporter activities, caused by genetic polymorphism and drug-drug interactions (DDIs), alter the systemic exposure of substrate drugs, leading to differences in drug responses. Recently, some endogenous substrates of drug transporters, particularly the solute carrier family transporters such as OATP1B1, OATP1B3, OAT1, OAT3, OCT1, OCT2, MATE1, and MATE2-K, have been identified to investigate variations in drug transporters in humans. Clinical data obtained support their performance as surrogate probes in terms of specificity and reproducibility. Pharmacokinetic parameters of the endogenous biomarkers depend on the genotypes of drug transporters and the systemic exposure to perpetrator drugs. Furthermore, the development of physiologically based pharmacokinetic models for the endogenous biomarkers has enabled a top-down approach to obtain insights into the effect of perpetrators on drug transporters and to more precisely simulate the DDI with victim drugs, including probe drugs. The endogenous biomarkers can address the uncertainty in the DDI prediction in the preclinical and early phases of clinical development and have the potential to fulfill regulatory requirements. Therefore, the endogenous biomarkers should be able to predict disease effects on the variations in drug transporter activities observed in patients. This mini-review focuses on recent progress in the identification and use of the endogenous drug transporter substrate biomarkers and their application in drug development. SIGNIFICANCE STATEMENT: Advances in analytical methods have enabled the identification of endogenous substrates of drug transporters. Changes in the pharmacokinetic parameters (Cmax, AUC, or CLR) of these endogenous biomarkers relative to baseline values can serve as a quantitative index to assess variations in drug transporter activities during clinical studies and thereby provide more precise DDI predictions.

Characterization of endogenous markers of hepatic function in patients receiving itraconazole treatment for prophylaxis of deep mycosis

BACKGROUND

Long-term use of itraconazole (ITZ) is associated with a risk of inducing hepatotoxicity. This study aimed to evaluate the associations of plasma concentrations of ITZ and its hydroxylated metabolite (OH-ITZ) with endogenous markers of hepatic function.

METHODS

Thirty six patients treated with oral ITZ solution for prophylaxis of deep mycosis were enrolled. Plasma concentrations of ITZ and OH-ITZ were determined on the 14th day or later after administration of ITZ. Their associations with endogenous marker levels of hepatic function including plasma coproporphyrin (CP)-I and OATP1B1 genotypes were assessed.

RESULTS

The serum level of total bilirubin (T-Bil) was moderately correlated with the plasma concentration of total ITZ (tITZ) and OH-ITZ (tOH-ITZ). T-Bil elevation above 0.3 mg/dL was observed in 19% of patients, although statistically significant difference was not identified. The plasma concentration of tITZ had no correlation with other endogenous markers levels including AST, ALT, albumin, and plasma CP-I. The serum AST and plasma CP-I levels were correlated with the plasma concentration of free OH-ITZ (fOH-ITZ). T-Bil and plasma CP-I, a marker of OATP1B1 activity, were not correlated with each other, and neither was associated with the OATP1B1 genotypes.

CONCLUSIONS

Plasma ITZ and OH-ITZ had a positive association with T-Bil. The patients with a higher fOH-ITZ level had lower OATP1B1 activity on the basis of plasma CP-I level. ITZ and OH-ITZ have the potential to slightly increase endogenous marker levels of hepatic function, although most likely by different mechanisms.

Positive correlation between organic anion transporter 1B function indicated by plasma concentration of coproporphyrin-I and blood concentration of cyclosporin A in real-world patients

AIMS

Cyclosporin A (CyA) has potent inhibitory activity on organic anion transporting polypeptide 1B (OATP1B), causing drug-drug interactions with its substrate drugs. 3-carboxy-4-methyl-5-propyl-2-furanpropionate (CMPF), a uraemic toxin, has also been suggested to inhibit OATP1B activity. Recent study has identified coproporphyrin-I (CP-I) as a specific endogenous substrate for OATP1B, which is useful to indicate OATP1B activity. We investigated the relationship of CP-I with CyA and CMPF concentrations in patients taking CyA.

METHODS

In total, 121 blood samples from 74 patients who took CyA and underwent routine therapeutic drug monitoring were divided into trough and peak samples.

RESULTS

CyA and CP-I concentrations were significantly higher in peak samples than in trough samples. A positive correlation between CP-I and CyA concentrations was found in all samples and in trough and peak samples, while no correlation was observed between CP-I and CMPF concentrations. Multiple regression analysis identified CyA and C-reactive protein concentrations as independent factors affecting CP-I concentration, with blood CyA concentration having markedly greater contribution to plasma CP-I concentration.

CONCLUSION

The present study suggests that CyA inhibits OATP1B activity in a concentration-dependent manner in clinical setting, and that dose adjustment of OATP1B substrate drugs coadministered with CyA according to plasma CMPF concentration may not be necessary.

The role of OATP1B1 and OATP1B3 transporter polymorphisms in drug disposition and response to anticancer drugs: a review of the recent literature

INTRODUCTION

Members of the solute carrier family of organic anion transporting polypeptides are responsible for the cellular uptake of a broad range of endogenous compounds and xenobiotics in multiple tissues. In particular, the polymorphic transporters OATP1B1 and OATP1B3 are highly expressed in the liver and have been identified as critical regulators of hepatic eliminaton. As these transporters are also expressed in cancer cells, the function alteration of these proteins have important consequences for an individual's susceptibility to certain drug-induced side effects, drug-drug interactions, and treatment efficacy.

AREAS COVERED

In this mini-review, we provide an update of this rapidly emerging field, with specific emphasis on the direct contribution of genetic variants in OATP1B1 and OATP1B3 to the transport of anticancer drugs, the role of these carriers in regulation of their disposition and toxicity profiles, and recent advances in attempts to integrate information on transport function in patients to derive individualized treatment strategies.

EXPERT OPINION

Based on currently available data, it appears imperative that different aspects of disease, physiology, and drugs of relevance should be evaluated along with an individual's genetic signature, and that tools such as biomarker levels can be implemented to achieve the most reliable prediction of clinically relevant pharmacodynamic endpoints.

Recent advances in the translation of drug metabolism and pharmacokinetics science for drug discovery and development

Drug metabolism and pharmacokinetics (DMPK) is an important branch of pharmaceutical sciences. The nature of ADME (absorption, distribution, metabolism, excretion) and PK (pharmacokinetics) inquiries during drug discovery and development has evolved in recent years from being largely descriptive to seeking a more quantitative and mechanistic understanding of the fate of drug candidates in biological systems. Tremendous progress has been made in the past decade, not only in the characterization of physiochemical properties of drugs that influence their ADME, target organ exposure, and toxicity, but also in the identification of design principles that can minimize drug-drug interaction (DDI) potentials and reduce the attritions. The importance of membrane transporters in drug disposition, efficacy, and safety, as well as the interplay with metabolic processes, has been increasingly recognized. Dramatic increases in investments on new modalities beyond traditional small and large molecule drugs, such as peptides, oligonucleotides, and antibody-drug conjugates, necessitated further innovations in bioanalytical and experimental tools for the characterization of their ADME properties. In this review, we highlight some of the most notable advances in the last decade, and provide future perspectives on potential major breakthroughs and innovations in the translation of DMPK science in various stages of drug discovery and development.

Physiologically-based pharmacokinetic model-based translation of OATP1B-mediated drug-drug interactions from coproporphyrin I to probe drugs

The accurate prediction of OATP1B-mediated drug-drug interactions (DDIs) is challenging for drug development. Here, we report physiologically-based pharmacokinetic (PBPK) model analysis for clinical DDI data generated in heathy subjects who received oral doses of cyclosporin A (CysA; 20 and 75 mg) as an OATP1B inhibitor, and the probe drugs (pitavastatin, rosuvastatin and valsartan). PBPK models of CysA and probe compounds were combined assuming inhibition of hepatic uptake of endogenous coproporphyrin I (CP-I) by CysA. In vivo K_i of unbound CysA for OATP1B (K_i,OATP1B ), and the overall intrinsic hepatic clearance per body weight of CP-I (CL_int,all,unit ) were optimized to account for the CP-I data (K_i,OATP1B , 0.657 ± 0.048 nM; CL_int,all,unit , 57.0 ± 6.3 L/h/kg). DDI simulation using K_i,OATP1B reproduced the dose-dependent effect of CysA (20 and 75 mg) and the dosing interval (1 h and 3 h) on the time profiles of blood concentrations of pitavastatin and rosuvastatin, but DDI simulation using in vitro K_i,OATP1B failed. The Cluster Gauss-Newton method was used to conduct parameter optimization using 1,000 initial parameter sets for the seven pharmacokinetic parameters of CP-I (β, CL_int,all , F_a F_g , R_dif , f_bile , f_syn , and v_syn ), and K_i,OATP1B , and K_i,MRP2 of CysA. Based on the accepted 498 parameter sets, the range of CL_int,all and K_i,OATP1B was narrowed, with coefficients of variation (CVs) of 9.3% and 11.1%, respectively, indicating that these parameters were practically identifiable. These results suggest that PBPK model analysis of CP-I is a promising translational approach to predict OATP1B-mediated DDIs in drug development.

Quantitative prediction of pharmacokinetic properties of drugs in humans: Recent advance in in vitro models to predict the impact of efflux transporters in the small intestine and blood-brain barrier

Efflux transport systems are essential to suppress the absorption of xenobiotics from the intestinal lumen and protect the critical tissues at the blood-tissue barriers, such as the blood-brain barrier. The function of drug efflux transport is dominated by various transporters. Accumulated clinical evidences have revealed that genetic variations of the transporters, together with coadministered drugs, affect the expression and/or function of transporters and subsequently the pharmacokinetics of substrate drugs. Thus, in the preclinical stage of drug development, quantitative prediction of the impact of efflux transporters as well as that of uptake transporters and metabolic enzymes on the pharmacokinetics of drugs in humans has been performed using various in vitro experimental tools. Various kinds of human-derived cell systems can be applied to the precise prediction of drug transport in humans. Mathematical modeling consisting of each intrinsic metabolic or transport process enables us to understand the disposition of drugs both at the organ level and at the level of the whole body by integrating a variety of experimental results into model parameters. This review focuses on the role of efflux transporters in the intestinal absorption and brain distribution of drugs, in addition to recent advances in predictive tools and methodologies.

Organic Anion Transporting Polypeptide 2B1 (OATP2B1) Genetic Variants: In Vitro Functional Characterization and Association With Circulating Concentrations of Endogenous Substrates

Organic anion transporting polypeptide 2B1 (OATP2B1, gene SLCO2B1) is an uptake transporter that is thought to determine drug disposition and in particular, the oral absorption of medications. At present, the clinical relevance of SLCO2B1 genetic variation on pharmacokinetics is poorly understood. We sought to determine the functional activity of 5 of the most common missense OATP2B1 variants (c.76_84del, c.601G>A, c.917G>A, c.935G>A, and c.1457C>T) and a predicted dysfunctional variant (c.332G>A) in vitro. Furthermore, we measured the basal plasma concentrations of endogenous OATP2B1 substrates, namely estrone sulfate, dehydroepiandrosterone sulfate (DHEAS), pregnenolone sulfate, coproporphyrin I (CPI), and CPIII, and assessed their relationships with SLCO2B1 genotypes in 93 healthy participants. Compared to reference OATP2B1, the transport activities of the c.332G>A, c.601G>A and c.1457C>T variants were reduced among the substrates examined (estrone sulfate, DHEAS, CPI, CPIII and rosuvastatin), although there were substrate-dependent effects. Lower transport function of OATP2B1 variants could be explained by diminished cell surface expression. Other OATP2B1 variants (c.76-84del, c.917G>A and c.935G>A) had similar activity to the reference transporter. In the clinical cohort, the SLCO2B1 c.935G>A allele was associated with both higher plasma CPI (42%) and CPIII (31%) concentrations, while SLCO2B1 c.917G>A was linked to lower plasma CPIII by 28% after accounting for the effects of age, sex, and SLCO1B1 genotypes. No association was observed between SLCO2B1 variant alleles and estrone sulfate or DHEAS plasma concentrations, however 45% higher plasma pregnenolone sulfate level was associated with SLCO2B1 c.1457C>T. Taken together, we found that the impacts of OATP2B1 variants on transport activities in vitro were not fully aligned with their associations to plasma concentrations of endogenous substrates in vivo. Additional studies are required to determine whether circulating endogenous substrates reflect OATP2B1 activity.

Performance of Plasma Coproporphyrin I and III as OATP1B1 Biomarkers in Humans

A previous study in 356 healthy Finnish volunteers showed that glycochenodeoxycholate 3‐O‐glucuronide (GCDCA‐3G) and glycodeoxycholate 3‐O‐glucuronide (GDCA‐3G) are promising biomarkers of organic anion transporting polypeptide 1B1 (OATP1B1). In the same cohort, we now evaluated the performances of two other OATP1B1 biomarkers, coproporphyrin I (CPI) and III (CPIII), and compared them with GCDCA‐3G and GDCA‐3G. Based on decreased (*5 and *15) and increased (*14 and *20) function SLCO1B1 haplotypes, we stratified the participants to poor, decreased, normal, increased, and highly increased OATP1B1 function groups. Fasting plasma CPI concentration was 68% higher in the poor (95% confidence interval, 44%, 97%; P = 1.74 × 10⁻¹⁰), 7% higher in the decreased (0%, 15%; P = 0.0385), 10% lower in the increased (3%, 18%; P = 0.0087), and 23% lower in the highly increased (1%, 40%; P = 0.0387) function group than in the normal function group. CPIII concentration was 27% higher (7%, 51%; P = 0.0071) in the poor function group than in the normal function group. CPI and CPIII detected poor OATP1B1 function with areas under the precision‐recall curve (AUPRC) of 0.388 (95% confidence interval, 0.197, 0.689) and 0.0798 (0.0485, 0.203), and receiver operating characteristic curve (AUROC) of 0.888 (0.851, 0.919) and 0.731 (0.682, 0.776). The AUPRC and AUROC of GCDCA‐3G were, however, 0.389 (0.258, 0.563) and 0.100 (−0.0046, 0.204; P = 0.0610) larger than those of CPI, and 0.697 (0.555, 0.831) and 0.257 (0.141, 0.373; P < 0.0001) larger than those of CPIII. In conclusion, these data indicate that plasma CPI outperforms CPIII in detecting altered OATP1B1 function, but GCDCA‐3G is an even more sensitive OATP1B1 biomarker than CPI.

Endogenous Biomarkers for SLC Transporter-Mediated Drug-Drug Interaction Evaluation

Membrane transporters play an important role in the absorption, distribution, metabolism, and excretion of xenobiotic substrates, as well as endogenous compounds. The evaluation of transporter-mediated drug-drug interactions (DDIs) is an important consideration during the drug development process and can guide the safe use of polypharmacy regimens in clinical practice. In recent years, several endogenous substrates of drug transporters have been identified as potential biomarkers for predicting changes in drug transport function and the potential for DDIs associated with drug candidates in early phases of drug development. These biomarker-driven investigations have been applied in both preclinical and clinical studies and proposed as a predictive strategy that can be supplanted in order to conduct prospective DDIs trials. Here we provide an overview of this rapidly emerging field, with particular emphasis on endogenous biomarkers recently proposed for clinically relevant uptake transporters.

Factors Influencing Plasma Coproporphyrin-I Concentration as Biomarker of OATP1B Activity in Patients With Rheumatoid Arthritis

Organic anion transporting polypeptides (OATPs) 1B are drug transporters mainly expressed in the sinusoidal membrane. In previous reports, genetic factor, 3-carboxy-4-methyl-5-propyl-2-furanpropanoic acid (CMPF), which is one of the uremic toxins, inflammatory cytokines such as tumor necrosis factor-α (TNF-α) and interleukin-6 (IL-6) decreased OATP1B1 activity in vitro, but in vivo effects of these factors have not been elucidated. Plasma coproporphyrin-I (CP-I) is spotlighted as a highly accurate endogenous substrate of OATP1B. This study focused on patients with rheumatoid arthritis (RA) and evaluated the influence of several factors comprising gene polymorphisms, uremic toxins, and inflammatory cytokines on OATP1B activity using plasma CP-I concentration. Thirty-seven outpatients with RA who satisfied the selection criteria were analyzed at the time of recruitment (baseline) and at the next visit. OATP1B1*15 carriers tended to have higher CP-I concentration compared with noncarriers. Plasma CP-I correlated positively with CMPF concentration, but did not correlate with IL-6 or TNF-α concentration. Multiple logistic regression analysis by stepwise selection identified plasma CMPF concentration and OATP1B1*15 allele as significant factors independently affecting plasma CP-I concentration at baseline and at the next visit, respectively. In conclusion, the present results suggest that inflammatory cytokines do not have clinically significant effects on OATP1B activity, whereas the effects of genetic polymorphisms and uremic toxins should be considered.

Minimal contribution of the hepatic uptake transporter OATP1B1 to the inter-individual variability in SN-38 pharmacokinetics in cancer patients without severe renal failure

PURPOSE

SN-38, a pharmacologically active metabolite of irinotecan, is taken up into hepatocytes by organic anion transporting polypeptide (OATP) 1B1. The effects of functional OATP1B1 521T>C on the pharmacokinetics of SN-38 remain controversial. Here, we prospectively examined the effects of OATP1B1 function on the area under the plasma total or unbound concentration-time curve (tAUC or uAUC) of SN-38 by assessing OATP1B1 521T>C and the plasma levels of endogenous OATP1B1 substrates, coproporphyrin (CP)-I and III, in cancer patients treated with irinotecan.

METHODS

We enrolled cancer patients who were treated with an irinotecan-containing regimen and did not have severe renal failure. The total and unbound concentrations of SN-38 in the plasma were measured by high-performance liquid chromatography. AUC values were calculated and normalized to the actual irinotecan dose (AUC/dose). The OATP1B1 521T>C was analyzed by direct sequencing. Concentrations of the endogenous substrates in plasma before irinotecan treatment (baseline) were determined by liquid chromatography with tandem mass spectrometry.

RESULTS

Twenty-two patients with a median estimated glomerular filtration rate of 74.8 mL/min (range 32.6-99.6) were examined. Both tAUC/dose and uAUC/dose were associated with the grade of neutropenia; however, they were not associated with OATP1B1 521T>C or baseline CP-I and III levels. It is worth noting that these baseline concentrations were significantly higher in patients with OATP1B1 521C, supporting functional changes in OATP1B1.

CONCLUSION

The contribution of OATP1B1 activity to inter-patient variability in the systemic exposure to SN-38 is likely minimal in patients without severe renal failure.

Relationship of hemoglobin level and plasma coproporphyrin-I concentrations as an endogenous probe for phenotyping OATP1B

Plasma coproporphyrin‐I (CP‐I) concentration is used as a sensitive and selective endogenous probe for phenotyping organic anion transporting polypeptides 1B (OATP1B) activity in many studies. CP‐I is produced in the process of heme synthesis, but the relationship between plasma CP‐I concentrations and heme synthesis activity is unknown. In this study, we evaluated the relationship between plasma CP‐I concentration and hemoglobin level as a biomarker of heme synthesis activity. The data of 391 subjects selected from the Japanese general population were analyzed. One hundred twenty‐six participants had OATP1B1*15 allele, 11 of whom were homozygous (OATP1B1*15/*15). Multiple regression analysis identified hemoglobin level as an independent variable associated with plasma CP‐I concentration (p < 0.0001). A significant positive correlation was observed between hemoglobin level and plasma CP‐I concentration in participants without OATP1B1*15 allele (n = 265; r_s = 0.35, p < 0.0001) and with OATP1B1*15 allele (n = 126; r_s =0.27, p = 0.0022). However, Kruskal–Wallis test showed no large difference in Kruskal–Wallis statistics between the distribution of plasma CP‐I concentrations and that of ratio of plasma CP‐I to hemoglobin among six OATP1B1 polymorphism groups. These findings suggest that the hemoglobin level seems to reflect biosynthesis of CP‐I. However, correction by hemoglobin level is not required when using basal plasma CP‐I concentration for phenotyping OATP1B activity.

Effects of Probenecid on Hepatic and Renal Disposition of Hexadecanedioate, an Endogenous Substrate of Organic Anion Transporting Polypeptide 1B in Rats

The aim of the present study was to investigate changes in plasma concentrations and tissue distribution of endogenous substrates of organic anion transporting polypeptide (OATP) 1B, hexadecanedioate (HDA), octadecanedioate (ODA), tetradecanedioate (TDA), and coproporphyrin-III, induced by its weak inhibitor, probenecid (PBD), in rats. PBD increased the plasma concentrations of these four compounds regardless of bile duct cannulation, whereas liver-to-plasma (K_p,liver) and kidney-to-plasma concentration ratios of HDA and TDA were reduced. Similar effects of PBD on plasma concentrations and K_p,liver of HDA, ODA, and TDA were observed in kidney-ligated rats, suggesting a minor contribution of renal disposition to the overall distribution of these three compounds. Tissue uptake clearance of deuterium-labeled HDA (d-HDA) in liver was 16-fold higher than that in kidney, and was reduced by 80% by PBD. This was compatible with inhibition by PBD of d-HDA uptake in isolated rat hepatocytes. Such inhibitory effects of PBD were also observed in the human OATP1B1-mediated uptake of d-HDA. Overall, the disposition of HDA is mainly determined by hepatic OATP-mediated uptake, which is inhibited by PBD. HDA might, thus, be a biomarker for OATPs minimally affected by urinary and biliary elimination in rats.