Next Article in Journal
The Role of ctDNA and Liquid Biopsy in the Diagnosis and Monitoring of Head and Neck Cancer: Towards Precision Medicine
Next Article in Special Issue
The Potential and Challenges of Proton FLASH in Head and Neck Cancer Reirradiation
Previous Article in Journal
Feasibility Study of Nivolumab in Combination with Carboplatin Plus Paclitaxel and Concurrent Thoracic Radiation in Patients with Untreated Unresectable Locally Advanced Non-Small Cell Lung Cancer
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Proton Therapy Adaptation of Perisinusoidal and Brain Areas in the Cyclotron Centre Bronowice in Krakow: A Dosimetric Analysis

1
Cyclotron Centre Bronowice, Institute of Nuclear Physics Polish Academy of Sciences, PL31342 Kraków, Poland
2
National Oncology Institute—National Research Institute, Krakow Branch, PL31115 Kraków, Poland
*
Author to whom correspondence should be addressed.
Cancers 2024, 16(18), 3128; https://doi.org/10.3390/cancers16183128
Submission received: 19 July 2024 / Revised: 6 September 2024 / Accepted: 9 September 2024 / Published: 11 September 2024

Abstract

:

Simple Summary

Adaptive proton therapy (APT) is an evolving approach to proton beam scanning treatment planning. We performed dosimetric study on two groups of head and neck (H&N) patients to evaluate the influence of plan adaptation on planning target volume (PTV) and organs at risk (OARs) doses, resulting from the changes in patient anatomy observed in control computed tomography (CT). The adapted treatment plans, which incorporated the changes observed in the control CT images, statistically improved mostly PTV coverage compared to initial plan. Study shows that applying adaptive procedures in clinical workflow may increased efficiency by controlling the proper irradiation of the treated area for H&N cancer patients.

Abstract

Applying a proton beam in radiotherapy enables precise irradiation of the tumor volume, but only for continuous assessment of changes in patient anatomy. Proton beam range uncertainties in the treatment process may originate not only from physical beam properties but also from patient-specific factors such as tumor shrinkage, edema formation and sinus filling, which are not incorporated in tumor volume safety margins. In this paper, we evaluated variations in dose distribution in proton therapy resulting from the differences observed in the control tomographic images and the dosimetric influence of applied adaptive treatment. The data from weekly computed tomography (CT) control scans of 21 patients, which serve as the basis for adaptive radiotherapy, were used for this study. Dosimetric analysis of adaptive proton therapy (APT) was performed on patients with head and neck (H&N) area tumors who were divided into two groups: patients with tumors in the sinus/nasal area and patients with tumors in the brain area. For this analysis, the reference treatment plans were forward-calculated using weekly control CT scans. A comparative evaluation of organ at risk (OAR) dose-volume histogram (DVH) parameters, as well as conformity and homogeneity indices, was conducted between the initial and recalculated dose distributions to assess the necessity of the adaptation process in terms of dosimetric parameters. Changes in PTV volume after replanning were observed in seventeen patient cases, showing a discrepancy of over 1 cm 3 in ten cases. In these cases, tumor progression occurred in 30% of patients, while regression was observed in 70%. The statistical analysis indicates that the use of the adaptive planning procedure results in a statistically significant improvement in dose distribution, particularly in the PTV area. The findings led to the conclusion that the adaptive procedure provides significant advantages in terms of dose distribution within the treated volume. However, when considering the entire patient group, APT did not result in a statistically significant dose reduction in OARs ( α = 0.05).

1. Introduction

Radiotherapy is typically planned based on a single computed tomography (CT) scan, assuming that the patient’s anatomy and its position relative to the radiation delivery system remain unchanged within predefined limits during the treatment course. Most random uncertainties are addressed by applying safety margins around the target volume. However, certain factors are not entirely considered in the planning target volume (PTV), such as substantial changes in patient anatomy due to variations in tumor volume or fluctuations in organ filling [1,2]. Some of these anatomical changes result from the treatment itself, such as tumor shrinkage, while others may be unrelated to radiation therapy. In general, changes in patient anatomy during proton therapy are likely the most significant source of range uncertainty, as the dose distribution from the proton beam is sensitive to variations in range caused by density changes along the beam path [3].
The goal of adaptive proton therapy (APT) is to ensure the precise delivery of the prescribed dose while adhering to constraints on healthy tissues, particularly if patient geometry undergoes changes throughout a fractionated treatment. There are generally two strategies in adaptive treatment: online and offline adaptation [4]. Online APT can reduce the safety margins by conducting a daily observation of patient anatomy and organ localization. This requires well-verified implementation of online plan reoptimization and validation protocol, which is a clinically difficult approach. The offline adaptation can be implemented using patient control imaging, but the new plan preparation needs additional workload and does not account for the daily patient variations. In general, the APT aims to improve the precision of radiation therapy by systematically monitoring variations in patient anatomy and integrating them into the treatment plan through reoptimization during the course of treatment. A significant challenge in offline plan adaptation is the demanding clinical workflow and the need of additional clinical resources. This process involves imaging, contour definition, plan evaluation, plan adaptation and, ultimately, plan verification [2]. Therefore, the workflow of APT must be optimized to balance potential benefits for the patient against the required resources engaged in this process. Recently, simulation studies were also presented for partial APT, which is an example of the time-optimized solution, but still with the need of clinically approved software for plan preparation and verification, which is not yet available [5].
Different solutions are proposed for APT by research groups and consortiums around the world, such as ProtOnART [6] and RAPTOR [7]. They work on solutions to bring near-real-time APT utilizing Monte Carlo methods, deep learning- or model-based solutions, such as normal tissue complication probability [8,9,10].
In this paper, we present the clinical APT workflow for sinus and brain patients treated at the Cyclotron Centre Bronowice (CCB) Krakow proton therapy center. We hypothesize that APT improves dose distributions delivered to the patient throughout the treatment course. To verify this, we conduct a retrospective dosimetric analysis comparing the dose distributions in the target and organs at risk (OARs), using initial treatment plans recalculated on control CT scans and reoptimized plans after APT.

2. Materials and Methods

