Translational pharmacokinetic-pharmacodynamic modeling from nonclinical to clinical development: a case study of anticancer drug, crizotinib

Exploring inter-ethnic and inter-patient variability and optimal dosing of osimertinib: a physiologically based pharmacokinetic modeling approach

Purpose: This study aimed to develop and validate a physiologically based pharmacokinetic (PBPK) model for osimertinib (OSI) to predict plasma trough concentration (Ctrough) and pulmonary EGFRm+ (T790M and L858R mutants) inhibition in Caucasian, Japanese, and Chinese populations. The PBPK model was also utilized to investigate inter-ethnic and inter-patient differences in OSI pharmacokinetics (PK) and determine optimal dosing regimens. Methods: Population PBPK models of OSI for healthy and disease populations were developed using physicochemical and biochemical properties of OSI and physiological parameters of different groups. And then the PBPK models were validated using the multiple clinical PK and drug-drug interaction (DDI) study data. Results: The model demonstrated good consistency with the observed data, with most of prediction-to-observation ratios of 0.8-1.25 for AUC, Cmax, and Ctrough. The PBPK model revealed that plasma exposure of OSI was approximately 2-fold higher in patients compared to healthy individuals, and higher exposure observed in Caucasians compared to other ethnic groups. This was primarily attributed to a lower CL/F of OSI in patients and Caucasian. The PBPK model displayed that key factors influencing PK and EGFRm+ inhibition differences included genetic polymorphism of CYP3A4, CYP1A2 expression, plasma free concentration (fup), albumin level, and auto-inhibition/induction on CYP3A4. Inter-patient PK variability was most influenced by CYP3A4 variants, fup, and albumin level. The PBPK simulations indicated that the optimal dosing regimen for patients across the three populations of European, Japanese, and Chinese ancestry was OSI 80 mg once daily (OD) to achieve the desired range of plasma Ctrough (328-677 nmol/L), as well as 80 mg and 160 mg OD for desirable pulmonary EGFRm+ inhibition (>80%). Conclusion: In conclusion, this study's PBPK simulations highlighted potential ethnic and inter-patient variability in OSI PK and EGFRm+ inhibition between Caucasian, Japanese, and Chinese populations, while also providing insights into optimal dosing regimens of OSI.

Prediction of ROS1 and TRKA/B/C occupancy in plasma and cerebrospinal fluid for entrectinib alone and in DDIs using physiologically based pharmacokinetic (PBPK) modeling approach

PURPOSE

Entrectinib (ENT) is a potent c-ros oncogene 1(ROS1) and neurotrophic tyrosine receptor kinase (NTRKA/B/C) inhibitor. To determine the optimum dosage of ENT using ROS1 and NTRKA/B/C occupancy in plasma and cerebrospinal fluid (CSF) in drug-drug interactions (DDIs), physiologically-based pharmacokinetic (PBPK) models for healthy subjects and cancer population were developed for ENT and M5 (active metabolite).

METHODS

The PBPK models were built using the modeling parameters of ENT and M5 that were mainly derived from the published paper on the ENT PBPK model, and then validated by the observed pharmacokinetics (PK) in plasma and CSF from healthy subjects and patients.

RESULTS

The PBPK model showed that AUC, Cmax, and Ctrough ratios between predictions and observations are within the range of 0.5-2.0, except that the M5 AUC ratio is slightly above 2.0 (2.34). Based on the efficacy (> 75% occupancy for ROS1 and NTRKA/B/C) and safety (AUC < 160 μM·h and Cmax < 8.9 μM), the appropriate dosing regimens were identified. The appropriate dosage is 600 mg once daily (OD) when administered alone, reduced to 200 mg and 400 mg OD with itraconazole and fluconazole, respectively. ENT is not recommended for co-administration with rifampicin or efavirenz, but is permitted with fluvoxamine or dexamethasone.

CONCLUSION

The PBPK models can serve as a powerful approach to predict ENT concentration as well as ROS1 and NTRKA/B/C occupancy in plasma and CSF.

Catalytic Degraders Effectively Address Kinase Site Mutations in EML4-ALK Oncogenic Fusions

Heterobifunctional degraders, known as proteolysis targeting chimeras (PROTACs), theoretically possess a catalytic mode-of-action, yet few studies have either confirmed or exploited this potential advantage of event-driven pharmacology. Degraders of oncogenic EML4-ALK fusions were developed by conjugating ALK inhibitors to cereblon ligands. Simultaneous optimization of pharmacology and compound properties using ternary complex modeling and physicochemical considerations yielded multiple catalytic degraders that were more resilient to clinically relevant ATP-binding site mutations than kinase inhibitor drugs. Our strategy culminated in the design of the orally bioavailable derivative CPD-1224 that avoided hemolysis (a feature of detergent-like PROTACs), degraded the otherwise recalcitrant mutant L1196M/G1202R in vivo, and commensurately slowed tumor growth, while the third generation ALK inhibitor drug lorlatinib had no effect. These results validate our original therapeutic hypothesis by exemplifying opportunities for catalytic degraders to proactively address binding site resistant mutations in cancer.

In vivo-in vitro correlation of antitumor activity of heat shock protein 90 (HSP90) inhibitors with a pharmacokinetics/pharmacodynamics analysis using NCI-N87 xenograft mice

