A kinetic model of continuous radiation damage to populations of cells: comparison to the LQ model and application to molecular radiotherapy

Designing combination therapies for cancer treatment: application of a mathematical framework combining CAR T-cell immunotherapy and targeted radionuclide therapy

Introduction

Cancer combination treatments involving immunotherapies with targeted radiation therapy are at the forefront of treating cancers. However, dosing and scheduling of these therapies pose a challenge. Mathematical models provide a unique way of optimizing these therapies.

Methods

Using a preclinical model of multiple myeloma as an example, we demonstrate the capability of a mathematical model to combine these therapies to achieve maximum response, defined as delay in tumor growth. Data from mice studies with targeted radionuclide therapy (TRT) and chimeric antigen receptor (CAR)-T cell monotherapies and combinations with different intervals between them was used to calibrate mathematical model parameters. The dependence of progression-free survival (PFS), overall survival (OS), and the time to minimum tumor burden on dosing and scheduling was evaluated. Different dosing and scheduling schemes were evaluated to maximize the PFS and optimize timings of TRT and CAR-T cell therapies.

Results

Therapy intervals that were too close or too far apart are shown to be detrimental to the therapeutic efficacy, as TRT too close to CAR-T cell therapy results in radiation related CAR-T cell killing while the therapies being too far apart result in tumor regrowth, negatively impacting tumor control and survival. We show that splitting a dose of TRT or CAR-T cells when administered in combination is advantageous only if the first therapy delivered can produce a significant benefit as a monotherapy.

Discussion

Mathematical models are crucial tools for optimizing the delivery of cancer combination therapy regimens with application along the lines of achieving cure, maximizing survival or minimizing toxicity.

Primary study of the relative and compound biological effectiveness model for boron neutron capture therapy based on nanodosimetry

BACKGROUND

The current radiobiological model employed for boron neutron capture therapy (BNCT) treatment planning, which relies on microdosimetry, fails to provide an accurate representation the biological effects of BNCT. The precision in calculating the relative biological effectiveness (RBE) and compound biological effectiveness (CBE) plays a pivotal role in determining the therapeutic efficacy of BNCT. Therefore, this study focuses on how to improve the accuracy of the biological effects of BNCT.

PURPOSE

The purpose of this study is to propose new radiation biology models based on nanodosimetry to accurately assess RBE and CBE for BNCT.

METHODS

Nanodosimetry, rooted in ionization cluster size distributions (ICSD), introduces a novel approach to characterize radiation quality by effectively delineating RBE through the ion track structure at the nanoscale. In the context of prior research, this study presents a computational model for the nanoscale assessment of RBE and CBE. We establish a simplified model of DNA chromatin fiber using the Monte Carlo code TOPAS-nBio to evaluate the applicability of ICSD to BNCT and compute nanodosimetric parameters.

RESULTS

Our investigation reveals that both homogeneous and heterogeneous nanodosimetric parameters, as well as the corresponding biological model coefficients α and β, along with RBE values, exhibit variations in response to varying intracellular 10B concentrations. Notably, the nanodosimetric parameterM 1 C 2 $M_1^{{{\mathrm{C}}}_2}$ effectively captures the fluctuations in model coefficients α and RBE.

CONCLUSION

Our model facilitates a nanoscale analysis of BNCT, enabling predictions of nanodosimetric quantities for secondary ions as well as RBE, CBE, and other essential biological metrics related to the distribution of boron. This contribution significantly enhances the precision of RBE calculations and holds substantial promise for future applications in treatment planning.

Investigating the Influence of Growth Arrest Mechanisms on Tumour Responses to Radiotherapy

Cancer is a heterogeneous disease and tumours of the same type can differ greatly at the genetic and phenotypic levels. Understanding how these differences impact sensitivity to treatment is an essential step towards patient-specific treatment design. In this paper, we investigate how two different mechanisms for growth control may affect tumour cell responses to fractionated radiotherapy (RT) by extending an existing ordinary differential equation model of tumour growth. In the absence of treatment, this model distinguishes between growth arrest due to nutrient insufficiency and competition for space and exhibits three growth regimes: nutrient limited, space limited (SL) and bistable (BS), where both mechanisms for growth arrest coexist. We study the effect of RT for tumours in each regime, finding that tumours in the SL regime typically respond best to RT, while tumours in the BS regime typically respond worst to RT. For tumours in each regime, we also identify the biological processes that may explain positive and negative treatment outcomes and the dosing regimen which maximises the reduction in tumour burden.