At the CCB, protons are generated and accelerated to therapeutic energies up to 230 MeV in the isochronous cyclotron C 230 (IBA, Louvain-la-Neuve, Belgium). Additionally, the CCB is equipped with two rotating gantries featuring modern proton pencil beam scanning (PBS) nozzles. Proton radiotherapy procedures performed at the CCB are conducted in collaboration with the National Oncology Institute (Kraków, Poland).
Approximately 12% of patients treated with protons at the CCB between 2020 and mid-2024 underwent APT. The overall statistics of patients treated at CCB with proton beam scanning by year are included in Table 1. The table also provides information on how many patients underwent the APT procedure each year. The majority of APT was performed for tumors located in head and neck (H&N) area (>90%). Usually, control CT scans are obtained weekly or in response to specific changes in anatomy due to drug administration or weight loss. Follow-up CT scans are further used for APT if necessary. Typically, it takes a few days from acquiring a control CT image to delivering the clinically approved adapted plan. This adaptive approach enhances dose conformity in the case of interfractional anatomical variation. However, this method proves insufficient for non-anatomical patient variations, such as organ cavity fillings appearing after the scan and not being detected in the control CT image.

2.1. Data

In the CCB Krakow proton center, two types of imaging are conducted. Two-dimensional X-ray imaging is performed daily for the purpose of patient positioning based on bony anatomy. However, these 2D images do not provide information on soft tissue and possible anatomy changes that may be an indication for APT. Therefore, the clinical imaging protocol includes weekly 3D kV CT imaging that provides comprehensive information about both bony structures and soft tissues.
Examination of a patient’s anatomical conditions using the CT scanner is conducted by radiation therapists (RTTs) under the supervision of the physician. The CCB employs a 64-slice Siemens Somatom Definition AS Open (Siemens Healthcare, Erlangen, Germany) for performing CT scans. For H&N localization, the CT scan parameters include a 120 kV X-ray tube voltage, a slice thickness of 1.2 mm and a reconstruction diameter, i.e., the field of view (FoV), of 350 mm.
The subsequent step in preparing the treatment plan involves defining the contours of the PTV with appropriate margins, as well as delineating the OARs. The PTV is created by extending the clinical treatment volume (CTV) by 3 mm to accommodate variations in patient anatomy and positioning. It also considers the uncertainty of the CT calibration curve, equivalent to a 3.5% safety margin of the maximal beam range plus 1 mm.
For each patient treated at the CCB, a planning CT is available with corresponding contours outlined by a radiation oncologist. An intensity modulated proton therapy (IMPT) plan is then optimized for each patient using the Eclipse (v16.1, Varian, Palo Alto, CA, USA) treatment planning system. The proton beam specifications are based on the IBA dedicated scanning nozzle beam model, with the spot sigma ( σ ) in air ranging from 2.5 mm to 6.4 mm for the applicable nominal beam energies (between 226.1 MeV and 70 MeV, respectively). Dose calculations and optimization during treatment planning are conducted using proton convolution superposition (PCS) or nonlinear universal proton optimization (NUPO), assuming a constant relative biological effectiveness (RBE) equal to 1.1.
A cohort of 21 H&N patients was chosen for analysis (refer to Table A1 and Table A2 from Appendix A). The selection criteria were based on the PTV’s location, requiring replanning due to dose distribution irregularities in the PTV and/or OAR areas. Patients with new immobilizations prepared prior to control CT were not included. Among the 21 patients, 10 had tumors within the sinus or nasal cavity areas, and 11 had brain tumors. The total prescribed dose ranged from 54 to 70 Gy(RBE), depending on the diagnosis. The plans comprised 2–6 fields.
According to the aforementioned tables, the number of days that elapsed from the CT scan to the start of irradiation using the adaptive plan was determined. The average replanning time was 5 (±2) days for patients with tumors in the sinus area and 6 (±4) days for patients with brain tumors. For the entire patient cohort, the average was 5 (±3) days.

2.2. Adaptation Protocol

Patients may experience various sources of geometric variation between fractions throughout the course of radiotherapy. Changes in patient anatomy, whether related to treatment or not, may have a dosimetric impact on OARs or lead to dosimetric discrepancies within treatment volume coverage, highlighting the importance of systematic anatomy monitoring [2].
If patient anatomy variations occur, anatomy monitoring may be followed by APT that, however, introduces an additional workload for the staff, mainly because of the necessity of replanning and patient-specific plan quality assurance (QA), and for radiation oncologists who review and approve the plans. Generally, offline adaptation has recently been employed in proton therapy facilities [11].
An offline adaptation is also implemented at the CCB. To preliminarily assess the necessity of APT, a rigid fusion of images from the control scans with the initial scan used for treatment planning is conducted. Changes between the specified scans are analyzed by radiation oncologists and medical physicists. The visual inspection of fused images includes possible variations in the quality of coverage of the treated volume and the exposure of critical organs that may be followed by dose recalculation and evaluation. Exemplary image and contour comparisons are shown in Figure 1.

2.3. Plan Recalculation and Data Analysis

To evaluate the dosimetric benefit of APT, the initial treatment plans were forward-calculated on the control CT scans. The analysis involved comparing the originally calculated dose distribution on the control CT scans with the corresponding APT dose distribution. The assessment of all IMPT plans relied on the calculated dose-volume histograms (DVHs) and their associated metrics, i.e., DVH indices. The procedure is illustrated schematically in Figure 2.
The APT aims to achieve the clinical goals established during the development of the initial treatment plan by replanning on the control CT. These goals can be specified for any designated volume, including PTVs and OARs. The comparison of DVH parameters was conducted between adapted plans (control CT, new plan/replan) and non-adapted plans (control CT, initial plan) to illustrate the dosimetric benefits of the adaptations performed in Krakow.
Typically, various constraints for different organs are defined, such as the dose received by percentage of the organ volume (i.e., D 98 % ) or the volume of the organ receiving at least a specific value of the absorbed dose expressed in gray (Gy) or percentage of the prescribed dose (i.e., V 95 % or V 54 Gy ). Given the location of the tumors in the selected patient group, structures such as the brain, brainstem, chiasm, optic nerves and spinal cord were chosen for analysis. As these structures correspond to serial organs (organs where disabling any subunit causes the entire organ to fail), parameters such as D max and D near max (dose delivered to 0.1 cm 3 of a structure volume) were selected for analysis.
Additionally, parameters of the PTV and metrics describing the quality of coverage of the planned volume were analyzed, including conformity index (CI), homogeneity index (HI) and conformation number (CN):
  • The homogeneity index (HI) signifies the slope of the DVH, reflecting the uniformity of the dose distribution within the PTV. The HI is calculated as follows:
    HI = D 2 % D 98 % D p × 100
    where D 2 % and D 98 % represent the doses received by 2% and 98%, respectively, of the PTV, and D p is the prescribed dose [12]. The ideal value is zero, indicating equality between D 2 % and D 98 % .
  • The conformity index (CI) describes the coverage of the PTV with the prescribed dose. The CI is expressed by the following Formula (2):
    CI = TV RI TV
    where TV RI denotes the target volume covered by the reference isodose (e.g., at least 98% of the prescribed dose D p ), and TV is the target volume. This index varies from 0 (when the entire target volume is located outside 98% D p isodose) to 1 (when the entire target volume receives at least 98% D p ).
  • The conformation number (CN) simultaneously considers the irradiation of the target volume and healthy tissues and is defined as follows:
    CN = TV RI TV · TV RI V RI
    where TV is the target volume, V RI is the volume of the reference isodose (98%). The CN varies from 0 (indicating either a complete lack of conformation or a very large volume of irradiation compared to the target volume) to 1 (in cases when dose coverage is ideal) [13].
