Search

Hepatic impairment (HI) is known to modulate drug disposition and may lead to elevated plasma exposure. The aim of this study was to quantitate the in vivo OATP1B-mediated hepatic uptake activity in populations with varying degrees of HI. First, we measured baseline levels of plasma coproporphyrin-I, an endogenous OATP1B biomarker, in an open-label, parallel cohort study in adult subjects with normal liver function and mild, moderate and severe HI (n=24, 6/cohort). The geometric mean plasma concentrations of coproporphyrin-I were 1.66-fold, 2.81-fold (P0.05) and 7.78-fold (P0.0001) higher in mild, moderate and severe impairment than those in healthy control. Second, we developed a dataset of 21 OATP1B substrate drugs with HI data extracted from literature. Median disease-to-healthy plasma area under the curve (AUC) ratios for substrate drugs were approximately 1.4, 3.0 and 6.4 for mild, moderate and severe HI, respectively. Additionally, significant linear relationship was noted between AUC ratios of substrate drugs without and with coadministration of rifampin, a prototypic OATP1B inhibitor, and AUC ratios in moderate (P0.01) and severe (P0.0001) hepatic impairment. Third, a physiologically-based pharmacokinetic model analysis was conducted with ten substrate drugs following estimation of relative OATP1B functional activity in virtual disease population models using coproporphyrin-I data and verification of substrate models with rifampin drug-drug interaction data. This approach adequately predicted plasma AUC change particularly in moderate (9 of 10 within 2-fold) and severe (5 of 5 within 2-fold) hepatic impairment. Collective findings indicate progressive reduction, by as much as 90-92%, in OATP1B activity in the HI population.

PDF file, 141K, Photographs of canine muzzles showing skin flushing following PF-04880594 treatment.

Supplemental Methods

ElasticNet Tables for PF-06939999 in NSCLC

Supplemental Figures S1-S8

PDF file - 44K

PDF file, 75K, Body weight changes in nude mice administered RAF inhibitor (PF-04880594) are maintained in mice co-treated with MEK inhibitor (PD-0325901).

Clinical trials of selective RAF inhibitors in patients with melanoma tumors harboring activated BRAFV600E have produced very promising results, and a RAF inhibitor has been approved for treatment of advanced melanoma. However, about a third of patients developed resectable skin tumors during the course of trials. This is likely related to observations that RAF inhibitors activate extracellular signal–regulated kinase (ERK) signaling, stimulate proliferation, and induce epithelial hyperplasia in preclinical models. Because these findings raise safety concerns about RAF inhibitor development, we further investigated the underlying mechanisms. We showed that the RAF inhibitor PF-04880594 induces ERK phosphorylation and RAF dimerization in those epithelial tissues that undergo hyperplasia. Hyperplasia and ERK hyperphosphorylation are prevented by treatment with the mitogen-activated protein/extracellular signal–regulated kinase (MEK) inhibitor PD-0325901 at exposures that extrapolate to clinically well-tolerated doses. To facilitate mechanistic and toxicologic studies, we developed a three-dimensional cell culture model of epithelial layering that recapitulated the RAF inhibitor–induced hyperplasia and reversal by MEK inhibitor in vitro. We also showed that PF-04880594 stimulates production of the inflammatory cytokine interleukin 8 in HL-60 cells, suggesting a possible mechanism for the skin flushing observed in dogs. The complete inhibition of hyperplasia by MEK inhibitor in epithelial tissues does not seem to reduce RAF inhibitor efficacy and, in fact, allows doubling of the PF-04880594 dose without toxicity usually associated with such doses. These findings indicated that combination treatment with MEK inhibitors might greatly increase the safety and therapeutic index of RAF inhibitors for the treatment of melanoma and other cancers. Mol Cancer Ther; 11(10); 2274–83. ©2012 AACR.

PDF file, 39K, Relationship between oral doses and average plasma concentrations of PD-0325901 in athymic nu/nu mice.

PDF file, 66K.

XLS file - 76K

PDF file - 348K

PDF file - 11K

