Lung Stereotactic Body Radiotherapy (SBRT) Using Spot-Scanning Proton Arc (SPArc) Therapy: A Feasibility Study

Purpose: We developed a 4D interplay effect model to quantitatively evaluate breathing-induced interplay effects and assess the feasibility of utilizing spot-scanning proton arc (SPArc) therapy for hypo-fractionated lung stereotactic body radiotherapy (SBRT). The model was then validated by retrospective application to clinical cases.

Materials and Methods: A digital lung 4DCT phantoms was used to mimic targets in diameter of 3cm with breathing motion amplitudes: 5, 10, 15, and 20 mm, respectively. Two planning groups based on robust optimization were generated: (1) Two-field Intensity Modulated Proton Therapy (IMPT) plans and (2) SPArc plans via a partial arc. 5,000 cGy relative biological effectiveness (RBE) was prescribed to the internal target volume (ITV) in five fractions. To quantitatively assess the breathing induced interplay effect, the 4D dynamic dose was calculated by synchronizing the breathing pattern with the simulated proton machine delivery sequence, including IMPT, Volumetric repainting (IMPT_volumetric), iso-layered repainting (IMPT_layer) and SPArc. Ten lung patients’ 4DCT previously treated with VMAT SBRT, were used to validate the digital lung tumor model. Normal tissue complicated probability (NTCP) of chestwall toxicity was calculated.

Result: Target dose were degraded as the tumor motion amplitude increased. The 4D interplay effect phantom model indicated that motion mitigation effectiveness using SPArc was about five times of IMPT_volumetric or IMPT_layer using maximum MU/spot as 0.5 MU at 20 mm motion amplitude. The retrospective study showed that SPArc has an advantage in normal tissue sparing. The probability of chestwall’s toxicity were significantly improved from 40.2 ± 29.0% (VMAT) (p = 0.01) and 16.3 ± 12.0% (IMPT) (p = 0.01) to 10.1 ± 5.4% (SPArc). SPArc could play a significant role in the interplay effect mitigation with breathing-induced motion more than 20 mm, where the target D99 of 4D dynamic dose for patient #10 was improved from 4,514 ± 138 cGy [RBE] (IMPT) vs. 4,755 ± 129 cGy [RBE] (SPArc) (p = 0.01).

Conclusion: SPArc effectively mitigated the interplay effect for proton lung SBRT compared to IMPT with repainting and was associated with normal tissue sparing. This technology may make delivery of proton SBRT more technically feasible and less complex with fewer concerns over underdosing the target compared to other proton therapy techniques.

Introduction

Lung cancer remains a leading cause of cancer mortality in the world (1). Compared to conventional radiotherapy, hypo-fractionated stereotactic body radiotherapy (SBRT) has been proved to improve local tumor control and survival rate for stage I non-small cell lung cancer (NSCLC) patients (2–6). Taking advantage of the unique beam characteristics, Bragg Peak, proton beam therapy could offer a superior dose distribution compared to photon radiotherapy technique in treating locally advanced lung cancer (7). Recently, with the development of pencil beam scanning (PBS) technology, intensity modulated proton therapy (IMPT) offers the potentials to spare the adjacent normal tissues further while maintaining similar or superior target coverage in a more efficient way without using beam specific blocks or compensators compared to passive scatter proton therapy (PSPT) (8–11). However, such scanning technique is susceptible to the interplay effect between proton spot scanning and respiratory induced motion during dose delivery. It eventually leads to an inaccurate dose delivery such as underdose targe or overdoses of the healthy tissue during lung cancer treatment (12, 13). Several motion management strategies were introduced to mitigate the interplay effect, such as repainting, gating, and tracking (14–16), in which volumetric or layer repainting technique has been widely adopted by the proton clinic. With volumetric repainting, the dose delivered during one full volume is equal to 1/N of the prescribed dose, where N was the number of rescans (14). An alternative approach is called iso-layered repainting, in which first delivered several rescans within one energy plane before switching to the next plane with the dose per spot being limited by a maximal MU value (17).

The concept of spot-scanning arc therapy (SPArc) technique was introduced in 2016 to improve the dosimetric plan quality, robustness, and delivery efficiency of proton beam therapy. The technique demonstrated potential clinical benefits in several disease sites or indications (18–23). Whether this novel technique has any potential clinical benefits in the management of stage I non-small cell lung cancer and whether it is robust enough to be implemented in the hypo-fractionated lung SBRT has yet to be explored. Therefore, we proposed a comprehensive study is to 1) to build a lung SBRT model to evaluate the effectiveness of motion interplay mitigation via SPArc quantitatively; 2) to validate the model using clinical data sets and exploit the potential benefits.

Materials and Methods

In Silico 4D Interplay Phantom Model

Due to the target deformation, motion, and imaging artifact in the four-dimensional computed tomography (4DCT), it is challenging to analyze the interplay effect quantitatively using the patient dataset directly. Previous studies have suggested using a digital lung cancer phantom as a surrogate (24, 25). By introducing a digital phantom employed in a prior study, we built an in silico 4D interplay phantom model to mimic patient’s 4DCT datasets while eliminating the artifact and target deformation uncertainties (25). Since most lung tumor motion happens in the superior-inferior (SI) direction (26), a set of digital lung tumor phantoms 4DCT with different breathing induced motion amplitudes (5, 10, 15, and 20 mm in SI direction) were created (25, 27). The target was simulated using a sphere 3 cm in diameter with 1.0 g/cc density (28), close to the average target size measured in the patient group for this study. The gross tumor volume (GTV) was contoured on the lung window through the HU (Hounsfield unit) threshold at each phase image. The internal target volume (ITV) was generated by union GTVs at each phase.

Treatment Planning on the Phantom Model

