The substantial spatial and temporal heterogeneity observed in patient tumors poses considerable challenges for the design of effective drug combinations with predictable outcomes. Currently, the implications of tissue heterogeneity and sampling bias during diagnosis are unclear for selection and subsequent performance of potential combination therapies. Here, we apply a multiobjective computational optimization approach integrated with empirical information on efficacy and toxicity for individual drugs with respect to a spectrum of genetic perturbations, enabling derivation of optimal drug combinations for heterogeneous tumors comprising distributions of subpopulations possessing these perturbations. Analysis across probabilistic samplings from the spectrum of various possible distributions reveals that the most beneficial (considering both efficacy and toxicity) set of drugs changes as the complexity of genetic heterogeneity increases. Importantly, a significant likelihood arises that a drug selected as the most beneficial single agent with respect to the predominant subpopulation in fact does not reside within the most broadly useful drug combinations for heterogeneous tumors. The underlying explanation appears to be that heterogeneity essentially homogenizes the benefit of drug combinations, reducing the special advantage of a particular drug on a specific subpopulation. Thus, this study underscores the importance of considering heterogeneity in choosing drug combinations and offers a principled approach toward designing the most likely beneficial set, even if the subpopulation distribution is not precisely known.systems biology | cancer | combination therapy G enetic intratumor heterogeneity has long been appreciated as present in cancer patients (1). Recent sequencing studies further revealed the extent of this tumor diversity, arising from highly complex clonal evolutionary processes (2). This phenomenon has been observed in many solid (3-5) and hematopoietic cancers (6-8). Moreover, treatments also may have dramatic effects on tumor composition-with a preexisting subclone at diagnosis often becoming dominant at relapse (9-12). To meet this challenge, rational drug combination treatments must be designed so as to account for intratumor heterogeneity to better predict reoccurrence.The history of theoretical studies aimed at drug optimization in cancer therapy is long. Some studies used differential equation models formulated as deterministic optimal control problems (13-16), whereas others used stochastic birth-death process models (17-21), to examine treatment regimens for tumors comprising a sensitive population along with a few resistant subpopulations. However, these studies, many of which dealt only with generic scenarios, focused primarily on drug scheduling, investigating the effects of frequency and dose intensity of drugs on resistance potential. Practical applications of such results have been limited, although some of these strategies recently were combined with experimental validation to examin...