A translational model-based approach to inform the choice of the dose in phase 1 oncology trials: the case study of erdafitinib

PURPOSE

Erdafitinib (JNJ-42756493, BALVERSA) is a tyrosine kinase inhibitor indicated for the treatment of advanced urothelial carcinoma. In this work, a translational model-based approach to inform the choice of the doses in phase 1 trials is illustrated.

METHODS

A pharmacokinetic (PK) model was developed to describe the time course of erdafitinib plasma concentrations in mice and rats. Data from multiple xenograft studies in mice and rats were analyzed using the Simeoni tumor growth inhibition (TGI) model. The model parameters were used to derive a range of erdafitinib exposures that might inform the choice of the doses in oncology phase 1 trials. Conversion of exposures to doses was based on preliminary PK assessments from the first-in human (FIH) study.

RESULTS

A one-compartment PK disposition model, with linear absorption and dose-dependent clearance, adequately described the PK data in both mice and rats via an allometric scaling approach. The TGI model was able to describe tumor growth dynamics, providing quantitative measurements of erdafitinib antitumor potency in mice and rats. Based on these estimates, ranges of efficacious unbound concentration were identified for erdafitinib in mice (0.642-5.364 μg/L) and rats (0.782-2.565 μg/L). Based on the FIH data, it was possible to transpose exposures into doses and doses of above 4 mg/day provided erdafitinib exposures associated with significant TGI in animals. The findings were in agreement with the results of the FIH trial, in which the first hints of clinical activities were observed at 6 mg.

CONCLUSION

The successful modeling exercise of erdafitinib preclinical data showed how translational PK-PD modeling might be a tool to help to inform the choice of the doses in FIH studies.

Population pharmacokinetics of the rilpivirine long-acting formulation after intramuscular dosing in healthy subjects and people living with HIV

OBJECTIVES

To characterize the population pharmacokinetics of the rilpivirine long-acting (LA) formulation after intramuscular administration.

METHODS

Rich and sparse rilpivirine plasma concentration data were obtained from seven clinical studies. In total, 18 261 rilpivirine samples were collected from 986 subjects (131 healthy subjects from Phase I studies and 855 people living with HIV from Phase IIb/III studies). Doses ranged from 300 to 1200 mg, as single-dose or multiple-dose regimens (every 4 or 8 weeks). In Phase III studies, an initiation injection of 900 mg followed by continuation injections of 600 mg every 4 weeks was used. Non-linear mixed-effects modelling was performed using NONMEM® software.

RESULTS

A one-compartment model with linear elimination and two parallel absorption pathways (fast and slow) with sequential zero-first-order processes adequately captured rilpivirine flip-flop pharmacokinetics after intramuscular administration of the LA formulation. The estimated apparent elimination half-life of rilpivirine LA was 200 days. None of the evaluated covariates (age, body weight, BMI, sex, race, health status and needle length) had a clinically relevant impact on rilpivirine pharmacokinetics.

CONCLUSIONS

The population pharmacokinetic model suitably describes the time course and associated variability of rilpivirine plasma concentrations after rilpivirine LA intramuscular administration. The monthly regimen consists of an oral lead-in period (rilpivirine 25 mg tablets once daily for 4 weeks), followed by an initiation injection of 900 mg rilpivirine LA, then 600 mg rilpivirine LA continuation injections monthly. The absence of a clinically relevant effect of covariates on rilpivirine pharmacokinetics suggests that rilpivirine LA dose adjustments for specific subgroups are not warranted.

Assessment of the effect of erdafitinib on cardiac safety: analysis of ECGs and exposure-QTc in patients with advanced or refractory solid tumors

PURPOSE

To characterize the effect of erdafitinib on electrocardiogram (ECG) parameters and the relationship between erdafitinib plasma concentrations and QTc interval changes in patients with advanced or refractory solid tumors.

METHODS

Triplicate ECGs and continuous 12-lead Holter data were collected in the dose escalation part (Part 1) of the first-in-human study, with doses ranging from 0.5 to 12 mg. Triplicate ECG monitoring continued in Parts 2-4 where 2 dose regimens selected from Part 1 were expanded in prespecified tumor types. Analyses of ECG data included central tendency analyses, identification of categorical outliers and morphological assessment. A concentration-QTc analysis was conducted using a linear mixed-effect model based on extracted time matching Holter data.