The in vitro antitumor activity (e.g. IC₅₀) of anticancer drugs is important for selecting candidate compounds for in vivo drug efficacy study in the early stage of drug discovery. In this study, we investigated the relationship between in vitro IC₅₀ and in vivo EC₅₀ using six heat shock protein 90 (HSP90) inhibitors.IC₅₀ of each compound was calculated from in vitro cell proliferation assays using the NCI-N87 cancer cell line. Each compound was administered to NCI-N87 xenograft mice, and EC₅₀ and the maximum tumour-killing rate constant were calculated from pharmacokinetics/pharmacodynamics analyses using plasma concentrations and tumour volumes.IC₅₀ obtained in vitro was poorly correlated with EC₅₀ obtained in vivo, while a good correlation (r = 0.856) was observed between them when corrected with the unbound fraction ratio.The results of this study using of HSP90 inhibitors as model compounds suggest importance of the consideration of an unbound fraction to evaluate the relationship between IC₅₀ and EC₅₀. These results will contribute to improvement in the prediction accuracy of in vivo drug efficacy from in vitro activity and the efficiency of drug discovery research.

A pharmacokinetic-pharmacodynamic model for the MET tyrosine kinase inhibitor, savolitinib, to explore target inhibition requirements for anti-tumour activity

BACKGROUND AND PURPOSE

Savolitinib (AZD6094, HMPL-504, volitinib) is an oral, potent, and highly MET receptor TK inhibitor. This series of studies aimed to develop a pharmacokinetic-pharmacodynamic (PK/PD) model to link inhibition of MET phosphorylation (pMET) by savolitinib with anti-tumour activity.

EXPERIMENTAL APPROACH

Cell line-derived xenograft (CDX) experiments using human lung cancer (EBC-1) and gastric cancer (MKN-45) cells were conducted in athymic nude mice using a variety of doses and schedules of savolitinib. Tumour pMET changes and growth inhibition were calculated after 28 days. Population PK/PD techniques were used to construct a PK/PD model for savolitinib.

KEY RESULTS

Savolitinib showed dose- and dose frequency-dependent anti-tumour activity in the CDX models, with more frequent, lower dosing schedules (e.g., twice daily) being more effective than intermittent, higher dosing schedules (e.g., 4 days on/3 days off or 2 days on/5 days off). There was a clear exposure-response relationship, with maximal suppression of pMET of >90%. Data from additional CDX and patient-derived xenograft (PDX) models overlapped, allowing calculation of a single EC₅₀ of 0.38 ng·ml^-1 . Tumour growth modelling demonstrated that prolonged, high levels of pMET inhibition (>90%) were required for tumour stasis and regression in the models.

CONCLUSION AND IMPLICATIONS

High and persistent levels of MET inhibition by savolitinib were needed for optimal monotherapy anti-tumour activity in preclinical models. The modelling framework developed here can be used to translate tumour growth inhibition from the mouse to human and thus guide choice of clinical dose and schedule.

Application of Pharmacokinetic-Pharmacodynamic Modeling in Drug Delivery: Development and Challenges

With the advancement of technology, drug delivery systems and molecules with more complex architecture are developed. As a result, the drug absorption and disposition processes after administration of these drug delivery systems and engineered molecules become exceedingly complex. As the pharmacokinetic and pharmacodynamic (PK-PD) modeling allows for the separation of the drug-, carrier- and pharmacological system-specific parameters, it has been widely used to improve understanding of the in vivo behavior of these complex delivery systems and help their development. In this review, we summarized the basic PK-PD modeling theory in drug delivery and demonstrated how it had been applied to help the development of new delivery systems and modified large molecules. The linkage between PK and PD was highlighted. In particular, we exemplified the application of PK-PD modeling in the development of extended-release formulations, liposomal drugs, modified proteins, and antibody-drug conjugates. Furthermore, the model-based simulation using primary PD models for direct and indirect PD responses was conducted to explain the assertion of hypothetical minimal effective concentration or threshold in the exposure-response relationship of many drugs and its misconception. The limitations and challenges of the mechanism-based PK-PD model were also discussed.

Does Crizotinib Auto-Inhibit CYP3A in vivo?

Crizotinib is a tyrosine kinase inhibitor used to treat anaplastic lymphoma kinase-positive lung cancer. There is in vitro evidence that crizotinib may auto-inhibit cytochrome P450 3A (CYP3A) activity, with important implications for crizotinib pharmacokinetics. In order to test whether crizotinib treatment alters CYP3A activity in vivo, mice were treated with 5 and 25 mg/kg crizotinib (p.o.) daily for 14 days. Results showed that crizotinib treatment did not alter CYP3A activity as determined by erythromycin N-demethylation. In addition, CYP3A polypeptide expression as measured by Western blot was unchanged. Therefore, our results do not support CYP3A inhibition by crizotinib in vivo.

Translational prediction of first-in-human pharmacokinetics and pharmacodynamics of janagliflozin, a selective SGLT2 inhibitor, using allometric scaling, dedrick and PK/PD modeling methods

AIM

Janagliflozin is an orally selective SGLT2 inhibitor. To predict human pharmacokinetics/pharmacodynamics (PK/PD) characteristics of janagliflozin. To design optimal starting dose and effective dose for janagliflozin first-in-human (FIH) study.

