Open Access
Open Peer Review

This article has Open Peer Review reports available.

How does Open Peer Review work?

Rapid initial decline in BCR-ABL1 is associated with superior responses to second-line nilotinib in patients with chronic-phase chronic myeloid leukemia

  • Andrew M Stein1Email author,
  • Giovanni Martinelli2,
  • Timothy P Hughes3,
  • Martin C Müller4,
  • Lan Beppu5,
  • Enrico Gottardi6,
  • Susan Branford3,
  • Simona Soverini2,
  • Richard C Woodman7,
  • Andreas Hochhaus8,
  • Dong-Wook Kim9,
  • Giuseppe Saglio6 and
  • Jerald P Radich5
BMC Cancer201313:173

DOI: 10.1186/1471-2407-13-173

Received: 23 August 2012

Accepted: 31 January 2013

Published: 2 April 2013



We evaluated BCR-ABL1 kinetics in patients treated with nilotinib and analyzed whether a dynamic model of changes in BCR-ABL1 levels over time could be used to predict long-term responses.


Patients from the nilotinib registration trial (CAMN107A2101; registered at as NCT00109707) who had imatinib-resistant or -intolerant Philadelphia chromosome–positive (Ph+) chronic myeloid leukemia (CML) in chronic phase (CP) or accelerated phase with BCR-ABL1 > 10% (on the international scale [IS]) at baseline and, in the first 6 months, had at least three BCR-ABL1 transcript measurements and an average daily dose of at least 720 mg were included in this analysis (N = 123).


More than half of patients (65/123; 53%) had a slow monophasic response and the remainder (58/123; 47%) had a biphasic response, in which patients had a rapid initial decrease in BCR-ABL1 transcripts followed by a more gradual response. The biphasic response type strongly correlated with improved event-free survival (EFS). Data in the first 6 months of follow-up were sufficient to predict EFS at 24 months.


Unlike newly diagnosed patients with Ph+ CML-CP—in whom the majority had a biphasic response—approximately half of patients with imatinib-resistant or -intolerant CML had a slower, monophasic response. Second-line patients who did have a biphasic response had an EFS outlook similar to that of newly diagnosed patients treated with imatinib. Our model was comparable to using BCR-ABL1 (IS) ≤ 10% at 6 months as a threshold for predicting EFS.


Chronic myeloid leukemia Nilotinib Mathematical modeling BCR-ABL Molecular response


Chronic myeloid leukemia (CML) is a clonal myeloproliferative disorder characterized by the expansion of hematopoietic cells carrying the oncogenic BCR-ABL1 fusion gene. The BCR-ABL1 fusion gene is formed by a reciprocal translocation of a fragment of the Abelson gene (ABL1), situated at the breakpoint on chromosome 9, with a fragment of the breakpoint cluster region (BCR) gene on chromosome 22. This translocation forms the Philadelphia chromosome (Ph), which encodes the constitutively active BCR-ABL1 protein tyrosine kinase [1, 2]. Deregulated phosphorylation of protein tyrosine residues by this constitutively active tyrosine kinase leads to activation of downstream effectors that enable growth factor–independent proliferation and neoplastic transformation in the transformed hematopoietic cells [3].

Imatinib is a BCR-ABL1 tyrosine kinase inhibitor (TKI) that demonstrated superior response rates and improved tolerability over interferon-α plus cytarabine in the phase III open-label International Randomized Study of Interferon and STI571 (IRIS) trial [4]. Despite the success of imatinib, a significant proportion of patients discontinue therapy (34% in the IRIS trial) [5] for reasons that include intolerance or resistance to the drug.

Nilotinib was rationally designed after the discovery of resistance to imatinib in early clinical trials [68] and laboratory evidence showing that point mutations in the BCR-ABL1 kinase domain mutation were the most common cause of such resistance [9]. A more potent BCR-ABL1 kinase inhibitor was hypothesized to reduce the reservoir of leukemic cells in patients, thereby impeding the emergence of drug resistance. Nilotinib is approved in more than 60 countries worldwide for the treatment of newly diagnosed patients with Ph+ CML in the chronic phase (CP) and in patients with Ph+ CML-CP and in the accelerated phase (AP) who have failed prior therapy, including imatinib. Approval as second-line treatment was based on results from a phase II open-label registration study [10] that showed durable responses and overall survival of 87% with 24 months of follow-up [11].

Normalization of blood counts, as measured by hematologic responses; reduction and elimination of the Ph chromosome, as measured by cytogenetic monitoring; and reduction and elimination of BCR-ABL1 gene expression, as measured by real-time quantitative polymerase chain reaction (RQ-PCR) are important measures of treatment success. Before the routine use of BCR-ABL1 TKIs, the evaluation of hematologic and cytogenetic responses was sufficient to determine treatment efficacy, as improvement in these outcomes was the limit of response in most patients. However, with more potent BCR-ABL1–targeting therapies, deeper responses are achieved in most patients, necessitating more sensitive methods of disease detection.

An abundance of literature explores the use of molecular response to predict long-term outcome [1218]. A common approach is the landmark analysis, whereby a particular threshold (e.g., BCR-ABL1 ratio ≤ 0.1% on the international scale [IS]) and time (e.g., 3 months) are chosen, and longer-term outcomes (e.g., progression-free survival or event-free survival [EFS] at 2 years) are evaluated according to whether patients attained the threshold at the landmark time point [1823]. While this approach is reasonable in controlled clinical trials, challenges may arise in practice, including missing samples or those inconsistent with prior or subsequent samples. For example, a patient may not have a BCR-ABL1 assessment at month 6 but may have assessments at months 4 and 8. Or, a patient may achieve major molecular response (MMR) at month 6, lose MMR at month 12, and then regain MMR at month 18. Rather than assessing a patient based on his or her response at a single time point, using all of the information from a patient’s full clinical history may provide a more robust method for classifying patient response.