RESULTS

Central tendency, categorical outlier, and ECG morphologic analyses from 187 patients revealed no clinically significant effect of erdafitinib on heart rate, atrioventricular conduction or cardiac depolarization (PR and QRS), and no effect on cardiac repolarization (QTc). Concentration-QTc analysis from 62 patients indicated that the slopes of relationship between total and free erdafitinib plasma concentrations and QTcI (mean exponent of 0.395) were estimated as - 0.00269 ms/(ng/mL) and - 1.138 ms/(ng/mL), respectively. The predicted change in QTcI at the observed geometric mean of total and free concentration at the highest therapeutic erdafitinib dose (9 mg daily) was < 10 ms at the upper bound of the two-sided 90% confidence interval.

CONCLUSIONS

ECG data and the concentration-QTc relationships demonstrate that erdafitinib does not prolong QTc interval and has no effects on cardiac repolarization or other ECG parameters. Clinical trial registration numbers NCT01703481, EudraCT: 2012-000697-34.

Immune network for viral hepatitis B: Topological representation

The liver is a well-known immunotolerogenic environment, which provides the adequate setting for liver infectious pathogens persistence such as the hepatitis B virus (HBV). Consequently, HBV infection can derive in the development of chronic disease in a proportion of the patients. If this situation persists in time, chronic hepatitis B (CHB) would end in cirrhosis, hepatocellular carcinoma and eventually, the death of the patient. It is thought that this immunotolerogenic environment is the result of complex interactions between different elements of the immune system and the viral biology. Therefore, the purpose of this work is to unravel the mechanisms implied in the development of CHB and to design a tool able to help in the study of adequate therapies. Firstly, a conceptual framework with the main components of the immune system and viral dynamics was constructed providing an overall insight on the pathways and interactions implied in this disease. Secondly, a review of the literature was performed in a modular fashion: (i) viral dynamics, (ii) innate immune response, (iii) humoral and (iv) cellular adaptive immune responses and (v) tolerogenic aspects. Finally, the information collected was integrated into a single topological representation that could serve as the plan for the systems pharmacology model architecture. This representation can be considered as the previous unavoidable step to the construction of a quantitative model that could assist in biomarker and target identification, drug design and development, dosing optimization and disease progression analysis.

Population Pharmacokinetics of AL-335 and Its Two Main Metabolites (ALS-022399, ALS-022227) in Monotherapy and in Combination with Odalasvir and/or Simeprevir

The aim of the current study was to characterize the time course of plasma concentrations of AL-335 and its main metabolites (ALS-022399 and ALS-022227) after oral administration in healthy and hepatitis C virus (HCV)-infected subjects, in monotherapy as well as in combination with simeprevir and/or odalasvir. AL-335, ALS-022399, and ALS-022227 plasma concentrations from subjects receiving 800 mg of AL-335 orally once daily (qd) as monotherapy or in combination were pooled and analyzed using a nonlinear mixed effect modeling approach. The typical values (between subject variability) of AL-335 and ALS-022399 apparent linear clearances were 3300 L/h (33.9%) and 1910 L/h (30.0%), respectively. ALS-022227 elimination was characterized as a nonlinear process, with typical values of V_{max,ALS-022227} and K_m,ALS-022227 estimated to be 84,799 ng/h (14.9%) and 450.2 ng/mL, respectively. AL-335 and ALS-022399 plasma concentrations were increased more than 2-fold in presence of simeprevir and/or odalasvir, while the effect on ALS-022227 plasma concentrations was limited. The effect of simeprevir and/or odalasvir might be explained by their capacity to inhibit P-glycoprotein. Internal evaluation confirmed that the population pharmacokinetic model developed was deemed appropriate to describe the time course of AL-335, ALS-022399, and ALS-022227 plasma concentrations and their associated variability in both healthy and HCV-infected subjects, as well as the interaction effect of simeprevir and/or odalasvir over AL-335 and its metabolites in healthy subjects. This model can be used as a starting point to evaluate drug-drug interaction processes in HCV-infected patients and support the development of a direct-acting antiviral (DAA) combination.