A Mathematical Modeling Approach for Targeted Radionuclide and Chimeric Antigen Receptor T Cell Combination Therapy

Simple Summary

Targeted radionuclide therapy (TRT) and immunotherapy, an example being chimeric antigen receptor T cells (CAR-Ts), represent two potent means of eradicating systemic cancers. Although each one as a monotherapy might have a limited effect, the potency can be increased with a combination of the two therapies. The complications involved in the dosing and scheduling of these therapies make the mathematical modeling of these therapies a suitable solution for designing combination treatment approaches. Here, we investigate a mathematical model for TRT and CAR-T cell combination therapies. Through an analysis of the mathematical model, we find that the tumor proliferation rate is the most important factor affecting the scheduling of TRT and CAR-T cell treatments with faster proliferating tumors requiring a shorter interval between the two therapies.

Abstract

Targeted radionuclide therapy (TRT) has recently seen a surge in popularity with the use of radionuclides conjugated to small molecules and antibodies. Similarly, immunotherapy also has shown promising results, an example being chimeric antigen receptor T cell (CAR-T) therapy in hematologic malignancies. Moreover, TRT and CAR-T therapies possess unique features that require special consideration when determining how to dose as well as the timing and sequence of combination treatments including the distribution of the TRT dose in the body, the decay rate of the radionuclide, and the proliferation and persistence of the CAR-T cells. These characteristics complicate the additive or synergistic effects of combination therapies and warrant a mathematical treatment that includes these dynamics in relation to the proliferation and clearance rates of the target tumor cells. Here, we combine two previously published mathematical models to explore the effects of dose, timing, and sequencing of TRT and CAR-T cell-based therapies in a multiple myeloma setting. We find that, for a fixed TRT and CAR-T cell dose, the tumor proliferation rate is the most important parameter in determining the best timing of TRT and CAR-T therapies.

Classification of tolerable/intolerable mucosal toxicity of head-and-neck radiotherapy schedules with a biomathematical model of cell dynamics

Purpose

The purpose of this study is to present a biomathematical model based on the dynamics of cell populations to predict the tolerability/intolerability of mucosal toxicity in head‐and‐neck radiotherapy.

Methods and Materials

Our model is based on the dynamics of proliferative and functional cell populations in irradiated mucosa, and incorporates the three As: Accelerated proliferation, loss of Asymmetric proliferation, and Abortive divisions. The model consists of a set of delay differential equations, and tolerability is based on the depletion of functional cells during treatment. We calculate the sensitivity (sen) and specificity (spe) of the model in a dataset of 108 radiotherapy schedules, and compare the results with those obtained with three phenomenological classification models, two based on a biologically effective dose (BED) function describing the tolerability boundary (Fowler and Fenwick) and one based on an equivalent dose in 2 Gy fractions (EQD₂) boundary (Strigari). We also perform a machine learning‐like cross‐validation of all the models, splitting the database in two, one for training and one for validation.

Results

When fitting our model to the whole dataset, we obtain predictive values (sen + spe) up to 1.824. The predictive value of our model is very similar to that of the phenomenological models of Fowler (1.785), Fenwick (1.806), and Strigari (1.774). When performing a k = 2 cross‐validation, the specificity and sensitivity in the validation dataset decrease for all models, from ˜1.82 to ˜1.55–1.63. For Fowler, the worsening is higher, down to 1.49.

Conclusions

Our model has proved useful to predict the tolerability/intolerability of a dataset of 108 schedules. As the model is more mechanistic than other available models, it could prove helpful when designing unconventional dose fractionations, schedules not covered by datasets to which phenomenological models of toxicity have been fitted.