Modeling Patient-Specific Dose-Function Response for Enhanced Characterization of Personalized Functional Damage

Integrating local and distant radiation-induced lung injury: Development and validation of a predictive model for ventilation loss

BACKGROUND

Investigations on radiation-induced lung injury (RILI) have predominantly focused on local effects, primarily those associated with radiation damage to lung parenchyma. However, recent studies from our group and others have revealed that radiation-induced damage to branching serial structures such as airways and vessels may also have a substantial impact on post-radiotherapy (RT) lung function. Furthermore, recent results from multiple functional lung avoidance RT trials, although promising, have demonstrated only modest toxicity reduction, likely because they were primarily focused on dose avoidance to lung parenchyma. These observations emphasize the critical need for predictive dose-response models that effectively incorporate both local and distant RILI effects.

PURPOSE

We develop and validate a predictive model for ventilation loss after lung RT. This model, referred to as P+A, integrates local (parenchyma [P]) and distant (central and peripheral airways [A]) radiation-induced damage, modeling partial (narrowing) and complete (collapse) obstruction of airways.

METHODS

In an IRB-approved prospective study, pre-RT breath-hold CTs (BHCTs) and pre- and one-year post-RT 4DCTs were acquired from lung cancer patients treated with definitive RT. Up to 13 generations of airways were automatically segmented on the BHCTs using a research virtual bronchoscopy software. Ventilation maps derived from the 4DCT scans were utilized to quantify pre- and post-RT ventilation, serving, respectively, as input data and reference standard (RS) in model validation. To predict ventilation loss solely due to parenchymal damage (referred to as P model), we used a normal tissue complication probability (NTCP) model. Our model used this NTCP-based estimate and predicted additional loss due radiation-induced partial or complete occlusion of individual airways, applying fluid dynamics principles and a refined version of our previously developed airway radiosensitivity model. Predictions of post-RT ventilation were estimated in the sublobar volumes (SLVs) connected to the terminal airways. To validate the model, we conducted a k-fold cross-validation. Model parameters were optimized as the values that provided the lowest root mean square error (RMSE) between predicted post-RT ventilation and the RS for all SLVs in the training data. The performance of the P+A and the P models was evaluated by comparing their respective post-RT ventilation values with the RS predictions. Additional evaluation using various receiver operating characteristic (ROC) metrics was also performed.

RESULTS

We extracted a dataset of 560 SLVs from four enrolled patients. Our results demonstrated that the P+A model consistently outperformed the P model, exhibiting RMSEs that were nearly half as low across all patients (13 ± 3 percentile for the P+A model vs. 24 ± 3 percentile for the P model on average). Notably, the P+A model aligned closely with the RS in ventilation loss distributions per lobe, particularly in regions exposed to doses ≥13.5 Gy. The ROC analysis further supported the superior performance of the P+A model compared to the P model in sensitivity (0.98 vs. 0.07), accuracy (0.87 vs. 0.25), and balanced predictions.

CONCLUSIONS

These early findings indicate that airway damage is a crucial factor in RILI that should be included in dose-response modeling to enhance predictions of post-RT lung function.

Evaluation of core-shell Fe3O4@Au nanoparticles as radioenhancer in A549 cell lung cancer model

In radiotherapy, metallic nanoparticles are of high interest in the fight against cancer for their radiosensitizing effects. This study aimed to evaluate the ability of core-shell Fe3O4@Au nanoparticles to potentiate the irradiation effects on redox-, pro-inflammatory markers, and cell death of A549 human pulmonary cancer cells. The hybrid Fe3O4@Au nanoparticles were synthesized using green chemistry principles by the sonochemistry method. Their characterization by transmission electron microscopy demonstrated an average size of 8 nm and a homogeneous distribution of gold. The decreased hydrodynamic size of these hybrid nanoparticles compared to magnetite (Fe3O4) nanoparticles showed that gold coating significantly reduced the aggregation of Fe3O4 particles. The internalization and accumulation of the Fe3O4@Au nanoparticles within the cells were demonstrated by Prussian Blue staining. The reactive oxygen species (ROS) levels measured by the fluorescent probe DCFH-DA were up-regulated, as well as mRNA expression of SOD, catalase, GPx antioxidant enzymes, redox-dependent transcription factor Nrf2, and ROS-producing enzymes (Nox2 and Nox4), quantified by RT-qPCR. Furthermore, irradiation coupled with Fe3O4@Au nanoparticles increased the expression of canonical pro-inflammatory cytokines and chemokines (TNF-α, IL-1β, IL-6, CXCL8, and CCL5) assessed by RT-qPCR and ELISA. Hybrid nanoparticles did not potentiate the increased DNA damage detected by immunofluorescence following the irradiation. Nevertheless, Fe3O4@Au caused cellular damage, leading to apoptosis through activation of caspase 3/7, secondary necrosis quantified by LDH release, and cell growth arrest evaluated by clonogenic-like assay. This study demonstrated the potential of Fe3O4@Au nanoparticles to potentiate the radiosensitivity of cancerous cells.