Purpose: Axitinib (AG-013736) is a potent and selective inhibitor of vascular endothelial growth factor (VEGF) receptor tyrosine kinases 1 to 3 that is in clinical development for the treatment of solid tumors. We provide a comprehensive description of its in vitro characteristics and activities, in vivo antiangiogenesis, and antitumor efficacy and translational pharmacology data.Experimental Design: The potency, kinase selectivity, pharmacologic activity, and antitumor efficacy of axitinib were assessed in various nonclinical models.Results: Axitinib inhibits cellular autophosphorylation of VEGF receptors (VEGFR) with picomolar IC50 values. Counterscreening across multiple kinase and protein panels shows it is selective for VEGFRs. Axitinib blocks VEGF-mediated endothelial cell survival, tube formation, and downstream signaling through endothelial nitric oxide synthase, Akt and extracellular signal-regulated kinase. Following twice daily oral administration, axitinib produces consistent and dose-dependent antitumor efficacy that is associated with blocking VEGFR-2 phosphorylation, vascular permeability, angiogenesis, and concomitant induction of tumor cell apoptosis. Axitinib in combination with chemotherapeutic or targeted agents enhances antitumor efficacy in many tumor models compared with single agent alone. Dose scheduling studies in a human pancreatic tumor xenograft model show that simultaneous administration of axitinib and gemcitabine without prolonged dose interruption or truncation of axitinib produces the greatest antitumor efficacy. The efficacious drug concentrations predicted in nonclinical studies are consistent with the range achieved in the clinic. Although axitinib inhibits platelet-derived growth factor receptors and KIT with nanomolar in vitro potencies, based on pharmacokinetic/pharmacodynamic analysis, axitinib acts primarily as a VEGFR tyrosine kinase inhibitor at the current clinical exposure.Conclusions: The selectivity, potency for VEGFRs, and robust nonclinical activity may afford broad opportunities for axitinib to improve cancer therapy.

PDF file - 481K

PDF file - 49K

Trastuzumab CRD Calculation Details

Description of common pharmacokinetic parameters used in the CRD calculations.

Dasatinib CRD Calculation Details

Vismodegib CRD Calculation Details

PDF file - 16K

Supplementary Figures S1-S7 from Nonclinical Antiangiogenesis and Antitumor Activities of Axitinib (AG-013736), an Oral, Potent, and Selective Inhibitor of Vascular Endothelial Growth Factor Receptor Tyrosine Kinases 1, 2, 3

Purpose: The translation of nonclinical oncology studies is a subject of continuous debate. We propose that translational oncology studies need to optimize both pharmacokinetic (drug exposure) and pharmacodynamic (xenograft model) aspects. While improvements in pharmacodynamic translatability can be obtained by choosing cell lines or patient-derived xenograft models closer to the clinical indication, significant ambiguity and variability exists when optimizing the pharmacokinetic translation of small molecule and biotherapeutic agents.Experimental Design and Results: In this work, we propose a pharmacokinetic-based strategy to select nonclinical doses for approved drug molecules. We define a clinically relevant dose (CRD) as the dosing regimen in mice that most closely approximates the relevant pharmacokinetic metric in humans. Such metrics include area under the time–concentration curve and maximal or minimal concentrations within the dosing interval. The methodology is applied to six drugs, including targeted agents and chemotherapeutics, small and large molecules (erlotinib, dasatinib, vismodegib, trastuzumab, irinotecan, and capecitabine). The resulting efficacy response at the CRD is compared with clinical responses.Conclusions: We conclude that nonclinical studies designed with the appropriate CRDs of approved drug molecules will maximize the translatability of efficacy results, which is critical when testing approved and investigational agents in combination. Clin Cancer Res; 23(4); 1080–90. ©2016 AACR.

Full Description with references for the data included in Figure 4.

Purpose: Triple-negative breast cancer (TNBC) patients have poor prognoses and survival outcomes such that the development of new targeted therapies is in strong demand. Mechanisms associated with high proliferation and aggressive tumor progression, such as PI3K/PTEN aberration, epidermal growth factor receptor (EGFR) overexpression, and cell-cycle upregulation, play important roles in TNBC. The molecular chaperone Hsp90 is required for the conformational maturation and stability of a variety of proteins in multiple pathways, such as EGFR, AKT, Raf, cdk4, etc. Therefore, an Hsp90 inhibitor may show therapeutic benefit in TNBC by targeting multiple pathways.Experimental Design: The novel oral Hsp90 inhibitor PF-4942847 was characterized in multiple in vitro and in vivo assays to determine its antitumor activity in TNBC cell lines. In addition, the correlation of AKT degradation and Hsp70 induction in host peripheral blood lymphocytes (PBL) and xenograft tumors was determined.Results: PF-4942847 induces degradation of multiple client proteins, cell-cycle block, apoptosis, and inhibits cell proliferation in TNBC lines, subsequently leading to tumor growth inhibition in mouse xenograft models. The correlation of AKT degradation and Hsp70 induction between PBLs and xenograft tumors reveals a differential modulation of Hsp90 activity between host and tumor tissues, and suggests that AKT degradation in PBLs may serve as a pharmacodynamic biomarker in future clinical development.Conclusions: The novel oral Hsp90 inhibitor, PF-4942847, is a candidate for clinical development in TNBC by collaboratively targeting multiple signaling pathways. In addition, AKT degradation in PBLs may serve as a biomarker in clinical development. Clin Cancer Res; 17(16); 5432–42. ©2011 AACR.

