Mitigation of the Interplay Effects of Combining 4D Robust With Layer Repainting Techniques in Proton-Based SBRT for Patients With Early-Stage Non-small Cell Lung Cancer

Abstract

Objective:

The objective of this study was to evaluate the interplay effects in proton-based stereotactic body radiotherapy (SBRT) using 4D robust optimization combined with iso-energy layer repainting techniques for non-small cell lung cancer (NSCLC).

Materials and Methods:

Twelve patients with early-stage NSCLC who underwent 4DCT were retrospectively selected. A robust CTV-based 4D plan was generated for each based on commercial Treatment planning system (TPS), considering patient setup errors, range uncertainties, and organ motion. The 4D static dose (4DSD) and 4D dynamic dose (4DDD) were calculated using a hybrid deformable algorithm and simulated proton delivery system. An index was developed to quantitatively evaluate the interplay effects. The interplay effects of the 4D robust plan and multiple iso-energy layers (3, 4, 5, 6, and 7) of the robust repainting 4D plan were calculated based on to select the optimal times for layer repainting.

Results:

Due to the interplay effects, the mean target values D₂ and D₅ increased to 1.28 and 1.01%, and the target values D₉₈ and D₉₅ decreased to 2.01 and 1.77%, respectively, for the 4D Robust SBRT plan. After multiple iso-energy repainting times, the interplay effects of the target values D₉₈ and D₉₅ tended to migrate, from 2.01 to 0.92% in target value D₉₈ and from 1.77 to 0.89% in target value D₉₅, respectively. Moreover, a positive linear correlation was observed between the optimal interplay effect mitigation and target range of motion.

Conclusion:

In proton-based 4D Robust SBRT, the interplay effects degraded the target dose distribution but were mitigated using iso-energy layer repainting techniques. However, there was no significant correlation between the number of repainting layers and improvements in the dose distributions. The optimal layer repainting times based on the interplay effect index were ascertained according to the patient characteristics.

Introduction

Proton spot scanning-based stereotactic body radiotherapy (SBRT) has been shown to outperform photon-based SBRT in patients with early-stage non-small cell lung cancer (NSCLC). Proton has a superior Bragg peak increasing the dose delivered to tumors and sparing healthy organs, such as the lungs, esophagus, and spinal cord (1, 2). However, patient setup errors and range and target motion uncertainties must be addressed to make full use of the advantages of proton-based SBRT for lung cancer (3).

The motion of lung tumors has been evaluated based on 4DCT imaging. Most lung tumors have limited motion, less than 5 mm, especially in locally advanced NSCLC (4). Patient setup and range uncertainties are the primary factors to consider during IMPT planning. In early-stage NSCLC, approximately 50% of lung tumors move more than 5 mm and some even move more than 2 cm in the superior–inferior direction (5), where interplay effects caused by the interference between the beam spot and intra-fractional respiratory motion is dominant and should be taken into account. Otherwise, the quality of the dose distribution can be severely degraded. The effects should be minimized as much as possible. Robust optimization combining with 4DCT imaging (4DRP) (6, 7) mitigates the interplay effects. The appropriate repainting strategy (8, 9) can also manage the interplay effects. However, there are few clinical reports on the effectiveness of combining 4DRP and repainting with proton-based SBRT for IMPT in patients with early-stage NSCLC. Moreover, an interplay index for quantitative evaluation is needed to assess the effectiveness of combining 4DRP and repainting. The interplay effects can be estimated by calculating the single-fraction 4D static dose (4DSD) and single-fraction 4D dynamic dose (4DDD) based on 4DCT images (10, 11). The 4DSD is calculated based on the assumption that the tumor moves in 4DCT images without considering the delivery system’s time dependence. The 4DDD is calculated by taking into account the delivery system’s time dependence. However, this method is not very intuitive and cannot provide changes in the target volume or organs at risk (OARs), such as target conformity, homogeneity, and OAR volume dosimetry.