In this analysis, eight parameters were assessed [14], with two parameters focusing on OARs in the H&N region ( D near _ max (0.1 cm 3 ), D max ). The remaining six parameters characterized the dose distribution and conformity within the PTV area ( D mean , D 98 % , V 95 % , CI, HI and CN). A summary of the parameter descriptions is provided in Table 2.
A nonparametric statistical hypothesis test, specifically the Wilcoxon signed-rank test, was conducted to compare the parameters of the initial clinical plans before and after the adaptive procedure.

3. Results

Table A1 and Table A2 in Appendix A contain a list of patients selected for the study. According to this table, changes in PTV after replanning were observed in seventeen cases, showing a discrepancy of over 1 cm 3 in ten cases. Each of the investigated patients underwent adaptation only one time. For a change in PTV volume of less than 1 cm 3 , it was assumed that it was a consequence of copying structures onto other CT scans. In cases where the PTV volume exceeded 1 cm 3 , tumor progression occurred in 30% of patients, while regression was observed in 70%. The most significant change in PTV volume was observed in patient #6 of the sinus tumor group (Table A1), reaching 21.4%. In 81% of cases, the motivation behind APT preparation was to enhance the quality of the dose distribution, providing patients with dosimetric benefits such as limiting the dose to critical organs or achieving more optimized PTV dose coverage. Interestingly, in one case, an adaptive plan was devised despite meeting the treatment plan constraints, indicating that APT still improved PTV coverage. The analysis of the results for the patient group reveals that 33.3% of cases necessitated an adaptive plan in the first half of the irradiation cycle, and the remaining 66.7% necessitated an adaptive plan in the second half of the cycle. Table 3 summarizes the timing of the start of irradiation using the adaptive plan during the radiotherapy course for two tumor location groups: sinus/nasal cavity and brain.
According to Table 3 in the sinus/nasal cavity tumor location patient group, irradiation of the adapted plan started at a random week during the treatment course. The median replanning week for this group is week 4, with a range spanning from week 1 to week 7.
In the brain tumor location patient group, there were no patients who began the irradiation of the new plan beyond week 5. The median replanning week for this group is week 5, ranging from week 1 to week 5. Notably, 6 out of 11 brain tumor patients, i.e., over 50%, were replanned during week 5. The results indicate that the more critical period from the point of view of adaptive radiotherapy in CCB is the second half of the treatment cycle, regardless of the group of patients analyzed.
Table 4 and Table 5 present the results regarding the occurrence of improvements or deteriorations in dose distribution parameters. It is noteworthy that a significant improvement in dose distribution was observed in the PTV area, with enhancements seen in 80% or more of cases for each PTV parameter. The results are presented in two separate tables, distinguishing between patients with tumors in the sinus area and those with brain tumors. For the group with brain tumors, the dose assessment within the brain structure was omitted. This was due to the fact that the medical physicist did not include a brain structure in the excision of the PTV area during the planning process, meaning that any assessment would also have included the PTV area. Note that some of the OAR were not delineated during the plan preparation. For brain tumor irradiation, the spinal cord was not delineated in most cases (54.5%), resulting in lower quality statistical analysis. In the case of proton radiotherapy, we are dealing with a much lower level of scattered radiation. Irradiation is more conformal, and consequently, organs located several centimeters from the target area are practically unaffected (if they are not in the path of the proton beam). Consequently, the physician decides not to contour some of the structures.
Table A3 in Appendix A shows a catalog of OAR and PTV parameters assessed in the analysis for the entire group of patients (without separating sinus and brain tumors), along with the number of patients exhibiting improvement, worsening or no change, considering the uncertainty in parameter determination. An analysis of Table A3 shows a clear statistically significant improvement in the PTV area, for all parameters analyzed.
The results obtained for the CI, HI and CN indicate improvement resulting from the APT. A statistical hypothesis test for the parameters specified for the PTV demonstrates that the preparation of an adaptive plan enhances the dose distribution both in the treated area (PTV) with the statistical significance at α = 0.05. However, the statistical hypothesis test for the parameters specified for selected OARs leads to the conclusion that there is no statistically significant evidence, at α = 0.05. A qualitative evaluation of the enhancement in dose distributions within the PTV area is most conveniently observed through the block diagrams presented in Figure 3, Figure 4 and Figure 5.
Based on the results presented in the graphs in Figure 3, it is apparent that the adaptive procedure has played a pivotal role in elevating the quality of the treatment plan. The mean CI parameter exhibited an improvement of approximately 15% across the entire patient group post the adaptive procedure. Particularly noteworthy is the more substantial improvement observed in patients with tumors in the sinus and pharynx region, amounting to approximately 28%. In contrast, patients with brain tumors demonstrated an improvement of about 4%.
According to Figure 4, the average HI parameter for the entire patient group showed an improvement of approximately 22% following the adaptive procedure. Similarly, patients with tumors in the sinus and pharynx area demonstrated a more significant improvement, reaching about 35%. An improvement of approximately 10% was in a group of patients with brain tumors.
Based on Figure 5, the average value of the CN parameter for the entire examined patient group increased by approximately 24%. In the case of patients with tumors in the sinuses and nasal cavities, the improvement in the abovementioned average value reached 44%, while for patients with brain tumors, there was an improvement of about 5%. A significant reduction (>30%) in the standard deviation was clearly observed for the group of patients with tumors in the sinus and nasal area. This indicates that the quality of the dose distribution was similar for each of the analyzed cases.
Figure 6 shows the feature boxplots for various parameters within the PTV range, encompassing the mean dose, D 98 % , and V 95 % , specifically for the two distinct groups of patients.
The obtained results lead to the conclusion that in a significant majority of the analyzed cases, the adaptive procedure enables the restoration of dose parameters to levels closely resembling the nominal values for the PTV. As depicted in Figure 6, the recalculation of the original plan to the control CT leads to a deterioration in dose distribution in both analyzed patient groups. Significantly greater differences were noted in the group of patients with tumors located in the sinus and nasal cavities. Within this particular patient cohort, anatomical changes and variations in sinus filling had a more pronounced impact on the quality of the dose distribution.