5,000 cGy relative biological effectiveness [RBE] was prescribed to ITV in five fractions SBRT with RBE = 1.1 for proton plans and RBE = 1.0 for photon plans (29). Two-field IMPT plans were generated using the single field optimization (SFO) technique via lateral and posterior beams. SPArc plans were regenerated using a partial arc from 180 to 30° clockwise with a sampling frequency of 2.5°. Both planning strategies used the same robust optimization on average CT with ±5% range and 5 mm setup uncertainties corresponding to 21 scenarios in total with a 3 mm dose grid. The minimum monitor unit (MU) threshold per spot was 0.02 MU based on the IBA proton system (19, 23, 30, 31). Similar objective constraints for organs at risk (OARs) were used in both planning groups. All plans were normalized to guarantee 99% ITV was covered by the prescription dose. The SPArc optimization algorithm starts from a multi-field IMPT with coarse sampling frequency using the worst-case scenario robust optimization and gradually resample the control point to achieve a proton arc plan (18). The algorithm integrated the iterative approaches includes (A) control point re-sampling; (B) control point energy layers re-distribution; (C) energy layers filtration; (D) energy layers re-sampling; and (E) spot number reduction by filtration. Details of the algorithm are described by Ding and Li et al. in 2016 (18).

Interplay Effect Evaluation

The 4D dynamic dose was calculated to assess the interplay effect by synchronizing the breathing pattern with the simulated proton machine delivery sequence (19, 32). To calculate a single fraction 4D dynamic dose, the dose calculated on each phase image was accumulated via the deformable image registration to the reference phase (expiration end, phase 50%) (19, 32). Ten different starting phases were simulated based on a clinical 360-degree gantry machine parameter with one revolution per minute (RPM) gantry rotation speed, 2 ms spot position switching time, energy layer switching time (ELST) of 1 s, as well as a respiratory motion period of 4 s (19). GTV D99 was assessed along with target motion amplitude variation.

A Quantitative Interplay Effect Mitigation Evaluation

In this study, IMPT treatment delivery simulation using a volumetric repainting technique was denoted as IMPT_volumetric, and IMPT using an iso-layered repainting technique denoted as IMPT_layer. To quantitatively evaluate the effectiveness of interplay effect mitigation in lung SBRT, the single fraction 4D dynamic was compared between SPArc without repainting and IMPT_volumetric, with different numbers of volumetric repainting (rescanning three, five, seven, and nine times respectively) and IMPT_layer with a series of maximum MU per spot (from 0.1 to 1.3 MU per spot).

A Retrospective Dosimetric Planning Study

Ten patients with stage I NSCLC previously treated with volumetric modulated Arc therapy (VMAT) based SBRT at our institution were selected. All patients received 4DCT simulation using a helical CT scanner (Philips Brilliance Big Bore, Philips Healthcare System, Cleveland, OH). The GTVs and ITV were generated through the strategies described above as well. The patient characteristics, including tumor location, tumor size in diameter, and tumor motion, are listed in Table 1.

TABLE 1

Table 1 Patient characteristics.

Treatment Planning in the Patient Dataset

The VMAT plans were generated using two to four partial arcs (control point frequency as 4°) based on the Elekta HD with 6 MV. The VMAT plan optimization starts from a coarse sampling of gantry position. New sample was added to achieve the desired sampling frequency, in which the Multileaf collimator (MLC) was linearly and gradually interpolated by the adjacent samples (33).

Two-field IMPT and SPArc with partial arc plans were generated, respectively. The prescription dose 5,000 cGy [RBE] was prescribed to 99% of the ITV. For a fair comparison between photon technique and proton technique, robust optimization-based ITV was used considering the same setup uncertainties as 5 mm, but ±5% range uncertainty was considered in proton planning.

Dosimetric Plan Quality Evaluation

The plan quality was evaluated based on the dose–volume histograms (DVHs) of target volume and OARs in the nominal plans. More specifically, all plans were compared for target coverage using conformity index (CI, the target volume covered by RX/the volume covered by RX). Dosimetric index for normal tissue sparing or OARs such as the Dmax or D0.1cc for the spinal cord, ribs and esophagus, Dmean for heart, ipsilateral lung (excluding ITV) as well as chest wall (CW) V30 (the volume received 3,000 cGy [RBE]) were evaluated by comparing SPArc planning group to VMAT and IMPT group.

The integral dose (ID) of radiation delivered to the whole patient body structure or external contour was analyzed (34). The ID definition was as following:

$\begin{array}{ll}\text{ID}(\text{Gy}·\text{L})=\overline{D}\left(\text{Gy}\right)·\text{V}\left(\text{L}\right)& (1)\end{array}$

where (Gy) is the mean dose delivered to volume V (L) (where L—liter).

Patient-Specific Interplay Effect Evaluation

Each case’s interplay effect was evaluated based on the 4D dynamic dose accumulation method mentioned in Interplay Effect Evaluation (19, 32).

Potential Clinical Benefit in Chestwall and Ipsilateral Lung Protection

Late chest wall toxicity after SBRT has been evaluated among three treatment modalities in this study. The probability of chest wall toxicity was calculated based on the odds ratios using the dosimetric parameter chest wall V30 (the volume of chest wall receiving 30 Gy) (35):

$\begin{array}{ll}Probability\hspace{0.5em}of\hspace{0.5em}CW\hspace{0.5em}toxicity\hspace{0.5em}b\text{y}\hspace{0.5em}V30\left(cc\right)=\frac{e[-3.151+(0.042*V30)]}{1+e[-3.151+(0.042*V30)]}& (2)\end{array}$

The incident of radiation pneumonitis for an ipsilateral lung was calculated based on the Lyman–Kutcher–Burman (LKB) model as following (36):

$\begin{array}{ll}\text{NTCP}=\frac{1}{\sqrt{2\pi }}{\int }_{-\infty }^{\text{t}}{e}^{-\frac{x^2}{2}}dx& (3)\end{array}$