Characterizing the Pharmacokinetic Interaction Between Simeprevir and Odalasvir in Healthy Volunteers Using a Population Modeling Approach

The aim of this study was to characterize the pharmacokinetic drug-drug interaction (DDI) between simeprevir (NS3/4A protease inhibitor) and odalasvir (NS5A inhibitor) after oral administration to support the design and dose selection of clinical studies with this combination for the treatment of chronic hepatitis C infection (HCV). Simeprevir and odalasvir plasma concentrations from 30 healthy subjects receiving these drugs in monotherapy as well as in combination were pooled and analyzed using a population pharmacokinetic modeling approach. Previous pharmacokinetic models developed to characterize the pharmacokinetics for each drug were used as starting point. The dual effect of simeprevir and odalasvir on their pharmacokinetic parameters was explored. Simulations were performed to assess the impact of the DDI on exposure parameters. In presence of odalasvir, the relative bioavailability of simeprevir increased by 26% and the apparent clearance was reduced following competitive inhibition depending on odalasvir plasma concentrations, with an inhibitory constant (K_i) estimated to be 1610 ng/mL. The apparent odalasvir clearance was reduced by simeprevir plasma concentrations following an I_max model, characterized by a maximum inhibitory effect of 46.7% and an IC₅₀ of 257 ng/mL. Model-based simulations indicated that both C_max and AUC_24h increased for both drugs, when co-administered. The pharmacokinetic model adequately describes the time course of plasma concentrations and their variability when simeprevir and/or odalasvir were orally administered. This model can be used as a first step to predict the exposures of concomitant administration of simeprevir and odalasvir in HCV-infected subjects. Data from study AL355-602 (NCT02512562) were used for this analysis.

Population Pharmacokinetics and Pharmacodynamics of the Calcimimetic Etelcalcetide in Chronic Kidney Disease and Secondary Hyperparathyroidism Receiving Hemodialysis

Etelcalcetide is a novel calcimimetic in development for the treatment of secondary hyperparathyroidism (SHPT). A population pharmacokinetic/pharmacodynamic (PK/PD) model was developed relating etelcalcetide exposures to markers of efficacy (parathyroid hormone [PTH]) and safety (calcium) using data from three clinical studies. The semimechanistic model was developed that included allosteric activation pharmacology and understanding of calcium homeostasis. The temporal profiles for all biomarkers were well described by the model. The cooperativity constant was 4.94, confirming allosteric activation mechanism. Subjects with more severe disease (higher PTH baseline) were predicted to experience less pronounced reduction in PTH (percentage change from baseline), but more reduction in calcium (Ca; percentage change from baseline). There was no evidence that dose adjustment by any covariate was needed. Model‐based simulations provided quantitative support to several elements of dosing, such as starting dose, monitoring, and titration timing for registration trials.

Romiplostim dose response in patients with immune thrombocytopenia

A pharmacodynamic model was developed for platelet counts in 52 patients with immune thrombocytopenia (ITP) receiving subcutaneous romiplostim in 3 phase I/II studies (dose range, 0.2-10 µg/kg). The model consisted of a drug-sensitive progenitor cell compartment linked to a peripheral blood compartment through 4 transition compartments. The baseline platelet count, mean transit time, and kinetics of drug effect constant were 11.1 × 10(9)/L, 170 hours, and 0.6 day(-1), respectively. The ITP patients had a shorter platelet life span and lower progenitor cell production rates than healthy volunteers. Romiplostim response was described for 2 subpopulations. The romiplostim stimulatory effect in ITP patients was 351%/100 µg/wk and 12%/100 µg/wk in 68% and 32% of patients, respectively. Visual and numerical predictive checks suggested accurate prediction of platelet time course and durable response rate in ITP patients. Model-based simulations confirmed the effectiveness of dose reduction to prevent platelet counts >400 × 10(9)/L.