In the current study, a 4D Robust plan was generated for selected patients. The interplay effects were quantitatively evaluated using a target index and OAR index. A multiple iso-energy layer repainting strategy was also used to further mitigate the interplay effects to explore the optimal mitigation outcomes.

Materials and Methods

Patient Characteristics

A total of 12 patients with early-stage (IA/IB) NSCLC were selected for this study, which was approved by the local institutional research review board. The clinical target volumes (CTVs) were contoured by the attending radiation oncologists at each 4DCT phase. The internal target volumes (ITVs) were created by encompassing the extent of 10 CTV motions in 10 4DCT phases. The patients’ information is summarized in Table 1. The inclusion criteria were as follows: (i) the tumor was small with no distant metastasis, each <5 cm in diameter, (ii) the patients’ body surface was more than 4 cm away from the tumor and did not require the use of a range shifter, and (iii) the patients underwent specialized respiratory training to maintain a stable respiratory cycle for approximately 3–5 s.

Target Range of Motion

The target range of motion was obtained by calculating the maximum deformation vector lengths (DVLs) in the target area (9, 12). The maximum inhale phase T₀ and maximum exhale phase T₅₀ were used for deformable image registration (DIR) (13) to obtain the DVL. The DIR algorithm was developed by RaySearch and performs well in lung applications (14). Voxels in the CTVs were selected and the target range of motion was calculated according to the DVL formula:

where x_T0,i−x_T50,i,y_T0,i−y_T50,ii,z_T0,i−z_T50,i are the components in voxel i of the deformation vector field between the T₀ and T₅₀ phase images in the 4DCT images.

4D Robust Treatment Planning

Spot Scanning SBRT

Pencil beam scanning proton plans were generated for the patients via RayStation (RaySearch Laboratories, Version 6.1 sp1, Stockholm, Sweden) using proton energies between 70 and 225 MeV with beam data from a typical pencil beam scanning dedicated nozzle manufactured by IBA (Ion Beam Applications S.A., Louvain-la-Neuve, Belgium). The spot sizes at the iso-center in air varied between 2.5 mm at 225 MeV and 6.8 mm at 70 MeV. At least greater than 95% of the CTV received a prescribed dose of 60 Gy [RBE] in five fractions. Two or three suitable oblique coplanar beams were used for the plan with the beam direction according to the tumor location. Then, a 4D Robust optimization algorithm (7) was used based on the 4DCT images considering the patient setup, range uncertainties, and target motion. Before optimization, the minimum and maximum spot weights were 0.02 and 4 MU, respectively.

Uncertainty Modeling

Proton-based SRBT is sensitive to the patient setup, range uncertainties, and organ motion, so all of the uncertainties should be considered in the model. Inter-fractional patient setup uncertainties were simulated by shifting the patient iso-center in the antero–posterior (A–P), superior–inferior (S–I), and right–left (R–L) directions by 5 mm, yielding six dose distributions and the corresponding influence matrices (the beamlet dose distributions per unit intensity). Range uncertainties were simulated by scaling the stopping power ratios by ±3.5% to generate two additional dose distributions and influence matrices corresponding to the minimum and maximum proton ranges, respectively. The organ motion uncertainty was considered using 4DCT images consisting of 10 respiratory cycle phases to generate 10 dose distributions with each respiratory cycle phase.

4D Robust Optimization

The 4D Robust optimization plans (4DRP) were generated by optimizing the CTV dose in 10 4DCT phases considering the modeling uncertainties. Robust optimization taking into account the set S of scenarios was implemented using minimax optimization (15). The objective function was

where X is the set of feasible variables (spot weights for spot scanning IMPT), d(x;s) is the dose distribution as a function of variable x, scenario s, and f_iis the i structure’s penalty function. The robust optimization objective were used in CTV with a minimum dose objective of 60 Gy [RBE] (weight 100) and a maximum dose objective of 60 Gy [RBE] (weight 60). A dose fall-off function from 60 Gy [RBE] to 10 Gy [RBE] over 1 cm (weight 10) was used to lower the dose to the normal tissue as much as possible. The equivalent uniform dose (EUD) (16) on the normal lung was approximately 5 Gy [RBE] with a dose volume effect parameter of 1 (A = 1). In the final dose calculation, the Monte Carlo dose engine was used with 0.5% statistical uncertainty and a 3 mm× 3 mm× 3 mm dose grid resolution.