PDF file - 356K

Supplementary Tables 1-2, Figures 1-3 from An Orally Available Small-Molecule Inhibitor of c-Met, PF-2341066, Exhibits Cytoreductive Antitumor Efficacy through Antiproliferative and Antiangiogenic Mechanisms

The c-Met receptor tyrosine kinase and its ligand, hepatocyte growth factor (HGF), have been implicated in the progression of several human cancers and are attractive therapeutic targets. PF-2341066 was identified as a potent, orally bioavailable, ATP-competitive small-molecule inhibitor of the catalytic activity of c-Met kinase. PF-2341066 was selective for c-Met (and anaplastic lymphoma kinase) compared with a panel of >120 diverse tyrosine and serine-threonine kinases. PF-2341066 potently inhibited c-Met phosphorylation and c-Met–dependent proliferation, migration, or invasion of human tumor cells in vitro (IC50 values, 5–20 nmol/L). In addition, PF-2341066 potently inhibited HGF-stimulated endothelial cell survival or invasion and serum-stimulated tubulogenesis in vitro, suggesting that this agent also exhibits antiangiogenic properties. PF-2341066 showed efficacy at well-tolerated doses, including marked cytoreductive antitumor activity, in several tumor models that expressed activated c-Met. The antitumor efficacy of PF-2341066 was dose dependent and showed a strong correlation to inhibition of c-Met phosphorylation in vivo. Near-maximal inhibition of c-Met activity for the full dosing interval was necessary to maximize the efficacy of PF-2341066. Additional mechanism-of-action studies showed dose-dependent inhibition of c-Met–dependent signal transduction, tumor cell proliferation (Ki67), induction of apoptosis (caspase-3), and reduction of microvessel density (CD31). These results indicated that the antitumor activity of PF-2341066 may be mediated by direct effects on tumor cell growth or survival as well as antiangiogenic mechanisms. Collectively, these results show the therapeutic potential of targeting c-Met with selective small-molecule inhibitors for the treatment of human cancers. [Cancer Res 2007;67(9):4408–17]

Protein arginine methyltransferase 5 (PRMT5) overexpression in hematologic and solid tumors methylates arginine residues on cellular proteins involved in important cancer functions including cell-cycle regulation, mRNA splicing, cell differentiation, cell signaling, and apoptosis. PRMT5 methyltransferase function has been linked with high rates of tumor cell proliferation and decreased overall survival, and PRMT5 inhibitors are currently being explored as an approach for targeting cancer-specific dependencies due to PRMT5 catalytic function. Here, we describe the discovery of potent and selective S-adenosylmethionine (SAM) competitive PRMT5 inhibitors, with in vitro and in vivo characterization of clinical candidate PF-06939999. Acquired resistance mechanisms were explored through the development of drug resistant cell lines. Our data highlight compound-specific resistance mutations in the PRMT5 enzyme that demonstrate structural constraints in the cofactor binding site that prevent emergence of complete resistance to SAM site inhibitors. PRMT5 inhibition by PF-06939999 treatment reduced proliferation of non–small cell lung cancer (NSCLC) cells, with dose-dependent decreases in symmetric dimethyl arginine (SDMA) levels and changes in alternative splicing of numerous pre-mRNAs. Drug sensitivity to PF-06939999 in NSCLC cells associates with cancer pathways including MYC, cell cycle and spliceosome, and with mutations in splicing factors such as RBM10. Translation of efficacy in mouse tumor xenograft models with splicing mutations provides rationale for therapeutic use of PF-06939999 in the treatment of splicing dysregulated NSCLC.