Mechanism-based Pharmacokinetic/Pharmacodynamic Meta-analysis of Trabectedin (ET-743, Yondelis) Induced Neutropenia

Myelosuppression was found to be one of the main toxicities of trabectedin (ET-743, Yondelis) during phase I/II studies. Our objective was to develop a pharmacokinetic-pharmacodynamic (PK/PD) model that describes the time course of the absolute neutrophil counts (ANCs) in cancer patients receiving trabectedin. Data from 699 patients who received intravenous trabectedin as monotherapy (dose range: 0.006-1.8 mg/m2) as a 1-, 3-, or 24-h infusion every 21 days; 1- or 3-h infusion on days 1, 8, and 15 every 28 days; or a 1-h infusion daily for 5 consecutive days every 21 days were used to develop (N=405; ANCs=7,291) and validate (N=294; ANCs=5,029) the model. The PK/PD model comprised a trabectedin-sensitive progenitor cell compartment, linked to the peripheral blood compartment, through three transition compartments representing the maturation chain in the bone marrow. To capture the rebound effect due to endogenous growth factors, the model included a feedback mechanism. The model estimated three system-related parameters: ANC at baseline (Circ0), mean transit time in bone marrow (MTT), and a feedback parameter (gamma). A first-order process quantified by the rate constant k(e0) described the trabectedin concentrations at the effect compartment (C(e)), which were assumed to reduce the proliferation rate and/or to increase the killing rate of the progenitor cells according to the function alphaC(e)beta. The model was qualified and simulations were undertaken to evaluate the neutropenia schedule dependency and the effects of selected covariates. NONMEM software was used to perform the modeling and simulation analyses. For a typical man of 70 kg, the mean values (between-subject variability; %) of the Circ0, MTT, gamma, k(e0), alpha, and beta were estimated to be 4.46 x 10(9)/l (37.9%), 4.0 days (37.5%), 0.218 (41.8%), 2.09 h(-1) (77.9%), 2.00 l/microg (85.1%), and 1.26, respectively. Although in women, k(e0) was reduced by 29% and a 25% increase in body weight resulted in a 12.6% reduction in the beta parameter, the clinical relevance of these effects is limited. The model evaluation procedure indicated accurate prediction of the observed incidence of neutropenia grades 3 and 4 across the dosing regimens evaluated. Simulations indicated that trabectedin dose and interdose interval, but not infusion duration, are the main determinants of the neutropenia severity. The model-predicted time course of the ANC and its variability confirmed that neutropenia is reversible, of short duration, and non-cumulative. The extent and time course of neutropenia following six different dosing regimens of trabectedin were well predicted by the semiphysiological PK/PD model.

Phase III double-blind placebo-controlled study of farnesyl transferase inhibitor R115777 in patients with refractory advanced colorectal cancer

PURPOSE

To determine whether R115777 improves survival in patients with refractory advanced colorectal cancer (CRC) in a multicenter, double-blind, prospective randomized study.

PATIENTS AND METHODS

Three hundred sixty-eight patients were randomly assigned to R115777 (300 mg twice daily) orally for 21 days every 28 days or placebo in a 2:1 ratio. All patients received best supportive care. The primary end point was overall survival; secondary end points were progression free survival, tumor response, toxicity, and quality of life.

RESULTS

The two treatment groups were well balanced for baseline demographics, including previous chemotherapy for advanced CRC. The median overall survival for R115777 was 174 days (95% CI, 157 to 198 days), and 185 days (95% CI, 158 to 238 days) for those patients receiving placebo (P =.376). One patient achieved a partial response in the R115777 arm. Stable disease (> 3 months) was observed in 24.3% patients in the R115777 group compared to 12.8% in the placebo arm. This did not translate into a statistically significant increase in progression-free survival. Overall, treatment was well tolerated. There was an increased incidence of reversible myelosuppression (neutropenia, thrombocytopenia), rash, and grade 1 to 2 diarrhea in the R115777 arm. There was no statistically significant difference in quality of life between arms.

CONCLUSION

Single agent R115777, given at this dose and schedule, has an acceptable toxicity profile, but does not improve overall survival compared to best supportive care alone in refractory advanced CRC.