METHODS

Animal PK/PD properties of janagliflozin were obtained from preclinical in vivo and in vitro study. Pharmacologically effective level of same class SGLT2 inhibitors were assessed through preclinical and clinical efficacy data of dapagliflozin, empagliflozin and canagliflozin. Human PK parameters and profiles of janagliflozin were predicted by various methods such as allometric scaling (AS), dedrick and PK/PD modeling analysis. Mechanistic PK/PD model was developed to describe janagliflozin-mediated impact on urinary glucose excretion (UGE). Human IC₅₀ was scaled from rat model-estimated IC₅₀ by correcting interspecies difference of in vitro IC₅₀ and plasma f_u of rat and human. The quantitative PK/PD prediction of janagliflozin was evaluated via observed PK/PD profiles of healthy subjects. Predicted PK/PD characteristics of janagliflozin were applied in FIH dose design. Optimal starting dose was suggested by considering preclinical PD and toxicity data of janagliflozin. Effective dose was suggested by considering pharmacologically effective level of same class drugs.

RESULTS

PK/PD characteristics of janagliflozin in preclinical species were summarized. Pharmacologically effective level for SGLT2 inhibitors was defined as 25~30% ΔUGE (ΔUGE=--(PG*GFR)_within24h) based on efficacy data of three same class drugs. Human predicted CL, V_ss and F were 1.04 L/h, 77.5 L and 0.80. Predicted AUC and C_max of janagliflozin of 10 and 50 mg were within 0.47~2.08 fold of observed values. Predicted human UGE_0-24 h and UGE_0-144 h of 10 and 50 mg dose range were within 0.66~1.41 fold of observed values. Optimal starting dose and pharmacologically active dose (PAD) were suggested as 10 mg and 50 mg. Dose range for FIH study was designed as 10-450 mg.

CONCLUSIONS

This study predicted human PK/PD characteristics of janagliflozin based on preclinical data and provide optimal dose design for janagliflozin FIH study based on pharmacologically effective level of same class drugs.

Pharmacokinetic-Pharmacodynamic-Efficacy Modeling of ONO-7579, a Novel Pan-Tropomyosin Receptor Kinase Inhibitor, in a Murine Xenograft Tumor Model

The orally available and novel small molecule ONO-7579 (N-{2-[4-(2-amino-5-chloropyridin-3-yl)phenoxy]pyrimidin-5-yl}-N'-[2-(methanesulfonyl)-5-(trifluoromethyl)phenyl]urea) is a highly potent and selective pan-tropomyosin receptor kinase (TRK) inhibitor. The objective of the present study was to characterize the pharmacokinetic (PK), pharmacodynamic (PD), and antitumor efficacy relationships of ONO-7579 in mice xenografted with a human colorectal cancer cell line, KM12 (harboring the tropomyosin 3 (TPM3) -neurotrophic tyrosine receptor kinase 1 fusion gene), via a PK/PD modeling approach. Plasma and tumor concentrations of ONO-7579, tumor levels of phosphorylated TPM3-TRKA (pTRKA), and tumor volumes in the murine model were measured with a single or multiple dose of ONO-7579 (0.06-0.60 mg/kg) administered once daily. The PK/PD/efficacy models were developed in a sequential manner. Changes in plasma concentrations of ONO-7579 were described with an oral one-compartment model. Tumor concentrations of ONO-7579 were higher than plasma concentrations, and changes in ONO-7579 tumor concentrations were described with an additional tumor compartment that had no influence on plasma concentrations. pTRKA in tumors was described with a direct E_max model, and the tumor ONO-7579 concentration causing 50% of the maximum effect was estimated to be 17.6 ng/g. In addition, a pTRKA-driven tumor growth inhibition model indicated that ONO-7579 started to sharply increase the antitumor effect at pTRKA inhibition rates >60% and required >91.5% to reduce tumors. In conclusion, the developed PK/PD/efficacy models revealed a "switch-like" relationship between pTRKA inhibition rate and antitumor effect in a murine KM12 xenograft model, demonstrating that pTRKA in tumors could serve as an effective biomarker for scheduling the dose regimen in early-stage clinical studies. SIGNIFICANCE STATEMENT: In recent years, clinical development of TRK inhibitors in patients with neurotrophic tyrosine receptor kinase fusion-positive solid tumors has been accelerated. This research found that phosphorylated TRKA was a useful biomarker for explaining the antitumor efficacy of TRK inhibitors using a pharmacokinetic/pharmacodynamic modeling approach in xenograft mice. This finding suggests a rational dosing regimen in early-stage clinical studies for ONO-7579 (N-{2-[4-(2-amino-5-chloropyridin-3-yl)phenoxy]pyrimidin-5-yl}-N'-[2-(methanesulfonyl)-5-(trifluoromethyl)phenyl]urea), a novel pan-TRK inhibitor.

A Microfluidic Perfusion Platform for In Vitro Analysis of Drug Pharmacokinetic-Pharmacodynamic (PK-PD) Relationships

