Clinical and Translational Science (Jun 2022)

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

  • Tatsuki Mochizuki,
  • Yasunori Aoki,
  • Takashi Yoshikado,
  • Kenta Yoshida,
  • Yurong Lai,
  • Hideki Hirabayashi,
  • Yoshiyuki Yamaura,
  • Kevin Rockich,
  • Kunal Taskar,
  • Tadayuki Takashima,
  • Xiaoyan Chu,
  • Maciej J. Zamek‐Gliszczynski,
  • Jialin Mao,
  • Kazuya Maeda,
  • Kenichi Furihata,
  • Yuichi Sugiyama,
  • Hiroyuki Kusuhara

DOI
https://doi.org/10.1111/cts.13272
Journal volume & issue
Vol. 15, no. 6
pp. 1519 – 1531

Abstract

Read online

Abstract The accurate prediction of OATP1B‐mediated drug–drug interactions (DDIs) is challenging for drug development. Here, we report a 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 Ki of unbound CysA for OATP1B (Ki,OATP1B), and the overall intrinsic hepatic clearance per body weight of CP‐I (CLint,all,unit) were optimized to account for the CP‐I data (Ki,OATP1B, 0.536 ± 0.041 nM; CLint,all,unit, 41.9 ± 4.3 L/h/kg). DDI simulation using Ki,OATP1B reproduced the dose‐dependent effect of CysA (20 and 75 mg) and the dosing interval (1 and 3 h) on the time profiles of blood concentrations of pitavastatin and rosuvastatin, but DDI simulation using in vitro Ki,OATP1B failed. The Cluster Gauss–Newton method was used to conduct parameter optimization using 1000 initial parameter sets for the seven pharmacokinetic parameters of CP‐I (β, CLint, all, FaFg, Rdif, fbile, fsyn, and vsyn), and Ki,OATP1B and Ki,MRP2 of CysA. Based on the accepted 546 parameter sets, the range of CLint, all and Ki,OATP1B was narrowed, with coefficients of variation of 12.4% and 11.5%, 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.