Mathematical modeling identifies LAG3 and HAVCR2 as biomarkers of T cell exhaustion in melanoma.

Cytotoxic T lymphocytes (CTLs) control tumors via lysis of antigen-presenting targets or through secretion of cytokines such as interferon-γ (IFNG), which inhibit tumor cell proliferation. Improved understanding of CTL interactions within solid tumors will aid the development of immunotherapeutic strategies against cancer. In this study, we take a systems biology approach to compare the importance of cytolytic versus IFNG-mediated cytostatic effects in a murine melanoma model (B16F10) and to dissect the contribution of immune checkpoints HAVCR2, LAG3, and PDCD1/CD274 to CTL exhaustion. We integrated multimodal data to inform an ordinary differential equation (ODE) model of CTL activities inside the tumor. Our model predicted that CTL cytotoxicity played only a minor role in tumor control relative to the cytostatic effects of IFNG. Furthermore, our analysis revealed that within B16F10 melanomas HAVCR2 and LAG3 better characterize the development of a dysfunctional CTL phenotype than does the PDCD1/CD274 axis.

Intra-tumor heterogeneity, turnover rate and karyotype space shape susceptibility to missegregation-induced extinction.

The phenotypic efficacy of somatic copy number alterations (SCNAs) stems from their incidence per base pair of the genome, which is orders of magnitudes greater than that of point mutations. One mitotic event stands out in its potential to significantly change a cell's SCNA burden-a chromosome missegregation. A stochastic model of chromosome mis-segregations has been previously developed to describe the evolution of SCNAs of a single chromosome type. Building upon this work, we derive a general deterministic framework for modeling missegregations of multiple chromosome types. The framework offers flexibility to model intra-tumor heterogeneity in the SCNAs of all chromosomes, as well as in missegregation- and turnover rates. The model can be used to test how selection acts upon coexisting karyotypes over hundreds of generations. We use the model to calculate missegregation-induced population extinction (MIE) curves, that separate viable from non-viable populations as a function of their turnover- and missegregation rates. Turnover- and missegregation rates estimated from scRNA-seq data are then compared to theoretical predictions. We find convergence of theoretical and empirical results in both the location of MIE curves and the necessary conditions for MIE. When a dependency of missegregation rate on karyotype is introduced, karyotypes associated with low missegregation rates act as a stabilizing refuge, rendering MIE impossible unless turnover rates are exceedingly high. Intra-tumor heterogeneity, including heterogeneity in missegregation rates, increases as tumors progress, rendering MIE unlikely.

Mathematical Modelling Based on In Vivo Imaging Suggests CD137-Stimulated Cytotoxic T Lymphocytes Exert Superior Tumour Control Due to an Enhanced Antimitotic Effect on Tumour Cells.

Several immunotherapeutic strategies for the treatment of cancer are under development. Two prominent strategies are adoptive cell transfer (ACT) of CTLs and modulation of CTL function with immune checkpoint inhibitors or with costimulatory antibodies. Despite some success with these approaches, there remains a lack of detailed and quantitative descriptions of the events following CTL transfer and the impact of immunomodulation. Here, we have applied ordinary differential equation models to two photon imaging data derived from a B16F10 murine melanoma. Models were parameterised with data from two different treatment conditions: either ACT-only, or ACT with intratumoural costimulation using a CD137 targeted antibody. Model dynamics and best fitting parameters were compared, in order to assess the mode of action of the CTLs and examine how the CD137 antibody influenced their activities. We found that the cytolytic activity of the transferred CTLs was minimal without CD137 costimulation, and that the CD137 targeted antibody did not enhance the per-capita killing ability of the transferred CTLs. Instead, the results of our modelling study suggest that an antiproliferative effect of CTLs exerted upon the tumour likely accounted for the majority of the reduction in tumour growth after CTL transfer. Moreover, we found that CD137 most likely improved tumour control via enhancement of this antiproliferative effect, as well as prolonging the period in which CTLs were inside the tumour, leading to a sustained duration of their antitumour effects following CD137 stimulation.

Heterogeneous, delayed-onset killing by multiple-hitting T cells: Stochastic simulations to assess methods for analysis of imaging data.

