- 1Experimental and Clinical Pharmacology, Centro di Riferimento Oncologico di Aviano (CRO), IRCCS, Aviano, Italy
- 2Department of Medical Oncology, Centro di Riferimento Oncologico di Aviano (CRO), IRCCS, Aviano, Italy
- 3Department of Medicine (DMED), University of Udine, Udine, Italy
- 4Doctoral School in Pharmacological Sciences, University of Padua, Padova, Italy
Palbociclib, an oral inhibitor of cyclin-dependent kinase 4 and 6, is approved for the treatment of metastatic breast cancer. This study investigated the influence of diverse clinical and biological factors—age, renal function, genetic variations, and concomitant medications (pharmacokinetic covariates)—on palbociclib pharmacokinetics. Employing a validated LC-MS/MS method, we analyzed the minimum plasma concentrations (Ctrough) of palbociclib in 68 women and determined the percentage deviations from the median Ctrough for each dosage group. Variations in a panel of absorption, distribution, metabolism, and excretion (ADME) genes were assessed using end-point allele-specific fluorescence detection and pyrosequencing. Two distinct patient cohorts were defined based on median values of age, creatinine, and eGFR, which exhibited statistically significant differences in percentage deviations (p = 0.0095, p = 0.0288, and p = 0.0005, respectively). Homozygous carriers of the PPARA variants displayed larger positive percentage deviations than the other group (p = 0.0292). Similarly, patients concurrently taking CYP3A and P-glycoprotein inhibitors alongside anticancer therapy exhibited significant variations (p = 0.0285 and p = 0.0334, respectively). Furthermore, exploring the drug–drug–gene interactions between inhibitors of CYP3A and P-glycoprotein with their respective genetic variants revealed two patient groups with statistically different percentage deviations (p = 0.0075, p = 0.0012, and p = 0.0191, respectively). These results could help address cases where pharmacokinetic covariates or subclinical conditions impair palbociclib adherence or response, aiming to offer tailored dosing strategies or monitoring for individual patients.
1 Introduction
Hormone receptor-positive (HR+) and human epidermal growth factor receptor 2-negative (HER2–) metastatic breast cancers account for a substantial proportion of all breast cancer cases. The evolution of targeted therapies has significantly altered the treatment paradigm for this subgroup, with a notable shift toward preferring endocrine therapy, integrated with cyclin-dependent kinase 4 and 6 (CDK4/6) inhibitors, over traditional chemotherapy (Bonotto et al., 2017; Giuliano et al., 2019).
Palbociclib, a selective CDK4/6 inhibitor, has emerged as a pivotal treatment in managing HR+/HER2– metastatic breast cancer. Approved by the Food and Drug Administration (FDA) and the European Medicines Agency (EMA) in 2015, palbociclib is administered at the starting dose of 125 mg once daily for 21 days, followed by 7 days off for a 28-day cycle (EMA, 2024; FDA, 2024b). Its pharmacokinetics reveal a low bioavailability (approximately 46%), a maximum concentration reached between 6- and 12-h post-ingestion, and a steady state achieved within 8 days of daily dosing. Its metabolism primarily occurs through cytochrome P450 (CYP) 3A isoform and sulfotransferase (SULT) 2A1. The major route of excretion is through feces, accounting for 74.1% of the dose, while approximately 17.5% of the dose is recovered in urine (Center for Drug Evaluation and Research, 2014). Additionally, palbociclib has been described, in in vivo and in vitro studies, as a substrate of the breast cancer resistance protein (BCRP) and the efflux transporter P-glycoprotein (P-gp) (Braal et al., 2021).
The pharmacokinetic profile of palbociclib, like other oral targeted therapies, is characterized by intra- and inter-individual variability, as indicated by a coefficient of variation of 41%–59% in trough concentrations (Groenland et al., 2020). This variability, influenced by covariates such as patient-specific characteristics (age, gender, race, weight, hormonal status, renal and hepatic function, and pharmacogenetics), diet, comorbidities, treatment adherence, and concomitant medications, underscores the need for a deeper understanding.
Previous research has shed light on several aspects of palbociclib pharmacokinetics; however, gaps remain in comprehensively understanding the influence of various covariates (Klümpen et al., 2011; Leenhardt et al., 2022; Roncato et al., 2022; Roncato et al., 2023; Peruzzi et al., 2023a).
This study aims to bridge these knowledge gaps by examining the clinical and biological factors that significantly impact the pharmacokinetics of palbociclib in patients with HR+/HER2- metastatic breast cancer. Additionally, given the higher likelihood of drug–drug interactions (DDIs) in cancer patients on multiple treatments, understanding the interactions resulting from concomitant medication is crucial. This includes pharmacokinetic interactions from the co-administration of inhibitors or inducers of enzymes and transporters involved in palbociclib’s pharmacokinetic profile, a critical aspect of all oral targeted therapies. Furthermore, this study also seeks to analyze the complex interplay between drug–gene interactions (DGIs) (Peruzzi et al., 2023b) and DDIs by assessing drug–drug–gene interactions (DDGIs), which allows us to evaluate the phenotype of absorption, distribution, metabolism, and excretion (ADME) enzymes and transporters based on both genetic variations and concomitant medications use.
Therefore, in a cohort of women with metastatic breast cancer treated with palbociclib, we determined the minimum palbociclib plasma concentration at the steady state (Ctrough) and investigated the influence of the aforementioned factors to highlight the clinical relevance of these factors in optimizing the therapeutic effect and improving patient safety by minimizing variability.
2 Materials and methods
2.1 Study design and patients
This study includes patients treated with palbociclib (days 1–21 of each 28-day cycle) affected by HR+/HER2- metastatic breast cancer and retrospectively selected from the cohort of the CRO–Aviano integrated pharmacological counseling program (Roncato et al., 2022) approved by the local ethics committee (Comitato Etico Unico Regionale del Friuli Venezia Giulia—CEUR) and conducted in accordance with the principles of the Declaration of Helsinki in its latest version (protocol ID: CRO 2022–14, approval date 12 April 2022; Parere-CEUR-2022-Os-65). Written informed consent was obtained from the patients.
The eligibility criteria were as follows: (i) patients with a diagnosis of HR+/HER2- metastatic breast cancer; (ii) patients who underwent treatment with palbociclib as either first- or second-line therapy, combined with aromatase inhibitors or fulvestrant; (iii) patients with available clinical data; (iv) patients with at least one blood sample collected at any palbociclib dose level (125 mg, 100 mg, or 75 mg) for evaluating the Ctrough.
2.2 Palbociclib plasma concentration
Palbociclib Ctrough was retrieved by plasma sample centrifuged at 16,200 g for 10 min at 4°C in whole blood EDTA tubes and analyzed with a validated LC-MS/MS method (Posocco et al., 2020; Posocco et al., 2024). To assess Ctrough, samples should be taken at a steady state and 24 h after the last dose. Samples collected within a 22–26 h window after the last intake of palbociclib were included. The median Ctrough value and the coefficient of variation (CV, in percentage) were calculated for each palbociclib dose (125 mg, 100 mg, and 75 mg). For each patient with more than one sample at the same dosage, the mean of the measured Ctrough values was considered.
To comprehensively analyze the Ctrough values associated with various doses of palbociclib, we calculated the percentage deviation of each Ctrough from the specific median Ctrough of each dosage level. This newly derived parameter, hereafter referred to as the “percentage deviations,” allowed for a detailed and comparative understanding of the pharmacokinetic behavior of palbociclib across its different therapeutic dosages. Positive percentage deviations represent Ctrough values higher than the dose-specific median, while negative values denote a Ctrough below the corresponding median. The magnitude of the absolute percentage deviations reflects the degree of deviation from the specific median Ctrough.
2.3 Pharmacokinetic covariates
Several factors were considered to analyze the palbociclib pharmacokinetic variability found in this selected case study. Such factors were classified as “pharmacokinetic covariates” and include patient clinical and pharmacological characteristics: age, hormonal status, treatment setting, hepatic and renal function, hematologic values, polymorphisms of ADME gene, and concomitant medications. All these factors were retrieved from the electronic medical record and/or directly from the patients at the time of sampling. Except for the pharmacogenetics, clinical characteristics were collected at baseline and at the time of sampling or, if not available, at the earliest time of sampling (within 2 weeks).
2.3.1 Clinical and biological characteristics
Only clinical–biological data associated to the sampling were considered for the evaluation of pharmacokinetic interactions. In instances where mean Ctrough values were considered, the average of the corresponding clinical and biological characteristics was also computed, provided that multiple data points were available. For hematologic data, white blood cell (WBC), absolute neutrophil count (ANC), and hemoglobin (Hgb) were evaluated. Aspartate and alanine aminotransferase (AST and ALT) and total bilirubin (BILT) were used to assess liver function, while creatinine and the calculated estimated glomerular filtration rate (eGFR) were used to assess renal function. The eGFR was calculated using the following formula derived from the Modification of Diet in Renal Disease (MDRD) study (Beck, 2023):
2.3.2 ADME gene polymorphisms
Polymorphism analyses of the ADME gene were performed by the SNPline PCR Genotyping System platform using Kompetitive allele–specific assays (LGC Genomics, Hoddesdon, UK) following DNA extraction from blood samples. The DNA extraction was done using the BioRobot EZ1 system workstation (Qiagen N.V., Germany) and the EZ1 DNA Blood 200 µL Kit. The analyzed polymorphisms were CYP3A (e.g., CYP3A4*22, rs35599367; CYP3A5*3, rs776746; CYP3A5*6, rs10264272; CYP3A5*7, rs41303343), nuclear receptor influencing CYP450 (e.g., POR*28, rs1057868; PPARA rs4253728; PPARA rs4823613; PXR rs2472677; PXR rs3814055; VDR rs11568820; VDR rs1544410; VDR rs4516035), and transporters (e.g., ABCB1 c.1236C>T, rs1128503; ABCB1 c.3435C>T, rs1045642; ABCG2 c.421C>A, rs2231142). In addition, the triallelic polymorphisms of ABCB1 2677G>T/A (rs2032582) were analyzed using pyrosequencing technology from PyroMark Q48 (Qiagen, Hilden, Germany) to achieve triallelic discrimination. Positive and negative control samples were included in each analysis.
2.3.2.1 Metabolizer status attribution
Each patient was given a CYP3A metabolizer status. Poor metabolizers (PM) were carriers of the CYP3A4*1/*22 or *22/*22 and CYP3A5*3/*3 diplotypes. Intermediate metabolizers (IM) were carriers of the CYP3A4*1/*1 and CYP3A5*3/*3 diplotypes as well as the CYP3A4*1/*22 or *22/*22 and CYP3A5*1/*3 or *1/*1 diplotypes. Lastly, normal metabolizers (NM) were carriers CYP3A4*1/*1 and CYP3A5*1/*3 or *1/*1 diplotypes (Lloberas et al., 2017; Sanchez Spitman et al., 2017).
2.3.3 Concomitant medication
All recorded concomitant medications were classified according to the therapeutic class and the anatomical–therapeutic–chemical (ATC) classification (last update 23 January 2023). To analyze the influence of concomitant medications on the pharmacokinetic variability of palbociclib, we first examined the influence of the number of concomitant medications each patient was taking concomitantly with oncology therapy and at the time of sampling. Patients with zero or one co-medication were compared with patients with 2–4 concomitant medications and with patients who had five or more concomitant medications. Next, the influence of co-administering at least one drug that acts as an inhibitor or inducer of CYP3A and P-gp, and BCRP was analyzed. Such inhibitors or inducers were classified based on four different sources: US Food and Drug Administration Pharmacological Review and Drug Interactions Flockhart Table and online sources “UptoDate” and “Cancer Liverpool Interactions” (Clinical Pharmacology School of Medicine, 2024; FDA, 2024a; UpToDate Free, 2023). Patients taking any type of inhibitor or inducer (weak, moderate, or strong) were compared with those not taking any them. Finally, the impact of the most common concomitant drugs recorded in the selected case study was also investigated.
2.3.4 Drug–drug–gene interactions
To evaluate the impact of DDGI on palbociclib Ctrough, patients with genetic variations that reduce the functionality of palbociclib ADME enzymes and transporters and/or those taking inhibitors of these enzymes and transporters were grouped together and compared to others.
2.4 Statistical analysis
Quantitative variables were reported as the number of observations (n), frequency (%), median, and interquartile range (IQR), while qualitative variables were described as the number of observations (n) and frequency (%). Then, the impact of each pharmacokinetic covariate was evaluated by a non-parametric Mann–Whitney or Kruskal–Wallis test. The plotted graphs report the median and the IQR for each compared group. Regarding ADME gene polymorphisms, the best-fitting genetic model (dominant or recessive) was used to perform the association analyses. All statistical analyses and graphs were performed using GraphPad Prism 9.2.0 version, and the significance level was set at 0.05.
3 Results
3.1 Patients and palbociclib plasma concentration
In this study, a total of 68 female patients diagnosed with HR+/HER2- metastatic breast cancer were enrolled, receiving treatment with palbociclib at the National Cancer Institute, IRCCS, Aviano. The baseline clinical characteristics of these patients are delineated in Table 1. Predominantly, the cohort comprised individuals undergoing first-line treatment with letrozole as the primary anti-hormonal therapy. Among the 68 patients, 45 (66%) contributed a single sample for the assessment of palbociclib Ctrough, while 23 (34%) provided multiple samples (two or three). Figure 1 shows the collected samples categorized based on the dosage. The median Ctrough values observed were 68 ng/mL (IQR 57–75; CV 30%) for the 125 mg dosage, 57 ng/mL (IQR: 47–76; CV: 40%) for the 100 mg dosage, and 33 ng/mL (IQR: 28–46; CV: 40%) for the 75 mg dosage.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1420174/fphar-15-1420174-HTML-r1/image_m/fphar-15-1420174-g001.jpg)
Figure 1. Distribution of samples based on the different dosages. The horizontal bar represents the median with the IQR, while the dots represent values of minimum plasma concentrations measured.
3.2 Clinical and biological data and palbociclib plasma concentration
Table 2 comprises the available clinical and biological parameters recorded for each patient’s sample. Figure 2 shows the results of the Kruskal–Wallis tests, examining the association between the percentage deviation and selected clinical–biological factors. Patients aged 65 years and older demonstrated a significantly more pronounced positive percentage deviation than those under 65 years, with median values for the older group being 13 percentage points higher (p = 0.0095). Although a similar trend was observed in relation to the menopausal status, it did not reach statistical significance (p = 0.0581). Patients with creatinine levels above 0.75 mg/dL had 20 percentage-point higher deviations than those with lower levels (p = 0.0288). Similarly, patients with eGFR below 82.02 mL/min/1.73 m2 had 33 percentage-point higher deviations than those with higher eGFR (p = 0.0005). Conversely, no significant differences were observed for other factors such as ANC, WBC, hemoglobin, ALT, AST, and bilirubin.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1420174/fphar-15-1420174-HTML-r1/image_m/fphar-15-1420174-g002.jpg)
Figure 2. Associations between the clinical–biological factors recorded and the percentage deviations. The box plot depicts the median (horizontal bar) and the IQR. The dots represent the percentage deviations. Abbreviations: ns, not significant; eGFR, estimated glomerular filtration rate.
3.3 Pharmacogenetic data and palbociclib plasma concentration
Sixteen single-nucleotide polymorphisms (SNPs) were tested, and their distributions across genotypes are reported in Supplementary Table S1, while Figure 3 illustrates the associations investigated. Although complete SNP data were not obtainable for all patients, those analyzed did not show significant deviation from the Hardy–Weinberg equilibrium. Notably, significant associations were observed only for two polymorphisms in the ADME gene, PPARA (rs4253728 and rs4823613), which are in linkage disequilibrium. Patients with homozygous variant genotypes for both polymorphisms exhibited significantly greater positive percentage deviations, with a median deviation 34 percentage points higher than those with wild-type or heterozygous genotypes (p = 0.0292).
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1420174/fphar-15-1420174-HTML-r1/image_m/fphar-15-1420174-g003.jpg)
Figure 3. Associations between pharmacogenetic data detected and the percentage deviations calculated. Regarding CYP3A, the status of the three metabolizers was evaluated. For ABCB1, the 1236T-3435T-2677T/A haplotype was considered, and the homozygous carriers were compared with all others. Finally, homozygous carriers of both variant polymorphisms of PPARA were compared with others. The box plot depicts the median (horizontal bar) and the IQR. The dots represent the percentage deviations. Abbreviations: ns, not significant;, PM, poor metabolizer; IM, intermediate metabolizer; NM, normal metabolizer.
3.4 Concomitant medication and palbociclib plasma concentration
Regarding concomitant medications, all patients were on at least one additional drug, with an average of four per patient. Ten patients (14.7%) were on a single concomitant medication, while 37 (54.4%) were taking between two and four medications, and 19 (27.9%) were on five or more. Figure 4 presents these data.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1420174/fphar-15-1420174-HTML-r1/image_m/fphar-15-1420174-g004.jpg)
Figure 4. Distribution of percentage deviation values based on the number of concomitant medications. The box plot depicts the median (horizontal bar) and the IQR. The dots represent the percentage deviations.
Our analysis specifically focused on inhibitors and inducers of CYP3A, P-gp, and BCRP. Only a limited number of weak and moderate inhibitors of CYP3A and P-gp were recorded, with just one patient taking a P-gp inducer. No BCRP inducers or inhibitors were found. Associations identified with these elements are shown in Figure 5. Patients taking weak or moderate inhibitors were grouped in the “yes” column. Notably, a statistically significant difference was observed with CYP3A and P-gp inhibitors, with median values for the group taking at least one CYP3A inhibitor being 21 percentage points higher and those for the group taking at least one P-gp inhibitor being 19 percentage points higher than those not taking these inhibitors (p = 0.0285 and p = 0.0334, respectively).
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1420174/fphar-15-1420174-HTML-r1/image_m/fphar-15-1420174-g005.jpg)
Figure 5. Associations between inhibitors detected and the percentage deviations. The box plot depicts the median (horizontal bar) and the IQR. The dots represent the percentage deviations.
Additionally, Figure 6 illustrates the associations with the most frequently recorded concomitant medications in our study: antihypertensives, anxiolytics/hypnotics (especially benzodiazepines), and proton pump inhibitors (PPIs). No statistically significant difference was observed in palbociclib exposure among patients taking these types of concomitant medications.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1420174/fphar-15-1420174-HTML-r1/image_m/fphar-15-1420174-g006.jpg)
Figure 6. Associations between concomitant medications detected and the percentage deviation. The box plot depicts the median (horizontal bar) and the IQR. The dots represent the percentage deviations. Abbreviations: ns, not significant; PPI, proton pump inhibitor.
3.5 Drug–drug–gene interactions and palbociclib plasma concentration
The DDGI analyses are shown in Figure 7. The first analysis, “CYP3A DDGI,” compares patients exhibiting a PM CYP3A phenotype and/or those taking CYP3A inhibitors as concomitant medications (grouped in the “yes” column) with other patients. In the second analysis, the “yes” column includes patients with a PM CYP3A phenotype, carriers of homozygous variants in both PPARA polymorphisms, and/or those using CYP3A inhibitors. Finally, a similar analysis was conducted for ABCB1, comparing patients carrying mutated ABCB1 haplotypes and/or those taking P-gp inhibitors as concomitant medications (grouped in the “yes” column) with other patients.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1420174/fphar-15-1420174-HTML-r1/image_m/fphar-15-1420174-g007.jpg)
Figure 7. Associations between DDGI and the percentage deviations. The box plot depicts the median (horizontal bar) and the IQR. The dots represent the percentage deviations. Abbreviation: DDGI, drug–drug–gene interactions.
Specifically, the median values for the group with CYP3A DDGI were 20 percentage points higher than those without. For the group with CYP3A-PPARA DDGI, the median values were 24.5 percentage points higher, and for the group with ABCB1 DDGI, the median values were 18 percentage points higher than those without these DDGIs. All three associations were statistically significant (p = 0.0075, p = 0.0012, and p = 0.0191, respectively).
4 Discussion
The findings of this study offer significant insights into the pharmacokinetic variability of palbociclib in patients with HR+/HER2- metastatic breast cancer, particularly in light of its complex metabolism and excretion pathways. The study’s detailed examination of palbociclib pharmacokinetic variability highlights the significance of patient-specific factors like age, renal function, genetic polymorphisms, DDIs, and DDGIs in influencing drug exposure, offering findings that might be crucial for a deeper understanding of pharmacokinetic variability.
The median Ctrough, calculated within our study group treated with the standard 125 mg dose of palbociclib aligns with the mean value of this parameter, as reported in the FDA’s Clinical Pharmacology Review. The reported palbociclib Ctrough value stands at 61 ng/mL with a CV of 42%. This characterization was conducted in both healthy subjects and patients with solid tumors, including those with advanced breast cancer (Center for Drug Evaluation and Research, 2014). Additionally, we provided median Ctrough values for two-step dose reductions, specifically 100 mg and 75 mg, showcasing decreases of 16% and 51%, respectively, from the median Ctrough calculated at the standard dosage.
Regarding the covariates analyzed, our analysis revealed a notable age-related difference in palbociclib exposure. Patients older than 65 years exhibited a higher median palbociclib exposure (expressed as a percentage deviation from the median Ctrough). This finding aligns with the understanding that physiological changes associated with aging can affect drug metabolism and distribution, suggesting a potential need for age-adjusted dosing strategies in older patients.
The study also identified significant associations between palbociclib exposure and renal function markers such as creatinine levels and eGFR. Given that eGFR is calculated using age and creatinine levels, the significant association of eGFR with percentage deviation may be influenced by these underlying factors. This is in line with findings by Leenhardt F. et al., who reported that palbociclib Ctrough was higher than the median values in older patients and those with reduced kidney function (Leenhardt et al., 2022). These findings suggest that even subclinical alterations in renal function, potentially associated with age, may still be important factors to consider in the dose adjustment and monitoring of palbociclib-treated patients. This is a noteworthy consideration, especially given the drug’s excretion profile, which indicates that palbociclib is predominantly excreted in the feces.
Our genetic analysis revealed a noteworthy association between specific PPARA gene polymorphisms and palbociclib exposure, underscoring the potential role of genetic factors in individual variability in drug response. Carriers of homozygous variant alleles for both PPARA gene polymorphisms exhibited a higher median palbociclib exposure (expressed as a percentage deviation from the median Ctrough) than patients who were either heterozygous or carriers of the wild-type alleles for both variants. Therefore, this insight brings to light the key role of nuclear receptors, which are known to be involved in the regulation and transcription of the cytochrome gene and to explain a percentage of the CYP3A4 activity variation.
Our study further explored the impact of concomitant medications on palbociclib exposure. While the drugs assessed in our study were not potent cytochrome inhibitors, the significant increase in exposure observed in patients receiving CYP3A and P-gp inhibitors underscores the necessity of incorporating DDIs into patient evaluations. This is essential to accurately ascertain the metabolizer status for ADME enzymes in conjunction with DGIs. This consideration becomes increasingly relevant due to the prevalent use of CYP3A substrates among cancer patients, both as part of concomitant medications and within oral anticancer agents. Leenhardt F. et al. also provided evidence of a higher palbociclib Ctrough in patients taking at least one CYP3A4 or P-gp inhibitor (Leenhardt et al., 2022). Furthermore, our data reveal that when considering various factors in combination, particularly the use of CYP3A inhibitors and a reduced-function phenotype (PM) of the CYP3A gene, there is a notable increase in palbociclib exposure (expressed as a percentage deviation from the median Ctrough) compared to other groups. This was also distinctly observed in individuals with homozygous variants of the PPARA polymorphisms.
Figure 8 summarizes our findings on the crucial roles of P-gp and CYP3A in the ADME processes involving palbociclib. Reduced function of P-gp and CYP3A, as assessed by evaluating DDGIs, may contribute to increased concentrations of palbociclib in the systemic circulation.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1420174/fphar-15-1420174-HTML-r1/image_m/fphar-15-1420174-g008.jpg)
Figure 8. Illustration of the roles of P-gp and CYP3A in the palbociclib ADME processes. Reduced function of P-gp and CYP3A, as assessed through DDGI, can lead to increased concentrations of palbociclib in the bloodstream. Created with BioRender.com.
Our study faced limitations due to its relatively small sample size and descriptive nature, which restricted our capacity to fully elucidate the clinical significance of the observed palbociclib trough concentration variations. A more comprehensive analysis, potentially incorporating clinical outcomes, would be necessary to fully understand the implications of these pharmacokinetic variations. From a clinical perspective, higher palbociclib Ctrough values, as indicated by the higher percentage deviations in our results, are expected to increase toxicity, particularly neutropenia. Although increased exposure seems to be associated with a greater reduction in ANC, the data on efficacy, such as progression-free survival, remain inconclusive (Flaherty et al., 2012; Center for Drug Evaluation and Research, 2014). Therefore, future research should prioritize not only validating our findings in larger cohorts but also exploring whether the patient-specific factors we investigated are implicated in the mechanisms underlying these exposure–toxicity associations. This will help to refine palbociclib dosing guidelines to manage toxicity more effectively.
Furthermore, distinct from numerous studies that rely on PK/PD models and population pharmacokinetics from healthy volunteers, our plasma concentration data originate from a real-world cohort of women diagnosed with metastatic breast cancer and undergoing treatment with palbociclib.
5 Conclusion
In conclusion, our study offers an in-depth perspective of the pharmacokinetic variability of palbociclib, highlighting the importance of considering a spectrum of factors such as age, renal function, genetic polymorphisms, and concomitant medication in tailoring treatment for patients with HR+/HER2- metastatic breast cancer. It notably highlights the impact of subclinical factors, including variations in creatinine levels and eGFR, or the involvement of moderate to weak inhibitors, on the patient-specific deviations from the median Ctrough values. This research contributes to the development of more personalized treatment strategies, with the goal of optimizing therapeutic outcomes for each patient undergoing palbociclib treatment for metastatic breast cancer.
Data availability statement
The original contributions presented in the study are included in the article and Supplementary Material, further inquiries can be directed to the corresponding authors.
Ethics statement
The studies involving humans were approved by the Comitato Etico Unico Regionale del Friuli Venezia Giulia—CEUR and conducted in accordance with the principles of the Declaration of Helsinki. The studies were conducted in accordance with the local legislation and institutional requirements. The participants provided their written informed consent to participate in this study.
Author contributions
EP: formal analysis, writing–original draft, writing–review and editing. BP: formal analysis, writing–original draft, writing–review and editing. LG: resources, writing–review and editing. MN: formal analysis, writing–review and editing. MO: formal analysis, writing–review and editing. SG: formal analysis, writing–review and editing. ED: resources, writing–review and editing. FP: writing–original draft, writing–review and editing. EC: supervision, writing–review and editing. GT: supervision, writing–review and editing. RR: supervision, writing–original draft, writing–review and editing.
Funding
The author(s) declare that financial support was received for the research, authorship, and/or publication of this article. This work was supported by the Italian Ministry of Health (Ricerca Corrente); 5 x 1000 “Seed Grant–Ricerca Sanitaria” J35F21000650007.
Conflict of interest
LG reports consulting or advisory roles with AstraZeneca, Daiichi Sankyo, Eli Lilly, GlaxoSmithKline, Incyte, Novartis, Pfizer, Menarini Stemmline, and AbbVie; and research funding for Menarini Silicon Biosystems. FP reports consulting or advisory roles with Roche, MSD, AstraZeneca, Novartis, Eli Lilly, Pfizer, Pierre Fabre, Daiichi Sankyo, Eisai, and Amgen, in addition to research funding from Eisai, AstraZeneca, and Roche.
The remaining authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
The author(s) declared that they were an editorial board member of Frontiers, at the time of submission. This had no impact on the peer review process and the final decision.
Publisher’s note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors, and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Supplementary material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fphar.2024.1420174/full#supplementary-material
References
Bonotto, M., Gerratana, L., Di Maio, M., De Angelis, C., Cinausero, M., Moroso, S., et al. (2017). Chemotherapy versus endocrine therapy as first-line treatment in patients with luminal-like HER2-negative metastatic breast cancer: a propensity score analysis. Breast 31, 114–120. doi:10.1016/j.breast.2016.10.021
Braal, C. L., Jongbloed, E. M., Wilting, S. M., Mathijssen, R. H. J., Koolen, S. L. W., and Jager, A. (2021). Inhibiting CDK4/6 in breast cancer with palbociclib, ribociclib, and Abemaciclib: Similarities and differences. Drugs 81, 317–331. doi:10.1007/s40265-020-01461-2
Center for Drug Evaluation and Research (2014). Clinical pharmacology and biopharmaceutics review. Available at: https://www.accessdata.fda.gov/drugsatfda_docs/nda/2015/207103orig1s000clinpharmr.pdf (Accessed February 10, 2024).
Clinical Pharmacology School of Medicine (2024). Cytochrome P450 Drug Interaction Table TM. Available at: https://drug-interactions.medicine.iu.edu/MainTable.aspx (Accessed February 10, 2024).
EMA (2024). Ibrance epar product information. Available at: https://www.ema.europa.eu/en/documents/product-information/ibrance-epar-product-information_en.pdf (Accessed February 10, 2024).
FDA (2024a). Drugs. Available at: https://www.fda.gov/drugs (Accessed February 13, 2024).
FDA (2024b). Palbociclib FDA prescribing information. Available at: https://www.accessdata.fda.gov/drugsatfda_docs/label/2016/207103s002lbl.pdf (Accessed February 10, 2024).
Flaherty, K. T., Lorusso, P. M., Demichele, A., Abramson, V. G., Courtney, R., Randolph, S. S., et al. (2012). Phase I, dose-escalation trial of the oral cyclin-dependent kinase 4/6 inhibitor PD 0332991, administered using a 21-day schedule in patients with advanced cancer. Clin. Cancer Res. 18, 568–576. doi:10.1158/1078-0432.CCR-11-0509
Giuliano, M., Schettini, F., Rognoni, C., Milani, M., Jerusalem, G., Bachelot, T., et al. (2019). Endocrine treatment versus chemotherapy in postmenopausal women with hormone receptor-positive, HER2-negative, metastatic breast cancer: a systematic review and network meta-analysis. Lancet Oncol. 20, 1360–1369. doi:10.1016/S1470-2045(19)30420-6
Groenland, S. L., Martínez-Chávez, A., van Dongen, M. G. J., Beijnen, J. H., Schinkel, A. H., Huitema, A. D. R., et al. (2020). Clinical pharmacokinetics and Pharmacodynamics of the cyclin-dependent kinase 4 and 6 inhibitors palbociclib, ribociclib, and Abemaciclib. Clin. Pharmacokinet. 59, 1501–1520. doi:10.1007/s40262-020-00930-x
Klümpen, H.-J., Samer, C. F., Mathijssen, R. H. J., Schellens, J. H. M., and Gurney, H. (2011). Moving towards dose individualization of tyrosine kinase inhibitors. Cancer Treat. Rev. 37, 251–260. doi:10.1016/j.ctrv.2010.08.006
Leenhardt, F., Fiteni, F., Gauthier, L., Alexandre, M., Guiu, S., Firmin, N., et al. (2022). Pharmacokinetic variability Drives palbociclib-Induced Neutropenia in metastatic breast cancer patients: drug-drug interactions are the Usual Suspects. Pharmaceutics 14, 841. doi:10.3390/pharmaceutics14040841
Lloberas, N., Elens, L., Llaudó, I., Padullés, A., van Gelder, T., Hesselink, D. A., et al. (2017). The combination of CYP3A4*22 and CYP3A5*3 single-nucleotide polymorphisms determines tacrolimus dose requirement after kidney transplantation. Pharmacogenet Genomics 27, 313–322. doi:10.1097/FPC.0000000000000296
Peruzzi, E., Gerratana, L., Montico, M., Posocco, B., Corsetti, S., Bartoletti, M., et al. (2023a). Association of ADME gene polymorphisms on toxicity to CDK4/6 inhibitors in patients with HR+ HER2-metastatic breast cancer. Biomed. Pharmacother. 167, 115479. doi:10.1016/j.biopha.2023.115479
Peruzzi, E., Roncato, R., De Mattia, E., Bignucolo, A., Swen, J. J., Guchelaar, H.-J., et al. (2023b). Implementation of pre-emptive testing of a pharmacogenomic panel in clinical practice: where do we stand? Br. J. Clin. Pharmacol. doi:10.1111/bcp.15956
Posocco, B., Buzzo, M., Poetto, A. S., Orleni, M., Gagno, S., Zanchetta, M., et al. (2020). Simultaneous quantification of palbociclib, ribociclib and letrozole in human plasma by a new LC-MS/MS method for clinical application. PLOS ONE 15, e0228822. doi:10.1371/journal.pone.0228822
Posocco, B., Zanchetta, M., Orleni, M., Gagno, S., Montico, M., Peruzzi, E., et al. (2024). Therapeutic monitoring of palbociclib, ribociclib, Abemaciclib, M2, M20, and letrozole in human plasma: a Novel LC-MS/MS method. Ther. Drug Monit. 46, 485–493. doi:10.1097/FTD.0000000000001174
Roncato, R., Gerratana, L., Palmero, L., Gagno, S., Poetto, A. S., Peruzzi, E., et al. (2022). An integrated pharmacological counselling Approach to Guide decision-Making in the treatment with CDK4/6 inhibitors for metastatic breast cancer. Front. Pharmacol. 13, 897951. doi:10.3389/fphar.2022.897951
Roncato, R., Peruzzi, E., Gerratana, L., Posocco, B., Nuzzo, S., Montico, M., et al. (2023). Clinical impact of body mass index on palbociclib treatment outcomes and effect on exposure. Biomed. Pharmacother. 164, 114906. doi:10.1016/j.biopha.2023.114906
Sanchez Spitman, A. B., Moes, D. J. R., Gelderblom, H., Dezentje, V. O., Swen, J. J., and Guchelaar, H. J. (2017). Effect of CYP3A4*22, CYP3A5*3, and CYP3A combined genotypes on tamoxifen metabolism. Eur. J. Clin. Pharmacol. 73, 1589–1598. doi:10.1007/s00228-017-2323-2
UpToDate Free (2023). Drugs Interactions 2023. Available at: https://pro.uptodatefree.ir/Drugs-interactions/(Accessed February 10, 2024).
Keywords: palbociclib, metastatic breast cancer, minimum plasma concentration, pharmacokinetic variability, pharmacokinetic covariate, drug-drug interactions, drug-gene interaction, drug-drug-gene interactions
Citation: Peruzzi E, Posocco B, Gerratana L, Nuti M, Orleni M, Gagno S, De Mattia E, Puglisi F, Cecchin E, Toffoli G and Roncato R (2024) Exploring pharmacokinetic variability of palbociclib in HR+/HER2- metastatic breast cancer: a focus on age, renal function, and drug–gene interactions. Front. Pharmacol. 15:1420174. doi: 10.3389/fphar.2024.1420174
Received: 19 April 2024; Accepted: 15 July 2024;
Published: 06 September 2024.
Edited by:
Yurong Lai, Gilead, United StatesReviewed by:
James M Rae, University of Michigan, United StatesTanja Blagus, University of Ljubljana, Slovenia
Copyright © 2024 Peruzzi, Posocco, Gerratana, Nuti, Orleni, Gagno, De Mattia, Puglisi, Cecchin, Toffoli and Roncato. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Bianca Posocco, bposocco@cro.it; Giuseppe Toffoli, gtoffoli@cro.it