$\begin{array}{ll}t=\frac{D-{TD}_{50}}{m\times {TD}_{50}}& (4)\end{array}$

where TD₅₀ is the tolerance dose for a 50% complication probability for uniform doses to the organ, and m is a dimensionless parameter for determining the slope of the complication probability according to the dose curve. And D is the equivalent dose(EUD), which is the DVH to a single dose value, representing the uniform dose that results in the survival of an equal number of clonogens in a non–homogeneously irradiated tumor. It is defined with the formula as (36):

$\begin{array}{ll}\text{EUD}={\left({\sum }_{i=1}^N{v}_i{D}_i^a\right)}^{\frac{1}{a}}& (5)\end{array}$

where Di is the dose for each bin in a differential DVH, vi is the volume in a specific dose bin i, and N is the unequal fractional sub-volume. The ‘a’ value is a parameter equal to 1/n, in which n represents the volume dependence of the complication probability. The parameter set for the lung tissue were taken from Burman et al. (TD50 = 24.5 Gy, m = 0.18, and a = 0.87) (37, 38). According to the LQ model, the dose axis of the DVH was re-scaled to the equivalent dose in 2 Gy per fractions and the method described by Van den Heuvel with assuming α/β ratio of 3 Gy for the ipsilateral normal lung tissue (39, 40).

The dosimetric index from SPArc was utilized as a reference. By comparing with two other treatment technologies (IMPT and VMAT), the differences were assessed with a paired, 2-tailed non-parametric Wilcoxon signed-rank test via SPSS 21.0 software (International Business Machines, Armonk, New York), respectively, and p values less than 0.05 were considered statistically significant.

Result

In Silico 4D Interplay Phantom Model

The Interplay Effect Evaluation

The study showed that the dose to target degraded as the tumor motion amplitude increased in the IMPT planning, which was agreed with the previous reports (16, 41, 42). Figure 1 displayed the target D99 in relationship with various motion amplitudes. SPArc could significantly improve the target coverage compared to IMPT through all the different motion amplitudes, even without any repainting. More specially, the average relative target D99 degradation via single fraction 4D dynamic dose accumulation were 2.51 vs 0.00% (p <0.01), 4.01 vs 0.10% (p <0.01), 6.61 vs 1.29% (p <0.01), 8.40 vs 1.70% (p <0.01) for IMPT vs SPArc at different breathing amplitude (5, 10, 15 and 20 mm), respectively.

FIGURE 1

Figure 1 Single fraction dynamic dose for target D99 along with different motion amplitude from 5 to 20 mm.

Comparison of Mitigation Effectiveness in the Interplay Effect With Repainting IMPT

The motion interplay effect can be compensated by increasing the total number of volumetric repainting or constraining the maximum MU per spot using iso-layered repainting. The single fraction 4D dynamic dose accumulation for target D99 with the number of volumetric repainting times and maximum MU per spot for the target motion 10 and 20 mm was displayed in Figures 2A, B, respectively. At 10 mm target motion amplitude, GTV D99 was 4,767 ± 63 cGy [RBE] (p <0.01) in IMPT without repainting, which is less than 4,950 ± 41 cGy [RBE] in SPArc. IMPT_volumetric increased GTV D99 to 4,959 ± 76 cGy [RBE] (p = 0.51) and 4,985 ± 66 cGy [RBE] (p = 0.33) with three and five times of volumetric repainting (Figure 2A); IMPT_layer increased GTV D99 to 4,931 ± 78 cGy [RBE] (p = 0.39) and 4,981 ± 66cGy [RBE] (p = 0.11) with maximum MU per spot as 0.75 and 0.50 MU respectively (Figure 2B), compared to SPArc. It is interesting to find that SPArc is as effective as three to five times of volumetric repainting IMPT or iso-layered repainting with maximum MU per spot as 0.75 to 0.5 MU at 10 mm target motion amplitude.

FIGURE 2

Figure 2 Single fraction 4D dynamic dose comparison between SPArc (green) and (A) different number of volumetric repainting times IMPT_volumetric, and (B) IMPTlayer with different maximum MU per spot.

Moreover, in the target motion with 20 mm amplitude, GTV D99 was 4,532 ± 180 cGy [RBE] (IMPT without repainting) vs 4,902 ± 94 cGy [RBE] (SPArc) (p = 0.01). SPArc was as effective as five to seven times of volumetric repainting or iso-layered repainting with maximum MU per spot as 0.7 to 0.4 MU in IMPT for where GTV D99 received 4,896 ± 75 cGy [RBE] (p = 0.96) and 4,912 ± 26 cGy [RBE] (p = 0.65) for IMPT with volumetric repainting five and seven times, respectively (Figure 2A). Meanwhile, GTV D99 reached as 4,841 ± 102 cGy [RBE] (p = 0.09) and 4,929 ± 71 cGy [RBE] (p = 0.65) during IMPT_layer with maximum MU per spot of 0.7 and 0.4 MU (Figure 2B), respectively.

Retrospective Study Using Patient Dataset

Plan Quality Evaluation

Taking advantage of more degrees of freedom in the optimization through arc(s) trajectory, VMAT and SPArc planning groups demonstrated superior dose conformity to the target. Isodose distributions of patient #10 for treatment plan using VMAT (first column), IMPT (the second column), and SPArc (the third column) was displayed in Figure 3.

FIGURE 3

Figure 3 Isodose distributions of patient #10 for treatment plan using VMAT (first column), IMPT (the second column), and SPArc (the third column). The green contour represents ITV. 100% dose is equal to prescription dose.