A mathematical model is a natural tool for integrating information from the full BCR-ABL1 transcript dynamic time course. Several groups have established such models for predicting BCR-ABL1 response in patients with CML [13, 15, 18]. We previously presented an approach using standard pharmacokinetic/pharmacodynamic methods to model the clinical data of newly diagnosed patients treated with imatinib in the IRIS trial [18]. Here, we applied this model to patients with CML who received imatinib treatment that had failed because of resistance or intolerance and who were subsequently treated with nilotinib in the registration study. The objectives of this analysis were to (1) determine whether BCR-ABL1 kinetics in patients treated with nilotinib were similar to those in patients treated with imatinib; (2) explore whether a model using BCR-ABL1 as a continuous measurement could potentially be used to predict long-term response based on short-term data; and (3) compare such methods with simpler methods.


Patients and samples

The nilotinib registration trial (CAMN107A2101; registered at as NCT00109707) was an open-label, multicenter, single-arm phase II study of nilotinib (400 mg twice daily) in patients with imatinib-resistant or -intolerant Ph+ CML-CP (N = 321; 226 imatinib-resistant, 95 imatinib-intolerant) with a median exposure to nilotinib of 561 days [10, 11]. This study was conducted in accordance with the Declaration of Helsinki, and the protocol was reviewed by the ethics committee or institutional review board at each participating institution. All patients gave written informed consent.

RQ-PCR measurements of BCR-ABL1 were taken monthly during the first 3 months of the trial and once every 3 months thereafter, as described previously [24]. For this analysis, patients were required to have at least three PCR measurements, an average daily dose of at least 720 mg (90% of the target dose) in the first 6 months of the study, and BCR-ABL1 (IS) > 10% at baseline (n = 123).

BCR-ABL1 (IS) ≤ 10% can be measured effectively [25]. The upper limit of 10% is determined by two sources of nonlinearity in the BCR-ABL1 ratio estimation: (1) depending on the primer design for the ABL1 control gene, both ABL1 and BCR-ABL1 are amplified when the BCR-ABL1 levels are high [26]; and (2) healthy cells have two copies of the BCR and ABL1 control genes, while leukemic cells have one control gene and one BCR-ABL1 gene. There is no simple formula that can be applied to remove this nonlinearity because expression rates of the control genes may be different at the time of diagnosis and at the time that a patient achieves a good response. Thus, while measurements above 10% are included in this analysis so that we can describe the full range of BCR-ABL1 dynamics for all patients, quantitative values of the parameters governing the dynamics at the start of treatment (A, α, μ; see model below) should be considered qualitatively.

Model selection and parameter fitting