4. Discussion

In this study, we evaluated the dosimetric impact of a clinically applied APT protocol on 21 H&N tumor patients, who represent 11.7% of the 116 patients treated with APT at the CCB proton therapy center between 2020 and mid-2024. Most patients with APT have tumors located in the H&N area (∼90%). APT led to improved target coverage and either similar or slightly reduced doses to OARs. When no plan adaptation was performed, target coverage fell below clinical goals, but this was restored after the application of APT. A single APT was, in most cases, sufficient to restore target coverage. Based on the available data, the second half of the treatment cycle is the more critical period for adaptive radiotherapy at CCB, regardless of the patient group analyzed.
The findings obtained in this research are in line with results from other institutions [15,16,17]. One of the strengths of this study is the number of patients that have undergone retrospective analysis, as in the majority of publications, the number of analyzed cases did not exceed 10 patients [16]. However, gathering a sufficient number of replanning cases for specific indications and treatment planning protocols—such as tumor size, location and number of fields—remains challenging. To confirm the findings presented here, we are actively extending the existing database, prospectively collecting data of patients that have undergone APT. Unifying the APT protocols and combining data from multiple institutions would further reinforce the evidence of APT. Moreover, this study indicates that the motivation for preparing an adaptive plan originated from both the disturbance in the dose distribution in the target and OARs [1,18]. Our data show that adapting the dose distribution in the target area typically improves dose distribution in the OARs region. This is contradictory to the majority of publications on APT, where the primary trigger for APT was the deterioration of the dose in the target volume only [15,17].
The great advantage of the presented analysis is that the presented data were used for APT at an operating proton therapy facility. Our findings show that the replanning procedure typically takes about five days. The current protocol involves multiple steps: medical decision-making, contour delineation, treatment plan reoptimization and dosimetric verification. While the protocol could be enhanced by incorporating AI to support medical decisions, automating contouring [19,20], replanning in the TPS [21] and using Monte Carlo-based dosimetric verification [22,23], it reflects legal requirements that do not differentiate between original treatment planning and treatment plan adaptation. Currently, the collaborative effort of the RAPTOR consortium is dedicated to advancing proton therapy by developing and implementing automative, online treatment plan adaptation strategies to improve precision and outcomes in cancer treatment. Online APT could offer substantial advantages. As indicated in [16], online daily APT allows for the quickest response to changes in patient anatomy. As it stands, APT at the CCB Krakow proton center demands significant manpower and technical resources, which must be weighed against the clinical benefits for the patient. Notably, even with offline APT, the protocol still leads to dosimetric improvements over no APT, validating the current clinical approach.
One crucial aspect of APT is the frequency of control imaging. More frequent control imaging would help monitor daily anatomical changes and assess target dose deterioration, aiding in the optimal timing for offline APT. This approach could also determine if and when online APT is advantageous, considering the patient’s exposure to daily imaging doses. Significant anatomical changes with high variability are most likely in tumors located in perisinusoidal areas. In contrast, tumors in locations such as the spine, brain or sacrum may experience less dynamic changes due to the radiation therapy process itself. Therefore, a thorough analysis of the benefits and drawbacks of both online and offline APT, especially in relation to patient imaging dose, is essential.
In summary, the shortcomings of this study, which are planned to be considered in the future, relate either to patient data collection and analysis or the clinical APT protocol that is currently applied in the CCB Krakow proton center. In terms of this study, increasing patient statistics would reduce the uncertainty of the results and could provide much more sound conclusions that may be the basis for correction of the existing APT protocol. As discussed above, increased patient statistics would enable the division of the patient cohort into smaller, more concise sub-groups, e.g., according to tumor size, location, tumor type, treatment dose and the number of irradiation fields, frequency of control imaging or patient age, etc. Addressing these shortcomings and expanding the patient dataset will be crucial for refining the APT protocol, leading to more precise and individualized treatment strategies in proton therapy.

5. Conclusions

We studied the treatment plan adaptation protocol applied clinically at the CCB Krakow proton therapy center for a cohort of 21 selected patients treated between 2020 and mid-2024. We evaluated dose distributions before and after plan adaptation using original and control CT scans. The examination of individual DVH parameters and the conducted statistical tests collectively affirm that the quality of the adapted plans has improved compared to the initial plans, especially regarding the coverage of the PTV area. While the reduction in the dose delivered to critical organs was a common objective of the adaptation procedure, the statistical tests indicate no significant differences in critical organ exposure between the initial and adapted plans. This study justifies the application of the offline APT in the CCB Krakow proton therapy center.

Author Contributions

Conceptualization, M.R., K.K., A.S. and J.G.; Data curation, T.S. and R.K.; Formal analysis, M.R., A.S. and M.L.; Funding acquisition, A.R.; Investigation, M.R., K.K., M.G., M.L., G.F., E.G. and D.K.; Methodology, M.R., T.S. and K.K.; Project administration, A.R.; Resources, T.S., K.K., D.K. and R.K.; Software, M.R. and A.S.; Supervision, T.S., D.K., R.K. and A.R.; Validation, M.R. and A.S.; Visualization, M.R. and M.G.; Writing—original draft, M.R.; Writing—review and editing, T.S., K.K., A.S., M.G., M.L., G.F., E.G., J.G., D.K., R.K. and A.R. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the Polish Ministry of Education and Science under grant No. NdS/544748/2021/2021.

Institutional Review Board Statement

Human Ethics and Consent to Participate declarations: not applicable. Ethical review and approval were waived for this study, due to the fact that analyzed patient data were collected from standard treatment datasets and anonymized. At the moment of admission to our institute, every patient gives written consent for the usage of his medical records for scientific studies.

