Multiobjective optimization identifies cancer-selective combination therapies

Abstract

Author summaryCancer is diagnosed in nearly 40% of people in the U.S at some point during their lifetimes. Despite decades of research to lower cancer incidence and mortality, cancer remains a leading cause of deaths worldwide. Therefore, new targeted therapies are required to further reduce the death rates and toxic effects of treatments. Here we developed a mathematical optimization framework for finding cancer-selective treatments that optimally use drugs and their combinations. The method uses multiobjective optimization to identify drug combinations that simultaneously show maximal therapeutic potential and minimal non-selectivity, to avoid severe side effects. Our systematic search approach is applicable to various cancer types and it enables optimization of combinations involving both targeted therapies as well as standard chemotherapies.Combinatorial therapies are required to treat patients with advanced cancers that have become resistant to monotherapies through rewiring of redundant pathways. Due to a massive number of potential drug combinations, there is a need for systematic approaches to identify safe and effective combinations for each patient, using cost-effective methods. Here, we developed an exact multiobjective optimization method for identifying pairwise or higher-order combinations that show maximal cancer-selectivity. The prioritization of patient-specific combinations is based on Pareto-optimization in the search space spanned by the therapeutic and nonselective effects of combinations. We demonstrate the performance of the method in the context of BRAF-V600E melanoma treatment, where the optimal solutions predicted a number of co-inhibition partners for vemurafenib, a selective BRAF-V600E inhibitor, approved for advanced melanoma. We experimentally validated many of the predictions in BRAF-V600E melanoma cell line, and the results suggest that one can improve selective inhibition of BRAF-V600E melanoma cells by combinatorial targeting of MAPK/ERK and other compensatory pathways using pairwise and third-order drug combinations. Our mechanism-agnostic optimization method is widely applicable to various cancer types, and it takes as input only measurements of a subset of pairwise drug combinations, without requiring target information or genomic profiles. Such data-driven approaches may become useful for functional precision oncology applications that go beyond the cancer genetic dependency paradigm to optimize cancer-selective combinatorial treatments.