Integrating local and distant radiation‐induced lung injury: Development and validation of a predictive model for ventilation loss
Authors: Vicente Esther M. , Grande Gutierrez Noelia , Oakes Jessica M. , Cammin Jochen , Gopal Arun , Kipritidis John , Modiri Arezoo , Mossahebi Sina , Mohindra Pranshu , Citron Wendla K. , Matuszak Martha M. , Timmerman Robert , Sawant Amit 
Published: May 2024

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.