Interplay Effect Calculation

An overview of the interplay effect process is shown in Figure 1. The black flowchart shows how the 4DSD was produced, while the blue flowchart shows how the 4DDD was produced.

FIGURE 1

The quantitative interplay effect index was represented by the differences in the DVH metrics between 4DSD and 4DDD over 4DSD for each region of interest (ROI).

where and are specific DVH metrics for one ROI in the 4DSD distribution or 4DDD distribution. represents the 4DSD at 95% of the CTV volume. If the value of was positive, the DVH metrics increased in the 4DDD distribution and vice versa.

We used the CTV and lung minus ITV (lung – ITV) as the ROIs to obtain the quantitative interplay effects. The DVH metrics in the CTV included the minimum target dose (D₉₈[cGy(RBE)]: dose at 98% of the target volume), prescription dose (D₉₅[cGy(RBE)]: dose at 95% of the target volume), and maximum dose (D₂[cGy(RBE)]: dose at 2% of the target volume). The target EUD (A = 10) metrics were also included for interplay evaluation. For lung – ITV, the DVH metrics of the lungs included V₅, V₂₀, and V₃₀, which were the percentage volume of lungs receiving 500 cGy [RBE], 2000 cGy [RBE], and 3000 cGy [RBE], respectively. The lung – ITV average dose and EUD (A = 1.2) (17) were also calculated for evaluation.

Iso-Energy Layer Repainting

The repainting strategy used for the 4D robust plans was layered repainting, where each energy layer is rescanned several times before switching to the next energy level. Zenklusen et al. (8) proposed two different methods to divide a plan into layers: scaled and iso-layered repainting. Scaled repainting involves simply dividing each layer in a present number of layers. In contrast, in iso-layered repainting, the MU per spot is limited by the maximum value. In this study, we segmented 4DRP by scaling the repainting 3, 4, 5, 6, or 7 times and obtained new repainting plans called 4DRP-SN3, 4DRP-SN4, 4DRP-SN5, 4DRP-SN6, and 4DRP-SN7, respectively. The layers were divided with respect to the minimum MU per spot, which means that some spots with the same energy levels consecutively contained fewer and fewer spots, but the minimum number of spots was not smaller than the limited minimum MU per spot. For the current study, the machine’s minimum MU (0.02 MU) was used during the division to ensure that all of the spots in the rescans were directly deliverable.

The optimal interplay effect mitigation in the 12 patients was assessed by comparing the value of interplay effect index and to obtain the minimum index of and over five different iso-energy repainting plans. In other words, the smaller the interplay effect index, the better the interplay effect mitigation outcome.

Statistical Analysis

Student’s t-test was used to compare the following results between 4DRP(SN3), 4DRP(SN4), 4DRP(SN5), 4DRP(SN6), and 4DRP(SN7), respectively: (1) the interplay effect in the target DVH metrics (D₂, D₅, D₉₅, D₉₈, and Target EUD) and (2) the interplay effect in lung minus ITV DVH metrics (V₅, V₂₀, V₃₀, and EUD); p < 0.05 was considered statistically significant. We compared (1) and (2) to investigate whether increasing the number repainting layers mitigated the interplay effects in the target and lungs. A linear regression model was used to evaluate the correlation between the optimal interplay mitigation vs. the tumor range of motion. This study was conducted using a linear model created with Excel software version (v.2016, Microsoft, Redmond, WA, United States)¹ to predict the correction.

Results

Interplay Effects of the 4DRP

As shown in Figure 2, the DVH target and OAR metrics changed in the 4D robust plan due to the interplay effects. Targets D₂ and D₅ increased while D₉₈ and D₉₅ obviously decreased.

FIGURE 2

The mean value of and increased to 1.28 and 1.01% according to the interplay index, and and decreased to 2.01 and 1.77%, respectively. For the OARs, the mean value of , , and increased less than 1%.

