<figure><img src="image_1.png" /></figure>

Article

# From Personalized to Precision Medicine in Oncology:
A Model-Based Dosing Approach to Optimize Achievement of Imatinib Target Exposure

Sylvain Goutelle ¹,²,³,* ID, Monia Guidi ⁴,⁵,⁶, Verena Gotta ⁷, Chantal Csajka ⁵,⁶,⁸ ID, Thierry Buclin ⁴ ID
and Nicolas Widmer ⁴,⁶,⁹ ID

1 Service de Pharmacie, GH Nord, Hospices Civils de Lyon, 69002 Lyon, France

2 Univ. Lyon, Université Claude Bernard Lyon 1, UMR CNRS 5558, LBBE—Laboratoire de Biométrie et Biologie Évolutive, 69100 Villeurbanne, France

3 Univ. Lyon, Université Claude Bernard Lyon 1, ISPB—Faculté de Pharmacie de Lyon, 69008 Lyon, France

4 Service of Clinical Pharmacology, Lausanne University Hospital and University of Lausanne, 1011 Lausanne, Switzerland; monia.guidi@chuv.ch (M.G.); nicolas.widmer@chuv.ch (N.W.)

5 Center for Research and Innovation in Clinical Pharmaceutical Sciences, University Hospital and University of Lausanne, 1011 Lausanne, Switzerland

6 Institute of Pharmaceutical Sciences of Western Switzerland, University of Geneva and University of Lausanne, 1211 Geneva, Switzerland

7 Pediatric Pharmacology and Pharmacometrics, University of Basel Children's Hospital, 4056 Basel, Switzerland

8 School of Pharmaceutical Sciences, University of Geneva, 1205 Geneva, Switzerland

9 Pharmacy of the Eastern Vaud Hospitals, 1847 Rennaz, Switzerland

* Correspondence: sylvain.goutelle@chu-lyon.fr; Tel.: +33-4-72-16-80-99

<figure><img src="image_2.png" /></figure>

check for
updates

**Citation:** Goutelle, S.; Guidi, M.; Gotta, V.; Csajka, C.; Buclin, T.; Widmer, N. From Personalized to Precision Medicine in Oncology: A Model-Based Dosing Approach to Optimize Achievement of Imatinib Target Exposure. *Pharmaceutics* **2023**, *15*, 1081. https://doi.org/10.3390/pharmaceutics15041081

Academic Editor: Federico Pea

Received: 28 January 2023

Revised: 23 March 2023

Accepted: 27 March 2023

Published: 28 March 2023

<figure><img src="image_3.png" /></figure>

**Copyright:** © 2023 by the authors.
Licensee MDPI, Basel, Switzerland.
This article is an open access article
distributed under the terms and
conditions of the Creative Commons
Attribution (CC BY) license (https://
creativecommons.org/licenses/by/
4.0/).

**Abstract:** Imatinib is a targeted cancer therapy that has significantly improved the care of patients with chronic myeloid leukemia (CML) and gastrointestinal stromal tumor (GIST). However, it has been shown that the recommended dosages of imatinib are associated with trough plasma concentration (Cmin) lower than the target value in many patients. The aims of this study were to design a novel model-based dosing approach for imatinib and to compare the performance of this method with that of other dosing methods. Three target interval dosing (TID) methods were developed based on a previously published PK model to optimize the achievement of a target Cmin interval or minimize underexposure. We compared the performance of those methods to that of traditional model-based target concentration dosing (TCD) as well as fixed-dose regimen using simulated patients (n = 800) as well as real patients' data (n = 85). Both TID and TCD model-based approaches were effective with about 65% of Cmin achieving the target imatinib Cmin interval of 1000–2000 ng/mL in 800 simulated patients and more than 75% using real data. The TID approach could also minimize underexposure. The standard 400 mg/24 h dosage of imatinib was associated with only 29% and 16.5% of target attainment in simulated and real conditions, respectively. Some other fixed-dose regimens performed better but could not minimize over- or underexposure. Model-based, goal-oriented methods can improve initial dosing of imatinib. Combined with subsequent TDM, these approaches are a rational basis for precision dosing of imatinib and other drugs with exposure-response relationships in oncology.

**Keywords:** imatinib; pharmacokinetics; model-informed precision dosing; oncology; antineoplastic agents

## 1. Introduction

Imatinib is a selective and potent tyrosine kinases inhibitor approved notably for the treatment of chronic myeloid leukemia (CML) and gastrointestinal stromal tumor (GIST). It has been a major breakthrough in cancer therapy [1], as it specifically inhibits
the genetically altered molecular targets (i.e., tyrosine kinase proteins) involved in the pathogenesis of both those two conditions and has been associated with significant increase in response rates [2,3]. Owing to its innovative mechanism of action and its clinical success, imatinib has been considered as a prototypical example of targeted drug therapy [4] and a successful application of the concept of personalized medicine [5]. Since imatinib's discovery, significant advances in this area have been made, especially in oncology [6], although the actual public health benefit expected from a generalization of this approach is still questionable [7].

In 2015, the precision medicine initiative was launched in the USA [8]. Precision medicine includes the concepts of personalized medicine (i.e., reliance on genetic or acquired biomarkers for therapeutic decisions), but it is also a step beyond as it encompasses identification, quantification and clinical applications of all potential determinants of individual response in prevention and treatment strategies. The initiative also promotes the utilization of big data in health as well as the deployment of computational tools to analyze them.