SPArc improved CI from 0.31 ± 0.08 in IMPT to 0.38 ± 0.10 (p = 0.01). This feature allows SPArc to spare more OARs such as spinal cord and ribs than IMPT. In addition, SPArc plans significantly reduced the Dmean of the ipsilateral lung from 503 ± 1 76 cGy [RBE] to 418 ± 140 cGy [RBE] (p = 0.01) and Dmax of ribs from 4,369 ± 978 cGy [RBE] to 4,151 ± 1,015 cGy [RBE] (p = 0.02) compared to IMPT respectively. V30 of the chest wall was significantly reduced from 30 ± 22 cc to 20 ± 14 cc (p = 0.02) compared to IMPT (Table 2).

TABLE 2

Table 2 Dosimetry results for the three planning modalities.

Compared to VMAT, SPArc significantly reduced the dose to OARs listed in Table 2. More specially, SPArc significantly reduced maximum dose to spinal cord: 1,026 ± 494 cGy [RBE] (VMAT) vs 300 ± 530 cGy [RBE] (SPArc) (p = 0.01), esophagus: 1,611 ± 1,361 cGy [RBE] (VMAT) vs 501 ± 1,565 cGy [RBE] (SPArc) (p <0.01), ribs: 4,770 ± 1,059 cGy [RBE] (VMAT) vs 4,151 ± 1,015 cGy [RBE] (SPArc) (p = 0.01). In addition, SPArc significantly reduced the mean dose of the ipsilateral lung from 659 ± 200 cGy [RBE] (VMAT) to 418 ± 140 cGy [RBE] (SPArc) (p = 0.01), and the mean dose of heart from 288 ± 253 cGy [RBE] (VMAT) to 8 ± 11 cGy [RBE] (SPArc) (p = 0.01). The chest wall V30 was also significantly reduced via SPArc plan from 60 ± 37 cc to 20 ± 14 cc (p <0.01). The study also found that SPArc (16.13 ± 5.36 Gy·L) significantly reduced ID comparing to both IMPT (18.40 ± 5.79 Gy·L, p = 0.01) and VMAT (38.38 ± 17.10 Gy·L, p = 0.01) planning group.

Interplay Effect Evaluation in the Patient Population

The single fraction 4D dynamic dose accumulation showed that target D99 was degraded due to the interplay effect in IMPT without repainting (Figure 4) among the ten patients. Similar to the 4D interplay phantom model, the retrospective study showed a similar trend of target dose coverage degradation with increased breathing-induced motion amplitudes. SPArc significantly mitigated interplay effect compared to IMPT among cases where breathing-induced motion cannot be ignored. Even though the target’s motion and shape were complicated as these parameters are patient-specific, the patient cases’ trend was consistent with the 4D interplay phantom study (Comparison of Mitigation Effectiveness in the Interplay Effect With Repainting IMPT) when the amplitudes of breathing induced motion is large (Figure 4). More specifically, the target D99 of 4D dynamic dose was 4,514 ± 138 cGy [RBE] (IMPT without repainting) vs 4,755 ± 129 cGy [RBE] (SPArc) (p = 0.01) fort the patient #10.

FIGURE 4

Figure 4 Single fraction dynamic dose for target coverage D99 for ten patients.

The Probability of Chest Wall Toxicity and Radiation-Induced Pneumonitis

Due to SPArc significantly spared chest wall better than both IMPT and VMAT. The consequence clinic benefit was obvious, where the probability of CW toxicity were improved from 40.2 ± 29.0% (VMAT) (p = 0.01) and 16.3 ± 12.0 (IMPT) (p = 0.01) to 10.1 ± 5.4% (SPArc).

All three treatment modalities were able to spare ipsilateral lung tissue well. The corresponding incidence of radiation pneumonitis was fairly low, all of them were approximate to 0% on average.

Discussion

This is the first study to explore the feasibility of using SPArc in hypo-fractionated treatment in mobile targets and its potential clinical benefits in lung SBRT. We quantify the effectiveness of using SPArc to mitigate the motion interplay effect using digital 4D lung cancer phantom and validate the model through a retrospective dosimetric study. The result confirmed the previous report that the interplay effect led to a deterioration of the dose distribution (16, 41, 42). The larger amplitude in motion, the more deterioration in the target coverage. Such a trend was consistent within the retrospective dosimetric study findings in the ten patient cases (Figure 4), even though the tumor shape, size, 3-dimentional tumor motion, and density variation during the breathing cycle were complicated for the patient population. For the patient group with 6 mm <target offset ≤10 mm, the target D99 degradation was 3.36 ± 0.55% on average for IMPT, and 2.05 ± 0.53% SPArc. Thus, it indicated that SPArc or IMPT with repainting was preferred to lung SBRT rather than IMPT alone without repainting. When the target offset >10 mm, the relative target dose D99 degradation was greater than 3% in both phantom and patient cases, indicating that IMPT poses a potential risk in missing part of the target in lung SBRT. In contrast, SPArc could mitigate the target dose degradation caused by interplay effect well in the 4D phantom model. The phantom study indicated that SPArc is as effective as five times of volumetric repainting IMPT in terms of interplay effect mitigation at 10 mm target motion amplitude. Only SI directional and rigid motion being considered in 4D phantom model, it required further investigation since 3D motion and complicated shape changes occurred for the clinical patients. Despite the effectiveness of motion interplay mitigation SPArc offered was compromised in the patient group, it was significantly superior to IMPT without repainting.

The previous study indicated that repainting might not be needed when multiple filed are applied with target amplitude up to 6 mm, because dose blurring effects appear negligible between standard delivery and repainting technique (43). A similar phenomenon was observed in our study as well, where the interplay effect between SPArc and IMPT is very close with target motion less than 6 mm. Additionally, Knopf et al. also demonstrated that IMPT with multiple beams was able to mitigate the interplay effect for targets with large motion amplitude (43), which was confirmed in our study. Our result indicated that SPArc plan could offer superior interplay effect mitigation through applying many beam angles via arc trajectory for the target motion more than10mm, compared to IMPT with two beams.