Interplay Effects After Repainting

Using patient #2’s treatment plan as an example, the 4DSD distribution, 4DDD distribution of the 4DRP, and 4DDD distributions of the 4DRP(SN3), 4DRP(SN4), 4DRP(SN5), 4DRP(SN6), and 4DRP(SN7) were investigated. In the 4DSD distribution, the isodose line of the prescription dose (PD) basically covered the tumor volume (Figure 3a). Considering the interplay effects, the distribution of the dose lines markedly deteriorated, and the PD isodose lines failed to cover the target area (Figure 3b). Further executing the repainting at different layer repainting times, the target area was again covered by the isodose lines of the PD dose lines (Figures 3c–g). No significant changes occurred in the isodose lines of the lungs.

FIGURE 3

Figure 4 shows the static interplay effect index in the target and lung DVH metrics. The interplay effect index and decreased as the number of layers increased (Figure 4A). Specifically, the mean values of were 2.01, 1.48, 1.21, 1.03, 1.01, and 0.92% and the mean values of were 1.77, 1.42, 1.13, 1.01, 0.91, and 0.89% for 4DRP(SN3), 4DRP(SN4), 4DRP(SN5), 4DRP(SN6), and 4DRP(SN7), respectively. Compared to and in 4DRP, 4DRP(SN3) was lower, with a significant difference (p < 0.05), whereas no significant differences in the other metrics were observed in the target. In the normal tissue lung, as shown in Figure 4B, the interplay effect index and increased when three layers repainting was conducted on 4DRP. The average value was 3% for 4DRP vs. 4.5% for 4DRP(SN3) in and 3.4% for 4DRP vs. 5.1% for 4DRP(SN3) in , but no significant difference was observed when the number of iso-energy layers was more than three.

FIGURE 4

Relationship Between the Optimal Interplay Effect Mitigation and Target Range of Motion

Figure 5 shows the process of exploring optimal interplay effect mitigation via multiple iso-energy layer repainting for target coverage metrics. The optimal interplay effect mitigation was attained by comparing the index among 4DRP(SN3), 4DRP(SN4), 4DRP(SN5), 4DRP(SN6), and 4DRP(SN7) with the smallest absolute value. Then, a statistical correlation study was conducted to evaluate the correlation between the optimal interplay mitigation and target range of motion that showed a linear relationship between the two (Figure 6). For , the expression was , while for , it was , where DVF is the maximum deformation vector field or target range of motion. A positive correlation was observed between the optimal interplay effect index and target range of motion.

FIGURE 5

FIGURE 6

Discussion

4D Robust optimization has been proven to mitigate the interplay effects of IMPT planning and can make the difference between the planned dose and delivery dose less than the established criterion (for example, 3%) (18). However, compared to conventional IMPT planning, proton-based SBRT has led to more consideration of the interplay effects and should have strict standards (19). Therefore, it might be ineffective to rely only on 4D Robust optimization alone to mitigate the interplay effects. Thus, we proposed 4D Robust optimization combined with an iso-energy layer repainting strategy in early-stage NSCLC patients.

In the current study, 4D robust plans were evaluated based on the quantitative interplay index for patients with early-stage NSCLC. The results showed that the target coverage decreased due to the interplay effects (Figures 2, 3). Multiple iso-energy layer repainting was subsequently used to mitigate the interplay effects. The target coverage increased as the number of iso-energy repainting layers increased (Figure 4), but this was not applicable for specific patients. Therefore, the optimal iso-energy repainting times were explored based on the quantitative interplay index. A positive linear relationship occurred between the optimal interplay effect mitigation and target range of motion. Our results confirmed that 4D Robust optimization combined with iso-energy layer repainting technology further mitigated the interplay effects, which has important clinical significance.