Static in vitro cell culture studies cannot capture the dynamic concentration profiles of drugs, nutrients, and other factors that cells experience in physiological systems. This limits the confidence in the translational relevance of in vitro experiments and increases the reliance on empirical testing of exposure-response relationships and dose optimization in animal models during preclinical drug development, introducing additional challenges owing to species-specific differences in drug pharmacokinetics (PK) and pharmacodynamics (PD). Here, we describe the development of a microfluidic cell culture device that enables perfusion of cells under 2D or 3D culture conditions with temporally programmable concentration profiles. Proof-of-concept studies using doxorubicin and gemcitabine demonstrated the ability of the microfluidic PK-PD device to examine dose- and time-dependent effects of doxorubicin as well as schedule-dependent effects of doxorubicin and gemcitabine combination therapy on cell viability using both step-wise drug concentration profiles and species-specific (i.e., mouse, human) drug PK profiles. The results demonstrate the importance of including physiologically relevant dynamic drug exposure profiles during in vitro drug testing to more accurately mimic in vivo drug effects, thereby improving translatability across nonclinical studies and reducing the reliance on animal models during drug development.

Translational Pharmacokinetic-Pharmacodynamic Modeling for an Orally Available Novel Inhibitor of Epigenetic Regulator Enhancer of Zeste Homolog 2

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 EC ₅₀ of 76 nM, whereas a signal-transduction model sufficiently fit dose-dependent disease responses (tumor growth inhibition) with an unbound tumor stasis concentration (T _sc ) of 168 nM. Thus, effective concentration for 70% of maximal effect (EC₇₀) for H3K27me3 inhibition was roughly comparable to T _sc , 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 EC₇₀ 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.

Preclinical pharmacokinetics, disposition, and translational pharmacokinetic/pharmacodynamic modeling of savolitinib, a novel selective cMet inhibitor

Savolitinib is a novel small-molecule selective cMet inhibitor. This work characterized its pharmacokinetics in preclinical phase, established the preclinical relationships between PK, cMet modulation and anti-tumor efficacy. In vitro and in vivo animal studies were performed for PK characterization. Savolitinib showed good absorption, moderate tissue distribution, low to intermediate clearance, and low accumulation. Hepatic oxidative metabolism followed by urinary and biliary excretions was the major elimination pathway. Based on preclinical PK data, human PK profiles were predicted using empirical methods. Pharmacodynamic studies for evaluating cMet inhibition and anti-tumor efficacy were conducted in nude mice bearing Hs746t xenograft. PK/PD models were built to link the PD measurements to nude mouse PK. The established integrated preclinical PK/PD model contained a two-compartment non-linear PK model, a biomarker link model and a tumor growth transit model. The IC₅₀ of cMet inhibition and the concentration achieving half of the maximal Hs746t tumor reduction by savolitinib were equal to 12.5 and 3.7 nM (free drug), respectively. Based on the predicted human PK data, as well as the established PK/PD model in nude mouse, the human PD (cMet inhibition) profiles were also simulated. This research supported clinical development of savolitinib. Understanding the preclinical PK/PD relationship of savolitinib provides translational insights into the cMet-targeted drug development.

Current status and future perspective on preclinical pharmacokinetic and pharmacodynamic (PK/PD) analysis: Survey in Japan pharmaceutical manufacturers association (JPMA)

Preclinical pharmacokinetic/pharmacodynamic (PK/PD) analysis is an efficient tool for the translational research and proof of mechanism/concept in animals. The questionnaire survey on the practice of preclinical PK/PD analysis was conducted in the member companies of the Japan Pharmaceutical Manufacturers Association (JPMA). According to the survey, 60% of companies conducted preclinical PK/PD analysis and its impact for drug development was different between each of the companies. The frequently analyzed therapeutic areas of preclinical PK/PD analysis were neurology, inflammation and metabolic disease, and those are different from the therapeutic area (infectious disease and oncology) in which PK/PD analysis was considered as effective by the present survey. Many companies which have used preclinical PK/PD analysis for the translation to human PK/PD and for the prediction of dose/regimen had good communication with other research & development (R&D) departments (e.g. pharmacology/clinical pharmacology). The increase in resources for preclinical PK/PD analysis including education was highly demanded. As a future perspective, the closer collaboration between pharmacokinetics scientists, pharmacologists, toxicologists and clinical pharmacologists and the increase in resources including upskilling and the comprehension of preclinical PK/PD analysis by the project team are considered to lead to efficient contributions to improve the success ratio of drug discovery and development.

Molecular Pharmacodynamics-Guided Scheduling of Biologically Effective Doses: A Drug Development Paradigm Applied to MET Tyrosine Kinase Inhibitors

