Population Pharmacokinetics of Ibrutinib and Its Dihydrodiol Metabolite in Patients with Lymphoid Malignancies

Lymphocyte migration and retention properties affected by ibrutinib in chronic lymphocytic leukemia

The Bruton tyrosine kinase (BTK) inhibitor ibrutinib is widely used for treatment of patients with relapsed/refractory or treatment-naïve chronic lymphocytic leukemia (CLL). A prominent effect of ibrutinib is to disrupt the retention of CLL cells from supportive lymphoid tissues, by altering BTK-dependent adhesion and migration. To further explore the mechanism of action of ibrutinib and its potential impact on non-leukemic cells, we quantified multiple motility and adhesion parameters of human primary CLL cells and non-leukemic lymphoid cells. In vitro, ibrutinib affected CCL19-, CXCL12- and CXCL13-evoked migration behavior of CLL cells and non-neoplastic lymphocytes, by reducing both motility speed and directionality. De-phosphorylation of BTK induced by ibrutinib in CLL cells was associated with defective polarization over fibronectin and inability to assemble the immunological synapse upon B-cell receptor engagement. In patients' samples collected during a 6-month monitoring of therapy, chemokine-evoked migration was repressed in CLL cells and marginally reduced in T cells. This was accompanied by profound modulation of the expression of chemokine receptors and adhesion molecules. Remarkably, the relative expression of the receptors governing lymph node entry (CCR7) versus exit (S1PR1) stood out as a reliable predictive marker of the clinically relevant treatment-induced lymphocytosis. Together, our data reveal a multifaceted modulation of motility and adhesive properties of ibrutinib on both CLL leukemic cell and T-cell populations and point to intrinsic differences in CLL recirculation properties as an underlying cause for variability in treatment response.

Therapeutic Drug Monitoring of Kinase Inhibitors in Oncology

Although kinase inhibitors (KI) frequently portray large interpatient variability, a 'one size fits all' regimen is still often used. In the meantime, relationships between exposure-response and exposure-toxicity have been established for several KIs, so this regimen could lead to unnecessary toxicity and suboptimal efficacy. Dose adjustments based on measured systemic pharmacokinetic levels-i.e., therapeutic drug monitoring (TDM)-could therefore improve treatment efficacy and reduce the incidence of toxicities. Therefore, the aim of this comprehensive review is to give an overview of the available evidence for TDM for the 77 FDA/EMA kinase inhibitors currently approved (as of July 1st, 2023) used in hematology and oncology. We elaborate on exposure-response and exposure-toxicity relationships for these kinase inhibitors and provide practical recommendations for TDM and discuss corresponding pharmacokinetic targets when possible.

Co-delivery of ibrutinib and hydroxychloroquine by albumin nanoparticles for enhanced chemotherapy of glioma

Ibrutinib (IBR) is an oral covalent inhibitor of Bruton's tyrosine kinase (BTK) that has been approved for the treatment of hematological malignancies. It was reported that IBR exhibited great therapeutic potential for glioma. However, the poor water solubility and high hepatic first-pass effect restrict its anti-glioma application. Meanwhile, IBR induces cytoprotective autophagy through Akt/mTOR signaling pathway, thus leading to a compromised antitumor effect. Herein, we aimed to develop a human serum albumin (HSA) based co-delivery system (IBR&HCQ HSA NPs) encapsulating IBR and hydroxychloroquine (HCQ). The bioavailability of IBR was largely improved, and enhanced sensitivity of glioma to IBR was achieved due to inhibition effect of HCQ on IBR-induced pro-survival autophagy. The physicochemical properties of IBR&HCQ HSA NPs were characterized to optimize the formulation. Biodistribution investigation revealed that HSA NPs (20 mg/kg, i.v.) dramatically increased the accumulation of IBR in glioma, which was 5.59 times higher than that of free IBR (100 mg/kg, i.g.). CCK-8 and apoptosis assays demonstrated that IBR&HCQ HSA NPs showed maximal cytotoxicity to C6 cells. In vivo studies indicated that the survival time was significantly prolonged in IBR&HCQ HSA NPs treated mice compared to those treated with IBR HSA NPs. Taken together, the HSA-based drug delivery system of IBR and HCQ opens a new avenue for efficient treatment of glioma.

