Investigation of the robustness of two models for assessing synergy in pre‐clinical drug combination studies

Pre‐clinical studies may be used to screen for synergistic combinations of drugs. The types of in vitro assays used for this purpose will depend upon the disease area of interest. In oncology, one frequently used study measures cell line viability: cells placed into wells on a plate are treated with doses of two compounds, and cell viability is assessed from an optical density measurement corrected for blank well values. These measurements are often transformed and analysed as cell survival relative to untreated wells. The monotherapies are assumed to follow the Hill equation with lower and upper asymptotes at 0 and 1, respectively. Additionally, a common variance about the dose–response curve may be assumed. In this paper, we consider two models for incorporating synergy parameters. We investigate the effect of different models of biological variation on the assessment of synergy from both of these models. We show that estimates of the synergy parameters appear to be robust, even when estimates of the other model parameters are biased. Using untransformed measurements provides better coverage of the 95% confidence intervals for the synergy parameters than using transformed measurements, and the requirement to fit the upper asymptote does not cause difficulties. Assuming homoscedastic variances appears to be robust. The added complexity of determining and fitting an appropriate heteroscedastic model does not seem to be justified. Copyright © 2013 John Wiley & Sons, Ltd.     

Experimental designs for detecting synergy and antagonism between two drugs in a pre‐clinical study

The identification of synergistic interactions between combinations of drugs is an important area within drug discovery and development. Pre‐clinically, large numbers of screening studies to identify synergistic pairs of compounds can often be ran, necessitating efficient and robust experimental designs. We consider experimental designs for detecting interaction between two drugs in a pre‐clinical in vitro assay in the presence of uncertainty of the monotherapy response. The monotherapies are assumed to follow the Hill equation with common lower and upper asymptotes, and a common variance. The optimality criterion used is the variance of the interaction parameter. We focus on ray designs and investigate two algorithms for selecting the optimum set of dose combinations. The first is a forward algorithm in which design points are added sequentially. This is found to give useful solutions in simple cases but can lack robustness when knowledge about the monotherapy parameters is insufficient. The second algorithm is a more pragmatic approach where the design points are constrained to be distributed log‐normally along the rays and monotherapy doses. We find that the pragmatic algorithm is more stable than the forward algorithm, and even when the forward algorithm has converged, the pragmatic algorithm can still out‐perform it. Practically, we find that good designs for detecting an interaction have equal numbers of points on monotherapies and combination therapies, with those points typically placed in positions where a 50% response is expected. More uncertainty in monotherapy parameters leads to an optimal design with design points that are more spread out. Copyright © 2015 John Wiley & Sons, Ltd.     

Adaptive phase I/II clinical trials for drug combination assessment in oncology using the outcomes of each cycle

Many new anticancer agents can be combined with existing drugs, as combining a number of drugs may be expected to have a better therapeutic effect than monotherapy owing to synergistic effects. Furthermore, to drive drug development and to reduce the associated cost, there has been a growing tendency to combine these as phase I/II trials. With respect to phase I/II oncology trials for the assessment of dose combinations, in the existing methodologies in which efficacy based on tumor response and safety based on toxicity are modeled as binary outcomes, it is not possible to enroll and treat the next cohort of patients unless the best overall response has been determined in the current cohort. Thus, the trial duration might be potentially extended to an unacceptable degree. In this study, we proposed a method that randomizes the next cohort of patients in the phase II part to the dose combination based on the estimated response rate using all the available observed data upon determination of the overall response in the current cohort. We compared the proposed method to the existing method using simulation studies. These demonstrated that the percentage of optimal dose combinations selected in the proposed method is not less than that in the existing method and that the trial duration in the proposed method is shortened compared to that in the existing method. The proposed method meets both ethical and financial requirements, and we believe it has the potential to contribute to expedite drug development.