Skip to main content

Radiobiological model-based approach to determine the potential of dose-escalated robust intensity-modulated proton radiotherapy in reducing gastrointestinal toxicity in the treatment of locally advanced unresectable pancreatic cancer of the head



The purpose of this study was to determine the potential of escalated dose radiation (EDR) robust intensity-modulated proton radiotherapy (ro-IMPT) in reducing GI toxicity risk in locally advanced unresectable pancreatic cancer (LAUPC) of the head in term of normal tissue complication probability (NTCP) predictive model.


For 9 patients, intensity-modulated radiotherapy (IMRT) was compared with ro-IMPT. For all plans, the prescription dose was 59.4GyE (Gray equivalent) in 33 fractions with an equivalent organ at risk (OAR) constraints. Physical dose distribution was evaluated. GI toxicity risk for different endpoints was estimated using published NTCP Lyman Kutcher Burman (LKB) models for stomach, duodenum, small bowel, and combine stomach and duodenum (Stoduo). A Wilcoxon signed-rank test was used for dosimetry parameters and NTCP values comparison.


The dosimetric results have shown that, with similar target coverage, ro-IMPT achieves a significant dose-volume reduction in the stomach, small bowel, and stoduo in low to high dose range in comparison to IMRT. NTCP evaluation for the endpoint gastric bleeding of stomach (10.55% vs. 13.97%, P = 0.007), duodenum (1.87% vs. 5.02%, P = 0.004), and stoduo (5.67% vs. 7.81%, P = 0.008) suggest reduced toxicity by ro-IMPT compared to IMRT. ∆NTCP IMRT – ro-IMPT (using parameter from Pan et al. for gastric bleed) of ≥5 to < 10% was seen in 3 patients (33%) for stomach and 2 patients (22%) for stoduo. An overall GI toxicity relative risk (NTCPro-IMPT/NTCPIMRT) reduction was noted (0.16–0.81) for all GI-OARs except for duodenum (> 1) with endpoint grade ≥ 3 GI toxicity (using parameters from Holyoake et al.).


With similar target coverage and better conformity, ro-IMPT has the potential to substantially reduce the risk of GI toxicity compared to IMRT in EDR of LAUPC of the head. This result needs to be further evaluated in future clinical studies.


Pancreatic cancer is the fourth primary cause of cancer-related death in Japan [1]. Locally advanced unresectable pancreatic cancer (LAUPC) has a 5-year survival of < 5% [2]. The main treatment option for LAUPC is chemotherapy, while definitive chemoradiotherapy (CRT) plays a pivotal role in therapeutic management to enhance survival and quality of life [3].

Despite the use of modern radiotherapy (RT) delivery techniques, potentially gastrointestinal (GI), severe toxicity had been noted [4, 5]. Studies using escalated dose radiation (EDR) intensity-modulated radiotherapy (IMRT), PTV-based double scatter (DS), and pencil beam scanning (PBS) proton beam therapy (PBT) have reported improved local control and survival with the main limiting GI toxicity [6,7,8]. A study by Kelly et al. and Ben-Josef et al., in LAUPC using EDR-IMRT have reported ≥ grade 2 GI toxicity in approximately 15 to 20% patients [6, 9]. The study by Takatori et al., using hypo-fractionated concurrent gemcitabine PBT (GPT) for LAUPC 67.5 Gray equivalent (GyE) in 25 fractions) have reported 49.4% rate of gastric/duodenal ulcer [7]. The study by Terashima et al. treated 45 patients with hypo-fractionated GPT have reported 10% of grade ≥ 3 late gastric ulcer and hemorrhage [8]. However, studies using conventional fractionated 1.8 GyE/fraction concurrent PBT with a dosage of 59.4GyE for LAUPC have resulted in a modest decrease in GI toxicity with no grade 3 toxicity during treatment, or during follow-up. These studies had a limitation of small sample size and short median follow up [10, 11].

The clinical target volume (CTV) to PTV set-up margin (SM) alone cannot guarantee the adequate dose coverage of the CTV in PTV-based DS, PBS, and intensity-modulated proton therapy (IMPT) plans [12, 13]. In PTV-based IMPT plans, under- or over dosage inside the PTV can occur in the patient from deviation in the position of high in-field dose gradients from spot to spot due to set-up errors or range uncertainties. Hence, the comparison of PTV-based proton and photon treatment is certainly not precise [14]. The robust IMPT (ro-IMPT) plan can result in even dose gradients per field across the target volume and can reduce the risk of pencil beams ceasing directly in front of an adjoining normal tissue [13]. However, for pancreatic cancer, a dosimetric and radiobiological model-based comparative treatment planning study between ro-IMPT and IMRT has not yet been reported.

