BackgroundThyroid cancer affects over ½ million people in the U.S. and the incidence of thyroid cancer has increased worldwide at a rate higher than any other cancer, while survival has remained largely unchanged. The aim of this research was to develop, calibrate and verify a mathematical disease model to simulate the natural history of papillary thyroid cancer, which will serve as a platform to assess the effectiveness of clinical and cancer control interventions.MethodsHerein, we modeled the natural pre-clinical course of both benign and malignant thyroid nodules with biologically relevant health states from normal to detected nodule. Using established calibration techniques, optimal parameter sets for tumor growth characteristics, development rate, and detection rate were used to fit Surveillance Epidemiology and End Results (SEER) incidence data and other calibration targets.ResultsModel outputs compared to calibration targets demonstrating sufficient calibration fit and model validation are presented including primary targets of SEER incidence data and size distribution at detection of malignancy. Additionally, we show the predicted underlying benign and malignant prevalence of nodules in the population, the probability of detection based on size of nodule, and estimates of growth over time in both benign and malignant nodules.ConclusionsThis comprehensive model provides a dynamic platform employable for future comparative effectiveness research. Future model analyses will test and assess various clinical management strategies to improve patient outcomes related to thyroid cancer and optimize resource utilization for patients with thyroid nodules.