Drugs bind to their target proteins, which interact with downstream effectors and ultimately perturb the transcriptome of a cancer cell. These perturbations reveal information about their source, i.e., drugs’ targets. Here, we investigate whether these perturbations and protein interaction networks can uncover drug targets and key pathways. We performed the first systematic analysis of over 500 drugs from the Connectivity Map. First, we show that the gene expression of drug targets is usually not significantly affected by the drug perturbation. Hence, expression changes after drug treatment on their own are not sufficient to identify drug targets. However, ranking of candidate drug targets by network topological measures prioritizes the targets. We introduce a novel measure, local radiality, which combines perturbed genes and functional interaction network information. The new measure outperforms other methods in target prioritization and proposes cancer-specific pathways from drugs to affected genes for the first time. Local radiality identifies more diverse targets with fewer neighbors and possibly less side effects.

Continuing tyrosine kinase inhibitor (TKI)-mediated targeting of the BCR-ABL1 oncoprotein is the standard therapy for chronic myeloid leukemia (CML) and allows for a sustained disease control in the majority of patients. While therapy cessation for patients appeared as a safe option for about half of those patients with optimal response, no systematic assessment of long-term TKI dose de-escalation has been made. We use a mathematical model to analyze and consistently describe biphasic treatment responses from TKI-treated patients from two independent clinical phase III trials. Scale estimates reveal that drug efficiency determines the initial response while the long-term behavior is limited by the rare activation of leukemic stem cells. We use this mathematical framework to investigate the influence of different dosing regimens on the treatment outcome. We provide strong evidence to suggest that TKI dose de-escalation (at least 50%) does not lead to a reduction of long-term treatment efficiency for most patients, who have already achieved sustained remission, and maintains the secondary decline of BCR-ABL1 levels. We demonstrate that continuous BCR-ABL1 monitoring provides patient-specific predictions of an optimal reduced dose without decreasing the anti-leukemic effect on residual leukemic stem cells. Our results are consistent with the interim results of the DESTINY trial and provide clinically testable predictions. Our results suggest that dose-halving should be considered as a long-term treatment option for CML patients with good response under continuing maintenance therapy with TKIs. We emphasize the clinical potential of this approach to reduce treatment-related side-effects and treatment costs.

Recent clinical findings in patients with chronic myeloid leukemia (CML) suggest that the risk of molecular recurrence after stopping tyrosine kinase inhibitor (TKI) treatment substantially depends on an individual's leukemia-specific immune response. However, it is still not possible to prospectively identify patients that will remain in treatment-free remission (TFR). Here, we used an ordinary differential equation model for CML, which explicitly includes an antileukemic immunologic effect, and applied it to 21 patients with CML for whom BCR-ABL1/ABL1 time courses had been quantified before and after TKI cessation. Immunologic control was conceptually necessary to explain TFR as observed in about half of the patients. Fitting the model simulations to data, we identified patientspecific parameters and classified patients into three different groups according to their predicted immune system configuration ("immunologic landscapes"). While one class of patients required complete CML eradication to achieve TFR, other patients were able to control residual leukemia levels after treatment cessation. Among them were a third class of patients that maintained TFR only if an optimal balance between leukemia abundance and immunologic activation was achieved before treatment cessation. Model simulations further suggested that changes in the BCR-ABL1 dynamics resulting from TKI dose reduction convey information about the patient-specific immune system and allow prediction of outcome after treatment cessation. This inference of individual immunologic configurations based on treatment alterations can also be applied to other cancer types in which the endogenous immune system supports maintenance therapy, long-term disease control, or even cure.Significance: This mathematical modeling approach provides strong evidence that different immunologic configurations in patients with CML determine their response to therapy cessation and that dose reductions can help to prospectively infer different risk groups.See related commentary by Triche Jr, p. 2083

Longitudinal monitoring of BCR-ABL transcript levels in peripheral blood of CML patients treated with tyrosine kinase inhibitors (TKI) revealed a typical biphasic response. Although second generation TKIs like dasatinib proved more efficient in achieving molecular remission compared to first generation TKI imatinib, it is unclear how individual responses differ between the drugs and whether mechanisms of drug action can be deduced from the dynamic data. We use time courses from the DASISION trial to address statistical differences in the dynamic response between first line imatinib vs. dasatinib treatment cohorts and we analyze differences between the cohorts by fitting an established mathematical model of functional CML treatment to individual time courses. On average, dasatinib-treated patients show a steeper initial response, while the long-term response only marginally differed between the treatments. Supplementing each patient time course with a corresponding confidence region, we illustrate the consequences of the uncertainty estimate for the underlying mechanisms of CML remission. Our model suggests that the observed BCR-ABL dynamics may result from different, underlying stem cell dynamics. These results illustrate that the perception and description of CML treatment response as a dynamic process on the level of individual patients is a prerequisite for reliable patient-specific response predictions and treatment optimizations.