For the pancreatic head cancer surrounded circumferentially by gastrointestinal OARs (GI-OARs), IMRT was stated to be superior compared to DS proton therapy [15]. The dosimetric study by Thompson et al., in pancreatic head cancer reported no dosimetric evidence that DS and PBS proton therapy facilitates EDR more readily in comparison to IMRT, as surrounding GI-OARs receive incrementally higher doses using DS and PBS proton therapy [16]. The ro-IMPT with a spot-scanning technique would offer a more-fair comparison with IMRT plans. In-silico study by Stefaowicz et al., using EDR ro-IMPT in advanced pancreatic cancer, have reported a better target homogenous dose distribution and minimized dose to the OARs with a 3 beam design configuration with at least one non-coplanar beam [14].

The dose-volume analysis study is usually restricted to just certain specific DVH parameters that might not always correspond directly to a clinical outcome. The radiobiological normal tissue complication probability (NTCP) model using parameters emanate from reported toxicity rates in clinical trials, and it assesses the treatment plans by analyzing the information from the entire DVH. However, each toxicity endpoint has a specific NTCP parameter set, and besides, it depends on the cohort of the patient and treatment technique used. It is essential to use a more accurate predictor while comparing treatment plans and taking a clinical decision based on dosimetric benefit and absolute NTCP reduction (∆NTCP) [17].

Hence, the research questions of the present in-silico planning comparison study were: 1) In comparison with IMRT, can the dose delivered to GI-OAR for LAUPC of the head be lowered using ro-IMPT? 2) what is the anticipated clinical advantage of this GI-OARs sparing? To answer these questions, we performed an NTCP radiobiological model-based comparison study between IMRT and ro-IMPT for LAUPC of the head with EDR, and we hypothesized that ro-IMPT could reduce GI-OAR toxicity.


The clinicopathological data of patients were reviewed from the hospital’s medical records. With the approval from the Institute Research Ethics Committee (Reference number: 2017–440), for this study, we identified nine locally-advanced pancreatic ductal adenocarcinoma (LAUPC) of the head patients with T4 disease (encasing superior mesenteric artery or celiac axis) as per the 7th edition of the American Joint Committee on Cancer (AJCC) staging manual from 2015 to 2018. The treatment planning computed tomography (CT) in the supine position was obtained for these 9 patients with a 3 mm slice thickness. Each patient was re-planned for IMRT and ro-IMPT.

Target volume and OAR delineation

Target volume and OARs were contoured on IV contrast CT simulation scans. Gross disease and clinically apparent nodes were included as gross tumor volume (GTV); however, the elective nodal region was not included in the target volume [18, 19]. For this study, the GTV to clinical target volume (CTV) was given margin of 0.5 cm, and the CTV was edited at the interface of the GI-OARs, and as per the previously published phase I/II dose-escalation studies, CTV to planning treatment volume (PTV) was given isotropic expansion margin of 0.5 cm as shown in Additional Figure 1 [6]. In this study for reducing the motion to estimate maximum potential benefit, it was supposed that all patients would be treated using breath-hold technique [16, 20].

The OARs were contoured for all patients, and it includes the whole stomach, the duodenum was from pylorus till ligament of Treitz, bilateral kidney, small bowel loops, liver, and spinal cord. The small bowel loops were contoured 2 cm superior-inferiorly to PTV [15]. The whole liver was contoured, including the vessels and intraductal biliary system. The organ contour “Stoduo” was created, which combines stomach and duodenum for comparison with previously published studies [21].

Dose prescription and OARs constraints

The prescription dose was 59.4GyE at 1.8GyE/ fraction in 33 fractions [10, 11]. The proton beam output was modulated with relative biological effectiveness (RBE) of 1.1 [22]. As all tissues are presumed to have nearly the same RBE, the dose stated in GyE is directly in comparison with the photon doses. The planning goal for IMRT and ro-IMPT was at least 100% of GTV receives ≥95% of the dose, at least ≥98% of CTV receives ≥95% of the dose, and 0% volume of CTV receives < 107% of the prescribed dose. Besides, our goal during IMRT was also to provide adequate PTV coverage of at least 95% of PTV receiving 95% of the dose. One physicist designed all IMRT plans, and all ro-IMPT plans were created by another physicist and were checked by two physicians.

The OARs constraints were

For stomach wall, ≤16 cc receive 50GyE, ≤ 10%volume receive 50GyE, ≤15% volume receive 45GyE, and 0.1 cc receive ≤60GyE [21, 23]. For duodenum, ≤45% volume receive 25GyE, 1 cc receive ≤55GyE, and 0.1 cc receive ≤60GyE [9, 24]. For small bowel loops, ≤ 10%volume receive 50GyE, ≤15% volume receive 45GyE, ≤5% volume receive 54GyE, and 0.1 cc volume receive ≤60GyE [23]. For Kidneys, mean dose ≤18GyE and V23GyE < 30%. The mean liver goal was ≤30GyE, V30GyE ≤ 50%, V35GyE ≤ 33%, and 0.1 cc of spinal cord receive <45GyE.

