Elucidating nonlinear pharmacokinetics of telmisartan: Integration of target-mediated drug disposition and OATP1B3-mediated hepatic uptake in a physiologically based model.

Telmisartan, a selective inhibitor of angiotensin II receptor type 1 (AT1), demonstrates nonlinear pharmacokinetics (PK) when orally administered in ascending doses to healthy volunteers, but the underlying mechanisms remain unclear. This study presents a physiologically based pharmacokinetic model integrated with target-mediated drug disposition (TMDD-PBPK model) to explore the mechanism of its nonlinear PK. We employed the Cluster-Gauss Newton method for top-down analysis, estimating the in vivo Km,OATP1B3 (Michaelis-Menten constant for telmisartan hepatic uptake via Organic Anion Transporting Polypeptide 1B3) to be 2.0-5.7 nM. This range is significantly lower than the reported in vitro value of 810 nM, obtained in 0.3% human serum albumin (HSA) conditions. Further validation was achieved through in vitro assessment in plated human hepatocytes with 4.5% HSA, showing a Km of 4.5 nM. These results underscore the importance of albumin-mediated uptake effect for the hepatic uptake of telmisartan. Our TMDD-PBPK model, developed through a "middle-out" approach, underwent sensitivity analysis to identify key factors in the nonlinear PK of telmisartan. We found that the nonlinearity in the area under the concentration-time curve (AUC) and/or maximum concentration (Cmax) of telmisartan is sensitive to Km,OATP1B3 across all dosages. Additionally, the dissociation constant (Kd) for telmisartan binding to the AT1 receptor, along with its receptor abundance, notably influences PK at lower doses (below 20 mg). In conclusion, the nonlinear PK of telmisartan appears primarily driven by hepatic uptake saturation across all dose ranges and by AT1-receptor binding saturation, notably at lower doses.

Cluster Gauss-Newton method for a quick approximation of profile likelihood: With application to physiologically-based pharmacokinetic models.

Physiologically-based pharmacokinetic (PBPK) models can be challenging to work with because they can have too many parameters to identify from observable data. The profile likelihood method can help solve this issue by determining parameter identifiability and confidence intervals, but it involves repetitive parameter optimizations that can be time-consuming. The Cluster Gauss-Newton method (CGNM) is a parameter estimation method that efficiently searches through a wide range of parameter space. In this study, we propose a method that approximates the profile likelihood by reusing intermediate computation results from CGNM, allowing us to obtain the upper bounds of the profile likelihood without conducting additional model evaluation. This method allows us to quickly draw approximate profile likelihoods for all unknown parameters. Additionally, the same approach can be used to draw two-dimensional profile likelihoods for all parameter combinations within seconds. We demonstrate the effectiveness of this method on three PBPK models.

A case of perforation of the prostatic abscess into the rectum resulting in a rectoprostatic fistula.

Introduction: A rectoprostatic fistula complicating a prostatic abscess is extremely rare, and there are many uncertainties regarding its treatment and prognosis. Case presentation: We present the case of a 68-year-old male who presented with rectal bleeding. Computed tomography reveals a prostatic abscess in contact with the rectum. As conservative treatment with antibiotics resulted in the aggravation of symptoms, we placed a transperineal catheter placement in the prostate, followed by drainage of the abscess through the catheter. The abscess communicated with the rectum, forming a rectoprostatic fistula. Subsequently, a transverse colostomy was performed. The abscess resolved, and no recurrence was observed even after 6 months of drainage. However, the rectal ulceration that resulted in the formation of the fistula persisted, and stoma closure was not achieved. Conclusion: In cases of prostatic abscess complicated by rectoprostatic fistula, drainage of the abscess along with colostomy seems essential for a curative approach.

Physiologically-based pharmacokinetic modeling for investigating the effect of simeprevir on concomitant drugs and an endogenous biomarker of OATP1B.