In his allocution in 2015, President Obama also stressed the importance of "figuring out the right dose of medicine" in personalized medicine (https://obamawhitehouse.archives.gov/node/333101, accessed on 12 December 2022). In line with other researchers [9,10], we believe that precision dosage design of drugs has been overlooked in personalized medicine so far, and that major progress can be made in this area. While selecting the right drug for a given individual based on his/her characteristics is important, prescribing the best dose of such drug also matters to optimize drug response. Precision dosing notably includes drug dosing based on models to better consider PK-PD variability and better achieve target objectives in a goal-oriented approach. This can be performed both a priori (without any concentration measurement) and a posteriori, with a further readjustment of the dosing schedule based on the patient's actual exposure. In this paper, we aim at showing how a priori precision dosing can be applied using imatinib as an illustrative example.

The approved dosages of imatinib in adults are 400 mg once daily in GIST and chronic phase of CML, and 600 mg in CML accelerated phase and blast crisis. The dosage may be increased up to a maximum of 800 mg per day in patients with CML. The recommended dosage regimen of imatinib in chronic phase CML, 400 mg/day, has been based on in vitro activity and a mathematical model of hematological dose-response indicating that 400 mg was sufficient to provide optimal response in most patients [1,11]. Since then, it has been shown that a higher dose may be necessary to optimize cytogenetic response and survival [12]. Indeed, a higher dosage of imatinib has been associated with earlier and higher rates of cytogenetic and molecular response in patients with CML [13–15].

Since imatinib's approval in the early 2000s, there has been mounting evidence of concentration-response relationships for both CML and GIST therapy. In patients with CML, several observational studies reported that achieving an imatinib trough concentration (Cmin) above a target of about 1000 ng/mL was associated with improved rates of hematological and molecular response [16–18]. In patients with GIST, it has been shown that Cmin values above 1100 mg/L were associated with significantly longer time to disease progression [19]. Although some adverse reactions were more frequent in patients with the highest Cmin, in the study from Larson and colleagues [16], data supporting an upper bound for imatinib concentration are limited [20].

Because of this concentration-dependent activity and large interindividual PK variability [16], it has been suggested that patients could benefit from therapeutic drug monitoring of imatinib (TDM), and real-life data appear to confirm the clinical value of TDM [20–25]. However, TDM services may not be available for all patients treated with imatinib. In addition, because of the relatively long half-life of imatinib (mean value of 19 h [11]), it may take many days to adequately measure concentrations at the steady-state, then adjust dosage and control the achievement of target Cmin. As a result, slow response to imatinib may occur because of initial suboptimal exposure.
The objectives of the present study were to design a new model-based approach for optimal targeting of imatinib Cmin and evaluate the ability of this novel approach to achieve a target Cmin interval compared with fixed-dose regimens and another model-based approach.

## 2. Materials and Methods

### 2.1. Pharmacokinetic Model

The model published by Widmer et al. [26] was used to build the new dosing method. Briefly, this model was based on the nonlinear mixed-effects population analysis of 321 plasma concentration of imatinib measured in 59 patients with CML or GIST. The final model was a one-compartment model with linear absorption and elimination. The absorption rate constant Ka was a non-random parameter, with an estimated value of 0.609 h⁻¹ in all individuals. Population median and coefficient of variation of apparent oral clearance (CL/F) and volume of distribution (V/F) were 14.3 L/h (36%) and 347 L (63%), respectively. These two parameters were assumed to have a lognormal distribution. Significant covariance was found between CL/F and V/F, with an estimated correlation coefficient of 0.8. Because imatinib oral bioavailability is close to 100% [11], we assumed that bioavailability was 100%, so CL/F and V/F will be simply denoted as CL and V in the rest of the manuscript.

Several covariates were found to influence imatinib CL and V. The relationships between typical values (TV) of CL and V and covariates were modeled as follows:

$$TV_{CL} = 14.3 + 5.42 \cdot \frac{BW - 70}{70} + 1.49 \cdot male - 1.49(1 - male) - 5.81 \cdot \frac{AGE - 50}{50} - 0.806 \cdot path + 0.806(1 - path) \quad (1)$$

$$TV_{V} = 347 + 46.2 \cdot male - 46.2(1 - male) \quad (2)$$

where TVCL is the typical value of CL, TVV is the typical value of V, BW is patient's body weight (in kg), AGE is the patient's age (in years). Male is a binary variable indicating male gender (1 if male, 0 if female), and path is a binary variable indicating the disease (1 if GIST, 0 if CML). In summary, a priori imatinib clearance increases with increasing BW, decreases with increasing age, and is greater in males and patients with CML compared with females and GIST patients, respectively. Typical imatinib V is greater in male than in female subjects.

The performance of this model in predicting imatinib concentration collected during routine TDM was assessed in a separate study, with two datasets (65 and 20 subjects) not used for model building. The results showed good overall predictive performance of the model estimates based on a Bayesian Maximum A Posteriori (MAP) calculation of individual PK parameters [27].

### 2.2. Model Discretization

Our approach for optimal targeting of a concentration interval has been described in detail in a previous publication [28]. It is inspired by previous works from D'Argenio and colleagues who used a discrete prior distribution of PK parameters to optimize model-based control of drug exposure [29,30].

The lognormal distribution of individual CL and Vd from Widmer's model can be written as follows:

$$CL_{i} = TV_{CL} \times \exp(\eta_{CLi}) \quad (3)$$

and

$$V_{i} = TV_{V} \times \exp(\eta_{Vi}) \quad (4)$$

where TVCL and TVV are as defined in Equations (1) and (2), and ηCLi and ηVi are random variables that follow a normal distribution with zero mean and variance-covariance matrix Ω² defined as follows:

$$\Omega^2 = \begin{bmatrix} \omega^2(CL) & \omega^2(CL, V) \\ \omega^2(CL, V) & \omega^2(V) \end{bmatrix} = \begin{bmatrix} 0.127 & 0.179 \\ 0.179 & 0.396 \end{bmatrix} \quad (5)$$
The joint distribution of η parameters was discretized into a grid of 81 discrete values. Those 81 values were obtained by random sampling from the corresponding bivariate normal distribution, using the multivariate normal random numbers routine of Matlab (version 2018b, Mathworks, Natick, MA, USA), and their probabilities were obtained from the corresponding probability function values. The choice of 81 as the number of grid points was inspired by a previous work from Katz and D'Argenio which showed that this number was sufficient to obtain a correct approximation of a bivariate continuous distribution [31]. Indeed, this number may be considered large, considering that a nonparametric population model of amikacin developed in 634 patients yielded a grid of 76 support points [32]. The median and variance of the discrete distribution were -0.024 (0.114) and -0.057 (0.390) for ηCLi and ηVi, respectively, with a covariance of 0.167, very similar to values of the original continuous distribution. The prior value of Ka was fixed to 0.609, in accordance with the model from Widmer et al. [26].

This grid of 81 values of random effects allowed deriving a discrete prior distribution of CLi and Vi for each patient. Indeed, Equations (1) and (2) can be used to calculate typical values of CL and V based on the patients' covariates and discrete prior distributions of CL and V can be obtained by incorporating the vectors of 81 values of random effects (ηCLi and ηVi) into Equations (3) and (4).

An example of such bivariate prior distribution of CLi and Vi is shown in Figure 1.

<figure><img src="image_4.png" /><figcaption>**Figure 1.** Discrete prior distribution of imatinib clearance and volume of distribution in a typical patient. The patient is a male with CML, his age and body weight are 50 years and 70 kg, respectively.</figcaption></figure>

The nonparametric, discrete PK model can be used to calculate a priori Cmin distribution for various imatinib dosage regimens. Let us take the example of a virtual patient, a 60-year-old woman weighing 60 kg, with CML. Based on Equations (1) and (2), her median a priori values of CL and V are 11.7 L/h and 300.8 L, respectively. However, using the discrete distribution, there is an a priori distribution of 81 values of CL and V similar to that shown in Figure 1.
Then, one can simulate the PK profile of imatinib under a given dosage regimen in this patient. The trough concentration of imatinib after multiple oral doses at the steady-state ($C_{ss_{min}}$, further denoted as $C_{min}$ for ease of reading) is given by the following equation:

$$C_{ss_{min}} = \frac{F.Dose.Ka}{V.(Ka - Ke)} \cdot \left( \frac{1}{1 - e^{-Ke.tau}} - \frac{1}{1 - e^{-Ka.tau}} \right) \quad (6)$$

where $F$ is the oral bioavailability (here assumed to be equal to 1, as previously mentioned), $Dose$ is the dose (in mg), $Ka$ and $V$ are as defined above, $tau$ is the dosing interval (in hours) and $Ke$ is the elimination rate constant, which is equal to $CL/V$.

The discrete distribution of $CL$ and $V$ parameters results in a discrete a priori distribution of $C_{min}$ for a given dosage regimen in this patient. Such a distribution is shown in Figure 2 for two dosage regimens, 400 mg/24 h and 600 mg/24 h.

<figure><img src="image_5.png" /><figcaption>Figure 2. A priori discrete distribution of imatinib trough concentration at the steady state in a virtual patient. Upper panel, 400 mg/24 h; lower panel, 600 mg/24 h. The patient is a 60-year-old woman, her weight is 60 kg and she has CML.</figcaption></figure>

## 2.3. The Target Interval Dosing (TID) Approach

For a patient who received a given dosage regimen, there are $n$ (here $n = 81$) possible $C_{min}$ values, corresponding to the $n$ vectors of PK parameters ($CL_i$, $Vd_i$) in the discrete collection, and each $C_{min_i}$ has an a priori probability, $prob_i$. The sum of probi is equal to one.

If the target interval of $C_{min}$ is defined by a lower bound $L$ and an upper bound $U$, the a priori probability that $C_{min}$ is within the interval $[L, U]$ for a given dosage regimen is given by

$$\sum_{i=0}^{n} prob_i(C_{min_i} \in [L; U]) \quad (7)$$
The optimal dosage (ODint) for achieving the target Cmin interval [L, U] is the combination of dose and dosing interval that maximizes the a priori probability as defined by Equation (7). Mathematically, this translates into

$$OD_{int}(Dose, tau) = \arg\left\{\max_{Dose,tau} [7]\right\} \qquad (8)$$

In the case of imatinib, we evaluated 14 possible initial dosage regimens of increasing intensity, with daily doses ranging from 100 to 800 mg with a 100 mg increment, and 3 dosing intervals, q8h, q12h and q24h. These 14 regimens were: 100 mg/24 h, 200 mg/24 h, 100 mg/12 h, 300 mg/24 h, 100 mg/8 h, 400 mg/24 h, 200 mg/12 h, 500 mg/24 h, 600 mg/24 h, 300 mg/12 h, 200 mg/8 h, 700 mg/24 h, 800 mg/24 h and 400 mg/12 h. The daily dose was limited to 800 mg, which is the maximum recommended daily dose for imatinib.

Instead of targeting an interval of concentration, one may be interested in minimizing the a priori probability of underexposure or overexposure. For imatinib, it is relevant to minimize underexposure, i.e., to avoid low Cmin. In such case, the probability of interest is

$$\sum_{i=0}^{n} prob_i(Cmin_i < L) \qquad (9)$$

The optimal dosage regimen (ODLmin) is the one that minimizes this probability:

$$OD_{Lmin}(Dose, tau) = \arg\left\{\min_{Dose,tau} [9]\right\} \qquad (10)$$

Because a very high dose may be necessary to achieve a zero probability of underdos- ing, it may be more clinically relevant to identify the dosage associated with a specified low but non-zero probability of underdosing, for example 5%:

$$OD_{L5\%}(Dose, tau) = \arg\left\{[9]_{Dose,tau} < 0.05\right\} \qquad (11)$$

ODint, ODLmin and OD5% are three ways to perform a general dosing approach that we have called target interval dosing (TID).

Figure 3 summarizes the various steps of the TID approach for imatinib dosing.

<figure><img src="image_6.png" /><figcaption>Figure 3. Graphical abstract of the target interval dosing (TID) approach for imatinib.</figcaption></figure>
## 2.4. Dosing Based on the Traditional Model-Based Approach

The traditional model-based dosing approach uses a single vector of PK parameters, the typical values provided by Equations (1) and (2) as a priori values and considers a single concentration value as the target. This approach will be denoted as target concentration dosing (TCD), in reference to Target Concentration Intervention (TCI) which has been used to qualify the activity consisting in model-based dose adjustment based on TDM results [33]. Most of the time, the target is fixed at the midpoint of the therapeutic interval. In our example, the dose targeting the center of the Cmin interval can be calculated by rearranging Equation (3) as follows (it is assumed that F = 1):

$$Dose^{target} = \frac{C_{min}^{target} \cdot Vd \cdot (Ka - Ke)}{Ka} \cdot \frac{1}{\frac{1}{1 - e^{-Ke.tau}} - \frac{1}{1 - e^{-Ka.tau}}} \quad (12)$$

## 2.5. Evaluation of TID and Other Dosing Methods in Simulated Patients

Virtual populations of patients, defined by PK parameters and covariate values of the imatinib model were created to assess performance the TID and alternative dosing methods.

We considered four groups for categorical covariates: females with CML, females with GIST, males with CML and males with GIST.

For each group, 200 values of body weight and age were randomly sampled from two independent truncated normal distributions, with mean ± standard deviation (minimum-maximum) values of 64 ± 11 kg (40–110) in females and 79 ± 11 kg (40–110) in males for body weight, and 62 ± 13 years (20–90) for age. These values were representative of individual data from 85 patients with GIST treated with imatinib [27]. Then, typical values of CL and V were calculated for each patient using Equations (1) and (2). True simulated parameter values were calculated using Equations (3) and (4), based on calculated typical values and random effect ($\eta_{CLi}$ and $\eta_{Vi}$) values obtained by random sampling from the bivariate joint distribution.

We evaluated the ability of the TID and TCD dosing methods, as well as standard fixed-dose approaches to achieve a common imatinib target interval of 1000–2000 ng/mL in the 4 simulated populations defined above.

For the TID method, three optimal dosages (OD) were derived for each target interval:

- ODint, the dosage maximizing the attainment of the target interval as defined in Equation (8);

- ODLmin, the dosage minimizing underexposure (Cmin < L, see Equation (10));

- ODL5%, the dosage associated with a priori probability of underexposure less than 5% (Equation (11)).

For the traditional model-based dosing approach (TCD), we used the mid-interval as the target Cmin, i.e., 1500 mg/L. The dose was calculated using Equation (12) for both 12 h and 24 h dosing intervals, rounded to the nearest hundred, and capped to a maximum of 800 mg and 400 mg for once and twice daily dosing, respectively.

We also evaluated four fixed-dose regimens in each simulated patient: 400 mg/24 h, 600 mg/24 h, 200 mg/12 h and 300 mg/12 h. Once daily dosages of 400 mg and 600 mg are dosages mentioned in imatinib's European summary of product characteristics.

For each calculated dosage, the corresponding steady-state Cmin was calculated using Equation (6), with the true simulated parameter values and the calculated dose and dosing interval. Percentages of Cmin within the target interval, below the lower bound and above the upper bound, were calculated for each dosing method.

## 2.6. Evaluation of TID and Other Dosing Methods Based on Data from Real Patients

Individual estimates of imatinib CL and V, as well as covariate values were available from two groups of 65 and 20 patients, all with GIST. These patients' characteristics as well as the Bayesian method used to estimate imatinib PK parameters and ethics committee approval have been described in a previous publication [27]. This study was approved by
the ethics committee of the Lausanne Faculty of Medicine (approval number 31/09, and informed written consent was obtained from the participants.

For each patient, at least one measured concentration of imatinib was available from routine TDM. In 20 patients, a pair of concentrations measured over the same dosing interval was available on at least one occasion. Individual values of CL and Vd were estimated for each patient on each occasion using a (MAP) method based on the model from Widmer et al. [26]. Ka was fixed at the population value of 0.609 in all patients. When multiple CL and Vd values were available from several occasions in the same patient, the mean estimate was used in subsequent calculations.

We used these MAP estimates as the “true value” to assess the ability of the various imatinib dosing methods to achieve a target Cmin interval. The target interval was fixed at 1100–2000 ng/mL in this case, as Cmin ≥ 1100 ng/mL is the target value recommended in patients with GIST [30].

A procedure similar to that previously described for simulated patients was carried out. For the TID and TCD approaches, several dosage regimens of imatinib were calculated, based on the patients’ individual covariates (body weight, age, gender and disease, here GIST), as described for simulated patients. Fixed-dose regimens were also evaluated.

We then computed the predicted steady-state Cmin corresponding to each calculated dosage, and the percentage of Cmin within and outside the desired target interval, as described previously, using Equation (6) and the individual MAP estimates of CL and Vd.

## 2.7. *Model-Based Dosing Recommendations*

We used the TID-ODint method to derive dosage recommendations according to age, body weight, sex and disease (CML or GIST). We set the Cmin target interval to 1000–2000 ng/mL to accommodate both CML and GIST. We only considered once daily and twice daily dosage, with daily dose ranging from 100 to 800 mg per day. The selected dosage was that with the lowest dose and the largest dosing interval which maximized the a priori probability to achieve the target interval, as less frequent dosing regimens are desirable for patient’s adherence [34].

## 2.8. *Software Tools*

The model discretization, random sampling for the 200-subject simulations and the calculation of the optimal dosage based on the TID approach were performed using Matlab software (version 2018b, Mathworks, Natick, MA, USA).

# 3. Results

## 3.1. *Imatinib Dosing and Target Attainment in Simulated Patients*

Table 1 shows the performance of the various dosing methods in achieving the target intervals of 1000–2000 ng/mL in all simulated patients (n = 800). Figure 4 shows the distributions of Cmin achieved with the various dosing methods evaluated.

The TID-ODint, TID-ODL5%, as well as the two TCD-based dosages (q24h and q12h) provided the highest percentages of target attainment, around 65%. Regimens based on those methods were associated with lower variability in Cmin and a better distribution of values within the target interval, compared with fixed-dose regimens. Those fixed-dose regimens provided poor rates of target attainment, especially the 400 mg/24 h regimen showing only 29% of Cmin within the target range. This regimen was also associated with a high proportion of underexposure (68.2%). The fixed-dose regimen with 600 mg/24 h or 300 mg/12 h of imatinib provided better achievement of the target interval (>50%) but were associated with large proportions of underexposure and overexposure, respectively. This confirms that any fixed-dose regimen is unlikely to be optimal for achieving an imatinib target concentration in all patients, because of PK variability.
<table><caption>**Table 1.** Imatinib Cmin and attainment of target interval (1000–2000 ng/mL) for the various dosing methods in all simulated patients.</caption><thead><tr><th></th><th>TID-OD<sub>int</sub></th><th>TID-OD<sub>Lmin</sub></th><th>TID-OD<sub>L5%</sub></th><th>TCD q24h</th><th>TCD q12h</th><th>400 mg/24 h</th><th>600 mg/24 h</th><th>200 mg/12 h</th><th>300 mg/12 g</th></tr></thead><tbody><tr><td>Number of predictions</td><td>800</td><td>800</td><td>791 <sup>a</sup></td><td>800</td><td>800</td><td>800</td><td>800</td><td>800</td><td>800</td></tr><tr><td>Mean imatinib dose per dosing interval (mg)</td><td>NA</td><td>NA</td><td>NA</td><td>683 (131)</td><td>278 (68)</td><td>400</td><td>600</td><td>200</td><td>300</td></tr><tr><td>Mean Cmin (ng/mL)</td><td>1420</td><td>2252</td><td>1624</td><td>1458</td><td>1553</td><td>902</td><td>1353</td><td>1182</td><td>1773</td></tr><tr><td>Coefficient of variation of Cmin (%)</td><td>35%</td><td>39%</td><td>35%</td><td>37%</td><td>36%</td><td>48%</td><td>48%</td><td>44%</td><td>44%</td></tr><tr><td>Median Cmin (ng/mL)</td><td>1340</td><td>2087</td><td>1515</td><td>1365</td><td>1454</td><td>804</td><td>1205</td><td>1069</td><td>1603</td></tr><tr><td>5th and 95th percentiles of Cmin (ng/mL)</td><td>755–2325</td><td>1118–3959</td><td>850–2727</td><td>739–2433</td><td>807–2626</td><td>398–1694</td><td>597–2541</td><td>559–2164</td><td>839–3245</td></tr><tr><td>Cmin within target interval (%)</td><td>66.0%</td><td>43.1%</td><td>66.2%</td><td>64.5%</td><td>64.6%</td><td>29%</td><td>54.7%</td><td>49.3%</td><td>56.6%</td></tr><tr><td>Cmin < 1000 ng/mL (%)</td><td>20.2%</td><td>2.8%</td><td>10.9%</td><td>18.9%</td><td>15.6%</td><td>68.2%</td><td>32%</td><td>43.2%</td><td>13.4%</td></tr><tr><td>Cmin > 2000 ng/mL (%)</td><td>13.8%</td><td>54.1%</td><td>22.9%</td><td>16.6%</td><td>19.8%</td><td>2.8%</td><td>13.3%</td><td>7.3%</td><td>30%</td></tr></tbody></table>

<footnote><sup>a</sup> The dosage calculation with OD<sub>L5%</sub> method was impossible in 9 simulated patients because the calculated probability of underexposure was > 5%. Abbreviations: NA, not applicable.</footnote>
<figure><img src="image_7.png" /><figcaption>**Figure 4.** Boxplot of imatinib trough concentrations achieved in simulated patients with various dosing methods targeting the interval 1000–2000 ng/mL. The box edges represent the 25th and 75th percentiles of Cmin distributions. The central mark is the median. The length of whiskers is equal to 1.5 times the interquartile range. The plus signs are outlier values. For ease of graphical display, Cmin values greater than 4500 ng/mL are all compressed in the upper area between the two grey solid lines. Of note, the proportions of Cmin values within and outside the target interval are slightly different from those shown in Table 1 because the data represented graphically are restricted to the 791 Cmin values that could be computed with all methods. Abbreviations: Cmin, trough concentration; od, once daily administration (every 24 h); td, twice daily administration (every 12 h); TCD, target concentration dosing.</figcaption></figure>

As predicted by theory, the TID-OD<sub>Lmin</sub> and TID-OD<sub>L5%</sub> methods were effective in minimizing underexposure, but this logically resulted in higher rates of overexposure, especially for the OD<sub>Lmin</sub> method. Compared with the latter method, OD<sub>L5%</sub> appears as a more reasonable compromise to minimize underexposure while limiting overexposure. Of note, the observed frequency of underexposure with this method was higher than expected (10.9% versus 5%). This discrepancy is likely due to a proportion of simulated patients with high imatinib CL and V values greater than expected based on the discrete prior distribution.

There were some differences in the results between the four patient subgroups. While the percentages of target interval achievement were relatively consistent across gender and disease groups for model-based dosages, they markedly varied for the fixed-dose regimens. For example, for the ODint, TCD q24h, 400 mg q24h and 300 mg q12h dosages, the respective proportions of Cmin within 1000–2000 ng/mL were as follows: 65%, 66%, 31% and 59% in women with CML; 68%, 65%, 49% and 41% in women with GIST; 60%, 55%, 14% and 61% in men with CML; 72%, 70%, 22% and 67% in men with GIST. The 400 mg q24h dosage was associated with higher rates of underexposure in men than in women, because of increased clearance predicted by the model.
Table 2 presents the characteristics of dosages based in the TID method. While twice daily dosing was the most frequent schedule for the methods minimizing underexposure, thrice daily and once daily dosing were found as the optimal schedule in a significant proportion of patients with each method. As expected, minimizing underexposure with the OD<sub>Lmin</sub> and OD<sub>L5%</sub> methods required larger mean doses than that maximizing attainment of the target interval (OD<sub>int</sub>).

<table><caption>**Table 2.** Characteristics of dosing regimens based on the TID approach in simulated patients.</caption><thead><tr><th></th><th></th><th>OD<sub>int</sub></th><th>OD<sub>Lmin</sub></th><th>OD<sub>L5%</sub></th></tr></thead><tbody><tr><td rowspan="3">Dosing schedule <sup>a</sup></td><td>q8h</td><td>12%</td><td>2.5%</td><td>16.2%</td></tr><tr><td>q12h</td><td>21%</td><td>95.3%</td><td>67.5%</td></tr><tr><td>q24h</td><td>67%</td><td>2.2%</td><td>16.3%</td></tr><tr><td rowspan="3">Mean dose per dosing interval</td><td>q8h</td><td>180 mg</td><td>200 mg</td><td>188 mg</td></tr><tr><td>q12h</td><td>266 mg</td><td>394 mg</td><td>299 mg</td></tr><tr><td>q24h</td><td>639 mg</td><td>611 mg</td><td>581 mg</td></tr></tbody></table>

## 3.2. Performance of Imatinib Dosing Methods Based on Real Patients' Data

Table 3 shows the demographic and PK data of the 85 patients with GIST. Table 4 summarizes the simulated imatinib Cmin and attainment of the target interval for each dosing method in those 85 patients. The model-based TID-OD<sub>int</sub> and TCD q24h methods performed best in terms of target interval attainment, but the fixed regimen with 300 mg/12 h also performed remarkably well in this setting and was associated with a very low proportion of Cmin < 1100 ng/mL.

<table><caption>**Table 3.** Characteristics of the 85 patients with GIST.</caption><thead><tr><th>Variable</th><th>Value</th></tr></thead><tbody><tr><td>Number of females/males</td><td>44/41</td></tr><tr><td>Age (years)</td><td>62 ± 13 (23–85)</td></tr><tr><td>Body weight (kg)</td><td>71.2 ± 13.2 (48–100)</td></tr><tr><td>Imatinib clearance (L/h)</td><td>12.9 ± 3.7 (2.6–24.4)</td></tr><tr><td>Imatinib volume of distribution (L)</td><td>376 ± 146 (51–717)</td></tr></tbody></table>

Table 5 provides dosage recommendations based on the TID-OD<sub>int</sub> to achieve imatinib Cmin within the 1000–2000 ng/mL interval. As predicted by the population PK model, the youngest and heaviest patients would require the larger daily dose and twice-daily administration to achieve the target Cmin, while lower doses administered once daily would be sufficient for older patients with low body weight.

<footnote><sup>a</sup> The percentages indicate the proportion of q8h, q12h, and q24h schedule for each method.</footnote>
<table><caption>**Table 4.** Predicted imatinib C<sub>min</sub> and attainment of target interval (1100–2000 ng/mL) in 85 patients with GIST.</caption><thead><tr><th></th><th>TID-OD<sub>int</sub></th><th>TID-OD<sub>Lmin</sub></th><th>TID-OD<sub>L5%</sub></th><th>TCD q24h</th><th>TCD q12h</th><th>400 mg/24 h</th><th>600 mg/24 h</th><th>200 mg/12 h</th><th>300 mg/12 h</th></tr></thead><tbody><tr><td>Mean Cmin (ng/mL)</td><td>1462</td><td>2369</td><td>1754</td><td>1514</td><td>1543</td><td>945</td><td>1417</td><td>1216</td><td>1824</td></tr><tr><td>Coefficient of variation of Cmin (%)</td><td>39%</td><td>44%</td><td>49%</td><td>37%</td><td>39%</td><td>40%</td><td>40%</td><td>44%</td><td>44%</td></tr><tr><td>Median Cmin (ng/mL)</td><td>1303</td><td>2208</td><td>1547</td><td>1413</td><td>1344</td><td>886</td><td>1329</td><td>1130</td><td>1695</td></tr><tr><td>5th and 95th percentiles of Cmin (ng/mL)</td><td>980–2144</td><td>1475–3744</td><td>1040–2944</td><td>1001–2109</td><td>966–2390</td><td>567–1488</td><td>851–2233</td><td>753–1896</td><td>1130–2845</td></tr><tr><td>Cmin within target interval (%)</td><td>75.3%</td><td>36.5%</td><td>64.7%</td><td>80.0%</td><td>72.9%</td><td>16.5%</td><td>65.9%</td><td>51.8%</td><td>72.9%</td></tr><tr><td>Cmin < 1100 ng/mL (%)</td><td>16.5%</td><td>0%</td><td>9.4%</td><td>11.8%</td><td>14.1%</td><td>82.3%</td><td>24.7%</td><td>45.9%</td><td>3.5%</td></tr><tr><td>Cmin > 2000 ng/mL (%)</td><td>8.2%</td><td>63.5%</td><td>25.9%</td><td>8.2%</td><td>13.0%</td><td>1.2%</td><td>9.4%</td><td>2.3%</td><td>23.6%</td></tr></tbody></table>

<table><caption>**Table 5.** Imatinib dosage recommendations based on the TID approach according to age, body weight, sex and disease. The suggested dosage regimens (dose in mg/dosing interval in hours) are those calculated with the TID-ODint approach to maximize the probability to achieve a Cmin target interval of 1000–2000 ng/mL at the steady-state. The color code is as follows.</caption><thead><tr><th rowspan="2">Body Weight (kg)</th><th colspan="15">Age (Years)</th></tr><tr><th>20</th><th>30</th><th>40</th><th>50</th><th>60</th><th>70</th><th>80</th><th>90</th></tr></thead><tbody><tr><td rowspan="2">40</td><td>400/12</td><td>300/12</td><td>300/12</td><td>800/24</td><td>300/12</td><td>700/24</td><td>700/24</td><td>600/24</td><td>600/24</td><td>200/12</td><td>600/24</td><td>400/24</td><td>200/12</td><td>400/24</td><td>400/24</td><td>300/24</td></tr><tr><td>300/12</td><td>700/24</td><td>800/24</td><td>600/24</td><td>700/24</td><td>200/12</td><td>600/24</td><td>500/24</td><td>500/24</td><td>400/24</td><td>500/24</td><td>300/24</td><td>400/24</td><td>100/12</td><td>300/24</td><td>200/24</td></tr><tr><td rowspan="2">50</td><td>400/12</td><td>300/12</td><td>400/12</td><td>300/12</td><td>300/12</td><td>700/24</td><td>800/24</td><td>600/24</td><td>700/24</td><td>600/24</td><td>600/24</td><td>500/24</td><td>500/24</td><td>400/24</td><td>200/12</td><td>300/24</td></tr><tr><td>300/12</td><td>800/24</td><td>800/24</td><td>700/24</td><td>800/24</td><td>600/24</td><td>700/24</td><td>500/24</td><td>600/24</td><td>400/24</td><td>500/24</td><td>400/24</td><td>400/24</td><td>300/24</td><td>400/24</td><td>100/12</td></tr><tr><td rowspan="2">60</td><td>400/12</td><td>400/12</td><td>400/12</td><td>300/12</td><td>300/12</td><td>800/24</td><td>800/24</td><td>700/24</td><td>800/24</td><td>600/24</td><td>700/24</td><td>500/24</td><td>600/24</td><td>400/24</td><td>200/12</td><td>300/24</td><td>400/24</td></tr><tr><td>400/12</td><td>300/12</td><td>300/12</td><td>800/24</td><td>300/12</td><td>700/24</td><td>700/24</td><td>200/12</td><td>600/24</td><td>500/24</td><td>600/24</td><td>400/24</td><td>200/12</td><td>300/24</td><td>400/24</td><td>300/24</td></tr><tr><td rowspan="2">70</td><td>400/12</td><td>400/12</td><td>400/12</td><td>300/12</td><td>400/12</td><td>300/12</td><td>300/12</td><td>800/24</td><td>300/12</td><td>700/24</td><td>700/24</td><td>600/24</td><td>600/24</td><td>200/12</td><td>600/24</td><td>400/24</td></tr><tr><td>400/12</td><td>300/12</td><td>400/12</td><td>300/12</td><td>300/12</td><td>700/24</td><td>800/24</td><td>600/24</td><td>700/24</td><td>200/12</td><td>600/24</td><td>500/24</td><td>500/24</td><td>400/24</td><td>500/24</td><td>300/24</td></tr><tr><td rowspan="2">80</td><td>400/12</td><td>400/12</td><td>400/12</td><td>400/12</td><td>400/12</td><td>400/12</td><td>300/12</td><td>400/12</td><td>300/12</td><td>300/12</td><td>700/24</td><td>800/24</td><td>600/24</td><td>700/24</td><td>600/24</td><td>600/24</td></tr><tr><td>400/12</td><td>400/12</td><td>400/12</td><td>300/12</td><td>300/12</td><td>800/24</td><td>800/24</td><td>700/24</td><td>800/24</td><td>600/24</td><td>700/24</td><td>500/24</td><td>600/24</td><td>400/24</td><td>500/24</td><td>400/24</td></tr><tr><td rowspan="2">90</td><td>400/12</td><td>400/12</td><td>400/12</td><td>400/12</td><td>400/12</td><td>400/12</td><td>400/12</td><td>300/12</td><td>300/12</td><td>800/24</td><td>800/24</td><td>700/24</td><td>800/24</td><td>600/24</td><td>700/24</td><td>500/24</td></tr><tr><td>400/12</td><td>400/12</td><td>400/12</td><td>400/12</td><td>400/12</td><td>400/12</td><td>400/12</td><td>300/12</td><td>400/12</td><td>300/12</td><td>300/12</td><td>800/24</td><td>300/12</td><td>700/24</td><td>200/12</td><td>600/24</td><td>500/24</td></tr><tr><td rowspan="2">100</td><td>400/12</td><td>400/12</td><td>400/12</td><td>400/12</td><td>400/12</td><td>400/12</td><td>400/12</td><td>400/12</td><td>400/12</td><td>400/12</td><td>400/12</td><td>300/12</td><td>300/12</td><td>300/12</td><td>700/24</td><td>700/24</td><td>600/24</td></tr><tr><td>400/12</td><td>400/12</td><td>400/12</td><td>400/12</td><td>400/12</td><td>400/12</td><td>400/12</td><td>400/12</td><td>400/12</td><td>400/12</td><td>400/12</td><td>400/12</td><td>400/12</td><td>300/12</td><td>300/12</td><td>700/24</td><td>800/24</td></tr><tr><td rowspan="2">110</td><td>400/12</td><td>400/12</td><td>400/12</td><td>400/12</td><td>400/12</td><td>400/12</td><td>400/12</td><td>400/12</td><td>400/12</td><td>400/12</td><td>400/12</td><td>400/12</td><td>400/12</td><td>400/12</td><td>400/12</td><td>400/12</td></tr><tr><td>400/12</td><td>400/12</td><td>400/12</td><td>400/12</td><td>400/12</td><td>400/12</td><td>400/12</td><td>400/12</td><td>400/12</td><td>400/12</td><td>400/12</td><td>400/12</td><td>400/12</td><td>400/12</td><td>400/12</td><td>400/12</td></tr><tr><td colspan="2">Man with CML</td><td colspan="2">Woman with CML</td><td colspan="2"></td><td colspan="2"></td><td colspan="2"></td><td colspan="2"></td><td colspan="2"></td><td colspan="2"></td><td colspan="2"></td></tr><tr><td colspan="2">Man with GIST</td><td colspan="2">Woman with GIST</td><td colspan="2"></td><td colspan="2"></td><td colspan="2"></td><td colspan="2"></td><td colspan="2"></td><td colspan="2"></td><td colspan="2"></td></tr></tbody></table>
## 4. Discussion

The suboptimal exposure associated with dosage recommendations approved at the time of commercialization of anticancer agents is increasingly recognized and justifies efforts towards dosage optimization after drug approval [35].

Because of increasing evidence of concentration-effect relationships, achieving therapeutic concentration of tyrosine kinases inhibitors has become a relevant objective in clinical practice [36,37]. As imatinib pharmacokinetic variability is large, with a 20-fold variability in steady-state trough concentrations in patients receiving the same dosage [16], imatinib dose is poorly predictive of drug exposure. The standard initial dose of imatinib of 400 mg has not been designed to achieve a target Cmin ≥ 1000 ng/mL. Larson et al. reported steady-state Cmin 25th, 50th and 75th percentile values of 647, 879 and 1170 ng/mL in 351 CML patients who received 400 mg of imatinib [16]. In another study in 108 patients with GIST, with 80% of them receiving 400 mg/24 h, 44.4% of patients had their first imatinib Cmin at steady-state < 1000 ng/mL [38]. This means that many patients under this dosage do not achieve the target concentration in both conditions and that precision dosing of imatinib is desirable.

The use of PK models is the state-of-the-art approach for designing drug dosage regimen targeting a specific exposure [39,40]. The traditional, target concentration dosing (TCD), approach for model-based initial drug dosing consists in using a single set of PK parameters, usually based on a population PK study, to calculate the dose necessary to achieve a single target value of the exposure index of interest (e.g., trough concentration, maximal concentration or area under the curve) within the therapeutic interval. Recently, we have developed an alternative model-based dosing approach that specifically optimizes the probability to achieve a concentration or exposure interval. We have previously evaluated this target interval dosing (TID) method for busulfan dosing in children [28], and this approach outperformed other model-based approaches in terms of target attainment.

In this study, we applied the TID method for imatinib dosing and compared its performance with that of other dosing approaches in simulated patients and using real patients' data. When targeting the imatinib Cmin interval of 1000–2000 ng/mL, the TID approach was effective in maximizing the rate of target attainment. However, in this example, it did not perform better than the traditional model-based approach targeting the mid-interval (TCD). This result may not be generalizable as the ability to achieve a target interval can depend on several features including the variability of PK parameters and the width of the target interval. However, the TCD method should be effective when the condition of a relatively symmetrical distribution of the exposure around the median is fulfilled, while the TID approach may be superior in case of heavy-tailed distribution of the exposure index [28]. The latter also has the unique characteristic of accommodating single bound intervals, and so can minimize under- or overexposure. In this application to imatinib, the TCI-OD<sub>Lmin</sub> and -OD<sub>L5%</sub> methods consistently minimized the proportion of Cmin < 1000 mg/L. However, these methods were associated with the highest dosage requirements and higher rates of Cmin > 2000 ng/mL, which may raise safety concerns. This illustrates the limitations of any dosing method based on population data only and the need to perform TDM for individual control of drug exposure during treatment [27,41].

Imatinib TDM has been officially recommended by the International Association of Therapeutic Drug Monitoring and Clinical Toxicology for patients with CML and GIST in 2021 [24]. The dosages based on the TID approach presented in Table 5 could be used when no TDM is possible, or before TDM can be performed, to optimize the achievement of the Cmin target interval of 1000 to 2000 ng/mL.

Our study results also illustrate the limitations of the “one size fits all” approach for imatinib dosing. The standard 400 mg/24 h dosage was associated with very low rates of simulated Cmin ≥ 1000 ng/mL, especially in male patients, which confirms the results of several observational studies [16,21]. Based on our simulations, twice-daily dosing with 200 mg or 300 mg of imatinib every 12 h would enhance the achievement of Cmin target. However, a large variability in exposure and significant proportions of
under or overexposure would still occur because fixed-dose approaches do not consider PK variability and are not goal-oriented. Only model-based approaches can consistently reduce the variability and control the drug exposure optimally in all patients' group identified by population PK. Cost savings may also be expected by optimizing the dosage and hence the efficacy of imatinib, now that generic copies are available, while second-line treatments remain expensive [25].

This study is based on several assumptions and has a number of limitations that should be acknowledged. First, all model-based dosages (with both the TID and TCD approaches) and dosing recommendations were based on the model from Widmer et al. [26]. So, it is assumed that this model adequately describes imatinib concentration, at least Cmin. While other population PK models of imatinib have been published [42], this model has been validated for Bayesian analysis of imatinib concentration in routine TDM conditions [27]. In addition, the concentrations simulated with this model are in close agreement with real data. The 25th, 50th and 75th percentiles of Cmin values simulated in 800 patients for the 400 mg/24 h dosage were 611, 804 and 1099 ng/mL and compared well with the values from 351 patients in the study from Larson et al. presented above [16]. However, as true parameter values of both simulated and real patients were generated using the Widmer model, these conditions may favor the performance of the model-based approaches. This is especially true for the evaluation based on real patient data, as the individual Bayesian estimates of PK parameters are likely to shrink towards the typical values and so favor the TCD approach. A prospective clinical evaluation of such model-based approaches is necessary to confirm their interest in controlling imatinib Cmin. The dosages suggested based on this modeling approach (Table 5) were defined for specific age (20 to 90 years) and weight (40 to 110 kg) ranges and cannot be extrapolated out of those ranges. They were based on the TID-ODint approach that should minimize both under- and overexposure. The other TID methods (ODLmin and ODLS%) could lead to different dosages.

We did not directly consider safety in the dosing. The upper bound of the target interval was arbitrarily fixed at 2000 ng/mL, i.e., twice the lower bound value. Although rates of adverse reactions to imatinib such as fluid retention, rash and anemia were found to rise with increasing Cmin [16,23,42], there is currently no clear recommendation regarding the maximum Cmin that should not be exceeded [24]. In their meta-analysis, Garcia-Ferrer et al. discussed this upper limit of imatinib Cmin, mentioning both 3000 ng/mL and 1500 ng/mL, but recognized the lack of evidence. The upper bound used in our study was in between and appears reasonable. As the TID-ODint approach is based on targeting a Cmin interval, not only exceeding the efficacy threshold, we believe it is relevant to prevent high exposure that might be more toxic. Finally, our dosing approach does not handle the influence of other covariates such as drug-drug interaction. Cytochrome P450 3A4-mediated drug interactions may significantly alter imatinib exposure [43]. This could be considered in further development of the approach. TDM is probably the best way to deal with this risk at the individual level.

**5. Conclusions**

This study confirms that fixed-dose regimens of imatinib currently recommended are associated with low Cmin in most patients with CML or GIST. Precision dosing of imatinib based on PK models can improve the achievement of target imatinib Cmin in patients treated for CML or GIST. Our target interval dosing approach has the ability to optimize the achievement of the imatinib target Cmin interval, as well as minimizing underexposure. Dosing recommendations based on this approach are suggested. While model-based dosing can improve initial dosing, TDM remains necessary to control the individual exposure and adjust the dosage during therapy. The combination of initial model-based dosing and subsequent dose individualization based on TDM are sound basis for precision medicine with imatinib therapy and other drugs with exposure-response relationships in oncology.
**Author Contributions:** Conceptualization, S.G.; methodology, S.G., M.G., V.G., T.B., C.C. and N.W.; software, S.G.; validation, S.G.; formal analysis, S.G.; investigation, V.G., T.B. and N.W.; data curation, V.G., N.W., T.B., C.C. and M.G.; writing—original draft preparation, S.G.; writing—review and editing, all authors; supervision, S.G., N.W., C.C. and T.B.; project administration, S.G. All authors have read and agreed to the published version of the manuscript.

**Funding:** This research received no external funding. It was performed as part of routine activity that is supported by the authors' institutions.

**Institutional Review Board Statement:** The model used in this study was developed in the population pharmacokinetic study from Widmer et al. [26]. This original study was approved by ethics committee of the Lausanne Faculty of Medicine (Lausanne, Switzerland), with approval number 100/02 and informed written consent was obtained from the participants. Part of the present study was based on patient's TDM data from another previously published study [27]. This study was also approved by the ethics committee of the Lausanne Faculty of Medicine (Lausanne, Switzerland), with approval number 31/09 and informed written consent was obtained from the participants.

**Informed Consent Statement:** Informed consent was obtained from all subjects involved in the study (see above).

**Data Availability Statement:** Data are available upon request.

**Conflicts of Interest:** The authors declare no conflict of interest.

## References

1. Capdeville, R.; Buchdunger, E.; Zimmermann, J.; Matter, A. Glivec (STI571, imatinib), a rationally developed, targeted anticancer drug. *Nat. Rev. Drug Discov.* **2002**, *1*, 493–502. [CrossRef] [PubMed]

2. Demetri, G.D.; von Mehren, M.; Blanke, C.D.; Van den Abbeele, A.D.; Eisenberg, B.; Roberts, P.J.; Heinrich, M.C.; Tuveson, D.A.; Singer, S.; Janicek, M.; et al. Efficacy and safety of imatinib mesylate in advanced gastrointestinal stromal tumors. *N. Engl. J. Med.* **2002**, *347*, 472–480. [CrossRef]

3. Druker, B.J.; Guilhot, F.; O'Brien, S.G.; Gathmann, I.; Kantarjian, H.; Gattermann, N.; Deininger, M.W.; Silver, R.T.; Goldman, J.M.; Stone, R.M.; et al. Five-year follow-up of patients receiving imatinib for chronic myeloid leukemia. *N. Engl. J. Med.* **2006**, *355*, 2408–2417. [CrossRef] [PubMed]

4. Maggiora, G. Is imatinib a prototypical example of targeted drug therapy? *Future Med. Chem.* **2016**, *8*, 1907–1911. [CrossRef] [PubMed]

5. Schilsky, R.L. Personalized medicine in oncology: The future is now. *Nat. Rev. Drug Discov.* **2010**, *9*, 363–366. [CrossRef]

6. Kalia, M. Personalized oncology: Recent advances and future challenges. *Metab. Clin. Exp.* **2013**, *62* (Suppl. 1), S11–S14. [CrossRef]

7. Joyner, M.J.; Paneth, N. Seven Questions for Personalized Medicine. *JAMA* **2015**, *314*, 999–1000. [CrossRef]

8. Collins, F.S.; Varmus, H. A new initiative on precision medicine. *N. Engl. J. Med.* **2015**, *372*, 793–795. [CrossRef]

9. Darwich, A.S.; Ogungbenro, K.; Vinks, A.A.; Powell, J.R.; Reny, J.L.; Marsousi, N.; Daali, Y.; Fairman, D.; Cook, J.; Lesko, L.J.; et al. Why has model-informed precision dosing not yet become common clinical reality? lessons from the past and a roadmap for the future. *Clin. Pharmacol. Ther.* **2017**, *101*, 646–656. [CrossRef]

10. Neely, M. Scalpels not hammers: The way forward for precision drug prescription. *Clin. Pharmacol. Ther.* **2017**, *101*, 368–372. [CrossRef]

11. Peng, B.; Lloyd, P.; Schran, H. Clinical pharmacokinetics of imatinib. *Clin. Pharmacokinet.* **2005**, *44*, 879–894. [CrossRef] [PubMed]

12. Stegmeier, F.; Warmuth, M.; Sellers, W.R.; Dorsch, M. Targeted cancer therapies in the twenty-first century: Lessons from imatinib. *Clin. Pharmacol. Ther.* **2010**, *87*, 543–552. [CrossRef] [PubMed]

13. Kantarjian, H.; Talpaz, M.; O'Brien, S.; Garcia-Manero, G.; Verstovsek, S.; Giles, F.; Rios, M.B.; Shan, J.; Letvak, L.; Thomas, D.; et al. High-dose imatinib mesylate therapy in newly diagnosed Philadelphia chromosome-positive chronic phase chronic myeloid leukemia. *Blood* **2004**, *103*, 2873–2878. [CrossRef] [PubMed]

14. Cortes, J.E.; Baccarani, M.; Guilhot, F.; Druker, B.J.; Branford, S.; Kim, D.W.; Pane, F.; Pasquini, R.; Goldberg, S.L.; Kalaycio, M.; et al. Phase III, randomized, open-label study of daily imatinib mesylate 400 mg versus 800 mg in patients with newly diagnosed, previously untreated chronic myeloid leukemia in chronic phase using molecular end points: Tyrosine kinase inhibitor optimization and selectivity study. *J. Clin. Oncol. Off. J. Am. Soc. Clin. Oncol.* **2010**, *28*, 424–430. [CrossRef]

15. Hoffmann, V.S.; Hasford, J.; Deininger, M.; Cortes, J.; Baccarani, M.; Hehlmann, R. Systematic review and meta-analysis of standard-dose imatinib vs. high-dose imatinib and second generation tyrosine kinase inhibitors for chronic myeloid leukemia. *J. Cancer Res. Clin. Oncol.* **2017**, *143*, 1311–1318. [CrossRef]

16. Larson, R.A.; Druker, B.J.; Guilhot, F.; O'Brien, S.G.; Riviere, G.J.; Krahnke, T.; Gathmann, I.; Wang, Y. Imatinib pharmacokinetics and its correlation with response and safety in chronic-phase chronic myeloid leukemia: A subanalysis of the IRIS study. *Blood* **2008**, *111*, 4022–4028. [CrossRef]
17. Picard, S.; Titier, K.; Etienne, G.; Teilhet, E.; Ducint, D.; Bernard, M.A.; Lassalle, R.; Marit, G.; Reiffers, J.; Begaud, B.; et al. Trough imatinib plasma levels are associated with both cytogenetic and molecular responses to standard-dose imatinib in chronic myeloid leukemia. *Blood* **2007**, *109*, 3496-3499. [CrossRef]

18. Takahashi, N.; Wakita, H.; Miura, M.; Scott, S.A.; Nishii, K.; Masuko, M.; Sakai, M.; Maeda, Y.; Ishige, K.; Kashimura, M.; et al. Correlation between imatinib pharmacokinetics and clinical response in Japanese patients with chronic-phase chronic myeloid leukemia. *Clin. Pharmacol. Ther.* **2010**, *88*, 809-813. [CrossRef]

19. Demetri, G.D.; Wang, Y.; Wehrle, E.; Racine, A.; Nikolova, Z.; Blanke, C.D.; Joensuu, H.; von Mehren, M. Imatinib plasma levels are correlated with clinical benefit in patients with unresectable/metastatic gastrointestinal stromal tumors. *J. Clin. Oncol. Off. J. Am. Soc. Clin. Oncol.* **2009**, *27*, 3141-3147. [CrossRef]

20. Garcia-Ferrer, M.; Wojnicz, A.; Mejia, G.; Koller, D.; Zubiaur, P.; Abad-Santos, F. Utility of Therapeutic Drug Monitoring of Imatinib, Nilotinib, and Dasatinib in Chronic Myeloid Leukemia: A Systematic Review and Meta-analysis. *Clin. Ther.* **2019**, *41*, 2558-2570.e2557. [CrossRef]

21. Bouchet, S.; Titier, K.; Moore, N.; Lassalle, R.; Ambrosino, B.; Poulette, S.; Schuld, P.; Belanger, C.; Mahon, F.X.; Molimard, M. Therapeutic drug monitoring of imatinib in chronic myeloid leukemia: Experience from 1216 patients at a centralized laboratory. *Fundam. Clin. Pharmacol.* **2013**, *27*, 690-697. [CrossRef] [PubMed]

22. Bouchet, S.; Poulette, S.; Titier, K.; Moore, N.; Lassalle, R.; Abouelfath, A.; Italiano, A.; Chevreau, C.; Bompas, E.; Collard, O.; et al. Relationship between imatinib trough concentration and outcomes in the treatment of advanced gastrointestinal stromal tumours in a real-life setting. *Eur. J. Cancer* **2016**, *57*, 31-38. [CrossRef] [PubMed]

23. Gotta, V.; Bouchet, S.; Widmer, N.; Schuld, P.; Decosterd, L.A.; Buclin, T.; Mahon, F.X.; Csajka, C.; Molimard, M. Large-scale imatinib dose-concentration-effect study in CML patients under routine care conditions. *Leuk. Res.* **2014**, *38*, 764-772. [CrossRef] [PubMed]

24. Clarke, W.A.; Chatelut, E.; Fotoohi, A.K.; Larson, R.A.; Martin, J.H.; Mathijssen, R.H.J.; Salamone, S.J. Therapeutic Drug Monitoring in Oncology: International Association of Therapeutic Drug Monitoring and Clinical Toxicology Consensus Guidelines for Imatinib Therapy. *Eur. J. Cancer* **2021**, *157*, 428-440. [CrossRef] [PubMed]

25. Zuidema, S.; Desar, I.M.E.; van Erp, N.P.; Kievit, W. Optimizing the dose in patients treated with imatinib as first line treatment for gastrointestinal stromal tumours: A cost-effectiveness study. *Br. J. Clin. Pharmacol.* **2019**, *85*, 1994-2001. [CrossRef] [PubMed]

26. Widmer, N.; Decosterd, L.A.; Csajka, C.; Leyvraz, S.; Duchosal, M.A.; Rosselet, A.; Rochat, B.; Eap, C.B.; Henry, H.; Biollaz, J.; et al. Population pharmacokinetics of imatinib and the role of alpha-acid glycoprotein. *Br. J. Clin. Pharmacol.* **2006**, *62*, 97-112. [CrossRef] [PubMed]

27. Gotta, V.; Widmer, N.; Montemurro, M.; Leyvraz, S.; Haouala, A.; Decosterd, L.A.; Csajka, C.; Buclin, T. Therapeutic drug monitoring of imatinib: Bayesian and alternative methods to predict trough levels. *Clin. Pharmacokinet.* **2012**, *51*, 187-201. [CrossRef]

28. Philippe, M.; Neely, M.; Bertrand, Y.; Bleyzac, N.; Goutelle, S. A Nonparametric Method to Optimize Initial Drug Dosing and Attainment of a Target Exposure Interval: Concepts and Application to Busulfan in Pediatrics. *Clin. Pharmacokinet.* **2017**, *56*, 435-447. [CrossRef]

29. Katz, D.; D'Argenio, D.Z. Implementation and evaluation of control strategies for individualizing dosage regimens, with application to the aminoglycoside antibiotics. *J. Pharmacokinet. Biopharm.* **1986**, *14*, 523-537. [CrossRef]

30. D'Argenio, D.Z.; Rodman, J.H. Targeting the systemic exposure of teniposide in the population and the individual using a stochastic therapeutic objective. *J. Pharmacokinet. Biopharm.* **1993**, *21*, 223-251. [CrossRef]

31. Katz, D.; D'Argenio, D.Z. Discrete approximation of multivariate densities with application to Bayesian estimation. *Comput. Stat. Data Anal.* **1984**, *2*, 27-36. [CrossRef]

32. Macdonald, I.; Staatz, C.E.; Jelliffe, R.W.; Thomson, A.H. Evaluation and comparison of simple multiple model, richer data multiple model, and sequential interacting multiple model (IMM) Bayesian analyses of gentamicin and vancomycin data collected from patients undergoing cardiothoracic surgery. *Ther. Drug Monit.* **2008**, *30*, 67-74. [CrossRef] [PubMed]

33. Holford, N.H. Target concentration intervention: Beyond Y2K. *Br. J. Clin. Pharmacol.* **1999**, *48*, 9-13. [CrossRef] [PubMed]

34. Claxton, A.J.; Cramer, J.; Pierce, C. A systematic review of the associations between dose regimens and medication compliance. *Clin. Ther.* **2001**, *23*, 1296-1310. [CrossRef]

35. Overbeek, J.K.; Heine, R.T.; Verheul, H.M.W.; Chatelut, E.; Rudek, M.A.; Gurney, H.; Plummer, R.; Gilbert, D.C.; Buclin, T.; Burger, D.M.; et al. Off-label, but on target: The evidence needed to implement alternative dosing regimens of anticancer drugs. *ESMO Open* **2023**, *8*, 100749. [CrossRef]

36. Yu, H.; Steeghs, N.; Nijenhuis, C.M.; Schellens, J.H.; Beijnen, J.H.; Huitema, A.D. Practical guidelines for therapeutic drug monitoring of anticancer tyrosine kinase inhibitors: Focus on the pharmacokinetic targets. *Clin. Pharmacokinet.* **2014**, *53*, 305-325. [CrossRef]

37. Verheijen, R.B.; Yu, H.; Schellens, J.H.M.; Beijnen, J.H.; Steeghs, N.; Huitema, A.D.R. Practical Recommendations for Therapeutic Drug Monitoring of Kinase Inhibitors in Oncology. *Clin. Pharmacol. Ther.* **2017**, *102*, 765-776. [CrossRef]

38. Farag, S.; Verheijen, R.B.; Martijn Kerst, J.; Cats, A.; Huitema, A.D.; Steeghs, N. Imatinib Pharmacokinetics in a Large Observational Cohort of Gastrointestinal Stromal Tumour Patients. *Clin. Pharmacokinet.* **2017**, *56*, 287-292. [CrossRef]
39. Jelliffe, R.W.; Schumitzky, A.; Bayard, D.; Milman, M.; Van Guilder, M.; Wang, X.; Jiang, F.; Barbaut, X.; Maire, P. Model-based, goal-oriented, individualised drug therapy. Linkage of population modelling, new 'multiple model' dosage design, bayesian feedback and individualised target goals. *Clin. Pharmacokinet.* **1998**, *34*, 57–77. [CrossRef]
40. Patel, J.N.; Papachristos, A. Personalizing chemotherapy dosing using pharmacological methods. *Cancer Chemother. Pharmacol.* **2015**, *76*, 879–896. [CrossRef]
41. Buclin, T.; Thoma, Y.; Widmer, N.; Andre, P.; Guidi, M.; Csajka, C.; Decosterd, L.A. The Steps to Therapeutic Drug Monitoring: A Structured Approach Illustrated with Imatinib. *Front. Pharmacol.* **2020**, *11*, 177. [CrossRef] [PubMed]
42. Gotta, V.; Buclin, T.; Csajka, C.; Widmer, N. Systematic review of population pharmacokinetic analyses of imatinib and relationships with treatment outcomes. *Ther. Drug Monit.* **2013**, *35*, 150–167. [CrossRef] [PubMed]
43. Haouala, A.; Widmer, N.; Duchosal, M.A.; Montemurro, M.; Buclin, T.; Decosterd, L.A. Drug interactions with the tyrosine kinase inhibitors imatinib, dasatinib, and nilotinib. *Blood* **2011**, *117*, e75–e87. [CrossRef] [PubMed]

**Disclaimer/Publisher's Note:** The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.