Lung SBRT for Stage I NSCLC is a highly effective treatment that is being increasingly utilized (4, 44, 45). Lung SBRT is characterized by using a hypo-fractionated treatment course with a biological equivalent dose of at least 100 Gy. Proton lung SBRT offers increased conformality compared to photon lung SBRT; however, there is uncertainty in tumor coverage mainly due to the interplay effect. Thus, proton lung SBRT commonly has been described with passive scattering techniques typically using at least ten fractions (46, 47). Chen et al. reported on using lung SBRT with IMPT with a few patients receiving eight fractions, although the majority had at least ten fractions (48). This study indicated that SPArc’s ability to mitigate the interplay effect could improve the normal tissue toxicity while also providing the means to use three to five fraction regimens commonly used with photon SBRT. Even single fractions of photon lung SBRT were shown in RTOG 0915 to have comparable efficacy and toxicity to four fractions (48, 49). The logistical and financial benefits of hypofractionation are attractive in this population of patients with multiple medical comorbidities. SPArc’s increased robustness enables the use.

Other motion management strategies such as passive pressure technique, gating, the breath-hold approach could be implemented in the clinical practice, but these procedures prolong treatment time, which causes additional intra-fractionation motion or setup uncertainties (50). This study demonstrated that SPArc could effectively mitigate motion interplay. This new finding of utilizing SPArc to mitigate the interplay effect opens a new direction of motion management strategy by increase the degree of freedom such as arc(s) trajectory to effectively reduce the dosimetric impact from each beam’s direction. More specifically, for patient #10, where the breathing-induced motion exceeds 2 cm, the effectiveness of interplay effect mitigation of using SPArc technique reached five times of IMPT with volumetric repainting or IMPT using iso-layered repainting with maximum MU per spot as 0.5 MU at most, in which the corresponding GTV D99 were 4,822 ± 98 cGy [RBE] (p = 0.39) and 4,854 ± 86 cGy [RBE] (p = 0.14) (Figures 5A, B).

FIGURE 5

Figure 5 Single fraction 4D dynamic dose of patient #10 for target D99 along with (A) volumetric repainting times and (B) maximum MU per spot. The top value is the p value for the comparison of target D99 between IMPT with repainting and SPArc without repainting.

In general, SBRT is associated with a low incidence of acute and late toxicity. However, late chest well toxicity such as chest pain has been reported, typically mild to moderate. Moreover, chest wall pain commonly occurs with a median time of onset of greater than six months after the treatment. This study shows that SPArc has significantly spared the chest wall V30. Consequently, the probability of CW toxicity was improved considerably from 40.2 ± 29.0% (VMAT) (p = 0.01) and 16.3 ± 12.0% (IMPT) (P = 0.01) to 10.1 ± 5.4% (SPArc), which would improve the probability of chest wall toxicity and patient’s life quality.

Conclusion

A 4D interplay digital phantom model for mobile lung target was established to evaluate the effectiveness of interplay effect mitigation quantitatively. SPArc, as a novel proton treatment technique, could significantly reduce the dosimetric impact from the interplay effect and potentially reduce the Chestwall pain in lung SBRT.

Data Availability Statement

The raw data supporting the conclusions of this article will be made available by the authors, without undue reservation.

Ethics Statement

The studies involving human participants were reviewed and approved by Beaumont Health Institutional Review Board. Written informed consent for participation was not required for this study in accordance with the national legislation and the institutional requirements.

Author Contributions

GL contributed to the acquisition, analysis, and interpretation of data and drafted and designed the paper. XL provided technical support. IG and RD contributed to revising the paper and providing clinical inputs. AQ, DY, CS, and SZ provided support in imaging acquisition. LZ contributed to the statistical analysis. XD contributed to the design of the study, revised the draft, and led the research direction. All authors contributed to the article and approved the submitted version.

Funding

The study is supported by Ion Beam Application S.A. (IBA, Belgium). The funder was not involved in the study design, collection, analysis, interpretation of data, the writing of this article or the decision to submit it for publication.

Conflict of Interest

XD, XL, and DY have a patent related to spot-scanning arc therapy.

The remaining authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

References

1. Grills IS, Hugo G, Kestin LL, Galerani AP, Chao KK, Wloch J, et al. Image-Guided Radiotherapy via Daily Online Cone-Beam CT Substantially Reduces Margin Requirements for Stereotactic Lung Radiotherapy. Int J Radiat Oncol Biol Phys (2008) 70:1045–56. doi: 10.1016/j.ijrobp.2007.07.2352

PubMed Abstract | CrossRef Full Text | Google Scholar

2. RTOG 0618. Stereotactic body radiation therapy (SBRT) to treat operable early-stage lung cancer patients. Available at: https://ascopubs.org/doi/abs/10.1200/jco.2013.31.15_suppl.7523 (Accessed March 12, 2021).

Google Scholar

3. Timmerman RD, Hu C, Michalski JM, Bradley JC, Galvin J, Johnstone DW, et al. Long-term Results of Stereotactic Body Radiation Therapy in Medically Inoperable Stage I Non–Small Cell Lung Cancer. JAMA Oncol (2018) 4:1287. doi: 10.1001/jamaoncol.2018.1258

PubMed Abstract | CrossRef Full Text | Google Scholar

4. Fakiris AJ, McGarry RC, Yiannoutsos CT, Papiez L, Williams M, Henderson MA, et al. Stereotactic body radiation therapy for early-stage non-small-cell lung carcinoma: four-year results of a prospective phase II study. Int J Radiat Oncol Biol Phys (2009) 75:677–82. doi: 10.1016/j.ijrobp.2008.11.042

PubMed Abstract | CrossRef Full Text | Google Scholar