The development of molecularly targeted agents has benefited from use of pharmacodynamic markers to identify "biologically effective doses" (BED) below MTDs, yet this knowledge remains underutilized in selecting dosage regimens and in comparing the effectiveness of targeted agents within a class. We sought to establish preclinical proof-of-concept for such pharmacodynamics-based BED regimens and effectiveness comparisons using MET kinase small-molecule inhibitors. Utilizing pharmacodynamic biomarker measurements of MET signaling (tumor pY1234/1235MET/total MET ratio) in a phase 0-like preclinical setting, we developed optimal dosage regimens for several MET kinase inhibitors and compared their antitumor efficacy in a MET-amplified gastric cancer xenograft model (SNU-5). Reductions in tumor pY1234/1235MET/total MET of 95%-99% were achievable with tolerable doses of EMD1214063/MSC2156119J (tepotinib), XL184 (cabozantinib), and XL880/GSK1363089 (foretinib), but not ARQ197 (tivantinib), which did not alter the pharmacodynamic biomarker. Duration of kinase suppression and rate of kinase recovery were specific to each agent, emphasizing the importance of developing customized dosage regimens to achieve continuous suppression of the pharmacodynamic biomarker at the required level (here, ≥90% MET kinase suppression). The customized dosage regimen of each inhibitor yielded substantial and sustained tumor regression; the equivalent effectiveness of customized dosage regimens that achieve the same level of continuous molecular target control represents preclinical proof-of-concept and illustrates the importance of proper scheduling of targeted agent BEDs. Pharmacodynamics-guided biologically effective dosage regimens (PD-BEDR) potentially offer a superior alternative to pharmacokinetic guidance (e.g., drug concentrations in surrogate tissues) for developing and making head-to-head comparisons of targeted agents. Mol Cancer Ther; 17(3); 698-709. ©2018 AACR.

Prediction of a Therapeutic Dose for Buagafuran, a Potent Anxiolytic Agent by Physiologically Based Pharmacokinetic/Pharmacodynamic Modeling Starting from Pharmacokinetics in Rats and Human

Physiologically based pharmacokinetic (PBPK)/pharmacodynamic (PD) models can contribute to animal-to-human extrapolation and therapeutic dose predictions. Buagafuran is a novel anxiolytic agent and phase I clinical trials of buagafuran have been completed. In this paper, a potentially effective dose for buagafuran of 30 mg t.i.d. in human was estimated based on the human brain concentration predicted by a PBPK/PD modeling. The software GastroPlus^TM was used to build the PBPK/PD model for buagafuran in rat which related the brain tissue concentrations of buagafuran and the times of animals entering the open arms in the pharmacological model of elevated plus-maze. Buagafuran concentrations in human plasma were fitted and brain tissue concentrations were predicted by using a human PBPK model in which the predicted plasma profiles were in good agreement with observations. The results provided supportive data for the rational use of buagafuran in clinic.

Practical Recommendations for Therapeutic Drug Monitoring of Kinase Inhibitors in Oncology

Despite the fact that pharmacokinetic exposure of kinase inhibitors (KIs) is highly variable and clear relationships exist between exposure and treatment outcomes, fixed dosing is still standard practice. This review aims to summarize the available clinical pharmacokinetic and pharmacodynamic data into practical guidelines for individualized dosing of KIs through therapeutic drug monitoring (TDM). Additionally, we provide an overview of prospective TDM trials and discuss the future steps needed for further implementation of TDM of KIs.

Current mathematical models for cancer drug discovery

INTRODUCTION

Pharmacometric models represent the most comprehensive approaches for extracting, summarizing and integrating information obtained in the often sparse, limited, and less-than-optimally designed experiments performed in the early phases of oncology drug discovery. Whilst empirical methodologies may be enough for screening and ranking candidate drugs, modeling approaches are needed for optimizing and making economically viable the learn-confirm cycles within an oncology research program and anticipating the dose regimens to be investigated in the subsequent clinical development. Areas covered: Papers appearing in the literature of approximately the last decade reporting modeling approaches applicable to anticancer drug discovery have been listed and commented. Papers were selected based on the interest in the proposed methodology or in its application. Expert opinion: The number of modeling approaches used in the discovery of anticancer drugs is consistently increasing and new models are developed based on the current directions of research of new candidate drugs. These approaches have contributed to a better understanding of new oncological targets and have allowed for the exploitation of the relatively sparse information generated by preclinical experiments. In addition, they are used in translational approaches for guiding and supporting the choice of dosing regimens in early clinical development.

Enhancing value of clinical pharmacodynamics in oncology drug development: An alliance between quantitative pharmacology and translational science

Clinical pharmacodynamic evaluation is a key component of the "pharmacologic audit trail" in oncology drug development. We posit that its value can and should be greatly enhanced via application of a robust quantitative pharmacology framework informed by biologically mechanistic considerations. Herein, we illustrate examples of intersectional blindspots across the disciplines of quantitative pharmacology and translational science and offer a roadmap aimed at enhancing the caliber of clinical pharmacodynamic research in the development of oncology therapeutics.

Translational modeling and simulation approaches for molecularly targeted small molecule anticancer agents from bench to bedside

INTRODUCTION

Recent advances in molecular biology have enabled personalized cancer therapies with molecularly targeted agents (MTAs), which offer a promising future for cancer therapy. Dynamic modeling and simulation (M&S) is a powerful mathematical approach linking drug exposures to pharmacological responses, providing a quantitative assessment of in vivo drug potency. Accordingly, a growing emphasis is being placed upon M&S to quantitatively understand therapeutic exposure-response relationships of MTAs in nonclinical models.

AREAS COVERED

An overview of M&S approaches for MTAs in nonclinical models is presented with discussion about mechanistic extrapolation of antitumor efficacy from bench to bedside. Emphasis is placed upon recent advances in M&S approaches linking drug exposures, biomarker responses (e.g. target modulation) and pharmacological outcomes (e.g. antitumor efficacy).

EXPERT OPINION