The orally available anti-hepatitis C virus (HCV) drug simeprevir exhibits nonlinear pharmacokinetics at the clinical doses due to saturation of cytochrome P450 (CYP) 3A4 metabolism and organic anion transporting peptide (OATP) 1B mediated hepatic uptake. Additionally, simeprevir increases exposures of concomitant drugs by CYP3A4 and OATP1B inhibition. The objective of this study was to develop physiologically-based pharmacokinetic (PBPK) models that could describe drug-drug interactions (DDIs) of simeprevir with concomitant drugs via CYP3A4 and OATP1B inhibition, and also to capture the effects on coproporphyrin-I (CP-I), an endogenous biomarker of OATP1B. PBPK modeling estimated unbound simeprevir inhibitory constant (Ki ) of 2.89 µM against CYP3A4 in the DDI results between simeprevir and midazolam in healthy volunteers. Then, we analyzed the DDIs between simeprevir and atorvastatin, a dual substrate of CYP3A4 and OATP1B, in healthy volunteers, and unbound Ki against OATP1B was estimated to be 0.00347 µM. Finally, we analyzed the increase in the blood level of CP-I by simeprevir to verify the Ki,OATP1B . Because CP-I was measured in subjects with HCV with various hepatic fibrosis state, Monte Carlo simulation was performed to involve the decreases in expression levels of hepatic CYP3A4 and OATP1B and their interindividual variabilities. The PBPK modeling coupled with Monte Carlo simulation using the Ki,OATP1B value obtained from atorvastatin study reasonably recovered the observed relationship between CP-I and simeprevir blood levels. In conclusion, the simeprevir PBPK model developed in this study can quantitatively describe the increase in exposures of concomitant drugs and an endogenous biomarker via inhibition of CYP3A4 and OATP1B.

A 20-Year Research Overview: Quantitative Prediction of Hepatic Clearance Using the In Vitro-In Vivo Extrapolation Approach Based on Physiologically Based Pharmacokinetic Modeling and Extended Clearance Concept.

Understanding the extended clearance concept and establishing a physiologically based pharmacokinetic (PBPK) model are crucial for investigating the impact of changes in transporter and metabolizing enzyme abundance/functions on drug pharmacokinetics in blood and tissues. This mini-review provides an overview of the extended clearance concept and a PBPK model that includes transporter-mediated uptake processes in the liver. In general, complete in vitro and in vivo extrapolation (IVIVE) poses challenges due to missing factors that bridge the gap between in vitro and in vivo systems. By considering key in vitro parameters, we can capture in vivo pharmacokinetics, a strategy known as the top-down or middle-out approach. We present the latest progress, theory, and practice of the Cluster Gauss-Newton method, which is used for middle-out analyses. As examples of poor IVIVE, we discuss "albumin-mediated hepatic uptake" and "time-dependent inhibition" of OATP1Bs. The hepatic uptake of highly plasma-bound drugs is more efficient than what can be accounted for by their unbound concentration alone. This phenomenon is referred to as "albumin-mediated" hepatic uptake. IVIVE was improved by measuring hepatic uptake clearance in vitro in the presence of physiologic albumin concentrations. Lastly, we demonstrate the application of Cluster Gauss-Newton method-based analysis to the target-mediated drug disposition of bosentan. Incorporating saturable target binding and OATP1B-mediated hepatic uptake into the PBPK model enables the consideration of nonlinear kinetics across a wide dose range and the prediction of receptor occupancy over time. SIGNIFICANCE STATEMENT: There have been multiple instances where researchers' endeavors to unravel the underlying mechanism of poor in vitro-in vivo extrapolation have led to the discovery of previously undisclosed truths. These include 1) albumin-mediated hepatic uptake, 2) the target-mediated drug disposition in small molecules, and 3) the existence of a trans-inhibition mechanism by inhibitors for OATP1B-mediated hepatic uptake of drugs. Consequently, poor in vitro-in vivo extrapolation and the subsequent inquisitiveness of scientists may serve as a pivotal gateway to uncover hidden mechanisms.

Predicting In Vivo Target Occupancy (TO) Profiles via Physiologically Based Pharmacokinetic-TO Modeling of Warfarin Pharmacokinetics in Blood: Importance of Low Dose Data and Prediction of Stereoselective Target Interactions.