5. Onishi H, Shirato H, Nagata Y, Hiraoka M, Fujino M, Gomi K, et al. Hypofractionated stereotactic radiotherapy (HypoFXSRT) for stage I non-small cell lung cancer: updated results of 257 patients in a Japanese multi-institutional study. J Thorac Oncol (2007) 2:S94–100. doi: 10.1097/JTO.0b013e318074de34

PubMed Abstract | CrossRef Full Text | Google Scholar

6. Uematsu M, Shioda A, Suda A, Fukui T, Ozeki Y, Hama Y, et al. Computed tomography-guided frameless stereotactic radiotherapy for stage I non-small cell lung cancer: a 5-year experience. Int J Radiat Oncol Biol Phys (2001) 51:666–70.

PubMed Abstract | Google Scholar

7. Newhauser WD, Zhang R. The physics of proton therapy. Phys Med Biol (2015) 60:R155–209. doi: 10.1088/0031-9155/60/8/R155

PubMed Abstract | CrossRef Full Text | Google Scholar

8. Georg D, Hillbrand M, Stock M, Dieckmann K, Pötter R. Can protons improve SBRT for lung lesions? Dosimetric considerations. Radiother Oncol (2008) 88:368–75. doi: 10.1016/j.radonc.2008.03.007

PubMed Abstract | CrossRef Full Text | Google Scholar

9. Zhang X, Li Y, Pan X, Xiaoqiang L, Mohan R, Komaki R, et al. Intensity-Modulated Proton Therapy Reduces the Dose to Normal Tissue Compared With Intensity-Modulated Radiation Therapy or Passive Scattering Proton Therapy and Enables Individualized Radical Radiotherapy for Extensive Stage IIIB Non-Small-Cell Lung Cancer: A Virtual Clinical Study. Int J Radiat Oncol Biol Phys (2010) 77:357–66. doi: 10.1016/j.ijrobp.2009.04.028

PubMed Abstract | CrossRef Full Text | Google Scholar

10. Ding X, Dionisi F, Tang S, Ingram M, Hung C-Y, Prionas E, et al. A comprehensive dosimetric study of pancreatic cancer treatment using three-dimensional conformal radiation therapy (3DCRT), intensity-modulated radiation therapy (IMRT), volumetric-modulated radiation therapy (VMAT), and passive-scattering and modulated-scanning proton therapy (PT). Med Dosimetry (2014) 39:139–45. doi: 10.1016/j.meddos.2013.11.005

CrossRef Full Text | Google Scholar

11. Chang JY, Zhang X, Wang X, Kang Y, Riley B, Bilton S, et al. Significant reduction of normal tissue dose by proton radiotherapy compared with three-dimensional conformal or intensity-modulated radiation therapy in Stage I or Stage III non–small-cell lung cancer. Int J Radiat Oncol Biol Phys (2006) 65:1087–96. doi: 10.1016/j.ijrobp.2006.01.052

PubMed Abstract | CrossRef Full Text | Google Scholar

12. Engwall E, Glimelius L, Hynning E. Effectiveness of different rescanning techniques for scanned proton radiotherapy in lung cancer patients. Phys Med Biol (2018) 63:095006. doi: 10.1088/1361-6560/aabb7b

PubMed Abstract | CrossRef Full Text | Google Scholar

13. Boria AJ, Uh J, Pirlepesov F, Stuckey JC, Axente M, Gargone MA, et al. Interplay Effect of Target Motion and Pencil-Beam Scanning in Proton Therapy for Pediatric Patients. Int J Particle Ther (2018) 5:1–10. doi: 10.14338/IJPT-17-00030.1

CrossRef Full Text | Google Scholar

14. Bernatowicz K, Lomax AJ, Knopf A. Comparative study of layered and volumetric rescanning for different scanning speeds of proton beam in liver patients. Phys Med Biol (2013) 58:7905–20. doi: 10.1088/0031-9155/58/22/7905

PubMed Abstract | CrossRef Full Text | Google Scholar

15. Court LE, Wagar M, Ionascu D, Berbeco R, Chin L. Management of the interplay effect when using dynamic MLC sequences to treat moving targets: Managing the interplay effect. Med Phys (2008) 35:1926–31. doi: 10.1118/1.2896083

PubMed Abstract | CrossRef Full Text | Google Scholar

16. Schätti A, Zakova M, Meer D, Lomax AJ. Experimental verification of motion mitigation of discrete proton spot scanning by re-scanning. Phys Med Biol (2013) 58:8555–72. doi: 10.1088/0031-9155/58/23/8555

PubMed Abstract | CrossRef Full Text | Google Scholar

17. Zenklusen SM, Pedroni E, Meer D. 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–21. doi: 10.1088/0031-9155/55/17/014

PubMed Abstract | CrossRef Full Text | Google Scholar

18. Ding X, Li X, Zhang JM, Kabolizadeh P, Stevens C, Yan D. Spot-Scanning Proton Arc (SPArc) Therapy: The First Robust and Delivery-Efficient Spot-Scanning Proton Arc Therapy. Int J Radiat Oncol Biol Phys (2016) 96:1107–16. doi: 10.1016/j.ijrobp.2016.08.049

PubMed Abstract | CrossRef Full Text | Google Scholar

19. Li X, Kabolizadeh P, Yan D, Qin A, Zhou J, Hong Y, et al. Improve dosimetric outcome in stage III non-small-cell lung cancer treatment using spot-scanning proton arc (SPArc) therapy. Radiat Oncol (2018) 13:35. doi: 10.1186/s13014-018-0981-6

PubMed Abstract | CrossRef Full Text | Google Scholar

20. Chang S, Liu G, Zhao L, Dilworth JT, Zheng W, Jawad S, et al. Feasibility study: spot-scanning proton arc therapy (SPArc) for left-sided whole breast radiotherapy. Radiat Oncol (2020) 15:232. doi: 10.1186/s13014-020-01676-3