As one of the key components in model‐informed drug discovery and development, physiologically‐based pharmacokinetic (PBPK) modeling linked with in vitro‐to‐in vivo extrapolation (IVIVE) is widely applied to quantitatively predict drug–drug interactions (DDIs) on drug‐metabolizing enzymes and transporters. This study aimed to investigate an IVIVE for intestinal P‐glycoprotein (Pgp, ABCB1)‐mediated DDIs among three Pgp substrates, digoxin, dabigatran etexilate, and quinidine, and two Pgp inhibitors, itraconazole and verapamil, via PBPK modeling. For Pgp substrates, assuming unbound Michaelis‐Menten constant (K m) to be intrinsic, in vitro‐to‐in vivo scaling factors for maximal Pgp‐mediated efflux rate (J max) were optimized based on the clinically observed results without co‐administration of Pgp inhibitors. For Pgp inhibitors, PBPK models utilized the reported in vitro values of Pgp inhibition constants (K i), 1.0 μM for itraconazole and 2.0 μM for verapamil. Overall, the PBPK modeling sufficiently described Pgp‐mediated DDIs between these substrates and inhibitors with the prediction errors of less than or equal to ±25% in most cases, suggesting a reasonable IVIVE for Pgp kinetics in the clinical DDI results. The modeling results also suggest that Pgp kinetic parameters of both the substrates (K m and J max) and the inhibitors (K i) are sensitive to Pgp‐mediated DDIs, thus being key for successful DDI prediction. It would also be critical to incorporate appropriate unbound inhibitor concentrations at the site of action into PBPK models. The present results support a quantitative prediction of Pgp‐mediated DDIs using in vitro parameters, which will significantly increase the value of in vitro studies to design and run clinical DDI studies safely and effectively.

PF06821497 has been identified as an orally available small-molecule enhancer of zeste homolog 2 inhibitor. The objectives of the present study were to characterize pharmacokinetic-pharmacodynamic-disease relationships of PF06821497 in xenograft mouse models with diffuse large B-cell lymphoma (Karpas422). An indirect-response model reasonably fit dose-dependent pharmacodynamic responses [histone H3 on lysine 27 (H3K27) me3 inhibition] with an unbound EC50 of 76 nM, whereas a signal-transduction model sufficiently fit dose-dependent disease responses (tumor growth inhibition) with an unbound tumor stasis concentration (Tsc) of 168 nM. Thus, effective concentration for 70% of maximal effect (EC70) for H3K27me3 inhibition was roughly comparable to Tsc, suggesting that 70% H3K27me3 inhibition could be required for tumor stasis. Consistently, an integrated pharmacokinetic-pharmacodynamic-disease model adequately describing tumor growth inhibition also suggested that ∼70% H3K27me3 inhibition was associated with tumor stasis. Based on these results, we would propose that an EC70 estimate for H3K27me3 inhibition corresponding to tumor stasis could be considered a minimum target efficacious concentration of PF06821497 in cancer patients. SIGNIFICANCE STATEMENT Using a mathematical modeling approach, the quantitative relationships of an orally available anticancer small-molecule enhancer of zeste homolog 2 inhibitor, PF06821497, were characterized among pharmacokinetics, pharmacodynamic biomarker inhibition, and disease responses in nonclinical xenograft models with diffuse large B-cell lymphoma. The modeling results suggest that >70% histone H3 on lysine 27 (H3K27) me3 inhibition would be required for tumor stasis (i.e., 100% tumor growth inhibition). Accordingly, we would propose that an effective concentration for 70% of maximal effect estimate for H3K27me3 inhibition could be considered a minimum target efficacious concentration of PF06821497 in cancer patients.

Physiologically-based pharmacokinetic (PBPK) modeling is a powerful tool to quantitatively describe drug disposition profiles in vivo, thereby providing an alternative to predict drug-drug interactions (DDIs) that have not been tested clinically. This study aimed to predict effects of rifampin-mediated intestinal P-glycoprotein (Pgp) induction on pharmacokinetics of Pgp substrates via PBPK modeling. First, we selected four Pgp substrates (digoxin, talinolol, quinidine, and dabigatran etexilate) to derive in vitro to in vivo scaling factors for intestinal Pgp kinetics. Assuming unbound Michaelis-Menten constant (Km ) to be intrinsic, we focused on the scaling factors for maximal efflux rate (Jmax ) to adequately recover clinically observed results. Next, we predicted rifampin-mediated fold increases in intestinal Pgp abundances to reasonably recover clinically observed DDI results. The modeling results suggested that threefold to fourfold increases in intestinal Pgp abundances could sufficiently reproduce the DDI results of these Pgp substrates with rifampin. Hence, the obtained fold increases can potentially be applicable to DDI prediction with other Pgp substrates.