Informed Consent Statement

Informed consent was obtained from all subjects involved in the study.

Data Availability Statement

The datasets generated and/or analyzed during the current study are accessible from the main author upon reasonable request, contingent on the consent of all co-authors, and given that there are no additional costs associated with data sharing.

Conflicts of Interest

The authors declare no conflicts of interest.

Abbreviations

The following abbreviations are used in this manuscript:
APTadaptive proton therapy
CCBCyclotron Centre Bronowice
CTcomputed tomography
CTVclinical target volume
DVHdose-volume histogram
H&Nhead and neck
IMPTintensity modulated proton therapy
OARorgan at risk
PBSpencil beam scanning
PTVplanning target volume
RBErelative biological effectiveness
RTTradiation therapist

Appendix A

Table A3 shows a catalog of OAR and PTV parameters assessed in the analysis for the entire group of patients (without separating sinus and brain tumors), along with the number of patients exhibiting improvement, worsening, or no change, considering the uncertainty in parameter determination. The final column presents the results of the nonparametric statistical hypothesis test.
Table A1 and Table A2 listing patients chosen for the study on proton therapy adaptation of perisinusoidal and brain areas at the Cyclotron Centre Bronowice in Krakow divided into two groups, along with tumour locations and treatment plan details.
Table A1. List of patients with tumors located in the sinus/nasal area selected for the study, including treatment plan specifications. In the column listing the number of fractions completed for the initial plan, information was added in parentheses indicating the week of treatment when the adaptive plan was introduced. Radiation therapy is performed daily, 5 days a week.
Table A1. List of patients with tumors located in the sinus/nasal area selected for the study, including treatment plan specifications. In the column listing the number of fractions completed for the initial plan, information was added in parentheses indicating the week of treatment when the adaptive plan was introduced. Radiation therapy is performed daily, 5 days a week.
Patient No.AgeTumour
Location
Days after
CT Reskan
to Adapted Plan
Original PTV
[cm3]
Adapted PTV
[cm3]
No.
of Fields
Fraction Dose
[GyRBE]
No of Fractions
with Original
Plan/All Delivered
Fractions
Replan
Motivation
Radiotherapy
Combined with:
138Sinus cavity431.6631.2931.819/30 (week 4)tumor regressionrecurrence
276Nasal cavity7211.11211.556224/35 (week 5)dose
increase in OARs
immunotherapy
341Nasal cavity360.4860.116216/33 (week 4)dose
increase in OARs
surgery
435H&N
Nasal cavity
+lymph nodes
6274.95235.354226/35 (week 6)tumor
regression
chemotherapy
543Sinus cavitythe same CT305.45370.886214/35 (week 3)tumor
progression
chemotherapy
669H&N
Sinus cavity
7850.81923.88628/30 (week 2)tumor
progression
chemotherapy
759Sinus cavity5165.84165.3561.831/37 (week 7)dose decrease
in PTV
chemotherapy
870H&N
Sinus cavity
+lymph nodes
4107.66107.936227/35 (week 6)disturbance
of dose distribution
in PTV nad OARs
chemotherapy
957H&N.
Sinus cavity
the same CT453.15453.154/5 *22/27 (week 1)dhange of the
dose constrains
for OAR
surgery
and chemotherapy
1045H&N
Nasal cavity
+lymph nodes
the same CT197.32197.32622/35 (week 1)dhange of the
dose constrains
for OAR
chemotherapy
* 4 fields in initial plan, 5 fields in APT.
Table A2. List of patients with tumors located in the brain selected for the study, including treatment plan specifications. In the column listing the number of fractions completed for the initial plan, information was added in parentheses indicating the week of treatment when the adaptive plan was introduced. Radiation therapy is performed daily, 5 days a week.
Table A2. List of patients with tumors located in the brain selected for the study, including treatment plan specifications. In the column listing the number of fractions completed for the initial plan, information was added in parentheses indicating the week of treatment when the adaptive plan was introduced. Radiation therapy is performed daily, 5 days a week.
Patient No.AgeTumour LocationDays after
CT Reskan
to Adapted Plan
Original PTV
[cm3]
Adapted PTV
[cm3]
No.
of Fields
Fraction Dose
[GyRBE]
No of Fractions
with Original
Plan/All Delivered
Fractions
Replan
Motivation
Radiotherapy
Combined with:
137Brainthe same CT148.46148.4641.82/30 (week 1)Dose increase
in OARs
surgery
221Brain6588.55589.2161.82/33 (week 1)weight gaining
—disturbance of dose
distribution in OARs
surgery
360Brain6352.13350.1731.824/30 (week 5)tumor
regression
surgery
438Brain1152.24150.7721.822/30 (week 5)tumor
regression
surgery
530Brainthe same CT663.92663.9231.86/30 (week 2)dhange of the
dose constrains
for OAR
surgery
648Brain14232.24231.2031.822/30 (week 5)disturbance of dose
distribution in OARs
surgery
744Brain6365.21363.8121.823/33 (week 5)increased Dmax
to Brainstem
-
871Brain2132.45131.7731.822/30 (week 5)disturbance of dose
distribution in CTV
-
940Brain7818.47834.1341.819/33 (week 4)tumor
progression
surgery
1034Brain6539.8537.3131.822/30 (week 5)increased Dmax
to Brainstem
surgery
1119Brain3673.84671.5131.822/30 (week 4)increased Dmax
to Brainstem and Chiasm
surgery
Table A3. List of OAR and PTV parameters evaluated simultaneously for both groups of patients.
Table A3. List of OAR and PTV parameters evaluated simultaneously for both groups of patients.
OAR or PTV
Parameter
No. of PatientsWilcoxon
Signed-Rank Test
ImprovingWorseningNo
Significant
Change
p-ValueTest
Result *
Brainstem D max 9660.179
D near _ max 10560.038
Chiasm D max 7950.355
D near _ max 10650.313
Optic Nerve L D max 8740.277
D near _ max 8650.421
Optic Nerve R D max 11530.395
D near _ max 8650.595
PTV D mean 15150.003
D 98 % 17400.007
V 95 % 17310.001
CI16230.001
CN15240.002
HI18300.001
* ✓—statistically significant improvement, ✗—no statistical evidence of improvement.