The motion of lung tumors was evaluated based on 4DCT images. Only 35 to 39% of the tumors moved more than 5 mm in locally advanced NSCLC, but the percentage increased to 50% for early-stage NSCLC (5). In the current study, the target range of motion was assessed using the maximum DVLs of the DIR within the target. The results showed that all of the patients had a more than 5-mm range of motion, which was also evaluated by the difference between the ITV volume and mean CTV volume over all phases (the CTV change rate in the 4DCT images). Figure 6 shows that there was a positive linear relationship between the tumor range of motion and CTV change rate. Thus, the target range of motion was obtained based on the differences in the ITV and CTV contours in the 4DCT images. DIR errors should be considered. Anaconda, implemented in RayStation, is a hybrid method utilizing a combination of image intensities and controlling structures from contoured image sets (13). Anaconda demonstrates a good performance in the thoracic region compared to other commercially available algorithms based on previous studies (14).

In practice, when using spot scanning, the number of spot rescans should be proportional to the patients’ spot weights in a range of motion up to 5 mm. However, the optimal number of repainting layers remains debatable due to patient- and machine-specific parameters, such as the patient breathing cycle, energy switch times, and other factors (20). Seco et al. (21) investigated a phase-controlled repainting and breathing-sampled strategy in which the number of rescans was decided using a motion-monitoring system. Engwall et al. (9) investigated offline breath-sampled layered repainting methods in which the number of iso-energy layers was spread uniformly throughout the breathing cycle, the optimal method of mitigating interplay effects. In the current study, multiple iso-energy layer repainting times were used to explore the optimal number of interplay mitigation rescans based on the quantitative metrics (Figure 5). A positive linear correlation occurred between the optimal interplay effect mitigation and target range of motion, demonstrating that breathing motions are dominated by the interplay effect.

One limitation of this study was associated with the method of simply dividing each layer in the present number of layers, which resulted in small weighted spots that might have been deleted after multiple iso-energy rescans. The dose distribution before and after the iso-energy rescans was compared to avoid this scenario in this study.

Conclusion

In proton-based SBRT, interplay effects degrade the target dose distribution and can be mitigated using iso-energy layer repainting techniques. However, in this study, there was no significant correlation between the number of repainting layers and improvement in the dose distributions. We recommend using the optimal layer repainting times based on the interplay effect index according to the patient characteristics.

References

• 1.

  BushDACheekGZaheerSWallenJMirshahidiHKaterelosAet alHigh-dose hypo-fractionated proton beam radiation therapy is safe and effective for central and peripheral early-stage non-small cell lung cancer: results of a 12-year experience at Loma Linda University Medical Center.Int J Radiat Oncol Biol Phys. (2013) 86:964–8. 10.1016/j.ijrobp.2013.05.002

  • CrossRef
  • Google Scholar

• 2.

  MakitaCNakamuraTTakadaATakayamaKSuzukiMAzamiYet alHigh-dose proton beam therapy for stage I non-small cell lung cancer: clinical outcomes and prognostic factors.Acta Oncologica. (2015) 54:307–14. 10.3109/0284186X.2014.948060

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 3.

  PaganettiH.Range uncertainties in proton therapy and the role of Monte Carlo simulations.Phys Med Biol. (2012) 57:R99. 10.1088/0031-9155/57/11/R99

  • CrossRef
  • Google Scholar

• 4.

  LiuHHBalterPTuttTChoiBZhangJWangCet alAssessing respiration-induced tumor motion and internal target volume using four-dimensional computed tomography for radiotherapy of lung cancer.Int J Radiat Oncol Biol Phys. (2007) 68:531–40. 10.1016/j.ijrobp.2006.12.066

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 5.

  YuZHLinSHBalterPZhangLDongL.A comparison of tumor motion characteristics between early stage and locally advanced stage lung cancers.Radiother Oncol. (2012) 104:33–8. 10.1016/j.radonc.2012.04.010

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 6.

  WeiLSchildSEChangJYLiaoZLiuZLiHet alImpact of respiratory motion on worst-case scenario optimized intensity modulated proton therapy for lung cancers.Int J Radiat Oncol Biol Phys. (2016) 95:523–33.

  • Google Scholar