Quantitative assessment of drug‐drug interactions (DDIs) involving breast cancer resistance protein (BCRP) inhibition is challenged by overlapping substrate/inhibitor specificity. This study used physiologically‐based pharmacokinetic (PBPK) modeling to delineate the effects of inhibitor drugs on BCRP‐ and organic anion transporting polypeptide (OATP)1B‐mediated disposition of rosuvastatin, which is a recommended BCRP clinical probe. Initial static model analysis using in vitro inhibition data suggested BCRP/OATP1B DDI risk while considering regulatory cutoff criteria for a majority of inhibitors assessed (25 of 27), which increased rosuvastatin plasma exposure to varying degree (~ 0–600%). However, rosuvastatin area under plasma concentration‐time curve (AUC) was minimally impacted by BCRP inhibitors with calculated G‐value (= gut concentration/inhibition potency) below 100. A comprehensive PBPK model accounting for intestinal (OATP2B1 and BCRP), hepatic (OATP1B, BCRP, and MRP4), and renal (OAT3) transport mechanisms was developed for rosuvastatin. Adopting in vitro inhibition data, rosuvastatin plasma AUC changes were predicted within 25% error for 9 of 12 inhibitors evaluated via PBPK modeling. This study illustrates the adequacy and utility of a mechanistic model‐informed approach in quantitatively assessing BCRP‐mediated DDIs.

A pyridone-derived phosphate prodrug of an enhancer of zeste homolog 2 (EZH2) inhibitor was designed and synthesized to improve the inhibitor's aqueous solubility. This prodrug (compound 5) was profiled in pharmacokinetic experiments to assess its ability to deliver the corresponding parent compound (compound 2) to animals in vivo following oral administration. Results from these studies showed that the prodrug was efficiently converted to its parent compound in vivo. In separate experiments, the prodrug demonstrated impressive in vivo tumor growth inhibition in a diffuse large B-cell lymphoma Karpas-422 cell line-derived xenograft model. The described prodrug strategy is expected to be generally applicable to poorly soluble pyridone-containing EZH2 inhibitors and provides a new option to enable such compounds to achieve sufficiently high exposures in vivo.

Volume of distribution at steady state (Vss) is an important pharmacokinetic parameter of a drug candidate. In this study, Vss prediction accuracy was evaluated by using: (1) seven methods for rat with 56 compounds, (2) four methods for human with 1276 compounds, and (3) four in vivo methods and three Kp (partition coefficient) scalar methods from scaling of three preclinical species with 125 compounds. The results showed that the global QSAR models outperformed the PBPK methods. Tissue fraction unbound (fu,t) method with adipose and muscle also provided high Vss prediction accuracy. Overall, the high performing methods for human Vss prediction are the global QSAR models, Oie-Tozer and equivalency methods from scaling of preclinical species, as well as PBPK methods with Kp scalar from preclinical species. Certain input parameter ranges rendered PBPK models inaccurate due to mass balance issues. These were addressed using appropriate theoretical limit checks. Prediction accuracy of tissue Kp were also examined. The fu,t method predicted Kp values more accurately than the PBPK methods for adipose, heart and muscle. All the methods overpredicted brain Kp and underpredicted liver Kp due to transporter effects. Successful Vss prediction involves strategic integration of in silico, in vitro and in vivo approaches.

Bosutinib is an orally available Src/Abl tyrosine kinase inhibitor indicated for the treatment of patients with Ph+ chronic myelogenous leukemia at a clinically recommended dose of 500 mg once daily. Clinical results indicated that increases in bosutinib oral exposures were supraproportional at the lower doses (50-200 mg) and approximately dose-proportional at the higher doses (200-600 mg). Bosutinib is a substrate of CYP3A4 and P-glycoprotein and exhibits pH-dependent solubility with moderate intestinal permeability. These findings led us to investigate the factors influencing the underlying pharmacokinetic mechanisms of bosutinib with physiologically based pharmacokinetic (PBPK) models. Our primary objectives were to: 1) refine the previously developed bosutinib PBPK model on the basis of the latest oral bioavailability data and 2) verify the refined PBPK model with P-glycoprotein kinetics on the basis of the bosutinib drug-drug interaction (DDI) results with ketoconazole and rifampin. Additionally, the verified PBPK model was applied to predict bosutinib DDIs with dual CYP3A/P-glycoprotein inhibitors. The results indicated that 1) the refined PBPK model adequately described the observed plasma concentration-time profiles of bosutinib and 2) the verified PBPK model reasonably predicted the effects of ketoconazole and rifampin on bosutinib exposures by accounting for intestinal P-glycoprotein inhibition/induction. These results suggested that bosutinib DDI mechanism could involve not only CYP3A4-mediated metabolism but also P-glycoprotein-mediated efflux on absorption. In summary, P-glycoprotein kinetics could constitute an element in the PBPK models critical to understanding the pharmacokinetic mechanism of dual CYP3A/P-glycoprotein substrates, such as bosutinib, that exhibit nonlinear pharmacokinetics owing largely to a saturation of intestinal P-glycoprotein-mediated efflux.