References

  1. Stützer, K.; Jakobi, A.; Bandurska-Luque, A.; Barczyk, S.; Arnsmeyer, C.; Löck, S.; Richter, C. Potential proton and photon dose degradation in advanced head and neck cancer patients by intratherapy changes. J. Appl. Clin. Med. Phys. 2017, 18, 104–113. [Google Scholar] [CrossRef] [PubMed]
  2. Paganetti, H.; Botas, P.; Sharp, G.C.; Winey, B. Adaptive proton therapy. Phys. Med. Biol. 2021, 66, 22TR01. [Google Scholar] [CrossRef] [PubMed]
  3. Lomax, A.J. Myths and realities of range uncertainty. Br. J. Radiol. 2020, 93, 20190582. [Google Scholar] [CrossRef] [PubMed]
  4. Bobić, M.; Lalonde, A.; Nesteruk, K.P.; Lee, H.; Nenoff, L.; Gorissen, B.L.; Bertolet, A.; Busse, P.M.; Chan, A.W.; Winey, B.A.; et al. Large anatomical changes in head-and-neck cancers–A dosimetric comparison of online and offline adaptive proton therapy. Clin. Transl. Radiat. Oncol. 2023, 40, 100625. [Google Scholar] [CrossRef]
  5. Gambetta, V.; Fredriksson, A.; Menkel, S.; Richter, C.; Stützer, K. The partial adaptation strategy for online-adaptive proton therapy: A proof of concept study in head and neck cancer patients. Med. Phys. 2024, 51, 5572–5581. [Google Scholar] [CrossRef]
  6. ProtOnART—A New Consortium for Proton Online Adaptive Radiation Therapy. Available online: https://www.raysearchlabs.com/media/press-releases/2022/protonart–a-new-consortium-for-proton-online-adaptive-radiation-therapy/ (accessed on 24 August 2024).
  7. Real-Time Adaptive Particle Therapy of Cancer (RAPTOR). Available online: https://raptor-consortium.com/ (accessed on 24 August 2024).
  8. Galapon, A.V., Jr.; Thummerer, A.; Langendijk, J.A.; Wagenaar, D.; Both, S. Feasibility of Monte Carlo Dropout-Based Uncertainty Maps to Evaluate Deep Learning-Based Synthetic CTs for Adaptive Proton Therapy. Large anatomical changes in head-and-neck cancers–A dosimetric comparison of online and offline adaptive proton therapy. Med. Phys. 2024, 51, 2499–2509. [Google Scholar] [CrossRef]
  9. Smolders, A.; Lomax, A.; Weber, D.C.; Albertini, F. Deep Learning Based Uncertainty Prediction of Deformable Image Registration for Contour Propagation and Dose Accumulation in Online Adaptive Radiotherapy. Med. Phys. 2023, 68, 245027. [Google Scholar] [CrossRef]
  10. Smolders, A.; Hengeveld, A.C.; Both, S.; Wijsman, R.; Langendijk, J.A.; Weber, D.C.; Lomax, A.J.; Albertini, F.; Guterres Marmitt, G. Inter- and Intrafractional 4D Dose Accumulation for Evaluating ΔNTCP Robustness in Lung Cancer. Radiother. Oncol. 2023, 182, 109488. [Google Scholar] [CrossRef]
  11. Zhang, Y.; Alshaikhi, J.; Amos, R.A.; Lowe, M.; Tan, W.; Bär, E.; Royle, G. Improving workflow for adaptive proton therapy with predictive anatomical modelling: A proof of concept. Radiother. Oncol. 2022, 173, 93–101. [Google Scholar] [CrossRef]
  12. Kataria, T.; Sharma, K.; Subramani, V.; Karrthick, K.P.; Bisht, S.S. Homogeneity Index: An objective tool for assessment of conformal radiation treatments. J. Med. Phys. 2012, 37, 207–213. [Google Scholar] [CrossRef]
  13. Feuvret, L.; Noël, G.; Mazeron, J.J.; Bey, P. Conformity index: A review. Int. J. Radiat. Oncol. Biol. Phys. 2006, 64, 333–342. [Google Scholar] [CrossRef] [PubMed]
  14. Noël, G.; Antoni, D. Organs at risk radiation dose constraints. Cancer/Radiothérapie 2022, 26, 59–75. [Google Scholar] [CrossRef]
  15. Lalonde, A.; Bobić, M.; Winey, B.; Verburg, J.; Sharp, G.C.; Paganetti, H. Anatomic changes in head and neck intensity-modulated proton therapy: Comparison between robust optimization and online adaptation. Radiother. Oncol. 2021, 159, 39–47. [Google Scholar] [CrossRef]
  16. Huiskes, M.; Astreinidou, E.; Kong, W.; Breedveld, S.; Heijmen, B.; Rasch, C. Dosimetric impact of adaptive proton therapy in head and neck cancer—A review. Clin. Transl. Radiat. Oncol. 2023, 39, 100598. [Google Scholar] [CrossRef]
  17. Bobić, M.; Lalonde, A.; Sharp, G.C.; Grassberger, C.; Verburg, J.M.; Winey, B.A.; Lomax, A.J.; Paganetti, H. Comparison of weekly and daily online adaptation for head and neck intensity-modulated proton therapy. Phys. Med. Biol. 2021, 66, 055023. [Google Scholar] [CrossRef] [PubMed]
  18. De Ornelas, M.; Xu, Y.; Padgett, K.; Schmidt, R.M.; Butkus, M.; Diwanji, T.; Luciani, G.; Lambiase, J.; Samuels, S.; Samuels, M.; et al. CBCT-based adaptive assessment workflow for intensity modulated proton therapy for head and neck cancer. Int. J. Part. Ther. 2021, 7, 29–41. [Google Scholar] [CrossRef] [PubMed]
  19. Taku, N.; Wahid, K.A.; van Dijk, L.V.; Sahlsten, J.; Jaskari, J.; Kaski, K.; Fuller, C.D.; Naser, M.A. Auto-detection and segmentation of involved lymph nodes in HPV-associated oropharyngeal cancer using a convolutional deep learning neural network. Clin. Transl. Radiat. Oncol. 2022, 36, 47–55. [Google Scholar] [CrossRef]
  20. Nikolov, S.; Blackwell, S.; Zverovitch, A.; Mendes, R.; Livne, M.; De Fauw, J.; Patel, Y.; Meyer, C.; Askham, H.; Romera-Paredes, B.; et al. Clinically applicable segmentation of head and neck anatomy for radiotherapy: Deep learning algorithm development and validation study. J. Med. Internet Res. 2021, 23, e26151. [Google Scholar] [CrossRef]
  21. Bernatowicz, K.; Geets, X.; Barragan, A.; Janssens, G.; Souris, K.; Sterpin, E. Feasibility of online IMPT adaptation using fast, automatic and robust dose restoration. Phys. Med. Biol. 2018, 63, 085018. [Google Scholar] [CrossRef]
  22. Ates, O.; Pirlepesov, F.; Zhao, L.; Hua, C.H.; Merchant, T.E. Development of a log file analysis tool for proton patient QA, system performance tracking, and delivered dose reconstruction. J. Appl. Clin. Med. Phys. 2023, 24, e13972. [Google Scholar] [CrossRef]
  23. Gajewski, J.; Garbacz, M.; Chang, C.W.; Czerska, K.; Durante, M.; Krah, N.; Krzempek, K.; Kopeć, R.; Lin, L.; Mojżeszek, N.; et al. Commissioning of GPU–accelerated Monte Carlo code FRED for clinical applications in proton therapy. Front. Phys. 2021, 8, 567300. [Google Scholar] [CrossRef]
