Molecular and cellular effects of radiotherapy on tumor microenvironment (TME) can help prime and propagate antitumor immunity. We hypothesized that delivering radiation to all tumor sites could augment response to immunotherapies. We tested an approach to enhance response to immune checkpoint inhibitors (ICIs) by using targeted radionuclide therapy (TRT) to deliver radiation semiselectively to tumors. NM600, an alkylphosphocholine analog that preferentially accumulates in most tumor types, chelates a radioisotope and semiselectively delivers it to the TME for therapeutic or diagnostic applications. Using serial ⁸⁶Y-NM600 positron emission tomography (PET) imaging, we estimated the dosimetry of ⁹⁰Y-NM600 in immunologically cold syngeneic murine models that do not respond to ICIs alone. We observed strong therapeutic efficacy and reported optimal dose (2.5 to 5 gray) and sequence for ⁹⁰Y-NM600 in combination with ICIs. After combined treatment, 45 to 66% of mice exhibited complete response and tumor-specific T cell memory, compared to 0% with ⁹⁰Y-NM600 or ICI alone. This required expression of STING in tumor cells. Combined TRT and ICI activated production of proinflammatory cytokines in the TME, promoted tumor infiltration by and clonal expansion of CD8⁺ T cells, and reduced metastases. In mice bearing multiple tumors, combining TRT with moderate-dose (12 gray) external beam radiotherapy (EBRT) targeting a single tumor augmented response to ICIs compared to combination of ICIs with either TRT or EBRT alone. The safety of TRT was confirmed in a companion canine study. Low-dose TRT represents a translatable approach to promote response to ICIs for many tumor types, regardless of location.

Uncertainties in the estimated mean excitation energies (I-values) needed for calculating proton stopping powers can be in the order of 10-15%, which introduces a fundamental limitation in the accuracy of proton range determination. Previous efforts have quantified shifts in proton depth dose distributions due to I-value uncertainties in water and homogenous tissue phantoms. This study is the first to quantify the clinical impact of I-value uncertainties on proton dose distributions within patient geometries. A previously developed Geant4 based Monte Carlo code was used to simulate a proton treatment plan for three patients (prostate, pancreases, and liver) with varying tissue I-values. A uniform variation study was conducted in which the tissue I-values were varied by ±5% and ±10% of the nominal values as well as a probabilistic variation study in which the I-values were randomly sampled according to a normal distribution with the mean equal to the nominal I-value and a standard deviation of 5 and 10% of the nominal values. Modification of tissue I-values impacted both the proton range and SOBP width. R(90) range shifts up to 7.7 mm (4.4.%) and R(80) range shifts up to 4.8 mm (1.9%) from the nominal range were recorded. Modulating the tissue I-values by 10% the nominal value resulted in up to a 3.5% difference mean dose in the target volumes and organs at risk compared to the nominal case. The range and dose differences were the largest for the deeper-seated prostate and pancreas cases. The treatments that were simulated with randomly sampled I-values resulted in range and dose differences that were generally within the upper and lower bounds set by the 10% uniform variations. This study demonstrated the impact of I-value uncertainties on patient dose distributions. Clearly, sub-millimeter precision in proton therapy would necessitate a reduction in I-value uncertainties to ensure an efficacious clinical outcome.

Nowadays, Monte Carlo models of proton therapy treatment heads are being used to improve beam delivery systems and to calculate the radiation field for patient dose calculations. The achievable accuracy of the model depends on the exact knowledge of the treatment head geometry and time structure, the material characteristics, and the underlying physics. This work aimed at studying the uncertainties in treatment head simulations for passive scattering proton therapy. The sensitivities of spread-out Bragg peak (SOBP) dose distributions on material densities, mean ionization potentials, initial proton beam energy spread and spot size were investigated. An improved understanding of the nature of these parameters may help to improve agreement between calculated and measured SOBP dose distributions and to ensure that the range, modulation width, and uniformity are within clinical tolerance levels. Furthermore, we present a method to make small corrections to the uniformity of spread-out Bragg peaks by utilizing the time structure of the beam delivery. In addition, we re-commissioned the models of the two proton treatment heads located at our facility using the aforementioned correction methods presented in this paper.