• 7.

  EngwallEFredrikssonAGlimeliusL.4D robust optimization including uncertainties in time structures can reduce the interplay effect in proton pencil beam scanning radiation therapy.Med Phys. (2018a) 45:4020–29. 10.1002/mp.13094

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 8.

  ZenklusenSMPedroniEMeerD.A study on repainting strategies for treating moderately moving targets with proton pencil beam scanning at the new gantry 2 at PSI.Phys Med Biol. (2010) 55:5103–512. 10.1088/0031-9155/55/17/014

  • CrossRef
  • Google Scholar

• 9.

  EngwallEGlimeliusLHynningE.Effectiveness of different rescanning techniques for scanned proton radiotherapy in lung cancer patients.Phys Med Biol. (2018b) 63:095006. 10.1088/1361-6560/aabb7b

  • CrossRef
  • Google Scholar

• 10.

  KardarLLiYLiXLiHCaoWChangJYet alEvaluation and mitigation of the interplay effects of intensity modulated proton therapy for lung cancer in a clinical setting.Pract Radiat Oncol. (2014) 4:e259–68. 10.1016/j.prro.2014.06.010

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 11.

  YuJZhangXLiaoLLiHZhuRParkPCet alMotion-robust intensity-modulated proton therapy for distal esophageal cancer.Med Phys. (2016) 43:1111. 10.1118/1.4940789

  • CrossRef
  • Google Scholar

• 12.

  InoueTWidderJvan DijkLVTakegawaHLiHZhuRParkPCet alLimited impact of setup and range uncertainties, breathing motion, and Interplay effects in robustly optimized intensity modulated proton therapy for stage III non-small cell lung cancer.Int J Radiation Oncol Biol Phys. (2016) 96:661–9. 10.1016/j.ijrobp.2016.06.2454

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 13.

  WeistrandOSvenssonS.The anaconda algorithm for deformable image registration in radiotherapy.Med Phys. (2015) 42:40–53. 10.1118/1.4894702

  • CrossRef
  • Google Scholar

• 14.

  KadoyaNNakajimaYSaitoMMiyabeYKurookaMKitoSet alMulti-institutional validation study of commercially available deformable image registration software for thoracic images.Int J Radiat Oncol Biol Phys. (2016) 96:422–31. 10.1016/j.ijrobp.2016.05.012

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 15.

  FredrikssonAForsgrenAHardemarkB.Minimax optimization for handling range and setup uncertainties in proton therapy.Med Phys. (2011) 38:1672–84. 10.1118/1.3556559

  • CrossRef
  • Google Scholar

• 16.

  NiemierkoA.A generalized concept of equivalent uniform dose (EUD).Med Phys. (1999) 26:1100.

  • Google Scholar

• 17.

  D’AndreaMStrolinSUnganiaSCacciatoreABruzzanitiVMarconiRet alRadiobiological optimization in lung stereotactic body radiation therapy: are we ready to apply radiobiological models.Front Oncol. (2018) 7:321. 10.3389/fonc.2017.00321

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 18.

  ChangJYZhangXKnopfALiHMoriSDongLet alConsensus guidelines for implementing pencil-beam scanning proton therapy for thoracic malignancies on behalf of the PTCOG thoracic and lymphoma subcommittee.Int J Radiat Oncol Biol Phys. (2017) 99:41–50. 10.1016/j.ijrobp.2017.05.014

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 19.

  MadhavaRKNeilDShankarPG.Stereotactic body radiation therapy in lung.Transl Cancer Res. (2015) 4:435–48.

  • Google Scholar

• 20.

  SchattiAZakovaMMeerDLomaxAJ.Experimental verification of motion mitigation of discrete proton spot scanning by re-scanning.Phys Med Biol. (2013) 58:8555–72. 10.1088/0031-9155/58/23/8555

  • CrossRef
  • Google Scholar

• 21.

  SecoJRobertsonDTrofimovAPaganettiH.Breathing interplay effects during proton beam scanning: simulation and statistical analysis.Phys Med Biol. (2009) 54:N283. 10.1088/0031-9155/54/14/N01

  • CrossRef
  • Google Scholar