A High-Throughput Clinical Laboratory Methodology for the Therapeutic Monitoring of Ibrutinib and Dihydrodiol Ibrutinib

Ibrutinib (IBR) is an oral anticancer medication that inhibits Bruton tyrosine kinase irreversibly. Due to the high risk of adverse effects and its pharmacokinetic variability, the safe and effective use of IBR is expected to be facilitated by precision dosing. Delivering suitable clinical laboratory information on IBR is a prerequisite of constructing fit-for-purpose population and individual pharmacokinetic models. The validation of a dedicated high-throughput method using liquid chromatography-mass spectrometry is presented for the simultaneous analysis of IBR and its pharmacologically active metabolite dihydrodiol ibrutinib (DIB) in human plasma. The 6 h benchtop stability of IBR, DIB, and the active moiety (IBR+DIB) was assessed in whole blood and in plasma to identify any risk of degradation before samples reach the laboratory. In addition, four regression algorithms were tested to determine the optimal assay error equations of IBR, DIB, and the active moiety, which are essential for the correct estimation of the error of their future nonparametric pharmacokinetic models. The noncompartmental pharmacokinetic properties of IBR and the active moiety were evaluated in three patients diagnosed with chronic lymphocytic leukemia to provide a proof of concept. The presented methodology allows clinical laboratories to efficiently support pharmacokinetics-based precision pharmacotherapy with IBR.

Intentional Modulation of Ibrutinib Pharmacokinetics through CYP3A Inhibition

Ibrutinib (Imbruvica; PCI-32765) is an orally administered inhibitor of Bruton's tyrosine kinase that has transformed the treatment of B-cell malignancies. However, ibrutinib has very low oral bioavailability that contributes to significant variability in systemic exposure between patients, and this has the potential to affect both efficacy and toxicity. We hypothesized that the oral bioavailability of ibrutinib is limited by CYP3A isoform-mediated metabolism, and that this pathway can be inhibited to improve the pharmacokinetic properties of ibrutinib. Pharmacokinetic studies were performed in wild-type mice and mice genetically engineered to lack all CYP3A isoforms [CYP3A(-/-)] that received ibrutinib alone or in combination with CYP3A inhibitors cobicistat or ketoconazole. Computational modeling was performed to derive doses of ibrutinib that, when given after a CYP3A inhibitor, results in therapeutically-relevant drug levels. Deficiency of CYP3A in mice was associated with a ~10-fold increase in the area under the curve of ibrutinib. This result could be phenocopied by administration of cobicistat before ibrutinib in wild-type mice, but cobicistat did not influence levels of ibrutinib in CYP3A(-/-) mice. Population pharmacokinetic and prospectively validated physiologically-based pharmacokinetic models established preclinical and clinical doses of ibrutinib that could be given safely in combination with cobicistat without negatively affecting anti-leukemic properties. These findings signify a dominant role for CYP3A-mediated metabolism in the elimination of ibrutinib, and suggest a role for pharmacological inhibitors of this pathway to intentionally modulate the plasma levels and improve the therapeutic use of this clinically important agent.

Exposure-response analysis of acalabrutinib and its active metabolite, ACP-5862, in patients with B-cell malignancies

AIMS

Examine relationships between the systemic exposure of acalabrutinib, a highly selective, next-generation Bruton tyrosine kinase inhibitor, and its active metabolite (ACP-5862) vs. efficacy and safety responses in patients with B-cell malignancies who received acalabrutinib as monotherapy or in combination with obinutuzumab.

METHODS

For exposure-efficacy analyses, patients with untreated chronic lymphocytic leukaemia were assessed for best overall response, progression-free survival and tumour regression. For exposure-safety analyses, incidences of grade ≥2 adverse events (AEs), grade ≥3 AEs and grade ≥2 events of clinical interest were assessed in patients with B-cell malignancies. Acalabrutinib and ACP-5862 pharmacokinetic (PK) parameter estimates were obtained from population PK modelling. Exposure calculations were based on study dosing regimens. Total active moieties were calculated to account for contributions of ACP-5862 to overall efficacy/safety.

RESULTS