PubMed Abstract | CrossRef Full Text | Google Scholar

21. Liu G, Li X, Qin A, Zheng W, Yan D, Zhang S, et al. Improve the dosimetric outcome in bilateral head and neck cancer (HNC) treatment using spot-scanning proton arc (SPArc) therapy: a feasibility study. Radiat Oncol (2020) 15:21. doi: 10.1186/s13014-020-1476-9

PubMed Abstract | CrossRef Full Text | Google Scholar

22. Ding X, Li X, Qin A, Zhou J, Yan D, Stevens C, et al. Have we reached proton beam therapy dosimetric limitations? – A novel robust, delivery-efficient and continuous spot-scanning proton arc (SPArc) therapy is to improve the dosimetric outcome in treating prostate cancer. Acta Oncol (2018) 57:435–7. doi: 10.1080/0284186X.2017.1358463

PubMed Abstract | CrossRef Full Text | Google Scholar

23. Ding X, Zhou J, Li X, Blas K, Liu G, Wang Y, et al. Improving dosimetric outcome for hippocampus and cochlea sparing whole brain radiotherapy using spot-scanning proton arc therapy. Acta Oncol (2019) 0:1–8. doi: 10.1080/0284186X.2018.1555374

CrossRef Full Text | Google Scholar

24. Zhang F, Kelsey CR, Yoo D, Yin F-F, Cai J. Uncertainties of 4-dimensional computed tomography-based tumor motion measurement for lung stereotactic body radiation therapy. Pract Radiat Oncol (2014) 4:e59–65. doi: 10.1016/j.prro.2013.02.009

PubMed Abstract | CrossRef Full Text | Google Scholar

25. Liu G, Hu F, Ding X, Li X, Shao Q, Wang Y, et al. Simulation of dosimetry impact of 4DCT uncertainty in 4D dose calculation for lung SBRT. Radiat Oncol (2019) 14:1–25. doi: 10.1186/s13014-018-1191-y

PubMed Abstract | CrossRef Full Text | Google Scholar

26. Glide-Hurst CK, Hugo GD, Liang J, Yan D. A simplified method of four-dimensional dose accumulation using the mean patient density representation: 4D dose accumulation using mean patient density representation. Med Phys (2008) 35:5269–77. doi: 10.1118/1.3002304

PubMed Abstract | CrossRef Full Text | Google Scholar

27. Bissonnette J-P, Franks KN, Purdie TG, Moseley DJ, Sonke J-J, Jaffray DA, et al. Quantifying Interfraction and Intrafraction Tumor Motion in Lung Stereotactic Body Radiotherapy Using Respiration-Correlated Cone Beam Computed Tomography. Int J Radiat Oncol Biol Phys (2009) 75:688–95. doi: 10.1016/j.ijrobp.2008.11.066

PubMed Abstract | CrossRef Full Text | Google Scholar

28. Suryanto A, Herlambang K, Rachmatullah P. Comparison of tumor density by CT scan based on histologic type in lung cancer patients. Acta Med Indones (2005) 37:195–8.

PubMed Abstract | Google Scholar

29. Paganetti H, Niemierko A, Ancukiewicz M, Gerweck LE, Goitein M, Loeffler JS, et al. Relative biological effectiveness (RBE) values for proton beam therapy. Int J Radiat Oncol Biol Phys (2002) 53:407–21. doi: 10.1016/S0360-3016(02)02754-2

PubMed Abstract | CrossRef Full Text | Google Scholar

30. Liu W, Frank SJ, Li X, Li Y, Park PC, Dong L, et al. Effectiveness of robust optimization in intensity-modulated proton therapy planning for head and neck cancers: Robust optimization for IMPT for H&N cancer. Med Phys (2013) 40:051711. doi: 10.1118/1.4801899

PubMed Abstract | CrossRef Full Text | Google Scholar

31. Liu W, Zhang X, Li Y, Mohan R. Robust optimization of intensity modulated proton therapy: Robust optimization of IMPT. Med Phys (2012) 39:1079–91. doi: 10.1118/1.3679340

PubMed Abstract | CrossRef Full Text | Google Scholar

32. Li H, Li Y, Zhang X, Li X, Liu W, Gillin MT, et al. Dynamically accumulated dose and 4D accumulated dose for moving tumors: Dynamic dose and 4D dose. Med Phys (2012) 39:7359–67. doi: 10.1118/1.4766434

PubMed Abstract | CrossRef Full Text | Google Scholar

33. Otto K. Volumetric modulated arc therapy: IMRT in a single gantry arc: Single arc radiation therapy. Med Phys (2007) 35:310–7. doi: 10.1118/1.2818738

CrossRef Full Text | Google Scholar

34. Aoyama H, Westerly DC, Mackie TR, Olivera GH, Bentzen SM, Patel RR, et al. Integral radiation dose to normal structures with conformal external beam radiation. Int J Radiat Oncol Biol Phys (2006) 64:962–7. doi: 10.1016/j.ijrobp.2005.11.005

PubMed Abstract | CrossRef Full Text | Google Scholar

35. Stephans KL, Djemil T, Tendulkar RD, Robinson CG, Reddy CA, Videtic GMM. Prediction of Chest Wall Toxicity From Lung Stereotactic Body Radiotherapy (SBRT). Int J Radiat Oncol Biol Phys (2012) 82:974–80. doi: 10.1016/j.ijrobp.2010.12.002

PubMed Abstract | CrossRef Full Text | Google Scholar

36. Daly ME, Luxton G, Choi CYH, Gibbs IC, Chang SD, Adler JR, et al. Normal Tissue Complication Probability Estimation by the Lyman-Kutcher-Burman Method Does Not Accurately Predict Spinal Cord Tolerance to Stereotactic Radiosurgery. Int J Radiat Oncol Biol Phys (2012) 82:2025–32. doi: 10.1016/j.ijrobp.2011.03.004

