Early prediction of clinical response to checkpoint inhibitor therapy in human solid tumors through mathematical modeling

Immune-checkpoint inhibitor therapy response evaluation using oncophysics-based mathematical models

The field of oncology has transformed with the advent of immunotherapies. The standard of care for multiple cancers now includes novel drugs that target key checkpoints that function to modulate immune responses, enabling the patient's immune system to elicit an effective anti-tumor response. While these immune-based approaches can have dramatic effects in terms of significantly reducing tumor burden and prolonging survival for patients, the therapeutic approach remains active only in a minority of patients and is often not durable. Multiple biological investigations have identified key markers that predict response to the most common form of immunotherapy-immune checkpoint inhibitors (ICI). These biomarkers help enrich patients for ICI but are not 100% predictive. Understanding the complex interactions of these biomarkers with other pathways and factors that lead to ICI resistance remains a major goal. Principles of oncophysics-the idea that cancer can be described as a multiscale physical aberration-have shown promise in recent years in terms of capturing the essence of the complexities of ICI interactions. Here, we review the biological knowledge of mechanisms of ICI action and how these are incorporated into modern oncophysics-based mathematical models. Building on the success of oncophysics-based mathematical models may help to discover new, rational methods to engineer immunotherapy for patients in the future. This article is categorized under: Therapeutic Approaches and Drug Discovery > Nanomedicine for Oncologic Disease.

SMoRe ParS: A novel methodology for bridging modeling modalities and experimental data applied to 3D vascular tumor growth

Multiscale systems biology is having an increasingly powerful impact on our understanding of the interconnected molecular, cellular, and microenvironmental drivers of tumor growth and the effects of novel drugs and drug combinations for cancer therapy. Agent-based models (ABMs) that treat cells as autonomous decision-makers, each with their own intrinsic characteristics, are a natural platform for capturing intratumoral heterogeneity. Agent-based models are also useful for integrating the multiple time and spatial scales associated with vascular tumor growth and response to treatment. Despite all their benefits, the computational costs of solving agent-based models escalate and become prohibitive when simulating millions of cells, making parameter exploration and model parameterization from experimental data very challenging. Moreover, such data are typically limited, coarse-grained and may lack any spatial resolution, compounding these challenges. We address these issues by developing a first-of-its-kind method that leverages explicitly formulated surrogate models (SMs) to bridge the current computational divide between agent-based models and experimental data. In our approach, Surrogate Modeling for Reconstructing Parameter Surfaces (SMoRe ParS), we quantify the uncertainty in the relationship between agent-based model inputs and surrogate model parameters, and between surrogate model parameters and experimental data. In this way, surrogate model parameters serve as intermediaries between agent-based model input and data, making it possible to use them for calibration and uncertainty quantification of agent-based model parameters that map directly onto an experimental data set. We illustrate the functionality and novelty of Surrogate Modeling for Reconstructing Parameter Surfaces by applying it to an agent-based model of 3D vascular tumor growth, and experimental data in the form of tumor volume time-courses. Our method is broadly applicable to situations where preserving underlying mechanistic information is of interest, and where computational complexity and sparse, noisy calibration data hinder model parameterization.

Mathematical modeling of cancer immunotherapy for personalized clinical translation

Encouraging advances are being made in cancer immunotherapy modeling, especially in the key areas of developing personalized treatment strategies based on individual patient parameters, predicting treatment outcomes and optimizing immunotherapy synergy when used in combination with other treatment approaches. Here we present a focused review of the most recent mathematical modeling work on cancer immunotherapy with a focus on clinical translatability. It can be seen that this field is transitioning from pure basic science to applications that can make impactful differences in patients' lives. We discuss how researchers are integrating experimental and clinical data to fully inform models so that they can be applied for clinical predictions, and present the challenges that remain to be overcome if widespread clinical adaptation is to be realized. Lastly, we discuss the most promising future applications and areas that are expected to be the focus of extensive upcoming modeling studies.

Protocol for mathematical prediction of patient response and survival to immune checkpoint inhibitor immunotherapy

This protocol describes the application of a mechanistic mathematical model of immune checkpoint inhibitor (ICI) immunotherapy to patient tumor imaging data for predicting solid tumor response and patient survival under ICI intervention. We describe steps for data collection and processing, data pipelines, and approaches to increase precision. The protocol is highly predictive as early as the first restaging after treatment start and can be used with standard-of-care imaging measures. For complete details on the use and execution of this protocol, please refer to Butner et al. (2020)¹ and Butner et al. (2021).².

Distinct Dynamics of Migratory Response to PD-1 and CTLA-4 Blockade Reveals New Mechanistic Insights for Potential T-Cell Reinvigoration following Immune Checkpoint Blockade

Cytotoxic T-lymphocyte-associated antigen 4 (CTLA-4) and programmed cell death protein 1 (PD-1), two clinically relevant targets for the immunotherapy of cancer, are negative regulators of T-cell activation and migration. Optimizing the therapeutic response to CTLA-4 and PD-1 blockade calls for a more comprehensive insight into the coordinated function of these immune regulators. Mathematical modeling can be used to elucidate nonlinear tumor-immune interactions and highlight the underlying mechanisms to tackle the problem. Here, we investigated and statistically characterized the dynamics of T-cell migration as a measure of the functional response to these pathways. We used a previously developed three-dimensional organotypic culture of patient-derived tumor spheroids treated with anti-CTLA-4 and anti-PD-1 antibodies for this purpose. Experiment-based dynamical modeling revealed the delayed kinetics of PD-1 activation, which originates from the distinct characteristics of PD-1 and CTLA-4 regulation, and followed through with the modification of their contributions to immune modulation. The simulation results show good agreement with the tumor cell reduction and active immune cell count in each experiment. Our findings demonstrate that while PD-1 activation provokes a more exhaustive intracellular cascade within a mature tumor environment, the time-delayed kinetics of PD-1 activation outweighs its preeminence at the individual cell level and consequently confers a functional dominance to the CTLA-4 checkpoint. The proposed model explains the distinct immunostimulatory pattern of PD-1 and CTLA-4 blockade based on mechanisms involved in the regulation of their expression and may be useful for planning effective treatment schemes targeting PD-1 and CTLA-4 functions.