For successful personalized cancer therapies with MTAs, it is critical to mechanistically and quantitatively understand their exposure-response relationships in nonclinical models, and to logically and properly apply such knowledge to the clinic. Particularly, M&S approaches to predict pharmacologically active concentrations of MTAs in patients based upon nonclinical data would be highly valuable in guiding the design and execution of clinical trials. Proactive approaches to understand their exposure-response relationships could substantially increase probability of achieving a positive proof-of-concept in the clinic.

Efficacy Screening of Gloriosa Superba Extracts in a Murine Pancreatic Cancer Model Using (18)F-FDG PET/CT for Monitoring Treatment Response

PURPOSE

In vivo efficacy of two herbal extracts of Gloriosa superba L. (Colchicaceae) was investigated in a murine pancreatic tumor model by tumor volume measurements and Positron Emission Tomography (PET) imaging using 2-deoxy-2-[(18)F]fluoro-d-glucose ((18)F-FDG).

MATERIALS AND METHODS

A crude extract of G. superba (GS) seeds rich in colchicine and a colchicine-poor extract (GS2B) containing mostly colchicoside as a putative prodrug were prepared. PANC02-bearing C57BL/6 mice were treated with either placebo, gemcitabine, or one of the extracts (three different doses) for 10 days. Tumor volume measurements were performed daily during treatment and additionally (18)F-FDG Positron emission tomography/computed tomography was acquired at baseline and after 7 days of treatment. Ki-67 and cleaved caspase-3 immunostaining was performed on the resected tumors.

RESULTS

After 7 days of treatment, a dose-dependent tumor growth inhibition of both extracts was observed with the highest in vivo response at the highest dose of GS and GS2B and gemcitabine. A positive significant correlation was found between Ki-67 scores and relative tumor volumes (RTV), and a negative significant correlation between caspase-3 staining scores and RTV. A decrease in (18)F-FDG uptake was clearly observed in all treatment groups.

CONCLUSIONS

The therapeutic efficacy of the two different herbal extracts was demonstrated in an in vivo pancreatic tumor model. (18)F-FDG PET was able to detect an early response as overall lower (18)F-FDG uptake was measured in the treated groups.

Mixture dynamics: Combination therapy in oncology

In recent years combination therapies have become increasingly popular in most therapeutic areas. We present a qualitative and quantitative approach and elucidate some of the challenges and solutions to a more optimal therapy. For tumor growth this involves the study of semi-mechanistic cell-growth/kill models with multiple sites of action. We introduce such models and analyze their dynamic properties using simulations and mathematical analysis. This is done for two specific case studies, one involving a single compound and one a combination of two compounds. We generalize the notion of Tumor Static Concentration to cases when two compounds are involved and develop a graphical method for determining the optimal combination of the two compounds, using ideas akin to those used in studies employing isobolograms. In studying the dynamics of the second case study we focus, not only on the different concentrations, but also on the different dosing regimens and pharmacokinetics of the two compounds.

Pharmacokinetics in Drug Discovery: An Exposure-Centred Approach to Optimising and Predicting Drug Efficacy and Safety

The role of pharmacokinetics (PK) in drug discovery is to support the optimisation of the absorption, distribution, metabolism and excretion (ADME) properties of lead compounds with the ultimate goal to attain a clinical candidate which achieves a concentration-time profile in the body that is adequate for the desired efficacy and safety profile. A thorough characterisation of the lead compounds aiming at the identification of the inherent PK liabilities also includes an early generation of PK/PD relationships linking in vitro potency and target exposure/engagement with expression of pharmacological activity (mode-of-action) and efficacy in animal studies. The chapter describes an exposure-centred approach to lead generation, lead optimisation and candidate selection and profiling that focuses on a stepwise generation of an understanding between PK/exposure and PD/efficacy relationships by capturing target exposure or surrogates thereof and cellular mode-of-action readouts in vivo. Once robust PK/PD relationship in animal PD models has been constructed, it is translated to anticipate the pharmacologically active plasma concentrations in patients and the human therapeutic dose and dosing schedule which is also based on the prediction of the PK behaviour in human as described herein. The chapter outlines how the level of confidence in the predictions increases with the level of understanding of both the PK and the PK/PD of the new chemical entities (NCE) in relation to the disease hypothesis and the ability to propose safe and efficacious doses and dosing schedules in responsive patient populations. A sound identification of potential drug metabolism and pharmacokinetics (DMPK)-related development risks allows proposing of an effective de-risking strategy for the progression of the project that is able to reduce uncertainties and to increase the probability of success during preclinical and clinical development.

Deciphering the In Vivo Performance of a Monoclonal Antibody to Neutralize Its Soluble Target at the Site of Action in a Mouse Collagen-Induced Arthritis Model

PURPOSE

Study the disposition and target-neutralization capability of an anti-interleukin-6 (IL-6) monoclonal antibody (mAb) at the joint in a mouse collagen-induced arthritis (CIA) model.

METHODS

A mechanistic pharmacokinetic/pharmacodynamic study was conducted in a mouse CIA model using CNTO 345, a rat anti-mouse IL-6 mAb, as model compound. The drug, total/free IL-6 concentrations in both serum and joint lavage fluid were quantitatively assessed and compared to those in the normal control mice.