Warfarin is well recognized for its high-affinity and capacity-limited binding to the pharmacological target and undergoes target-mediated drug disposition. Here, we developed a physiologically based pharmacokinetic (PBPK) model that incorporated saturable target binding and other reported hepatic disposition components of warfarin. The PBPK model parameters were optimized by fitting to the reported blood pharmacokinetic (PK) profiles of warfarin with no stereoisomeric separation after oral dosing of racemic warfarin (0.1, 2, 5, or 10 mg) using the Cluster Gauss-Newton method (CGNM). The CGNM-based analysis yielded multiple "accepted" sets for six optimized parameters, which were then used to simulate the warfarin blood PK and in vivo target occupancy (TO) profiles. When further analyses examined the impact of dose selection on uncertainty in parameter estimation by the PBPK modeling, the PK data from 0.1 mg dose (well below target saturation) was important in practically identifying the target binding-related parameters in vivo. When stereoselective differences were incorporated for both hepatic disposition and target interactions, our PBPK modeling predicted that R-warfarin (of slower clearance and lower target affinity than S-warfarin) contributes to TO prolongation after oral dosing of racemic warfarin. Our results extend the validity of the approach by which the PBPK-TO modeling of blood PK profiles can yield TO prediction in vivo (applicable to the drugs with targets of high affinity and abundance and limited distribution volume via nontarget interactions). Our findings support that model-informed dose selection and PBPK-TO modeling may aid in TO and efficacy assessment in preclinical and clinical phase 1 studies. SIGNIFICANCE STATEMENT: The current physiologically based pharmacokinetic modeling incorporated the reported hepatic disposition components and target binding of warfarin and analyzed the blood pharmacokinetic (PK) profiles from varying warfarin doses, practically identifying target binding-related parameters in vivo. By implementing the stereoselective differences between R- and S-warfarin, our analysis predicted the role of R-warfarin in prolonging overall target occupancy. Our results extend the validity of analyzing blood PK profiles to predict target occupancy in vivo, which may guide efficacy assessment.

Cluster Gauss-Newton method analyses of PBPK model parameter combinations of coproporphyrin-I based on OATP1B-mediated rifampicin interaction studies.

Coproporphyrin I (CP-I) is an endogenous biomarker supporting the prediction of drug-drug interactions (DDIs) involving hepatic organic anion transporting polypeptide 1B (OATP1B). We previously constructed a physiologically-based pharmacokinetic (PBPK) model for CP-I using clinical DDI data with an OATP1B inhibitor, rifampicin (RIF). In this study, PBPK model parameters for CP-I were estimated using the cluster Gauss-Newton method (CGNM), an algorithm used to find multiple approximate solutions for nonlinear least-squares problems. Eight unknown parameters including the hepatic overall intrinsic clearance (CLint,all ), the rate of biosynthesis (vsyn ), and the OATP1B inhibition constant of RIF(Ki,u,OATP ) were estimated by fitting to the observed CP-I blood concentrations in two different clinical studies involving changing the RIF dose. Multiple parameter combinations were obtained by CGNM that could well capture the clinical data. Among those, CLint,all , Ki,u,OATP , and vsyn were sensitive parameters. The obtained Ki,u,OATP for CP-I was 5.0- and 2.8-fold lower than that obtained for statins, confirming our previous findings describing substrate-dependent Ki,u,OATP values. In conclusion, CGNM analyses of PBPK model parameter combinations enables estimation of the three essential parameters for CP-I to capture the DDI profiles, even if the other parameters remain unidentified. The CGNM also clarified the importance of appropriate combinations of other unidentified parameters to enable capture of the CP-I concentration time course under the influence of RIF. The described CGNM approach may also support the construction of robust PBPK models for additional transporter biomarkers beyond CP-I.

Frequency of null genotypes of glutathione S-transferase M1 and T1 in Japanese patients with drug-induced liver injury.