Incorporation of Functional Lung Imaging Into Radiation Therapy Planning in Patients With Lung Cancer: A Systematic Review and Meta-Analysis

Our purpose was to provide an understanding of current functional lung imaging (FLI) techniques and their potential to improve dosimetry and outcomes for patients with lung cancer receiving radiation therapy (RT). Excerpta Medica dataBASE (EMBASE), PubMed, and Cochrane Library were searched from 1990 until April 2023. Articles were included if they reported on FLI in one of: techniques, incorporation into RT planning for lung cancer, or quantification of RT-related outcomes for patients with lung cancer. Studies involving all RT modalities, including stereotactic body RT and particle therapy, were included. Meta-analyses were conducted to investigate differences in dose-function parameters between anatomic and functional RT planning techniques, as well as to investigate correlations of dose-function parameters with grade 2+ radiation pneumonitis (RP). One hundred seventy-eight studies were included in the narrative synthesis. We report on FLI modalities, dose-response quantification, functional lung (FL) definitions, FL avoidance techniques, and correlations between FL irradiation and toxicity. Meta-analysis results show that FL avoidance planning gives statistically significant absolute reductions of 3.22% to the fraction of well-ventilated lung receiving 20 Gy or more, 3.52% to the fraction of well-perfused lung receiving 20 Gy or more, 1.3 Gy to the mean dose to the well-ventilated lung, and 2.41 Gy to the mean dose to the well-perfused lung. Increases in the threshold value for defining FL are associated with decreases in functional parameters. For intensity modulated RT and volumetric modulated arc therapy, avoidance planning results in a 13% rate of grade 2+ RP, which is reduced compared with results from conventional planning cohorts. A trend of increased predictive ability for grade 2+ RP was seen in models using FL information but was not statistically significant. FLI shows promise as a method to spare FL during thoracic RT, but interventional trials related to FL avoidance planning are sparse. Such trials are critical to understanding the effect of FL avoidance planning on toxicity reduction and patient outcomes.

Shrinkage priors for isotonic probability vectors and binary data modeling, with applications to dose-response modeling

Motivated by the need to model dose-response or dose-toxicity curves in clinical trials, we develop a new horseshoe-based prior for Bayesian isotonic regression modeling a binary outcome against an ordered categorical predictor, where the probability of the outcome is assumed to be monotonically non-decreasing with the predictor. The set of differences between outcome probabilities in consecutive categories of the predictor is equipped with a multivariate prior having support over simplex. The Dirichlet distribution, which can be derived from a normalized sum of independent gamma-distributed random variables, is a natural choice of prior, but using mathematical and simulation-based arguments, we show that the resulting posterior is prone to underflow and other numerical instabilities, even under simple data configurations. We propose an alternative prior based on horseshoe-type shrinkage that is numerically more stable. We show that this horseshoe-based prior is not subject to the numerical instability seen in the Dirichlet/gamma-based prior and that the horseshoe-based posterior can estimate the underlying true curve more efficiently than the Dirichlet-based one. We demonstrate the use of this prior in a model predicting the occurrence of radiation-induced lung toxicity in lung cancer patients as a function of dose delivered to normal lung tissue. Our methodology is implemented in the R package isotonicBayes and therefore suitable for use in the design of dose-finding studies or other dose-response modeling contexts.

Combining Serial and Parallel Functionality in Functional Lung Avoidance Radiation Therapy

PURPOSE

Functional lung avoidance (FLA) radiation therapy (RT) aims to minimize post-RT pulmonary toxicity by preferentially avoiding dose to high-functioning lung (HFL) regions. A common limitation is that FLA approaches do not consider the conducting architecture for gas exchange. We previously proposed the functionally weighted airway sparing (FWAS) method to spare airways connected to HFL regions, showing that it is possible to substantially reduce risk of radiation-induced airway injury. Here, we compare the performance of FLA and FWAS and propose a novel method combining both approaches.

METHODS