A total of 573 patients were included (exposure-efficacy analyses, n = 274; exposure-safety analyses, n = 573). Most patients (93%) received acalabrutinib 100 mg twice daily. Median total active area under the concentration-time curve (AUC_24h,ss ) and total active maximal concentration at steady-state (C_max,ss ) were similar for patients who received acalabrutinib as monotherapy or in combination with obinutuzumab, and for responders and nonresponders. No relationship was observed between AUC_24h,ss /C_max,ss and progression-free survival or tumour regression. Acalabrutinib AUC_24h,ss and C_max,ss were generally comparable across groups regardless of AE incidence.

CONCLUSION

No clinically meaningful correlations between acalabrutinib PK exposure and efficacy and safety outcomes were observed. These data support the fixed acalabrutinib dose of 100 mg twice daily in the treatment of patients with B-cell malignancies.

Easy and reliable maximum a posteriori Bayesian estimation of pharmacokinetic parameters with the open-source R package mapbayr

Pharmacokinetic (PK) parameter estimation is a critical and complex step in the model‐informed precision dosing (MIPD) approach. The mapbayr package was developed to perform maximum a posteriori Bayesian estimation (MAP‐BE) in R from any population PK model coded in mrgsolve. The performances of mapbayr were assessed using two approaches. First, “test” models with different features were coded, for example, first‐order and zero‐order absorption, lag time, time‐varying covariates, Michaelis–Menten elimination, combined and exponential residual error, parent drug and metabolite, and small or large inter‐individual variability (IIV). A total of 4000 PK profiles (combining single/multiple dosing and rich/sparse sampling) were simulated from each test model, and MAP‐BE of parameters was performed in both mapbayr and NONMEM. Second, a similar procedure was conducted with seven “real” previously published models to compare mapbayr and NONMEM on a PK outcome used in MIPD. For the test models, 98% of mapbayr estimations were identical to those given by NONMEM. Some discordances could be observed when dose‐related parameters were estimated or when models with large IIV were used. The exploration of objective function values suggested that mapbayr might outdo NONMEM in specific cases. For the real models, a concordance close to 100% on PK outcomes was observed. The mapbayr package provides a reliable solution to perform MAP‐BE of PK parameters in R. It also includes functions dedicated to data formatting and reporting and enables the creation of standalone Shiny web applications dedicated to MIPD, whatever the model or the clinical protocol and without additional software other than R.

Estimation of drug exposure by machine learning based on simulations from published pharmacokinetic models: The example of tacrolimus

We previously demonstrated that Machine learning (ML) algorithms can accurately estimate drug area under the curve (AUC) of tacrolimus or mycophenolate mofetil (MMF) based on limited information, as well as or even better than maximum a posteriori Bayesian estimation (MAP-BE). However, the major limitation in the development of such ML algorithms is the limited availability of large databases of concentration vs. time profiles for such drugs. The objectives of this study were: (i) to develop a Xgboost model to estimate tacrolimus inter-dose AUC based on concentration-time profiles obtained from a literature population pharmacokinetic (POPPK) model using Monte Carlo simulation; and (ii) to compare its performance with that of MAP-BE in external datasets of rich concentration-time profiles. The population parameters of a previously published PK model were used in the mrgsolve R package to simulate 9000 rich interdose tacrolimus profiles (one concentration simulated every 30 min) at steady-state. Data splitting was performed to obtain a training set (75%) and a test set (25%). Xgboost algorithms able to estimate tacrolimus AUC based on 2 or 3 concentrations were developed in the training set and the model with the lowest RMSE in a ten-fold cross-validation experiment was evaluated in the test set, as well as in 4 independent, rich PK datasets from transplant patients. ML algorithms based on 2 or 3 concentrations and a few covariates yielded excellent AUC estimation in the external validation datasets (relative bias < 5% and relative RMSE < 10%), comparable to those obtained with MAP-BE. In conclusion, Xgboost machine learning models trained on concentration-time profiles simulated using literature POPPK models allow accurate tacrolimus AUC estimation based on sparse concentration data. This study paves the way to the development of artificial intelligence at the service of precision therapeutic drug monitoring in different therapeutic areas.

Population Pharmacokinetics of Ibrutinib in Healthy Adults

BACKGROUND AND OBJECTIVES