A general objective of drug-drug interaction (DDI) studies is to determine whether potential interactions of new molecular entities with concomitantly administered other drugs exist and, if DDIs occur, whether dosage adjustments are required. A typical end point for DDI evaluations is the ratio of area under the plasma concentration-time curve (AUC) of substrate drugs (AUCR), whereas the ratios of maximal plasma concentration (Cmax ) and terminal half-life (t1/2 ) are also important to understand DDI mechanisms (Cmax R and t1/2 R, respectively). Because changes in substrate AUC by precipitant drugs ultimately result from alterations of Cmax and t1/2 , AUCR can be considered a hybrid parameter of Cmax R and t1/2 R, for example, AUCR ≈ Cmax R × t1/2 R. The primary objective of this study was to investigate the relationships between AUCR, Cmax R, and t1/2 R in physiologically based pharmacokinetic model-predicted and clinically observed DDI results. First, the model-predicted results showed the excellent proportional relationship between AUCR and (Cmax R × t1/2 R) in DDI results of virtual substrates having a wide range of oral bioavailability with coadministration of ketoconazole, ritonavir, and rifampin. Second, the reasonable proportional relationships were also observed in the clinically observed DDI results of midazolam and statins (atorvastatin, cerivastatin, fluvastatin, lovastatin, pitavastatin, pravastatin, rosuvastatin, and simvastatin) with various inhibitors and inducers. Finally, these results suggest that utilization of the proportional relationship between AUCR and (Cmax R × t1/2 R) can provide an additional framework to further interpret DDI results reasonably and clearly. Furthermore, the proportional relationship can be purposely used to assess study design and pharmacokinetic analyses in DDI studies.

A new series of lactam-derived EZH2 inhibitors was designed via ligand-based and physicochemical-property-based strategies to address metabolic stability and thermodynamic solubility issues associated with previous lead compound 1. The new inhibitors incorporated an sp3 hybridized carbon atom at the 7-position of the lactam moiety present in lead compound 1 as a replacement for a dimethylisoxazole group. This transformation enabled optimization of the physicochemical properties and potency compared to compound 1. Analysis of relationships between calculated log D (clogD) values and in vitro metabolic stability and permeability parameters identified a clogD range that afforded an increased probability of achieving favorable ADME data in a single molecule. Compound 23a exhibited the best overlap of potency and pharmaceutical properties as well as robust tumor growth inhibition in vivo and was therefore advanced as a development candidate (PF-06821497). A crystal structure of 23a in complex with the three-pr...

Bosutinib is an orally available Src/Abl tyrosine kinase inhibitor indicated for the treatment of patients with Philadelphia chromosome-positive chronic myelogenous leukemia. Bosutinib is predominantly metabolized by CYP3A4 as the primary clearance mechanism. The main objectives of this study were to 1) develop physiologically based pharmacokinetic (PBPK) models of bosutinib; 2) verify and refine the PBPK models based on clinical study results of bosutinib single-dose drug-drug interaction (DDI) with ketoconazole and rifampin, as well as single-dose drug-disease interaction (DDZI) in patients with renal and hepatic impairment; 3) apply the PBPK models to predict DDI outcomes in patients with weak and moderate CYP3A inhibitors; and 4) apply the PBPK models to predict DDZI outcomes in renally and hepatically impaired patients after multiple-dose administration. Results showed that the PBPK models adequately predicted bosutinib oral exposures in patients after single- and multiple-dose administrations. The PBPK models also reasonably predicted changes in bosutinib exposures in the single-dose DDI and DDZI results, suggesting that the PBPK models were sufficiently developed and verified based on the currently available data. Finally, the PBPK models predicted 2- to 4-fold increases in bosutinib exposures by moderate CYP3A inhibitors, as well as comparable increases in bosutinib exposures in renally and hepatically impaired patients between single- and multiple-dose administrations. Given the challenges in conducting numerous DDI and DDZI studies of anticancer drugs in patients, we believe that the PBPK models verified in our study would be valuable to reasonably predict bosutinib exposures under various scenarios that have not been tested clinically.