We used breath-hold computed tomography (BHCT) and simulation 4-dimensional computed tomography (4DCT) from 12 lung stereotactic ablative radiation therapy patients. Four planning strategies were examined: (1) Conventional: no sparing other than clinical dose-volume constraints; (2) FLA: using a 4DCT-based ventilation map to delineate the HFL, plans were optimized to reduce mean dose and V_13.50 in HFL; (3) FWAS: we autosegemented 11 to 13 generations of individual airways from each patient's BHCT and assigned priorities based on the relative contribution of each airway to total ventilation. We used these priorities in the optimization along with airway dose constraints, estimated as a function of airway diameter and 5% probability of collapse; and (4) FLA + FWAS: we combined information from the 2 strategies. We prioritized clinical dose constraints for organs at risk and planning target volume in all plans. We performed the evaluation in terms of ventilation preservation accounting for radiation-induced damage to both lung parenchyma and airways.

RESULTS

We observed average ventilation preservation for FLA, FWAS, and FLA + FWAS as 3%, 8.5%, and 14.5% higher, respectively, than for Conventional plans for patients with ventilation preservation in Conventional plans <90%. Generalized estimated equations showed that all improvements were statistically significant (P ≤ .036). We observed no clinically relevant improvements in outcomes of the sparing techniques in patients with ventilation preservation in Conventional plans ≥90%.

CONCLUSIONS

These initial results suggest that it is crucial to consider the parallel and the serial nature of the lung to improve post-radiation therapy lung function and, consequently, quality of life for patients.

Investigating the SPECT Dose-Function Metrics Associated With Radiation-Induced Lung Toxicity Risk in Patients With Non-small Cell Lung Cancer Undergoing Radiation Therapy

Purpose

Dose to normal lung has commonly been linked with radiation-induced lung toxicity (RILT) risk, but incorporating functional lung metrics in treatment planning may help further optimize dose delivery and reduce RILT incidence. The purpose of this study was to investigate the impact of the dose delivered to functional lung regions by analyzing perfusion (Q), ventilation (V), and combined V/Q single-photon-emission computed tomography (SPECT) dose-function metrics with regard to RILT risk in patients with non-small cell lung cancer (NSCLC) patients who received radiation therapy (RT).

Methods and Materials

SPECT images acquired from 88 patients with locally advanced NSCLC before undergoing conventionally fractionated RT were retrospectively analyzed. Dose was converted to the nominal dose equivalent per 2 Gy fraction, and SPECT intensities were normalized. Regional lung segments were defined, and the average dose delivered to each lung region was quantified. Three functional categorizations were defined to represent low-, normal-, and high-functioning lungs. The percent of functional lung category receiving ≥20 Gy and mean functional intensity receiving ≥20 Gy (iV₂₀) were calculated. RILT was defined as grade 2+ radiation pneumonitis and/or clinical radiation fibrosis. A logistic regression was used to evaluate the association between dose-function metrics and risk of RILT.

Results

By analyzing V/Q normalized intensities and functional distributions across the population, a wide range in functional capability (especially in the ipsilateral lung) was observed in patients with NSCLC before RT. Through multivariable regression models, global lung average dose to the lower lung was found to be significantly associated with RILT, and Q and V iV₂₀ were correlated with RILT when using ipsilateral lung metrics. Through a receiver operating characteristic analysis, combined V/Q low-function receiving ≥20 Gy (low-functioning V/Q₂₀) in the ipsilateral lung was found to be the best predictor (area under the curce: 0.79) of RILT risk.

Conclusions

Irradiation of the inferior lung appears to be a locational sensitivity for RILT risk. The multivariable correlation between ipsilateral lung iV₂₀ and RILT, as well as the association of low-functioning V/Q₂₀ and RILT, suggest that irradiating low-functioning regions in the lung may lead to higher toxicity rates.

Functionally weighted airway sparing (FWAS): a functional avoidance method for preserving post-treatment ventilation in lung radiotherapy