Although quantitative insights into the killing behaviour of Cytotoxic T Lymphocytes (CTLs) are necessary for the rational design of immune-based therapies, CTL killing function remains insufficiently characterised. One established model of CTL killing treats CTL cytotoxicity as a Poisson process, based on the assumption that CTLs serially kill antigen-presenting target cells via delivery of lethal hits, each lethal hit corresponding to a single injection of cytotoxic proteins into the target cell cytoplasm. Contradicting this model, a recent in vitro study of individual CTLs killing targets over a 12-hour period found significantly greater heterogeneity in CTL killing performance than predicted by Poisson-based killing. The observed killing process was dynamic and varied between CTLs, with the best performing CTLs exhibiting a marked increase in killing during the final hours of the experiments, along with a "burst killing" kinetic. Despite a search for potential differences between CTLs, no mechanistic explanation for the heterogeneous killing kinetics was found. Here we have used stochastic simulations to assess whether target cells might require multiple hits from CTLs before undergoing apoptosis, in order to verify whether multiple-hitting could explain the late onset, burst killing dynamics observed in vitro. We found that multiple-hitting from CTLs was entirely consistent with the observed killing kinetics. Moreover, the number of available targets and the spatiotemporal kinetics of CTL:target interactions influenced the realised CTL killing rate. We subsequently used realistic, spatial simulations to assess methods for estimating the hitting rate and the number of hits required for target death, to be applied to microscopy data of individual CTLs killing targets. We found that measuring the cumulative duration of individual contacts that targets have with CTLs would substantially improve accuracy when estimating the killing kinetics of CTLs.

Contact-Dependent Killing by Cytotoxic T Lymphocytes Is Insufficient for EL4 Tumor Regression In Vivo.

Immunotherapies are an emerging strategy for treatment of solid tumors. Improved understanding of the mechanisms employed by cytotoxic T lymphocytes (CTL) to control tumors will aid in the development of immunotherapies. CTLs can directly kill tumor cells in a contact-dependent manner or may exert indirect effects on tumor cells via secretion of cytokines. Here, we aim to quantify the importance of these mechanisms in murine thymoma EL4/EG7 cells. We developed an agent-based model (ABM) and an ordinary differential equation model of tumor regression after adoptive transfer of a population of CTLs. Models were parameterized based on in vivo measurements of CTL infiltration and killing rates applied to EL4/EG7 tumors and OTI T cells. We quantified whether infiltrating CTLs are capable of controlling tumors through only direct, contact-dependent killing. Both models agreed that the low measured killing rate of CTLs in vivo was insufficient to cause tumor regression. In our ABM, we also simulated CTL production of the cytokine IFNγ in order to explore how an antiproliferative effect of IFNγ might aid CTLs in tumor control. In this model, IFNγ substantially reduced tumor growth compared with direct killing alone. Collectively, these data demonstrate that contact-dependent killing is insufficient for EL4 regression in vivo and highlight the potential importance of cytokine-induced antiproliferative effects in T-cell-mediated tumor control. SIGNIFICANCE: Computational modeling highlights the importance of cytokine-induced antiproliferative effects in T-cell-mediated control of tumor progression.

Single-cell imaging of CAR T cell activity in vivo reveals extensive functional and anatomical heterogeneity.

CAR T cells represent a potentially curative strategy for B cell malignancies. However, the outcome and dynamics of CAR T cell interactions in distinct anatomical sites are poorly understood. Using intravital imaging, we tracked interactions established by anti-CD19 CAR T cells in B cell lymphoma-bearing mice. Circulating targets trapped CAR T cells in the lungs, reducing their access to lymphoid organs. In the bone marrow, tumor apoptosis was largely due to CAR T cells that engaged, killed, and detached from their targets within 25 min. Notably, not all CAR T cell contacts elicited calcium signaling or killing while interacting with tumors, uncovering extensive functional heterogeneity. Mathematical modeling revealed that direct killing was sufficient for tumor regression. Finally, antigen-loss variants emerged in the bone marrow, but not in lymph nodes, where CAR T cell cytotoxic activity was reduced. Our results identify a previously unappreciated level of diversity in the outcomes of CAR T cell interactions in vivo, with important clinical implications.

Simulation of lower limb axial arterial length change during locomotion.

The effect of external forces on axial arterial wall mechanics has conventionally been regarded as secondary to hemodynamic influences. However, arteries are similar to muscles in terms of the manner in which they traverse joints, and their three-dimensional geometrical requirements for joint motion. This study considers axial arterial shortening and elongation due to motion of the lower extremity during gait, ascending stairs, and sitting-to-standing motion. Arterial length change was simulated by means of a graphics based anatomic and kinematic model of the lower extremity. This model estimated the axial shortening to be as much as 23% for the femoropopliteal arterial region and as much as 21% for the iliac artery. A strong correlation was observed between femoropopliteal artery shortening and maximum knee flexion angle (r²=0.8) as well as iliac artery shortening and maximum hip angle flexion (r²=0.9). This implies a significant mechanical influence of locomotion on arterial behavior in addition to hemodynamics factors. Vascular tissue has high demands for axial compliance that should be considered in the pathology of atherosclerosis and the design of vascular implants.