IMRT and ro-IMPT planning, beam configuration, and optimization

For each patient, two plans were generated (IMRT and ro-IMPT). The non-coplanar 6 beam IMRT plan was made using Raystation v6.2 (Raysearch Laboratories, Stockholm, Sweden) treatment planning system (TPS).

Non-coplanar CTV-based robust multifield optimization IMPT plan was made using Eclipse (v15.1) TPS (Varian Medical System, Inc., Palo Alto, CA). All the ro-IMPT plan was delivered using 3 beams, 2 co-planar beams (1350 and 2200), and one non-coplanar beam (gantry at 2700 with couch at 50) as shown in Fig. 1. The two posterior oblique fields were used so that the proton beam minimally intersects with high uncertainty tissue, in particular, the diaphragm and bowel, and is deemed to be more robust against intra-fractional motions. Also, to reduce the dosimetric impact of organ filling and motion uncertainty, a right lateral field was added through the liver [10, 14, 25]. A more detailed description of IMRT and ro_IMPT planning, beam configuration, and optimization is presented in Additional file 1.

Fig. 1
figure 1

Overview of the beam configuration (a and b) and axial, sagittal, and coronal CT slices showing dose distribution of IMRT (c) and ro-IMPT (d) treatment plan. The range of beam direction and couch angle in all patients is given per beam direction. Non-coplanar beam direction is marked in green, and co-planar beam direction in white

Plan evaluation