Figure 1. Illustration of changes in sinus filling: (a) treatment planning CT overlaid with the initial dose distribution and delineated PTV (red line), (b) fusion of the control CT images with the initial PTV outline (red line), (c) illustration of a fusion of the initial and control CT images with a color scale highlighting differences: red indicates the greatest difference, while blue signifies no differences.
Figure 1. Illustration of changes in sinus filling: (a) treatment planning CT overlaid with the initial dose distribution and delineated PTV (red line), (b) fusion of the control CT images with the initial PTV outline (red line), (c) illustration of a fusion of the initial and control CT images with a color scale highlighting differences: red indicates the greatest difference, while blue signifies no differences.
Cancers 16 03128 g001
Figure 2. Schematic illustration of the adaptive plan preparation procedure. The yellow box represents the initial plan developed based on the initial CT scans. The blue box illustrates the fusion of the initial and control CT scans. The purple box represents an initial dose distribution recalculated on the control CT scans. The green box represents the newly reoptimized dose calculated on the control CT scans.
Figure 2. Schematic illustration of the adaptive plan preparation procedure. The yellow box represents the initial plan developed based on the initial CT scans. The blue box illustrates the fusion of the initial and control CT scans. The purple box represents an initial dose distribution recalculated on the control CT scans. The green box represents the newly reoptimized dose calculated on the control CT scans.
Cancers 16 03128 g002
Figure 3. Boxplots illustrating the CI following the implementation of the adaptation procedure. The boxes extend from the Q1 to Q3 quartile values of the data, with a line at the median value. Red boxes signify the outcomes derived from recalculating the initial plan on control scans. Green boxes portray the outcomes achieved after adapting the plan to the current anatomical conditions of the patient (replan); (a) CI results for the entire patient group, (b) CI results for patients exclusively with brain tumors, (c) CI results for patients with tumors in the sinus and nasal cavities.
Figure 3. Boxplots illustrating the CI following the implementation of the adaptation procedure. The boxes extend from the Q1 to Q3 quartile values of the data, with a line at the median value. Red boxes signify the outcomes derived from recalculating the initial plan on control scans. Green boxes portray the outcomes achieved after adapting the plan to the current anatomical conditions of the patient (replan); (a) CI results for the entire patient group, (b) CI results for patients exclusively with brain tumors, (c) CI results for patients with tumors in the sinus and nasal cavities.
Cancers 16 03128 g003
Figure 4. Boxplots illustrating the changes in HI following the implementation of the adaptation procedure. The boxes extend from the Q1 to Q3 quartile values of the data, with a line at the median value. Red boxes represent the outcomes obtained after recalculating the original plan on control scans. Green boxes represent the results obtained after adapting the plan to the patient’s current anatomical conditions (replan); (a) HI results for the entire patient group, (b) HI results for a group of patients with brain tumors only, (c) HI results for a group of patients with sinus and nasal cavities tumors.
Figure 4. Boxplots illustrating the changes in HI following the implementation of the adaptation procedure. The boxes extend from the Q1 to Q3 quartile values of the data, with a line at the median value. Red boxes represent the outcomes obtained after recalculating the original plan on control scans. Green boxes represent the results obtained after adapting the plan to the patient’s current anatomical conditions (replan); (a) HI results for the entire patient group, (b) HI results for a group of patients with brain tumors only, (c) HI results for a group of patients with sinus and nasal cavities tumors.
Cancers 16 03128 g004
Figure 5. Boxplots illustrating the changes in the CN after applying the adaptation procedure was applied. The boxes extend from the Q1 to Q3 quartile values of the data, with a line at the median value. Red boxes represent the results obtained after recalculating the initial plan into control scans. Green boxes represent the results obtained after adapting the plan to the patient’s current anatomical conditions (replan); (a) the results of the CN for the whole group of patients, (b) the results of the CN for a group of patients with brain tumors only, (c) the results of the CN for a group of patients with sinus or nasal cavity tumors.
Figure 5. Boxplots illustrating the changes in the CN after applying the adaptation procedure was applied. The boxes extend from the Q1 to Q3 quartile values of the data, with a line at the median value. Red boxes represent the results obtained after recalculating the initial plan into control scans. Green boxes represent the results obtained after adapting the plan to the patient’s current anatomical conditions (replan); (a) the results of the CN for the whole group of patients, (b) the results of the CN for a group of patients with brain tumors only, (c) the results of the CN for a group of patients with sinus or nasal cavity tumors.
Cancers 16 03128 g005
Figure 6. Boxplots illustrating the DVH parameters ( D m e a n , D 98 % , V 95 % ) within the PTV for patients with tumors located in the brain (a) and patients with cancer within the sinus and nasal areas (b). The boxes extend from the Q1 to Q3 quartile values of the data, with a line at the median value.
Figure 6. Boxplots illustrating the DVH parameters ( D m e a n , D 98 % , V 95 % ) within the PTV for patients with tumors located in the brain (a) and patients with cancer within the sinus and nasal areas (b). The boxes extend from the Q1 to Q3 quartile values of the data, with a line at the median value.
Cancers 16 03128 g006
Table 1. The overall statistics of patients treated at CCB with proton beam scanning by year.
Table 1. The overall statistics of patients treated at CCB with proton beam scanning by year.
YearNumber of Patients
with APT
Total Number
of Patients
Treated in CCB
Fraction of
APT Patients
2020 14 133 10.5%
2021 8 171 4.7%
2022 20 177 11.3%
2023 38 312 12.2%
mid.2024 36 200 18.0%
sum:11699311.7%
Table 2. Description of the parameters employed in the evaluation of the plans.
Table 2. Description of the parameters employed in the evaluation of the plans.
ParameterDescriptionEvaluated for
D max , D mean [Gy]Maximal/mean dose within the organ structureBrain, chiasm, optic nerves, brainstem, spinal cord, pituitary, cochleas and PTV
D near _ max [Gy]Dose delivered to 0.1 cm 3 of a structure volumeBrain, chiasm, optic nerves, brainstem, spinal cord, pituitary, cochleas
V 95 % [%]Volume of organ receiving at least 95% of prescribed dosePTV
D 98 % [Gy]Dose received by 98% of the organ volumePTV
Table 3. Statistics on the initiation of adaptive plan irradiation for sinus/nasal cavity and brain tumor patient groups in proton radiotherapy.
Table 3. Statistics on the initiation of adaptive plan irradiation for sinus/nasal cavity and brain tumor patient groups in proton radiotherapy.
Tumor
Location
Group
No of Patients Started Irradiation of the Adapted Plan at: Sum
of
Patients
Median
Week
1
Week
2
Week
3
Week
4
Week
5
Week
6
Week
7
Sinus/nasal
cavity
2 1 1 2 1 2 1 10 week 4
Brain 2 1 0 2 6 0 0 11 week 5
Table 4. A catalog of OARs and PTV parameters assessed in the analysis, along with the number of sinus/nasal area cancer patients exhibiting improvement, worsening or no change, considering the uncertainty of accuracy in parameter determination. The final column includes the outcomes from the nonparametric statistical hypothesis test.
Table 4. A catalog of OARs and PTV parameters assessed in the analysis, along with the number of sinus/nasal area cancer patients exhibiting improvement, worsening or no change, considering the uncertainty of accuracy in parameter determination. The final column includes the outcomes from the nonparametric statistical hypothesis test.
OAR or PTV
Parameter
No. of PatientsWilcoxon
Signed-Rank Test
ImprovingWorseningNo
Significant
Change
p-ValueTest
Result *
Brain D max 7300.064
D near _ max 4420.432
Brainstem D max 4511.000
D near _ max 5321.000
Chiasm D max 4420.846
D near _ max 5320.375
Optic Nerve L D max 6300.084
D near _ max 5400.432
Optic Nerve R D max 6301.000
D near _ max 5220.359
Spinal Cord D max 6210.301
D near _ max 4410.820
PTV D mean 8020.006
D 98 % 8200.020
V 95 % 8200.014
CI8200.027
CN8110.049
HI9100.004
* ✓—statistically significant improvement, ✗—no statistical evidence of improvement.
Table 5. A catalog of OARs and PTV parameters assessed in the analysis, along with the number of brain area tumor patients exhibiting improvement, worsening or no change, considering the uncertainty of accuracy in parameter determination. The final column includes the outcomes from the nonparametric statistical hypothesis test.
Table 5. A catalog of OARs and PTV parameters assessed in the analysis, along with the number of brain area tumor patients exhibiting improvement, worsening or no change, considering the uncertainty of accuracy in parameter determination. The final column includes the outcomes from the nonparametric statistical hypothesis test.
OAR or PTV
Parameter
No. of PatientsWilcoxon
Signed-Rank Test
ImprovingWorseningNo
Significant
Change
p-ValueTest
Result *
Brainstem D max 5150.054
D near _ max 5240.320
Chiasm D max 3530.278
D near _ max 5330.508
Optic Nerve L D max 2440.492
D near _ max 3250.767
Optic Nerve R D max 5230.275
D near _ max 3430.846
Pituitary D max 5400.496
D near _ max 6300.164
Cochlea L D max 4320.250
D near _ max 6030.012
Cochlea R D max 6300.301
D near _ max 6300.250
PTV D mean 7130.365
D 98 % 9200.275
V 95 % 9110.019
CI8030.007
CN7130.017
HI9200.175
* ✓—statistically significant improvement, ✗—no statistical evidence of improvement.
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Rydygier, M.; Skóra, T.; Kisielewicz, K.; Spaleniak, A.; Garbacz, M.; Lipa, M.; Foltyńska, G.; Góra, E.; Gajewski, J.; Krzempek, D.; et al. Proton Therapy Adaptation of Perisinusoidal and Brain Areas in the Cyclotron Centre Bronowice in Krakow: A Dosimetric Analysis. Cancers 2024, 16, 3128. https://doi.org/10.3390/cancers16183128