PubMed Abstract | CrossRef Full Text | Google Scholar

37. Willner J, Jost A, Baier K, Flentje M. A Little to a Lot or a Lot to a Little?: An Analysis of Pneumonitis Risk from Dose-Volume Histogram Parameters of theLung in Patients with Lung Cancer Treated with 3-D Conformal Radiotherapy. Strahlenther Onkol (2003) 179:548–56. doi: 10.1007/s00066-003-1078-0

PubMed Abstract | CrossRef Full Text | Google Scholar

38. Burman C, Kutcher GJ, Emami B, Goitein M. Fitting of normal tissue tolerance data to an analytic function. Int J Radiat Oncol Biol Phys (1991) 21:123–35. doi: 10.1016/0360-3016(91)90172-Z

PubMed Abstract | CrossRef Full Text | Google Scholar

39. Ricardi U, Filippi AR, Guarneri A, Giglioli FR, Mantovani C, Fiandra C, et al. Dosimetric predictors of radiation-induced lung injury in stereotactic body radiation therapy. Acta Oncol (2009) 48:571–7. doi: 10.1080/02841860802520821

PubMed Abstract | CrossRef Full Text | Google Scholar

40. Zimmermann FB, Geinitz H, Schill S, Thamm R, Nieder C, Schratzenstaller U, et al. Stereotactic hypofractionated radiotherapy in stage I (T1-2 N0 M0) non-small-cell lung cancer (NSCLC). Acta Oncol (2006) 45:796–801. doi: 10.1080/02841860600913210

PubMed Abstract | CrossRef Full Text | Google Scholar

41. Grassberger C, Dowdell S, Lomax A, Sharp G, Shackleford J, Choi N, et al. Motion Interplay as a Function of Patient Parameters and Spot Size in Spot Scanning Proton Therapy for Lung Cancer. Int J Radiat Oncol Biol Phys (2013) 86:380–6. doi: 10.1016/j.ijrobp.2013.01.024

PubMed Abstract | CrossRef Full Text | Google Scholar

42. Seco J, Robertson D, Trofimov A, Paganetti H. Breathing interplay effects during proton beam scanning: simulation and statistical analysis. Phys Med Biol (2009) 54:N283–94. doi: 10.1088/0031-9155/54/14/N01

PubMed Abstract | CrossRef Full Text | Google Scholar

43. Knopf A-C, Hong TS, Lomax A. Scanned proton radiotherapy for mobile targets—the effectiveness of re-scanning in the context of different treatment planning approaches and for different motion characteristics. Phys Med Biol (2011) 56:7257–71. doi: 10.1088/0031-9155/56/22/016

PubMed Abstract | CrossRef Full Text | Google Scholar

44. Zhang Y, Chen Y, Qiu J, Yang J. Dosimetric Comparisons of Lung SBRT with Multiple Metastases by Two Advanced Planning Systems. Int J Med Phys Clin Eng Radiat Oncol (2014) 03:252–61. doi: 10.4236/ijmpcero.2014.34032

CrossRef Full Text | Google Scholar

45. Nagata Y, Takayama K, Matsuo Y, Norihisa Y, Mizowaki T, Sakamoto T, et al. Clinical outcomes of a phase I/II study of 48 Gy of stereotactic body radiotherapy in 4 fractions for primary lung cancer using a stereotactic body frame. Int J Radiat Oncol Biol Phys (2005) 63:1427–31. doi: 10.1016/j.ijrobp.2005.05.034

PubMed Abstract | CrossRef Full Text | Google Scholar

46. Liu Q, Zhu Z, Chen Y, Deng J, Ai D, Liu Q, et al. Phase 2 Study of Stereotactic Body Radiation Therapy for Patients with Oligometastatic Esophageal Squamous Cell Carcinoma. Int J Radiat Oncol Biol Phys (2020) 108:707–15. doi: 10.1016/j.ijrobp.2020.05.003

PubMed Abstract | CrossRef Full Text | Google Scholar

47. Westover KD, Seco J, Adams JA, Lanuti M, Choi NC, Engelsman M, et al. Proton SBRT for medically inoperable stage I NSCLC. J Thorac Oncol (2012) 7:1021–5. doi: 10.1097/JTO.0b013e31824de0bf

PubMed Abstract | CrossRef Full Text | Google Scholar

48. Videtic GM, Paulus R, Singh AK, Chang JY, Parker W, Olivier KR, et al. Long-term Follow-up on NRG Oncology RTOG 0915 (NCCTG N0927): A Randomized Phase 2 Study Comparing 2 Stereotactic Body Radiation Therapy Schedules for Medically Inoperable Patients With Stage I Peripheral Non-Small Cell Lung Cancer. Int J Radiat Oncol Biol Phys (2019) 103:1077–84. doi: 10.1016/j.ijrobp.2018.11.051

PubMed Abstract | CrossRef Full Text | Google Scholar

49. Videtic GM, Paulus R, Singh AK, Chang JY, Parker W, Olivier K, et al. Long-Term Follow-Up on NRG Oncology RTOG 0915 (NCCTG N0927): A Randomized Phase 2 Study Comparing 2 Stereotactic Body Radiation Therapy Schedules for Medically Inoperable Patients with Stage I Peripheral Non–small Cell Lung Cancer. Int J Radiat Oncol Biol Phys (2017) 99:S15–6. doi: 10.1016/j.ijrobp.2017.06.052

CrossRef Full Text | Google Scholar

50. Engwall E, Glimelius L, Hynning E. Effectiveness of different rescanning techniques for scanned proton radiotherapy in lung cancer patients. Phys Med Biol (2018) 63:095006. doi: 10.1088/1361-6560/aabb7b

PubMed Abstract | CrossRef Full Text | Google Scholar