The time course of BCR-ABL1 transcript levels was modeled using nonlinear mixed effects (NLME) methodology. This is an accepted approach for deriving population models from clinical data [27, 28]. We used maximum likelihood methods, as implemented in Monolix ( NLME software in MATLAB (The MathWorks, Inc., Natick, MA), to fit the following between-patient mixture model to the log10 of the BCR-ABL1 ratios measured at time t, R(t), where treatment begins at t = 0:
R t = { l o g 10 A e μt + ε Slow BCR - ABL respons e with probability p l o g 10 A e αt + B e βt + ε Fast BCR - ABL response with probability 1 - p

The model represents two typical profiles: (1) slow responders with a monophasic decline, who typically had little if any decline in BCR-ABL1 (IS)% (at rate μ), and (2) fast responders with a biphasic decline, in whom a rapid initial decrease in the BCR-ABL1 (IS)% (at rate α) was followed by a more gradual response (at rate β) that could be either increasing or decreasing. A similar version of the above model was used to fit patients from the IRIS study [18], although we excluded the triphasic mixture population due to limited follow-up time (2 years) in the present study. The constant terms (A, B) were restricted to be greater than zero because BCR-ABL1 (IS)% is always positive. The rate constants (μ, α, β) have units of 1/year. Before treatment starts, the log10 of the BCR-ABL1 (IS)% is assumed to be log(A) or log(A + B) with additive measurement error. The model has an additive residual error term (ε) in log space that describes variability in the measurement due to assay and interoccasion variability and model error. All mono- or biexponential parameters (A, B, α, β, μ) were treated as mixed effects with both a fixed component that represents the typical value for the population and a random component that represents the interpatient variability. The mixture probability p was a fixed effect. The constant terms (A, B) were chosen from a log-normal distribution because BCR-ABL1 (IS)% is always greater than zero (e.g., A = θΑ*exp(ηΑ), where ηΑ is normally distributed with mean zero and variance ωA). Similarly, both α and μ are always less than or equal to zero and so are chosen from a log-normal distribution with a negative fixed-effect term (e.g., θα < 0). Thus, we have α = θα*exp(ηα). Finally, we chose β from a normal distribution, allowing it to be either negative or positive (i.e., β = θβ + ηβ).

We used 0.0032% (a 4.5-log10 decrease in the BCR-ABL1 [IS]%) as the minimal level of detection of the PCR assay for all laboratories and treated these measurements as left-censored data. BCR-ABL1 transcript levels below 0.0032% were established as “undetectable” in the IRIS trial [29], and BCR-ABL (IS) ≤ 0.0032% has subsequently been defined as molecular remission at 4.5 logs (MR4.5) [30, 31]. Rather than remove all measurements below the limit of quantitation or fix them at a certain value, we used the M3 maximum-likelihood approach [32], which treated any measurement falling below the limit of detection as having an upper boundary of 0.0032% to make optimal use of all measurements.

The maximum-likelihood approach yields information about both the population and each individual patient. For the population, the typical model parameters (A, B, μ, α, β, p) and the interpatient variability (of A, B, μ, α, β) were estimated. In addition, for each individual, the most likely set of parameters that describe that individual was estimated, along with the patient’s BCR-ABL1 response classification (slow vs. fast). In particular, each patient was categorized as a fast or slow responder based on the best description of that patient as calculated by the maximum-likelihood method. For the biphasic patients, we also computed the transition time at which the patient transitioned from the α to the β phase of the decline. We defined the transition time as the point of maximum curvature of the fitted BCR-ABL1 transcript profile, where curvature was defined in the differential geometrical sense as κ(t) = |R”(t)|/(1 + R’(t) 2 ) 3/2 .

Response prediction and model comparison

We then tested the ability of the model to predict the event status of individual patients at year 2 from data in the first 1, 3, 6, 9, 12, and 18 months. EFS was defined as the time between randomization and any of the following events on treatment: death due to any cause, progression to AP or blast crisis, loss of complete hematologic response, or loss of major cytogenetic response. For each landmark time, patients who progressed or died before that time were excluded from the analysis.

Categorization was performed by fixing the average model parameters for the population (A, B, μ, α, β, p) and then using the data from baseline to the landmark time of interest to estimate the individual parameters and response classification (fast vs. slow responder) for each patient. When fewer than 6 months of data were used, the long-term response parameters B and β were not estimable, in which case the model naturally assumed the population estimates for each patient. Because the purpose of this analysis was to categorize each patient as a fast or slow responder, the estimation of the long-term response was not essential.

Finally, we compared the models’ ability to predict event status at 2 years based on a single BCR-ABL1 transcript measurement at 3 or 6 months, respectively. We predicted each patient’s event status by whether the BCR-ABL1 levels were above or below a particular threshold at a particular landmark time. The thresholds used were BCR-ABL1 (IS) > 10%, > 1%, and > 0.1%, and the landmark times were 3 and 6 months. The models were compared according to the percentage of patients who were correctly classified, the positive predictive value (PPV; proportion of nonprogressing patients correctly classified as nonprogressing), and the negative predictive value (NPV; proportion of progressing patients correctly classified as progressing). In this case, we restricted the analysis to the 76 patients who did not progress by month 6 and whose event status by year 2 was known.


Of the total 123 patients analyzed, 65 (53%) were categorized as slow monophasic responders and 58 (47%) as fast biphasic responders (Figures 1A and 1B). Individual fits to 16 representative patients are shown in Figure 2. Note that some patients (patients 14 and 15) obtained an undetectable level of transcripts before transitioning to the β phase. The model accounts for these patients by treating them as biphasic responders with a β phase below the limit of quantification of the assay, although the true dynamics below 0.0032% BCR-ABL1 are unknown.
Figure 1

BCR-ABL1 (international scale [IS])% for all 123 patients in this analysis, divided into two groups, with an example patient shown in black within each group: (A) slow monophasic responders (n = 65) with μ = −0.4/year for the example patient and (B) fast biphasic responders (n = 58) with α = −15.1/year and β = −0.9/year for the example patient. The μ*, α*, β* parameters are shown for the example patient, where μ* = μlog10e, α* = αlog10e, and β* = βlog10e. The asterisks (*) indicate when a BCR-ABL1 (IS)% measurement fell below 0.0032% (= 100 × 10-4.5), the approximate limit of quantitation of the assay for all laboratories in the International Randomized Study of Interferon and STI571 trial; (C) histogram of the initial reduction rates μ and α. Note that there is no overlap between the two populations and that the vertical dashed line at −3/year demonstrates where the separation occurs.
Figure 2

Individual fits to 16 patients. Patients 1 to 8 are slow monophasic, and patients 9 to 16 are fast biphasic. The asterisk (*) represents measurements below the limit of quantification of the assay.

The parameter p was estimated accordingly to be 51% ± 10%, differing significantly from IRIS, where the probability of a monophasic response was only 14%. However, the other model parameters were all comparable, with typical values of: θA,IRIS = 34%, θA,2101 = 45%; θB,IRIS = 0.09%, θB,2101 = 0.05%; θμ,IRIS = −0.1/year, θμ,2101 = −0.2/year; θα,IRIS = −14/year, θα,2101 = −17/year; θβ,IRIS = −0.5/year, θβ,2101 = −0.6/year [18]. In addition, the variances on the parameters from 2101 were ωA = 0.4, ωB = 8.8, ωμ = 1.4, ωα = 0.3, ωβ = 0.3. The large variability in B indicates that the β-phase transition for the biphasic patients can occur over a large range of BCR-ABL1 values, from below 0.0032% (undetectable) to 7%. The initial response rate—which was equal to μ for the slow responders and α for the fast responders—had a bimodal distribution, with a threshold of approximately −2 to −6/year separating the two distributions (Figure 1C).

A fast biphasic response, according to the model, was associated with the best hematologic, cytogenetic, and molecular responses by 24 months, and patients classified as fast responders were more likely to be considered optimal responders according to European LeukemiaNet criteria [33] (Table 1). Furthermore, slow monophasic responders were more likely to be imatinib-resistant (as opposed to imatinib-intolerant) and have nilotinib-insensitive mutations (half maximal inhibitory concentration [IC50] > 150 nm; Y253H, E255V/K, F359V/C) [24] in BCR-ABL1 at baseline (Table 2).
Table 1

Patient response by model response type


Best response achieved over 24 months, n (%)

ELN response at 12 months, n (%)

Response type

All patients, n






Warning or failure

Not assessable

Slow (monophasic)


47 (72.3)

5 (7.7)

1 (1.5)

2 (3.1)

2 (3.1)

39 (60.0)

22 (33.8)

Fast (biphasic)


57 (98.3)

51 (87.9)

41 (70.7)

40 (69.0)

7 (12.1)

4 (6.9)

7 (12.1)

Abbreviations: CCyR = complete cytogenetic response; CHR = complete hematologic response; ELN = European LeukemiaNet; MMR = major molecular response.

Table 2

Patient baseline characteristics by model response type

Response type

All patients, n

Sensitive baseline mutation,*n (%)

Insensitive baseline mutation,n (%)

Imatinib-resistant, n (%)

Imatinib-intolerant, n (%)

Slow (monophasic)


21 (32.3)

18 (27.7)

58 (89.2)

7 (10.8)

Fast (biphasic)


18 (31.0)

2 (3.4)

42 (72.4)

16 (27.6)

* Half maximal inhibitory concentration (IC50) ≤ 150 nm.

IC50 > 150 nm; Y253H, E255V/K, F359V/C.

A small number of patients achieved a fast response without a detectable β phase. These patients had small estimated B values, indicating that the transition from the α to the β phase occurred below MR4.5. The true BCR-ABL1 dynamics below MR4.5 cannot be accurately determined at present.

The EFS outcomes of the two responder groups are shown in Figure 3, overlaid with the EFS for newly diagnosed patients with CML in the IRIS trial. Fast responders on second-line nilotinib had higher EFS rates than did slow responders and had a comparable outlook to that of a typical patient on frontline imatinib.
Figure 3

Event-free survival (EFS) according to patient group. The Kaplan-Meier plot of the imatinib arm of the International Randomized Study of Interferon and STI571 (IRIS) trial for frontline patients with chronic myeloid leukemia is also shown, to illustrate that biphasic responders do approximately as well as patients receiving frontline imatinib treatment.

Assuming the same average population parameters with different landmark times, we assessed the ability of the model to predict EFS at 2 years. At a landmark time of 6 months, the model showed a statistically significant difference in EFS between patients classified as slow (monophasic) and fast (biphasic) responders (Figure 4A). We changed the follow-up time in the model from 1 month to 18 months and found that the classification improved up to month 6, but after that there was no improvement in the percentage of patients correctly classified as having progressed (Figure 4B). The median time to transition from the α to the β phase for biphasic patients was 6.1 months, with 66% of these patients having a transition point between months 4 and 8 (Figures 4C and 4D). The mean BCR-ABL1 ratio at the time of the transition was 0.035%.
Figure 4

(A) Event-free survival (EFS) at 2 years based on patient classification using only the first 6 months of data. The difference between the EFS curves was found to be statistically significant at p = 0.0001 by a log-rank test; (B) the percentage of patients correctly classified as having an event at year 2 as a function of the amount of data used to classify each patient; (C) example patient: the ✩ indicates the transition point, defined as the point of maximum curvature where fast responders transition from the α to the β phase; the * corresponds to measurements below the limit of quantification of the polymerase chain reaction assay, assumed to be 0.0032% for all labs; (D) histogram of transition times between the α and β phases. Abbreviation: IS = international scale.

Finally, we compared model classification to other predictive measures of patient response, in particular, the molecular response at 3 and 6 months (Table 3). We found that the model classification had the highest percentage of correctly classified patients as well as the highest NPV. However, the improvement shown in the model-based classification over a simpler classification method—one categorizing patients based on whether their BCR-ABL1 (IS) was above or below 10% at 6 months—was not substantial (2.6%). We also found that the model using 6 months of follow-up was more accurate at predicting EFS at 2 years than the model using 3 months of follow-up.
Table 3

Comparison of different categorization methods for predicting event-free survival at 2 years in patients alive and event free at 6 months (n=76)

Follow-up time for subpopulation predicted to progress

N (above, below) threshold or (slow, fast) response

Correctly classified,%

Positive predictive value,%

Negative predictive value,%

3 months


  > 10% BCR-ABL

32, 44




  > 1% BCR-ABL

51, 25




  > 0.1% BCR-ABL (MMR)

63, 13


92.3 *


6 months


  > 10% BCR-ABL

26, 50




  > 1% BCR-ABL

38, 38




  > 0.1% BCR-ABL (MMR)

50, 26




 Slow (monophasic) responders

27, 49

77.6 *


61.5 *

Abbreviation: MMR = major molecular response.

* The number in bold indicates the highest score for that column.


A model-based analysis was previously developed to classify patients from the imatinib arm of the IRIS trial into three main categories: (1) slow monophasic patients who either do not respond to therapy or exhibit a very gradual BCR-ABL1 decline; (2) fast biphasic patients with an initial rapid decline in BCR-ABL1 during the first 6 months of therapy, followed by a more gradual decline; and (3) fast triphasic patients who follow the biphasic trend but relapse 2 to 4 years into treatment, exhibiting a rapid rise in BCR-ABL1 [18]. Because second-line patients treated with nilotinib in the registration trial have only been followed for 2 to 3 years thus far, the triphasic patients have been difficult to detect. Thus, we applied a scaled-down version of the model that categorized patients as either slow or fast BCR-ABL1 responders.

Patients were nearly evenly divided between the slow monophasic population (65/123; 53%) and the biphasic population (58/123; 47%). This was in contrast to the IRIS trial [18], where the vast majority of patients exhibited a biphasic response. However, the biphasic patients in the nilotinib registration trial had a similar EFS outlook to that of newly diagnosed biphasic patients in the imatinib arm of the IRIS trial, and the model parameters describing the monophasic and biphasic parameters were similar for IRIS and the nilotinib registration trial. This could be because the resistance to imatinib was overwhelmed by the more potent BCR-ABL inhibition by nilotinib and/or because "resistance" was actually intolerance or nonadherence, whereby patients did not take imatinib daily. Thus, some of the nilotinib biphasic responders may have actually been patients intolerant to or nonadherent with imatinib who might have been responders had they remained adherent to their medication.

The biphasic response type strongly correlated with improved EFS, as in IRIS. In addition, for this patient subpopulation, a 6-month follow-up time was optimal for predicting EFS at 24 months. This coincides with the observation that the transition time from the α to the β phase in biphasic patients typically occurs between months 4 and 8. Additional measurements at later times were of limited value for prediction based on model categorization. The median transition time at 6 months, at a median transition depth of 0.035% BCR-ABL1 (IS), also provides insight into why using a single measurement at 6 months outperforms using a measurement at 3 months. This level of response—nearly a 4-log reduction in BCR-ABL1 from standardized baseline—is a deeper response than MMR (BCR-ABL1 [IS] ≤ 0.1%), and achieving such a response has been associated with excellent long-term outcomes [22, 23, 34].

These findings are consistent with other studies showing that early molecular responses were associated with superior EFS and overall survival [22, 3537]. These landmark analyses, most in frontline patients treated with imatinib or imatinib-based combinations, typically focused on responses at 3 months as a predictor of later responses. Recent data from patients who received second-generation TKIs as second-line treatment at the Hammersmith Hospital in London revealed that molecular responses at 3 months were predictive of overall survival [38]. Another landmark analysis of imatinib-resistant or -intolerant patients enrolled in the nilotinib registration trial showed that molecular responses at 3 months predicted EFS at 24 months [21], but did not compare the utility of 3- versus 6-month BCR-ABL levels. A recent evaluation of frontline dasatinib suggested that the predictive value of the 3-month landmark could be improved with drug-specific BCR-ABL cutoffs [39].

We note that it may seem counterintuitive to use this model to categorize a patient as a fast biphasic responder when only 3 to 6 months of data is available, and the β phase is not yet observable. In using the biphasic model even in the presence of limited information, we are making use of prior information from both IRIS as well as this study, where almost all fast responders ultimately have a β phase. Furthermore, in some cases, the β phase does occur before month 6 (Figure 4D), so accounting for this second slope may be important for these patients.

Because the definition of the transition point as the point of maximum curvature was somewhat complex, we initially tested simpler definitions for the transition point, such as the point of maximum second derivative (R″(t)), or the point where Ae αt = Be βt . We found that the point of maximum curvature definition most consistently agreed with our estimation of where the transition point should lie based on clinical experience.

The model used in this study was not superior at predicting the longer-term event status of individual patients versus using BCR-ABL1 (IS) above or below 10% at 6 months. Furthermore, the NPV of all the classifiers in the model was relatively low (maximum of 62%). This may be due in part to 2 years of follow-up being too short and that, over longer times, the monophasic patients may continue to be more likely to progress than the biphasic patients. The NPVs and PPVs depend on the prevalence of second-line patients who do not respond well to nilotinib therapy; this prevalence may be lower in a controlled trial than in the general population. Nevertheless, PPV and NPV are preliminary tools for comparing the predictive power of different methodologies. Our modeling approach represents one of many classification schemes that can be used for dividing a population into two groups. In the future, we could consider exploring other classification approaches, such as neural networks or support vector machines.

Key benefits of the modeling approach (classification of a second-line patient as either a fast or slow responder) are that (1) the model helps to identify 6 months as an ideal follow-up time; and (2) the model provides a natural means for incorporating all measurements and treating the BCR-ABL1 dynamics as a continuous variable up to the time at which classification takes place. In the clinic, samples may be missed or multiple conflicting measurements may be obtained in the same time period. In these instances, a simple algorithm may not be able to adequately determine patient response based on one data point at a particular time. The model, however, provides a natural technique for integrating missing data and multiple samples at the same time, and produces a dichotomous classification associated with EFS. While we required patients to have at least three data points in the first 6 months, the NLME approach used here can work with fewer data points such that a patient with only one assessment at month 6 could be classified as a fast responder if he or she has achieved BCR-ABL1 (IS) ≤ 0.1% (at least MMR) or a slow responder if his or her BCR-ABL1 levels are above 10%. In the future, it will be useful to investigate how the model performance changes depending on the amount of data collected.

A related modeling approach that also treated BCR-ABL1 response as a continuous parameter was recently reported [40]. In this model, BCR-ABL1 doubling time was used to characterize relapse rate. Patients entering blast crisis had similar BCR-ABL1 doubling times to patients discontinuing therapy, whereas patients with BCR-ABL–resistant mutations who remained in CP had a longer time to relapse. These differences in the velocity of leukemic cell growth were not always evident by the BCR-ABL1 fold rise, which was dependent on the time between measurements.

The nilotinib registration trial demonstrated high rates of major cytogenetic response (59%) and complete cytogenetic response (CCyR; 44%) [11]. Most patients who achieved CCyR also achieved MMR (56%), and cytogenetic responses were durable, with 84% of patients who achieved CCyR maintaining response at 24 months. The overall survival at 24 months was 87%. Here, we showed that the EFS rates of fast responders to nilotinib were comparable to those observed with newly diagnosed patients treated with imatinib, and that BCR-ABL1 dynamics in the first 6 months of therapy were useful in predicting such responses.

This work represents an exploratory step in using a model to measure response based on BCR-ABL1 as a continuous variable and transform it into a dichotomous classification associated with EFS. Because the patient subpopulation in this study was highly heterogeneous and the subpopulation used to compare different predictors was small (n = 76), more testing and validation will be necessary before this approach can be directly applied in the clinical setting.


Landmark analyses, in which longer-term outcomes are analyzed according to achievement of molecular milestones at early time points (3 and 6 months), have become an increasingly common means of evaluating responses to both front- and second-line TKI therapies used for the treatment of CML. However, this emphasis on achieving specific molecular targets may be challenging in routine clinical practice, where samples are frequently missing. Here, we described a mathematical model that utilizes the full BCR-ABL1 transcript profile to predict longer-term responses among patients treated with nilotinib after imatinib failure. Evaluating BCR-ABL1 transcript dynamics over time in an individual patient may provide a more robust and practical method for classifying patient response.



Chronic myeloid leukemia


Abelson gene


Breakpoint cluster region


Philadelphia chromosome


International Randomized Study of Interferon and STI571


Chronic phase


Accelerated phase


Real-time quantitative polymerase chain reaction


International scale


Event-free survival


Major molecular response


Nonlinear mixed effects


Molecular remission at 4.5 logs


Positive predictive value


Negative predictive value


Half maximal inhibitory concentration


Complete cytogenetic response.



We thank Dean Bottino, Alison Carter, Anny Chen, Yen Lin Chia, Ovidiu Chiparus, Ariful Haque, Gabriel Helmlinger, Mariya Mazlina, Jerry Nedelman, John Reynolds, Yaping Shou, and Tomasz Szczudlo for their contributions to this manuscript. Financial support for medical editorial assistance was provided by Novartis Pharmaceuticals Corporation. We thank Erinn Goldman for medical editorial assistance with this manuscript.

Authors’ Affiliations

Novartis Institutes for BioMedical Research, Inc.
Department of Hematology and Oncological Sciences “L. e A. Seràgnoli, University of Bologna
Centre for Cancer Biology/SA Pathology, Royal Adelaide Hospital
Medizinische Fakultät Mannheim der Universität Heidelberg
Clinical Research Division, Fred Hutchinson Cancer Research Center
University of Turin
Novartis Pharmaceuticals Corporation
Universitätsklinikum Jena
Seoul St. Mary’s Hospital, The Catholic University of Korea


  1. Faderl S, Talpaz M, Estrov Z, O’Brien S, Kurzrock R, Kantarjian HM: The biology of chronic myeloid leukemia. N Engl J Med. 1999, 341: 164-172. 10.1056/NEJM199907153410306.View ArticlePubMedGoogle Scholar
  2. Sawyers CL: Chronic myeloid leukemia. N Engl J Med. 1999, 340: 1330-1340. 10.1056/NEJM199904293401706.View ArticlePubMedGoogle Scholar
  3. Kurzrock R, Talpaz M: The molecular pathology of chronic myelogenous leukaemia. Br J Haematol. 1991, 79 (Suppl 1): 34-37.View ArticlePubMedGoogle Scholar
  4. O’Brien SG, Guilhot F, Larson RA, Gathmann I, Baccarani M, Cervantes F, Cornelissen JJ, Fischer T, Hochhaus A, Hughes T, Lechner K, Nielsen JL, Rousselot P, Reiffers J, Saglio G, Shepherd J, Simonsson B, Gratwohl A, Goldman JM, Kantarjian H, Taylor K, Verhoef G, Bolton AE, Capdeville R, Druker BJ: Imatinib compared with interferon and low-dose cytarabine for newly diagnosed chronic-phase chronic myeloid leukemia. N Engl J Med. 2003, 348: 994-1004. 10.1056/NEJMoa022457.View ArticlePubMedGoogle Scholar
  5. Hochhaus A, O’Brien SG, Guilhot F, Druker BJ, Branford S, Foroni L, Goldman JM, Muller MC, Radich JP, Rudoltz M, Mone M, Gathmann I, Hughes TP, Larson RA: Six-year follow-up of patients receiving imatinib for the first-line treatment of chronic myeloid leukemia. Leukemia. 2009, 23: 1054-1061. 10.1038/leu.2009.38.View ArticlePubMedGoogle Scholar
  6. Ottmann OG, Druker BJ, Sawyers CL, Goldman JM, Reiffers J, Silver RT, Tura S, Fischer T, Deininger MW, Schiffer CA, Baccarani M, Gratwohl A, Hochhaus A, Hoelzer D, Fernandes-Reese S, Gathmann I, Capdeville R, O’Brien SG: A phase 2 study of imatinib in patients with relapsed or refractory Philadelphia chromosome-positive acute lymphoid leukemias. Blood. 2002, 100: 1965-1971. 10.1182/blood-2001-12-0181.View ArticlePubMedGoogle Scholar
  7. Sawyers CL, Hochhaus A, Feldman E, Goldman JM, Miller CB, Ottmann OG, Schiffer CA, Talpaz M, Guilhot F, Deininger MW, Fischer T, O’Brien SG, Stone RM, Gambacorti-Passerini CB, Russell NH, Reiffers JJ, Shea TC, Chapuis B, Coutre S, Tura S, Morra E, Larson RA, Saven A, Peschel C, Gratwohl A, Mandelli F, Ben-Am M, Gathmann I, Capdeville R, Paquette RL, Druker BJ: Imatinib induces hematologic and cytogenetic responses in patients with chronic myelogenous leukemia in myeloid blast crisis: results of a phase II study. Blood. 2002, 99: 3530-3539. 10.1182/blood.V99.10.3530.View ArticlePubMedGoogle Scholar
  8. Talpaz M, Silver RT, Druker BJ, Goldman JM, Gambacorti-Passerini C, Guilhot F, Schiffer CA, Fischer T, Deininger MW, Lennard AL, Hochhaus A, Ottmann OG, Gratwohl A, Baccarani M, Stone R, Tura S, Mahon FX, Fernandes-Reese S, Gathmann I, Capdeville R, Kantarjian HM, Sawyers CL: Imatinib induces durable hematologic and cytogenetic responses in patients with accelerated phase chronic myeloid leukemia: results of a phase 2 study. Blood. 2002, 99: 1928-1937. 10.1182/blood.V99.6.1928.View ArticlePubMedGoogle Scholar
  9. Gorre ME, Mohammed M, Ellwood K, Hsu N, Paquette R, Rao PN, Sawyers CL: Clinical resistance to STI-571 cancer therapy caused by BCR-ABL gene mutation or amplification. Science. 2001, 293: 876-880. 10.1126/science.1062538.View ArticlePubMedGoogle Scholar
  10. Kantarjian HM, Giles F, Gattermann N, Bhalla K, Alimena G, Palandri F, Ossenkoppele GJ, Nicolini FE, O’Brien SG, Litzow M, Bhatia R, Cervantes F, Haque A, Shou Y, Resta DJ, Weitzman A, Hochhaus A, le Coutre P: Nilotinib (formerly AMN107), a highly selective BCR-ABL tyrosine kinase inhibitor, is effective in patients with Philadelphia chromosome-positive chronic myelogenous leukemia in chronic phase following imatinib resistance and intolerance. Blood. 2007, 110: 3540-3546. 10.1182/blood-2007-03-080689.View ArticlePubMedGoogle Scholar
  11. Kantarjian HM, Giles FJ, Bhalla KN, Pinilla-Ibarz JA, Larson RA, Gattermann N, Ottmann OG, Hochhaus A, Radich JP, Saglio G, Hughes TP, Martinelli G, Kim DW, Shou Y, Gallagher NJ, Blakesley R, Baccarani M, Cortes J, le Coutre PD: Nilotinib is effective in patients with chronic myeloid leukemia in chronic phase following imatinib resistance or intolerance: 24-month follow-up results. Blood. 2011, 117: 1141-1145. 10.1182/blood-2010-03-277152.View ArticlePubMedGoogle Scholar
  12. Dingli D, Michor F: Successful therapy must eradicate cancer stem cells. Stem Cells. 2006, 24: 2603-2610. 10.1634/stemcells.2006-0136.View ArticlePubMedGoogle Scholar
  13. Michor F, Hughes TP, Iwasa Y, Branford S, Shah NP, Sawyers CL, Nowak MA: Dynamics of chronic myeloid leukaemia. Nature. 2005, 435: 1267-1270. 10.1038/nature03669.View ArticlePubMedGoogle Scholar
  14. Komarova NL, Wodarz D: Effect of cellular quiescence on the success of targeted CML therapy. PLoS One. 2007, 2: e990-10.1371/journal.pone.0000990.View ArticlePubMedPubMed CentralGoogle Scholar
  15. Roeder I, Horn M, Glauche I, Hochhaus A, Mueller MC, Loeffler M: Dynamic modeling of imatinib-treated chronic myeloid leukemia: functional insights and clinical implications. Nat Med. 2006, 12: 1181-1184. 10.1038/nm1487.View ArticlePubMedGoogle Scholar
  16. Wodarz D: Heterogeneity in chronic myeloid leukaemia dynamics during imatinib treatment: role of immune responses. Proc Biol Sci. 2010, 277: 1875-1880. 10.1098/rspb.2009.2179.View ArticlePubMedPubMed CentralGoogle Scholar
  17. Lenaerts T, Pacheco JM, Traulsen A, Dingli D: Tyrosine kinase inhibitor therapy can cure chronic myeloid leukemia without hitting leukemic stem cells. Haematologica. 2010, 95: 900-907. 10.3324/haematol.2009.015271.View ArticlePubMedGoogle Scholar
  18. Stein AM, Bottino D, Modur V, Branford S, Kaeda J, Goldman JM, Hughes TP, Radich JP, Hochhaus A: BCR-ABL transcript dynamics support the hypothesis that leukemic stem cells are reduced during imatinib treatment. Clin Cancer Res. 2011, 17: 6812-6821. 10.1158/1078-0432.CCR-11-0396.View ArticlePubMedGoogle Scholar
  19. Druker BJ, Guilhot F, O’Brien SG, Gathmann I, Kantarjian H, Gattermann N, Deininger MW, Silver RT, Goldman JM, Stone RM, Cervantes F, Hochhaus A, Powell BL, Gabrilove JL, Rousselot P, Reiffers J, Cornelissen JJ, Hughes T, Agis H, Fischer T, Verhoef G, Shepherd J, Saglio G, Gratwohl A, Nielsen JL, Radich JP, Simonsson B, Taylor K, Baccarani M, So C: Five-year follow-up of patients receiving imatinib for chronic myeloid leukemia. N Engl J Med. 2006, 355: 2408-2417. 10.1056/NEJMoa062867.View ArticlePubMedGoogle Scholar
  20. Branford S, Martinelli G, Saglio G, Kim D, Shou Y, Reynolds J, Woodman RC, Kantarjian H, Hochhaus A, Radich JP: Association of early molecular response to nilotinib with probability of cytogenetic response in chronic myeloid leukemia patients (pts) who fail imatinib [abstract]. J Clin Oncol. 2010, 28: 490s-10.1200/JCO.2010.29.6293.View ArticleGoogle Scholar
  21. Branford S, Kim DW, Soverini S, Haque A, Shou Y, Woodman RC, Kantarjian HM, Martinelli G, Radich JP, Saglio G, Hochhaus A, Hughes TP, Muller MC: Initial molecular response at 3 months may predict both response and event-free survival at 24 months in imatinib-resistant or -intolerant patients with Philadelphia chromosome -positive chronic myeloid leukemia in chronic phase treated with nilotinib. J Clin Oncol. 2012, 30: 4323-4329. 10.1200/JCO.2011.40.5217.View ArticlePubMedGoogle Scholar
  22. Hughes TP, Hochhaus A, Branford S, Muller MC, Kaeda JS, Foroni L, Druker BJ, Guilhot F, Larson RA, O’Brien SG, Rudoltz MS, Mone M, Wehrle E, Modur V, Goldman JM, Radich JP: Long-term prognostic significance of early molecular response to imatinib in newly diagnosed chronic myeloid leukemia: an analysis from the International Randomized Study of Interferon versus STI571 (IRIS). Blood. 2010, 116: 3758-3765. 10.1182/blood-2010-03-273979.View ArticlePubMedPubMed CentralGoogle Scholar
  23. Hehlmann R, Lauseker M, Jung-Munkwitz S, Leitner A, Muller MC, Pletsch N, Proetel U, Haferlach C, Schlegelberger B, Balleisen L, Hanel M, Pfirrmann M, Krause SW, Nerl C, Pralle H, Gratwohl A, Hossfeld DK, Hasford J, Hochhaus A, Saussele S: Tolerability-adapted imatinib 800 mg/d versus 400 mg/d versus 400 mg/d plus interferon-α in newly diagnosed chronic myeloid leukemia. J Clin Oncol. 2011, 29: 1634-1642. 10.1200/JCO.2010.32.0598.View ArticlePubMedGoogle Scholar
  24. Hughes T, Saglio G, Branford S, Soverini S, Kim DW, Müller MC, Martinelli G, Cortes J, Beppu L, Gottardi E, Kim D, Erben P, Shou Y, Haque A, Gallagher N, Radich J, Hochhaus A: Impact of baseline BCR-ABL mutations on response to nilotinib in patients with chronic myeloid leukemia in chronic phase. J Clin Oncol. 2009, 27: 4204-4210. 10.1200/JCO.2009.21.8230.View ArticlePubMedGoogle Scholar
  25. Branford S, Fletcher L, Cross NC, Muller MC, Hochhaus A, Kim DW, Radich JP, Saglio G, Pane F, Kamel-Reid S, Wang YL, Press RD, Lynch K, Rudzki Z, Goldman JM, Hughes T: Desirable performance characteristics for BCR-ABL measurement on an international reporting scale to allow consistent interpretation of individual patient response and comparison of response rates between clinical trials. Blood. 2008, 112: 3330-3338. 10.1182/blood-2008-04-150680.View ArticlePubMedGoogle Scholar
  26. Zhang T, Grenier S, Nwachukwu B, Wei C, Lipton JH, Kamel-Reid S, Association for Molecular Pathology Hematopathology Subdivision: Inter-laboratory comparison of chronic myeloid leukemia minimal residual disease monitoring: summary and recommendations. J Mol Diagn. 2007, 9: 421-430. 10.2353/jmoldx.2007.060134.View ArticlePubMedPubMed CentralGoogle Scholar
  27. Davidian M, Giltinan DM: Nonlinear models for repeated measurement data: an overview and update. JABES. 2003, 8: 387-419. 10.1198/1085711032697.View ArticleGoogle Scholar
  28. Sheiner BL, Beal SL: Evaluation of methods for estimating population pharmacokinetic parameters. II. Biexponential model and experimental pharmacokinetic data. J Pharmacokinet Biopharm. 1981, 9: 635-651. 10.1007/BF01061030.View ArticlePubMedGoogle Scholar
  29. Hughes TP, Kaeda J, Branford S, Rudzki Z, Hochhaus A, Hensley ML, Gathmann I, Bolton AE, van Hoomissen IC, Goldman JM, Radich JP: Frequency of major molecular responses to imatinib or interferon alfa plus cytarabine in newly diagnosed chronic myeloid leukemia. N Engl J Med. 2003, 349: 1423-1432. 10.1056/NEJMoa030513.View ArticlePubMedGoogle Scholar
  30. Kantarjian HM, Hochhaus A, Saglio G, De Souza C, Flinn IW, Stenke L, Goh Y, Rosti G, Nakamae H, Gallagher NJ, Hoenekopp A, Blakesley RE, Larson RA, Hughes TP: Nilotinib versus imatinib for the treatment of patients with newly diagnosed chronic phase, Philadelphia chromosome-positive, chronic myeloid leukaemia: 24-month minimum follow-up of the phase 3 randomised ENESTnd trial. Lancet Oncol. 2011, 12: 841-851. 10.1016/S1470-2045(11)70201-7.View ArticlePubMedGoogle Scholar
  31. Kantarjian HM, Shah NP, Cortes JE, Baccarani M, Agarwal MB, Undurraga MS, Wang J, Kassack Ipiña JJ, Kim D, Ogura M, Pavlovsky C, Junghanss C, Milone JH, Nicolini FE, Robak T, Van Droogenbroeck J, Vellenga E, Bradley-Garelik MB, Zhu C, Hochhaus A: Dasatinib or imatinib in newly diagnosed chronic phase chronic myeloid leukemia: 2-year follow-up from a randomized phase 3 trial (DASISION). Blood. 2012, 119: 1123-1129. 10.1182/blood-2011-08-376087.View ArticlePubMedGoogle Scholar
  32. Beal SL: Ways to fit a PK model with some data below the quantification limit. J Pharmacokinet Pharmacodyn. 2001, 28: 481-504. 10.1023/A:1012299115260.View ArticlePubMedGoogle Scholar
  33. Baccarani M, Cortes J, Pane F, Niederwieser D, Saglio G, Apperley J, Cervantes F, Deininger M, Gratwohl A, Guilhot F, Hochhaus A, Horowitz M, Hughes T, Kantarjian H, Larson R, Radich J, Simonsson B, Silver RT, Goldman J, Hehlmann R: Chronic myeloid leukemia: an update of concepts and management recommendations of European LeukemiaNet. J Clin Oncol. 2009, 27: 6041-6051. 10.1200/JCO.2009.25.0779.View ArticlePubMedGoogle Scholar
  34. Hochhaus A, Muller MC, Radich J, Branford S, Kantarjian HM, Hanfstein B, Rousselot P, Kim DW, Lipton JH, Bleickardt E, Lambert A, Hughes TP: Dasatinib-associated major molecular responses in patients with chronic myeloid leukemia in chronic phase following imatinib failure: response dynamics and predictive value. Leukemia. 2009, 23: 1628-1633. 10.1038/leu.2009.156.View ArticlePubMedGoogle Scholar
  35. Quintas-Cardama A, Kantarjian H, Jones D, Shan J, Borthakur G, Thomas D, Kornblau S, O’Brien S, Cortes J: Delayed achievement of cytogenetic and molecular response is associated with increased risk of progression among patients with chronic myeloid leukemia in early chronic phase receiving high-dose or standard-dose imatinib therapy. Blood. 2009, 113: 6315-6321. 10.1182/blood-2008-07-166694.View ArticlePubMedPubMed CentralGoogle Scholar
  36. Marin D, Ibrahim AR, Lucas C, Gerrard G, Wang L, Szydlo RM, Clark RE, Apperley JF, Milojkovic D, Bua M, Pavlu J, Paliompeis C, Reid A, Rezvani K, Goldman JM, Foroni L: Assessment of BCR-ABL1 transcript levels at 3 months is the only requirement for predicting outcome for patients with chronic myeloid leukemia treated with tyrosine kinase inhibitors. J Clin Oncol. 2012, 30: 232-238. 10.1200/JCO.2011.38.6565.View ArticlePubMedGoogle Scholar
  37. Hochhaus A, Saglio G, Chuah C, Pavlovsky C, Garelick MBB, Lambert A, Shah N: Dasatinib and imatinib-induced reductions in BCR-ABL transcript levels below 10% at 3 months are associated with improved responses in patients with newly diagnosed chronic myeloid leukemia in chronic phase (CML-CP): analysis of molecular response kinetics in the DASISION trial [abstract]. Blood. 2011, 118: 1190-1191. 10.1182/blood-2011-06-359836.View ArticleGoogle Scholar
  38. Milojkovic D, Apperley JF, Gerrard G, Ibrahim AR, Szydlo R, Bua M, Reid A, Rezvani K, Foroni L, Goldman J, Marin D: Responses to second-line tyrosine kinase inhibitors are durable: an intention-to-treat analysis in chronic myeloid leukemia patients. Blood. 2012, 119: 1838-1843. 10.1182/blood-2011-10-383000.View ArticlePubMedGoogle Scholar
  39. Marin D, Hedgley C, Clark RE, Apperley J, Foroni L, Milojkovic D, Pocock C, Goldman JM, O’Brien S: Predictive value of early molecular response in patients with chronic myeloid leukemia treated with first line dasatinib. Blood. 2012, EPub Ahead of Print May 29Google Scholar
  40. Branford S, Yeung DT, Prime JA, Choi SY, Bang JH, Park JE, Kim DW, Ross DM, Hughes TP: BCR-ABL1 doubling-times more reliably assess the dynamics of CML relapse compared with the BCR-ABL1 fold-rise: implications for monitoring and management. Blood. 2012, 119: 4264-4271. 10.1182/blood-2011-11-393041.View ArticlePubMedGoogle Scholar
  41. Pre-publication history

    1. The pre-publication history for this paper can be accessed here:


© Stein et al.; licensee BioMed Central Ltd. 2013

This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.