Rapid optimization of drug combinations for the optimal angiostatic treatment of cancer

Drug combinations can improve angiostatic cancer treatment efficacy and enable the reduction of side effects and drug resistance. Combining drugs is non-trivial due to the high number of possibilities. We applied a feedback system control (FSC) technique with a population-based stochastic search algorithm to navigate through the large parametric space of nine angiostatic drugs at four concentrations to identify optimal low-dose drug combinations. This implied an iterative approach of in vitro testing of endothelial cell viability and algorithm-based analysis. The optimal synergistic drug combination, containing erlotinib, BEZ-235 and RAPTA-C, was reached in a small number of iterations. Final drug combinations showed enhanced endothelial cell specificity and synergistically inhibited proliferation (p < 0.001), but not migration of endothelial cells, and forced enhanced numbers of endothelial cells to undergo apoptosis (p < 0.01). Successful translation of this drug combination was achieved in two preclinical in vivo tumor models. Tumor growth was inhibited synergistically and significantly (p < 0.05 and p < 0.01, respectively) using reduced drug doses as compared to optimal single-drug concentrations. At the applied conditions, single-drug monotherapies had no or negligible activity in these models. We suggest that FSC can be used for rapid identification of effective, reduced dose, multi-drug combinations for the treatment of cancer and other diseases. Electronic supplementary material The online version of this article (doi:10.1007/s10456-015-9462-9) contains supplementary material, which is available to authorized users.