Recent changes to the guidelines for screening and early diagnosis of lung cancer have increased the interest in preserving post-radiotherapy lung function. Current investigational approaches are based on spatially mapping functional regions and generating regional avoidance plans that preferentially spare highly ventilated/perfused lung. A potentially critical, yet overlooked, aspect of functional avoidance is radiation injury to peripheral airways, which serve as gas conduits to and from functional lung regions. Dose redistribution based solely on regional function may cause irreparable damage to the 'supply chain'. To address this deficiency, we propose the functionally weighted airway sparing (FWAS) method. FWAS (i) maps the bronchial pathways to each functional sub-lobar lung volume; (ii) assigns a weighting factor to each airway based on the relative contribution of the sub-volume to overall lung function; and (iii) creates a treatment plan that aims to preserve these functional pathways. To evaluate it, we used four cases from a retrospective cohort of SAbR patients treated for lung cancer. Each patient's airways were auto-segmented from a diagnostic-quality breath-hold CT using a research virtual bronchoscopy software. A ventilation map was generated from the planning 4DCT to map regional lung function. For each terminal airway, as resolved by the segmentation software, the total ventilation within the sub-lobar volume supported by that airway was estimated and used as a function-based weighting factor. Upstream airways were weighted based on the cumulative volumetric ventilation supported by corresponding downstream airways. Using a previously developed model for airway radiosensitivity, dose constraints were determined for each airway corresponding to a <5% probability of airway collapse. Airway dose constraints, ventilation scores, and clinical dose constraints were input to a swarm optimization-based inverse planning engine to create a 3D conformal SAbR plan (CRT). The FWAS plans were compared to the patients' prescribed CRT clinical plans and the inverse-optimized clinical plans. Depending on the size and location of the tumour, the FWAS plan showed superior preservation of ventilation due to airflow preservation through open pathways (i.e. cumulative ventilation score from the sub-lobar volumes of open pathways). Improvements ranged between 3% and 23%, when comparing to the prescribed clinical plans, and between 3% and 35%, when comparing to the inverse-optimized clinical plans. The three plans satisfied clinical requirements for PTV coverage and OAR dose constraints. These initial results suggest that by sparing pathways to high-functioning lung subregions it is possible to reduce post-SAbR loss of respiratory function.

Comparison of regional lung perfusion response on longitudinal MAA SPECT/CT in lung cancer patients treated with and without functional tissue-avoidance radiation therapy

OBJECTIVE

The effect of functional lung avoidance planning on radiation dose-dependent changes in regional lung perfusion is unknown. We characterized dose-perfusion response on longitudinal perfusion single photon emission computed tomography (SPECT)/CT in two cohorts of lung cancer patients treated with and without functional lung avoidance techniques.

METHODS

The study included 28 primary lung cancer patients: 20 from interventional (NCT02773238) (FLARE-RT) and eight from observational (NCT01982123) (LUNG-RT) clinical trials. FLARE-RT treatment plans included perfused lung dose constraints while LUNG-RT plans adhered to clinical standards. Pre- and 3 month post-treatment macro-aggregated albumin (MAA) SPECT/CT scans were rigidly co-registered to planning four-dimensional CT scans. Tumour-subtracted lung dose was converted to EQD2 and sorted into 5 Gy bins. Mean dose and percent change between pre/post-RT MAA-SPECT uptake (%ΔPERF), normalized to total tumour-subtracted lung uptake, were calculated in each binned dose region. Perfusion frequency histograms of pre/post-RT MAA-SPECT were analyzed. Dose-response data were parameterized by sigmoid logistic functions to estimate maximum perfusion increase (%ΔPERFmaxincrease), maximum perfusion decrease (%ΔPERFmaxdecrease), dose midpoint (Dmid), and dose-response slope (k).

RESULTS

Differences in MAA perfusion frequency distribution shape between time points were observed in 11/20 (55%) FLARE-RT and 2/8 (25%) LUNG-RT patients (p < 0.05). FLARE-RT dose response was characterized by >10% perfusion increase in the 0-5 Gy dose bin for 8/20 patients (%ΔPERFmaxincrease = 10-40%), which was not observed in any LUNG-RT patients (p = 0.03). The dose midpoint Dmid at which relative perfusion declined by 50% trended higher in FLARE-RT compared to LUNG-RT cohorts (35 GyEQD2 vs 21 GyEQD2, p = 0.09), while the dose-response slope k was similar between FLARE-RT and LUNG-RT cohorts (3.1-3.2, p = 0.86).

CONCLUSION

Functional lung avoidance planning may promote increased post-treatment perfusion in low dose regions for select patients, though inter-patient variability remains high in unbalanced cohorts. These preliminary findings form testable hypotheses that warrant subsequent validation in larger cohorts within randomized or case-matched control investigations.

ADVANCES IN KNOWLEDGE

This novel preliminary study reports differences in dose-response relationships between patients receiving functional lung avoidance radiation therapy (FLARE-RT) and those receiving conventionally planned radiation therapy (LUNG-RT). Following further validation and testing of these effects in larger patient populations, individualized estimation of regional lung perfusion dose-response may help refine future risk-adaptive strategies to minimize lung function deficits and toxicity incidence.