Abstract
Optimal integration of molecularly targeted therapies, such as tyrosine kinase inhibitors (TKI), with concurrent chemotherapy and radiation (CRT) to improve outcomes in genotype-defined cancers remains a current challenge in clinical settings. Important questions regarding optimal scheduling and length of induction period for neoadjuvant use of targeted agents remain unsolved and vary among clinical trial protocols. Here, we develop and validate a biomathematical framework encompassing drug resistance and radiobiology to simulate patterns of local versus distant recurrences in a non-small cell lung cancer (NSCLC) population with mutated EGFR receiving TKIs and CRT. Our model predicted that targeted induction before CRT, an approach currently being tested in clinical trials, may render adjuvant targeted therapy less effective due to proliferation of drug-resistant cancer cells when using very long induction periods. Furthermore, simulations not only demonstrated the competing effects of drug-resistant cell expansion versus overall tumor regression as a function of induction length, but also directly estimated the probability of observing an improvement in progression-free survival at a given cohort size. We thus demonstrate that such stochastic biological simulations have the potential to quantitatively inform the design of multimodality clinical trials in genotype-defined cancers. SIGNIFICANCE: A biomathematical framework based on fundamental principles of evolution and radiobiology for in silico clinical trial design allows clinicians to optimize administration of TKIs before chemoradiotherapy in oncogene-driven NSCLC.