Introduction
Anti-angiogenic therapies are routinely used in the treatment of various cancers [1][2][3]. Their contribution to the prolongation of patient survival, however, is often limited mainly due to disease and patient heterogeneity [4,5], toxicity [6], induction of metastasis [7] and drug resistance [8,9]. Redundancy of growth factor signaling pathways makes angiogenesis a robust physiological function [10,11], where targeting multiple pathways with drug combinations may be necessary for efficient therapy [12]. Although difficult to predict, in such drug combinations one may encounter synergistic, additive or antagonistic interactions between drugs. Synergistic interactions can Electronic supplementary material The online version of this article (doi:10.1007/s10456-015-9462-9) contains supplementary material, which is available to authorized users. lead to effective angiogenesis inhibition at reduced doses as compared to single-drug therapies. Combination strategies may thus lead to enhanced efficacy [13,14] with limited side effects [15] and reduced probability of developing drug resistance [16,17].
Combinations of anti-angiogenic drugs have often resulted in significant clinical toxicity [18], even when designed to target complementary pathways [19]. This is because drugs to be combined are frequently selected based on their success as single agents [20] and tend to be used in combination at their maximum tolerated single agent doses, thus increasing the risk of toxicity and resistance [21]. When trying to identify an optimal combination starting from, for instance, 10 drugs at 5 doses, one will have to test nearly 10 million (5 10 ) combinations. To overcome this challenge, we employed a feedback system control (FSC) technique to rapidly identify the most powerful drug combinations with minimal experimental effort [22] (Supplementary Methods). In combination with the differential evolution (DE) algorithm [23], an iterative approach of experimental testing in an endothelial cell viability assay and mathematical analysis (a process of selection, where only the permutations which improve the system's response are maintained) drove the system to converge toward an optimal solution, i.e., maximal inhibition of endothelial cell growth. Although others have tried to optimize drug combinations [24][25][26][27], see Supplementary Methods, the advantage of our approach is that FSC is phenotypically driven, i.e., no mechanistic information is required in order to rapidly identify experimentally verifiable optimal drug combinations [22].
The aim of the present study was to find an optimal lowdose, synergistic anti-angiogenic drug combination using the FSC technology, and to validate this drug combination in preclinical tumor models. The FSC technique, together with a second-order linear regression model to allow for elimination of less effective drugs, resulted in the identification of the optimal low-dose combination containing erlotinib (EGFR inhibitor [28]), RAPTA-C (histone inactivator [29]) and BEZ-235 (a dual PI3K/mTOR inhibitor [30]). This final drug combination synergistically inhibited ECRF24 viability, while having minimal effects on nonendothelial cell types. We successfully translated this in vitro optimized drug combination to inhibit tumor growth in two preclinical tumor models.

Materials and methods
Cell viability, migration and apoptosis assay Cell viability and migration assays were performed as previously described. Cells were seeded in a 96-well culture plate at a density of 2.5-10 9 10 3 cells/well. Cells were incubated with drugs for 72 h (for drug acquisition and cells and culture conditions, see Supplementary Methods). Drugs were premixed in culture medium and applied at the doses provided in Table 1. Cell viability was assessed using the CellTiter-Glo luminescence assay (Promega, Madison, WI, USA). For migration assays, ECRF24 and 786-O were seeded in 96-well cell culture plates (3 9 10 4 cells/well) 24 h prior to making the scratch (Peira Scientific Instruments, Beerse, Belgium). Drugs were premixed in culture medium and applied at doses indicated in Supplementary Fig. S2A. Images were automatically captured on a Leica DMI3000 microscope (Leica, Rijswijk, Netherlands) at 59 magnification with Universal Grab 6.3 software (DCILabs, Keerbergen, Belgium). Scratch sizes were determined at t = 0 h and t = 7 h using Scratch Assay 6.2 (DCILabs), and values reported represented the absolute closure of the scratch (initially subtracting the final scratch area). Apoptosis was measured after drug exposure, trypsinization and incubation with propidium iodide (PI, 20 µg/ml) in DNA extraction buffer [31], by flow cytometry. Tip cells were flow cytometrically quantified by CD34 [32] staining, and morphology was studied in vivo using the CAM assay [33] (see Supplementary Methods).

The feedback system control (FSC) technique and data modeling
The FSC technique was employed as previously described [34,35]. FSC was implemented using the DE algorithm, and two separate optimizations were performed with the cellular outputs of ECRF24 cell viability (proliferation) and migration assays. Nineteen drug combinations were tested per iteration, and 10 iterations were performed in each optimization until a plateau in the best output value was reached. For dilutions and culture conditions, see Supplementary Methods. The cells were incubated in 50 µl of each combination for 72 h in the viability assay or for 7 h in the migration assay. Second-order linear regression models were generated using the data obtained from each optimization. Data were modeled using real concentration values, and both concentration values and cell viability output data were transformed using the z score function in MATLAB. For detailed description, see Supplementary Methods.

Human ovarian carcinoma grown on the chicken chorioallantoic membrane (CAM)
Human ovarian carcinoma tumors were implanted on the CAM as previously described [36]. On embryo development day (EDD) 7, 1 9 10 6 A2780 carcinoma cells were prepared as a spheroid in a 25-µl hanging drop and were transplanted onto the CAM surface 3 h after preparation. Treatment began 3 days after tumor implantation (EDD10) when vascularized tumors were visible. Drug combinations were freshly prepared and administered as a 20-µl intravenous injection. Treatment was performed twice, and tumor growth was monitored and measured daily (volume = width 2 9 length 9 0.52).

Colorectal carcinoma xenograft model
Female Swiss nu/nu mice aged 6-8 weeks were obtained from Charles River (weight 20-30 g). Mice were inoculated in the right flank with 100 µl DMEM with 1 million LS174T cells. LS174T cells were obtained from ECACC, Salisbury, UK (authentication by STR PCR), and were used within 6 months of resuscitation. Palpable tumors were present within 3-5 days, at which time treatment was initiated. Mice were treated daily by oral gavage and i.p. injection as indicated ( Table 2) and were monitored daily for tumor size and body weight (see Supplementary Methods).

Statistical analysis
Values are given as mean values ± SD. Statistical analysis was performed using a two-sided student's t test and the two-way ANOVA assay. *p \ 0.05 and **p \ 0.01 were considered statistically significant.

Selection of drug combinations by the FSC technique
Nine drugs targeting a broad spectrum of endothelial cell signaling pathways (Supplementary Fig. S1; Supplementary Methods) were selected for FSC-based screening ( Fig. 1a): anginex (1), bevacizumab (2), axitinib (3), erlotinib (4), anti-HMGB1 Ab (5), sunitinib (6), antivimentin Ab (7), RAPTA-C (8) and BEZ-235 (9). Singledrug dose-response curves were generated for both cell viability (example for sunitinib provided, Fig. 1b) and migration, using in vitro bioassays ( Supplementary Fig.  S2A). The optimization was carried out with each compound at four low doses. The highest concentration, dose 3 or ED 10 , was the dose where 10 % of the maximal response was observed, dose 2 or ED 5 , where 5 % of the maximal response was observed, dose 1 or ED 0 , represented half the maximal dose where no effect was observed, and dose 0, where no drug was present (Table 1). Starting from randomly selected drug combinations (Fig. 1a, yellow arrow), the FSC technique implements an algorithm-guided closedloop feedback search to iteratively optimize the results of an in vitro cell assay (blue arrows). The box plot in Fig. 1c provides the median and interquartile ranges of the output results of the drug combinations identified by the end of each iterative cycle of the FSC optimization. After 10 iterations of 19 drug mixtures, the optimization goal was reached, i.e., no further improvement of the lowest output efficacy could be achieved, indicating that the maximum activity (approx. 70 % inhibition) had been reached.
The data obtained from this optimization process were used to build a second-order stepwise linear regression model [37] (Supplementary Methods) to determine the relative importance of the individual drugs. This model generated regression coefficients (Fig. 1d) corresponding to single-drug linear effects (left panel), two-drug pair-wise interaction effects (middle panel) and single-drug quadratic effects (right panel). Compounds with the largest negative regression coefficients, i.e., axitinib, erlotinib, RAPTA-C and BEZ-235, inhibited ECRF24 viability most effectively (Fig. 1d, green arrows). A regression model containing all regression coefficients (i.e., a non-stepwise linear regression model) is provided in Supplementary Fig. S2B.
In a parallel approach, we also investigated the best drug combinations for ECRF24 migration inhibition. Even though single drugs generally showed a stronger response in the cell migration assay ( Supplementary Fig. S2A), the process of migration was less affected, reaching a maximum effect of 40 % inhibition in the given conditions ( Supplementary Fig. S3B). The optimization of EC migration inhibition was not further pursued. Yet, regression analysis also revealed strong single-drug linear and quadratic contributions for erlotinib, RAPTA-C and BEZ-235.

Refined search leads to further optimized synergistic drug combinations
Subsequently, a second FSC-based optimization was performed with the above-selected compounds, i.e., axitinib (3), erlotinib (4), RAPTA-C (8) and BEZ-235 (9), each now considered at five drug doses with a maximum activity of 25 % at the highest dose ( Fig. 2a; single-drug effects in Supplementary Fig. S3). The most effective combinations resulting in more than 50 % inhibition of ECRF24 cell viability identified in the second screen are provided in Linear regression modeling of data showed the singledrug linear contributions of all compounds, as well as the single-drug quadratic effect of 4, to be significant (Fig. 2b, green arrows). Response surfaces (Fig. 2c) provide a visual representation of the relationship between the system output (EC viability) and the varying dose of only two drugs in the combination. Interestingly, surfaces containing 3 (bottom row) show that increasing the dose of 3 does not enhance the combination efficacy (red). These response surfaces show a relatively smooth response when doses of the given drugs are varied. This "smoothness" indicates that a moderate change in the dose of a single drug in the range of the experimental conditions near the optimal output investigated will most likely not result in a significant change in the output response. This implies a certain amount of "stability" in the optimal drug mixture which may facilitate its translation to different models.  Selected optimized drug combinations exhibit enhanced endothelial cell specificity The optimized drug combinations I-VIII (Fig. 2a) and corresponding single drugs were tested for viability of different cell types and shown in comparison with ECRF24 ( Fig. 3a; Supplementary Fig. S5). The activity in ECRF24 was confirmed in primary ECs (HUVEC, Fig. 3a) and was much stronger than that of non-malignant cell types (adult human dermal fibroblasts (HDFa), human peripheral blood mononuclear cells ( (Fig. 3c).
Finally, the effect of combination VI on the inhibition of tip cells was assessed both in vitro and in vivo. FACS analysis shows a reduction in the number of CD34 + tip cells [32,38,39] after treating HUVEC cultures in vitro with VI. Quantification of these results is provided in the bar graph as compared to single-drug treatments (Fig. 3d,  top left). Additionally, CD34 + tip cells treated with VI present with a clearly different cellular organization of the actin fibers stained with phalloidin, as compared to control cells (Fig. 3d, top right), compatible with decreased migratory activity.
The inhibition of tip cells in vivo was shown in the chicken chorioallantoic membrane (CAM) model following vaso-occlusive Visudyne ® -photodynamic therapy (V-PDT) [33]. In the control treatment group, vascular sprouts led by tip cells can be seen growing into the treated area 24 h after V-PDT treatment starting the revascularization of the tissue (Fig. 3d, bottom left, green arrows). A significant reduction in the number of sprouting tip cells can be seen in the group treated with VI following V-PDT (Fig. 3d, bottom right, red arrows).
Successful translation of optimal drug combinations into in vivo cancer models A2780 cells were transplanted onto the chorioallantoic membrane of the chicken embryo and were allowed to grow forming vascularized tumors. Tumors were treated Fig. 1 Optimization of the inhibition of endothelial cell viability. a Schematic diagram of the FSC technique loop (blue arrow loop) and modeling (green arrow loop) used for in vitro drug optimization. FSC starts with randomly selected drug combinations (yellow arrow) and implements an algorithm-guided closed-loop feedback search to iteratively optimize the results of an in vitro cell assay (blue arrows). Once a plateau in the output is reached, the data obtained from the optimization are used to model the system, analyze drug interactions and eliminate certain drugs (green arrows). Using a refined set of drugs, the drug combination is again optimized with FSC (blue arrows). b Dose-response curve of sunitinib for cell viability bioassay. c Output results (in vitro EC viability, represented as a percentage of the control) for the 10 iterations of the FSC optimization performed. Box plots provide median and interquartile ranges of the cell response to the 19 best drug combinations identified by the end of each iterative cycle of the FSC optimization. Dotted lines, representing maximum and minimum (red) output values, showed no improvement in the best-optimized combination over iterations 8-10.  with combinations I, II, VI, VII and VIII by intravenous injection on treatment days 1 and 2 (Fig. 4a, red arrows). Doses (subsequently identified by a subscript) were translated to this model maintaining the drug dose ratios and taking into account the single-drug efficacy in this model ( Fig. 4b; Table 2, Supplementary Methods). Drug combination VII (3 3 +4 4 +8 1 +9 2 ) synergistically (CI 0.66) inhibited tumor growth by 87 % (* p \ 0.03, Fig. 4a). Based on results in Fig. 3a, this activity could be due to the dual action on both ECRF24 and A2780 cells. Combination VI (4 1 +8 2 +9 4 ) synergistically inhibited tumor growth by 68 % (Fig. 4a, c, **p \ 0.002, CI 0.34) through mainly anti-angiogenic activity (compare Fig. 3a). Of note, none of these doses inhibited tumor growth significantly when applied individually (Fig. 4b). As group VII experienced weight loss (Fig. 4d, **p\0.004), it was not further examined. Microvessel density (MVD) assessment (Fig. 4e) revealed that control tumors were well vascularized. MVD was 50-60 % lower in tumors treated with VI (**p\0.008) and VIII (*p = 0.01, Fig. 4f). Based on these data, VI was selected as the most promising combination.

Discussion
Next to the identification of novel targets as well as endogenous and synthetic novel angiogenesis inhibitors [40,41], the combination of therapies is globally seen as a b Fig. 3 Validation of the best drug combinations. The effects of the most promising drug combinations (I-VIII from Fig. 2)  promising strategy to improve cancer therapy. The FSC technique was used to navigate through the large parametric space of nine compounds, each considered at four doses, aiming for an optimal angiostatic drug combination. Using a simple in vitro endothelial cell (EC) viability bioassay as the output, an optimal low-dose drug combinations containing axitinib (3), erlotinib (4), RAPTA-C (8) and BEZ-235 (9) was found. The most efficient of these combinations was also effectively inhibiting cancer in two in vivo animal models. We observed that (1) while some drugs showed synergistic interactions, others showed additive or even antagonistic behavior, (2) the observed synergy was drug dose ratio dependent, (3) the combination of angiostatic drugs enhanced endothelial specificity, (4) screening on EC migration did not identify highly efficient drug combinations, and (5) in vitro optimized anti-angiogenic drug combinations translated to anti-angiogenic anticancer effects in vivo. We previously demonstrated that multi-drug effects can be expressed by a quadratic relationship of the drug-drug interactions [42], which was confirmed in bacterial systems [43]. Here, we have further demonstrated that the response surface for the whole range of drugs and drug doses applied can be expressed as a second-order equation that can be used to formulate optimal drug combinations. The results of this regression modeling (Supplementary Methods) permitted us to eliminate sunitinib (6) a compound which is known to have a similar target profile as axitinib (3) (note that both inhibit signaling of VEGF and PDGF [44]). The exclusion of sunitinib over axitinib appears justifiable, as axitinib is known to be a more selective TKI with stronger affinity for the same targets [44]. Similarly, the exclusion of bevacizumab (2) was expected, as it is known that EC does not use VEGF as an autocrine growth factor, and tumor angiogenesis is mainly driven by tumor produced VEGF [10].
The four drugs with significant inhibitory single-drug linear contributions to cell viability were compounds axitinib, erlotinib, RAPTA-C and BEZ-235 (Fig. 1d). In terms of intracellular signaling, this combination of drugs appears to make sense in retrospect. EGFR targeting by erlotinib and VEGFR targeting by axitinib result in inhibition of two largely synergistic and widely used cellular signaling pathways, i.e., the PI3K/AKT/mTOR and the ras/raf/MEK/ MAPK signaling pathway, respectively. Since mTORC1 and mTORC2 belong to the PI3K/AKT pathway, one would expect that both signaling pathways are inhibited by EGFR and VEGFR inhibitors. It is also expected that a drug that targets histone proteins, such as RAPTA-C [45], can reinforce the angiostatic effect, as intervention with histone-DNA interactions is known to be angiostatic from the many reports on histone deacetylase inhibitors [46][47][48].  Table 2. b Inhibition of tumor growth by single compounds at indicated doses. Data points represent the average tumor volume as a percentage of the final CTRL volume per experiment, and error bars represent the SEM; N = 3-9. *p \ 0.05 and **p \ 0.01 (two-way ANOVA). "S" indicates synergy (CI \ 1). c Representative images of vehicle-treated (CTRL) and tumors treated with drug combination VI on the last experiment day.
d Representative images of immunohistochemical staining for the EC marker CD31 and corresponding quantification of microvessel density, measured as the number of vessels per mm 2 and presented as a percentage of the CTRL. Results show significantly reduced microvessel density in tumors treated with VI. The bar in the lower panel image represents 0.2 mm and is valid for both images. The whole images were linearly adjusted for brightness and contrast. e Body weight change during the experiment. *p \ 0.05 and **p \ 0.01 (student t test). mTOR and EGFR inhibitors have already been identified as a synergistic combination in various cancer cell types [49,50], despite clinically observed toxicity [51,52]. Based on the analysis of the response surfaces of the second-order linear regression model generated from the four-drug optimization data (Fig. 2c) and embryo weight loss observed in the CAM model (Fig. 4d), axitinib was eliminated from further investigation. Thus, the optimal drug combination containing erlotinib, RAPTA-C and BEZ-235 was identified. It allowed for dose reductions of 5-, 11-and 6-fold, respectively, as compared to the equivalent single-drug dose efficiencies in vitro. Interestingly, enhanced EC specificity was observed for the combinations when compared to the individual compounds. This is another indication that the parallel blocking of multiple angiogenesis pathways can result in synergism for the angiostatic outcome. Simultaneous targeting of different signaling pathways may limit the probability of cells to develop acquired resistance [16].
The migration-based optimization screen failed to reach effective combinations ( Supplementary Fig. S3B). This may suggest that proliferation is more dominant in the process of angiogenesis than cell migration, which has also been proposed by others [53]. The same may also be reflected by clinical trials, where proliferation inhibitors (such as sunitinib and BEZ-235) were more successful than migration inhibitors (the α v β 3 inhibitor cilengitide [54] and the α 5 β 1 antibody volociximab [55]). Another possible explanation for enhanced success with the proliferation assay over the migration assay may be the selective nature of synergistic drug interactions. As seen by Lehar et al. [15], drug combinations could attain greater selectivity. They suggested that "synergistic combinations tend to be more specific to particular cellular phenotypes than are single drugs." This may explain the preferred success of the screen through selection on basis of proliferation, rather than cell migration.
The optimal drug mixture inhibited tumor growth by approximately 80 %, most likely by an inhibitory effect on angiogenesis. Although the detailed mechanism of combination therapy still needs to be understood, the induction of apoptosis as well as the inhibition of tip cells shows part of the effector mechanism. Targeting of tip cells may be another attractive strategy as these cells are indispensable for sprouting angiogenesis. The results provide a promising option for future clinical anti-angiogenic applications.
One might expect that the differences in pharmacokinetics between the components of the drug mixture may interfere with obtaining good results in vivo. Our results imply that (1) the best drug combinations found show smooth response surfaces (Fig. 2c, i.e., moderate changes in drug ratios do not significantly change the output), (2) response surfaces, giving a mathematical description of the magnitude of the interaction for all drug pairs, confirmed in vivo treatment outcome (compare Figs. 2c, 5a), and (3) EC viability observed in vitro seems to be a relatively good parameter for translation to vascular density reduction and tumor growth inhibition in vivo.
The current study shows that FSC applied in vitro can be used for the fast and reliable identification of potent, lowdose angiostatic drug combinations in vivo. It is likely that combining the optimal anti-angiogenic compounds with other treatment strategies may lead to even better cancer treatment outcomes. The capacity to normalize the tumor vasculature provides angiostatic strategies with outstanding combination therapy features [56,57]. The impact of the method also lies in the fact that it can be applied in a variety of situations, e.g., for finding drug mixtures directly targeting tumor cells. This strategy would then also offer the opportunity for a personalized approach, by performing a drug screen on freshly isolated tumor cells from a patient biopsy. Such strategy would also depend on faster screening methods. We are currently working on improvement of the FSC method in order to make the selection procedure faster. In conclusion, designing effective, synergistic and specific multi-component drug combinations may become a key approach in developing new therapies for cancer and other diseases.