Pharmacodynamic models provide inroads to understanding key mechanisms of action and may significantly improve patient outcomes in cancer with improved ability to determine therapeutic benefit. Additionally, these models may also lead to insights into potential biomarkers that can be utilized for prediction in prognosis and therapeutic decisions. As an example of this potential, here we present an advanced computational Ordinary Differential Equation (ODE) model of PARP1 signalling and downstream effects due to its inhibition. The model has been validated experimentally and further evaluated through a global sensitivity analysis. The sensitivity analysis uncovered two model parameters related to protein synthesis and degradation rates that were also found to contribute the most variability to the therapeutic prediction. Because this variability may define cancer patient subpopulations, we interrogated genomic, transcriptomic, and clinical databases, to uncover a biomarker that may correspond to patient outcomes in the model. In particular, GSPT2, a GTPase with translation function, was discovered and if mutations serve to alter catalytic activity, its presence may explain the variability in the model′s parameters. This work offers an analysis of ODE models, inclusive of model development, sensitivity analysis, and ensuing experimental data analysis, and demonstrates the utility of this methodology in uncovering biomarkers in cancer.