Ibrutinib is an antineoplastic agent that reduces B-cell proliferation by inhibiting Bruton's tyrosine kinase. We describes population pharmacokinetics of ibrutinib in healthy adults, and explores potential patient characteristics associated with ibrutinib pharmacokinetics.

METHODS

A population pharmacokinetic modeling approach was applied to 39 healthy subjects. Modeling was performed using Monolix (v.2019R2). Serial blood samples to measure the plasma ibrutinib concentration were collected following the oral administration of 140 mg ibrutinib on two different occasions under fasting conditions. Demographic and clinical information were evaluated as possible predictors of ibrutinib pharmacokinetics during model development. Simulations (using mlxR: R package v.4.0.2) following the administration of therapeutic doses were performed to explore the clinical implications of identified covariates on ibrutinib steady-state concentrations.

RESULTS

A two-compartment model with zero order absorption best fit the data. Inter-individual and inter-occasion variability were quantified by the proposed model. We identified smoking status as a significant covariate associated with ibrutinib clearance. Smoking was found to increase ibrutinib clearance by approximately 60%, which resulted in a reduction in simulated steady-state concentrations by around 40%.

CONCLUSION

The model can be used to simulate clinical trials or various dosing scenarios. The proposed model can be used to optimize ibrutinib dosing based on the smoking status.

Limited Sampling Strategy for Determination of Ibrutinib Plasma Exposure: Joint Analyses with Metabolite Data

Therapeutic drug monitoring of ibrutinib is based on the area under the curve of concentration vs. time (AUC_IBRU) instead of trough concentration (C_min,ss) because of a limited accumulation in plasma. Our objective was to identify a limited sampling strategy (LSS) to estimate AUC_IBRU associated with Bayesian estimation. The actual AUC_IBRU of 85 patients was determined by the Bayesian analysis of the full pharmacokinetic profile of ibrutinib concentrations (pre-dose T0 and 0.5, 1, 2, 4 and 6 h post-dose) and experimental AUC_IBRU were derived considering combinations of one to four sampling times. The T0–1–2–4 design was the most accurate LSS (root-mean-square error RMSE = 11.0%), and three-point strategies removing the 1 h or 2 h points (RMSE = 22.7% and 14.5%, respectively) also showed good accuracy. The correlation between the actual AUC_IBRU and C_min,ss was poor (r² = 0.25). The joint analysis of dihydrodiol-ibrutinib metabolite concentrations did not improve the predictive performance of AUC_IBRU. These results were confirmed in a prospective validation cohort (n = 27 patients). At least three samples, within the pre-dose and 4 h post-dose period, are necessary to estimate ibrutinib exposure accurately.

Population PK-PD Modeling of Circulating Lymphocyte Dynamics in Chronic Lymphocytic Leukemia Patients Under Ibrutinib Treatment

Ibrutinib is indicated for the treatment of chronic lymphocytic leukemia (CLL). Absolute lymphocyte count (ALC) is a clinical criterion used for the monitoring of CLL. Ibrutinib has several effects on lymphocytes, and has highly variable pharmacokinetics (PK). The objective of this work was to build a PK-pharmacodynamic (PD) model describing ALC dynamics under ibrutinib treatment in patients with CLL. ALC observations before and after ibrutinib treatment initiation in patients with CLL were included in the analysis. A population PK-PD model was developed based on physio-pharmacological knowledge. Individual PK concentrations at each hospital visit were included in the model. The association between PD parameters and lymphocytosis, and between PD parameters and response to treatment were assessed. A total of 94 patients, 658 ALC and 1,501 PK observations were included in model development. The final PK-PD model accurately described ALC dynamics for different patient profiles. It consisted in two compartments (tissues and blood circulation) with ibrutinib plasmatic concentration inducing two drug effects: stimulation of lymphocyte redistribution and death. Patients with hyperlymphocytosis had significantly higher tissues to circulation baseline lymphocyte count ratio, and lower death effect. Patients who progressed under ibrutinib had significantly lower baseline lymphocyte counts in tissues (2-fold lower) and blood (3-fold lower). The first PK-PD model for ALC in patients with CLL under ibrutinib treatment was developed. This model suggests that estimated lymphocyte counts in tissues and blood could be used as an early predictor of response in patients with CLL.