AIM: Previous reports suggest that the null genotype (*0/*0) of glutathione S-transferase (GST) M1 and/or GSTT1 could be risk factors for drug-induced liver injury (DILI). However, multi-institutional pharmacogenetic research with various suspected drugs has rarely been performed in Japan. Therefore, the aim of this study was to investigate the role of GSTM1 and GSTT1 null genotype in the occurrence of DILI in Japanese patients. METHODS: Blood samples of 270 DILI patients from 23 hospitals throughout Japan collected between 2010 and 2018 were subjected to genotyping of null genotypes of GSTM1 and GSTT1 using the SmartAmp-2 method. We also collected information on DILI types, time to onset of DILI, pharmacological classification of suspected drugs and Digestive Disease Week-Japan score, as well as genotypes of GSTM1 and GSTT1 in each patient with DILI. RESULTS: The distribution of a combination of null genotypes of GSTM1 and GSTT1 in Japanese patients with DILI was significantly different from that reported in the general Japanese population. Notably, the incidence of the GSTM1 null genotype in patients with DILI was significantly higher than that of the control population. A significant relationship between the frequency of GSTM1 and GSTT1 null genotypes and pharmacological classification of suspected drugs, clinical laboratory data for liver function, time to onset of DILI, and Digestive Disease Week-Japan scores was not observed. CONCLUSIONS: The GSTM1 null genotype was associated with an increased incidence of DILI in Japanese patients.

Clinical Relevance of Hepatic and Renal P-gp/BCRP Inhibition of Drugs: An International Transporter Consortium Perspective.

The role of P-glycoprotein (P-gp) and breast cancer resistance protein (BCRP) in drug-drug interactions (DDIs) and limiting drug absorption as well as restricting the brain penetration of drugs with certain physicochemical properties is well known. P-gp/BCRP inhibition by drugs in the gut has been reported to increase the systemic exposure to substrate drugs. A previous International Transporter Consortium (ITC) perspective discussed the feasibility of P-gp/BCRP inhibition at the blood-brain barrier and its implications. This ITC perspective elaborates and discusses specifically the hepatic and renal P-gp/BCRP (referred as systemic) inhibition of drugs and whether there is any consequence for substrate drug disposition. This perspective summarizes the clinical evidence-based recommendations regarding systemic P-gp and BCRP inhibition of drugs with a focus on biliary and active renal excretion pathways. Approaches to assess the clinical relevance of systemic P-gp and BCRP inhibition in the liver and kidneys included (i) curation of DDIs involving intravenously administered substrates or inhibitors; (ii) in vitro-to-in vivo extrapolation of P-gp-mediated DDIs at the systemic level; and (iii) curation of drugs with information available about the contribution of biliary excretion and related DDIs. Based on the totality of evidence reported to date, this perspective supports limited clinical DDI risk upon P-gp or BCRP inhibition in the liver or kidneys.

Robust physiologically based pharmacokinetic model of rifampicin for predicting drug-drug interactions via P-glycoprotein induction and inhibition in the intestine, liver, and kidney.

P-glycoprotein (P-gp) is an efflux transporter that plays an important role in the pharmacokinetics of its substrate, and P-gp activities can be altered by induction and inhibition effects of rifampicin. This study aimed to establish a physiologically based pharmacokinetic (PBPK) model of rifampicin to predict the P-gp-mediated drug-drug interactions (DDIs) and assess the DDI impact in the intestine, liver, and kidney. The induction and inhibition parameters of rifampicin for P-gp were estimated using two of seven DDI cases of rifampicin and digoxin and incorporated into our previously constructed PBPK model of rifampicin. The constructed rifampicin model was verified using the remaining five DDI cases with digoxin and five DDI cases with other P-gp substrates (talinolol and quinidine). Based on the established PBPK model, following repeated dosing of 600 mg rifampicin, the deduced net effect was an approximately threefold induction in P-gp activities in the intestine, liver, and kidney. Furthermore, in all 12 cases the predicted area under the plasma concentration-time curve ratios of the P-gp substrates were within the predefined acceptance criteria with various dosing regimens. Intestinal effects of P-gp-mediated DDIs had their greatest impact on the pharmacokinetics of digoxin and talinolol, with a minimal impact on the liver and kidney. For quinidine, predicted intestinal P-gp/cytochrome P450 3A-mediated DDIs were slightly underestimated because of the complexity of nonlinearity and transporter-enzyme interplay. These findings demonstrate that our rifampicin model can be applicable to quantitatively predict the net impact of P-gp induction and/or inhibition on diverse P-gp substrates and investigate the magnitude of DDIs in each tissue.

Transporters in Drug Development: International Transporter Consortium Update on Emerging Transporters of Clinical Importance.

During its fourth transporter workshop in 2021, the International Transporter Consortium (ITC) provided updates on emerging clinically relevant transporters for drug development. Previously highlighted and new transporters were considered based on up-to-date clinical evidence of their importance in drug-drug interactions and potential for altered drug efficacy and safety, including drug-nutrient interactions leading to nutrient deficiencies. For the first time, folate transport pathways (PCFT, RFC, and FRα) were examined in-depth as a potential mechanism of drug-induced folate deficiency and related toxicities (e.g., neural tube defects and megaloblastic anemia). However, routine toxicology studies conducted in support of drug development appear sufficient to flag such folate deficiency toxicities, whereas prospective prediction from in vitro folate metabolism and transport inhibition is not well enough established to inform drug development. Previous suggestion of a retrospective study of intestinal OATP2B1 inhibition to explain unexpected decreases in drug exposure were updated. Furthermore, when the absorption of a new molecular entity is more rapid and extensive than can be explained by passive permeability, evaluation of the OATP2B1 transport may be considered. Emerging research on hepatic and renal OAT2 is summarized, but current understanding of the importance of OAT2 was deemed insufficient to justify specific consideration for drug development. Hepatic, renal, and intestinal MRPs (MRP2, MRP3, and MRP4) were revisited. MRPs may be considered when they are suspected to be the major determinant of drug disposition (e.g., direct glucuronide conjugates); MRP2 inhibition as a mechanistic explanation for drug-induced hyperbilirubinemia remains justified. There were no major changes in recommendations from previous ITC whitepapers.

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 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 , Fa Fg , 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.

Clinical evaluation of [18F]pitavastatin for quantitative analysis of hepatobiliary transporter activity.

It is widely accepted that uptake and efflux transporters on clearance organs play crucial roles in drug disposition. Although in vitro transporter assay system can identify the intrinsic properties of the target transporters, it is not so easy to precisely predict in vivo pharmacokinetic parameters from in vitro data. Positron emission tomography (PET) imaging is a useful tool to directly assess the activity of drug transporters in humans. We recently developed a practical synthetic method for fluorine-18-labeled pitavastatin ([18F]PTV) as a PET probe for quantitative evaluation of hepatobiliary transport. In the present study, we conducted clinical PET imaging with [18F]PTV and compared the pharmacokinetic properties of the probe for healthy subjects with or without rifampicin pretreatment. Rifampicin pretreatment significantly suppressed the hepatic maximum concentration and biliary excretion of the probe to 52% and 34% of the control values, respectively. Rifampicin treatment markedly decreased hepatic uptake clearance (21% of the control), and moderately canalicular efflux clearance with regard to hepatic concentration (52% of the control). These results demonstrate that [18F]PTV is a useful probe for clinical investigation of the activities of hepatobiliary uptake/efflux transporters in humans.

Experimental and Modeling Evidence Supporting the Trans-Inhibition Mechanism for Preincubation Time-Dependent, Long-Lasting Inhibition of Organic Anion Transporting Polypeptide 1B1 by Cyclosporine A.

Cyclosporine A (CsA) and rifampin are potent inhibitors of organic anion transporting polypeptide (OATP) 1B1 and are widely used to assess the risk for drug-drug interactions. CsA displays preincubation time-dependent, long-lasting inhibition of OATP1B1 in vitro and in rats in vivo, and a proposed mechanism is the trans-inhibition by which CsA inhibits OATP1B1 from the inside of cells. The current study aimed to experimentally validate the proposed mechanism using human embryonic kidney 293 cells stably expressing OATP1B1. The uptake of CsA reached a plateau following an approximate 60-minute incubation, with the cell-to-buffer concentration ratio of 3930, reflective of the high-affinity, high-capacity intracellular binding of CsA. The time course of CsA uptake was analyzed to estimate the kinetic parameters for permeability clearance and intracellular binding. When the OATP1B1-mediated uptake of [3H]estradiol-17ß-glucuronide was measured following preincubation with CsA for 5 to 120 minutes, apparent Ki values became lower with longer preincubation. Our kinetic modeling incorporated the two reversible inhibition constants [Ki,trans and Ki,cis for the inhibition from inside (trans-inhibition) and outside (cis-inhibition) of cells, respectively] and estimated Ki,trans value of CsA was smaller by 48-fold than the estimated Ki,cis value. Rifampin also displayed preincubation time-dependent inhibition of OATP1B1, albeit the extent of enhancement was only twofold. The current study provides experimental evidence for the preincubation time-dependent shift of apparent Ki values and a mechanistic basis for physiologically based pharmacokinetic modeling that incorporates permeability clearance, extensive intracellular binding, and asymmetry of Ki values between the inside and outside of cells. SIGNIFICANCE STATEMENT: In vitro data and kinetic modeling support that preincubation time-dependent, long-lasting inhibition of OATP1B1 by CsA can be explained by the extensive intracellular binding and reversible OATP1B1 inhibition intracellularly (trans-inhibition) as well as extracellularly (cis-inhibition). For inhibitors to display time-dependency, the following factors were found important: time to reach a steady-state cellular concentration, trans-inhibition potency relative to cis-inhibition, and the degree of cellular inhibitor accumulation. This study would aid in the accurate prediction of drug-drug interactions mediated by OATP1B1 inhibition.

Effect of Cyclosporin A and Impact of Dose Staggering on OATP1B1/1B3 Endogenous Substrates and Drug Probes for Assessing Clinical Drug Interactions.

This study was designed to assess the quantitative performance of endogenous biomarkers for organic anion transporting polypeptide (OATP) 1B1/1B3-mediated drug-drug interactions (DDIs). Ten healthy volunteers orally received OATP1B1/1B3 probe cocktail (0.2 mg pitavastatin, 1 mg rosuvastatin, and 2 mg valsartan) and an oral dose of cyclosporin A (CysA, 20 mg and 75 mg) separated by a 1-hour interval (20 mg (-1 hour), and 75 mg (-1 hour)). CysA 75 mg was also given with a 3-hour interval (75 mg (-3 hours)) to examine the persistence of OATP1B1/1B3 inhibition. The area under the plasma concentration-time curve ratios (AUCRs) were 1.63, 3.46, and 2.38 (pitavastatin), 1.39, 2.16, and 1.81 (rosuvastatin), and 1.42, 1.77, and 1.85 (valsartan), at 20 mg, 75 mg (-1 hour) and 75 mg (-3 hours) of CysA, respectively. CysA effect on OATP1B1/1B3 was unlikely to persist at the dose examined. Among 26 putative OATP1B1/1B3 biomarkers evaluated, AUCR and maximum concentration ratio (Cmax R) of CP-I showed the highest Pearson's correlation coefficient with CysA AUC (0.94 and 0.93, respectively). Correlation between AUCR of pitavastatin, and Cmax R or AUCR of CP-I were consistent between this study and our previous study using rifampicin as an OATP1B1/1B3 inhibitor. Nonlinear regression analysis of AUCR-1 of pitavastatin and CP-I against CysA Cmax yielded Ki,OATP1B1/1B3,app (109 ± 35 and 176 ± 42 nM, respectively), similar to the Ki ,OATP1B1/1B3 estimated by our physiologically-based pharmacokinetic model analysis described previously (107 nM). The endogenous OATP1B1/1B3 biomarkers, particularly Cmax R and AUCR of CP-I, corroborates OATP1B1/1B3 inhibition and yields valuable information that improve accurate DDI predictions in drug development, and enhance our understanding of interindividual variability in the magnitude of DDIs.

Consideration of albumin-mediated hepatic uptake for highly protein-bound anionic drugs: Bridging the gap of hepatic uptake clearance between in vitro and in vivo.

The accuracy in predicting in vivo hepatic clearance of drugs from in vitro data (often termed as in vitro-to-in vivo extrapolation, IVIVE) has improved in part by applying the extended-clearance concept that considers the interplay between hepatic metabolism and uptake/efflux processes. However, the IVIVE-based prediction performs poorly in predicting the hepatic uptake clearance of highly albumin-bound anionic drugs. Their hepatic uptake clearances tend to be much higher than expected based on the free-drug theory. Such observation can be attributable to a phenomenon called albumin-mediated hepatic uptake, for which various models have been thus far proposed. Our group has been applying a facilitated-dissociation model, which assumes the enhanced dissociation of the drug-albumin complex upon interaction with the cell surface. By considering the albumin-mediated hepatic uptake (using the facilitated-dissociation model or alternative kinetic models), a number of investigations demonstrated the improvement in the prediction accuracy for the hepatic clearance of highly protein-bound anionic drugs that are substrates for hepatic uptake transporters. This review summarizes the reported kinetic analyses of the albumin-mediated hepatic uptake of highly albumin-bound drugs concerning the IVIVE and the clinical and physiological relevance.

Physiologically Based Pharmacokinetic/Pharmacodynamic Model for the Treatment of Dengue Infections Applied to the Broad Spectrum Antiviral Soraphen A.

While a drug treatment is unavailable, the global incidence of Dengue virus (DENV) infections and its associated severe manifestations continues to rise. We report the construction of the first physiologically based pharmacokinetic/pharmacodynamic (PBPK/PD) model that predicts viremia levels in relevant target organs based on preclinical data with the broad spectrum antiviral soraphen A (SorA), an inhibitor of the host cell target acetyl-CoA-carboxylase. SorA was highly effective against DENV in vitro (EC50 = 4.7 nM) and showed in vivo efficacy by inducing a significant reduction of viral load in the spleen and liver of IFNAR-/- mice infected with DENV-2. PBPK/PD predictions for SorA matched well with the experimental infection data. Transfer to a human PBPK/PD model for DENV to mimic a clinical scenario predicted a reduction in viremia by more than one log10 unit for an intravenous infusion regimen of SorA. The PBPK/PD model is applicable to any DENV drug lead and, thus, represents a valuable tool to accelerate and facilitate DENV drug discovery and development.

Revisiting Nonlinear Bosentan Pharmacokinetics by Physiologically Based Pharmacokinetic Modeling: Target Binding, Albeit Not a Major Contributor to Nonlinearity, Can Offer Prediction of Target Occupancy.

Bosentan is a high-affinity antagonist of endothelin receptors and one of the earliest examples for target-mediated drug disposition [a type of nonlinear pharmacokinetics (PKs) caused by saturable target binding]. The previous physiologically based PK (PBPK) modeling indicated that the nonlinear PKs of bosentan was explainable by considering saturable hepatic uptake. However, it remained unexamined to what extent the saturable target binding contributes to the nonlinear PKs of bosentan. Here, we developed a PBPK model incorporating saturable target binding and hepatic uptake and analyzed the clinical bosentan PK data using the cluster Gauss-Newton method (CGNM). The PBPK model without target binding fell short in capturing the bosentan concentrations below 100 nM, based on the PK profiles and the goodness-of-fit plot. Both global and local identifiability analyses (using the CGNM and Fisher information matrix, respectively) informed that the target binding parameters were identifiable only if the observations from the lowest dose (10 mg) were included. By analyzing blood PK profiles alone, the PBPK model with target binding yielded practically identifiable target binding parameters and predicted the maximum target occupancies of 0.6-0.8 at clinical bosentan doses. Our results indicate that target binding, albeit not a major contributor to the nonlinear bosentan PKs, may offer a prediction of target occupancy from blood PK profiles alone and potential guidance on achieving optimal efficacy outcomes, under the condition when the high-affinity drug target is responsible for the efficacy of interest and when the dose ranges cover varying degrees of target binding. SIGNIFICANCE STATEMENT: By incorporating saturable target binding, our physiologically based pharmacokinetic (PBPK) model predicted in vivo target occupancy of bosentan based only on the blood concentration-time profiles obtained from a wide range of doses. Our analysis highlights the potential utility of PBPK models that incorporate target binding in predicting target occupancy in vivo.

Intestinal P-gp and Putative Hepatic OATP1B Induction: International Transporter Consortium Perspective on Drug Development Implications.

There is an increasing interest in transporter induction (i.e., decreased systemic drug exposure due to increased efflux-limited absorption or transporter-mediated clearance) as a mechanism of drug-drug interactions (DDIs), although evidence of clinical relevance is still evolving. Intestinal P-glycoprotein (P-gp) and hepatic organic anion transporting polypeptides 1B (OATP1B) can be important determinants of drug absorption and disposition, as well as targets for DDIs. Current data indicate that intestinal P-gp protein levels can be induced up to threefold to fourfold in humans primarily with pregnane X receptor (PXR) activators, and that this induction can decrease the systemic exposure of drugs with P-gp efflux-limited absorption (e.g., ≤ 67% decrease in the exposure of total dabigatran following rifampin multiple oral dosing). Evaluation of the clinical relevance of P-gp induction as a DDI mechanism must consider the induction potential of the perpetrator drug for P-gp and attenuation of exposure of the victim drug in the context of its therapeutic window. Practical drug development recommendations are provided herein. Reports are contradictory on OATP1B induction by PXR activators in human hepatocytes and liver biopsies. Some clinical investigations demonstrated that rifampin pretreatment decreased exposure of OATP1B substrates, while other studies found no differences, and the potential involvement of other mechanisms in these observed DDIs cannot be definitively ruled out. Thus, further studies are needed to understand hepatic OATP1B induction and potential involvement of other mechanisms contributing to reduced exposure of OATP1B substrates. This review critically summarizes the state-of-the-art on intestinal P-gp and hepatic OATP1B induction, and highlights implications for drug development.