RESULTS

CNTO 345 exhibited higher clearance and significantly higher joint lavage/serum ratio in the CIA mice than in the normal control mice. The mAb concentrations in the joint lavage are approximately proportional to the serum concentrations at all the time points being examined. Dosing of CNTO 345 led to sustained free IL-6 suppression in both serum and joint lavage in a dose-dependent manner. A dose-dependent increase in total IL-6 was observed in serum, but not in the joint lavage fluid. Though no change in disease activity was observed following a single dose of anti-IL-6 mAb at peak of the disease, a dose-dependent decrease in serum amyloid A, a downstream biomarker of IL-6, was observed.

CONCLUSIONS

This study provided quantitative assessments of the distribution and target-neutralization capability of an anti-IL-6 mAb at the site of action in an animal disease model.

Optimization of human dose prediction by using quantitative and translational pharmacology in drug discovery

In this perspective article, we explain how quantitative and translational pharmacology, when well-implemented, is believed to lead to improved clinical candidates and drug targets that are differentiated from current treatment options. Quantitative and translational pharmacology aims to build and continuously improve the quantitative relationship between drug exposure, target engagement, efficacy, safety and its interspecies relationship at every phase of drug discovery. Drug hunters should consider and apply these concepts to develop compounds with a higher probability of interrogating the clinical biological hypothesis. We offer different approaches to set an initial effective concentration or pharmacokinetic-pharmacodynamic target in man and to predict human pharmacokinetics that determine together the predicted human dose and dose schedule. All concepts are illustrated with ample literature examples.

Model-Based Design of a Decision Tree for Treating HER2+ Cancers Based on Genetic and Protein Biomarkers

Human cancers are incredibly diverse with regard to molecular aberrations, dependence on oncogenic signaling pathways, and responses to pharmacological intervention. We wished to assess how cellular dependence on the canonical PI3K vs. MAPK pathways within HER2+ cancers affects responses to combinations of targeted therapies, and biomarkers predictive of their activity. Through an integrative analysis of mechanistic model simulations and in vitro cell line profiling, we designed a six-arm decision tree to stratify treatment of HER2+ cancers using combinations of targeted agents. Activating mutations in the PI3K and MAPK pathways (PIK3CA and KRAS), and expression of the HER3 ligand heregulin determined sensitivity to combinations of inhibitors against HER2 (lapatinib), HER3 (MM-111), AKT (MK-2206), and MEK (GSK-1120212; trametinib), in addition to the standard of care trastuzumab (Herceptin). The strategy used to identify effective combinations and predictive biomarkers in HER2-expressing tumors may be more broadly extendable to other human cancers.

Pharmacokinetics of crizotinib in NSCLC patients

INTRODUCTION

For a subpopulation of NSCLC patients genetic rearrangement of the anaplastic lymphoma kinase (ALK) was found as driver mutation, which can be targeted by the selective inhibitor crizotinib.

AREAS COVERED

This article presents an overview of the clinical studies that provided the characterization of the pharmacokinetic parameters for the administration of crizotinib to cancer patients and the factors influencing the clinical profiles of this drug.

EXPERT OPINION

Crizotinib is administered orally as a capsule and clinical studies indicated 250 mg crizotinib BID continuously as the maximal tolerated dose in cancer patients. Bioavailability is ∼ 40% and pharmacokinetic parameters are influenced by food only to a minor degree. This dose of the drug corresponds to a significant inhibition of the mutated ALK, retards tumor growth and achieves clinical responses in the majority of patients. Crizotinib lactam is the single metabolite with minor inhibitory activity for the ALK fusion protein. Metabolization is executed mainly by CYP3A4/5 and is modulated by other drugs interacting with this cytochrome oxidase. Despite the one-fits-all approach in administration of crizotinib at a fixed dose the pharmacokinetic parameters indicate a stable steady state upon continuous administration, which achieves sufficient inhibition of the ALK drug target.

Mechanistic understanding of translational pharmacokinetic-pharmacodynamic relationships in nonclinical tumor models: a case study of orally available novel inhibitors of anaplastic lymphoma kinase