AMA Style

Rydygier M, Skóra T, Kisielewicz K, Spaleniak A, Garbacz M, Lipa M, Foltyńska G, Góra E, Gajewski J, Krzempek D, et al. Proton Therapy Adaptation of Perisinusoidal and Brain Areas in the Cyclotron Centre Bronowice in Krakow: A Dosimetric Analysis. Cancers. 2024; 16(18):3128. https://doi.org/10.3390/cancers16183128

Chicago/Turabian Style

Rydygier, Marzena, Tomasz Skóra, Kamil Kisielewicz, Anna Spaleniak, Magdalena Garbacz, Monika Lipa, Gabriela Foltyńska, Eleonora Góra, Jan Gajewski, Dawid Krzempek, and et al. 2024. "Proton Therapy Adaptation of Perisinusoidal and Brain Areas in the Cyclotron Centre Bronowice in Krakow: A Dosimetric Analysis" Cancers 16, no. 18: 3128. https://doi.org/10.3390/cancers16183128

APA Style

Rydygier, M., Skóra, T., Kisielewicz, K., Spaleniak, A., Garbacz, M., Lipa, M., Foltyńska, G., Góra, E., Gajewski, J., Krzempek, D., Kopeć, R., & Ruciński, A. (2024). Proton Therapy Adaptation of Perisinusoidal and Brain Areas in the Cyclotron Centre Bronowice in Krakow: A Dosimetric Analysis. Cancers, 16(18), 3128. https://doi.org/10.3390/cancers16183128

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop