Abstract
Background: Bruton’s tyrosine kinase (BTK) inhibitors are commonly used in the targeted therapy of B-cell malignancies. It is reported that myelosuppression and fungal infections might occur during antitumor therapy of BTK inhibitors, therefore a combination therapy with triazole antifungals is usually required.
Objective: To evaluate the influence of different triazoles (voriconazole, fluconazole, itraconazole) on the pharmacokinetics of BTK inhibitors (zanubrutinib, acalabrutinib) and to quantify the drug-drug interactions (DDIs) between them.
Methods: The physiologically-based pharmacokinetic (PBPK) models were developed based on pharmacokinetic parameters and physicochemical data using Simcyp® software. These models were validated using clinically observed plasma concentrations data which based on existing published studies. The successfully validated PBPK models were used to evaluate and predict potential DDIs between BTK inhibitors and different triazoles. BTK inhibitors and triazole antifungal agents were simulated by oral administration.
Results: Simulated plasma concentration-time profiles of the zanubrutinib, acalabrutinib, voriconazole, fluconazole, and itraconazole are consistent with the clinically observed profiles which based on existing published studies, respectively. The exposures of BTK inhibitors increase by varying degrees when co-administered with different triazole antifungals. At multiple doses regimen, voriconazole, fluconazole and itraconazole may increase the area under plasma concentration-time curve (AUC) of zanubrutinib by 127%, 81%, and 48%, respectively, and may increase the AUC of acalabrutinib by 326%, 119%, and 264%, respectively.
Conclusion: The PBPK models sufficiently characterized the pharmacokinetics of BTK inhibitors and triazole antifungals, and were used to predict untested clinical scenarios. Voriconazole exhibited the greatest influence on the exposures of BTK inhibitors. The dosage of zanubrutinib or acalabrutinib need to be reduced when co-administered with moderate CYP3A inhibitors.
Introduction
Hematologic malignancies are severe hematopoietic diseases which often accompanied by invasive fungal infections (IFIs) (Neofytos et al., 2013; Zeng et al., 2021). This not only due to the malignancies, but also due to the antitumor treatment, such as cytotoxic chemotherapy (Hamalainen et al., 2008), targeted immunotherapies (Lanini et al., 2011), long-term intravenous catheters (Heidenreich et al., 2022), and chemo-radiotherapy (Martino et al., 1997). Hematological malignancies accompanied by IFIs may increase the tumor recurrence and mortality of the patients (Lewis et al., 2013), so it is necessary to start the antifungal treatment as soon as possible.
According to the clinical practice guidelines of Infectious Diseases Society of America (IDSA), triazole antifungal agents are recommended for the prevention and treatment of IFIs, such as voriconazole, fluconazole and itraconazole (Perfect et al., 2010; Pappas et al., 2016; Patterson et al., 2016). Triazole antifungals are mainly metabolized by cytochrome P450 enzymes (CYP450), including CYP2C19, CYP2C9, and CYP3A4 (Bellmann and Smuszkiewicz, 2017), meanwhile they strongly inhibit CYP3A enzymes (Bellmann and Smuszkiewicz, 2017; Han et al., 2021; Ou et al., 2021). In fact, it is difficult to avoid the long-term consolidation therapy for antitumor and antifungal. In this process, the drug-drug interactions (DDIs) may increase the risk of drug toxicities, sub-optimal therapy, and drug resistance.
Over the past decade, with the rapid development of targeted therapy, many tyrosine kinase inhibitors (TKIs) have been approved for the treatment of hematological malignancies. Bruton’s tyrosine kinase (BTK) inhibitors such as zanubrutinib and acalabrutinib are increasingly replacing chemotherapy-based regimens, especially for patients with mantle cell lymphoma (MCL), chronic lymphocytic leukemia (CLL) (Burger, 2019) and small lymphocytic lymphoma (SLL) (Abbas and Wierda, 2021; Tam et al., 2021). According to the pharmacokinetic studies, zanubrutinib and acalabrutinib are mainly metabolized by CYP3A in the liver. When BTK inhibitors are co-administered with triazoles, the exposures of BTK inhibitors tend to increase, which may result in serious adverse effects, such as hematological toxicity, dermatological toxicities and diarrhea (Lipsky and Lamanna, 2020). To the best of our knowledge, at present, only a few reports have suggested the empirical reduction of BTK inhibitors in combination with CYP inhibitors (Hardy-Abeloos et al., 2020; Bruggemann et al., 2022). Therefore, it is essential to evaluate the DDIs between triazoles and BTK inhibitors.
Physiologically-based pharmacokinetic (PBPK) model is a mathematical model that integrated knowledge of physiology, biochemistry and anatomy, in order to simulate the absorption, distribution, metabolism and excretion (ADME) characteristics of drugs in humans (Ellison, 2018). Recently PBPK model has been increasingly accepted by regulatory agencies as a method to inform clinical research strategies. And it has become a useful tool in the simulation of multiple inducers or inhibitors, relevant metabolites, and multiple mechanisms of interaction. Therefore, it has been allowed to predict the complex DDIs involving transporters, enzymes, and multiple interaction mechanisms (Sinha et al., 2014; Sager et al., 2015). The U.S. Food and Drug Administration (FDA) Office of Clinical Pharmacology has been tracking the use of PBPK models in regulatory submissions since 2008. According to 2013 submissions, the models included in regulatory files were most commonly used for DDI (60%), pediatric (21%), and absorption (6%) predictions (Sager et al., 2015). Simcyp (version 20, Certara, Sheffield, United Kingdom), a platform and database for “bottom-up” mechanistic modeling and simulation of the processes of oral absorption, tissue distribution, metabolism and excretion of drugs and drug candidates in healthy and disease populations, is often used to develop PBPK models and to predict the pharmacokinetics and DDIs (Jamei et al., 2009).
In this study, a PBPK model was used to investigate the influence of different triazoles on the pharmacokinetics of BTK inhibitors (zanubrutinib, acalabrutinib) by Simcyp, and the DDIs were quantified to provide a general guidance for the dosage adjustment of BTK inhibitors when co-administered with triazole antifungals.
Materials and methods
Physiologically-based pharmacokinetic model development and verification of bruton’s tyrosine kinase inhibitors
A basic framework of PBPK model development and verification is presented in Figure 1. The developments of zanubrutinib and acalabrutinib PBPK models were based on clinical pharmacokinetic parameters, physicochemical properties data, and in vitro experiments parameters. The essential physicochemical properties parameters for the development of PBPK models including molecular weight, the acid dissociation constant (pKa), solubility, octanol/water partition coefficient (logP), fraction unbound in plasma (fup), fraction unbound in gut (fu,gut), blood-to-plasma concentration ratio (Rbp), and effective permeability (Peff). These physicochemical properties parameters and the corresponding references (Zane and Thakker, 2014; Qi et al., 2017; Li et al., 2018; Zhou et al., 2019; Cai et al., 2020; Li et al., 2020; Wang et al., 2021) are listed in Table 1. The first order absorption model and advanced dissolution, absorption, and metabolism (ADAM) model were used to describe the absorption processes of acalabrutinib and zanubrutinib, respectively. The minimal PBPK model and full PBPK model were used to simulate the distribution processes of acalabrutinib and zanubrutinib, respectively. The selected distribution models are based on published literatures (Zhou et al., 2019; Wang et al., 2021), and the results of model validation showed that the models are reliable and robust. For zanubrutinib, according to the human liver microsome study, the intrinsic clearance value for CYP3A is 120 μL/(minmg); an additional clearance value of 60 µL/(minmg) was inputted to account for non-CYP3A mediated clearance. The renal clearance value is 0.5 L/h (Wang et al., 2021).
FIGURE 1
TABLE 1
| Parameter | Zanubrutinib Wang et al. (2021) | Acalabrutinib Zhou et al. (2019) | Voriconazole Zane and Thakker. (2014); Qi et al. (2017); Li et al. (2018); Li et al. (2020) | Fluconazole Cai et al. (2020) | Itraconazole Cai et al. (2020) |
|---|---|---|---|---|---|
| Base pKa | 3.3 | 3.54, 5.77 | 1.6 | 1.76a | 4.28a |
| Molecular weight (g/mol) | 471.55 | 465.5 | 349.3 | 306.3a | 705.6 |
| Solubility (mg/ml) | 3.2 | 1.39 | 0.00964 | ||
| Rbp | 0.804 | 0.787 | 1 | 1a | 0.58a |
| Peff (×10−4cm/s) | 0.9 | 4 | 3.8 | - | 0.28 |
| logP | 4.2 | 2.03 | 1.8 | 0.2a | 4.47a |
| Papp,caco-2 (×10−6cm/s) | - | - | - | 29.8a | - |
| fu,gut (%) | - | 2.6 | 89a | 1.6a | |
| fup (%) | 5.82 | 2.6 | 42 | 89a | 1.6a |
| CYP3A4 Km (μM) | - | - | - | - | 0.004 |
| CYP3A4 Vmax [pmol/(min·pmol)] | - | - | - | - | 0.065 |
| CLint [μL/(min·mg)] | 120 | 9.63μL/min/pmol | - | - | |
| Hepatic CLint [μL/(min·mg)] | - | - | 4.3 | - | - |
| Additional clearance HLM [μL/(min·mg)] | 60 | 289.5 | - | - | |
| CLR (L/h) | 0.5 | 1.33 | 0.096 | 0.86a | - |
| CYP3A4 Ki | - | - | 0.66 μM | 10.7 μMa | 0.001 μMa |
Physicochemical property values used for PBPK modeling of zanubrutinib, acalabrutinib, voriconazole, fluconazole and itraconazole.
From Simcyp Data Management system.
pKa, acid dissociation constant; Rbp, blood-to-plasma concentration ratio; Peff, effective permeability; logP, octanol/water partition coefficient; Papp,caco-2, apparent permeability of Caco-2 cell line; fu,gut, fraction unbound in gut; fup, fraction unbound in plasma; Km, MichaelisMenten constant; Vmax, maximum rate of metabolism formation; CLint, intrinsic clearance; CLR, renal clearance.
After the PBPK models were developed, simulations were performed at doses of 80 mg zanubrutinib capsule and 100 mg acalabrutinib capsule which were based on the conventional clinical administration regimens. The time-concentration curves were simulated by PBPK models and the maximum plasma concentration (Cmax) is calculated as the peak concentration in the curve and area under the plasma concentration-time curve (AUC) integrated from 0.00 to t is calculated using log-linear trapezoidal rule in Simcyp. Specifically, Simcyp calculates AUC from 0.00 to t as where n is the number of time points in which and . The rule for is as follows. If , the log-down formula is used to calculate . Otherwise, the linear-up formula is applied as . The developed PBPK models were verified by comparing the simulated plasma concentration curves and pharmacokinetic parameters with corresponding clinically observed plasma concentration curves and pharmacokinetic data in healthy adults which based on existing published studies (Podoll et al., 2019; Ou et al., 2020). The observed data was extracted by applying GetData Graph Digitizer (http://getdata-graph-digitizer.com/). GetData Graph Digitizer is software used to digitize and extract sufficient data (Giang et al., 2019; Shen et al., 2021). The fold-error was used to assess the credibility of the developed PBPK models. The developed PBPK models were considered credible only when the fold-error was less than 2 (Cai et al., 2020). If the observed value is greater than the predicted value, fold-error = observed/predicted; if the observed value is smaller than the predicted value, fold-error = predicted/observed (Fan et al., 2019).
Physiologically-based pharmacokinetic model development and verification of triazole antifungal agents
The PBPK models developed for triazole antifungal agents were similar to the BTK inhibitors. Voriconazole, fluconazole and itraconazole are all described as inhibitors of CYP3A4 (Bellmann and Smuszkiewicz, 2017). The physicochemical properties parameters used in PBPK models and the corresponding references (Zane and Thakker, 2014; Qi et al., 2017; Li et al., 2018; Zhou et al., 2019; Cai et al., 2020; Li et al., 2020; Wang et al., 2021) are listed in Table 1. The absorption processes of voriconazole, fluconazole and itraconazole were dscribed using the first order absorption models. The distribution processes of voriconazole, fluconazole and itraconazole were dscribed using full PBPK model, minimal PBPK model and minimal PBPK model, respectively. The recombinant enzyme and kinetic parameters [Michaelis-Menten constant (Km) and maximum reaction velocity (Vmax)] were used to describe the metabolic process of drugs. The apparent Km and Vmax values of itraconazole were 0.004 μM and 0.065 pmol/(minpmol) for CYP3A4, respectively. The essential parameters of voriconazole, fluconazole and itraconazole were listed in Table 1. The accuracy of developed PBPK models were verified by comparing the simulated plasma concentration curves and pharmacokinetic parameters with corresponding clinically observed data (Thorpe et al., 1990; Jaruratanasirikul and Sriwiriyajan, 1998; Purkins et al., 2002).
Drug-drug interactions simulations of bruton’s tyrosine kinase inhibitors and triazole antifungal agents
After the verification, the PBPK model was used to simulate clinical DDI scenarios to quantitatively evaluate the pharmacokinetic changes of zanubrutinib or acalabrutinib when co-administered with triazoles. For the simulation of single dose, all virtual volunteers were given zanubrutinib capsule 160 mg or acalabrutinib capsule 100 mg, combined with 200 mg voriconazole or 200 mg fluconazole or 200 mg itraconazole orally. For the simulation of multiple doses zanubrutinib, the virtual volunteers were given 160 mg zanubrutinib capsule twice daily concomitantly with 200 mg fluconazole once daily for 14 days or 200 mg itraconazole once-daily for 14 days or voriconazole at a loading dose of 400 mg twice-daily (day 1) and a subsequent dose of 200 mg twice-daily (days 2–14). For acalabrutinib group, the virtual volunteers were given 100 mg acalabrutinib capsule twice daily concomitantly with 200 mg fluconazole once daily for 7 days or 200 mg itraconazole once-daily for 7 days or voriconazole at a loading dose of 400 mg twice-daily (day 1) and a subsequent dose of 200 mg twice-daily (days 2–7). The inhibitory potency of triazole antifungals can be measured by the inhibition constant (Ki) value. The Ki values of triazole antifungals were entered into PBPK models to predict the potential DDIs. The Ki values of voriconazole, fluconazole and itraconazole were laid in Table 1.
Results
Physiologically-based pharmacokinetic model development and verification of bruton’s tyrosine kinase inhibitors and triazole antifungal agents
The robustness of the PBPK models were assessed by comparing predicted with corresponding clinically observed plasma concentration-time profiles and pharmacokinetic parameters (Thorpe et al., 1990; Jaruratanasirikul and Sriwiriyajan, 1998; Purkins et al., 2002; Podoll et al., 2019; Ou et al., 2020). As presented in Figure 2, the predicted plasma concentration curves of zanubrutinib, acalabrutinib, voricoanzole, fluconazole and itraconazole were consistent with the observed curves. Besides, the Cmax and AUC values were successfully predicted with fold-errors ≤ 2. The Cmax and AUC values of zanubrutinib, acalabrutinib, voricoanzole, fluconazole and itraconazole and the fold-error values are presented in Table 2. It is obvious that the developed PBPK models are credible.
FIGURE 2
TABLE 2
| Cmax (ng/ml) | Tmax (h) | AUC (ng·h/mL)* | ||
|---|---|---|---|---|
| Zanubrutinib 80 mg | Observed | 162.8 | 1.5 | 663 |
| Predicted | 108 | 1.68 | 1030 | |
| Fold-error | 1.51 | 1.12 | 1.55 | |
| Acalabrutinib 100 mg | Observed | 639 | 0.5 | 643 |
| Predicted | 390 | 0.56 | 491 | |
| Fold-error | 1.64 | 1.12 | 1.31 | |
| Voriconazole 300 mg | Observed | 2360 | 1.41 | 12650 |
| Predicted | 2300 | 0.99 | 21800 | |
| Fold-error | 1.03 | 1.42 | 1.72 | |
| Fluconazole 100 mg | Observed | 1700 | 4.29 | 93000 |
| Predicted | 1560 | 2.49 | 75200 | |
| Fold-error | 1.09 | 1.72 | 1.24 | |
| Itraconazole 200 mg | Observed | 280 | 4.36 | 1970 |
| Predicted | 201 | 3.24 | 1930 | |
| Fold-error | 1.39 | 1.35 | 1.02 | |
Observed and predicted PK parameters of zanubrutinib, acalabrutinib, voriconazole, fluconazole and itraconazole.
*AUClast for zanubrutinb, acalabrutinib and voriconazole; AUCinf for fluconazole; AUC24 for itraconazole (single dose).
PK, pharmacokinetics; AUC, area under the plasma concentration-time curve; Cmax, maximum plasma concentration; Tmax, time-to-maximum plasma concentration.
Drug-drug interactions simulations of bruton’s tyrosine kinase inhibitors and triazole antifungal agents
The developed PBPK model was applied to predict clinical DDI scenarios for zanubrutinib or acalabrutinib when co-administered with triazole antifungal agents. The simulated DDI results are presented in Table 3, Table 4, Figure 3, Figure 4, Figure 5 and Figure 6. The results indicate that exposures of zanubrutinib and acalabrutinib may increase when co-administered with triazole antifungals. The Cmax of zanubrutinib increased by 94%, 60%, and 34% and the AUC increased by 127%, 81%, and 48% when co-administered with voriconazole, fluconazole or itraconazole at multiple doses, respectively. The Cmax of acalabrutinib increased by 220%, 93%, and 200% and the AUC increased by 326%, 119% and 264% when co-administered with voriconazole, fluconazole or itraconazole at multiple doses, respectively. Compared with fluconazole and itraconazole, voriconazole exhibited the greatest influence on exposures of zanubrutinib and acalabrutinib.
TABLE 3
| Compound | Parameters | |||
|---|---|---|---|---|
| Cmax (ng/ml) | Tmax (h) | AUC (ng·h/mL) | ||
| Zanubrutinib | Alone (single dose) | 161 | 1.44 | 1290 |
| DDI with voriconazole (single dose) | 238 | 1.44 | 2200 | |
| Ratio with voriconazole (single dose) | 1.48 | 1.00 | 1.71 | |
| Alone (multiple doses) | 216 | 1.92 | 1580 | |
| DDI with voriconazole (multiple doses) | 419 | 1.92 | 3580 | |
| Ratio with voriconazole (multiple doses) | 1.94 | 1.00 | 2.27 | |
| Alone (single dose) | 161 | 1.44 | 1290 | |
| DDI with fluconazole (single dose) | 197 | 1.44 | 1760 | |
| Ratio with fluconazole (single dose) | 1.22 | 1.00 | 1.36 | |
| Alone (multiple doses) | 216 | 1.92 | 1580 | |
| DDI with fluconazole (multiple doses) | 345 | 1.92 | 2860 | |
| Ratio with fluconazole (multiple doses) | 1.60 | 1.00 | 1.81 | |
| Alone (single dose) | 164 | 1.44 | 1350 | |
| DDI with itraconazole (single dose) | 243 | 1.44 | 2250 | |
| Ratio with itraconazole (single dose) | 1.48 | 1.00 | 1.67 | |
| Alone (multiple doses) | 222 | 1.92 | 1640 | |
| DDI with itraconazole (multiple doses) | 299 | 1.92 | 2430 | |
| Ratio with itraconazole (multiple doses) | 1.34 | 1.00 | 1.48 | |
Model-predicted PK parameters and ratios of zanubrutinib given alone and with triazoles.
PK, pharmacokinetics; DDI, drug-drug interaction.
TABLE 4
| Compound | Parameters | |||
|---|---|---|---|---|
| Cmax (ng/ml) | Tmax (h) | AUC (ng·h/mL) | ||
| Acalabrutinib | Alone (single dose) | 385 | 0.6 | 513 |
| DDI with voriconazole (single dose) | 1170 | 0.6 | 1930 | |
| Ratio with voriconazole (single dose) | 3.04 | 1.00 | 3.76 | |
| Alone (multiple doses) | 402 | 1.08 | 513 | |
| DDI with voriconazole (multiple doses) | 1286 | 1.08 | 2184 | |
| Ratio with voriconazole (multiple doses) | 3.20 | 1.00 | 4.26 | |
| Alone (single dose) | 385 | 0.6 | 513 | |
| DDI with fluconazole (single dose) | 658 | 0.6 | 937 | |
| Ratio with fluconazole (single dose) | 1.71 | 1.00 | 1.83 | |
| Alone (multiple doses) | 402 | 1.08 | 513 | |
| DDI with fluconazole (multiple doses) | 776 | 1.08 | 1124 | |
| Ratio with fluconazole (multiple doses) | 1.93 | 1.00 | 2.19 | |
| Alone (single dose) | 387 | 0.6 | 512 | |
| DDI with itraconazole (single dose) | 1160 | 0.6 | 1790 | |
| Ratio with itraconazole (single dose) | 3.00 | 1.00 | 3.50 | |
| Alone (multiple doses) | 404 | 1.08 | 513 | |
| DDI with itraconazole (multiple doses) | 1213 | 1.08 | 1865 | |
| Ratio with itraconazole (multiple doses) | 3.00 | 1.00 | 3.64 | |
Model-predicted PK parameters and ratios of acalabrutinib given alone and with triazoles.
PK, pharmacokinetics; DDI, drug-drug interaction.
FIGURE 3
FIGURE 4
FIGURE 5
FIGURE 6
Discussion
The results of the DDI simulations showed that the pharmacokinetic exposures of zanubrutinib and acalabruitnib increased to varying degrees when combined with voriconazole, fluconazole, or itraconazole, respectively. In brief, compared with taking zanubrutinib alone, the AUC of zanubrutinib increased by 127%, 81%, and 48% when combined with voriconazole, fluconazole or itraconazole at multiple doses, respectively. Furthermore, compared with taking acalabrutinib alone, the AUC of acalabrutinib increased by 326%, 119%, and 264% when combined with voriconazole, fluconazole or itraconazole at multiple doses, respectively.
According to the results above, co-administered of BTK inhibitors and triazoles will increase the pharmacokinetic exposures of BTK inhibitors, and among the three triazoles, voriconazole exhibited the most significant effect on the pharmacokinetic exposures of zanubrutinib and acalabruitnib. Nonetheless, the degree of elevation was markedly different between zanubrutinib and acalatinib, especially co-administered with voriconazole and itraconazole. The reason may be related to the fact that zanubrutinib can decrease the systemic exposure of CYP3A and CYP2C19 substrates (Ou et al., 2021). Voriconazole, which happens to be a substrate for CYP2C19, CYP2C9 and CYP3A4, and itraconazole is a substrate for CYP3A4 (Bellmann and Smuszkiewicz, 2017). Therefore, zanubrutinib decreased the systemic exposures of voriconazole and itraconazole, resulting in less inhibitory effects on zanubrutinib caused by voriconazole and itraconazole compared with acalabrutinib. Whereas fluconazole’s metabolic pathways are not qualitatively or quantitatively significant, and its main route of elimination is renal excretion (Debruyne and Ryckelynck, 1993), which will not be influenced by zanubrutinib and acalabruitnib, so both of the pharmacokinetic exposures increased in similar degree.
Therapeutic drug monitoring (TDM) is the clinical practice of measuring drugs at specified time intervals to support individualized PK-based dose adjustments, thus maintaining consistent concentrations in patient’s blood, reducing regimen-related toxicities and improving treatment efficacy. TDM has been shown its advantage in optimization the dosing of voriconazole (Ashbee et al., 2014), vancomycin (Pai et al., 2014), valproic acid (Johannessen Landmark et al., 2020), cyclosporine (Jorga et al., 2004) and so on. Moreover, the exposure-response and/or exposure-toxicity relationships of several oral targeted antineoplastic drugs have been established, and TDM has been proven to be practical for individualized dosing of imatinib, sunitinib, abiraterone, everolimus, etc., (Verheijen et al., 2017; Mueller-Schoell et al., 2021). Even though there has not any recommendation for TDM of the BTK inhibitors to date, TDM can still be conducted to clarify the DDIs between BTK inhibitors and triazole antifungal agents, so as to guide individualized dosing, optimize therapy and prevent toxicity. Overall, our study indicated that in order to avoid the increased concentration of BTK inhibitors, we should reduce the dosage of BTK inhibitors when co-administered with triazoles, especially voriconazole.
Although the PBPK model is well-established, reasonably refined and validated, limitations still exist in the present study. Firstly, genetic polymorphisms of CYP3A4 may alter the metabolic enzyme activities of zanubrutinib and acalabrutinib. The inhibitory potency also varies among different variants when co-administered with a CYP inhibitor (Han et al., 2021). Secondly, the DDIs between zanubrutinib, acalabrutinib and triazoles were predicted in healthy subjects in our study. However, the enzyme activity of CYP3A4 may be different in disease state such as CLL, SLL, and MCL (Gao et al., 2022). Therefore, the DDIs between zanubrutinib, acalabrutinib and triazoles in patients with hematologic malignancies need to be studied in further research.
Conclusion
In conclusion, the developed and validated PBPK models were successfully used to predict the DDIs between zanubrutinib, acalabrutinib and different triazoles. Compared with taking zanubrutinib or acalabrutinib alone, the pharmacokinetic exposures of zanubrutinib and acalabruitnib increased to varying degrees when co-administered with voriconazole, fluconazole, or itraconazole, respectively. The dosage of zanubrutinib and acalabrutinib need to be reduced when co-administered with triazole antifungal agents.
Statements
Data availability statement
The original contributions presented in the study are included in the article/supplementary material, further inquiries can be directed to the corresponding author.
Author contributions
LC, CL, and WC contributed to conception and design of the study. CL collected the data. LC performed the statistical analysis. LC wrote the first draft of the manuscript. CL, HB, LL, WC and LC wrote sections of the manuscript. All authors contributed to manuscript revision, read, and approved the submitted version.
Funding
The authors appreciate the support of the Fundamental Research Funds for the Central Universities (2021CDJYGRH-014), the Natural Science Foundation of Chongqing, China (cstc2021jcyj-msxmX1154) and the Chongqing Key Specialty Construction Project of Clinical Pharmacy.
Acknowledgments
Certara United Kingdom (Simcyp Division) granted free access to the Simcyp® Simulators through an academic licence.
Conflict of interest
The 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.
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.
References
1
AbbasH. A.WierdaW. G. (2021). Acalabrutinib: A selective Bruton tyrosine kinase inhibitor for the treatment of B-cell malignancies. Front. Oncol.11, 668162. 10.3389/fonc.2021.668162
2
AshbeeH. R.BarnesR. A.JohnsonE. M.RichardsonM. D.GortonR.HopeW. W. (2014). Therapeutic drug monitoring (TDM) of antifungal agents: Guidelines from the British society for medical mycology. J. Antimicrob. Chemother.69, 1162–1176. 10.1093/jac/dkt508
3
BellmannR.SmuszkiewiczP. (2017). Pharmacokinetics of antifungal drugs: Practical implications for optimized treatment of patients. Infection45, 737–779. 10.1007/s15010-017-1042-z
4
BruggemannR. J.VerheggenR.BoerrigterE.StanzaniM.VerweijP. E.BlijlevensN. M. A.et al (2022). Management of drug-drug interactions of targeted therapies for haematological malignancies and triazole antifungal drugs. Lancet. Haematol.9, e58–e72. 10.1016/S2352-3026(21)00232-5
5
BurgerJ. A. (2019). Bruton tyrosine kinase inhibitors: Present and future. Cancer J.25, 386–393. 10.1097/PPO.0000000000000412
6
CaiT.LiaoY.ChenZ.ZhuY.QiuX. (2020). The influence of different triazole antifungal agents on the pharmacokinetics of cyclophosphamide. Ann. Pharmacother.54, 676–683. 10.1177/1060028019896894
7
DebruyneD.RyckelynckJ. P. (1993). Clinical pharmacokinetics of fluconazole. Clin. Pharmacokinet.24, 10–27. 10.2165/00003088-199324010-00002
8
EllisonC. A. (2018). Structural and functional pharmacokinetic analogs for physiologically based pharmacokinetic (PBPK) model evaluation. Regul. Toxicol. Pharmacol.99, 61–77. 10.1016/j.yrtph.2018.09.008
9
FanJ.ChenL.LuX.LiM.ZhuL. (2019). The pharmacokinetic prediction of cyclosporin A after coadministration with wuzhi capsule. Aaps Pharmscitech20, 247. 10.1208/s12249-019-1444-6
10
GaoN.ZhangX.HuX.KongQ.CaiJ.HuG.et al (2022). The influence of CYP3A4 genetic polymorphism and proton pump inhibitors on osimertinib metabolism. Front. Pharmacol.13, 794931. 10.3389/fphar.2022.794931
11
GiangH. T. N.AhmedA. M.FalaR. Y.KhattabM. M.OthmanM. H. A.AbdelrahmanS. A. M.et al (2019). Methodological steps used by authors of systematic reviews and meta-analyses of clinical trials: A cross-sectional study. BMC Med. Res. Methodol.19, 164. 10.1186/s12874-019-0780-2
12
HamalainenS.KuittinenT.MatinlauriI.NousiainenT.KoivulaI.JantunenE. (2008). Neutropenic fever and severe sepsis in adult acute myeloid leukemia (AML) patients receiving intensive chemotherapy: Causes and consequences. Leuk. Lymphoma49, 495–501. 10.1080/10428190701809172
13
HanM.QianJ.YeZ.XuR.ChenD.XieS.et al (2021). Functional assessment of the effects of CYP3A4 variants on acalabrutinib metabolism in vitro. Chem. Biol. Interact.345, 109559. 10.1016/j.cbi.2021.109559
14
Hardy-AbeloosC.PinottiR.GabriloveJ. (2020). Ibrutinib dose modifications in the management of CLL. J. Hematol. Oncol.13, 66. 10.1186/s13045-020-00870-w
15
HeidenreichD.HansenE.KreilS.NolteF.JawharM.HechtA.et al (2022). The insertion site is the main risk factor for central venous catheter-related complications in patients with hematologic malignancies. Am. J. Hematol.97, 303–310. 10.1002/ajh.26445
16
JameiM.MarciniakS.FengK.BarnettA.TuckerG.Rostami-HodjeganA. (2009). The Simcyp population-based ADME simulator. Expert Opin. Drug Metab. Toxicol.5, 211–223. 10.1517/17425250802691074
17
JaruratanasirikulS.SriwiriyajanS. (1998). Effect of omeprazole on the pharmacokinetics of itraconazole. Eur. J. Clin. Pharmacol.54, 159–161. 10.1007/s002280050438
18
Johannessen LandmarkC.JohannessenS. I.PatsalosP. N. (2020). Therapeutic drug monitoring of antiepileptic drugs: Current status and future prospects. Expert Opin. Drug Metab. Toxicol.16, 227–238. 10.1080/17425255.2020.1724956
19
JorgaA.HoltD. W.JohnstonA. (2004). Therapeutic drug monitoring of cyclosporine. Transpl. Proc.36, 396S–403S. 10.1016/j.transproceed.2004.01.013
20
LaniniS.MolloyA. C.FineP. E.PrenticeA. G.IppolitoG.KibblerC. C. (2011). Risk of infection in patients with lymphoma receiving rituximab: Systematic review and meta-analysis. BMC Med.9, 36. 10.1186/1741-7015-9-36
21
LewisR. E.Cahyame-ZunigaL.LeventakosK.ChamilosG.Ben-AmiR.TamboliP.et al (2013). Epidemiology and sites of involvement of invasive fungal infections in patients with haematological malignancies: A 20-year autopsy study. Mycoses56, 638–645. 10.1111/myc.12081
22
LiN.ZhuL.QiF.LiM.XuG.GeT. (2018). Prediction of the effect of voriconazole on the pharmacokinetics of non-steroidal anti-inflammatory drugs. J. Chemother.30, 240–246. 10.1080/1120009X.2018.1500197
23
LiX.FrechenS.MojD.LehrT.TaubertM.HsinC. H.et al (2020). A physiologically based pharmacokinetic model of voriconazole integrating time-dependent inhibition of CYP3A4, genetic polymorphisms of CYP2C19 and predictions of drug-drug interactions. Clin. Pharmacokinet.59, 781–808. 10.1007/s40262-019-00856-z
24
LipskyA.LamannaN. (2020). Managing toxicities of Bruton tyrosine kinase inhibitors. Hematol. Am. Soc. Hematol. Educ. Program2020, 336–345. 10.1182/hematology.2020000118
25
MartinoR.LopezR.SuredaA.BrunetS.Domingo-AlbosA. (1997). Risk of reactivation of a recent invasive fungal infection in patients with hematological malignancies undergoing further intensive chemo-radiotherapy. A single-center experience and review of the literature. Haematologica82, 297–304.
26
Mueller-SchoellA.GroenlandS. L.Scherf-ClavelO.van DykM.HuisingaW.MicheletR.et al (2021). Correction to: Therapeutic drug monitoring of oral targeted antineoplastic drugs. Eur. J. Clin. Pharmacol.77, 465. 10.1007/s00228-020-03067-9
27
NeofytosD.LuK.Hatfield-SeungA.BlackfordA.MarrK. A.TreadwayS.et al (2013). Epidemiology, outcomes, and risk factors of invasive fungal infections in adult patients with acute myelogenous leukemia after induction chemotherapy. Diagn. Microbiol. Infect. Dis.75, 144–149. 10.1016/j.diagmicrobio.2012.10.001
28
OuY. C.PrestonR. A.MarburyT. C.TangZ.NovotnyW.TawashiM.et al (2020). A phase 1, open-label, single-dose study of the pharmacokinetics of zanubrutinib in subjects with varying degrees of hepatic impairment. Leuk. Lymphoma61, 1355–1363. 10.1080/10428194.2020.1719097
29
OuY. C.TangZ.NovotnyW.TawashiM.LiT. K.ColemanH. A.et al (2021). Evaluation of drug interaction potential of zanubrutinib with cocktail probes representative of CYP3A4, CYP2C9, CYP2C19, P-gp and BCRP. Br. J. Clin. Pharmacol.87, 2926–2936. 10.1111/bcp.14707
30
PaiM. P.NeelyM.RodvoldK. A.LodiseT. P. (2014). Innovative approaches to optimizing the delivery of vancomycin in individual patients. Adv. Drug Deliv. Rev.77, 50–57. 10.1016/j.addr.2014.05.016
31
PappasP. G.KauffmanC. A.AndesD. R.ClancyC. J.MarrK. A.Ostrosky-ZeichnerL.et al (2016). Executive summary: Clinical practice guideline for the management of candidiasis: 2016 update by the infectious diseases society of America. Clin. Infect. Dis.62, 409–417. 10.1093/cid/civ1194
32
PattersonT. F.ThompsonG. R.DenningD. W.FishmanJ. A.HadleyS.HerbrechtR.et al (2016). Practice guidelines for the diagnosis and management of aspergillosis: 2016 update by the infectious diseases society of America. Clin. Infect. Dis.63, e1–e60. 10.1093/cid/ciw326
33
PerfectJ. R.DismukesW. E.DromerF.GoldmanD. L.GraybillJ. R.HamillR. J.et al (2010). Clinical practice guidelines for the management of cryptococcal disease: 2010 update by the infectious diseases society of America. Clin. Infect. Dis.50, 291–322. 10.1086/649858
34
PodollT.PearsonP. G.EvartsJ.IngallineraT.BibikovaE.SunH.et al (2019). Bioavailability, biotransformation, and excretion of the covalent Bruton tyrosine kinase inhibitor acalabrutinib in rats, dogs, and humans. Drug Metab. Dispos.47, 145–154. 10.1124/dmd.118.084459
35
PurkinsL.WoodN.GhahramaniP.GreenhalghK.AllenM. J.KleinermansD. (2002). Pharmacokinetics and safety of voriconazole following intravenous- to oral-dose escalation regimens. Antimicrob. Agents Chemother.46, 2546–2553. 10.1128/AAC.46.8.2546-2553.2002
36
QiF.ZhuL.LiN.GeT.XuG.LiaoS. (2017). Influence of different proton pump inhibitors on the pharmacokinetics of voriconazole. Int. J. Antimicrob. Agents49, 403–409. 10.1016/j.ijantimicag.2016.11.025
37
SagerJ. E.YuJ.Ragueneau-MajlessiI.IsoherranenN. (2015). Physiologically based pharmacokinetic (PBPK) modeling and simulation approaches: A systematic review of published models, applications, and model verification. Drug Metab. Dispos.43, 1823–1837. 10.1124/dmd.115.065920
38
ShenT.YueY.HeT.HuangC.QuB.LvW.et al (2021). The association between the gut microbiota and Parkinson's disease, a meta-analysis. Front. Aging Neurosci.13, 636545. 10.3389/fnagi.2021.636545
39
SinhaV.ZhaoP.HuangS. M.ZinehI. (2014). Physiologically based pharmacokinetic modeling: From regulatory science to regulatory policy. Clin. Pharmacol. Ther.95, 478–480. 10.1038/clpt.2014.46
40
TamC. S.OuY. C.TrotmanJ.OpatS. (2021). Clinical pharmacology and PK/PD translation of the second-generation Bruton's tyrosine kinase inhibitor, zanubrutinib. Expert Rev. Clin. Pharmacol.14, 1329–1344. 10.1080/17512433.2021.1978288
41
ThorpeJ. E.BakerN.Bromet-PetitM. (1990). Effect of oral antacid administration on the pharmacokinetics of oral fluconazole. Antimicrob. Agents Chemother.34, 2032–2033. 10.1128/AAC.34.10.2032
42
VerheijenR. B.YuH.SchellensJ. H. M.BeijnenJ. H.SteeghsN.HuitemaA. D. R. (2017). Practical recommendations for therapeutic drug monitoring of kinase inhibitors in oncology. Clin. Pharmacol. Ther.102, 765–776. 10.1002/cpt.787
43
WangK.YaoX.ZhangM.LiuD.GaoY.SahasranamanS.et al (2021). Comprehensive PBPK model to predict drug interaction potential of Zanubrutinib as a victim or perpetrator. CPT. Pharmacometrics Syst. Pharmacol.10, 441–454. 10.1002/psp4.12605
44
ZaneN. R.ThakkerD. R. (2014). A physiologically based pharmacokinetic model for voriconazole disposition predicts intestinal first-pass metabolism in children. Clin. Pharmacokinet.53, 1171–1182. 10.1007/s40262-014-0181-y
45
ZengH.WuZ.YuB.WangB.WuC.WuJ.et al (2021). Network meta-analysis of triazole, polyene, and echinocandin antifungal agents in invasive fungal infection prophylaxis in patients with hematological malignancies. BMC Cancer21, 404. 10.1186/s12885-021-07973-8
46
ZhouD.PodollT.XuY.MoorthyG.VishwanathanK.WareJ.et al (2019). Evaluation of the drug-drug interaction potential of acalabrutinib and its active metabolite, ACP-5862, using a physiologically-based pharmacokinetic modeling approach. CPT. Pharmacometrics Syst. Pharmacol.8, 489–499. 10.1002/psp4.12408
Summary
Keywords
BTK inhibitors, voriconazole, fluconazole, itraconazole, drug-drug interactions, physiologically-based pharmacokinetic
Citation
Chen L, Li C, Bai H, Li L and Chen W (2022) Use of modeling and simulation to predict the influence of triazole antifungal agents on the pharmacokinetics of zanubrutinib and acalabrutinib. Front. Pharmacol. 13:960186. doi: 10.3389/fphar.2022.960186
Received
02 June 2022
Accepted
26 September 2022
Published
10 October 2022
Volume
13 - 2022
Edited by
Oscar Garcia-Algar, Hospital Clinic of Barcelona, Spain
Reviewed by
Ronette Gehring, Utrecht University, Netherlands
Ning Ji, Tianjin Medical University Cancer Institute and Hospital, China
Updates
Copyright
© 2022 Chen, Li, Bai, Li and Chen.
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: Wanyi Chen, 32003371@qq.com
†These authors have contributed equally to this work
This article was submitted to Drug Metabolism and Transport, a section of the journal Frontiers in Pharmacology
Disclaimer
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.