For IMRT and ro-IMPT plan evaluation, the DVH of targets (GTV and CTV) and OARs were generated on nominal dose distributions. The IMRT and ro-IMPT plans were compared for target homogeneity and conformity. The target volume and OARs DVH physical dose parameters were documented. Homogeneity is defined by the dose distribution consistency of a plan throughout the target volume. The RTOG formula [(D2%-D98%/D50%)] was used to computed HI (Homogeneity index). Where D2%, D98%, and D50% are the dose received by 2, 98, and 50% of the target volume. The conformation number (CN) formula for CTV [(CTV95) [2]/ (CTV*V95) was used to define conformity around the CTV. Where CTV is a target volume, CTV95 is target volume covered by 95% of reference isodose, and V95 is a volume of 95% isodose. As CN value approaches 1, the plan is deemed to be more conformal, and plan with CN 0 indicates the total absence of conformity or a huge volume of irradiation compared to the target volume.

Dose-volume data and Normal tissue complication probability

Coverage of target volume and various dose-volume parameters were assessed. The Digital Imaging and Communications in Medicine (DICOM) standard RT doses from IMRT and ro-IMPT plan were transferred to MIM (v6.86, MIM Software Inc., Cleveland, United States). Before NTCP calculation, the linear-quadratic (LQ) equation with α/β=4 (for the stomach, duodenum, small bowel, and stoduo) was used to convert the cumulative physical dose into an equivalent dose of 2Gy (EQD2) per fraction.

The radiobiological Lyman-Kutcher-Butcher (LKB) model was used to computed NTCP for GI toxicity endpoints using parameters from Pan et al., Burman et al., and Holyoake et al., as shown in Table 1 [26,27,28,29]. Computed NTCP values were used in a relative sense for comparison between ro-IMPT and IMRT. The RADBIOMOD Visual Basic for Application (VBA) software was used to calculate NTCP values from EQD2 DVH’s ASCII files [30]. The absolute NTCP reduction (∆NTCPIMRT – ro-IMPT) and quantitative relative risk (RR = NTCPro-IMPT/NTCPIMRT) ratios for GI-OARs was also computed.

Table 1 Normal tissue complication probability (NTCP) LKB model parameters used in biological evaluation of IMRT and ro-IMPT plans

Statistical analysis

The mean and standard deviation (SD) was used to describe all continuous variables. The non-parametric Wilcoxon sign rank exact test provides an estimate of statistical significance between techniques. Two-sided P-value < 0.05 was considered to be statistically significant. R statistical software version 3.6.3 (R commander EZR version 2-6.2) was used for all statistical analysis.


Target dose parameters evaluation

The patient characteristics are shown in Table 2. The median GTV, CTV, and PTV volumes were 37.6 cc (range, 22.9 to 54.3 cc), 76.1 cc (range, 49.3 to 99.8 cc), and 135.5 cc (range, 93.9 to 178.0 cc) respectively. IMRT and ro-IMPT dose distribution for one representative patient is shown in Fig. 1. Target coverage for all plans met the required goal for the GTV (D100% ≥ 95%) and CTV (V95% ≥ 98%), and the result for target coverage is shown in Fig. 2 (a and b). CTV CN for ro-IMPT plan show significantly (P = 0.004) better conformation of the dose; as a result, a lower percentage of the body outside the CTV was irradiated to high doses with ro-IMPT than with IMRT.

Table 2 Patient characteristics
Fig. 2
figure 2

Comparisons of average cumulative DVH curve for target volume (GTV and CTV), stomach, duodenum, small bowel, and stoduo using IMRT and ro-IMPT plans. Average DVH are shown for each cohort of plans (IMRT and ro-IMPT). Radiation dose is shown along the X-axis and cumulative volume receiving at least dose is plotted on Y-axis

Dose delivered to OARs

Physical dosimetric OAR DVH parameters were significantly lower in the ro-IMPT plan in comparison to IMRT, in low to high dose range (V10GyE to V55GyE) for the stomach, small bowel, and stoduo (Table 3 and Fig. 2 c, e, and f). For stoduo, V50GyE was 7.91 ± 4.4 cc with ro-IMPT vs. 12.9 ± 5.6 cc with IMRT (P = 0.007; Table 3). For duodenum, ro-IMPT delivered a significantly lesser dose in low dose area (≤10Gy) in comparison to IMRT (Table 3 and Fig. 2d).

Table 3 Target volumes and OARs physical dosimetric parameters and comparative analysis between IMRT and ro-IMPT plans

For liver, V35GyE and V30GyE were significantly lower in ro-IMPT in comparison to IMRT. In contrast, the kidneys Dmean and V23GyE were significantly higher with ro-IMPT (Table 3).

For the spinal cord, D0.1cc was 7.27 ± 4.3GyE with ro-IMPT vs. 20.8 ± 1.9GyE with IMRT (P = 0.004; Table 3).

The change in dose to OARs with robustness on the CTV at its worst iteration compared to nominal doses is shown in Additional figures 3, 4, and 5. The dose constraint was met for all OARs with robustness on the CTV in worst-case iteration except for the duodenal constraints, V25GyE ≤ 45% for three patients, D0.1cc ≤ 60Gy for one patient, and D1cc ≤ 55Gy for six patients as shown in Additional figure 3.

NTCP analysis

As reported in Table 4 and Fig. 3, the probability of gastric ulceration/perforation and gastric bleed was significantly worse in the IMRT plans in comparison to ro-IMPT plans according to the model of Pan et al. (stomach, duodenum, and stoduo) and Burman et al. (stomach) [26, 28]. The NTCP value for small bowel was not significantly different in two irradiation techniques (Table 4).

Table 4 Normal tissue complication probability (NTCP), Relative risk (RR) ratio, number of patients with ∆NTCPIMRT – ro-IMPT in specific range for gastro-intestinal OARs toxicity
Fig. 3
figure 3

Box and whisker plot of NTCP (%) comparison for gastrointestinal OARs (Stomach, duodenum, small bowel, and stoduo) using Pan et al., Burman et al., Holyoake et al., LKB model parameters between IMRT and ro-IMPT treatment plans

The ∆NTCPIMRT – ro-IMPT of ≥5 to < 10% was seen for endpoint gastric bleeding of the stomach (3 patients) and stoduo (2 patients) as shown in Table 4 and Fig. 4. The NTCP mode application to GI OARs has demonstrated an overall GI toxicity relative risk reduction (RR < 1) except for endpoint grade ≥ 3 toxicity for duodenum (RR > 1) for all ro-IMPT plans in comparison to IMRT plans (Table 4). According to the considered toxicity endpoint for the stomach, small bowel, and stoduo, the RR values ranged from 0.16 to 0.81 (Table 4).

Fig. 4
figure 4

Bar graph of NTCP reduction (∆NTCPIMRT – ro-IMPT) for GI-OARs (Stomach, duodenum, small bowel, and stoduo) of each patient using Pan et al., Burman et al., Holyoake et al., LKB model parameters between IMRT and ro-IMPT treatment plans


Our study demonstrates a significant GI-OARs sparing benefit using ro-IMPT over IMRT plans in EDR for LAUPC of the head with better target conformation. Clinically acceptable plan with target coverage goal and OAR dose constraint were made for all patient with both the techniques. To our knowledge, this is the first radiobiological model-based comparative study in LAUPC of the head with EDR to assess the potential radiobiological-based clinical implication of ro-IMPT in reducing GI-OARs toxicity.

Dose to stomach (V50GyE ≤ 16 cc), stoduo (V50GyE ≤ 33 cc), and duodenum (V25GyE ≤ 45%, D1cc, and D0.1cc) were below the threshold predicting the low risk of grade ≥ 2 acute GI toxicity and upper GI bleeding in both the radiation delivery technique [9, 21, 24]. Also, in contrast to the study by Thompson et al. and Bouchard et al., the current study demonstrates that the ro-IMPT plan significantly decreases volume receiving intermediate and higher dose for stomach and small bowel, in addition to a decrease in doses < 30 Gy [15, 16]. Using multi-field optimization ro-IMPT in our study compared with single field optimization in Thompson et al., and reducing the beam penumbra by using a small PBS spot size, explains this dissimilarity in intermediate and high doses region. This result corresponds with the study by Jethwa et al., using ro-IMPT in pancreatic cancer, which might facilitate EDR for cases in which OAR are closely surrounding the GTV in almost all the directions [31].

The study published by Thompson et al., reported that the proton therapy, in comparison to IMRT, substantially reduces the dose in low-intermediate dose range [16]. However, the clinical implication of their result is uncertain. In contrast to their study, we performed an absolute NTCP reduction and relative risk (RR) assessment for GI-OARs. In our study, the ∆NTCPIMRT – ro-IMPT of ≥5 to < 10% was seen for the gastric bleeding endpoint of Pan et al., for GI-OARs stomach in 3 patients (33%) and stoduo in 2 patients (22%). The ro-IMPT plans reduce the relative risk of toxicity for the stomach (gastric bleeding, ulceration, and perforation), small bowel (obstruction and perforation), duodenum (gastric bleeding), and stoduo (gastric bleeding) for all patients.

For small bowel, the significant lesser dose volume in high, intermediate, and low dose range for ro-IMPT did not translate into a decrease in NTCP. This shows that, even though a reduction of dose to OAR been excellent, a statistically significant dosimetric difference may not interpret into clinically considerable differences. The use of a radiobiological NTCP model and NTCP-based quantitative relative risk assessment simplifies the task for different planning technique comparisons. It is more robust than DVH parameters for investigation of GI-OARs related toxicity in spite of the uncertainties in NTCP model parameters.

Our present study has several limitations, and the potential limitation is the use of photon-derived tissue NTCP models. To authenticate the results of this study, a large and reliable clinical outcome data are needed. Because such data are lacking, possibly a significant change in model-based toxicity and uncertainties are seen when these radiobiological parameters are used to define the advantage between radiation technique. This can impact the absolute NTCP values, and thus the ∆NTCP. Although clinical validation of these NTCP models was out of the scope of this study, the relative NTCP comparison must be meaningful. The NTCP model selected in our study was generated based on similar patient cohorts and treatment for upper gastrointestinal tumors. Cautious interpretation of these results is essential because it may be affected by model uncertainties.

In LAUPC of the head, while designing a proton beam with the posterior field, the GI-OARs may be positioned inside the distal Bragg peak. Caution should be taken when delivering a dose for LAUPC of the head given its proximity to the duodenum, and the end range uncertainty that must be taken into account is particularly crucial for dose escalation strategy. An increase in effective biological dose in these organs may result in a higher risk of adverse events. In most cases, this increase in dose can be accounted for with small alteration to the physical dose or treating at a lower total physical dose. In certain cases, it is advantageous to add a field to decrease the overall biological effective dose [32]. Our study had not taken into consideration the effectiveness of variable RBE for protons assuming interpatient variability of α/β [33]. A significant uncertainty with the NTCP values and ∆NTCP can occur, as a result of considerable uncertainty with the RBE variation. This uncertainty may cause a substantial increase in dose to the OARs if the OARs are close to the target volume [34].

The target overlap with GI-OARs could restrict EDR for pancreatic cancers. The ITV is generated to account for pancreatic tumor motion with respiration, which could hinder safe EDR, and respiratory gating could benefit such cases. In advanced pancreatic cancer, regional recurrence remains uncommon in comparison to the local and distant recurrence [19]. Hence, in our study, we opted for EDR to a limited CTV without elective nodal irradiation (ENI) to reduce the risk of an adverse event. Inclusion of ENI in our study would have resulted in an increase in the irradiated volume of the stomach, duodenum, small bowel, kidney, liver, and spinal cord using IMRT in comparison to ro-IMPT as demonstrated by Jethwa KR et al., comparing ro-IMPT with VMAT [31].

A comparison of ro-IMPT plans was carried out for nominal dose distributions supposing an idealized patient setup model based on a single CT scan, wherein the anatomical and geometric changes were not taken into consideration. The GI-OARs are an expansible and movable organ, and as a result, determining the accurate dose-volume constraints is quite challenging [35]. Therefore, well-defined image guidance protocol and adaptive treatment strategy are essential during the clinical implementation of ro-IMPT. Despite that, with this approach, further uncertainties on dose distributions are being introduced through deformable image registration [36]. The ro-IMPT plan optimized and evaluated considering the setup and range uncertainties is generally robust for non-rigid anatomical changes visualized on a repeat CT scans [37]. Nevertheless, according to our understanding, the biases of organ motion, positioning, and respiration tend to occur among patients who are treated using both RT techniques, and hence, should not undermine the comparison of the GI-OARs DVHs.

In future studies, for a better comparison of proton and photon plans, proton plans should be calculated, taking variable RBE into accounts [38]. Investigating the use of image registration and fusion algorithm for dose mapping may be necessary to precisely compute the dose delivered to GI-OARs to confirm the eminence of ro-IMPT plan during radiation.

In conclusion, given the smaller sample size and design of our study, ro-IMPT can potentially provide a substantial decrease in GI-toxicity risk for LAUPC of the head in EDR in comparison to IMRT. The quantitative risk evaluation also supports the potential clinical benefit of EDR IMPT for LAUPC of head due to the lower risk of GI morbidity. The result of our study using EDR ro-IMPT should be considered as hypothesis-generating for future clinical trials and research to verify the expected risk reduction in GI toxicity.

Availability of data and materials

The datasets supporting the conclusion of this article are included within the article.



Locally advanced unresectable pancreatic cancer;


Concurrent chemoradiotherapy


Intensity-modulated radiotherapy




Double scatter proton beam therapy


Pencil beam scanning proton therapy


Proton beam therapy


Clinical target volume


Planning treatment volume


Organ at risk


Gray equivalent


Intensity-modulated proton radiotherapy


Robust intensity-modulated proton radiotherapy


Relative biological equivalent


Normal tissue complication probability


NTCP reduction


Dose received by X% of the volume


Volume receiving X-Gy


Dose-volume histogram


Conformation number


Homogeneity index




Probability value


Relative risk


  1. Bray F, Ferlay J, Soerjomataram I, Siegel RL, Torre LA, Jemal A. Global cancer statistics 2018: GLOBOCAN estimates of incidence and mortality worldwide for 36 cancers in 185 countries. CA Cancer J Clin. 2018;68(6):394–424.

    Article  PubMed  Google Scholar 

  2. Li D, Xie K, Wolff R, Abbruzzese JL. Pancreatic cancer. Lancet. 2004;363(9414):1049–57.

  3. Cardenes HR, Moore AM, Johnson CS, et al. A phase II study of gemcitabine in combination with radiation therapy in patients with localized, unresectable, pancreatic cancer: a Hoosier oncology group study. Am J Clin Oncol. 2011;34(5):460–5.

    Article  CAS  PubMed  Google Scholar 

  4. Bittner M-I, Grosu A-L, Brunner TB. Comparison of toxicity after IMRT and 3D-conformal radiotherapy for patients with pancreatic cancer - a systematic review. Radiother Oncol. 2015;114(1):117–21.

    Article  PubMed  Google Scholar 

  5. Polistina F, Costantin G, Casamassima F, et al. Unresectable locally advanced pancreatic cancer: a multimodal treatment using neoadjuvant chemoradiotherapy (gemcitabine plus stereotactic radiosurgery) and subsequent surgical exploration. Ann Surg Oncol. 2010;17(8):2092–101.

    Article  PubMed  Google Scholar 

  6. Ben-Josef E, Schipper M, Francis IR, et al. A phase I/II trial of intensity modulated radiation (IMRT) dose escalation with concurrent fixed-dose rate gemcitabine (FDR-G) in patients with unresectable pancreatic cancer. Int J Radiat Oncol Biol Phys. 2012;84(5):1166–71.

    Article  PubMed  PubMed Central  Google Scholar 

  7. Takatori K, Terashima K, Yoshida R, et al. Upper gastrointestinal complications associated with gemcitabine-concurrent proton radiotherapy for inoperable pancreatic cancer. J Gastroenterol. 2014;49(6):1074–80.

    Article  CAS  PubMed  Google Scholar 

  8. Terashima K, Demizu Y, Hashimoto N, et al. A phase I/II study of gemcitabine-concurrent proton radiotherapy for locally advanced pancreatic cancer without distant metastasis. Radiother Oncol. 2012;103(1):25–31.

    Article  PubMed  Google Scholar 

  9. Kelly P, Das P, Pinnix CC, et al. Duodenal toxicity after fractionated Chemoradiation for Unresectable pancreatic Cancer. Int J Radiat Oncol. 2013;85(3):e143–9.

    Article  Google Scholar 

  10. Nichols RC, George TJ, Zaiden RA, et al. Proton therapy with concomitant capecitabine for pancreatic and ampullary cancers is associated with a low incidence of gastrointestinal toxicity. Acta Oncol (Madr). 2013;52(3):498–505.

    Article  CAS  Google Scholar 

  11. Sachsman S, Nichols RC, Morris CG, et al. Proton Therapy and Concomitant Capecitabine for Non-Metastatic Unresectable Pancreatic Adenocarcinoma. Int J Part Ther. 2014;1(3):692–701.

    Article  Google Scholar 

  12. Langen K, Zhu M. Concepts of PTV and robustness in passively scattered and pencil beam scanning proton therapy. Semin Radiat Oncol. 2018;28(3):248–55.

    Article  PubMed  Google Scholar 

  13. Stuschke M, Kaiser A, Jawad JA, Pöttgen C, Levegrün S, Farr J. Multi-scenario based robust intensity-modulated proton therapy (IMPT) plans can account for set-up errors more effectively in terms of normal tissue sparing than planning target volume (PTV) based intensity-modulated photon plans in the head and neck regi. Radiat Oncol. 2013;8(1):1–5.

    Article  Google Scholar 

  14. Stefanowicz S, Stützer K, Zschaeck S, Jakobi A, Troost EGC. Comparison of different treatment planning approaches for intensity-modulated proton therapy with simultaneous integrated boost for pancreatic cancer. Radiat Oncol. 2018;13(1):228.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  15. Bouchard M, Amos RA, Briere TM, Beddar S, Crane CH. Dose escalation with proton or photon radiation treatment for pancreatic cancer. Radiother Oncol. 2009;92(2):238–43.

    Article  PubMed  Google Scholar 

  16. Thompson RF, Mayekar SU, Zhai H, et al. A dosimetric comparison of proton and photon therapy in unresectable cancers of the head of pancreas. Med Phys. 2014;41(8Part1):081711.

    Article  PubMed  Google Scholar 

  17. Langendijk JA, Lambin P, De Ruysscher D, Widder J, Bos M, Verheij M. Selection of patients for radiotherapy with protons aiming at reduction of side effects: the model-based approach. Radiother Oncol. 2013;107(3):267–73.

    Article  CAS  PubMed  Google Scholar 

  18. Murphy JD, Adusumilli S, Griffith KA, et al. Full-dose gemcitabine and concurrent radiotherapy for Unresectable pancreatic Cancer. Int J Radiat Oncol Biol Phys. 2007;68(3):801–8.

    Article  CAS  PubMed  Google Scholar 

  19. McGinn CJ, Zalupski MM, Shureiqi I, et al. Phase I trial of radiation dose escalation with concurrent weekly full-dose gemcitabine in patients with advanced pancreatic cancer. J Clin Oncol. 2001;19(22):4202–8.

    Article  CAS  PubMed  Google Scholar 

  20. Wong JW, Sharpe MB, Jaffray DA, et al. The use of active breathing control (ABC) to reduce margin for breathing motion. Int J Radiat Oncol Biol Phys. 1999;44(4):911–9.

    Article  CAS  PubMed  Google Scholar 

  21. Nakamura A, Shibuya K, Matsuo Y, et al. Analysis of dosimetric parameters associated with acute gastrointestinal toxicity and upper gastrointestinal bleeding in locally advanced pancreatic cancer patients treated with gemcitabine-based concurrent chemoradiotherapy. Int J Radiat Oncol Biol Phys. 2012;84(2):369–75.

    Article  PubMed  Google Scholar 

  22. Uzawa A, Ando K, Furusawa Y, et al. Biological Intercomparison Using Gut Crypt Survivals for Proton and Carbon-Ion Beams. 2007;48 Accessed June 12, 2019.

  23. Chuong M, Badiyan SN, Yam M, et al. Pencil beam scanning versus passively scattered proton therapy for unresectable pancreatic cancer. J Gastrointest Oncol. 2018;9(4):687–93.

    Article  PubMed  PubMed Central  Google Scholar 

  24. Huang J, Robertson JM, Ye H, Margolis J, Nadeau L, Yan D. Dose–volume analysis of predictors for gastrointestinal toxicity after concurrent full-dose gemcitabine and radiotherapy for locally advanced pancreatic adenocarcinoma. Int J Radiat Oncol. 2012;83(4):1120–5.

    Article  CAS  Google Scholar 

  25. Nichols RC, Huh SN, Prado KL, et al. Protons offer reduced Normal-tissue exposure for patients receiving postoperative radiotherapy for resected pancreatic head Cancer. Int J Radiat Oncol. 2012;83(1):158–63.

    Article  Google Scholar 

  26. 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(1):123–35.

    Article  CAS  PubMed  Google Scholar 

  27. Lyman JT. Complication probability as assessed from dose-volume histograms. Radiat Res Suppl. 1985;8.

  28. Pan C, Dawson L, McGinn C, Lawrence T, Ten Haken R. Analysis of radiation-induced gastric and duodenal bleeds using the Lyman-Kutcher-Burman model. Int J Radiat Oncol. 2003;57(2):S217–8.

    Article  Google Scholar 

  29. Holyoake DLP, Aznar M, Mukherjee S, Partridge M, Hawkins MA. Modelling duodenum radiotherapy toxicity using cohort dose-volume-histogram data. Radiother Oncol. 2017;123(3):431–7.

    Article  PubMed  PubMed Central  Google Scholar 

  30. Chang JH, Gehrke C, Prabhakar R, et al. RADBIOMOD: a simple program for utilising biological modelling in radiotherapy plan evaluation. Phys Medica. 2016;32(1):248–54.

    Article  Google Scholar 

  31. Jethwa KR, Tryggestad EJ, Whitaker TJ, et al. Initial experience with intensity modulated proton therapy for intact, clinically localized pancreas cancer: clinical implementation, dosimetric analysis, acute treatment-related adverse events, and patient-reported outcomes. Adv Radiat Oncol. 2018;3(3):314–21.

    Article  PubMed  PubMed Central  Google Scholar 

  32. Beltran C, Schultz HL, Anand A, Merrell K. Radiation biology considerations of proton therapy for gastrointestinal cancers. J Gastrointest Oncol. 2020;11(1):225–30.

    Article  PubMed  PubMed Central  Google Scholar 

  33. Paganetti H. Relating the proton relative biological effectiveness to tumor control and normal tissue complication probabilities assuming interpatient variability in α/β. Acta Oncol (Madr). 2017;56(11):1379–86.

    Article  Google Scholar 

  34. Paganetti H, Niemierko A, Ancukiewicz M, et al. Relative biological effectiveness (RBE) values for proton beam therapy. Int J Radiat Oncol Biol Phys. 2002;53(2):407–21.

    Article  PubMed  Google Scholar 

  35. Kavanagh BD, Pan CC, Dawson LA, et al. Radiation dose–volume effects in the stomach and small bowel. Int J Radiat Oncol. 2010;76(3):S101–7.

    Article  Google Scholar 

  36. Stützer K, Jakobi A, Haase R, et al. Characterization of deformation vector fields for the registration of dose distribution in adaptive treatment planning and sequential boost protocols. In: Joint Conference of the SSRMP, DGMP, ÖGMP; 2014. p. 140–1. ISBN: 987-3-9816508-5-3.

    Google Scholar 

  37. Korevaar EW, Habraken SJM, Scandurra D, et al. Practical robustness evaluation in radiotherapy – a photon and proton-proof alternative to PTV-based plan evaluation. Radiother Oncol. 2019;141:267–74.

    Article  CAS  PubMed  Google Scholar 

  38. Paganetti H. Relating proton treatments to photon treatments via the relative biological effectiveness - should we revise current clinical practice? Int J Radiat Oncol Biol Phys. 2015;91(5):892–4.

    Article  PubMed  Google Scholar 

Download references


Not applicable.


This study was supported in part by the Japan Agency for Medical Research and Development (AMED) under Grant Number18ck0106210h0003 and the National Cancer Center Research and Development Fund [grant numbers 31-A-17].

Author information

Authors and Affiliations



The collection of the patient’s clinical data was performed by VR. The treatment planning was done by HB and KH. Data analysis was done was VR and HH. Correction and/or improvement were suggested by NN, SZ, AM, HH, HB, KH, RT, HT, KM, AA, MN, MO, YH, TR. Major revision were done by TA. All the authors have read and approved the manuscript and agree to its submission.

Corresponding author

Correspondence to Tetsuo Akimoto.

Ethics declarations

Ethics approval and consent to participate

The study was approved by National Cancer Center Hospital East (NCCHE) Institute Research Ethics Committee (Reference number: 2017–440).

Consent for publication

Not applicable.

Competing interests

The authors declare that they have no competing interest.

Additional information

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Supplementary information

Additional file 1: Figure 1.

Axial, sagittal, and coronal CT image of one representative patient with gross tumor volume (GTV), clinical target volume (CTV), planning treatment volume (PTV), and organs at risk (OARs) contours.

Additional file 2: Figure 2.

Plan uncertainty DVH for one representative patient to quantify robustness of the ro-IMPT treatment plan (Dose coverage in each worst-case scenario met the CTV criteria of V95% ≥ 98% and D0% < 107%). Radiation dose in Gy and % is shown along the X-axis and structure volume (%) along the Y-axis.

Additional file 3: Figure 3.

Duodenal dose with robustness analysis on the CTV, at its worst iteration. In this figure the nominal and worst-case scenario value for the duodenum are shown. Dashed lines represent the dose constraints. (Abbreviations: DXcc: dose received by Xcc volume; VxGy: volume receiving x-Gy).

Additional file 4: Figure 4.

Stomach and small bowel dose with robustness analysis on the CTV. In this figure the nominal and worst-case scenario value for stomach and small bowel are shown. Dashed lines represent the dose constraints. (Abbreviations: DXcc: dose received by Xcc volume; VxGy: volume receiving X-Gy).

Additional file 5: Figure 5.

Kidneys, liver, and spinal cord dose with robustness analysis on the CTV at its worst iteration. In this figure the nominal and worst-case scenario value for Kidneys, liver, and spinal cord are shown. Dashed lines represent the dose constraints. (Abbreviations: Dmean: mean dose; DXcc: dose received by X cc volume; VxGy: volume receiving X-Gy).

Additional file 6.

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit The Creative Commons Public Domain Dedication waiver ( applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Raturi, V.P., Hojo, H., Hotta, K. et al. Radiobiological model-based approach to determine the potential of dose-escalated robust intensity-modulated proton radiotherapy in reducing gastrointestinal toxicity in the treatment of locally advanced unresectable pancreatic cancer of the head. Radiat Oncol 15, 157 (2020).

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: