Use of mechanistic modeling to assess interindividual variability and interspecies differences in active uptake in human and rat hepatocytes

In Vitro-In Vivo Extrapolation and Scaling Factors for Clearance of Human and Preclinical Species with Liver Microsomes and Hepatocytes

In vitro-in vivo extrapolation ((IVIVE) and empirical scaling factors (SF) of human intrinsic clearance (CL_int) were developed using one of the largest dataset of 455 compounds with data from human liver microsomes (HLM) and human hepatocytes (HHEP). For extended clearance classification system (ECCS) class 2/4 compounds, linear SFs (SF_lin) are approximately 1, suggesting enzyme activities in HLM and HHEP are similar to those in vivo under physiological conditions. For ECCS class 1A/1B compounds, a unified set of SFs was developed for CL_int. These SFs contain both SF_lin and an exponential SF (SF_β) of fraction unbound in plasma (f_u,p). The unified SFs for class 1A/1B eliminate the need to identify the transporters involved prior to clearance prediction. The underlying mechanisms of these SFs are not entirely clear at this point, but they serve practical purposes to reduce biases and increase prediction accuracy. Similar SFs have also been developed for preclinical species. For HLM-HHEP disconnect (HLM > HHEP) ECCS class 2/4 compounds that are mainly metabolized by cytochrome P450s/FMO, HLM significantly overpredicted in vivo CL_int, while HHEP slightly underpredicted and geometric mean of HLM and HHEP slightly overpredicted in vivo CL_int. This observation is different than in rats, where rat liver microsomal CL_int correlates well with in vivo CL_int for compounds demonstrating permeability-limited metabolism. The good CL_int IVIVE developed using HLM and HHEP helps build confidence for prospective predictions of human clearance and supports the continued utilization of these assays to guide structure-activity relationships to improve metabolic stability.

Mathematical models for biomarker calculation of drug-induced liver injury in humans and experimental models based on gadoxetate enhanced magnetic resonance imaging

BACKGROUND

Drug induced liver injury (DILI) is a major concern when developing new drugs. A promising biomarker for DILI is the hepatic uptake rate of the contrast agent gadoxetate. This rate can be estimated using a novel approach combining magnetic resonance imaging and mathematical modeling. However, previous work has used different mathematical models to describe liver function in humans or rats, and no comparative study has assessed which model is most optimal to use, or focused on possible translatability between the two species.

AIMS

Our aim was therefore to do a comparison and assessment of models for DILI biomarker assessment, and to develop a conceptual basis for a translational framework between the species.

METHODS AND RESULTS

We first established which of the available pharmacokinetic models to use by identifying the most simple and identifiable model that can describe data from both human and rats. We then developed an extension of this model for how to estimate the effects of a hepatotoxic drug in rats. Finally, we illustrated how such a framework could be useful for drug dosage selection, and how it potentially can be applied in personalized treatments designed to avoid DILI.

CONCLUSION

Our analysis provides clear guidelines of which mathematical model to use for model-based assessment of biomarkers for liver function, and it also suggests a hypothetical path to a translational framework for DILI.

An orally available small molecule that targets soluble TNF to deliver anti-TNF biologic-like efficacy in rheumatoid arthritis

Tumor necrosis factor (TNF) is a pleiotropic cytokine belonging to a family of trimeric proteins with both proinflammatory and immunoregulatory functions. TNF is a key mediator in autoimmune diseases and during the last couple of decades several biologic drugs have delivered new therapeutic options for patients suffering from chronic autoimmune diseases such as rheumatoid arthritis and chronic inflammatory bowel disease. Attempts to design small molecule therapies directed to this cytokine have not led to approved products yet. Here we report the discovery and development of a potent small molecule inhibitor of TNF that was recently moved into phase 1 clinical trials. The molecule, SAR441566, stabilizes an asymmetrical form of the soluble TNF trimer, compromises downstream signaling and inhibits the functions of TNF in vitro and in vivo. With SAR441566 being studied in healthy volunteers we hope to deliver a more convenient orally bioavailable and effective treatment option for patients suffering with chronic autoimmune diseases compared to established biologic drugs targeting TNF.

Novel method for kinetic analysis applied to transport by the uniporter OCT2

The kinetics of solute transport shed light on the roles these processes play in cellular physiology, and the absolute values of the kinetic parameters that quantitatively describe transport are increasingly used to model its impact on drug clearance. However, accurate assessment of transport kinetics is challenging. Although most carrier-mediated transport is adequately described by the Michaelis-Menten equation, its use presupposes that the rates of uptake used in the analysis of maximal rates of transport (Jmax) and half-saturation constants (Kt) reflect true unidirectional rates of influx from known concentrations of substrate. Most experimental protocols estimate the initial rate of transport from net substrate accumulation determined at a single time point (typically between 0.5 and 5 min) and assume it reflects unidirectional influx. However, this approach generally results in systematic underestimates of Jmax and overestimates of Kt; the former primarily due to the unaccounted impact of efflux of accumulated substrate, and the latter due to the influence of unstirred water layers. Here, we describe the bases of these time-dependent effects and introduce a computational model that analyzes the time course of net substrate uptake at several concentrations to calculate Jmax and Kt for unidirectional influx, taking into account the influence of unstirred water layers and mediated efflux. This method was then applied to calculate the kinetics of transport of 1-methyl-4-phenylpryridinium and metformin by renal organic cation transporter 2 as expressed in cultured Chinese hamster ovary cells.NEW & NOTEWORTHY Here, we describe a mathematical model that uses the time course of net substrate uptake into cells from several increasing concentrations to calculate unique kinetic parameters [maximal rates of transport (Jmax) and half-saturation constants (Kt)] of the process. The method is the first to take into consideration the common complicating factors of unstirred layers and carrier-mediated efflux in the experimental determination of transport kinetics.

Coproporphyrin I as an Endogenous Biomarker to Detect Reduced OATP1B Activity and Shift in Elimination Route in Chronic Kidney Disease

Coproporphyrin I (CPI) is an endogenous biomarker of organic anion transporting polypeptide 1B transporter (OATP1B). CPI plasma baseline was reported to increase with severity of chronic kidney disease (CKD). Further, ratio of CPI area under the plasma concentration-time curve (AUCR) in the presence/absence of OATP1B inhibitor rifampin was higher in patients with CKD compared with healthy participants, in contrast to pitavastatin (a clinical OATP1B probe). This study investigated mechanism(s) contributing to altered CPI baseline in patients with CKD by extending a previously developed physiologically-based pharmacokinetic (PBPK) model to this patient population. CKD-related covariates were evaluated in a stepwise manner on CPI fraction unbound in plasma (f_u,p ), OATP1B-mediated hepatic uptake clearance (CL_active ), renal clearance (CL_R ), and endogenous synthesis (k_syn ). The CPI model successfully recovered increased baseline and rifampin-mediated AUCR in patients with CKD by accounting for the following disease-related changes: 13% increase in f_u,p , 29% and 39% decrease in CL_active in mild and moderate to severe CKD, respectively, decrease in CL_R proportional to decline in glomerular filtration rate, and 27% decrease in k_syn in severe CKD. Almost complete decline in CPI renal elimination in severe CKD increased its fraction transported by OATP1B, rationalizing differences in the CPI-rifampin interaction observed between healthy participants and patients with CKD. In conclusion, mechanistic modeling performed here supports CKD-related decrease in OATP1B function to inform prospective PBPK modeling of OATP1B-mediated drug-drug interaction in these patients. Monitoring of CPI allows detection of CKD-drug interaction risk for OATP1B drugs with combined hepatic and renal elimination which may be underestimated by extrapolating the interaction risk based on pitavastatin data in healthy participants.

Recent Advances in Drug Transporter Sciences: Highlights From the Year 2020

Drug Metabolism Reviews has an impressive track record of providing scientific reviews in the area of xenobiotic biotransformation over 47 years. It has consistently proved to be resourceful to many scientists from pharmaceutical industry, academia, regulatory agencies working in diverse areas including enzymology, pharmacology, pharmacokinetics and toxicology. Over the last 5 years Drug metabolism Reviews has annually published an industry commentary aimed to highlight novel insights and approaches that have made significant impacts on the field of biotransformation (led by Cyrus Khojasteh). We hope to continue this tradition by providing an overview of advances made in the field of drug transporters during 2020. The field of drug transporters is rapidly evolving as they play an essential role in drug absorption, distribution, clearance and elimination. In this review we have selected outstanding drug transporter articles that have significantly contributed to moving forward the field of transporter science with respect to translation and improved understanding of diverse aspects including uptake clearance, clinical biomarkers, induction, proteomics, emerging transporters and tissue targeting.The theme of this review consists of synopsis that summarizes each article followed by our commentary. The objective of this work is not to provide a comprehensive review but rather exemplify novel insights and state-of-the-art highlights of recent research that have advanced our understanding of drug transporters in drug disposition. We are hopeful that this effort will prove useful to the scientific community and as such request feedback, and further extend an invitation to anyone interested in contributing to future reviews.

Recent developments in in vitro and in vivo models for improved translation of preclinical pharmacokinetics and pharmacodynamics data

Improved pharmacokinetics/pharmacodynamics (PK/PD) prediction in the early stages of drug development is essential to inform lead optimization strategies and reduce attrition rates. Recently, there have been significant advancements in the development of new in vitro and in vivo strategies to better characterize pharmacokinetic properties and efficacy of drug leads. Herein, we review advances in experimental and mathematical models for clearance predictions, advancements in developing novel tools to capture slowly metabolized drugs, in vivo model developments to capture human etiology for supporting drug development, limitations and gaps in these efforts, and a perspective on the future in the field.

Interaction between Omeprazole and Gliclazide in Relation to CYP2C19 Phenotype

The antidiabetic drug gliclazide is partly metabolized by CYP2C19, the main enzyme involved in omeprazole metabolism. The aim of the study was to explore the interaction between omeprazole and gliclazide in relation to CYP2C19 phenotype using physiologically based pharmacokinetic (PBPK) modeling approach. Developed PBPK models were verified using in vivo pharmacokinetic profiles obtained from a clinical trial on omeprazole-gliclazide interaction in healthy volunteers, CYP2C19 normal/rapid/ultrarapid metabolizers (NM/RM/UM). In addition, the association of omeprazole cotreatment with gliclazide-induced hypoglycemia was explored in 267 patients with type 2 diabetes (T2D) from the GoDARTS cohort, Scotland. The PBPK simulations predicted 1.4–1.6-fold higher gliclazide area under the curve (AUC) after 5-day treatment with 20 mg omeprazole in all CYP2C19 phenotype groups except in poor metabolizers. The predicted gliclazide AUC increased 2.1 and 2.5-fold in intermediate metabolizers, and 2.6- and 3.8-fold in NM/RM/UM group, after simulated 20-day dosing with 40 mg omeprazole once and twice daily, respectively. The predicted results were corroborated by findings in patients with T2D which demonstrated 3.3-fold higher odds of severe gliclazide-induced hypoglycemia in NM/RM/UM patients concomitantly treated with omeprazole. Our results indicate that omeprazole may increase exposure to gliclazide and thus increase the risk of gliclazide-associated hypoglycemia in the majority of patients.

Comparison of Hepatic Transporter Tissue Expression in Rodents and Interspecies Hepatic OCT1 Activity

Hepatic clearance may be uptake rate limited by organic anion transporting polypeptides (OATPs) and organic cation transporter 1 (OCT1). While comparison of OATP activity has been investigated across species, little has been reported for OCT1. Additionally, while data on interspecies transporter expression in the liver exist, quantitative comparison of these transporters in multiple tissues is lacking. In the current research, the pharmacokinetics of OCT1 substrates (sumatriptan and metformin) were assessed in Oct knockout rats for comparison with previous Oct1/2-/- mice data and OCT1 pharmacogenetics in humans. Effect of OCT1 inhibitors verapamil and erlotinib on OCT1 substrate liver partitioning was also evaluated in rats. Expression of 18 transporters, including Oatps and Octs, in 9 tissues from mice and rats was quantitated using nanoLC/MS-MS, along with uptake transporters in hepatocytes from 5 species. Interspecies differences in OCT1 activity were further evaluated via uptake of OCT1 substrates in hepatocytes with corresponding in vivo liver partitioning in rodents and monkey. In Oct1-/- rats, sumatriptan hepatic clearance and liver partitioning decreased; however, metformin pharmacokinetics were unaffected. OCT1 inhibitor coadministration decreased sumatriptan liver partitioning. In rodents, Oatp expression was highest in the liver, although comparable expression of Oatps in other tissues was determined. Expression of Octs was highest in the kidney, with liver Oct1 expression comparably lower than Oatps. Liver partitioning of OCT1 substrates was lower in rodents than in monkey, in agreement with the highest OCT1 expression and uptake of OCT1 substrates in monkey hepatocytes. Species-dependent OCT1 activity requires consideration when translating preclinical data to the clinic.

A Systematic Study of the In Vitro Pharmacokinetics and Estimated Human In Vivo Clearance of Indole and Indazole-3-Carboxamide Synthetic Cannabinoid Receptor Agonists Detected on the Illicit Drug Market

In vitro pharmacokinetic studies were conducted on enantiomer pairs of twelve valinate or tert-leucinate indole and indazole-3-carboxamide synthetic cannabinoid receptor agonists (SCRAs) detected on the illicit drug market to investigate their physicochemical parameters and structure-metabolism relationships (SMRs). Experimentally derived Log D_7.4 ranged from 2.81 (AB-FUBINACA) to 4.95 (MDMB-4en-PINACA) and all SCRAs tested were highly protein bound, ranging from 88.9 ± 0.49% ((R)-4F-MDMB-BINACA) to 99.5 ± 0.08% ((S)-MDMB-FUBINACA). Most tested SCRAs were cleared rapidly in vitro in pooled human liver microsomes (pHLM) and pooled cryopreserved human hepatocytes (pHHeps). Intrinsic clearance (CL_int) ranged from 13.7 ± 4.06 ((R)-AB-FUBINACA) to 2944 ± 95.9 mL min⁻¹ kg⁻¹ ((S)-AMB-FUBINACA) in pHLM, and from 110 ± 34.5 ((S)-AB-FUBINACA) to 3216 ± 607 mL min⁻¹ kg⁻¹ ((S)-AMB-FUBINACA) in pHHeps. Predicted Human in vivo hepatic clearance (CL_H) ranged from 0.34 ± 0.09 ((S)-AB-FUBINACA) to 17.79 ± 0.20 mL min⁻¹ kg⁻¹ ((S)-5F-AMB-PINACA) in pHLM and 1.39 ± 0.27 ((S)-MDMB-FUBINACA) to 18.25 ± 0.12 mL min⁻¹ kg⁻¹ ((S)-5F-AMB-PINACA) in pHHeps. Valinate and tert-leucinate indole and indazole-3-carboxamide SCRAs are often rapidly metabolised in vitro but are highly protein bound in vivo and therefore predicted in vivo CL_H is much slower than CL_int. This is likely to give rise to longer detection windows of these substances and their metabolites in urine, possibly as a result of accumulation of parent drug in lipid-rich tissues, with redistribution into the circulatory system and subsequent metabolism.

A Critical Perspective on 3D Liver Models for Drug Metabolism and Toxicology Studies

The poor predictability of human liver toxicity is still causing high attrition rates of drug candidates in the pharmaceutical industry at the non-clinical, clinical, and post-marketing authorization stages. This is in part caused by animal models that fail to predict various human adverse drug reactions (ADRs), resulting in undetected hepatotoxicity at the non-clinical phase of drug development. In an effort to increase the prediction of human hepatotoxicity, different approaches to enhance the physiological relevance of hepatic in vitro systems are being pursued. Three-dimensional (3D) or microfluidic technologies allow to better recapitulate hepatocyte organization and cell-matrix contacts, to include additional cell types, to incorporate fluid flow and to create gradients of oxygen and nutrients, which have led to improved differentiated cell phenotype and functionality. This comprehensive review addresses the drug-induced hepatotoxicity mechanisms and the currently available 3D liver in vitro models, their characteristics, as well as their advantages and limitations for human hepatotoxicity assessment. In addition, since toxic responses are greatly dependent on the culture model, a comparative analysis of the toxicity studies performed using two-dimensional (2D) and 3D in vitro strategies with recognized hepatotoxic compounds, such as paracetamol, diclofenac, and troglitazone is performed, further highlighting the need for harmonization of the respective characterization methods. Finally, taking a step forward, we propose a roadmap for the assessment of drugs hepatotoxicity based on fully characterized fit-for-purpose in vitro models, taking advantage of the best of each model, which will ultimately contribute to more informed decision-making in the drug development and risk assessment fields.

Organic Anion-Transporting Polypeptide 1B1/1B3-Mediated Hepatic Uptake Determines the Pharmacokinetics of Large Lipophilic Acids: In Vitro-In Vivo Evaluation in Cynomolgus Monkey

It is generally presumed that uptake transport mechanisms are of limited significance in hepatic clearance for lipophilic or high passive-permeability drugs. In this study, we evaluated the mechanistic role of the hepato-selective organic anion-transporting polypeptides (OATPs) 1B1/1B3 in the pharmacokinetics of compounds representing large lipophilic acid space. Intravenous pharmacokinetics of 16 compounds with molecular mass ∼400-730 Da, logP ∼3.5-8, and acid pKa <6 were obtained in cynomolgus monkey after dosing without and with a single-dose rifampicin-OATP1B1/1B3 probe inhibitor. Rifampicin (30 mg/kg oral) significantly (P < 0.05) reduced monkey clearance and/or steady-state volume of distribution (VDss) for 15 of 16 acids evaluated. Additionally, clearance of danoprevir was reduced by about 35%, although statistical significance was not reached. A significant linear relationship was noted between the clearance ratio (i.e., ratio of control to treatment groups) and VDss ratio, suggesting hepatic uptake contributes to the systemic clearance and distribution simultaneously. In vitro transport studies using primary monkey and human hepatocytes showed uptake inhibition by rifampicin (100 µM) for compounds with logP ≤6.5 but not for the very lipophilic acids (logP > 6.5), which generally showed high nonspecific binding in hepatocyte incubations. In vitro uptake clearance and fraction transported by OATP1B1/1B3 (f_t,OATP1B) were found to be similar in monkey and human hepatocytes. Finally, for compounds with logP ≤6.5, good agreement was noted between in vitro f_t,OATP1B and clearance ratio (as well as VDss ratio) in cynomolgus monkey. In conclusion, this study provides mechanistic evidence for the pivotal role of OATP1B-mediated hepatic uptake in the pharmacokinetics across a wide, large lipophilic acid space. SIGNIFICANCE STATEMENT: This study provides mechanistic insight into the pharmacokinetics of a broad range of large lipophilic acids. Organic anion-transporting polypeptides 1B1/1B3-mediated hepatic uptake is of key importance in the pharmacokinetics and drug-drug interactions of almost all drugs and new molecular entities in this space. Diligent in vitro and in vivo transport characterization is needed to avoid the false negatives often noted because of general limitations in the in vitro assays while handling compounds with such physicochemical attributes.

Hepatic transporter-mediated pharmacokinetic drug-drug interactions: Recent studies and regulatory recommendations

Transporter-mediated drug-drug interactions are one of the major mechanisms in pharmacokinetic-based drug interactions and correspondingly affecting drugs' safety and efficacy. Regulatory bodies underlined the importance of the evaluation of transporter-mediated interactions as a part of the drug development process. The liver is responsible for the elimination of a wide range of endogenous and exogenous compounds via metabolism and biliary excretion. Therefore, hepatic uptake transporters, expressed on the sinusoidal membranes of hepatocytes, and efflux transporters mediating the transport from hepatocytes to the bile are determinant factors for pharmacokinetics of drugs, and hence, drug-drug interactions. In parallel with the growing research interest in this area, regulatory guidances have been updated with detailed assay models and criteria. According to well-established preclinical results, observed or expected hepatic transporter-mediated drug-drug interactions can be taken into account for clinical studies. In this paper, various methods including in vitro, in situ, in vivo, in silico approaches, and combinational concepts and several clinical studies on the assessment of transporter-mediated drug-drug interactions were reviewed. Informative and effective evaluation by preclinical tools together with the integration of pharmacokinetic modeling and simulation can reduce unexpected clinical outcomes and enhance the success rate in drug development.

Evaluation of the Utility of PXB Chimeric Mice for Predicting Human Liver Partitioning of Hepatic Organic Anion-Transporting Polypeptide Transporter Substrates

The ability to predict human liver-to-plasma unbound partition coefficient (K_puu) is important to estimate unbound liver concentration for drugs that are substrates of hepatic organic anion-transporting peptide (OATP) transporters with asymmetric distribution into the liver relative to plasma. Herein, we explored the utility of PXB chimeric mice with humanized liver that are highly repopulated with human hepatocytes to predict human hepatic disposition of OATP substrates, including rosuvastatin, pravastatin, pitavastatin, valsartan, and repaglinide. In vitro total uptake clearance and transporter-mediated active uptake clearance in C57 mouse hepatocytes were greater than in PXB chimeric mouse hepatocytes for rosuvastatin, pravastatin, pitavastatin, and valsartan. Consistent with in vitro uptake data, enhanced hepatic uptake and resulting total systemic clearance were observed with the above four compounds in severely compromised immune-deficient (SCID) control mice compared with the PXB chimeric mice, which suggest that mouse has a stronger transporter-mediated hepatic uptake than human. In vivo liver-to-plasma K_puu from PXB chimeric and SCID control mice were also compared, and rosuvastatin and pravastatin K_puu in SCID mice were more than 10-fold higher than that in PXB chimeric mice, whereas pitavastatin, valsartan, and repaglinide K_puu in SCID mice were comparable with K_puu in PXB chimeric mice. Finally, PXB chimeric mouse liver-to-plasma K_puu values were compared with the reported human K_puu, and a good correlation was observed as the PXB K_puu vales were within 3-fold of human K_puu Our results indicate that PXB mice could be a useful tool to delineate hepatic uptake and enable prediction of human liver-to-plasma K_puu of hepatic uptake transporter substrates. SIGNIFICANCE STATEMENT: We evaluated PXB mouse with humanized liver for its ability to predict human liver disposition of five organic anion-transporting polypeptide transporter substrates. Both in vitro and in vivo data suggest that mouse liver has a stronger transporter-mediated hepatic uptake than the humanized liver in PXB mouse. More importantly, PXB liver-to-plasma unbound partition coefficient (K_puu) values were compared with the reported human K_puu, and a good correlation was observed. PXB mice could be a useful tool to project human liver-to-plasma K_puu of hepatic uptake transporter substrates.

Evaluation of Hepatic Uptake of OATP1B Substrates by Short Term-Cultured Plated Human Hepatocytes: Comparison With Isolated Suspended Hepatocytes

Hepatic uptake clearance has been measured in suspended human hepatocytes (SHH). Plated human hepatocytes (PHH) after short-term culturing are increasingly employed to study hepatic transport driven mainly by its higher throughput. To know pros/cons of both systems, the hepatic uptake clearances of several organic anion transporting polypeptide 1B substrates were compared between PHH and SHH by determining the initial uptake velocities or through dynamic model-based analyses. For cerivastatin, pitavastatin and rosuvastatin, initial uptake clearances (PS_inf) obtained using PHH were comparable to those using SHH, while cell-to-medium concentration (C/M) ratios were 2.7- to 5.4-fold higher. For pravastatin and dehydropravastatin, hydrophilic compounds with low uptake/cellular binding, their PS_inf and C/M ratio in PHH were 1.8- to 3.2-fold lower than those in SHH. These hydrophilic substrates are more prone to wash-off during the uptake study using PHH, which may explain the apparently lower uptake than SHH. The C/M ratios obtained using PHH were stable over an extended time, making PHH suitable to estimate the C/M ratios and hepatocyte-to-medium unbound concentration ratios (K_p,uu). In conclusion, PHH is useful in evaluating hepatic uptake/efflux clearances and K_p,uu of OATP1B substrates in a high-throughput manner, however, a caution is warranted for hydrophilic drugs with low uptake/cellular binding.

Predicting the Human Hepatic Clearance of Acidic and Zwitterionic Drugs

Prospective predictions of human hepatic clearance for anionic/zwitterionic compounds, which are oftentimes subjected to transporter-mediated uptake, are challenging in drug discovery. We evaluated the utility of preclinical species, rats and cynomolgus monkeys [nonhuman primates (NHPs)], to predict the human hepatic clearance using a diverse set of acidic/zwitterionic drugs. Preclinical clearance data were generated following intravenous dosing in rats/NHPs and compared to the human clearance data (n = 18/27). Single-species scaling of NHP clearance with an allometric exponent of 0.50 allowed for good prediction of human clearance (fold error ∼2.1, bias ∼1.0), with ∼86% predictions within 3-fold. In comparison, rats underpredicted the clearance of lipophilic acids, while overprediction was noted for hydrophilic acids. Finally, an in vitro clearance assay based on human hepatocytes, which is routinely used in discovery setting, markedly underpredicted human clearance (bias ∼0.12). Collectively, this study provides insights into the usefulness of the preclinical models in enabling pharmacokinetic optimization for acid/zwitterionic drug candidates.

Physiologically Based Pharmacokinetic Models for Prediction of Complex CYP2C8 and OATP1B1 (SLCO1B1) Drug-Drug-Gene Interactions: A Modeling Network of Gemfibrozil, Repaglinide, Pioglitazone, Rifampicin, Clarithromycin and Itraconazole

Background

Drug–drug interactions (DDIs) and drug–gene interactions (DGIs) pose a serious health risk that can be avoided by dose adaptation. These interactions are investigated in strictly controlled setups, quantifying the effect of one perpetrator drug or polymorphism at a time, but in real life patients frequently take more than two medications and are very heterogenous regarding their genetic background.

Objectives

The first objective of this study was to provide whole-body physiologically based pharmacokinetic (PBPK) models of important cytochrome P450 (CYP) 2C8 perpetrator and victim drugs, built and evaluated for DDI and DGI studies. The second objective was to apply these models to describe complex interactions with more than two interacting partners.

Methods

PBPK models of the CYP2C8 and organic-anion-transporting polypeptide (OATP) 1B1 perpetrator drug gemfibrozil (parent–metabolite model) and the CYP2C8 victim drugs repaglinide (also an OATP1B1 substrate) and pioglitazone were developed using a total of 103 clinical studies. For evaluation, these models were applied to predict 34 different DDI studies, establishing a CYP2C8 and OATP1B1 PBPK DDI modeling network.

Results

The newly developed models show a good performance, accurately describing plasma concentration–time profiles, area under the plasma concentration–time curve (AUC) and maximum plasma concentration (C_max) values, DDI studies as well as DGI studies. All 34 of the modeled DDI AUC ratios (AUC during DDI/AUC control) and DDI C_max ratios (C_max during DDI/C_max control) are within twofold of the observed values.

Conclusions

Whole-body PBPK models of gemfibrozil, repaglinide, and pioglitazone have been built and qualified for DDI and DGI prediction. PBPK modeling is applicable to investigate complex interactions between multiple drugs and genetic polymorphisms.

Electronic supplementary material

The online version of this article (10.1007/s40262-019-00777-x) contains supplementary material, which is available to authorized users.

Predictability of human pharmacokinetics of drugs that undergo hepatic organic anion transporting polypeptide (OATP)-mediated transport using single-species allometric scaling in chimeric mice with humanized liver: integration with hepatic drug metabolism

We previously reported a prediction method for human pharmacokinetics (PK) using single species allometric scaling (SSS) and the complex Dedrick plot in chimeric mice with humanized liver to predict the total clearance (CL_t), distribution volumes in steady state (Vd_ss) and plasma concentration-time profiles of several drugs metabolized by cytochrome P450 (P450) and non-P450 enzymes. In the present study, we examined eight compounds (bosentan, cerivastatin, fluvastatin, pitavastatin, pravastatin, repaglinide, rosuvastatin, valsartan) as typical organic anion transporting polypeptide (OATP) substrates and six compounds metabolized by P450 and non-P450 enzymes to evaluate the predictability of CL_t, Vd_ss and plasma concentration-time profiles after intravenous administration to chimeric mice. The predicted CL_t and Vd_ss of drugs that undergo OATP-mediated uptake and P450/non-P450-mediated metabolism reflected the observed data from humans within a threefold error range. We also examined the possibility of predicting plasma concentration-time profiles of drugs that undergo OATP-mediated uptake using the complex Dedrick plot in chimeric mice. Most profiles could be superimposed with observed profiles from humans within a two- to threefold error range. PK prediction using SSS and the complex Dedrick plot in chimeric mice can be useful for evaluating drugs that undergo both OATP-mediated uptake and P450/non-P450-mediated metabolism.

Prediction of Clinical Transporter-Mediated Drug-Drug Interactions via Comeasurement of Pitavastatin and Eltrombopag in Human Hepatocyte Models

A structurally identifiable micro‐rate constant mechanistic model was used to describe the interaction between pitavastatin and eltrombopag, with improved goodness‐of‐fit values through comeasurement of pitavastatin and eltrombopag. Transporter association and dissociation rate constants and passive rates out of the cell were similar between pitavastatin and eltrombopag. Translocation into the cell through transporter‐mediated uptake was six times greater for pitavastatin, leading to pronounced inhibition of pitavastatin uptake by eltrombopag. The passive rate into the cell was 91 times smaller for pitavastatin compared with eltrombopag. A semimechanistic physiologically‐based pharmacokinetic (PBPK) model was developed to evaluate the potential for clinical drug–drug interactions (DDIs). The PBPK model predicted a twofold increase in the pitavastatin peak blood concentration and area under the concentration‐time curve in the presence of eltrombopag in simulated healthy volunteers. The use of structural identifiability supporting experimental design combined with robust micro‐rate constant parameter estimates and a semimechanistic PBPK model gave more informed predictions of transporter‐mediated DDIs.

Multiscale modelling of drug transport and metabolism in liver spheroids

In early preclinical drug development, potential candidates are tested in the laboratory using isolated cells. These in vitro experiments traditionally involve cells cultured in a two-dimensional monolayer environment. However, cells cultured in three-dimensional spheroid systems have been shown to more closely resemble the functionality and morphology of cells in vivo. While the increasing usage of hepatic spheroid cultures allows for more relevant experimentation in a more realistic biological environment, the underlying physical processes of drug transport, uptake and metabolism contributing to the spatial distribution of drugs in these spheroids remain poorly understood. The development of a multiscale mathematical modelling framework describing the spatio-temporal dynamics of drugs in multicellular environments enables mechanistic insight into the behaviour of these systems. Here, our analysis of cell membrane permeation and porosity throughout the spheroid reveals the impact of these properties on drug penetration, with maximal disparity between zonal metabolism rates occurring for drugs of intermediate lipophilicity. Our research shows how mathematical models can be used to simulate the activity and transport of drugs in hepatic spheroids and in principle any organoid, with the ultimate aim of better informing experimentalists on how to regulate dosing and culture conditions to more effectively optimize drug delivery.

Use of imaging to assess the activity of hepatic transporters

Introduction: Membrane transporters of the SLC and ABC families are abundantly expressed in the liver, where they control the transfer of drugs/drug metabolites across the sinusoidal and canalicular hepatocyte membranes and play a pivotal role in hepatic drug clearance. Noninvasive imaging methods, such as PET, SPECT or MRI, allow for measuring the activity of hepatic transporters in vivo, provided that suitable transporter imaging probes are available.Areas covered: We give an overview of the working principles of imaging-based assessment of hepatic transporter activity. We discuss different currently available PET/SPECT radiotracers and MRI contrast agents and their applications to measure hepatic transporter activity in health and disease. We cover mathematical modeling approaches to obtain quantitative parameters of transporter activity and provide a critical assessment of methodological limitations and challenges associated with this approach.Expert opinion: PET in combination with pharmacokinetic modeling can be potentially applied in drug development to study the distribution of new drug candidates to the liver and their clearance mechanisms. This approach bears potential to mechanistically assess transporter-mediated drug-drug interactions, to assess the influence of disease on hepatic drug disposition and to validate and refine currently available in vitro-in vivo extrapolation methods to predict hepatic clearance of drugs.

A hybrid model to evaluate the impact of active uptake transport on hepatic distribution of atorvastatin in rats

1. Mathematical modeling remains a useful tool to study the impact of transporters on overall and intracellular drug disposition. The impact of organic anion transporter protein mediated uptake on atorvastatin systemic and intracellular pharmacokinetics, specifically distribution volume, was studied in rats with mathematical modeling and conducting in vivo pharmacokinetic studies for atorvastatin in presence and absence of rifampicin. A previously developed 5-compartment explicit membrane model for the liver was combined with a compartmental model to develop a semi-physiological hybrid model for atorvastatin disposition. 2. Rifampicin treatment resulted in a decrease in systemic clearance and steady-state distribution volume, and an increase in half-life of atorvastatin. The hybrid model predicted higher unbound intracellular liver atorvastatin concentrations than unbound plasma concentrations in both rifampicin treated and untreated groups, indicating involvement of uptake transporters. The intracellular unbound concentrations during the distributive phase were unaffected by rifampicin. The dependence of clearance on blood flow as well as hepatic uptake for atorvastatin (a moderate-to-high extraction ratio drug) can explain this lack of change in intracellular concentrations if there is incomplete inhibition of transport at the tested rifampicin dose. 3. The hybrid model successfully allowed the evaluation of effect of active uptake on intracellular and plasma atorvastatin disposition.

Quantitative Translation of Microfluidic Transporter in Vitro Data to in Vivo Reveals Impaired Albumin-Facilitated Indoxyl Sulfate Secretion in Chronic Kidney Disease

Indoxyl sulfate (IxS), a highly albumin-bound uremic solute, accumulates in chronic kidney disease (CKD) due to reduced renal clearance. This study was designed to specifically investigate the role of human serum albumin (HSA) in IxS renal secretion via organic anion transporter 1 (OAT1) in a microfluidic system and subsequently apply quantitative translation of in vitro data to predict extent of change in IxS renal clearance in CKD stage IV relative to healthy. Conditionally immortalized human proximal tubule epithelial cells overexpressing OAT1 were incubated with IxS (5-200 μM) in the HSA-free medium or in the presence of either HSA or CKD-modified HSA. IxS uptake in the presence of HSA resulted in more than 20-fold decrease in OAT1 affinity (K_m,u) and 37-fold greater in vitro unbound intrinsic clearance (CL_int,u) versus albumin-free condition. In the presence of CKD-modified albumin, K_m,u increased four-fold and IxS CL_int,u decreased almost seven-fold relative to HSA. Fold-change in parameters exceeded differences in IxS binding between albumin conditions, indicating additional mechanism and facilitating role of albumin in IxS OAT1-mediated uptake. Quantitative translation of IxS in vitro OAT1-mediated CL_int,u predicted a 60% decrease in IxS renal elimination as a result of CKD, in agreement with the observed data (80%). The findings of the current study emphasize the role of albumin in IxS transport via OAT1 and explored the impact of modifications in albumin on renal excretion via active secretion in CKD. For the first time, this study performed quantitative translation of transporter kinetic data generated in a novel microfluidic in vitro system to a clinically relevant setting. Knowledge gaps and future directions in quantitative translation of renal drug disposition from microphysiological systems are discussed.

A mechanistic modelling approach for the determination of the mechanisms of inhibition by cyclosporine on the uptake and metabolism of atorvastatin in rat hepatocytes using a high throughput uptake method

Determine the inhibition mechanism through which cyclosporine inhibits the uptake and metabolism of atorvastatin in fresh rat hepatocytes using mechanistic models applied to data generated using a high throughput oil spin method.Atorvastatin was incubated in fresh rat hepatocytes (0.05-150 nmol/ml) with or without 20 min pre-incubation with 10 nmol/ml cyclosporine and sampled over 0.25-60 min using a high throughput oil spin method. Micro-rate constant and macro-rate constant mechanistic models were ranked based on goodness of fit values.The best fitting model to the data was a micro-rate constant mechanistic model including non-competitive inhibition of uptake and competitive inhibition of metabolism by cyclosporine (Model 2). The association rate constant for atorvastatin was 150-fold greater than the dissociation rate constant and 10-fold greater than the translocation into the cell. The association and dissociation rate constants for cyclosporine were 7-fold smaller and 10-fold greater, respectively, than atorvastatin. The simulated atorvastatin-transporter-cyclosporine complex derived using the micro-rate constant parameter estimates increased in line with the incubation concentration of atorvastatin.The increased amount of data generated with the high throughput oil spin method, combined with a micro-rate constant mechanistic model helps to explain the inhibition of uptake by cyclosporine following pre-incubation.

The clinical significance of statins-macrolides interaction: comprehensive review of in vivo studies, case reports, and population studies

The objectives of this article were to review the mechanism and clinical significance of statins-macrolides interaction, determine which combination has the highest risk for the interaction, and identify key patients' risk factors for the interaction in relation to the development of muscle toxicity. A literature review was conducted in PubMed and Embase (1946 to December 2018) using combined terms: statins - as group and individual agents, macrolides - as group and individual agents, drug interaction, muscle toxicity, rhabdomyolysis, CYP3A4 inhibitors, and OAT1B inhibitors, with forward and backward citation tracking. Relevant English language in vivo studies in healthy volunteers, case reports, and population studies were included. The interaction between statins and macrolides depends on the type of statin and macrolide used. The mechanism of the interaction is due to macrolides' inhibition of CYP3A4 isoenzyme and OAT1B transporter causing increased exposure to statins. The correlation of this increased statin's exposure to the development of muscle toxicity could not be established, unless the patient had other risk factors such as advanced age, cardiovascular diseases, renal impairment, diabetes, and the concomitant use of other CYP3A4 inhibitors. Simvastatin, lovastatin, and to lesser extent atorvastatin are the statins most affected by this interaction. Rosuvastatin, fluvastatin, and pravastatin are not significantly affected by this interaction. Telithromycin, clarithromycin, and erythromycin are the most "offending" macrolides, while azithromycin appears to be safe to use with statins. This review presented a clear description of the clinical significance of this interaction in real practice. Also, it provided health care professionals with clear suggestions and recommendations on how to overcome this interaction. In conclusion, understanding the different characteristics of each statin and macrolide, as well as key patients' risk factors, will enable health care providers to utilize both groups effectively without compromising patient safety.

Membrane transporter data to support kinetically-informed chemical risk assessment using non-animal methods: Scientific and regulatory perspectives

Humans are continuously exposed to low levels of thousands of industrial chemicals, most of which are poorly characterised in terms of their potential toxicity. The new paradigm in chemical risk assessment (CRA) aims to rely on animal-free testing, with kinetics being a key determinant of toxicity when moving from traditional animal studies to integrated in vitro-in silico approaches. In a kinetically informed CRA, membrane transporters, which have been intensively studied during drug development, are an essential piece of information. However, how existing knowledge on transporters gained in the drug field can be applied to CRA is not yet fully understood. This review outlines the opportunities, challenges and existing tools for investigating chemical-transporter interactions in kinetically informed CRA without animal studies. Various environmental chemicals acting as substrates, inhibitors or modulators of transporter activity or expression have been shown to impact TK, just as drugs do. However, because pollutant concentrations are often lower in humans than drugs and because exposure levels and internal chemical doses are not usually known in contrast to drugs, new approaches are required to translate transporter data and reasoning from the drug sector to CRA. Here, the generation of in vitro chemical-transporter interaction data and the development of transporter databases and classification systems trained on chemical datasets (and not only drugs) are proposed. Furtheremore, improving the use of human biomonitoring data to evaluate the in vitro-in silico transporter-related predicted values and developing means to assess uncertainties could also lead to increase confidence of scientists and regulators in animal-free CRA. Finally, a systematic characterisation of the transportome (quantitative monitoring of transporter abundance, activity and maintenance over time) would reinforce confidence in the use of experimental transporter/barrier systems as well as in established cell-based toxicological assays currently used for CRA.

Hepatic Organic Anion Transporting Polypeptide-Mediated Clearance in the Beagle Dog: Assessing In Vitro-In Vivo Relationships and Applying Cross-Species Empirical Scaling Factors to Improve Prediction of Human Clearance

In the present study, the beagle dog was evaluated as a preclinical model to investigate organic anion transporting polypeptide (OATP)-mediated hepatic clearance. In vitro studies were performed with nine OATP substrates in three lots of plated male dog hepatocytes ± OATP inhibitor cocktail to determine total uptake clearance (CL_uptake) and total and unbound cell-to-medium concentration ratio (Kp_uu). In vivo intrinsic hepatic clearances (CL_int,H) were determined following intravenous drug administration (0.1 mg/kg) in male beagle dogs. The in vitro parameters were compared with those previously reported in plated human, monkey, and rat hepatocytes; the ability of cross-species scaling factors to improve prediction of human in vivo clearance was assessed. CL_uptake in dog hepatocytes ranged from 9.4 to 135 µl/min/10⁶ cells for fexofenadine and telmisartan, respectively. Active process contributed >75% to CL_uptake for 5/9 drugs. Rosuvastatin and valsartan showed Kp_uu > 10, whereas cerivastatin, pitavastatin, repaglinide, and telmisartan had Kp_uu < 5. The extent of hepatocellular binding in dog was consistent with other preclinical species and humans. The bias (2.73-fold) obtained from comparison of predicted versus in vivo dog CL_int,H was applied as an average empirical scaling factor (ESF_av) for in vitro-in vivo extrapolation of human CL_int,H The ESF_av based on dog reduced underprediction of human CL_int,H for the same data set (geometric mean fold error = 2.1), highlighting its utility as a preclinical model to investigate OATP-mediated uptake. The ESF_av from all preclinical species resulted in comparable improvement of human clearance prediction, in contrast to drug-specific empirical scalars, rationalized by species differences in expression and/or relative contribution of particular transporters to drug hepatic uptake.

Advancing Predictions of Tissue and Intracellular Drug Concentrations Using In Vitro, Imaging and Physiologically Based Pharmacokinetic Modeling Approaches

This white paper examines recent progress, applications, and challenges in predicting unbound and total tissue and intra/subcellular drug concentrations using in vitro and preclinical models, imaging techniques, and physiologically based pharmacokinetic (PBPK) modeling. Published examples, regulatory submissions, and case studies illustrate the application of different types of data in drug development to support modeling and decision making for compounds with transporter-mediated disposition, and likely disconnects between tissue and systemic drug exposure. The goals of this article are to illustrate current best practices and outline practical strategies for selecting appropriate in vitro and in vivo experimental methods to estimate or predict tissue and plasma concentrations, and to use these data in the application of PBPK modeling for human pharmacokinetic (PK), efficacy, and safety assessment in drug development.

Comparison of uptake transporter functions in hepatocytes in different species to determine the optimal model for evaluating drug transporter activities in humans

A thorough understanding of species-dependent differences in hepatic uptake transporters is critical for predicting human pharmacokinetics (PKs) from preclinical data. In this study, the activities of organic anion transporting polypeptide (OATP/Oatp), organic cation transporter 1 (OCT1/Oct1), and sodium-taurocholate cotransporting polypeptide (NTCP/Ntcp) in cultured rat, dog, monkey and human hepatocytes were compared. The activities of hepatic uptake transporters were evaluated with respect to culture duration, substrate and species-dependent differences in hepatocytes. Longer culture duration reduced hepatic uptake transporter activities across species except for Oatp and Ntcp in rats. Comparable apparent Michaelis-Menten constant (K_m,app) values in hepatocytes were observed across species for atorvastatin, estradiol-17β-glucuronide and metformin. The K_m,app values for rosuvastatin and taurocholate were significantly different across species. Rat hepatocytes exhibited the highest Oatp percentage of uptake transporter-mediated permeation clearance (PS_inf,act) while no difference in %PS_inf,act of probe substrates were observed across species. The in vitro hepatocyte inhibition data in rats, monkeys and humans provided reasonable predictions of in vivo drug-drug interaction (DDIs) between atorvastatin/rosuvastatin and rifampin. These findings suggested that using human hepatocytes with a short culture time is the most robust preclinical model for predicting DDIs for compounds exhibiting active hepatic uptake in humans.

Cryopreserved Human Intestinal Mucosal Epithelium: A Novel In Vitro Experimental System for the Evaluation of Enteric Drug Metabolism, Cytochrome P450 Induction, and Enterotoxicity

We report here a novel in vitro enteric experimental system, cryopreserved human intestinal mucosa (CHIM), for the evaluation of enteric drug metabolism, drug-drug interaction, drug toxicity, and pharmacology. CHIM was isolated from the small intestines of four human donors. The small intestines were first dissected into the duodenum, jejunum, and ileum, followed by collagenase digestion of the intestinal lumen. The isolated mucosa was gently homogenized to yield multiple cellular fragments, which were then cryopreserved in a programmable liquid cell freezer and stored in liquid nitrogen. After thawing and recovery, CHIM retained robust cytochrome P450 (P450) and non-P450 drug-metabolizing enzyme activities and demonstrated dose-dependent induction of transcription of CYP24A1 (approximately 300-fold) and CYP3A4 (approximately 3-fold) by vitamin D₃ as well as induction of CYP3A4 (approximately 3-fold) by rifampin after 24 hours of treatment. Dose-dependent decreases in cell viability quantified by cellular ATP content were observed for naproxen and acetaminophen, with higher enterotoxicity observed for naproxen, consistent with that observed in humans in vivo. These results suggest that CHIM may be a useful in vitro experimental model for the evaluation of enteric drug properties, including drug metabolism, drug-drug interactions, and drug toxicity.

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.

Physiologically Based Pharmacokinetic Modeling of Bosentan Identifies the Saturable Hepatic Uptake As a Major Contributor to Its Nonlinear Pharmacokinetics

Bosentan is a substrate of hepatic uptake transporter organic anion-transporting polypeptides (OATPs), and undergoes extensive hepatic metabolism by cytochrome P450 (P450), namely, CYP3A4 and CYP2C9. Several clinical investigations have reported a nonlinear relationship between bosentan doses and its systemic exposure, which likely involves the saturation of OATP-mediated uptake, P450-mediated metabolism, or both in the liver. Yet, the underlying causes for the nonlinear bosentan pharmacokinetics are not fully delineated. To address this, we performed physiologically based pharmacokinetic (PBPK) modeling analyses for bosentan after its intravenous administration at different doses. As a bottom-up approach, PBPK modeling analyses were performed using in vitro kinetic parameters, other relevant parameters, and scaling factors. As top-down approaches, three different types of PBPK models that incorporate the saturation of hepatic uptake, metabolism, or both were compared. The prediction from the bottom-up approach (models 1 and 2) yielded blood bosentan concentration-time profiles and their systemic clearance values that were not in good agreement with the clinically observed data. From top-down approaches (models 3, 4, 5-1, and 5-2), the prediction accuracy was best only with the incorporation of the saturable hepatic uptake for bosentan. Taken together, the PBPK models for bosentan were successfully established, and the comparison of different PBPK models identified the saturation of the hepatic uptake process as a major contributing factor for the nonlinear pharmacokinetics of bosentan.

Prediction of the Pharmacokinetics of Pravastatin as an OATP Substrate Using Plateable Human Hepatocytes With Human Plasma Data and PBPK Modeling

Plateable human hepatocytes with human plasma were utilized to generate the uptake transporter kinetic data for pravastatin, an organic anion-transporting polypeptide (OATP) transporter substrate. The active hepatic uptake of pravastatin was determined with a J_max value of 134.4 pmol/min/million cells and K_m of 76.77 µM in plateable human hepatocytes with human plasma. The physiologically-based pharmacokinetic (PBPK) model with incorporation of these in vitro kinetic data successfully simulated the i.v. pharmacokinetic profile of pravastatin without applying scaling factor (the mean predicted area under the curve (AUC) is within 1.5-fold of the observed). Furthermore, the PBPK model also adequately described the oral plasma concentration-time profiles of pravastatin at different dose levels. The current investigation demonstrates an approach allowing us to build upon the translation of in vitro OATP uptake transporter data to in vivo, with a hope of utilizing the in vitro data for the prospective human pharmacokinetic (PK) prediction.

Successful Prediction of In Vivo Hepatobiliary Clearances and Hepatic Concentrations of Rosuvastatin Using Sandwich-Cultured Rat Hepatocytes, Transporter-Expressing Cell Lines, and Quantitative Proteomics

We determined whether in vivo transporter-mediated hepatobiliary clearance (CL) and hepatic concentrations of rosuvastatin (RSV) in the rat could be predicted by transport activity in sandwich-cultured rat hepatocytes (SCRHs) and/or transporter-expressing cell lines scaled by differences in transporter protein expression between SCRHs, cell lines, and rat liver. The predicted hepatobiliary CLs and hepatic concentrations of RSV were compared with our previously published positron emission tomography imaging data. Sinusoidal uptake CL ([Formula: see text]) and efflux (canalicular and sinusoidal) CLs of [³H]-RSV in SCRHs were evaluated in the presence and absence of Ca²⁺ and in the absence and presence of 1 mM unlabeled RSV (to estimate passive diffusion CL). [Formula: see text] of RSV into cells expressing organic anion transporting polypeptide (Oatp) 1a1, 1a4, and 1b2 was also determined. Protein expression of Oatps in SCRHs and Oatp-expressing cells was quantified by liquid chromatography tandem mass spectrometry. SCRHs well predicted the in vivo RSV sinusoidal and canalicular efflux CLs but significantly underestimated in vivo [Formula: see text]. Oatp expression in SCRHs was significantly lower than that in the rat liver. [Formula: see text], based on RSV [Formula: see text] into Oatp-expressing cells (active transport) plus passive diffusion CL in SCRHs, scaled by the difference in protein expression in Oatp cells versus SCRH versus rat liver, was within 2-fold of that observed in SCRHs or in vivo. In vivo hepatic RSV concentrations were well predicted by Oatp-expressing cells after correcting [Formula: see text] for Oatp protein expression. This is the first demonstration of the successful prediction of in vivo hepatobiliary CLs and hepatic concentrations of RSV using transporter-expressing cells and SCRHs.

Physiologically Based Pharmacokinetic Modeling of Transporter-Mediated Hepatic Clearance and Liver Partitioning of OATP and OCT Substrates in Cynomolgus Monkeys

In the present investigations, we evaluate in vitro hepatocyte uptake and partitioning for the prediction of in vivo clearance and liver partitioning. Monkeys were intravenously co-dosed with rosuvastatin and bosentan, substrates of the organic anion transporting polypeptides (OATPs), and metformin, a substrate of organic cation transporter 1 (OCT1). Serial plasma and liver samples were collected over time. Liver and plasma unbound fraction was determined using equilibrium dialysis. In vivo unbound partitioning (Kp_u,u) for rosuvastatin, bosentan, and metformin, calculated from total concentrations in the liver and plasma, were 243, 553, and 15, respectively. A physiologically based pharmacokinetic monkey model that incorporates active and passive hepatic uptake was developed to fit plasma and liver concentrations. In addition, a two-compartment model was used to fit in vitro hepatic uptake curves in suspended monkey hepatocyte to determine active uptake, passive diffusion, and intracellular unbound fraction parameters. At steady-state in the model, in vitro Kp_u,u was determined. The results demonstrated that in vitro values under-predicted in vivo active uptake for rosuvastatin, bosentan, and metformin by 6.7-, 28-, and 1.5-fold, respectively, while passive diffusion was over-predicted. In vivo Kp_u,u values were under-predicted from in vitro data by 30-, 79-, and 3-fold. In conclusion, active uptake and liver partitioning in monkeys for OATP substrates were greatly under-predicted from in vitro hepatocyte uptake, while OCT-mediated uptake and partitioning scaled reasonably well from in vitro, demonstrating substrate- and transporter-dependent scaling factors. The combination of in vitro experimental and modeling approaches proved useful for assessing prediction of in vivo intracellular partitioning.

Bosentan and Rifampin Interactions Modulate Influx Transporter and Cytochrome P450 Expression and Activities in Primary Human Hepatocytes

The incidence of polypharmacy-which can result in drug-drug interactions-has increased in recent years. Drug-metabolizing enzymes and drug transporters are important polypharmacy modulators. In this study, the effects of bosentan and rifampin on the expression and activities of organic anion-transporting peptide (OATP) and cytochrome P450 (CYP450) 2C9 and CYP3A4 were investigated in vitro. HEK293 cells and primary human hepatocytes overexpressing the target genes were treated with bosentan and various concentrations of rifampin, which decreased the uptake activities of OATP transporters in a dose-dependent manner. In primary human hepatocytes, CYP2C9 and CYP3A4 gene expression and activities decreased upon treatment with 20 μM bosentan+200 μM rifampin. Rifampin also reduced gene expression of OATP1B1, OATP1B3, and OATP2B1 transporter, and inhibited bosentan influx in human hepatocytes at increasing concentrations. These results confirm rifampin- and bosentan-induced interactions between OATP transporters and CYP450.

Physiologically Based Pharmacokinetic Modeling of Drug Transporters to Facilitate Individualized Dose Prediction

Physiologically based pharmacokinetic modeling is a commonly used strategy in the drug development and regulatory submissions. This commentary provides a critical overview of the current status of physiologically based pharmacokinetic methodologies to predict transporter-mediated pharmacokinetics, in addition to the impact of disease and genetics with respect to local and systemic concentration.

Investigating Transporter-Mediated Drug-Drug Interactions Using a Physiologically Based Pharmacokinetic Model of Rosuvastatin

Rosuvastatin is a frequently used probe in transporter-mediated drug-drug interaction (DDI) studies. This report describes the development of a physiologically based pharmacokinetic (PBPK) model of rosuvastatin for prediction of pharmacokinetic (PK) DDIs. The rosuvastatin model predicted the observed single (i.v. and oral) and multiple dose PK profiles, as well as the impact of coadministration with transporter inhibitors. The predicted effects of rifampin and cyclosporine (6.58-fold and 5.07-fold increase in rosuvastatin area under the curve (AUC), respectively) were mediated primarily via inhibition of hepatic organic anion-transporting polypeptide (OATP)1B1 (Inhibition constant (K_i ) ∼1.1 and 0.014 µM, respectively) and OATP1B3 (K_i ∼0.3 and 0.007 µM, respectively), with cyclosporine also inhibiting intestinal breast cancer resistance protein (BCRP; K_i ∼0.07 µM). The predicted effects of gemfibrozil and its metabolite were moderate (1.88-fold increase in rosuvastatin AUC) and mediated primarily via inhibition of hepatic OATP1B1 and renal organic cation transporter 3. This model of rosuvastatin will be useful in prospectively predicting transporter-mediated DDIs with novel pharmaceutical agents in development.