The orally available novel small molecules PF06463922 [(10R)-7-amino-12-fluoro-2,10,16-trimethyl-15-oxo-10,15,16,17-tetrahydro-2H-8,4-(metheno)pyrazolo[4,3-h][2,5,11]benzoxadiazacyclotetradecine-3-carbonitrile] and PF06471402 [(10R)-7-amino-12-fluoro-2,10,16-trimethyl-15-oxo-10,15,16,17-tetrahydro-2H-8,4-(azeno)pyrazolo[4,3-h][2,5,11]benzoxadiazacyclo-tetradecine-3-carbonitrile] are second-generation anaplastic lymphoma kinase (ALK) inhibitors targeted to both naïve and resistant patients with non-small cell lung cancer (NSCLC) to the first-generation ALK inhibitor crizotinib. The objectives of the present study were to characterize and compare the pharmacokinetic-pharmacodynamic (PKPD) relationships of PF06463922 and PF06471402 for target modulation in tumor and antitumor efficacy in athymic mice implanted with H3122 NSCLC cells expressing a crizotinib-resistant echinoderm microtubule-associated protein-like 4 (EML4)-ALK mutation, EML4-ALK(L1196M). Furthermore, the PKPD relationships for these ALK inhibitors were evaluated and compared between oral administration and subcutaneous constant infusion (i.e., between different pharmacokinetic [PK] profiles). Oral and subcutaneous PK profiles of these ALK inhibitors were adequately described by a one-compartment PK model. An indirect response model extended with a modulator fit the time courses of PF06463922- and PF06471402-mediated target modulation (i.e., ALK phosphorylation) with an estimated unbound EC50,in vivo of 36 and 20 nM, respectively, for oral administration, and 100 and 69 nM, respectively, for subcutaneous infusion. A drug-disease model based on the turnover concept fit tumor growth curves inhibited by PF06463922 and PF06471402 with estimated unbound tumor stasis concentrations of 51 and 27 nM, respectively, for oral administration, and 116 and 70 nM, respectively, for subcutaneous infusion. Thus, the EC50,in vivo to EC60,in vivo estimates for ALK inhibition corresponded to the concentrations required tumor stasis in all cases, suggesting that the pharmacodynamic relationships of target modulation to antitumor efficacy were consistent among the ALK inhibitors, even when the PK profiles with different administration routes were considerably different.

Translation of anticancer efficacy from nonclinical models to the clinic

Mouse cancer models have provided critical insights into tumor biology; however, clinical translation of these findings has been challenging. This perspective posits that factors impacting on successful translation start with limitations in capturing human cancer pathophysiology and end with challenges in generating robust translatable preclinical end points. A comprehensive approach that considers clinically relevant mouse models with both an integrated biomarker strategy and a complementary modeling and simulation effort will strengthen the current oncology drug development paradigm.

Translational pharmacokinetic-pharmacodynamic modeling for an orally available novel inhibitor of anaplastic lymphoma kinase and c-Ros oncogene 1

An orally available macrocyclic small molecule, PF06463922 [(10R)-7-amino-12-fluoro-2,10,16-trimethyl-15-oxo-10,15,16,17-tetrahydro-2H-8,4-(metheno)pyrazolo[4,3-h][2,5,11]benzoxadiazacyclotetradecine-3-carbonitrile], is a selective inhibitor of anaplastic lymphoma kinase (ALK) and c-Ros oncogene 1 (ROS1). The objectives of the present study were to characterize the pharmacokinetic-pharmacodynamic relationships of PF06463922 between its systemic exposures, pharmacodynamic biomarker (target modulation), and pharmacologic response (antitumor efficacy) in athymic mice implanted with H3122 non-small cell lung carcinomas expressing echinoderm microtubule-associated protein-like 4 (EML4)-ALK mutation (EML4-ALK(L1196M)) and with NIH3T3 cells expressing CD74-ROS1. In these nonclinical tumor models, PF06463922 was orally administered to animals with EML4-ALK(L1196M) and CD74-ROS1 at twice daily doses of 0.3-20 and 0.01-3 mg/kg per dose, respectively. Plasma concentration-time profiles of PF06463922 were adequately described by a one-compartment pharmacokinetic model. Using the model-simulated plasma concentrations, a pharmacodynamic indirect response model with a modulator sufficiently fit the time courses of target modulation (i.e., ALK phosphorylation) in tumors of EML4-ALK(L1196M)-driven models with EC50,in vivo of 36 nM free. A drug-disease model based on an indirect response model reasonably fit individual tumor growth curves in both EML4-ALK(L1196M)- and CD74-ROS1-driven models with the estimated tumor stasis concentrations of 51 and 6.2 nM free, respectively. Thus, the EC60,in vivo (52 nM free) for ALK inhibition roughly corresponded to the tumor stasis concentration in an EML4-ALK(L1196M)-driven model, suggesting that 60% ALK inhibition would be required for tumor stasis. Accordingly, we proposed that the EC60,in vivo for ALK inhibition corresponding to the tumor stasis could be considered a minimum target efficacious concentration of PF06463922 for cancer patients in a phase I trial.

Pharmacokinetics of Peptide-Fc fusion proteins

Peptide-Fc fusion proteins (or peptibodies) are chimeric proteins generated by fusing a biologically active peptide with the Fc-domain of immunoglobulin G. In this review, we describe recent studies that have evaluated the absorption, distribution, metabolism, and excretion characteristics of peptibodies. Key features of the pharmacokinetics of peptibodies include their extended half-life due to recycling by the neonatal Fc receptor (FcRn), a substantial contribution by renal excretion to total clearance and, for certain peptibodies, target-mediated drug disposition. The prolonged half-life of peptibodies permits less-frequent dose administration compared with small therapeutic peptides, thereby supporting patient convenience and compliance. Hence, a considerable number of peptibodies are currently in preclinical and clinical development. Investigation of the metabolism (biotransformation) of biologics is an evolving area of research: ligand-binding mass spectrometry techniques have been employed for the characterization of the peptibody romiplostim, providing a new approach to evaluation of the degradation products of biologics. Pharmacokinetic/pharmacodynamic modeling and simulation techniques have been used to predict the pharmacokinetics of peptibodies which can inform clinical decision-making, particularly selection of dosing regimens. This integrated review highlights the distinct pharmacokinetic characteristics of peptibodies and their influence on the drug development process for this emerging family of therapeutics.