Abstract Circulating tumor DNA assays are promising tools for the prediction of cancer treatment response. Here, we build a framework for the design of ctDNA biomarkers of therapy response that incorporate variations in ctDNA dynamics driven by specific treatment mechanisms. We develop mathematical models of ctDNA kinetics driven by tumor response to several therapy classes, and utilize them to simulate randomized virtual patient cohorts to test candidate biomarkers.