Wang, Yixuan, Bergman, Daniel R., Trujillo, Erica, Fernald, Anthony A., Li, Lie, Pearson, Alexander T., Sweis, Randy F., and Jackson, Trachette L.
Simple Summary: Immune checkpoint inhibitors (ICIs) are cancer immunotherapeutics that reinvigorate immune cells' ability to attack tumor cells. Despite remarkable results in some patients, ICIs do not demonstrate the same efficacy across all individuals. In this study, we present the first side-by-side comparison of an agent-based model (ABM) with an ordinary differential equation (ODE) model for ICIs targeting the PD-1/PD-L1 immune checkpoint. We consider tumor cells of high and low antigenicity and two distinct immune-cell kill mechanisms. Using key parameters calibrated from mouse bladder cancer studies, we simulate virtual tumors using both models. Our research identifies crucial tumor-immune characteristics that influence the efficacy of ICIs. By exploring the unique spatial insights provided by the ABM, we underscore the importance of considering the spatial complexity of the tumor microenvironment in mathematical models of ICIs, potentially paving the way for more effective cancer treatments. Since the introduction of the first immune checkpoint inhibitor (ICI), immunotherapy has changed the landscape of molecular therapeutics for cancers. However, ICIs do not work equally well on all cancers and for all patients. There has been a growing interest in using mathematical and computational models to optimize clinical responses. Ordinary differential equations (ODEs) have been widely used for mechanistic modeling in immuno-oncology and immunotherapy. They allow rapid simulations of temporal changes in the cellular and molecular populations involved. Nonetheless, ODEs cannot describe the spatial structure in the tumor microenvironment or quantify the influence of spatially-dependent characteristics of tumor-immune dynamics. For these reasons, agent-based models (ABMs) have gained popularity because they can model more detailed phenotypic and spatial heterogeneity that better reflect the complexity seen in vivo. In the context of anti-PD-1 ICIs, we compare treatment outcomes simulated from an ODE model and an ABM to show the importance of including spatial components in computational models of cancer immunotherapy. We consider tumor cells of high and low antigenicity and two distinct cytotoxic T lymphocyte (CTL) killing mechanisms. The preferred mechanism differs based on the antigenicity of tumor cells. Our ABM reveals varied phenotypic shifts within the tumor and spatial organization of tumor and CTLs despite similarities in key immune parameters, initial simulation conditions, and early temporal trajectories of the cell populations. [ABSTRACT FROM AUTHOR]